PythiaFilter Class Reference

#include <IOMC/PythiaFilter/src/PythiaFilter.cc>

Inheritance diagram for PythiaFilter:

Public Member Functions
virtual bool	filter (edm::Event &, const edm::EventSetup &)
	PythiaFilter (const edm::ParameterSet &)
	~PythiaFilter ()
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
double	maxetacut
double	maxpcut
double	maxphicut
double	maxptcut
double	maxrapcut
double	minetacut
double	minpcut
double	minphicut
double	minptcut
double	minrapcut
int	motherID
int	particleID
int	processID
double	rapidity
int	status
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

Description: <one line="" class="" summary>="">

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

Definition at line 42 of file PythiaFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 11 of file PythiaFilter.cc.

                                                         :
token_(consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared"))),
particleID(iConfig.getUntrackedParameter("ParticleID", 0)),
minpcut(iConfig.getUntrackedParameter("MinP", 0.)),
maxpcut(iConfig.getUntrackedParameter("MaxP", 10000.)),
minptcut(iConfig.getUntrackedParameter("MinPt", 0.)),
maxptcut(iConfig.getUntrackedParameter("MaxPt", 10000.)),
minetacut(iConfig.getUntrackedParameter("MinEta", -10.)),
maxetacut(iConfig.getUntrackedParameter("MaxEta", 10.)),
minrapcut(iConfig.getUntrackedParameter("MinRapidity", -20.)),
maxrapcut(iConfig.getUntrackedParameter("MaxRapidity", 20.)),
minphicut(iConfig.getUntrackedParameter("MinPhi", -3.5)),
maxphicut(iConfig.getUntrackedParameter("MaxPhi", 3.5)),
status(iConfig.getUntrackedParameter("Status", 0)),
motherID(iConfig.getUntrackedParameter("MotherID", 0)),
processID(iConfig.getUntrackedParameter("ProcessID", 0))
{
   //now do what ever initialization is needed

}

PythiaFilter::~PythiaFilter

(

)

Definition at line 33 of file PythiaFilter.cc.

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

}

Member Function Documentation

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

virtual

Implements edm::EDFilter.

Definition at line 47 of file PythiaFilter.cc.

References funct::abs(), GenParticle::GenParticle, edm::Event::getByToken(), cmsBatch::log, maxetacut, maxpcut, maxphicut, maxptcut, maxrapcut, minetacut, minpcut, minphicut, minptcut, minrapcut, motherID, AlCaHLTBitMon_ParallelJobs::p, particleID, processID, rapidity, status, and token_.

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

   const HepMC::GenEvent * myGenEvent = evt->GetEvent();
    
   if(processID == 0 || processID == myGenEvent->signal_process_id()) {
     
     for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
       p != myGenEvent->particles_end(); ++p ) {
       
       rapidity = 0.5*log( ((*p)->momentum().e()+(*p)->momentum().pz()) / ((*p)->momentum().e()-(*p)->momentum().pz()) );
       
       if ( abs((*p)->pdg_id()) == particleID 
        && (*p)->momentum().rho() > minpcut 
        && (*p)->momentum().rho() < maxpcut
        && (*p)->momentum().perp() > minptcut 
        && (*p)->momentum().perp() < maxptcut
        && (*p)->momentum().eta() > minetacut
        && (*p)->momentum().eta() < maxetacut 
        && rapidity > minrapcut
        && rapidity < maxrapcut 
        && (*p)->momentum().phi() > minphicut
        && (*p)->momentum().phi() < maxphicut ) {
     
         
     
     if (status == 0 && motherID == 0){
       accepted = true;
     }
     if (status != 0 && motherID == 0){
       if ((*p)->status() == status)   
         accepted = true;
     }
     
     HepMC::GenParticle* mother = (*((*p)->production_vertex()->particles_in_const_begin()));
     
     if (status == 0 && motherID != 0){    
       if (abs(mother->pdg_id()) == abs(motherID)) {
         accepted = true;
       }
     }
     if (status != 0 && motherID != 0){
       
       if ((*p)->status() == status && abs(mother->pdg_id()) == abs(motherID)){   
         accepted = true;
         
       }
     }
     
     /*
       if (status == 0 && motherID != 0){    
       if (abs(((*p)->mother())->pdg_id()) == abs(motherID)) {
       accepted = true;
       }
       }
       if (status != 0 && motherID != 0){
           
       if ((*p)->status() == status && abs(((*p)->mother())->pdg_id()) == abs(motherID)){   
       accepted = true;
       
       }
       }
     */
     
       }    
     }
     
   } else { accepted = true; }

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

}

Member Data Documentation

double PythiaFilter::maxetacut

private

Definition at line 59 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::maxpcut

private

Definition at line 55 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::maxphicut

private

Definition at line 63 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::maxptcut

private

Definition at line 57 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::maxrapcut

private

Definition at line 61 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::minetacut

private

Definition at line 58 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::minpcut

private

Definition at line 54 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::minphicut

private

Definition at line 62 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::minptcut

private

Definition at line 56 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::minrapcut

private

Definition at line 60 of file PythiaFilter.h.

Referenced by filter().

int PythiaFilter::motherID

private

Definition at line 68 of file PythiaFilter.h.

Referenced by filter().

int PythiaFilter::particleID

private

Definition at line 53 of file PythiaFilter.h.

Referenced by filter().

int PythiaFilter::processID

private

Definition at line 69 of file PythiaFilter.h.

Referenced by filter().

double PythiaFilter::rapidity

private

Definition at line 65 of file PythiaFilter.h.

Referenced by filter().

int PythiaFilter::status

private

Definition at line 67 of file PythiaFilter.h.

Referenced by dirstructure.Comparison::__make_image(), dirstructure.Comparison::__repr__(), and filter().

edm::EDGetTokenT<edm::HepMCProduct> PythiaFilter::token_

private

Definition at line 52 of file PythiaFilter.h.

Referenced by filter().