#include <L1TriggerDPG/L1Ntuples/src/L1GenTreeProducer.cc>

Inheritance diagram for L1GenTreeProducer:

Public Member Functions
	L1GenTreeProducer (const edm::ParameterSet &)

	~L1GenTreeProducer () override

Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

	EDAnalyzer ()

SerialTaskQueue *	globalLuminosityBlocksQueue ()

SerialTaskQueue *	globalRunsQueue ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

	~EDAnalyzer () 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 Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	beginJob (void) override

void	endJob () override

Private Attributes
edm::Service< TFileService >	fs_

edm::EDGetTokenT< GenEventInfoProduct >	genInfoToken_

edm::EDGetTokenT< reco::GenJetCollection >	genJetToken_

edm::EDGetTokenT< reco::GenParticleCollection >	genParticleToken_

std::unique_ptr< L1Analysis::L1AnalysisGeneratorDataFormat >	l1GenData_

unsigned	maxL1Upgrade_

edm::EDGetTokenT< std::vector< PileupSummaryInfo > >	pileupInfoToken_

TTree *	tree_

Additional Inherited Members
Public Types inherited from edm::EDAnalyzer
typedef EDAnalyzer	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::EDAnalyzer
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: Produce L1 Extra tree

Implementation:

Definition at line 54 of file L1GenTreeProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 83 of file L1GenTreeProducer.cc.

References fs_, genInfoToken_, genJetToken_, genParticleToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), l1GenData_, TFileService::make(), pileupInfoToken_, and tree_.

                                                                    {
   genJetToken_ = consumes<reco::GenJetCollection>(iConfig.getUntrackedParameter<edm::InputTag>("genJetToken"));
   genParticleToken_ =
       consumes<reco::GenParticleCollection>(iConfig.getUntrackedParameter<edm::InputTag>("genParticleToken"));
   pileupInfoToken_ =
       consumes<std::vector<PileupSummaryInfo>>(iConfig.getUntrackedParameter<edm::InputTag>("pileupInfoToken"));
   genInfoToken_ = consumes<GenEventInfoProduct>(iConfig.getParameter<edm::InputTag>("genInfoToken"));
 
   l1GenData_ = std::make_unique<L1Analysis::L1AnalysisGeneratorDataFormat>();
 
   // set up output
   tree_ = fs_->make<TTree>("L1GenTree", "L1GenTree");
   tree_->Branch("Generator", "L1Analysis::L1AnalysisGeneratorDataFormat", l1GenData_.get(), 32000, 3);
 }

L1GenTreeProducer::~L1GenTreeProducer ( )

override

Definition at line 98 of file L1GenTreeProducer.cc.

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

Member Function Documentation

void L1GenTreeProducer::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overrideprivate

Definition at line 108 of file L1GenTreeProducer.cc.

References funct::abs(), GenEventInfoProduct::binningValues(), l1GtPatternGenerator_cfi::bx, reco::LeafCandidate::charge(), reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), Exception, genInfoToken_, ttbarCategorization_cff::genJets, genJetToken_, genParticles2HepMC_cfi::genParticles, genParticleToken_, edm::Event::getByToken(), GenEventInfoProduct::hasBinningValues(), mps_fire::i, edm::HandleBase::isValid(), l1GenData_, reco::CompositeRefCandidateT< D >::mother(), nPart(), reco::CompositeRefCandidateT< D >::numberOfMothers(), AlCaHLTBitMon_ParallelJobs::p, EgammaValidation_cff::pdgId, reco::LeafCandidate::pdgId(), reco::LeafCandidate::phi(), pileupInfoToken_, reco::LeafCandidate::pt(), reco::LeafCandidate::status(), tree_, and GenEventInfoProduct::weight().

                                                                                    {
   edm::Handle<GenEventInfoProduct> genInfo;
   iEvent.getByToken(genInfoToken_, genInfo);
 
   l1GenData_->Reset();
 
   if (genInfo.isValid()) {
     l1GenData_->weight = genInfo->weight();
     l1GenData_->pthat = genInfo->hasBinningValues() ? (genInfo->binningValues())[0] : 0.0;
   }
 
   edm::Handle<reco::GenJetCollection> genJets;
   iEvent.getByToken(genJetToken_, genJets);
 
   if (genJets.isValid()) {
     reco::GenJetCollection::const_iterator jetItr = genJets->begin();
     reco::GenJetCollection::const_iterator jetEnd = genJets->end();
     for (; jetItr != jetEnd; ++jetItr) {
       l1GenData_->jetPt.push_back(jetItr->pt());
       l1GenData_->jetEta.push_back(jetItr->eta());
       l1GenData_->jetPhi.push_back(jetItr->phi());
       l1GenData_->nJet++;
     }
 
   } else {
     edm::LogWarning("MissingProduct") << "Gen jets not found. Branch will not be filled" << std::endl;
   }
 
   edm::Handle<reco::GenParticleCollection> genParticles;
   iEvent.getByToken(genParticleToken_, genParticles);
 
   if (genParticles.isValid()) {
     int nPart{0};
 
     for (size_t i = 0; i < genParticles->size(); ++i) {
       const reco::GenParticle& p = (*genParticles)[i];
       int id = p.pdgId();
 
       // See if the parent was interesting
       int parentID = -10000;
       unsigned int nMo = p.numberOfMothers();
       for (unsigned int i = 0; i < nMo; ++i) {
         int thisParentID = dynamic_cast<const reco::GenParticle*>(p.mother(i))->pdgId();
         //
         // Is this a bottom hadron?
         int hundredsIndex = abs(thisParentID) / 100;
         int thousandsIndex = abs(thisParentID) / 1000;
         if (((abs(thisParentID) >= 23) && (abs(thisParentID) <= 25)) || (abs(thisParentID) == 6) ||
             (hundredsIndex == 5) || (hundredsIndex == 4) || (thousandsIndex == 5) || (thousandsIndex == 4))
           parentID = thisParentID;
       }
       if ((parentID == -10000) && (nMo > 0))
         parentID = dynamic_cast<const reco::GenParticle*>(p.mother(0))->pdgId();
       //
       // If the parent of this particle is interesting, store all of the info
       if ((parentID != p.pdgId()) &&
           ((parentID > -9999) || (abs(id) == 11) || (abs(id) == 13) || (abs(id) == 23) || (abs(id) == 24) ||
            (abs(id) == 25) || (abs(id) == 4) || (abs(id) == 5) || (abs(id) == 6))) {
         l1GenData_->partId.push_back(p.pdgId());
         l1GenData_->partStat.push_back(p.status());
         l1GenData_->partPt.push_back(p.pt());
         l1GenData_->partEta.push_back(p.eta());
         l1GenData_->partPhi.push_back(p.phi());
         l1GenData_->partE.push_back(p.energy());
         l1GenData_->partParent.push_back(parentID);
         l1GenData_->partCh.push_back(p.charge());
         ++nPart;
       }
     }
     l1GenData_->nPart = nPart;
   }
 
   edm::Handle<std::vector<PileupSummaryInfo>> puInfoCollection;
   iEvent.getByToken(pileupInfoToken_, puInfoCollection);
 
   if (!puInfoCollection.isValid()) {
     throw cms::Exception("ProductNotValid") << "pileupInfoSource not valid";
   }
 
   // Loop over vector, find in-time entry, then store the relevant info
   std::vector<PileupSummaryInfo>::const_iterator puItr = puInfoCollection->begin();
   std::vector<PileupSummaryInfo>::const_iterator puEnd = puInfoCollection->end();
   for (; puItr != puEnd; ++puItr) {
     int bx = puItr->getBunchCrossing();
     if (bx == 0) {
       l1GenData_->nMeanPU = puItr->getTrueNumInteractions();
       l1GenData_->nVtx = puItr->getPU_NumInteractions();
       break;
     }
   }
 
   tree_->Fill();
 }