PythiaFilterMultiMother Class Reference

#include <PythiaFilterMultiMother.h>

Inheritance diagram for PythiaFilterMultiMother:

Public Member Functions
bool	filter (edm::Event &, const edm::EventSetup &) override
	PythiaFilterMultiMother (const edm::ParameterSet &)
	~PythiaFilterMultiMother () 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
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
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)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Attributes
double	betaBoost
double	maxetacut
double	maxpcut
double	maxphicut
double	maxptcut
double	maxrapcut
double	minetacut
double	minpcut
double	minphicut
double	minptcut
double	minrapcut
std::vector< int >	motherIDs
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
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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
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 24 of file PythiaFilterMultiMother.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 12 of file PythiaFilterMultiMother.cc.

References betaBoost, and processID.

                                                                               :
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)),
motherIDs(iConfig.getUntrackedParameter("MotherIDs", std::vector<int>{0})),
processID(iConfig.getUntrackedParameter("ProcessID", 0)),
betaBoost(iConfig.getUntrackedParameter("BetaBoost",0.))
{
   //now do what ever initialization is needed

}

PythiaFilterMultiMother::~PythiaFilterMultiMother ( )

override

Definition at line 35 of file PythiaFilterMultiMother.cc.

{

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

}

Member Function Documentation

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

override

Definition at line 49 of file PythiaFilterMultiMother.cc.

References funct::abs(), cms::dd::accepted(), betaBoost, GenParticle::GenParticle, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), cmsBatch::log, maxetacut, maxpcut, maxphicut, maxptcut, maxrapcut, minetacut, minpcut, minphicut, minptcut, minrapcut, motherIDs, AlCaHLTBitMon_ParallelJobs::p, particleID, processID, rapidity, status, token_, and MCFilterZboostHelper::zboost().

{
   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 ) {
       HepMC::FourVector mom = MCFilterZboostHelper::zboost((*p)->momentum(),betaBoost);
       rapidity = 0.5*log( (mom.e()+mom.pz()) / (mom.e()-mom.pz()) );

       if ( abs((*p)->pdg_id()) == particleID
        && mom.rho() > minpcut
        && mom.rho() < maxpcut
        && (*p)->momentum().perp() > minptcut
        && (*p)->momentum().perp() < maxptcut
        && mom.eta() > minetacut
        && mom.eta() < maxetacut
        && rapidity > minrapcut
        && rapidity < maxrapcut
        && (*p)->momentum().phi() > minphicut
        && (*p)->momentum().phi() < maxphicut ) {


     for(std::vector<int>::const_iterator motherID = motherIDs.begin(); motherID != motherIDs.end(); ++motherID) {
       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 PythiaFilterMultiMother::betaBoost

private

Definition at line 53 of file PythiaFilterMultiMother.h.

Referenced by filter(), and PythiaFilterMultiMother().

double PythiaFilterMultiMother::maxetacut

private

Definition at line 41 of file PythiaFilterMultiMother.h.

Referenced by filter().

double PythiaFilterMultiMother::maxpcut

private

Definition at line 37 of file PythiaFilterMultiMother.h.

Referenced by filter().

double PythiaFilterMultiMother::maxphicut

private

Definition at line 45 of file PythiaFilterMultiMother.h.

Referenced by filter().

double PythiaFilterMultiMother::maxptcut

private

Definition at line 39 of file PythiaFilterMultiMother.h.

Referenced by filter().

double PythiaFilterMultiMother::maxrapcut

private

Definition at line 43 of file PythiaFilterMultiMother.h.

Referenced by filter().

double PythiaFilterMultiMother::minetacut

private

Definition at line 40 of file PythiaFilterMultiMother.h.

Referenced by filter().

double PythiaFilterMultiMother::minpcut

private

Definition at line 36 of file PythiaFilterMultiMother.h.

Referenced by filter().

double PythiaFilterMultiMother::minphicut

private

Definition at line 44 of file PythiaFilterMultiMother.h.

Referenced by filter().

double PythiaFilterMultiMother::minptcut

private

Definition at line 38 of file PythiaFilterMultiMother.h.

Referenced by filter().

double PythiaFilterMultiMother::minrapcut

private

Definition at line 42 of file PythiaFilterMultiMother.h.

Referenced by filter().

std::vector<int> PythiaFilterMultiMother::motherIDs

private

Definition at line 50 of file PythiaFilterMultiMother.h.

Referenced by filter().

int PythiaFilterMultiMother::particleID

private

Definition at line 35 of file PythiaFilterMultiMother.h.

Referenced by filter().

int PythiaFilterMultiMother::processID

private

Definition at line 51 of file PythiaFilterMultiMother.h.

Referenced by filter(), and PythiaFilterMultiMother().

double PythiaFilterMultiMother::rapidity

private

Definition at line 47 of file PythiaFilterMultiMother.h.

Referenced by filter().

int PythiaFilterMultiMother::status

private

Definition at line 49 of file PythiaFilterMultiMother.h.

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

edm::EDGetTokenT<edm::HepMCProduct> PythiaFilterMultiMother::token_

private

Definition at line 34 of file PythiaFilterMultiMother.h.

Referenced by filter().