CMS 3D CMS Logo

List of all members | Public Member Functions | Private Attributes
PythiaFilter Class Reference

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

Inheritance diagram for PythiaFilter:
edm::global::EDFilter<> edm::global::EDFilterBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

bool filter (edm::StreamID, edm::Event &, const edm::EventSetup &) const override
 
 PythiaFilter (const edm::ParameterSet &)
 
 ~PythiaFilter () override
 
- Public Member Functions inherited from edm::global::EDFilter<>
 EDFilter ()=default
 
bool hasAbilityToProduceInLumis () const final
 
bool hasAbilityToProduceInRuns () const final
 
bool wantsGlobalLuminosityBlocks () const final
 
bool wantsGlobalRuns () const final
 
bool wantsStreamLuminosityBlocks () const final
 
bool wantsStreamRuns () const final
 
- Public Member Functions inherited from edm::global::EDFilterBase
 EDFilterBase ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDFilterBase () 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< ConsumesInfoconsumesInfo () 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
 
EDConsumerBaseoperator= (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

const double betaBoost
 
const double maxetacut
 
const double maxpcut
 
const double maxphicut
 
const double maxptcut
 
const double maxrapcut
 
const double minetacut
 
const double minpcut
 
const double minphicut
 
const double minptcut
 
const double minrapcut
 
const int motherID
 
const int particleID
 
const int processID
 
const int status
 
const edm::EDGetTokenT< edm::HepMCProducttoken_
 

Additional Inherited Members

- Public Types inherited from edm::global::EDFilterBase
typedef EDFilterBase 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::global::EDFilterBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 
- 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: <one line="" class="" summary>="">

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

Definition at line 41 of file PythiaFilter.h.

Constructor & Destructor Documentation

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

Definition at line 12 of file PythiaFilter.cc.

12  :
13 token_(consumes<edm::HepMCProduct>(edm::InputTag(iConfig.getUntrackedParameter("moduleLabel",std::string("generator")),"unsmeared"))),
14 particleID(iConfig.getUntrackedParameter("ParticleID", 0)),
15 minpcut(iConfig.getUntrackedParameter("MinP", 0.)),
16 maxpcut(iConfig.getUntrackedParameter("MaxP", 10000.)),
17 minptcut(iConfig.getUntrackedParameter("MinPt", 0.)),
18 maxptcut(iConfig.getUntrackedParameter("MaxPt", 10000.)),
19 minetacut(iConfig.getUntrackedParameter("MinEta", -10.)),
20 maxetacut(iConfig.getUntrackedParameter("MaxEta", 10.)),
21 minrapcut(iConfig.getUntrackedParameter("MinRapidity", -20.)),
22 maxrapcut(iConfig.getUntrackedParameter("MaxRapidity", 20.)),
23 minphicut(iConfig.getUntrackedParameter("MinPhi", -3.5)),
24 maxphicut(iConfig.getUntrackedParameter("MaxPhi", 3.5)),
25 status(iConfig.getUntrackedParameter("Status", 0)),
26 motherID(iConfig.getUntrackedParameter("MotherID", 0)),
27 processID(iConfig.getUntrackedParameter("ProcessID", 0)),
28 betaBoost(iConfig.getUntrackedParameter("BetaBoost",0.))
29 {
30  //now do what ever initialization is needed
31 
32 }
T getUntrackedParameter(std::string const &, T const &) const
const double maxetacut
Definition: PythiaFilter.h:58
const double minptcut
Definition: PythiaFilter.h:55
const int status
Definition: PythiaFilter.h:64
const double maxpcut
Definition: PythiaFilter.h:54
const double maxrapcut
Definition: PythiaFilter.h:60
const double minrapcut
Definition: PythiaFilter.h:59
const double maxptcut
Definition: PythiaFilter.h:56
const double minpcut
Definition: PythiaFilter.h:53
const int motherID
Definition: PythiaFilter.h:65
const int processID
Definition: PythiaFilter.h:66
const double maxphicut
Definition: PythiaFilter.h:62
const double minphicut
Definition: PythiaFilter.h:61
const double betaBoost
Definition: PythiaFilter.h:68
const double minetacut
Definition: PythiaFilter.h:57
const edm::EDGetTokenT< edm::HepMCProduct > token_
Definition: PythiaFilter.h:51
const int particleID
Definition: PythiaFilter.h:52
PythiaFilter::~PythiaFilter ( )
override

Definition at line 35 of file PythiaFilter.cc.

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

Member Function Documentation

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

Definition at line 49 of file PythiaFilter.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, motherID, AlCaHLTBitMon_ParallelJobs::p, particleID, processID, status, token_, and MCFilterZboostHelper::zboost().

50 {
51  using namespace edm;
52  bool accepted = false;
54  iEvent.getByToken(token_, evt);
55 
56  const HepMC::GenEvent * myGenEvent = evt->GetEvent();
57 
58  if(processID == 0 || processID == myGenEvent->signal_process_id()) {
59 
60  for ( HepMC::GenEvent::particle_const_iterator p = myGenEvent->particles_begin();
61  p != myGenEvent->particles_end(); ++p ) {
62  HepMC::FourVector mom = MCFilterZboostHelper::zboost((*p)->momentum(),betaBoost);
63  double rapidity = 0.5*log( (mom.e()+mom.pz()) / (mom.e()-mom.pz()) );
64 
65  if ( abs((*p)->pdg_id()) == particleID
66  && mom.rho() > minpcut
67  && mom.rho() < maxpcut
68  && (*p)->momentum().perp() > minptcut
69  && (*p)->momentum().perp() < maxptcut
70  && mom.eta() > minetacut
71  && mom.eta() < maxetacut
72  && rapidity > minrapcut
73  && rapidity < maxrapcut
74  && (*p)->momentum().phi() > minphicut
75  && (*p)->momentum().phi() < maxphicut ) {
76 
77 
78 
79  if (status == 0 && motherID == 0){
80  accepted = true;
81  }
82  if (status != 0 && motherID == 0){
83  if ((*p)->status() == status)
84  accepted = true;
85  }
86 
87  HepMC::GenParticle* mother = (*((*p)->production_vertex()->particles_in_const_begin()));
88 
89  if (status == 0 && motherID != 0){
90  if (abs(mother->pdg_id()) == abs(motherID)) {
91  accepted = true;
92  }
93  }
94  if (status != 0 && motherID != 0){
95 
96  if ((*p)->status() == status && abs(mother->pdg_id()) == abs(motherID)){
97  accepted = true;
98 
99  }
100  }
101 
102  /*
103  if (status == 0 && motherID != 0){
104  if (abs(((*p)->mother())->pdg_id()) == abs(motherID)) {
105  accepted = true;
106  }
107  }
108  if (status != 0 && motherID != 0){
109 
110  if ((*p)->status() == status && abs(((*p)->mother())->pdg_id()) == abs(motherID)){
111  accepted = true;
112 
113  }
114  }
115  */
116 
117  }
118  }
119 
120  } else { accepted = true; }
121 
122  if (accepted){
123  return true; } else {return false;}
124 
125 }
HepMC::FourVector zboost(const HepMC::FourVector &, double)
const double maxetacut
Definition: PythiaFilter.h:58
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:517
const double minptcut
Definition: PythiaFilter.h:55
const int status
Definition: PythiaFilter.h:64
const double maxpcut
Definition: PythiaFilter.h:54
const double maxrapcut
Definition: PythiaFilter.h:60
const double minrapcut
Definition: PythiaFilter.h:59
const double maxptcut
Definition: PythiaFilter.h:56
const double minpcut
Definition: PythiaFilter.h:53
const int motherID
Definition: PythiaFilter.h:65
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
const int processID
Definition: PythiaFilter.h:66
const double maxphicut
Definition: PythiaFilter.h:62
const HepMC::GenEvent * GetEvent() const
Definition: HepMCProduct.h:38
bool accepted(std::vector< std::string_view > const &, std::string_view)
const double minphicut
Definition: PythiaFilter.h:61
HLT enums.
const double betaBoost
Definition: PythiaFilter.h:68
const double minetacut
Definition: PythiaFilter.h:57
const edm::EDGetTokenT< edm::HepMCProduct > token_
Definition: PythiaFilter.h:51
const int particleID
Definition: PythiaFilter.h:52

Member Data Documentation

const double PythiaFilter::betaBoost
private

Definition at line 68 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::maxetacut
private

Definition at line 58 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::maxpcut
private

Definition at line 54 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::maxphicut
private

Definition at line 62 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::maxptcut
private

Definition at line 56 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::maxrapcut
private

Definition at line 60 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::minetacut
private

Definition at line 57 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::minpcut
private

Definition at line 53 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::minphicut
private

Definition at line 61 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::minptcut
private

Definition at line 55 of file PythiaFilter.h.

Referenced by filter().

const double PythiaFilter::minrapcut
private

Definition at line 59 of file PythiaFilter.h.

Referenced by filter().

const int PythiaFilter::motherID
private

Definition at line 65 of file PythiaFilter.h.

Referenced by filter().

const int PythiaFilter::particleID
private

Definition at line 52 of file PythiaFilter.h.

Referenced by filter().

const int PythiaFilter::processID
private

Definition at line 66 of file PythiaFilter.h.

Referenced by filter().

const int PythiaFilter::status
private
const edm::EDGetTokenT<edm::HepMCProduct> PythiaFilter::token_
private

Definition at line 51 of file PythiaFilter.h.

Referenced by filter().