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

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 &&)=default

	EDConsumerBase (EDConsumerBase const &)=delete

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::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase &	operator= (EDConsumerBase &&)=default

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

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)

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::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	wantsInputProcessBlocks ()

static bool	wantsProcessBlocks ()

static bool	wantsStreamLuminosityBlocks ()

static bool	wantsStreamRuns ()

Protected Member Functions inherited from edm::EDConsumerBase
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

template<BranchType B = InEvent>
EDConsumerBaseAdaptor< B >	consumes (edm::InputTag tag) noexcept

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<Transition Tr = Transition::Event>
constexpr auto	esConsumes () noexcept

template<typename ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag const &tag)

template<Transition Tr = Transition::Event>
auto	esConsumes (ESInputTag tag) noexcept

template<Transition Tr = Transition::Event>
ESGetTokenGeneric	esConsumes (eventsetup::EventSetupRecordKey const &iRecord, eventsetup::DataKey const &iKey)
	Used with EventSetupRecord::doGet. More...

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

Detailed Description

Definition at line 17 of file JetValidation.h.

Constructor & Destructor Documentation

◆ JetValidation()

JetValidation::JetValidation ( edm::ParameterSet const & cfg )

Definition at line 24 of file JetValidation.cc.

                                                        {
   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");
 }

References looper::cfg, mix_2012_lumiLevel_15_20_50ns_PoissonOOTPU_cfi::histoFileName, ECF_cff::Njets, and HLT_FULL_cff::PtMin.

Member Function Documentation

◆ analyze()

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

overrideprivatevirtual

Implements edm::EDAnalyzer.

Definition at line 82 of file JetValidation.cc.

                                                                             {
   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());
   }
 }

References gather_cfg::cout, PbPb_ZMuSkimMuonDPG_cff::deltaR, HGC3DClusterGenMatchSelector_cfi::dR, MillePedeFileConverter_cfg::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, EgHLTOffHistBins_cfi::mass, ECF_cff::Njets, BeamSpotPI::nTracks, DiDispStaMuonMonitor_cfi::pt, reco::LeafCandidate::pt(), HLT_FULL_cff::PtMin, edm::EventID::run(), reco::JetTracksAssociation::tracksNumber(), and reco::JetTracksAssociation::tracksP4().

◆ beginJob()

void JetValidation::beginJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

Definition at line 35 of file JetValidation.cc.

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

References mix_2012_lumiLevel_15_20_50ns_PoissonOOTPU_cfi::histoFileName, and M_PI.

◆ endJob()

void JetValidation::endJob ( void )

overrideprivatevirtual

Reimplemented from edm::EDAnalyzer.

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

◆ FillHist1D()

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

private

Definition at line 210 of file JetValidation.cc.

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

References gather_cfg::cout.

◆ FillHist2D()

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

private

Definition at line 218 of file JetValidation.cc.

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

References gather_cfg::cout.

◆ FillHistProfile()

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

private

Definition at line 226 of file JetValidation.cc.

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

References gather_cfg::cout.

Member Data Documentation

◆ calAlgo

std::string JetValidation::calAlgo

private

Definition at line 40 of file JetValidation.h.

◆ dRmatch

double JetValidation::dRmatch

private

Definition at line 35 of file JetValidation.h.

◆ genAlgo

std::string JetValidation::genAlgo

private

Definition at line 39 of file JetValidation.h.

◆ histoFileName

std::string JetValidation::histoFileName

private

Definition at line 38 of file JetValidation.h.

◆ jetTracksAssociator

std::string JetValidation::jetTracksAssociator

private

Definition at line 41 of file JetValidation.h.

◆ m_file

TFile* JetValidation::m_file

private

Definition at line 32 of file JetValidation.h.

◆ m_HistNames1D

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

private

Definition at line 29 of file JetValidation.h.

◆ m_HistNames2D

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

private

Definition at line 30 of file JetValidation.h.

◆ m_HistNamesProfile

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

private

Definition at line 31 of file JetValidation.h.

◆ MCarlo

bool JetValidation::MCarlo

private

Definition at line 37 of file JetValidation.h.

◆ Njets

int JetValidation::Njets

private

Definition at line 36 of file JetValidation.h.

◆ PtMin

double JetValidation::PtMin

private

Definition at line 34 of file JetValidation.h.

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ JetValidation()

Member Function Documentation

◆ analyze()

◆ beginJob()

◆ endJob()

◆ FillHist1D()

◆ FillHist2D()

◆ FillHistProfile()

Member Data Documentation

◆ calAlgo

◆ dRmatch

◆ genAlgo

◆ histoFileName

◆ jetTracksAssociator

◆ m_file

◆ m_HistNames1D

◆ m_HistNames2D

◆ m_HistNamesProfile

◆ MCarlo

◆ Njets

◆ PtMin