CMS 3D CMS Logo
List of all members | Public Member Functions | Static Public Member Functions | Protected Attributes
LHETablesProducer Class Reference
Inheritance diagram for LHETablesProducer:
edm::global::EDProducer<> edm::global::EDProducerBase edm::ProducerBase edm::EDConsumerBase edm::ProductRegistryHelper

Public Member Functions

void fillLHEObjectTable (const LHEEventProduct &lheProd, nanoaod::FlatTable &out) const
 
 LHETablesProducer (edm::ParameterSet const &params)
 
void produce (edm::StreamID id, edm::Event &iEvent, const edm::EventSetup &iSetup) const override
 
 ~LHETablesProducer () override
 
- Public Member Functions inherited from edm::global::EDProducer<>
 EDProducer ()=default
 
- Public Member Functions inherited from edm::global::EDProducerBase
 EDProducerBase ()
 
ModuleDescription const & moduleDescription () const
 
 ~EDProducerBase () 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)
 
virtual ~ProducerBase () noexcept(false)
 
- 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
 
void updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
 
virtual ~EDConsumerBase () noexcept(false)
 

Static Public Member Functions

static void fillDescriptions (edm::ConfigurationDescriptions &descriptions)
 
- Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string & baseType ()
 
static void fillDescriptions (ConfigurationDescriptions &descriptions)
 
static void prevalidate (ConfigurationDescriptions &descriptions)
 

Protected Attributes

const edm::EDGetTokenT< LHEEventProduct > lheTag_
 

Additional Inherited Members

- Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase 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
 
- 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 13 of file LHETablesProducer.cc.

Constructor & Destructor Documentation

LHETablesProducer::LHETablesProducer ( edm::ParameterSet const &  params)
inline

Definition at line 15 of file LHETablesProducer.cc.

15  :
16  lheTag_(consumes<LHEEventProduct>(params.getParameter<edm::InputTag>("lheInfo")))
17  {
18  produces<nanoaod::FlatTable>();
19  }
LHETablesProducer::lheTag_
const edm::EDGetTokenT< LHEEventProduct > lheTag_
Definition: LHETablesProducer.cc:95
LHETablesProducer::~LHETablesProducer ( )
inlineoverride

Definition at line 21 of file LHETablesProducer.cc.

21 {}

Member Function Documentation

static void LHETablesProducer::fillDescriptions ( edm::ConfigurationDescriptions &  descriptions)
inlinestatic

Definition at line 88 of file LHETablesProducer.cc.

References edm::ConfigurationDescriptions::add(), and edm::ParameterSetDescription::add().

88  {
89  edm::ParameterSetDescription desc;
90  desc.add<edm::InputTag>("lheInfo", edm::InputTag("externalLHEProducer"))->setComment("tag for the LHE information (LHEEventProduct)");
91  descriptions.add("lheInfoTable", desc);
92  }
edm::ParameterSetDescription::add
ParameterDescriptionBase * add(U const &iLabel, T const &value)
Definition: ParameterSetDescription.h:89
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:63
void LHETablesProducer::fillLHEObjectTable ( const LHEEventProduct &  lheProd,
nanoaod::FlatTable &  out 
) const
inline

Definition at line 34 of file LHETablesProducer.cc.

References funct::abs(), nanoaod::FlatTable::addColumnValue(), nanoaod::FlatTable::FloatColumn, LHEEventProduct::hepeup(), mps_fire::i, SiStripPI::max, gen::n, LHEEventProduct::npLO(), LHEEventProduct::npNLO(), EnergyCorrector::pt, lhef::HEPEUP::PUP, mps_update::status, nanoaod::FlatTable::UInt8Column, and findQualityFiles::v.

Referenced by produce().

34  {
35  double lheHT = 0, lheHTIncoming = 0;
36  unsigned int lheNj = 0, lheNb = 0, lheNc = 0, lheNuds = 0, lheNglu = 0;
37  double lheVpt = 0;
38 
39  const auto & hepeup = lheProd.hepeup();
40  const auto & pup = hepeup.PUP;
41  int lep = -1, lepBar = -1, nu = -1, nuBar = -1;
42  for (unsigned int i = 0, n = pup.size(); i < n; ++i) {
43  int status = hepeup.ISTUP[i];
44  int idabs = std::abs(hepeup.IDUP[i]);
45  if ( (status == 1) && ( ( idabs == 21 ) || (idabs > 0 && idabs < 7) ) ) { //# gluons and quarks
46  // object counters
47  lheNj++;
48  if (idabs==5) lheNb++;
49  else if (idabs==4) lheNc++;
50  else if (idabs <= 3 ) lheNuds++;
51  else if (idabs == 21) lheNglu++;
52  // HT
53  double pt = std::hypot( pup[i][0], pup[i][1] ); // first entry is px, second py
54  lheHT += pt;
55  int mothIdx = std::max(hepeup.MOTHUP[i].first-1,0); //first and last mother as pair; first entry has index 1 in LHE; incoming particles return motherindex 0
56  int mothIdxTwo = std::max(hepeup.MOTHUP[i].second-1,0);
57  int mothStatus = hepeup.ISTUP[mothIdx];
58  int mothStatusTwo = hepeup.ISTUP[mothIdxTwo];
59  bool hasIncomingAsMother = mothStatus<0 || mothStatusTwo<0;
60  if (hasIncomingAsMother) lheHTIncoming += pt;
61  } else if (idabs == 12 || idabs == 14 || idabs == 16) {
62  (hepeup.IDUP[i] > 0 ? nu : nuBar) = i;
63  } else if (idabs == 11 || idabs == 13 || idabs == 15) {
64  (hepeup.IDUP[i] > 0 ? lep : lepBar) = i;
65  }
66  }
67  std::pair<int,int> v(0,0);
68  if (lep != -1 && lepBar != -1) v = std::make_pair(lep,lepBar);
69  else if (lep != -1 && nuBar != -1) v = std::make_pair(lep, nuBar);
70  else if (nu != -1 && lepBar != -1) v = std::make_pair(nu ,lepBar);
71  else if (nu != -1 && nuBar != -1) v = std::make_pair(nu , nuBar);
72  if (v.first != -1 && v.second != -1) {
73  lheVpt = std::hypot( pup[v.first][0] + pup[v.second][0], pup[v.first][1] + pup[v.second][1] );
74  }
75 
76  out.addColumnValue<uint8_t>("Njets", lheNj, "Number of jets (partons) at LHE step", nanoaod::FlatTable::UInt8Column);
77  out.addColumnValue<uint8_t>("Nb", lheNb, "Number of b partons at LHE step", nanoaod::FlatTable::UInt8Column);
78  out.addColumnValue<uint8_t>("Nc", lheNc, "Number of c partons at LHE step", nanoaod::FlatTable::UInt8Column);
79  out.addColumnValue<uint8_t>("Nuds", lheNuds, "Number of u,d,s partons at LHE step", nanoaod::FlatTable::UInt8Column);
80  out.addColumnValue<uint8_t>("Nglu", lheNglu, "Number of gluon partons at LHE step", nanoaod::FlatTable::UInt8Column);
81  out.addColumnValue<float>("HT", lheHT, "HT, scalar sum of parton pTs at LHE step", nanoaod::FlatTable::FloatColumn);
82  out.addColumnValue<float>("HTIncoming", lheHTIncoming, "HT, scalar sum of parton pTs at LHE step, restricted to partons", nanoaod::FlatTable::FloatColumn);
83  out.addColumnValue<float>("Vpt", lheVpt, "pT of the W or Z boson at LHE step", nanoaod::FlatTable::FloatColumn);
84  out.addColumnValue<uint8_t>("NpNLO", lheProd.npNLO(), "number of partons at NLO", nanoaod::FlatTable::UInt8Column);
85  out.addColumnValue<uint8_t>("NpLO", lheProd.npLO(), "number of partons at LO", nanoaod::FlatTable::UInt8Column);
86  }
mps_fire.i
i
Definition: mps_fire.py:269
LHEEventProduct::hepeup
const lhef::HEPEUP & hepeup() const
Definition: LHEEventProduct.h:47
mps_update.status
status
Definition: mps_update.py:67
lhef::HEPEUP::PUP
std::vector< FiveVector > PUP
Definition: LesHouches.h:248
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
nanoaod::FlatTable::addColumnValue
void addColumnValue(const std::string &name, const C &value, const std::string &docString, ColumnType type=defaultColumnType< T >(), int mantissaBits=-1)
Definition: FlatTable.h:101
LHEEventProduct::npLO
int npLO() const
Definition: LHEEventProduct.h:41
LHEEventProduct::npNLO
int npNLO() const
Definition: LHEEventProduct.h:42
void LHETablesProducer::produce ( edm::StreamID  id,
edm::Event &  iEvent,
const edm::EventSetup &  iSetup 
) const
inlineoverride

Definition at line 23 of file LHETablesProducer.cc.

References fillLHEObjectTable(), edm::Event::getByToken(), nano_cff::lheInfo, lheTag_, eostools::move(), and edm::Event::put().

23  {
24  auto lheTab = std::make_unique<nanoaod::FlatTable>(1, "LHE", true);
25 
26  edm::Handle<LHEEventProduct> lheInfo;
27  if (iEvent.getByToken(lheTag_, lheInfo)) {
28  fillLHEObjectTable(*lheInfo, *lheTab);
29  }
30 
31  iEvent.put(std::move(lheTab));
32  }
edm::Event::put
OrphanHandle< PROD > put(std::unique_ptr< PROD > product)
Put a new product.
Definition: Event.h:127
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:508
LHETablesProducer::lheTag_
const edm::EDGetTokenT< LHEEventProduct > lheTag_
Definition: LHETablesProducer.cc:95
nano_cff.lheInfo
lheInfo
Definition: nano_cff.py:53
LHETablesProducer::fillLHEObjectTable
void fillLHEObjectTable(const LHEEventProduct &lheProd, nanoaod::FlatTable &out) const
Definition: LHETablesProducer.cc:34
eostools.move
def move(src, dest)
Definition: eostools.py:510

Member Data Documentation

const edm::EDGetTokenT<LHEEventProduct> LHETablesProducer::lheTag_
protected

Definition at line 95 of file LHETablesProducer.cc.

Referenced by produce().

Generated for CMSSW Reference Manual by doxygen 1.8.11