#include <PythiaMomDauFilter.h>

Inheritance diagram for PythiaMomDauFilter:

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

	PythiaMomDauFilter (const edm::ParameterSet &)

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

bool	chargeconju

int	daughterID

std::vector< int >	dauIDs

std::vector< int >	desIDs

edm::EDGetTokenT< edm::HepMCProduct >	label_

double	maxetacut

double	maxptcut

double	minetacut

double	minptcut

double	mom_maxetacut

double	mom_maxptcut

double	mom_minetacut

double	mom_minptcut

int	ndaughters

int	ndescendants

int	particleID

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

Description: Filter events using MotherId and ChildrenIds infos

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

Definition at line 44 of file PythiaMomDauFilter.h.

Constructor & Destructor Documentation

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

explicit

Definition at line 13 of file PythiaMomDauFilter.cc.

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

                                                                      :
 label_(consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared"))),
 particleID(iConfig.getUntrackedParameter<int>("ParticleID", 0)),
 daughterID(iConfig.getUntrackedParameter<int>("DaughterID", 0)),
 chargeconju(iConfig.getUntrackedParameter<bool>("ChargeConjugation", true)),
 ndaughters(iConfig.getUntrackedParameter<int>("NumberDaughters", 0)),
 ndescendants(iConfig.getUntrackedParameter<int>("NumberDescendants", 0)),
 minptcut(iConfig.getUntrackedParameter<double>("MinPt", 0.)),
 maxptcut(iConfig.getUntrackedParameter<double>("MaxPt", 14000.)),
 minetacut(iConfig.getUntrackedParameter<double>("MinEta", -10.)),
 maxetacut(iConfig.getUntrackedParameter<double>("MaxEta", 10.)),
 mom_minptcut(iConfig.getUntrackedParameter<double>("MomMinPt", 0.)),
 mom_maxptcut(iConfig.getUntrackedParameter<double>("MomMaxPt", 14000.)),
 mom_minetacut(iConfig.getUntrackedParameter<double>("MomMinEta", -10.)),
 mom_maxetacut(iConfig.getUntrackedParameter<double>("MomMaxEta", 10.)),
 betaBoost(iConfig.getUntrackedParameter("BetaBoost",0.))
 {
    //now do what ever initialization is needed
    vector<int> defdauID;
    defdauID.push_back(0);
    vector<int> defdesID;
    defdesID.push_back(0);
    dauIDs = iConfig.getUntrackedParameter< vector<int> >("DaughterIDs",    defdauID);
    desIDs = iConfig.getUntrackedParameter< vector<int> >("DescendantsIDs", defdesID);
 }

PythiaMomDauFilter::~PythiaMomDauFilter ( )

override

Definition at line 40 of file PythiaMomDauFilter.cc.

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

Member Function Documentation

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

override

Definition at line 54 of file PythiaMomDauFilter.cc.

References cms::dd::accepted(), betaBoost, chargeconju, class-composition::children, daughterID, dauIDs, desIDs, edm::Event::getByToken(), edm::HepMCProduct::GetEvent(), mps_fire::i, label_, maxetacut, maxptcut, minetacut, minptcut, mom_maxetacut, mom_maxptcut, mom_minetacut, mom_minptcut, ndaughters, ndescendants, AlCaHLTBitMon_ParallelJobs::p, particleID, PYCOMP, and MCFilterZboostHelper::zboost().

 {
  using namespace edm;
  bool accepted = false;
  bool mom_accepted = false;
  Handle<HepMCProduct> evt;
  iEvent.getByToken(label_, evt);
 
  const HepMC::GenEvent * myGenEvent = evt->GetEvent();
  int ndauac = 0;
  int ndau = 0;
  int ndesac = 0;
  int ndes   = 0;   
    
  for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();  p != myGenEvent->particles_end(); ++p ) {
      
   if( (*p)->pdg_id() != particleID ) continue ;
   HepMC::FourVector mom = MCFilterZboostHelper::zboost((*p)->momentum(),betaBoost);
   if( (*p)->momentum().perp() >  mom_minptcut  && (*p)->momentum().perp() <  mom_maxptcut  && mom.eta()  >  mom_minetacut &&  mom.eta()  <  mom_maxetacut ){ 
    mom_accepted = true;
    ndauac = 0;
    ndau   = 0;
    ndesac = 0;
    ndes   = 0;     
    if ( (*p)->end_vertex() ) {  
     for ( HepMC::GenVertex::particle_iterator dau  =(*p)->end_vertex()->particles_begin(HepMC::children); dau != (*p)->end_vertex()->particles_end(HepMC::children); ++dau ) {
       ++ndau;
       for( unsigned int i=0; i<dauIDs.size(); ++i) {
        if( (*dau)->pdg_id() != dauIDs[i] ) continue ;
         ++ndauac;
         break;
      }
       if((*dau)->pdg_id() == daughterID){
        for(HepMC::GenVertex::particle_iterator des  = (*dau)->end_vertex()->particles_begin(HepMC::children); des != (*des)->end_vertex()->particles_end(HepMC::children); ++des ){
         ++ndes;        
         for( unsigned int i=0; i<desIDs.size(); ++i) {
          if( (*des)->pdg_id() != desIDs[i] ) continue ;
          HepMC::FourVector dau_i = MCFilterZboostHelper::zboost((*des)->momentum(),betaBoost);
          if( (*des)->momentum().perp() >  minptcut  && (*des)->momentum().perp() <  maxptcut  && dau_i.eta()  >  minetacut &&  dau_i.eta()  <  maxetacut ) {
           ++ndesac;
           break;
       } 
      }
     }                    
        }
       } 
      }
     }  
     if( ndau ==  ndaughters && ndauac == ndaughters && mom_accepted && ndes == ndescendants && ndesac == ndescendants ) {
        accepted = true;
        break;
      }    
      
    }
    
    
    if( !accepted && chargeconju ) {
    
      mom_accepted = false;
      for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin(); p != myGenEvent->particles_end(); ++p ) {
        
        if( (*p)->pdg_id() != -particleID ) continue ;
        HepMC::FourVector mom = MCFilterZboostHelper::zboost((*p)->momentum(),betaBoost);
        if( (*p)->momentum().perp() >  mom_minptcut  && (*p)->momentum().perp() <  mom_maxptcut  &&  mom.eta()  >  mom_minetacut && mom.eta()  <  mom_maxetacut ){ 
        mom_accepted = true;
        ndauac = 0;
        ndau = 0;
        ndesac = 0;
        ndes   = 0;      
        if ( (*p)->end_vertex() ) {
      for ( HepMC::GenVertex::particle_iterator dau  =(*p)->end_vertex()->particles_begin(HepMC::children); dau != (*p)->end_vertex()->particles_end(HepMC::children);  ++dau ) {
        ++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( (*dau)->pdg_id() != IDanti ) continue ;
             ++ndauac;
             break;
            }
        int daughterIDanti = -daughterID;
        int pythiaCode = PYCOMP(daughterID);
        int has_antipart = pydat2.kchg[3-1][pythiaCode-1];
        if( has_antipart == 0 ) daughterIDanti = daughterID;
        if((*dau)->pdg_id() == daughterIDanti) {
             for( HepMC::GenVertex::particle_iterator des  = (*dau)->end_vertex()->particles_begin(HepMC::children); des != (*des)->end_vertex()->particles_end(HepMC::children);  ++des ){
                  ++ndes;
         for( unsigned int i=0; i<desIDs.size(); ++i) {
          int IDanti = -desIDs[i];
          int pythiaCode = PYCOMP(desIDs[i]);
          int has_antipart = pydat2.kchg[3-1][pythiaCode-1];
          if( has_antipart == 0 ) IDanti = desIDs[i];
          if( (*des)->pdg_id() != IDanti ) continue ;
              HepMC::FourVector dau_i = MCFilterZboostHelper::zboost((*des)->momentum(),betaBoost);
          if(   (*des)->momentum().perp() >  minptcut  && (*des)->momentum().perp() <  maxptcut  &&  dau_i.eta()  >  minetacut && dau_i.eta()  <  maxetacut ) {
            ++ndesac;
            break;
          } 
        }                     
           }  
      }
     } 
        }
       }
        if( ndau ==  ndaughters && ndauac == ndaughters && mom_accepted && ndes == ndescendants && ndesac == ndescendants ) {
      accepted = true;
      break;
        }    
      }
      
    }    
 
    if (accepted){
    return true; } else {return false;}
 
 }

Member Data Documentation

double PythiaMomDauFilter::betaBoost

private

Definition at line 70 of file PythiaMomDauFilter.h.

Referenced by filter().

bool PythiaMomDauFilter::chargeconju

private

Definition at line 59 of file PythiaMomDauFilter.h.

Referenced by filter().

int PythiaMomDauFilter::daughterID

private

Definition at line 58 of file PythiaMomDauFilter.h.

Referenced by filter().

std::vector<int> PythiaMomDauFilter::dauIDs

private

Definition at line 55 of file PythiaMomDauFilter.h.

Referenced by filter(), and PythiaMomDauFilter().

std::vector<int> PythiaMomDauFilter::desIDs

private

Definition at line 56 of file PythiaMomDauFilter.h.

Referenced by filter(), and PythiaMomDauFilter().

edm::EDGetTokenT<edm::HepMCProduct> PythiaMomDauFilter::label_

private

Definition at line 54 of file PythiaMomDauFilter.h.

Referenced by SequenceTypes.DummyModule::__repr__(), Mixins._Labelable::_findDependencies(), filter(), and Mixins._Labelable::setLabel().

double PythiaMomDauFilter::maxetacut

private

Definition at line 65 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::maxptcut

private

Definition at line 63 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::minetacut

private

Definition at line 64 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::minptcut

private

Definition at line 62 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::mom_maxetacut

private

Definition at line 69 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::mom_maxptcut

private

Definition at line 67 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::mom_minetacut

private

Definition at line 68 of file PythiaMomDauFilter.h.

Referenced by filter().

double PythiaMomDauFilter::mom_minptcut

private

Definition at line 66 of file PythiaMomDauFilter.h.

Referenced by filter().

int PythiaMomDauFilter::ndaughters

private

Definition at line 60 of file PythiaMomDauFilter.h.

Referenced by filter().

int PythiaMomDauFilter::ndescendants

private

Definition at line 61 of file PythiaMomDauFilter.h.

Referenced by filter().

int PythiaMomDauFilter::particleID

private

Definition at line 57 of file PythiaMomDauFilter.h.

Referenced by filter().

Public Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation