#include <MCSingleParticleYPt.h>

Inheritance diagram for MCSingleParticleYPt:

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

	MCSingleParticleYPt (const edm::ParameterSet &)

	~MCSingleParticleYPt () 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

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)

virtual	~EDConsumerBase () noexcept(false)

Private Attributes
bool	fchekantiparticle

int	fVerbose

std::vector< int >	particleID

std::vector< double >	ptMin

double	rapidity

std::vector< double >	rapMax

std::vector< double >	rapMin

std::vector< 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
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 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

Definition at line 39 of file MCSingleParticleYPt.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 11 of file MCSingleParticleYPt.cc.

References fchekantiparticle, fVerbose, edm::ParameterSet::getUntrackedParameter(), mps_fire::i, particleID, ptMin, rapMax, rapMin, and status.

                                                                        :
 token_(consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared")))
 {
    fVerbose = iConfig.getUntrackedParameter("verbose",0);
    fchekantiparticle = iConfig.getUntrackedParameter("CheckAntiparticle",true);
    //here do whatever other initialization is needed
    vector<int> defpid ;
    defpid.push_back(0) ;
    particleID = iConfig.getUntrackedParameter< vector<int> >("ParticleID",defpid);  
    vector<double> defptmin ;
    defptmin.push_back(0.);
    ptMin = iConfig.getUntrackedParameter< vector<double> >("MinPt", defptmin);
    vector<double> defrapmin ;
    defrapmin.push_back(-10.);
    rapMin = iConfig.getUntrackedParameter< vector<double> >("MinY", defrapmin);
    vector<double> defrapmax ;
    defrapmax.push_back(10.);
    rapMax = iConfig.getUntrackedParameter< vector<double> >("MaxY", defrapmax);
    vector<int> defstat ;
    defstat.push_back(0);
    status = iConfig.getUntrackedParameter< vector<int> >("Status", defstat);
 
     // check for same size
     if ( (ptMin.size() > 1 &&  particleID.size() != ptMin.size()) 
      ||  (rapMin.size() > 1 && particleID.size() != rapMin.size()) 
      ||  (rapMax.size() > 1 && particleID.size() != rapMax.size())
      ||  (status.size() > 1 && particleID.size() != status.size()) ) {
      edm::LogWarning("MCSingleParticleYPt") << "WARNING: MCSingleParticleYPt : size of vector cuts do not match!!" << endl;
     }
 
     // if ptMin size smaller than particleID , fill up further with defaults
     if (particleID.size() > ptMin.size() ){
        vector<double> defptmin2 ;
        for (unsigned int i = 0; i < particleID.size(); i++){ defptmin2.push_back(0.);}
        ptMin = defptmin2;   
     } 
     // if etaMin size smaller than particleID , fill up further with defaults
     if (particleID.size() > rapMin.size() ){
        vector<double> defrapmin2 ;
        for (unsigned int i = 0; i < particleID.size(); i++){ defrapmin2.push_back(-10.);}
        rapMin = defrapmin2;   
     } 
     // if etaMax size smaller than particleID , fill up further with defaults
     if (particleID.size() > rapMax.size() ){
        vector<double> defrapmax2 ;
        for (unsigned int i = 0; i < particleID.size(); i++){ defrapmax2.push_back(10.);}
        rapMax = defrapmax2;   
     }     
     // if status size smaller than particleID , fill up further with defaults
     if (particleID.size() > status.size() ){
        vector<int> defstat2 ;
        for (unsigned int i = 0; i < particleID.size(); i++){ defstat2.push_back(0);}
        status = defstat2;   
     } 
 
     if (fVerbose > 0) {
        edm::LogInfo("MCSingleParticleYPt") <<     "----------------------------------------------------------------------" << std::endl;
        edm::LogInfo("MCSingleParticleYPt") <<     "----- MCSingleParticleYPt" << std::endl;
        for (unsigned int i=0; i<particleID.size(); ++i) {
            edm::LogInfo("MCSingleParticleYPt") << " ID: " <<  particleID[i] << " pT > " << ptMin[i] << ",   " << rapMin[i] << " < y < " << rapMax[i] << ",   status = " << status[i] << std::endl;
        }
        if (fchekantiparticle) edm::LogInfo("MCSingleParticleYPt") << " anti-particles will be tested as well." << std::endl;
           edm::LogInfo("MCSingleParticleYPt") <<     "----------------------------------------------------------------------" << std::endl;
     }
 }

MCSingleParticleYPt::~MCSingleParticleYPt ( )

override

Definition at line 78 of file MCSingleParticleYPt.cc.

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

Member Function Documentation

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

override

Definition at line 86 of file MCSingleParticleYPt.cc.

References fchekantiparticle, fVerbose, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), mps_fire::i, cmsBatch::log, AlCaHLTBitMon_ParallelJobs::p, particleID, ptMin, rapidity, rapMax, rapMin, status, and token_.

 {
    using namespace edm;
    bool accepted = false;
    Handle<HepMCProduct> evt;
    iEvent.getByToken(token_, evt);
 
    const HepMC::GenEvent * myGenEvent = evt->GetEvent();
    for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
      p != myGenEvent->particles_end(); ++p ) {
      if (fVerbose > 3) edm::LogInfo("MCSingleParticleYPt") << "Looking at particle : " << (*p)->pdg_id() << " status : " << (*p)->status() << std::endl;
 
      for (unsigned int i = 0; i < particleID.size(); i++) {
        if (particleID[i] == (*p)->pdg_id() || (fchekantiparticle && (-particleID[i] == (*p)->pdg_id())) || particleID[i] == 0) {
          // calculate rapidity just for the desired particle and make sure, this particles has enough energy
          rapidity = ((*p)->momentum().e()-(*p)->momentum().pz()) > 0. ? 0.5*log( ((*p)->momentum().e()+(*p)->momentum().pz()) / ((*p)->momentum().e()-(*p)->momentum().pz()) ) : rapMax[i]+.1;    
          if (fVerbose > 2) edm::LogInfo("MCSingleParticleYPt") << "Testing particle : " << (*p)->pdg_id() << " pT: " << (*p)->momentum().perp() << " y: " << rapidity << " status : " << (*p)->status() << endl;
      if ( (*p)->momentum().perp() > ptMin[i] 
               && rapidity > rapMin[i] 
               && rapidity < rapMax[i] 
               && ((*p)->status() == status[i] || status[i] == 0) ) { 
            accepted = true;
            if (fVerbose > 1) 
               edm::LogInfo("MCSingleParticleYPt") << "Accepted particle : " << (*p)->pdg_id() << " pT: " << (*p)->momentum().perp() << " y: " << rapidity << " status : " << (*p)->status() << endl; 
            break;
      }  
      
        }
      }
      if (accepted) break;
    }
    
    if (accepted) { return true; } 
    else { return false; }
 }