#include <davidlw/HighMultiplicityGenFilter/src/HighMultiplicityGenFilter.cc>

Inheritance diagram for HighMultiplicityGenFilter:

Public Member Functions
	HighMultiplicityGenFilter (const edm::ParameterSet &)

	~HighMultiplicityGenFilter ()

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 Member Functions
virtual void	beginJob ()

virtual void	endJob ()

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

Private Attributes
double	etaMax

int	nAccepted

int	nMin

edm::ESHandle< ParticleDataTable >	pdt

double	ptMin

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: <one line="" class="" summary>="">

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

Definition at line 38 of file HighMultiplicityGenFilter.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 67 of file HighMultiplicityGenFilter.cc.

References nAccepted.

                                                                                    :
 etaMax(iConfig.getUntrackedParameter<double>("etaMax")),
 ptMin(iConfig.getUntrackedParameter<double>("ptMin")),
 nMin(iConfig.getUntrackedParameter<int>("nMin"))
 {
   //now do what ever initialization is needed
   nAccepted = 0; 
 }

HighMultiplicityGenFilter::~HighMultiplicityGenFilter ( )

Definition at line 77 of file HighMultiplicityGenFilter.cc.

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

Member Function Documentation

void HighMultiplicityGenFilter::beginJob ( void )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 125 of file HighMultiplicityGenFilter.cc.

126 {}

void HighMultiplicityGenFilter::endJob ( void )

privatevirtual

Reimplemented from edm::EDFilter.

Definition at line 130 of file HighMultiplicityGenFilter.cc.

References gather_cfg::cout, nAccepted, and nMin.

                                   {
   std::cout<<"There are "<<nAccepted<<" events with multiplicity greater than "<<nMin<<std::endl;
 }

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

privatevirtual

Implements edm::EDFilter.

Definition at line 92 of file HighMultiplicityGenFilter.cc.

References abs, DeDxDiscriminatorTools::charge(), gather_cfg::cout, etaMax, edm::Event::getByLabel(), edm::EventSetup::getData(), nAccepted, nMin, AlCaHLTBitMon_ParallelJobs::p, pdt, evf::utils::pid, and ptMin.

 {
 
   bool accepted = false;
   edm::Handle<edm::HepMCProduct> evt;
   iEvent.getByLabel("generator", evt);
 
   iSetup.getData(pdt);
 
   const HepMC::GenEvent * myGenEvent = evt->GetEvent();
 
   int nMult=0;
   for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();   p != myGenEvent->particles_end(); ++p ) {
 
     if((*p)->status()!=1) continue;
 
     double charge = 0;
     int pid = (*p)->pdg_id();
     if(abs(pid) > 100000) { std::cout<<"pid="<<pid<<" status="<<(*p)->status()<<std::endl; continue; }
     const ParticleData* part = pdt->particle(pid);
     if(part) charge = part->charge();
     if(charge == 0) continue;
 
     if ( 
      (*p)->momentum().perp() > ptMin 
      && fabs((*p)->momentum().eta()) < etaMax  ) nMult++;
   }
   if(nMult>=nMin) { nAccepted++; accepted = true; }
   return accepted;
 }