L1JetRecoTreeProducer Class Reference

#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
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
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
virtual void	analyze (const edm::Event &, const edm::EventSetup &)
virtual void	beginJob (void)
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)
virtual void	endJob ()
bool	jetId (const reco::PFJet &jet)

Private Attributes
bool	caloMetBEMissing_
edm::EDGetTokenT< reco::CaloMETCollection >	caloMetBEToken_
bool	caloMetMissing_
edm::EDGetTokenT< reco::CaloMETCollection >	caloMetToken_
edm::Service< TFileService >	fs_
edm::EDGetTokenT< reco::JetCorrector >	jecToken_
double	jetetaMax_
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 58 of file L1JetRecoTreeProducer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 119 of file L1JetRecoTreeProducer.cc.

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

                                                                          :
  pfJetsMissing_(false),
  pfMetMissing_(false),
  pfJetCorrMissing_(false),
  caloMetMissing_(false),
  caloMetBEMissing_(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("pfMetT1")));
  caloMetToken_ = consumes<reco::CaloMETCollection>(iConfig.getUntrackedParameter("caloMetToken",edm::InputTag("caloMet")));
  caloMetBEToken_ = consumes<reco::CaloMETCollection>(iConfig.getUntrackedParameter("caloMetBEToken",edm::InputTag("caloMetBE")));

  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

(

)

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

Definition at line 165 of file L1JetRecoTreeProducer.cc.

References CaloMET_cfi::caloMet, CaloMET_cfi::caloMetBE, caloMetBEMissing_, caloMetBEToken_, caloMetMissing_, caloMetToken_, doCaloMet(), doCaloMetBE(), doPFJetCorr(), doPFJets(), doPFMet(), edm::Event::getByToken(), edm::HandleBase::isValid(), jecToken_, jet_data, met_data, L1Analysis::L1AnalysisRecoJetDataFormat::nJets, pfJetCorrMissing_, pfJetBenchmark_cfi::pfJets, pfJetsMissing_, pfJetToken_, PFMET_cfi::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 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);

  edm::Handle<reco::CaloMETCollection> caloMet;
  iEvent.getByToken(caloMetToken_, caloMet);

  edm::Handle<reco::CaloMETCollection> caloMetBE;
  iEvent.getByToken(caloMetBEToken_, caloMetBE);

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

  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  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 456 of file L1JetRecoTreeProducer.cc.

{
}

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

private

Definition at line 250 of file L1JetRecoTreeProducer.cc.

References L1Analysis::L1AnalysisRecoJetDataFormat::e, L1Analysis::L1AnalysisRecoJetDataFormat::eEmEB, L1Analysis::L1AnalysisRecoJetDataFormat::eEmEE, L1Analysis::L1AnalysisRecoJetDataFormat::eEMF, L1Analysis::L1AnalysisRecoJetDataFormat::eEmHF, L1Analysis::L1AnalysisRecoJetDataFormat::eHadHB, L1Analysis::L1AnalysisRecoJetDataFormat::eHadHE, L1Analysis::L1AnalysisRecoJetDataFormat::eHadHF, L1Analysis::L1AnalysisRecoJetDataFormat::eHadHO, L1Analysis::L1AnalysisRecoJetDataFormat::eMaxEcalTow, L1Analysis::L1AnalysisRecoJetDataFormat::eMaxHcalTow, L1Analysis::L1AnalysisRecoJetDataFormat::et, L1Analysis::L1AnalysisRecoJetDataFormat::eta, L1Analysis::L1AnalysisRecoJetDataFormat::isPF, jet_data, maxJet_, L1Analysis::L1AnalysisRecoJetDataFormat::nJets, L1Analysis::L1AnalysisRecoJetDataFormat::phi, and L1Analysis::L1AnalysisRecoJetDataFormat::towerArea.

                                                                           {


  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::doCaloMet ( edm::Handle< reco::CaloMETCollection >  caloMet )

private

Definition at line 403 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 415 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 325 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 = new TVector2(mHx,mHy);
  met_data->mHt    = tv2->Mod();
  met_data->mHtPhi = tv2->Phi();

  // 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;
  // unsigned int 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 284 of file L1JetRecoTreeProducer.cc.

References L1Analysis::L1AnalysisRecoJetDataFormat::cemef, L1Analysis::L1AnalysisRecoJetDataFormat::chef, L1Analysis::L1AnalysisRecoJetDataFormat::chMult, L1Analysis::L1AnalysisRecoJetDataFormat::cmef, L1Analysis::L1AnalysisRecoJetDataFormat::cMult, L1Analysis::L1AnalysisRecoJetDataFormat::e, L1Analysis::L1AnalysisRecoJetDataFormat::eef, L1Analysis::L1AnalysisRecoJetDataFormat::elMult, L1Analysis::L1AnalysisRecoJetDataFormat::et, L1Analysis::L1AnalysisRecoJetDataFormat::eta, L1Analysis::L1AnalysisRecoJetDataFormat::hfemef, L1Analysis::L1AnalysisRecoJetDataFormat::hfemMult, L1Analysis::L1AnalysisRecoJetDataFormat::hfhef, L1Analysis::L1AnalysisRecoJetDataFormat::hfhMult, L1Analysis::L1AnalysisRecoJetDataFormat::isPF, jet_data, maxJet_, L1Analysis::L1AnalysisRecoJetDataFormat::mef, L1Analysis::L1AnalysisRecoJetDataFormat::muMult, L1Analysis::L1AnalysisRecoJetDataFormat::nemef, L1Analysis::L1AnalysisRecoJetDataFormat::nhef, L1Analysis::L1AnalysisRecoJetDataFormat::nhMult, L1Analysis::L1AnalysisRecoJetDataFormat::nJets, L1Analysis::L1AnalysisRecoJetDataFormat::nMult, L1Analysis::L1AnalysisRecoJetDataFormat::pef, L1Analysis::L1AnalysisRecoJetDataFormat::phi, and L1Analysis::L1AnalysisRecoJetDataFormat::phMult.

Referenced by analyze().

                                                                     {
  

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

  }

}

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

private

Definition at line 389 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(), 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();
  met_data->metPx   = theMet.px();
  met_data->metPy   = theMet.py();

}

void L1JetRecoTreeProducer::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 462 of file L1JetRecoTreeProducer.cc.

References DEFINE_FWK_MODULE.

                              {
}

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

private

Definition at line 427 of file L1JetRecoTreeProducer.cc.

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

                                                 {

  bool tmp = true;

  tmp &= jet.neutralHadronEnergyFraction() < 0.9 ;
  tmp &= jet.neutralEmEnergyFraction() < 0.9 ;
  tmp &= (jet.chargedMultiplicity() + jet.neutralMultiplicity()) > 1 ;
  tmp &= jet.muonEnergyFraction() < 0.8 ;
  if (fabs(jet.eta()) < 2.4) {
    tmp &= jet.chargedHadronEnergyFraction() > 0.0 ;
    tmp &= jet.chargedMultiplicity() > 0 ;
    tmp &= jet.chargedEmEnergyFraction() < 0.9 ;
  }
  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;

}

Member Data Documentation

bool L1JetRecoTreeProducer::caloMetBEMissing_

private

Definition at line 114 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

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

private

Definition at line 99 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

bool L1JetRecoTreeProducer::caloMetMissing_

private

Definition at line 113 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

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

private

Definition at line 98 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

edm::Service<TFileService> L1JetRecoTreeProducer::fs_

private

Definition at line 86 of file L1JetRecoTreeProducer.cc.

Referenced by L1JetRecoTreeProducer().

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

private

Definition at line 95 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

L1Analysis::L1AnalysisRecoJetDataFormat* L1JetRecoTreeProducer::jet_data

Definition at line 80 of file L1JetRecoTreeProducer.cc.

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

double L1JetRecoTreeProducer::jetetaMax_

private

Definition at line 105 of file L1JetRecoTreeProducer.cc.

Referenced by doPFJetCorr(), and L1JetRecoTreeProducer().

double L1JetRecoTreeProducer::jetptThreshold_

private

Definition at line 104 of file L1JetRecoTreeProducer.cc.

Referenced by doPFJetCorr(), and L1JetRecoTreeProducer().

unsigned int L1JetRecoTreeProducer::maxCl_

private

Definition at line 106 of file L1JetRecoTreeProducer.cc.

unsigned int L1JetRecoTreeProducer::maxJet_

private

Definition at line 107 of file L1JetRecoTreeProducer.cc.

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

unsigned int L1JetRecoTreeProducer::maxTrk_

private

Definition at line 109 of file L1JetRecoTreeProducer.cc.

unsigned int L1JetRecoTreeProducer::maxVtx_

private

Definition at line 108 of file L1JetRecoTreeProducer.cc.

L1Analysis::L1AnalysisRecoMetDataFormat* L1JetRecoTreeProducer::met_data

Definition at line 81 of file L1JetRecoTreeProducer.cc.

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

bool L1JetRecoTreeProducer::pfJetCorrMissing_

private

Definition at line 112 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

bool L1JetRecoTreeProducer::pfJetsMissing_

private

Definition at line 103 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

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

private

Definition at line 92 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

bool L1JetRecoTreeProducer::pfMetMissing_

private

Definition at line 111 of file L1JetRecoTreeProducer.cc.

Referenced by analyze().

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

private

Definition at line 97 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().

TTree* L1JetRecoTreeProducer::tree_

private

Definition at line 89 of file L1JetRecoTreeProducer.cc.

Referenced by analyze(), and L1JetRecoTreeProducer().