#include <L1Trigger/L1TNtuples/src/L1JetRecoTreeProducer.cc>

Inheritance diagram for L1JetRecoTreeProducer:

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

	~L1JetRecoTreeProducer ()

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

Public Attributes
L1Analysis::L1AnalysisRecoJetDataFormat *	jet_data

L1Analysis::L1AnalysisRecoMetDataFormat *	met_data

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

virtual void	beginJob (void)

void	doCaloJets (edm::Handle< reco::CaloJetCollection > caloJets)

void	doPFJetCorr (edm::Handle< reco::PFJetCollection > pfJets, edm::Handle< reco::JetCorrector > pfJetCorr)

void	doPFJets (edm::Handle< reco::PFJetCollection > pfJets)

void	doPFMet (edm::Handle< reco::PFMETCollection > pfMet)

virtual void	endJob ()

Private Attributes
edm::Service< TFileService >	fs_

edm::EDGetTokenT < reco::JetCorrector >	jecToken_

double	jetptThreshold_

unsigned int	maxCl_

unsigned int	maxJet_

unsigned int	maxTrk_

unsigned int	maxVtx_

bool	pfJetCorrMissing_

bool	pfJetsMissing_

edm::EDGetTokenT < reco::PFJetCollection >	pfJetToken_

bool	pfMetMissing_

edm::EDGetTokenT < reco::PFMETCollection >	pfMetToken_

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: Produces tree containing reco quantities

Definition at line 56 of file L1JetRecoTreeProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 109 of file L1JetRecoTreeProducer.cc.

References fs_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_FULL_cff::InputTag, jecToken_, jet_data, jetptThreshold_, TFileService::make(), maxJet_, met_data, pfJetToken_, pfMetToken_, and tree_.

                                                                           :
   pfJetsMissing_(false),
   pfMetMissing_(false),
   pfJetCorrMissing_(false)
 {
   
   //  caloJetToken_ = consumes<reco::CaloJetCollection>(iConfig.getUntrackedParameter("caloJetToken",edm::InputTag("ak4CaloJets")));
   pfJetToken_ = consumes<reco::PFJetCollection>(iConfig.getUntrackedParameter("pfJetToken",edm::InputTag("ak4PFJetsCHS")));
   //  caloJetIdToken_ = consumes<edm::ValueMap<reco::JetID> >(iConfig.getUntrackedParameter("jetIdToken",edm::InputTag("ak4JetID")));
   jecToken_ = consumes<reco::JetCorrector>(iConfig.getUntrackedParameter<edm::InputTag>("jecToken",edm::InputTag("ak4PFCHSL1FastL2L3ResidualCorrector")));
 
   pfMetToken_ = consumes<reco::PFMETCollection>(iConfig.getUntrackedParameter("pfMetToken",edm::InputTag("pfMet")));
 
   jetptThreshold_ = iConfig.getParameter<double>      ("jetptThreshold");
   maxJet_         = iConfig.getParameter<unsigned int>("maxJet");
 
   jet_data = new L1Analysis::L1AnalysisRecoJetDataFormat();
   met_data = new L1Analysis::L1AnalysisRecoMetDataFormat();
 
   // set up output
   tree_=fs_->make<TTree>("JetRecoTree", "JetRecoTree");
   tree_->Branch("Jet", "L1Analysis::L1AnalysisRecoJetDataFormat", &jet_data, 32000, 3);
   tree_->Branch("Sums", "L1Analysis::L1AnalysisRecoMetDataFormat", &met_data, 32000, 3);
 
 }

L1JetRecoTreeProducer::~L1JetRecoTreeProducer ( )

Definition at line 136 of file L1JetRecoTreeProducer.cc.

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

Member Function Documentation

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

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 150 of file L1JetRecoTreeProducer.cc.

References doPFJetCorr(), doPFJets(), doPFMet(), edm::Event::getByToken(), edm::HandleBase::isValid(), jecToken_, jet_data, met_data, L1Analysis::L1AnalysisRecoJetDataFormat::nJets, pfJetCorrMissing_, pfJets_cff::pfJets, pfJetsMissing_, pfJetToken_, PFMET_cfi::pfMet, pfMetMissing_, pfMetToken_, L1Analysis::L1AnalysisRecoJetDataFormat::Reset(), L1Analysis::L1AnalysisRecoMetDataFormat::Reset(), and tree_.

 {
 
   jet_data->Reset();
   met_data->Reset();
   
   // get jets
   edm::Handle<reco::PFJetCollection> pfJets;
   iEvent.getByToken(pfJetToken_, pfJets);
 
   //get sums
   edm::Handle<reco::PFMETCollection> pfMet;
   iEvent.getByToken(pfMetToken_, pfMet);
   
   // get jet ID
   //  edm::Handle<edm::ValueMap<reco::JetID> > jetsID;
   //iEvent.getByLabel(jetIdTag_,jetsID);
 
   edm::Handle<reco::JetCorrector> pfJetCorr;
   iEvent.getByToken(jecToken_, pfJetCorr);
 
 
   if (pfJets.isValid()) {
 
     jet_data->nJets=0;
 
     doPFJets(pfJets);
 
   }
   else {
     if (!pfJetsMissing_) {edm::LogWarning("MissingProduct") << "PFJets not found.  Branch will not be filled" << std::endl;}
     pfJetsMissing_ = true;
   }
  
   if (pfJetCorr.isValid()) {
  
     doPFJetCorr(pfJets,pfJetCorr);
  
   }
   else {
     if (!pfJetCorrMissing_)  {edm::LogWarning("MissingProduct") << "Jet Corrector not found.  Branch will not be filled" << std::endl;}
     pfJetCorrMissing_ = true;
   }
  
   if (pfMet.isValid()) {
  
     doPFMet(pfMet);
 
   }
   else {
     if (!pfMetMissing_) {edm::LogWarning("MissingProduct") << "PFMet not found.  Branch will not be filled" << std::endl;}
     pfMetMissing_ = true;
   }
 
 
   tree_->Fill();
 
 }

void L1JetRecoTreeProducer::beginJob ( void )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 346 of file L1JetRecoTreeProducer.cc.

347 {

348 }

void L1JetRecoTreeProducer::doCaloJets ( edm::Handle< reco::CaloJetCollection > caloJets )

private

Definition at line 211 of file L1JetRecoTreeProducer.cc.

                                                                            {
 
 
   for( auto it=caloJets->begin();
        it!=caloJets->end() && jet_data->nJets < maxJet_;
        ++it) {
     
     jet_data->et.push_back(it->et());
     jet_data->eta.push_back(it->eta());
     jet_data->phi.push_back(it->phi());
     jet_data->e.push_back(it->energy());
     jet_data->isPF.push_back(false);
     
     jet_data->eEMF.push_back(it->emEnergyFraction());
     jet_data->eEmEB.push_back(it->emEnergyInEB());
     jet_data->eEmEE.push_back(it->emEnergyInEE());
     jet_data->eEmHF.push_back(it->emEnergyInHF());
     jet_data->eHadHB.push_back(it->hadEnergyInHB());
     jet_data->eHadHE.push_back(it->hadEnergyInHE());
     jet_data->eHadHO.push_back(it->hadEnergyInHO());
     jet_data->eHadHF.push_back(it->hadEnergyInHF());
     jet_data->eMaxEcalTow.push_back(it->maxEInEmTowers());
     jet_data->eMaxHcalTow.push_back(it->maxEInHadTowers());
     jet_data->towerArea.push_back(it->towersArea());
 
     jet_data->nJets++;
 
   }
 
 
 }

void L1JetRecoTreeProducer::doPFJetCorr	(	edm::Handle< reco::PFJetCollection >	pfJets,
		edm::Handle< reco::JetCorrector >	pfJetCorr
	)

private

Definition at line 302 of file L1JetRecoTreeProducer.cc.

References corr, L1Analysis::L1AnalysisRecoJetDataFormat::corrFactor, L1Analysis::L1AnalysisRecoJetDataFormat::etCorr, jet_data, maxJet_, and edm::Handle< T >::product().

Referenced by analyze().

                                                                                                                {
 
 
   std::vector< std::pair<float,float> > corrJetEtsAndCorrs;
   
   //get jet correction and fill corrected jet ets and corrections
   for( auto it=pfJets->begin(); it!=pfJets->end(); ++it) 
     {
       float corr = pfJetCorr.product()->correction(*it);
       std::pair<float,float> corrJetEtAndCorr(corr*it->et(),corr);
       corrJetEtsAndCorrs.push_back(corrJetEtAndCorr);
     }
 
   // sort corrected jet ets and correction factors 
   // by corrected jet et
   std::sort(corrJetEtsAndCorrs.rbegin(),corrJetEtsAndCorrs.rend());
   
   //fill jet data array with sorted jet ets and corr factors
   std::vector<std::pair<float,float> >::iterator it;
   uint nJets = 0;
   
   for(it = corrJetEtsAndCorrs.begin(); it != corrJetEtsAndCorrs.end() && nJets < maxJet_; ++it){
     jet_data->etCorr.push_back(it->first);
     jet_data->corrFactor.push_back(it->second);
     nJets++;
   }
 }

void L1JetRecoTreeProducer::doPFJets ( edm::Handle< reco::PFJetCollection > pfJets )

private

Definition at line 245 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

                                                                      {
   
   float mHx = 0;
   float mHy = 0;
 
   met_data->Ht     = 0;
   met_data->mHt    = -999.;
   met_data->mHtPhi = -999.;
 
   for( auto it=pfJets->begin();
        it!=pfJets->end() && jet_data->nJets < maxJet_;
        ++it) {
     jet_data->et.push_back(it->et());
     jet_data->eta.push_back(it->eta());
     jet_data->phi.push_back(it->phi());
     jet_data->e.push_back(it->energy());
     jet_data->isPF.push_back(true);
 
     jet_data->chef.push_back(it->chargedHadronEnergyFraction());
     jet_data->nhef.push_back(it->neutralHadronEnergyFraction());
     jet_data->pef.push_back(it->photonEnergyFraction());
     jet_data->eef.push_back(it->electronEnergyFraction());
     jet_data->mef.push_back(it->muonEnergyFraction());
     jet_data->hfhef.push_back(it->HFHadronEnergyFraction());
     jet_data->hfemef.push_back(it->HFEMEnergyFraction());
     jet_data->chMult.push_back(it->chargedHadronMultiplicity());
     jet_data->nhMult.push_back(it->neutralHadronMultiplicity());
     jet_data->phMult.push_back(it->photonMultiplicity());
     jet_data->elMult.push_back(it->electronMultiplicity());
     jet_data->muMult.push_back(it->muonMultiplicity());
     jet_data->hfhMult.push_back(it->HFHadronMultiplicity());
     jet_data->hfemMult.push_back(it->HFEMMultiplicity());
 
     jet_data->cemef.push_back(it->chargedEmEnergyFraction());     
     jet_data->cmef.push_back(it->chargedMuEnergyFraction());      
     jet_data->nemef.push_back(it->neutralEmEnergyFraction());     
     jet_data->cMult.push_back(it->chargedMultiplicity());
     jet_data->nMult.push_back(it->neutralMultiplicity()); 
     
     jet_data->nJets++;
 
     
     if (it->pt()>jetptThreshold_){
       mHx += -1.*it->px();
       mHy += -1.*it->py();
       met_data->Ht  += it->pt();
     }
   }
 
   TVector2 *tv2 = new TVector2(mHx,mHy);
   met_data->mHt    = tv2->Mod();
   met_data->mHtPhi = tv2->Phi();
 
 }

void L1JetRecoTreeProducer::doPFMet ( edm::Handle< reco::PFMETCollection > pfMet )

private

Definition at line 331 of file L1JetRecoTreeProducer.cc.

References reco::LeafCandidate::et(), L1Analysis::L1AnalysisRecoMetDataFormat::met, met_data, L1Analysis::L1AnalysisRecoMetDataFormat::metPhi, reco::LeafCandidate::phi(), edm::Handle< T >::product(), L1Analysis::L1AnalysisRecoMetDataFormat::sumEt, and reco::MET::sumEt().

Referenced by analyze().

                                                                    {
 
   const reco::PFMETCollection *metCol = pfMet.product();
   const reco::PFMET theMet = metCol->front();
 
   met_data->met     = theMet.et();
   met_data->metPhi  = theMet.phi();
   met_data->sumEt   = theMet.sumEt();
 
 }