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

Inheritance diagram for L1JetRecoTreeProducer:

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

	~L1JetRecoTreeProducer () override

Public Member Functions inherited from edm::one::EDAnalyzer< edm::one::SharedResources >
	EDAnalyzer ()=default

	EDAnalyzer (const EDAnalyzer &)=delete

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

const EDAnalyzer &	operator= (const EDAnalyzer &)=delete

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

bool	wantsInputProcessBlocks () const final

bool	wantsProcessBlocks () const final

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

	EDAnalyzerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

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

std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const

std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (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::array< std::vector< ModuleDescription const > , NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, 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

void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)

ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)

virtual	~EDConsumerBase () noexcept(false)

Public Attributes
L1Analysis::L1AnalysisRecoJetDataFormat *	jet_data

L1Analysis::L1AnalysisRecoMetDataFormat *	met_data

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

void	beginJob (void) override

bool	caloJetID (const reco::CaloJet &jet)

void	doCaloJetCorr (edm::Handle< reco::CaloJetCollection > caloJets, edm::Handle< reco::JetCorrector > caloJetCorr)

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

void	doCaloMet (edm::Handle< reco::CaloMETCollection > caloMet)

void	doCaloMetBE (edm::Handle< reco::CaloMETCollection > caloMetBE)

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)

void	doPFMetNoMu (edm::Handle< reco::PFMETCollection > pfMet, edm::Handle< reco::MuonCollection >)

void	endJob () override

bool	pfJetID (const reco::PFJet &jet)

Private Attributes
edm::EDGetTokenT< reco::JetCorrector >	caloJECToken_

bool	caloJetCorrMissing_

bool	caloJetIDMissing_

edm::EDGetTokenT< edm::ValueMap< reco::JetID > >	caloJetIDToken_

bool	caloJetsMissing_

edm::EDGetTokenT< reco::CaloJetCollection >	caloJetToken_

bool	caloMetBEMissing_

edm::EDGetTokenT< reco::CaloMETCollection >	caloMetBEToken_

bool	caloMetMissing_

edm::EDGetTokenT< reco::CaloMETCollection >	caloMetToken_

double	jetetaMax_

double	jetptThreshold_

unsigned int	maxCl_

unsigned int	maxJet_

unsigned int	maxTrk_

unsigned int	maxVtx_

bool	muonsMissing_

edm::EDGetTokenT< reco::MuonCollection >	muonToken_

edm::EDGetTokenT< reco::JetCorrector >	pfJECToken_

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::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

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)

void	resetItemsToGetFrom (BranchType iType)

Detailed Description

Description: Produces tree containing reco quantities

Definition at line 58 of file L1JetRecoTreeProducer.cc.

Constructor & Destructor Documentation

◆ L1JetRecoTreeProducer()

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

explicit

Definition at line 122 of file L1JetRecoTreeProducer.cc.

References caloJECToken_, caloJetIDToken_, caloJetToken_, caloMetBEToken_, caloMetToken_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_2022v15_cff::InputTag, jet_data, jetetaMax_, jetptThreshold_, TFileService::kSharedResource, TFileService::make(), maxJet_, met_data, muonToken_, pfJECToken_, pfJetToken_, pfMetToken_, and tree_.

     : pfJetsMissing_(false),
       pfJetCorrMissing_(false),
       caloJetCorrMissing_(false),
       caloJetsMissing_(false),
       caloJetIDMissing_(false),
       pfMetMissing_(false),
       caloMetMissing_(false),
       caloMetBEMissing_(false),
       muonsMissing_(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("caloJetIDToken", edm::InputTag("ak4JetID")));
   pfJECToken_ = consumes<reco::JetCorrector>(
       iConfig.getUntrackedParameter<edm::InputTag>("pfJECToken", edm::InputTag("ak4PFCHSL1FastL2L3ResidualCorrector")));
   caloJECToken_ = consumes<reco::JetCorrector>(iConfig.getUntrackedParameter<edm::InputTag>(
       "caloJECToken", edm::InputTag("ak4CaloL1FastL2L3ResidualCorrector")));
 
   pfMetToken_ = consumes<reco::PFMETCollection>(iConfig.getUntrackedParameter("pfMetToken", edm::InputTag("pfMetT1")));
   caloMetToken_ =
       consumes<reco::CaloMETCollection>(iConfig.getUntrackedParameter("caloMetToken", edm::InputTag("caloMet")));
   caloMetBEToken_ =
       consumes<reco::CaloMETCollection>(iConfig.getUntrackedParameter("caloMetBEToken", edm::InputTag("caloMetBE")));
 
   muonToken_ = consumes<reco::MuonCollection>(iConfig.getUntrackedParameter("muonToken", edm::InputTag("muons")));
 
   usesResource(TFileService::kSharedResource);
 
   jetptThreshold_ = iConfig.getParameter<double>("jetptThreshold");
   jetetaMax_ = iConfig.getParameter<double>("jetetaMax");
   maxJet_ = iConfig.getParameter<unsigned int>("maxJet");
 
   jet_data = new L1Analysis::L1AnalysisRecoJetDataFormat();
   met_data = new L1Analysis::L1AnalysisRecoMetDataFormat();
 
   // set up output
   edm::Service<TFileService> fs_;
   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::~L1JetRecoTreeProducer ( )

override

Definition at line 167 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

◆ analyze()

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

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 177 of file L1JetRecoTreeProducer.cc.

References caloJECToken_, caloJetCorrMissing_, caloJetIDMissing_, caloJetIDToken_, isolatedTracks_cfi::caloJets, caloJetsMissing_, caloJetToken_, heavyionUCCDQM_cfi::caloMet, CaloMET_cfi::caloMetBE, caloMetBEMissing_, caloMetBEToken_, caloMetMissing_, caloMetToken_, doCaloJetCorr(), doCaloJets(), doCaloMet(), doCaloMetBE(), doPFJetCorr(), doPFJets(), doPFMet(), doPFMetNoMu(), iEvent, edm::HandleBase::isValid(), jet_data, met_data, PDWG_BPHSkim_cff::muons, muonsMissing_, muonToken_, L1Analysis::L1AnalysisRecoJetDataFormat::nCaloJets, L1Analysis::L1AnalysisRecoJetDataFormat::nJets, pfJECToken_, pfJetCorrMissing_, pfJetBenchmark_cfi::pfJets, pfJetsMissing_, pfJetToken_, RecoPFMET_cff::pfMet, pfMetMissing_, pfMetToken_, L1Analysis::L1AnalysisRecoMetDataFormat::Reset(), L1Analysis::L1AnalysisRecoJetDataFormat::Reset(), and tree_.

                                                                                        {
   jet_data->Reset();
   met_data->Reset();
 
   // get jets
   edm::Handle<reco::PFJetCollection> pfJets;
   iEvent.getByToken(pfJetToken_, pfJets);
 
   // get calo jets
   edm::Handle<reco::CaloJetCollection> caloJets;
   iEvent.getByToken(caloJetToken_, caloJets);
 
   //get sums
   edm::Handle<reco::PFMETCollection> pfMet;
   iEvent.getByToken(pfMetToken_, pfMet);
 
   // get jet ID
   edm::Handle<edm::ValueMap<reco::JetID> > jetsID;
   iEvent.getByToken(caloJetIDToken_, jetsID);
 
   edm::Handle<reco::JetCorrector> pfJetCorr;
   iEvent.getByToken(pfJECToken_, pfJetCorr);
 
   edm::Handle<reco::JetCorrector> caloJetCorr;
   iEvent.getByToken(caloJECToken_, caloJetCorr);
 
   edm::Handle<reco::CaloMETCollection> caloMet;
   iEvent.getByToken(caloMetToken_, caloMet);
 
   edm::Handle<reco::CaloMETCollection> caloMetBE;
   iEvent.getByToken(caloMetBEToken_, caloMetBE);
 
   // get muons
   edm::Handle<reco::MuonCollection> muons;
   iEvent.getByToken(muonToken_, muons);
 
   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") << "PF Jet Corrector not found.  Branch will not be filled" << std::endl;
     }
     pfJetCorrMissing_ = true;
   }
 
   if (caloJets.isValid()) {
     jet_data->nCaloJets = 0;
 
     doCaloJets(caloJets);
 
   } else {
     if (!caloJetsMissing_) {
       edm::LogWarning("MissingProduct") << "Calo Jets not found.  Branch will not be filled" << std::endl;
     }
     caloJetsMissing_ = true;
   }
 
   if (caloJetCorr.isValid()) {
     doCaloJetCorr(caloJets, caloJetCorr);
 
   } else {
     if (!caloJetCorrMissing_) {
       edm::LogWarning("MissingProduct") << "Calo Jet Corrector not found.  Branch will not be filled" << std::endl;
     }
     caloJetCorrMissing_ = true;
   }
 
   if (!jetsID.isValid()) {
     if (!caloJetIDMissing_) {
       edm::LogWarning("MissingProduct") << "Calo Jet ID not found.  Branch will not be filled" << std::endl;
     }
     caloJetIDMissing_ = true;
   }
 
   if (pfMet.isValid()) {
     doPFMet(pfMet);
 
     if (muons.isValid()) {
       doPFMetNoMu(pfMet, muons);
 
     } else {
       if (!muonsMissing_) {
         edm::LogWarning("MissingProduct") << "Muons not found.  PFMetNoMu branch will not be filled" << std::endl;
       }
       muonsMissing_ = true;
     }
   } else {
     if (!pfMetMissing_) {
       edm::LogWarning("MissingProduct") << "PFMet not found.  Branch will not be filled" << std::endl;
     }
     pfMetMissing_ = true;
   }
 
   if (caloMet.isValid()) {
     doCaloMet(caloMet);
 
   } else {
     if (!caloMetMissing_) {
       edm::LogWarning("MissingProduct") << "CaloMet not found. Branch will not be filled" << std::endl;
     }
     caloMetMissing_ = true;
   }
 
   if (caloMetBE.isValid()) {
     doCaloMetBE(caloMetBE);
 
   } else {
     if (!caloMetBEMissing_) {
       edm::LogWarning("MissingProduct") << "CaloMetBE not found. Branch will not be filled" << std::endl;
     }
     caloMetBEMissing_ = true;
   }
 
   tree_->Fill();
 }

◆ beginJob()

void L1JetRecoTreeProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 524 of file L1JetRecoTreeProducer.cc.

524 {}

◆ caloJetID()

bool L1JetRecoTreeProducer::caloJetID ( const reco::CaloJet & jet )

private

Definition at line 517 of file L1JetRecoTreeProducer.cc.

References createJobs::tmp.

Referenced by doCaloJetCorr(), and doCaloJets().

                                                             {
   bool tmp = true;
 
   return tmp;
 }

◆ doCaloJetCorr()

void L1JetRecoTreeProducer::doCaloJetCorr	(	edm::Handle< reco::CaloJetCollection >	caloJets,
		edm::Handle< reco::JetCorrector >	caloJetCorr
	)

private

Definition at line 404 of file L1JetRecoTreeProducer.cc.

References funct::abs(), L1Analysis::L1AnalysisRecoJetDataFormat::caloCorrFactor, L1Analysis::L1AnalysisRecoJetDataFormat::caloEtCorr, L1Analysis::L1AnalysisRecoMetDataFormat::caloHt, caloJetID(), caloJetIDMissing_, isolatedTracks_cfi::caloJets, reco::JetCorrector::correction(), jet_data, jetetaMax_, jetptThreshold_, maxJet_, met_data, and edm::Handle< T >::product().

Referenced by analyze().

                                                                                    {
   float caloCorrFactor = 1.;
   unsigned int nCaloJets = 0;
 
   met_data->caloHt = 0;
 
   for (auto it = caloJets->begin(); it != caloJets->end() && nCaloJets < maxJet_; ++it) {
     if (!caloJetIDMissing_)
       if (!caloJetID(*it))
         continue;
 
     caloCorrFactor = caloJetCorr.product()->correction(*it);
 
     jet_data->caloEtCorr.push_back(it->et() * caloCorrFactor);
     jet_data->caloCorrFactor.push_back(caloCorrFactor);
 
     nCaloJets++;
 
     if (it->pt() * caloCorrFactor > jetptThreshold_ && std::abs(it->eta()) < jetetaMax_) {
       met_data->caloHt += it->pt() * caloCorrFactor;
     }
   }
 }

◆ doCaloJets()

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

private

Definition at line 306 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

                                                                                 {
   for (auto it = caloJets->begin(); it != caloJets->end() && jet_data->nCaloJets < maxJet_; ++it) {
     if (!caloJetIDMissing_)
       if (!caloJetID(*it))
         continue;
 
     jet_data->caloEt.push_back(it->et());
     jet_data->caloEta.push_back(it->eta());
     jet_data->caloPhi.push_back(it->phi());
     jet_data->caloE.push_back(it->energy());
 
     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->n60.push_back(it->n60());
 
     jet_data->nCaloJets++;
   }
 }

◆ doCaloMet()

void L1JetRecoTreeProducer::doCaloMet ( edm::Handle< reco::CaloMETCollection > caloMet )

private

Definition at line 471 of file L1JetRecoTreeProducer.cc.

References heavyionUCCDQM_cfi::caloMet, L1Analysis::L1AnalysisRecoMetDataFormat::caloMet, L1Analysis::L1AnalysisRecoMetDataFormat::caloMetPhi, L1Analysis::L1AnalysisRecoMetDataFormat::caloSumEt, reco::LeafCandidate::et(), met_data, reco::LeafCandidate::phi(), and reco::MET::sumEt().

Referenced by analyze().

                                                                               {
   const reco::CaloMETCollection* metCol = caloMet.product();
   const reco::CaloMET theMet = metCol->front();
 
   met_data->caloMet = theMet.et();
   met_data->caloMetPhi = theMet.phi();
   met_data->caloSumEt = theMet.sumEt();
 }

◆ doCaloMetBE()

void L1JetRecoTreeProducer::doCaloMetBE ( edm::Handle< reco::CaloMETCollection > caloMetBE )

private

Definition at line 480 of file L1JetRecoTreeProducer.cc.

References CaloMET_cfi::caloMetBE, L1Analysis::L1AnalysisRecoMetDataFormat::caloMetBE, L1Analysis::L1AnalysisRecoMetDataFormat::caloMetPhiBE, L1Analysis::L1AnalysisRecoMetDataFormat::caloSumEtBE, reco::LeafCandidate::et(), met_data, reco::LeafCandidate::phi(), and reco::MET::sumEt().

Referenced by analyze().

                                                                                   {
   const reco::CaloMETCollection* metCol = caloMetBE.product();
   const reco::CaloMET theMet = metCol->front();
 
   met_data->caloMetBE = theMet.et();
   met_data->caloMetPhiBE = theMet.phi();
   met_data->caloSumEtBE = theMet.sumEt();
 }

◆ doPFJetCorr()

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

private

Definition at line 369 of file L1JetRecoTreeProducer.cc.

References funct::abs(), reco::JetCorrector::correction(), L1Analysis::L1AnalysisRecoJetDataFormat::corrFactor, L1Analysis::L1AnalysisRecoJetDataFormat::etCorr, L1Analysis::L1AnalysisRecoMetDataFormat::Ht, jet_data, jetetaMax_, jetptThreshold_, maxJet_, met_data, L1Analysis::L1AnalysisRecoMetDataFormat::mHt, L1Analysis::L1AnalysisRecoMetDataFormat::mHtPhi, l1ctJetFileWriter_cfi::nJets, pfJetID(), pfJetBenchmark_cfi::pfJets, and edm::Handle< T >::product().

Referenced by analyze().

                                                                                {
   float corrFactor = 1.;
   unsigned int nJets = 0;
 
   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() && nJets < maxJet_; ++it) {
     if (!pfJetID(*it))
       continue;
 
     corrFactor = pfJetCorr.product()->correction(*it);
 
     jet_data->etCorr.push_back(it->et() * corrFactor);
     jet_data->corrFactor.push_back(corrFactor);
 
     nJets++;
 
     if (it->pt() * corrFactor > jetptThreshold_ && std::abs(it->eta()) < jetetaMax_) {
       mHx += -1. * it->px() * corrFactor;
       mHy += -1. * it->py() * corrFactor;
       met_data->Ht += it->pt() * corrFactor;
     }
   }
 
   TVector2 tv2 = TVector2(mHx, mHy);
   met_data->mHt = tv2.Mod();
   met_data->mHtPhi = tv2.Phi();
 }

◆ doPFJets()

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

private

Definition at line 334 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

                                                                           {
   for (auto it = pfJets->begin(); it != pfJets->end() && jet_data->nJets < maxJet_; ++it) {
     if (!pfJetID(*it))
       continue;
 
     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->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++;
   }
 }

◆ doPFMet()

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

private

Definition at line 429 of file L1JetRecoTreeProducer.cc.

References reco::LeafCandidate::et(), L1Analysis::L1AnalysisRecoMetDataFormat::met, met_data, L1Analysis::L1AnalysisRecoMetDataFormat::metPhi, L1Analysis::L1AnalysisRecoMetDataFormat::metPx, L1Analysis::L1AnalysisRecoMetDataFormat::metPy, RecoPFMET_cff::pfMet, reco::LeafCandidate::phi(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), reco::MET::sumEt(), and L1Analysis::L1AnalysisRecoMetDataFormat::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();
   met_data->metPx = theMet.px();
   met_data->metPy = theMet.py();
 }

◆ doPFMetNoMu()

void L1JetRecoTreeProducer::doPFMetNoMu	(	edm::Handle< reco::PFMETCollection >	pfMet,
		edm::Handle< reco::MuonCollection >	muons
	)

private

Definition at line 440 of file L1JetRecoTreeProducer.cc.

References reco::LeafCandidate::et(), met_data, PDWG_BPHSkim_cff::muons, RecoPFMET_cff::pfMet, L1Analysis::L1AnalysisRecoMetDataFormat::pfMetNoMu, L1Analysis::L1AnalysisRecoMetDataFormat::pfMetNoMuPhi, L1Analysis::L1AnalysisRecoMetDataFormat::pfMetNoMuPx, L1Analysis::L1AnalysisRecoMetDataFormat::pfMetNoMuPy, reco::LeafCandidate::phi(), reco::LeafCandidate::px(), reco::LeafCandidate::py(), and reco::LeafCandidate::setP4().

Referenced by analyze().

                                                                              {
   const reco::PFMETCollection* metCol = pfMet.product();
   const reco::PFMET theMet = metCol->front();
   reco::PFMET thePFMetNoMu = metCol->front();
 
   double pfMetNoMuPx = theMet.px();
   double pfMetNoMuPy = theMet.py();
 
   double muPx(0.), muPy(0.);
 
   for (auto it = muons->begin(); it != muons->end(); ++it) {
     if (it->isPFMuon()) {
       muPx += it->px();
       muPy += it->py();
     }
   }
 
   pfMetNoMuPx += muPx;
   pfMetNoMuPy += muPy;
 
   math::XYZTLorentzVector pfMetNoMuP4(pfMetNoMuPx, pfMetNoMuPy, 0, hypot(pfMetNoMuPx, pfMetNoMuPy));
 
   thePFMetNoMu.setP4(pfMetNoMuP4);
 
   met_data->pfMetNoMu = thePFMetNoMu.et();
   met_data->pfMetNoMuPhi = thePFMetNoMu.phi();
   met_data->pfMetNoMuPx = thePFMetNoMu.px();
   met_data->pfMetNoMuPy = thePFMetNoMu.py();
 }

◆ endJob()

void L1JetRecoTreeProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 527 of file L1JetRecoTreeProducer.cc.

527 {}

◆ pfJetID()

bool L1JetRecoTreeProducer::pfJetID ( const reco::PFJet & jet )

private

Definition at line 489 of file L1JetRecoTreeProducer.cc.

References funct::abs(), metsig::jet, and createJobs::tmp.

Referenced by doPFJetCorr(), and doPFJets().

                                                         {
   bool tmp = true;
   if (std::abs(jet.eta()) < 2.7) {
     tmp &= jet.neutralHadronEnergyFraction() < 0.9;
     tmp &= jet.neutralEmEnergyFraction() < 0.9;
     tmp &= (jet.chargedMultiplicity() + jet.neutralMultiplicity()) > 1;
     tmp &= jet.muonEnergyFraction() < 0.8;
     tmp &= jet.chargedHadronEnergyFraction() > 0.0;
     tmp &= jet.chargedMultiplicity() > 0;
     tmp &= jet.chargedEmEnergyFraction() < 0.9;
   }
   if (std::abs(jet.eta()) > 2.7 && std::abs(jet.eta()) < 3.0) {
     tmp &= jet.neutralEmEnergyFraction() > 0.01;
     tmp &= jet.neutralHadronEnergyFraction() < 0.98;
     tmp &= jet.neutralMultiplicity() > 2;
   }
   if (std::abs(jet.eta()) > 3.0) {
     tmp &= jet.neutralEmEnergyFraction() < 0.9;
     tmp &= jet.neutralMultiplicity() > 10;
   }
 
   // our custom selection
   //tmp &= jet.muonMultiplicity() == 0;
   //tmp &= jet.electronMultiplicity() == 0;
 
   return tmp;
 }

Member Data Documentation

◆ caloJECToken_

edm::EDGetTokenT<reco::JetCorrector> L1JetRecoTreeProducer::caloJECToken_

private

Definition at line 94 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

◆ caloJetCorrMissing_

bool L1JetRecoTreeProducer::caloJetCorrMissing_

private

Definition at line 112 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

◆ caloJetIDMissing_

bool L1JetRecoTreeProducer::caloJetIDMissing_

private

Definition at line 114 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), doCaloJetCorr(), and doCaloJets().

◆ caloJetIDToken_

edm::EDGetTokenT<edm::ValueMap<reco::JetID> > L1JetRecoTreeProducer::caloJetIDToken_

private

Definition at line 92 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

◆ caloJetsMissing_

bool L1JetRecoTreeProducer::caloJetsMissing_

private

Definition at line 113 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

◆ caloJetToken_

edm::EDGetTokenT<reco::CaloJetCollection> L1JetRecoTreeProducer::caloJetToken_

private

Definition at line 91 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

◆ caloMetBEMissing_

bool L1JetRecoTreeProducer::caloMetBEMissing_

private

Definition at line 117 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

◆ caloMetBEToken_

edm::EDGetTokenT<reco::CaloMETCollection> L1JetRecoTreeProducer::caloMetBEToken_

private

Definition at line 98 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

◆ caloMetMissing_

bool L1JetRecoTreeProducer::caloMetMissing_

private

Definition at line 116 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

◆ caloMetToken_

edm::EDGetTokenT<reco::CaloMETCollection> L1JetRecoTreeProducer::caloMetToken_

private

Definition at line 97 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

◆ jet_data

L1Analysis::L1AnalysisRecoJetDataFormat* L1JetRecoTreeProducer::jet_data

Definition at line 82 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), doCaloJetCorr(), doCaloJets(), doPFJetCorr(), doPFJets(), and L1JetRecoTreeProducer().

◆ jetetaMax_

double L1JetRecoTreeProducer::jetetaMax_

private

Definition at line 105 of file L1JetRecoTreeProducer.cc.

Referenced by doCaloJetCorr(), doPFJetCorr(), and L1JetRecoTreeProducer().

◆ jetptThreshold_

double L1JetRecoTreeProducer::jetptThreshold_

private

Definition at line 104 of file L1JetRecoTreeProducer.cc.

Referenced by doCaloJetCorr(), doPFJetCorr(), and L1JetRecoTreeProducer().

◆ maxCl_

unsigned int L1JetRecoTreeProducer::maxCl_

private

Definition at line 106 of file L1JetRecoTreeProducer.cc.

◆ maxJet_

unsigned int L1JetRecoTreeProducer::maxJet_

private

Definition at line 107 of file L1JetRecoTreeProducer.cc.

Referenced by doCaloJetCorr(), doCaloJets(), doPFJetCorr(), doPFJets(), and L1JetRecoTreeProducer().

◆ maxTrk_

unsigned int L1JetRecoTreeProducer::maxTrk_

private

Definition at line 109 of file L1JetRecoTreeProducer.cc.

◆ maxVtx_

unsigned int L1JetRecoTreeProducer::maxVtx_

private

Definition at line 108 of file L1JetRecoTreeProducer.cc.

◆ met_data

L1Analysis::L1AnalysisRecoMetDataFormat* L1JetRecoTreeProducer::met_data

Definition at line 83 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), doCaloJetCorr(), doCaloMet(), doCaloMetBE(), doPFJetCorr(), doPFMet(), doPFMetNoMu(), and L1JetRecoTreeProducer().

◆ muonsMissing_

bool L1JetRecoTreeProducer::muonsMissing_

private

Definition at line 119 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

◆ muonToken_

edm::EDGetTokenT<reco::MuonCollection> L1JetRecoTreeProducer::muonToken_

private

Definition at line 100 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

◆ pfJECToken_

edm::EDGetTokenT<reco::JetCorrector> L1JetRecoTreeProducer::pfJECToken_

private

Definition at line 93 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

◆ pfJetCorrMissing_

bool L1JetRecoTreeProducer::pfJetCorrMissing_

private

Definition at line 111 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

◆ pfJetsMissing_

bool L1JetRecoTreeProducer::pfJetsMissing_

private

Definition at line 103 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

◆ pfJetToken_

edm::EDGetTokenT<reco::PFJetCollection> L1JetRecoTreeProducer::pfJetToken_

private

Definition at line 90 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

◆ pfMetMissing_

bool L1JetRecoTreeProducer::pfMetMissing_

private

Definition at line 115 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

◆ pfMetToken_

edm::EDGetTokenT<reco::PFMETCollection> L1JetRecoTreeProducer::pfMetToken_

private

Definition at line 96 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

◆ tree_

TTree* L1JetRecoTreeProducer::tree_

private

Definition at line 87 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

Public Member Functions

Public Attributes

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ L1JetRecoTreeProducer()

◆ ~L1JetRecoTreeProducer()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ caloJetID()

◆ doCaloJetCorr()

◆ doCaloJets()

◆ doCaloMet()

◆ doCaloMetBE()

◆ doPFJetCorr()

◆ doPFJets()

◆ doPFMet()

◆ doPFMetNoMu()

◆ endJob()

◆ pfJetID()

Member Data Documentation

◆ caloJECToken_

◆ caloJetCorrMissing_

◆ caloJetIDMissing_

◆ caloJetIDToken_

◆ caloJetsMissing_

◆ caloJetToken_

◆ caloMetBEMissing_

◆ caloMetBEToken_

◆ caloMetMissing_

◆ caloMetToken_

◆ jet_data

◆ jetetaMax_

◆ jetptThreshold_

◆ maxCl_

◆ maxJet_

◆ maxTrk_

◆ maxVtx_

◆ met_data

◆ muonsMissing_

◆ muonToken_

◆ pfJECToken_

◆ pfJetCorrMissing_

◆ pfJetsMissing_

◆ pfJetToken_

◆ pfMetMissing_

◆ pfMetToken_

◆ tree_