PythiaMomDauFilter Class Reference

#include <PythiaMomDauFilter.h>

Inheritance diagram for PythiaMomDauFilter:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)
	PythiaMomDauFilter (const edm::ParameterSet &)
	~PythiaMomDauFilter ()
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 Attributes
bool	chargeconju
int	daughterID
std::vector< int >	dauIDs
std::vector< int >	desIDs
std::string	label_
double	maxetacut
double	maxptcut
double	minetacut
double	minptcut
double	mom_maxetacut
double	mom_maxptcut
double	mom_minetacut
double	mom_minptcut
int	ndaughters
int	ndescendants
int	particleID

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

Description: Filter events using MotherId and ChildrenIds infos

Implementation: <Notes on="" implementation>="">

Definition at line 40 of file PythiaMomDauFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file PythiaMomDauFilter.cc.

References dauIDs, desIDs, and edm::ParameterSet::getUntrackedParameter().

                                                                     :
label_(iConfig.getUntrackedParameter("moduleLabel",std::string("generator"))),
particleID(iConfig.getUntrackedParameter("ParticleID", 0)),
daughterID(iConfig.getUntrackedParameter("DaughterID", 0)),
chargeconju(iConfig.getUntrackedParameter("ChargeConjugation", true)),
ndaughters(iConfig.getUntrackedParameter("NumberDaughters", 0)),
ndescendants(iConfig.getUntrackedParameter("NumberDescendants", 0)),
minptcut(iConfig.getUntrackedParameter("MinPt", 0.)),
maxptcut(iConfig.getUntrackedParameter("MaxPt", 14000.)),
minetacut(iConfig.getUntrackedParameter("MinEta", -10.)),
maxetacut(iConfig.getUntrackedParameter("MaxEta", 10.)),
mom_minptcut(iConfig.getUntrackedParameter("MomMinPt", 0.)),
mom_maxptcut(iConfig.getUntrackedParameter("MomMaxPt", 14000.)),
mom_minetacut(iConfig.getUntrackedParameter("MomMinEta", -10.)),
mom_maxetacut(iConfig.getUntrackedParameter("MomMaxEta", 10.))
{
   //now do what ever initialization is needed
   vector<int> defdauID;
   defdauID.push_back(0);
   vector<int> defdesID;
   defdesID.push_back(0);
   dauIDs = iConfig.getUntrackedParameter< vector<int> >("DaughterIDs",    defdauID);
   desIDs = iConfig.getUntrackedParameter< vector<int> >("DescendantsIDs", defdesID);
}

PythiaMomDauFilter::~PythiaMomDauFilter

(

)

Definition at line 39 of file PythiaMomDauFilter.cc.

{
 
   // do anything here that needs to be done at desctruction time
   // (e.g. close files, deallocate resources etc.)

}

Member Function Documentation

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

virtual

Implements edm::EDFilter.

Definition at line 53 of file PythiaMomDauFilter.cc.

References chargeconju, daughterID, dauIDs, desIDs, edm::Event::getByLabel(), i, label_, maxetacut, maxptcut, minetacut, minptcut, mom_maxetacut, mom_maxptcut, mom_minetacut, mom_minptcut, ndaughters, ndescendants, AlCaHLTBitMon_ParallelJobs::p, particleID, and PYCOMP.

{
 using namespace edm;
 bool accepted = false;
 bool mom_accepted = false;
 Handle<HepMCProduct> evt;
 iEvent.getByLabel(label_, evt);

 const HepMC::GenEvent * myGenEvent = evt->GetEvent();
 int ndauac = 0;
 int ndau = 0;
 int ndesac = 0;
 int ndes   = 0;   
   
 for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();  p != myGenEvent->particles_end(); ++p ) {
     
  if( (*p)->pdg_id() != particleID ) continue ;
  if( (*p)->momentum().perp() >  mom_minptcut  && (*p)->momentum().perp() <  mom_maxptcut  && (*p)->momentum().eta()  >  mom_minetacut &&  (*p)->momentum().eta()  <  mom_maxetacut ){ 
   mom_accepted = true;
   ndauac = 0;
   ndau   = 0;
   ndesac = 0;
   ndes   = 0;     
   if ( (*p)->end_vertex() ) {  
    for ( HepMC::GenVertex::particle_iterator dau  =(*p)->end_vertex()->particles_begin(HepMC::children); dau != (*p)->end_vertex()->particles_end(HepMC::children); ++dau ) {
      ++ndau;
      for( unsigned int i=0; i<dauIDs.size(); ++i) {
       if( (*dau)->pdg_id() != dauIDs[i] ) continue ;
        ++ndauac;
        break;
     }
      if((*dau)->pdg_id() == daughterID){
       for(HepMC::GenVertex::particle_iterator des  = (*dau)->end_vertex()->particles_begin(HepMC::children); des != (*des)->end_vertex()->particles_end(HepMC::children); ++des ){
        ++ndes;        
        for( unsigned int i=0; i<desIDs.size(); ++i) {
         if( (*des)->pdg_id() != desIDs[i] ) continue ;
         if( (*des)->momentum().perp() >  minptcut  && (*des)->momentum().perp() <  maxptcut  && (*des)->momentum().eta()  >  minetacut &&  (*des)->momentum().eta()  <  maxetacut ) {
          ++ndesac;
          break;
      } 
     }
    }                    
       }
      } 
     }
    }  
    if( ndau ==  ndaughters && ndauac == ndaughters && mom_accepted && ndes == ndescendants && ndesac == ndescendants ) {
       accepted = true;
       break;
     }    
     
   }
   
   
   if( !accepted && chargeconju ) {
   
     mom_accepted = false;
     for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end(); ++p ) {
       
       if( (*p)->pdg_id() != -particleID ) continue ;
       if( (*p)->momentum().perp() >  mom_minptcut  && (*p)->momentum().perp() <  mom_maxptcut  &&  (*p)->momentum().eta()  >  mom_minetacut && (*p)->momentum().eta()  <  mom_maxetacut ){ 
       mom_accepted = true;
       ndauac = 0;
       ndau = 0;
       ndesac = 0;
       ndes   = 0;      
       if ( (*p)->end_vertex() ) {
     for ( HepMC::GenVertex::particle_iterator dau  =(*p)->end_vertex()->particles_begin(HepMC::children); dau != (*p)->end_vertex()->particles_end(HepMC::children);  ++dau ) {
       ++ndau;
       for( unsigned int i=0; i<dauIDs.size(); ++i) {
            int IDanti = -dauIDs[i];
            int pythiaCode = PYCOMP(dauIDs[i]);
            int has_antipart = pydat2.kchg[3-1][pythiaCode-1];
            if( has_antipart == 0 ) IDanti = dauIDs[i];
            if( (*dau)->pdg_id() != IDanti ) continue ;
            ++ndauac;
            break;
           }
       if((*dau)->pdg_id() == -daughterID) {
            for( HepMC::GenVertex::particle_iterator des  = (*dau)->end_vertex()->particles_begin(HepMC::children); des != (*des)->end_vertex()->particles_end(HepMC::children);  ++des ){
                 ++ndes;
        for( unsigned int i=0; i<desIDs.size(); ++i) {
         int IDanti = -desIDs[i];
         int pythiaCode = PYCOMP(desIDs[i]);
         int has_antipart = pydat2.kchg[3-1][pythiaCode-1];
         if( has_antipart == 0 ) IDanti = desIDs[i];
         if( (*des)->pdg_id() != IDanti ) continue ;
         if(   (*des)->momentum().perp() >  minptcut  && (*des)->momentum().perp() <  maxptcut  &&  (*des)->momentum().eta()  >  minetacut && (*des)->momentum().eta()  <  maxetacut ) {
           ++ndesac;
           break;
         } 
       }                     
          }  
     }
    } 
       }
      }
       if( ndau ==  ndaughters && ndauac == ndaughters && mom_accepted && ndes == ndescendants && ndesac == ndescendants ) {
     accepted = true;
     break;
       }    
     }
     
   }    

   if (accepted){
   return true; } else {return false;}

}

Member Data Documentation

bool PythiaMomDauFilter::chargeconju

private

Definition at line 57 of file PythiaMomDauFilter.h.

Referenced by filter().

int PythiaMomDauFilter::daughterID

private

Definition at line 56 of file PythiaMomDauFilter.h.

Referenced by filter().

std::vector<int> PythiaMomDauFilter::dauIDs

private

Definition at line 53 of file PythiaMomDauFilter.h.

Referenced by filter(), and PythiaMomDauFilter().

std::vector<int> PythiaMomDauFilter::desIDs

private

Definition at line 54 of file PythiaMomDauFilter.h.

Referenced by filter(), and PythiaMomDauFilter().

std::string PythiaMomDauFilter::label_

private

Definition at line 52 of file PythiaMomDauFilter.h.

Referenced by Mixins._Labelable::_findDependencies(), filter(), and Mixins._Labelable::setLabel().

double PythiaMomDauFilter::maxetacut

private

Definition at line 63 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::maxptcut

private

Definition at line 61 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::minetacut

private

Definition at line 62 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::minptcut

private

Definition at line 60 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::mom_maxetacut

private

Definition at line 67 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::mom_maxptcut

private

Definition at line 65 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::mom_minetacut

private

Definition at line 66 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::mom_minptcut

private

Definition at line 64 of file PythiaMomDauFilter.h.

Referenced by filter().

int PythiaMomDauFilter::ndaughters

private

Definition at line 58 of file PythiaMomDauFilter.h.

Referenced by filter().

int PythiaMomDauFilter::ndescendants

private

Definition at line 59 of file PythiaMomDauFilter.h.

Referenced by filter().

int PythiaMomDauFilter::particleID

private

Definition at line 55 of file PythiaMomDauFilter.h.

Referenced by filter().