#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

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)

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 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 60 of file L1JetRecoTreeProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 130 of file L1JetRecoTreeProducer.cc.

References caloJECToken_, caloJetIDToken_, caloJetToken_, caloMetBEToken_, caloMetToken_, fs_, edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), 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 169 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 183 of file L1JetRecoTreeProducer.cc.

References caloJECToken_, caloJetCorrMissing_, caloJetIDMissing_, caloJetIDToken_, isolatedTracks_cfi::caloJets, caloJetsMissing_, caloJetToken_, CaloMET_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, nano_cff::muons, muonsMissing_, muonToken_, L1Analysis::L1AnalysisRecoJetDataFormat::nCaloJets, L1Analysis::L1AnalysisRecoJetDataFormat::nJets, pfJECToken_, pfJetCorrMissing_, pfJetBenchmark_cfi::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 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 585 of file L1JetRecoTreeProducer.cc.

586 {

587 }

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

private

Definition at line 573 of file L1JetRecoTreeProducer.cc.

References tmp.

Referenced by doCaloJets().

                                                        {
 
   bool tmp = true;
 
   return tmp;
 
 }

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

private

Definition at line 441 of file L1JetRecoTreeProducer.cc.

References L1Analysis::L1AnalysisRecoJetDataFormat::caloCorrFactor, L1Analysis::L1AnalysisRecoJetDataFormat::caloEtCorr, L1Analysis::L1AnalysisRecoMetDataFormat::caloHt, 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) {
 
     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 317 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 518 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 530 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 397 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, 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) {
     
     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 353 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 468 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 482 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 591 of file L1JetRecoTreeProducer.cc.

References DEFINE_FWK_MODULE.

591 {

592 }

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

private

Definition at line 542 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 tmp.

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