JetValidation Class Reference

#include <JetValidation.h>

Inheritance diagram for JetValidation:

Public Member Functions
	JetValidation (edm::ParameterSet const &cfg)
Public Member Functions inherited from edm::one::EDAnalyzer<>
	EDAnalyzer ()=default
	EDAnalyzer (const EDAnalyzer &)=delete
SerialTaskQueue *	globalLuminosityBlocksQueue () final
SerialTaskQueue *	globalRunsQueue () final
const EDAnalyzer &	operator= (const EDAnalyzer &)=delete
bool	wantsGlobalLuminosityBlocks () const final
bool	wantsGlobalRuns () const final
bool	wantsInputProcessBlocks () const final
bool	wantsProcessBlocks () const final
Public Member Functions inherited from edm::one::EDAnalyzerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
	EDAnalyzerBase ()
ModuleDescription const &	moduleDescription () const
bool	wantsStreamLuminosityBlocks () const
bool	wantsStreamRuns () const
	~EDAnalyzerBase () override
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ESProxyIndex const *	esGetTokenIndices (edm::Transition iTrans) const
std::vector< ESProxyIndex > const &	esGetTokenIndicesVector (edm::Transition iTrans) const
std::vector< ESRecordIndex > const &	esGetTokenRecordIndicesVector (edm::Transition iTrans) const
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector < ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::array< std::vector< ModuleDescription const * > *, NumBranchTypes > &modulesAll, std::vector< ModuleProcessName > &modulesInPreviousProcesses, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	selectInputProcessBlocks (ProductRegistry const &productRegistry, ProcessBlockHelperBase const &processBlockHelperBase)
ProductResolverIndexAndSkipBit	uncheckedIndexFrom (EDGetToken) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
void	updateLookup (eventsetup::ESRecordsToProxyIndices const &)
virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
void	analyze (edm::Event const &e, edm::EventSetup const &iSetup) override
void	beginJob () override
void	endJob () override
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::one::EDAnalyzerBase
typedef EDAnalyzerBase	ModuleType
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::one::EDAnalyzerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()
template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)
template<Transition Tr = Transition::Event>
constexpr auto	esConsumes ()
template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag)
template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
void	resetItemsToGetFrom (BranchType iType)

Detailed Description

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 edm::ParameterSet::getParameter(), diJetCalib::histoFileName, and HLT_FULL_cff::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  )

overrideprivatevirtual

Implements edm::one::EDAnalyzerBase.

Definition at line 82 of file JetValidation.cc.

References gather_cfg::cout, HLT_FULL_cff::deltaR, alignCSCRings::e, reco::CaloJet::emEnergyInEB(), reco::CaloJet::emEnergyInEE(), reco::CaloJet::emEnergyInHF(), reco::LeafCandidate::energy(), PVValHelper::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(), HLT_FULL_cff::jetTracks, ResonanceBuilder::mass, BeamSpotPI::nTracks, DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), HLT_FULL_cff::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->empty())
    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->empty())
      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  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 35 of file JetValidation.cc.

References diJetCalib::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  )

overrideprivatevirtual

Reimplemented from edm::one::EDAnalyzerBase.

Definition at line 194 of file JetValidation.cc.

                           {
  if (m_file != nullptr) {
    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 = nullptr;
  }
}

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

private

Definition at line 210 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 218 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 226 of file JetValidation.cc.

References gather_cfg::cout.

                                                                                                           {
  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.