#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_

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

L1GenTreeProducer::~L1GenTreeProducer ( )

Definition at line 102 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 117 of file L1GenTreeProducer.cc.

 {
   
   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);
 
   for(size_t i = 0; i < genParticles->size(); ++ i) {
     const reco::GenParticle & p = (*genParticles)[i];
     int id = p.pdgId();
     //int st = p.status();  
     if (abs(id) == 13) {
       unsigned int nMo=p.numberOfMothers();
       //std::cout << "id " << id << "; st " << st 
       //<< "; nMo " << nMo << std::endl;
       for(unsigned int i=0;i<nMo;++i){
     //int thisParentID = dynamic_cast<const reco::GenParticle*>(p.mother(i))->pdgId();
     //std::cout << "   mother ID " << thisParentID << std::endl;
       }
     }
     
     //
     // 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);
       }
   }
 
 
   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_->nPUPoissonMean = puItr->getTrueNumInteractions();
       l1GenData_->nVtx           = puItr->getPU_NumInteractions();
       break;
     }
   }
 
 
 
   tree_->Fill();
 
 }