Inheritance diagram for LHETablesProducer:

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

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsStreamLuminosityBlocks () const final

bool	wantsStreamRuns () const final

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

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) 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.

                                                             :
             lheTag_(consumes<LHEEventProduct>(params.getParameter<edm::InputTag>("lheInfo")))
         {
             produces<nanoaod::FlatTable>();
         }

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().

                                                                                   {
             edm::ParameterSetDescription desc;
             desc.add<edm::InputTag>("lheInfo", edm::InputTag("externalLHEProducer"))->setComment("tag for the LHE information (LHEEventProduct)");
             descriptions.add("lheInfoTable", desc);
         }

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().

                                                                                                {
             double lheHT = 0, lheHTIncoming = 0;
             unsigned int lheNj = 0, lheNb = 0, lheNc = 0, lheNuds = 0, lheNglu = 0;
             double lheVpt = 0;
 
             const auto & hepeup = lheProd.hepeup();
             const auto & pup = hepeup.PUP;
             int lep = -1, lepBar = -1, nu = -1, nuBar = -1;
             for (unsigned int i = 0, n = pup.size(); i  < n; ++i) {
                 int status = hepeup.ISTUP[i];
                 int idabs = std::abs(hepeup.IDUP[i]);
                 if ( (status == 1) && ( ( idabs == 21 ) || (idabs > 0 && idabs < 7) ) ) { //# gluons and quarks
                     // object counters
                     lheNj++;
                     if (idabs==5) lheNb++;
                     else if (idabs==4) lheNc++;
                     else if (idabs <= 3 ) lheNuds++;
                     else if (idabs == 21) lheNglu++;
                     // HT
                     double pt = std::hypot( pup[i][0], pup[i][1] ); // first entry is px, second py
                     lheHT += pt;
                     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
                     int mothIdxTwo = std::max(hepeup.MOTHUP[i].second-1,0);
                     int mothStatus  = hepeup.ISTUP[mothIdx];
                     int mothStatusTwo  = hepeup.ISTUP[mothIdxTwo];
                     bool hasIncomingAsMother = mothStatus<0 || mothStatusTwo<0;
                     if (hasIncomingAsMother) lheHTIncoming += pt;
                 } else if (idabs == 12 || idabs == 14 || idabs == 16) {
                     (hepeup.IDUP[i] > 0 ? nu : nuBar) = i;
                 } else if (idabs == 11 || idabs == 13 || idabs == 15) {
                     (hepeup.IDUP[i] > 0 ? lep : lepBar) = i;
                 }
             }
             std::pair<int,int> v(0,0);
             if      (lep != -1 && lepBar != -1) v = std::make_pair(lep,lepBar);
             else if (lep != -1 &&  nuBar != -1) v = std::make_pair(lep, nuBar);
             else if (nu  != -1 && lepBar != -1) v = std::make_pair(nu ,lepBar);
             else if (nu  != -1 &&  nuBar != -1) v = std::make_pair(nu , nuBar);
             if (v.first != -1 && v.second != -1) {
                 lheVpt = std::hypot( pup[v.first][0] + pup[v.second][0], pup[v.first][1] + pup[v.second][1] ); 
             }
 
             out.addColumnValue<uint8_t>("Njets", lheNj, "Number of jets (partons) at LHE step", nanoaod::FlatTable::UInt8Column);
             out.addColumnValue<uint8_t>("Nb",    lheNb, "Number of b partons at LHE step", nanoaod::FlatTable::UInt8Column);
             out.addColumnValue<uint8_t>("Nc",    lheNc, "Number of c partons at LHE step", nanoaod::FlatTable::UInt8Column);
             out.addColumnValue<uint8_t>("Nuds", lheNuds, "Number of u,d,s partons at LHE step", nanoaod::FlatTable::UInt8Column);
             out.addColumnValue<uint8_t>("Nglu", lheNglu, "Number of gluon partons at LHE step", nanoaod::FlatTable::UInt8Column);
             out.addColumnValue<float>("HT", lheHT, "HT, scalar sum of parton pTs at LHE step", nanoaod::FlatTable::FloatColumn);
             out.addColumnValue<float>("HTIncoming", lheHTIncoming, "HT, scalar sum of parton pTs at LHE step, restricted to partons", nanoaod::FlatTable::FloatColumn);
             out.addColumnValue<float>("Vpt", lheVpt, "pT of the W or Z boson at LHE step", nanoaod::FlatTable::FloatColumn);
             out.addColumnValue<uint8_t>("NpNLO", lheProd.npNLO(), "number of partons at NLO", nanoaod::FlatTable::UInt8Column);
             out.addColumnValue<uint8_t>("NpLO", lheProd.npLO(), "number of partons at LO", nanoaod::FlatTable::UInt8Column);
         }

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().

                                                                                                    {
             auto lheTab  = std::make_unique<nanoaod::FlatTable>(1, "LHE", true);
 
             edm::Handle<LHEEventProduct> lheInfo;
             if (iEvent.getByToken(lheTag_, lheInfo)) {
                 fillLHEObjectTable(*lheInfo, *lheTab);
             }
 
             iEvent.put(std::move(lheTab));
         }

Member Data Documentation

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

protected

Definition at line 95 of file LHETablesProducer.cc.

Referenced by produce().

Public Member Functions

Static Public Member Functions

Protected Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation