#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::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)

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_

edm::Service< TFileService >	fs_

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

Definition at line 58 of file L1JetRecoTreeProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 125 of file L1JetRecoTreeProducer.cc.

References caloJECToken_, caloJetIDToken_, caloJetToken_, caloMetBEToken_, caloMetToken_, fs_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), HLT_2018_cff::InputTag, jet_data, jetetaMax_, jetptThreshold_, 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")));
 
   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
   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 ( )

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

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

overrideprivate

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(), edm::Event::getByToken(), 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::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 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();
 }

void L1JetRecoTreeProducer::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 524 of file L1JetRecoTreeProducer.cc.

524 {}

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;
 }

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

private

Definition at line 404 of file L1JetRecoTreeProducer.cc.

References L1Analysis::L1AnalysisRecoJetDataFormat::caloCorrFactor, L1Analysis::L1AnalysisRecoJetDataFormat::caloEtCorr, L1Analysis::L1AnalysisRecoMetDataFormat::caloHt, caloJetID(), caloJetIDMissing_, 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_ && fabs(it->eta()) < jetetaMax_) {
       met_data->caloHt += it->pt() * caloCorrFactor;
     }
   }
 }

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++;
   }
 }

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

private

Definition at line 471 of file L1JetRecoTreeProducer.cc.

References L1Analysis::L1AnalysisRecoMetDataFormat::caloMet, L1Analysis::L1AnalysisRecoMetDataFormat::caloMetPhi, L1Analysis::L1AnalysisRecoMetDataFormat::caloSumEt, reco::LeafCandidate::et(), met_data, reco::LeafCandidate::phi(), edm::Handle< T >::product(), 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();
 }

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

private

Definition at line 480 of file L1JetRecoTreeProducer.cc.

References L1Analysis::L1AnalysisRecoMetDataFormat::caloMetBE, L1Analysis::L1AnalysisRecoMetDataFormat::caloMetPhiBE, L1Analysis::L1AnalysisRecoMetDataFormat::caloSumEtBE, reco::LeafCandidate::et(), met_data, reco::LeafCandidate::phi(), edm::Handle< T >::product(), 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();
 }

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

private

Definition at line 369 of file L1JetRecoTreeProducer.cc.

References 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, pfJetID(), 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_ && fabs(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();
 }

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++;
   }
 }

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, reco::LeafCandidate::phi(), edm::Handle< T >::product(), 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();
 }

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, L1Analysis::L1AnalysisRecoMetDataFormat::pfMetNoMu, L1Analysis::L1AnalysisRecoMetDataFormat::pfMetNoMuPhi, L1Analysis::L1AnalysisRecoMetDataFormat::pfMetNoMuPx, L1Analysis::L1AnalysisRecoMetDataFormat::pfMetNoMuPy, reco::LeafCandidate::phi(), edm::Handle< T >::product(), 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();
 }

void L1JetRecoTreeProducer::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 527 of file L1JetRecoTreeProducer.cc.

References DEFINE_FWK_MODULE.

527 {}

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

private

Definition at line 489 of file L1JetRecoTreeProducer.cc.

References reco::PFJet::chargedEmEnergyFraction(), reco::PFJet::chargedHadronEnergyFraction(), reco::PFJet::chargedMultiplicity(), reco::LeafCandidate::eta(), reco::PFJet::muonEnergyFraction(), reco::PFJet::neutralEmEnergyFraction(), reco::PFJet::neutralHadronEnergyFraction(), reco::PFJet::neutralMultiplicity(), and createJobs::tmp.

Referenced by doPFJetCorr(), and doPFJets().

                                                         {
   bool tmp = true;
   if (fabs(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 (fabs(jet.eta()) > 2.7 && fabs(jet.eta()) < 3.0) {
     tmp &= jet.neutralEmEnergyFraction() > 0.01;
     tmp &= jet.neutralHadronEnergyFraction() < 0.98;
     tmp &= jet.neutralMultiplicity() > 2;
   }
   if (fabs(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;
 }