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

Inheritance diagram for L1GenTreeProducer:

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

	~L1GenTreeProducer ()

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

	EDAnalyzer ()

ModuleDescription const &	moduleDescription () const

std::string	workerType () const

virtual	~EDAnalyzer ()

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

ProductHolderIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductHolderIndexAndSkipBit > &) const

std::vector < ProductHolderIndexAndSkipBit > const &	itemsToGetFromEvent () const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesDependentUpon (std::string const &iProcessName, std::string const &iModuleLabel, bool iPrint, std::vector< char const * > &oModuleLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

bool	registeredToConsume (ProductHolderIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductHolderIndexHelper const &)

virtual	~EDConsumerBase ()

Private Member Functions
virtual void	analyze (const edm::Event &, const edm::EventSetup &)

virtual void	beginJob (void)

virtual void	endJob ()

Private Attributes
edm::Service< TFileService >	fs_

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 &)

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

Description: Produce L1 Extra tree

Implementation:

Definition at line 55 of file L1GenTreeProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 88 of file L1GenTreeProducer.cc.

References fs_, genJetToken_, genParticleToken_, 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"));
   
   l1GenData_ = std::make_unique<L1Analysis::L1AnalysisGeneratorDataFormat>();
 
   // set up output
   tree_=fs_->make<TTree>("L1GenTree", "L1GenTree");
   tree_->Branch("Generator", "L1Analysis::L1AnalysisGeneratorDataFormat", &l1GenData_, 32000, 3);
 }

L1GenTreeProducer::~L1GenTreeProducer ( )

Definition at line 103 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
	)

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 118 of file L1GenTreeProducer.cc.

References funct::abs(), reco::LeafCandidate::charge(), reco::LeafCandidate::energy(), reco::LeafCandidate::eta(), Exception, genJetToken_, genParticleCandidates2GenParticles_cfi::genParticles, genParticleToken_, edm::Event::getByToken(), i, edm::HandleBase::isValid(), l1GenData_, reco::CompositeRefCandidateT< D >::mother(), nPart(), reco::CompositeRefCandidateT< D >::numberOfMothers(), AlCaHLTBitMon_ParallelJobs::p, SingleMuPt40Fwdv2_cfi_GEN_SIM::pdgId, reco::LeafCandidate::pdgId(), reco::LeafCandidate::phi(), pileupInfoToken_, reco::LeafCandidate::pt(), reco::LeafCandidate::status(), and tree_.

 {
   
   l1GenData_->Reset();
 
   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();
 
 }