JetValidation Class Reference

#include <JetValidation.h>

Inheritance diagram for JetValidation:

Public Member Functions
	JetValidation (edm::ParameterSet const &cfg)
Public Member Functions inherited from edm::EDAnalyzer
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzer ()
ModuleDescription const &	moduleDescription () const
std::string	workerType () const
virtual	~EDAnalyzer ()
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
	EDConsumerBase ()
	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 &	itemsToGetFromEvent () 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)

Private Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &iSetup)
void	beginJob ()
void	endJob ()
void	FillHist1D (const TString &histName, const Double_t &x)
void	FillHist2D (const TString &histName, const Double_t &x, const Double_t &y)
void	FillHistProfile (const TString &histName, const Double_t &x, const Double_t &y)

Private Attributes
std::string	calAlgo
double	dRmatch
std::string	genAlgo
std::string	histoFileName
std::string	jetTracksAssociator
TFile *	m_file
std::map< TString, TH1 * >	m_HistNames1D
std::map< TString, TH2 * >	m_HistNames2D
std::map< TString, TProfile * >	m_HistNamesProfile
bool	MCarlo
int	Njets
double	PtMin

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

Definition at line 17 of file JetValidation.h.

Constructor & Destructor Documentation

JetValidation::JetValidation

(

edm::ParameterSet const &

cfg

)

Definition at line 24 of file JetValidation.cc.

References CaloRecoTauDiscriminationAgainstMuon_cfi::dRmatch, edm::ParameterSet::getParameter(), mix_2012_lumiLevel_15_20_50ns_PoissonOOTPU_cfi::histoFileName, ECF_cfi::Njets, and HIPixelMedianVertex_cfi::PtMin.

{
  dRmatch             = cfg.getParameter<double> ("dRmatch");
  PtMin               = cfg.getParameter<double> ("PtMin");
  Njets               = cfg.getParameter<int> ("Njets");
  MCarlo              = cfg.getParameter<bool> ("MCarlo");
  genAlgo             = cfg.getParameter<string> ("genAlgo");
  calAlgo             = cfg.getParameter<string> ("calAlgo");
  jetTracksAssociator = cfg.getParameter<string> ("jetTracksAssociator"); 
  histoFileName       = cfg.getParameter<string> ("histoFileName");
}

Member Function Documentation

void JetValidation::analyze ( edm::Event const &  e, edm::EventSetup const &  iSetup  )

privatevirtual

Implements edm::EDAnalyzer.

Definition at line 81 of file JetValidation.cc.

References gather_cfg::cout, deltaR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, MillePedeFileConverter_cfg::e, reco::CaloJet::emEnergyInEB(), reco::CaloJet::emEnergyInEE(), reco::CaloJet::emEnergyInHF(), reco::LeafCandidate::energy(), stringResolutionProvider_cfi::eta, reco::LeafCandidate::eta(), edm::EventID::event(), edm::Event::getByLabel(), reco::CaloJet::hadEnergyInHB(), reco::CaloJet::hadEnergyInHE(), reco::CaloJet::hadEnergyInHF(), reco::CaloJet::hadEnergyInHO(), edm::EventBase::id(), impactParameterTagInfos_cfi::jetTracks, ResonanceBuilder::mass, ECF_cfi::Njets, nTracks(), EnergyCorrector::pt, reco::LeafCandidate::pt(), HIPixelMedianVertex_cfi::PtMin, edm::EventID::run(), reco::JetTracksAssociation::tracksNumber(), and reco::JetTracksAssociation::tracksP4().

{
  math::XYZTLorentzVector p4jet[2];
  int jetInd,jetCounter,nTracks;
  double dRmin,dR,e,eta,emEB,emEE,emHF,hadHB,hadHE,hadHO,hadHF,pt,phi,pthat,chf;
  Handle<CaloJetCollection> caljets;
  Handle<GenJetCollection> genjets;
  Handle<double> genEventScale;
  Handle<JetTracksAssociation::Container> jetTracks;
  CaloJetCollection::const_iterator i_caljet;
  GenJetCollection::const_iterator i_genjet;
  evt.getByLabel(calAlgo,caljets);
  evt.getByLabel(jetTracksAssociator,jetTracks);
  jetInd = 0;
  jetCounter = 0;
  if (caljets->size()==0)
    cout<<"WARNING: NO calo jets in event "<<evt.id().event()<<", Run "<<evt.id().run()<<" !!!!"<<endl;
  for(i_caljet = caljets->begin(); i_caljet != caljets->end() && jetInd<Njets; ++i_caljet) 
    {
      e = i_caljet->energy();
      pt = i_caljet->pt();
      phi = i_caljet->phi(); 
      eta = i_caljet->eta();
      emEB = i_caljet->emEnergyInEB();
      emEE = i_caljet->emEnergyInEE();
      emHF = i_caljet->emEnergyInHF();
      hadHB = i_caljet->hadEnergyInHB();
      hadHE = i_caljet->hadEnergyInHE(); 
      hadHO = i_caljet->hadEnergyInHO(); 
      hadHF = i_caljet->hadEnergyInHF();  
      nTracks = JetTracksAssociation::tracksNumber(*jetTracks,*i_caljet);
      chf = (JetTracksAssociation::tracksP4(*jetTracks,*i_caljet)).pt()/pt;
      if (jetInd<2)
        p4jet[jetInd] = i_caljet->p4();
      if (pt>PtMin)
        {  
          FillHist1D("ptCalo",pt);
          FillHist1D("etaCalo",eta);
          FillHist1D("phiCalo",phi);
          FillHist1D("emEnergyFraction",i_caljet->emEnergyFraction());
      FillHist1D("nTracks",nTracks);
      FillHist1D("chargeFraction",chf); 
          FillHist1D("emEnergyInEB",emEB); 
          FillHist1D("emEnergyInEE",emEE); 
          FillHist1D("emEnergyInHF",emHF); 
          FillHist1D("hadEnergyInHB",hadHB); 
          FillHist1D("hadEnergyInHE",hadHE); 
          FillHist1D("hadEnergyInHF",hadHF); 
          FillHist1D("hadEnergyInHO",hadHO);
          FillHistProfile("EBfractionVsEta",eta,emEB/e);
          FillHistProfile("EEfractionVsEta",eta,emEE/e);
          FillHistProfile("HBfractionVsEta",eta,hadHB/e);
          FillHistProfile("HOfractionVsEta",eta,hadHO/e);
          FillHistProfile("HEfractionVsEta",eta,hadHE/e);
          FillHistProfile("HFfractionVsEta",eta,(hadHF+emHF)/e);
          FillHistProfile("CaloEnergyVsEta",eta,e);
          FillHistProfile("emEnergyVsEta",eta,emEB+emEE+emHF);
          FillHistProfile("hadEnergyVsEta",eta,hadHB+hadHO+hadHE+hadHF);
          jetCounter++;
        }
      jetInd++;
    }
  FillHist1D("CaloJetMulti",jetCounter);
  if (jetInd>1)
    FillHist1D("m2jCalo",(p4jet[0]+p4jet[1]).mass());
  if (MCarlo)
    { 
      evt.getByLabel(genAlgo,genjets);  
      evt.getByLabel("genEventScale",genEventScale);
      pthat = *genEventScale;
      FillHist1D("ptHat",pthat);
      CaloJet MatchedJet;
      jetInd = 0;
      if (genjets->size()==0)
        cout<<"WARNING: NO gen jets in event "<<evt.id().event()<<", Run "<<evt.id().run()<<" !!!!"<<endl;
      for(i_genjet = genjets->begin(); i_genjet != genjets->end() && jetInd<Njets; ++i_genjet) 
        {
          if (jetInd<2)
            p4jet[jetInd] = i_genjet->p4();
          FillHist1D("ptGen",i_genjet->pt());
          FillHist1D("etaGen",i_genjet->eta());
          FillHist1D("phiGen",i_genjet->phi());
          FillHistProfile("GenEnergyVsEta",i_genjet->eta(),i_genjet->energy());
          dRmin=1000.0;
          for(i_caljet = caljets->begin(); i_caljet != caljets->end(); ++i_caljet)
            {
              dR = deltaR(i_caljet->eta(),i_caljet->phi(),i_genjet->eta(),i_genjet->phi());
              if (dR<dRmin)
                {
                  dRmin = dR;           
              MatchedJet = *i_caljet;       
                }
            }
          FillHist1D("dR",dRmin); 
          e = MatchedJet.energy();
          pt = MatchedJet.pt();
          eta = MatchedJet.eta();
          emEB = MatchedJet.emEnergyInEB();
          emEE = MatchedJet.emEnergyInEE();
          emHF = MatchedJet.emEnergyInHF();
          hadHB = MatchedJet.hadEnergyInHB();
          hadHE = MatchedJet.hadEnergyInHE(); 
          hadHO = MatchedJet.hadEnergyInHO(); 
          hadHF = MatchedJet.hadEnergyInHF();    
          if (dRmin<dRmatch && pt>PtMin)
            {
              FillHistProfile("CaloErespVsEta",eta,e/i_genjet->energy());
              FillHistProfile("emErespVsEta",eta,(emEB+emEE+emHF)/i_genjet->energy());
              FillHistProfile("hadErespVsEta",eta,(hadHB+hadHO+hadHE+hadHF)/i_genjet->energy());
              if (fabs(i_genjet->eta())<1.)
                FillHistProfile("respVsPtBarrel",i_genjet->pt(),pt/i_genjet->pt()); 
            }
          jetInd++;
        } 
      FillHist1D("GenJetMulti",jetInd);
      if (jetInd>1)
        FillHist1D("m2jGen",(p4jet[0]+p4jet[1]).mass()); 
    }
}

void JetValidation::beginJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 36 of file JetValidation.cc.

References mix_2012_lumiLevel_15_20_50ns_PoissonOOTPU_cfi::histoFileName, and M_PI.

{
  m_file = new TFile(histoFileName.c_str(),"RECREATE"); 
  
  m_HistNames1D["CaloJetMulti"] = new TH1F("CaloJetMulti","Multiplicity of CaloJets",100,0,100);
  m_HistNames1D["ptCalo"] = new TH1F("ptCalo","p_{T} of CaloJets",7000,0,7000);
  m_HistNames1D["etaCalo"] = new TH1F("etaCalo","#eta of CaloJets",100,-5.0,5.0);
  m_HistNames1D["phiCalo"] = new TH1F("phiCalo","#phi of CaloJets",72,-M_PI, M_PI);
  m_HistNames1D["m2jCalo"] = new TH1F("m2jCalo","Dijet Mass of leading CaloJets",7000,0,14000);
  m_HistNames1D["nTracks"] = new TH1F("nTracks","Number of tracks associated with a jet",100,0,100);
  m_HistNames1D["chargeFraction"] = new TH1F("chargeFraction","Fraction of charged tracks pt",500,0,5);
  m_HistNames1D["emEnergyFraction"] = new TH1F("emEnergyFraction","Jets EM Fraction",110,0,1.1);
  m_HistNames1D["emEnergyInEB"] = new TH1F("emEnergyInEB","Jets emEnergyInEB",7000,0,14000);
  m_HistNames1D["emEnergyInEE"] = new TH1F("emEnergyInEE","Jets emEnergyInEE",7000,0,14000);
  m_HistNames1D["emEnergyInHF"] = new TH1F("emEnergyInHF","Jets emEnergyInHF",7000,0,14000);
  m_HistNames1D["hadEnergyInHB"] = new TH1F("hadEnergyInHB","Jets hadEnergyInHB",7000,0,14000);
  m_HistNames1D["hadEnergyInHE"] = new TH1F("hadEnergyInHE","Jets hadEnergyInHE",7000,0,14000);
  m_HistNames1D["hadEnergyInHF"] = new TH1F("hadEnergyInHF","Jets hadEnergyInHF",7000,0,14000);
  m_HistNames1D["hadEnergyInHO"] = new TH1F("hadEnergyInHO","Jets hadEnergyInHO",7000,0,14000);
  m_HistNamesProfile["EBfractionVsEta"] = new TProfile("EBfractionVsEta","Jets EBfraction vs #eta",100,-5.0,5.0);
  m_HistNamesProfile["EEfractionVsEta"] = new TProfile("EEfractionVsEta","Jets EEfraction vs #eta",100,-5.0,5.0);
  m_HistNamesProfile["HBfractionVsEta"] = new TProfile("HBfractionVsEta","Jets HBfraction vs #eta",100,-5.0,5.0);
  m_HistNamesProfile["HOfractionVsEta"] = new TProfile("HOfractionVsEta","Jets HOfraction vs #eta",100,-5.0,5.0);
  m_HistNamesProfile["HEfractionVsEta"] = new TProfile("HEfractionVsEta","Jets HEfraction vs #eta",100,-5.0,5.0);
  m_HistNamesProfile["HFfractionVsEta"] = new TProfile("HFfractionVsEta","Jets HFfraction vs #eta",100,-5.0,5.0); 
  m_HistNamesProfile["CaloEnergyVsEta"] = new TProfile("CaloEnergyVsEta","CaloJets Energy Vs. Eta",100,-5.0,5.0);
  m_HistNamesProfile["emEnergyVsEta"] = new TProfile("emEnergyVsEta","Jets EM Energy Vs. Eta",100,-5.0,5.0);
  m_HistNamesProfile["hadEnergyVsEta"] = new TProfile("hadEnergyVsEta","Jets HAD Energy Vs. Eta",100,-5.0,5.0);
  if (MCarlo)
    {
      m_HistNames1D["GenJetMulti"] = new TH1F("GenJetMulti","Multiplicity of GenJets",100,0,100);
      m_HistNames1D["ptHat"] = new TH1F("ptHat","p_{T}hat",7000,0,7000);
      m_HistNames1D["ptGen"] = new TH1F("ptGen","p_{T} of GenJets",7000,0,7000);
      m_HistNames1D["etaGen"] = new TH1F("etaGen","#eta of GenJets",100,-5.0,5.0);
      m_HistNames1D["phiGen"] = new TH1F("phiGen","#phi of GenJets",72,-M_PI, M_PI);
      m_HistNames1D["m2jGen"] = new TH1F("m2jGen","Dijet Mass of leading GenJets",7000,0,14000);
      m_HistNames1D["dR"] = new TH1F("dR","GenJets dR with matched CaloJet",200,0,1);
      m_HistNamesProfile["GenEnergyVsEta"] = new TProfile("GenEnergyVsEta","GenJets Energy Vs. Eta",100,-5.0,5.0);
      m_HistNamesProfile["respVsPtBarrel"] = new TProfile("respVsPtBarrel","CaloJet Response of GenJets in Barrel",7000,0,7000);
      m_HistNamesProfile["CaloErespVsEta"] = new TProfile("CaloErespVsEta","Jets Energy Response Vs. Eta",100,-5.0,5.0);
      m_HistNamesProfile["emErespVsEta"] = new TProfile("emErespVsEta","Jets EM Energy Response Vs. Eta",100,-5.0,5.0);
      m_HistNamesProfile["hadErespVsEta"] = new TProfile("hadErespVsEta","Jets HAD Energy Response Vs. Eta",100,-5.0,5.0);
    } 
}

void JetValidation::endJob ( void  )

privatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 201 of file JetValidation.cc.

{
  if (m_file !=0) 
    {
      m_file->cd();
      for (std::map<TString, TH1*>::iterator hid = m_HistNames1D.begin(); hid != m_HistNames1D.end(); hid++)
        hid->second->Write();
      for (std::map<TString, TH2*>::iterator hid = m_HistNames2D.begin(); hid != m_HistNames2D.end(); hid++)
        hid->second->Write();
      for (std::map<TString, TProfile*>::iterator hid = m_HistNamesProfile.begin(); hid != m_HistNamesProfile.end(); hid++)
        hid->second->Write(); 
      delete m_file;
      m_file = 0;      
    }
}

void JetValidation::FillHist1D ( const TString &  histName, const Double_t &  x  )

private

Definition at line 218 of file JetValidation.cc.

References gather_cfg::cout.

{
  std::map<TString, TH1*>::iterator hid=m_HistNames1D.find(histName);
  if (hid==m_HistNames1D.end())
    std::cout << "%fillHist -- Could not find histogram with name: " << histName << std::endl;
  else
    hid->second->Fill(value);
}

void JetValidation::FillHist2D ( const TString &  histName, const Double_t &  x, const Double_t &  y  )

private

Definition at line 227 of file JetValidation.cc.

References gather_cfg::cout.

{
  std::map<TString, TH2*>::iterator hid=m_HistNames2D.find(histName);
  if (hid==m_HistNames2D.end())
    std::cout << "%fillHist -- Could not find histogram with name: " << histName << std::endl;
  else
    hid->second->Fill(valuex,valuey);
}

void JetValidation::FillHistProfile ( const TString &  histName, const Double_t &  x, const Double_t &  y  )

private

Definition at line 236 of file JetValidation.cc.

References gather_cfg::cout, and DEFINE_FWK_MODULE.

{
  std::map<TString, TProfile*>::iterator hid=m_HistNamesProfile.find(histName);
  if (hid==m_HistNamesProfile.end())
    std::cout << "%fillHist -- Could not find histogram with name: " << histName << std::endl;
  else
    hid->second->Fill(valuex,valuey);
} 

Member Data Documentation

std::string JetValidation::calAlgo

private

Definition at line 40 of file JetValidation.h.

double JetValidation::dRmatch

private

Definition at line 35 of file JetValidation.h.

std::string JetValidation::genAlgo

private

Definition at line 39 of file JetValidation.h.

std::string JetValidation::histoFileName

private

Definition at line 38 of file JetValidation.h.

std::string JetValidation::jetTracksAssociator

private

Definition at line 41 of file JetValidation.h.

TFile* JetValidation::m_file

private

Definition at line 32 of file JetValidation.h.

std::map<TString, TH1*> JetValidation::m_HistNames1D

private

Definition at line 29 of file JetValidation.h.

std::map<TString, TH2*> JetValidation::m_HistNames2D

private

Definition at line 30 of file JetValidation.h.

std::map<TString, TProfile*> JetValidation::m_HistNamesProfile

private

Definition at line 31 of file JetValidation.h.

bool JetValidation::MCarlo

private

Definition at line 37 of file JetValidation.h.

int JetValidation::Njets

private

Definition at line 36 of file JetValidation.h.

double JetValidation::PtMin

private

Definition at line 34 of file JetValidation.h.