#include <ZgammaMassFilter.h>

Inheritance diagram for ZgammaMassFilter:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)

	ZgammaMassFilter (const edm::ParameterSet &)

	~ZgammaMassFilter ()

Public Member Functions inherited from edm::EDFilter
	EDFilter ()

ModuleDescription const &	moduleDescription () const

virtual	~EDFilter ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

virtual	~ProducerBase ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (const std::string &iProcessName, std::vector< const char * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
int	charge (const int &Id)

Private Attributes
double	maxLeptonEta

double	maxPhotonEta

double	minDileptonMass

double	minLeptonEta

double	minLeptonPt

double	minPhotonEta

double	minPhotonPt

double	minZgMass

edm::EDGetTokenT < edm::HepMCProduct >	token_

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDFilter
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 43 of file ZgammaMassFilter.h.

Constructor & Destructor Documentation

ZgammaMassFilter::ZgammaMassFilter ( const edm::ParameterSet & iConfig )

explicit

Definition at line 22 of file ZgammaMassFilter.cc.

References edm::ParameterSet::getParameter().

 {
   token_          =consumes<edm::HepMCProduct>(iConfig.getParameter<InputTag>("HepMCProduct"));
   minPhotonPt     =iConfig.getParameter<double>("minPhotonPt");
   minLeptonPt     =iConfig.getParameter<double>("minLeptonPt");
   minPhotonEta    =iConfig.getParameter<double>("minPhotonEta");
   minLeptonEta    =iConfig.getParameter<double>("minLeptonEta");
   maxPhotonEta    =iConfig.getParameter<double>("maxPhotonEta");
   maxLeptonEta    =iConfig.getParameter<double>("maxLeptonEta");
   minDileptonMass =iConfig.getParameter<double>("minDileptonMass");
   minZgMass       =iConfig.getParameter<double>("minZgMass");
 }

ZgammaMassFilter::~ZgammaMassFilter ( )

Definition at line 35 of file ZgammaMassFilter.cc.

36 {

37 }

Member Function Documentation

int ZgammaMassFilter::charge ( const int & Id )

private

bool ZgammaMassFilter::filter	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

virtual

Implements edm::EDFilter.

Definition at line 40 of file ZgammaMassFilter.cc.

References funct::abs(), edm::Event::getByToken(), AlCaHLTBitMon_ParallelJobs::p, and configurableAnalysis::Photon.

 {
   using namespace edm;
   
   bool accepted = false;
   Handle<HepMCProduct> evt;
   iEvent.getByToken(token_, evt);
   const HepMC::GenEvent * myGenEvent = evt->GetEvent();
   
   vector<TLorentzVector> Lepton; Lepton.clear();
   vector<TLorentzVector> Photon; Photon.clear();   
   vector<float> Charge; Charge.clear();
   
   for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
     p != myGenEvent->particles_end(); ++p ) 
     {
       if ((*p)->status() == 1 && (abs((*p)->pdg_id()) == 11 || abs((*p)->pdg_id()) == 13 || abs((*p)->pdg_id()) == 15)) { 
     TLorentzVector LeptP((*p)->momentum().px(), (*p)->momentum().py(), (*p)->momentum().pz(), (*p)->momentum().e()); 
     if (LeptP.Pt() > minLeptonPt) { 
       Lepton.push_back(LeptP);
     } // if pt
       }// if lepton       
       
       if ( abs((*p)->pdg_id()) == 22 && (*p)->status() == 1) {
       TLorentzVector PhotP((*p)->momentum().px(), (*p)->momentum().py(), (*p)->momentum().pz(), (*p)->momentum().e()); 
       if (PhotP.Pt() > minPhotonPt) {
       Photon.push_back(PhotP);
       }// if pt
      }// if photon
 
     }// loop over particles  
   
   
   // std::cout << "\n" << "Photon size: " << Photon.size() << std::endl;
   // for (unsigned int u=0; u<Photon.size(); u++){
   //   std::cout << "BEF photon PT: " << Photon[u].Pt() << std::endl;
   // }
   // std::cout << "\n" << "Lepton size: " << Lepton.size() << std::endl;
   // for (unsigned int u=0; u<Lepton.size(); u++){
   //   std::cout << "BEF lepton PT: " << Lepton[u].Pt() << std::endl;
   // }
 
   // order Lepton and Photon according to Pt
   std::stable_sort(Photon.begin(), Photon.end(), orderByPt());
   std::stable_sort(Lepton.begin(), Lepton.end(), orderByPt());
 
   //  std::cout << "\n" << std::endl;
   //  std::cout << "\n" << "Photon size: " << Photon.size() << std::endl;
   //  for (unsigned int u=0; u<Photon.size(); u++){
   //    std::cout << "AFT photon PT: " << Photon[u].Pt() << std::endl;
   //  }
   //  std::cout << "\n" << "Lepton size: " << Lepton.size() << std::endl;
   //  for (unsigned int u=0; u<Lepton.size(); u++){
   //    std::cout << "AFT lepton PT: " << Lepton[u].Pt() << std::endl;
   //  }
   //  std::cout << "\n" << std::endl;
   
   if (
       Photon.size() > 0 && Lepton.size() > 1 &&  
       Photon[0].Pt()  > minPhotonPt  && 
       Lepton[0].Pt()  > minLeptonPt  && 
       Lepton[1].Pt()  > minLeptonPt  &&  
       Photon[0].Eta() > minPhotonEta && 
       Lepton[0].Eta() > minLeptonEta && 
       Lepton[1].Eta() > minLeptonEta &&  
       Photon[0].Eta() < maxPhotonEta && 
       Lepton[0].Eta() < maxLeptonEta && 
       Lepton[1].Eta() < maxLeptonEta &&  
       (Lepton[0]+Lepton[1]).M()           > minDileptonMass  &&
       (Lepton[0]+Lepton[1]+Photon[0]).M() > minZgMass 
       )
     { // satisfy molteplicity, kinematics, and ll llg minimum mass
       accepted = true; 
     }
 
   //  std::cout << "++ returning: " << accepted << "\n" << std::endl;
 
   return accepted;    
 }

Member Data Documentation

double ZgammaMassFilter::maxLeptonEta

private

Definition at line 65 of file ZgammaMassFilter.h.

double ZgammaMassFilter::maxPhotonEta

private

Definition at line 64 of file ZgammaMassFilter.h.

double ZgammaMassFilter::minDileptonMass

private

Definition at line 67 of file ZgammaMassFilter.h.

double ZgammaMassFilter::minLeptonEta

private

Definition at line 62 of file ZgammaMassFilter.h.

double ZgammaMassFilter::minLeptonPt

private

Definition at line 59 of file ZgammaMassFilter.h.

double ZgammaMassFilter::minPhotonEta

private

Definition at line 61 of file ZgammaMassFilter.h.

double ZgammaMassFilter::minPhotonPt

private

Definition at line 58 of file ZgammaMassFilter.h.

double ZgammaMassFilter::minZgMass

private

Definition at line 68 of file ZgammaMassFilter.h.

edm::EDGetTokenT<edm::HepMCProduct> ZgammaMassFilter::token_

private

Definition at line 56 of file ZgammaMassFilter.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation