PythiaDauFilter Class Reference

#include <PythiaDauFilter.h>

Inheritance diagram for PythiaDauFilter:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)
	PythiaDauFilter (const edm::ParameterSet &)
	~PythiaDauFilter ()
Public Member Functions inherited from edm::EDFilter
	EDFilter ()
virtual	~EDFilter ()
Public Member Functions inherited from edm::ProducerBase
	ProducerBase ()
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
boost::function< void(const BranchDescription &)>	registrationCallback () const
	used by the fwk to register list of products More...
virtual	~ProducerBase ()

Private Attributes
bool	chargeconju
std::vector< int >	dauIDs
std::string	label_
double	maxetacut
double	maxptcut
double	minetacut
double	minptcut
int	ndaughters
int	particleID

Additional Inherited Members
Public Types inherited from edm::EDFilter
typedef EDFilter	ModuleType
typedef WorkerT< EDFilter >	WorkerType
Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
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::EDFilter
CurrentProcessingContext const *	currentContext () const
Protected Member Functions inherited from edm::ProducerBase
template<class TProducer , class TMethod >
void	callWhenNewProductsRegistered (TProducer *iProd, TMethod iMethod)

Detailed Description

Description: Filter events using MotherId and ChildrenIds infos

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

Definition at line 40 of file PythiaDauFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file PythiaDauFilter.cc.

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

                                                               :
label_(iConfig.getUntrackedParameter("moduleLabel",std::string("generator"))),
particleID(iConfig.getUntrackedParameter("ParticleID", 0)),
chargeconju(iConfig.getUntrackedParameter("ChargeConjugation", true)),
ndaughters(iConfig.getUntrackedParameter("NumberDaughters", 0)),
minptcut(iConfig.getUntrackedParameter("MinPt", 0.)),
maxptcut(iConfig.getUntrackedParameter("MaxPt", 14000.)),
minetacut(iConfig.getUntrackedParameter("MinEta", -10.)),
maxetacut(iConfig.getUntrackedParameter("MaxEta", 10.))
{
   //now do what ever initialization is needed
   vector<int> defdauID;
   defdauID.push_back(0);
   dauIDs = iConfig.getUntrackedParameter< vector<int> >("DaughterIDs",defdauID);
}

PythiaDauFilter::~PythiaDauFilter

(

)

Definition at line 30 of file PythiaDauFilter.cc.

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

}

Member Function Documentation

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

virtual

Implements edm::EDFilter.

Definition at line 44 of file PythiaDauFilter.cc.

References chargeconju, dauIDs, edm::Event::getByLabel(), i, label_, maxetacut, maxptcut, minetacut, minptcut, ndaughters, AlCaHLTBitMon_ParallelJobs::p, particleID, and PYCOMP.

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

   const HepMC::GenEvent * myGenEvent = evt->GetEvent();
    
   
   for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
     p != myGenEvent->particles_end(); ++p ) {
     
     if( (*p)->pdg_id() != particleID ) continue ;
     int ndauac = 0;
     int ndau = 0;     
     if ( (*p)->end_vertex() ) {    
       for ( HepMC::GenVertex::particle_iterator 
           des=(*p)->end_vertex()->particles_begin(HepMC::children);
         des != (*p)->end_vertex()->particles_end(HepMC::children);
         ++des ) {
     ++ndau;       
     for( unsigned int i=0; i<dauIDs.size(); ++i) {
       if( (*des)->pdg_id() != dauIDs[i] ) continue ;
       if(   (*des)->momentum().perp() >  minptcut  &&
         (*des)->momentum().perp() <  maxptcut  &&
         (*des)->momentum().eta()  >  minetacut && 
         (*des)->momentum().eta()  <  maxetacut ) {
         ++ndauac;
         break;
       } 
     }                   
       }
     }  
     if( ndau ==  ndaughters && ndauac == ndaughters ) {
       accepted = true;
       break;
     }    
     
   }
   
   
   if( !accepted && chargeconju ) {
     
     for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
       p != myGenEvent->particles_end(); ++p ) {
       
       if( (*p)->pdg_id() != -particleID ) continue ;
       int ndauac = 0;
       int ndau = 0;     
       if ( (*p)->end_vertex() ) {
     for ( HepMC::GenVertex::particle_iterator 
         des=(*p)->end_vertex()->particles_begin(HepMC::children);
           des != (*p)->end_vertex()->particles_end(HepMC::children);
           ++des ) {
       ++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( (*des)->pdg_id() != IDanti ) continue ;
         if(   (*des)->momentum().perp() >  minptcut  &&
           (*des)->momentum().perp() <  maxptcut  &&
           (*des)->momentum().eta()  >  minetacut && 
           (*des)->momentum().eta()  <  maxetacut ) {
           ++ndauac;
           break;
         } 
       }                     
     }
       }
       if( ndau ==  ndaughters && ndauac == ndaughters ) {
     accepted = true;
     break;
       }    
     }
     
   }    

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

}

Member Data Documentation

bool PythiaDauFilter::chargeconju

private

Definition at line 55 of file PythiaDauFilter.h.

Referenced by filter().

std::vector<int> PythiaDauFilter::dauIDs

private

Definition at line 53 of file PythiaDauFilter.h.

Referenced by filter(), and PythiaDauFilter().

std::string PythiaDauFilter::label_

private

Definition at line 52 of file PythiaDauFilter.h.

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

double PythiaDauFilter::maxetacut

private

Definition at line 60 of file PythiaDauFilter.h.

Referenced by filter().

double PythiaDauFilter::maxptcut

private

Definition at line 58 of file PythiaDauFilter.h.

Referenced by filter().

double PythiaDauFilter::minetacut

private

Definition at line 59 of file PythiaDauFilter.h.

Referenced by filter().

double PythiaDauFilter::minptcut

private

Definition at line 57 of file PythiaDauFilter.h.

Referenced by filter().

int PythiaDauFilter::ndaughters

private

Definition at line 56 of file PythiaDauFilter.h.

Referenced by filter().

int PythiaDauFilter::particleID

private

Definition at line 54 of file PythiaDauFilter.h.

Referenced by filter().