#include <PythiaFilterZgamma.h>

Inheritance diagram for PythiaFilterZgamma:

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

	PythiaFilterZgamma (const edm::ParameterSet &)

	~PythiaFilterZgamma () override

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

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

	~EDFilter () override

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

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

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

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) override

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

void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const

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

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

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

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

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

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
double	etaElMax

double	etaMuMax

double	etaPhotonMax

double	ptElMin

double	ptMuMin

double	ptPhotonMin

int	selProc

int	theNumberOfSelected

edm::EDGetTokenT< edm::HepMCProduct >	token_

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

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

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 &)

static bool	wantsGlobalLuminosityBlocks ()

static bool	wantsGlobalRuns ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

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

PythiaFilterZgamma filter implements generator-level preselections for Z0 + photon like events to be used in TGC.

Author: A.Kyriakis, NCSR "Demokritos"

Definition at line 29 of file PythiaFilterZgamma.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 8 of file PythiaFilterZgamma.cc.

References theNumberOfSelected.

                                                                      :
 token_(consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared"))),
 selProc(iConfig.getUntrackedParameter<int>("SelectProcess")),
 ptElMin(iConfig.getUntrackedParameter<double>("MinElPt", 5.0)),
 ptMuMin(iConfig.getUntrackedParameter<double>("MinMuPt", 3.0)),
 ptPhotonMin(iConfig.getUntrackedParameter<double>("MinPhotPt", 5.0)),
 etaElMax(iConfig.getUntrackedParameter<double>("MaxElecEta", 2.7)),
 etaMuMax(iConfig.getUntrackedParameter<double>("MaxMuonEta", 2.4)),
 etaPhotonMax(iConfig.getUntrackedParameter<double>("MaxPhotEta", 2.7))
 {  
   theNumberOfSelected = 0;
 }

PythiaFilterZgamma::~PythiaFilterZgamma ( )

override

Definition at line 22 of file PythiaFilterZgamma.cc.

22 {}

Member Function Documentation

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

override

Definition at line 26 of file PythiaFilterZgamma.cc.

References funct::abs(), cms::dd::accepted(), etaElMax, etaMuMax, etaPhotonMax, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), RPCpg::mu, AlCaHLTBitMon_ParallelJobs::p, ptElMin, ptMuMin, ptPhotonMin, selProc, theNumberOfSelected, and token_.

                                                                             {
 
 
   bool accepted = false;
   edm::Handle<edm::HepMCProduct> evt;
   iEvent.getByToken(token_, evt);
 
   const HepMC::GenEvent * myGenEvent = evt->GetEvent();
 
   std::vector<const HepMC::GenParticle *> el;
   std::vector<const HepMC::GenParticle *> mu;
   std::vector<const HepMC::GenParticle *> gam;
 
   for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();  
                p != myGenEvent->particles_end(); ++p ) {
               
     if(selProc == 1 ) {
       if ( std::abs((*p)->pdg_id())==11 && (*p)->status()==1 )
         el.push_back(*p);
       if(el.size()>1) break;
     } else if(selProc == 2 ) {
       if ( std::abs((*p)->pdg_id())==13 && (*p)->status()==1 )
         mu.push_back(*p);
       if(mu.size()>1) break;        
     }
       
   } // end of first particle loop for finding Z0 decays
 
   for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();  
                p != myGenEvent->particles_end(); ++p ) {
 
     if ( std::abs((*p)->pdg_id())==22 && (*p)->status()==1 )
       gam.push_back(*p);
     if(gam.size()==1) break;
 
   } // end of second particle loop for finding gamma
 
   if(selProc == 1 ) {
     if(el.size()>1){
       double ptEl1 = el[0]->momentum().perp();
       double ptEl2 = el[1]->momentum().perp();
       double etaEl1 = el[0]->momentum().eta();
       double etaEl2 = el[1]->momentum().eta();
       if (ptEl1 > ptElMin && ptEl2 > ptElMin  && 
     std::abs(etaEl1) < etaElMax &&
     std::abs(etaEl2) < etaElMax) {
         if(gam.size()==1) {
       double ptGam = gam[0]->momentum().perp();
           double etaGam = gam[0]->momentum().eta();
       if(ptGam > ptPhotonMin && std::abs(etaGam) < etaPhotonMax)
                  accepted=true;
     }     
       } 
     } 
   } else if(selProc == 2 ) {
 
     if(mu.size()>1){
       double ptMu1 = mu[0]->momentum().perp();
       double ptMu2 = mu[1]->momentum().perp();
       double etaMu1 = mu[0]->momentum().eta();
       double etaMu2 = mu[1]->momentum().eta();
       if (ptMu1 > ptMuMin && ptMu2 > ptMuMin  && 
     std::abs(etaMu1) < etaMuMax &&
     std::abs(etaMu2) < etaMuMax) {
         if(gam.size()==1) {
       double ptGam = gam[0]->momentum().perp();
           double etaGam = gam[0]->momentum().eta();
       if(ptGam > ptPhotonMin && std::abs(etaGam) < etaPhotonMax)
                  accepted=true;
     }     
       } 
     } 
     
   }
 /*
   if(accepted) {
     std::cout << "Accepted event Number: " << theNumberOfSelected 
               << "  of category " << selProc << std::endl;
     if(selProc == 1 ) { 
       std::cout << "Electon 1 pt, eta = " << el[0]->momentum().perp() << ", " 
                 <<  el[0]->momentum().eta() << std::endl;
       std::cout << "Electon 2 pt, eta = " << el[1]->momentum().perp() << ", " 
                 <<  el[1]->momentum().eta() << std::endl;
       std::cout << "Photon pt, eta = " << gam[0]->momentum().perp() << ", " 
                 <<  gam[0]->momentum().eta() << std::endl;
     } else if(selProc == 2 ) {
       std::cout << "Muon 1 pt, eta = " << mu[0]->momentum().perp() << ", " 
                 <<  mu[0]->momentum().eta() << std::endl;
       std::cout << "Muon 2 pt, eta = " << mu[1]->momentum().perp() << ", " 
                 <<  mu[1]->momentum().eta() << std::endl;
       std::cout << "Photon pt, eta = " << gam[0]->momentum().perp() << ", " 
                 <<  gam[0]->momentum().eta() << std::endl;
 
     }  
   }
 */
   if (accepted) {
     theNumberOfSelected++;
     return true; 
   }
   else return false;
 
 }

Member Data Documentation

double PythiaFilterZgamma::etaElMax

private

Definition at line 46 of file PythiaFilterZgamma.h.

Referenced by filter().

double PythiaFilterZgamma::etaMuMax

private

Definition at line 47 of file PythiaFilterZgamma.h.

Referenced by filter().

double PythiaFilterZgamma::etaPhotonMax

private

Definition at line 48 of file PythiaFilterZgamma.h.

Referenced by filter().

double PythiaFilterZgamma::ptElMin

private

Definition at line 42 of file PythiaFilterZgamma.h.

Referenced by filter().

double PythiaFilterZgamma::ptMuMin

private

Definition at line 43 of file PythiaFilterZgamma.h.

Referenced by filter().

double PythiaFilterZgamma::ptPhotonMin

private

Definition at line 44 of file PythiaFilterZgamma.h.

Referenced by filter().

int PythiaFilterZgamma::selProc

private

Definition at line 40 of file PythiaFilterZgamma.h.

Referenced by filter().

int PythiaFilterZgamma::theNumberOfSelected

private

Definition at line 50 of file PythiaFilterZgamma.h.

Referenced by filter(), and PythiaFilterZgamma().

edm::EDGetTokenT<edm::HepMCProduct> PythiaFilterZgamma::token_

private

Definition at line 38 of file PythiaFilterZgamma.h.

Referenced by filter().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation