#include <GenWeightValidation.h>

Inheritance diagram for GenWeightValidation:

Public Member Functions
void	analyze (const edm::Event &, const edm::EventSetup &) override

void	bookHistograms (DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override

void	dqmBeginRun (const edm::Run &, const edm::EventSetup &) override

	GenWeightValidation (const edm::ParameterSet &)

	~GenWeightValidation () override=default

Public Member Functions inherited from DQMEDAnalyzer
void	accumulate (edm::Event const &event, edm::EventSetup const &setup) final

void	beginLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	beginRun (edm::Run const &run, edm::EventSetup const &setup) final

void	beginStream (edm::StreamID id) final

	DQMEDAnalyzer ()

void	endLuminosityBlock (edm::LuminosityBlock const &lumi, edm::EventSetup const &setup) final

void	endRun (edm::Run const &run, edm::EventSetup const &setup) final

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
	EDProducer ()=default

	EDProducer (const EDProducer &)=delete

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginProcessBlocks () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndProcessBlocks () const final

bool	hasAbilityToProduceInEndRuns () const final

const EDProducer &	operator= (const EDProducer &)=delete

Private Member Functions
void	bookTemplates (DQMHelper &aDqmHelper, std::vector< MonitorElement *> &tmps, const std::string &name, const std::string &title, int nbin, float low, float high, const std::string &xtitle, const std::string &ytitle)

void	fillTemplates (std::vector< MonitorElement *> &tmps, float obs)

Private Attributes
const edm::EDGetTokenT< reco::GenJetCollection >	genJetToken_

const edm::EDGetTokenT< reco::GenParticleCollection >	genParticleToken_

const int	idxFSRdown_

const int	idxFSRup_

const int	idxGenEvtInfo_

const int	idxISRdown_

const int	idxISRup_

int	idxMax_

const double	jetEtaCut_

std::vector< MonitorElement * >	jetMultTemp_

const double	jetPtCut_

const int	jetPtNbin_

const double	jetPtRange_

std::vector< MonitorElement * >	leadJetEtaTemp_

std::vector< MonitorElement * >	leadJetPtTemp_

std::vector< MonitorElement * >	leadLepEtaTemp_

const double	leadLepPtCut_

const int	leadLepPtNbin_

const double	leadLepPtRange_

std::vector< MonitorElement * >	leadLepPtTemp_

const double	lepEtaCut_

MonitorElement *	nEvt_

const int	nJetsNbin_

MonitorElement *	nlogWgt_

const int	rapidityNbin_

const double	rapidityRange_

double	weight_

std::vector< std::vector< double > >	weights_

MonitorElement *	wgtVal_

WeightManager	wmanager_

Additional Inherited Members
Public Types inherited from DQMEDAnalyzer
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >
using	CacheTypes = CacheContexts< T... >

using	GlobalCache = typename CacheTypes::GlobalCache

using	HasAbility = AbilityChecker< T... >

using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache

using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache

using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >

using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache

using	RunCache = typename CacheTypes::RunCache

using	RunContext = RunContextT< RunCache, GlobalCache >

using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Static Public Member Functions inherited from DQMEDAnalyzer
static void	globalEndJob (DQMEDAnalyzerGlobalCache const *)

static void	globalEndLuminosityBlockProduce (edm::LuminosityBlock &lumi, edm::EventSetup const &setup, LuminosityBlockContext const *context)

static void	globalEndRunProduce (edm::Run &run, edm::EventSetup const &setup, RunContext const *context)

static std::unique_ptr< DQMEDAnalyzerGlobalCache >	initializeGlobalCache (edm::ParameterSet const &)

Protected Member Functions inherited from DQMEDAnalyzer
uint64_t	meId () const

Protected Attributes inherited from DQMEDAnalyzer
edm::EDPutTokenT< DQMToken >	lumiToken_

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 26 of file GenWeightValidation.h.

Constructor & Destructor Documentation

◆ GenWeightValidation()

GenWeightValidation::GenWeightValidation ( const edm::ParameterSet & iPSet )

explicit

Definition at line 8 of file GenWeightValidation.cc.

References idxFSRdown_, idxFSRup_, idxISRdown_, idxISRup_, idxMax_, and jetsAK4_CHS_cff::sort.

     : wmanager_(iPSet, consumesCollector()),
       genParticleToken_(consumes<reco::GenParticleCollection>(iPSet.getParameter<edm::InputTag>("genParticles"))),
       genJetToken_(consumes<reco::GenJetCollection>(iPSet.getParameter<edm::InputTag>("genJets"))),
       idxGenEvtInfo_(iPSet.getParameter<int>("whichGenEventInfo")),
       idxFSRup_(iPSet.getParameter<int>("idxFSRup")),
       idxFSRdown_(iPSet.getParameter<int>("idxFSRdown")),
       idxISRup_(iPSet.getParameter<int>("idxISRup")),
       idxISRdown_(iPSet.getParameter<int>("idxISRdown")),
       leadLepPtNbin_(iPSet.getParameter<int>("leadLepPtNbin")),
       rapidityNbin_(iPSet.getParameter<int>("rapidityNbin")),
       leadLepPtRange_(iPSet.getParameter<double>("leadLepPtRange")),
       leadLepPtCut_(iPSet.getParameter<double>("leadLepPtCut")),
       lepEtaCut_(iPSet.getParameter<double>("lepEtaCut")),
       rapidityRange_(iPSet.getParameter<double>("rapidityRange")),
       nJetsNbin_(iPSet.getParameter<int>("nJetsNbin")),
       jetPtNbin_(iPSet.getParameter<int>("jetPtNbin")),
       jetPtCut_(iPSet.getParameter<double>("jetPtCut")),
       jetEtaCut_(iPSet.getParameter<double>("jetEtaCut")),
       jetPtRange_(iPSet.getParameter<double>("jetPtRange")) {
   std::vector<int> idxs = {idxFSRup_, idxFSRdown_, idxISRup_, idxISRdown_};
   std::sort(idxs.begin(), idxs.end(), std::greater<int>());
   idxMax_ = idxs.at(0);
 }

◆ ~GenWeightValidation()

GenWeightValidation::~GenWeightValidation ( )

overridedefault

Member Function Documentation

◆ analyze()

void GenWeightValidation::analyze	(	const edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 122 of file GenWeightValidation.cc.

References funct::abs(), dqm::impl::MonitorElement::Fill(), fillTemplates(), genJetToken_, genParticleToken_, heavyIonCSV_trainingSettings::idx, idxGenEvtInfo_, idxMax_, iEvent, jetEtaCut_, jetMultTemp_, jetPtCut_, leadJetEtaTemp_, leadJetPtTemp_, leadLepEtaTemp_, leadLepPtCut_, leadLepPtTemp_, lepEtaCut_, HLT_2022v15_cff::leptons, nEvt_, nlogWgt_, jetsAK4_CHS_cff::sort, weight_, weights_, WeightManager::weightsCollection(), wgtVal_, and wmanager_.

                                                                                      {
   weights_ = wmanager_.weightsCollection(iEvent);
 
   unsigned weightsSize = weights_.at(idxGenEvtInfo_).size();
   if (weightsSize < 2)
     return;  // no baseline weight in GenEventInfo
 
   weight_ = weights_.at(idxGenEvtInfo_)[0] / std::abs(weights_.at(idxGenEvtInfo_)[0]);
   nEvt_->Fill(0.5, weight_);
   nlogWgt_->Fill(std::log10(weightsSize), weight_);
 
   for (unsigned idx = 0; idx < weightsSize; idx++)
     wgtVal_->Fill(weights_.at(idxGenEvtInfo_)[idx] / weights_.at(idxGenEvtInfo_)[1], weight_);
 
   if ((int)weightsSize <= idxMax_)
     return;  // no PS weights in GenEventInfo
 
   edm::Handle<reco::GenParticleCollection> ptcls;
   iEvent.getByToken(genParticleToken_, ptcls);
   edm::Handle<reco::GenJetCollection> genjets;
   iEvent.getByToken(genJetToken_, genjets);
 
   std::vector<reco::GenParticleRef> leptons;
 
   for (unsigned iptc = 0; iptc < ptcls->size(); iptc++) {
     reco::GenParticleRef ptc(ptcls, iptc);
     if (ptc->status() == 1 && (std::abs(ptc->pdgId()) == 11 || std::abs(ptc->pdgId()) == 13)) {
       if (ptc->pt() > leadLepPtCut_ && std::abs(ptc->eta()) < lepEtaCut_)
         leptons.push_back(ptc);
     }
   }
 
   std::sort(leptons.begin(), leptons.end(), HepMCValidationHelper::sortByPtRef<reco::GenParticleRef>);
 
   if (!leptons.empty()) {
     reco::GenParticleRef leadLep = leptons.at(0);
     fillTemplates(leadLepPtTemp_, leadLep->pt());
     fillTemplates(leadLepEtaTemp_, leadLep->eta());
   }
 
   std::vector<reco::GenJetRef> genjetVec;
 
   for (unsigned igj = 0; igj < genjets->size(); igj++) {
     reco::GenJetRef genjet(genjets, igj);
 
     if (genjet->pt() > jetPtCut_ && std::abs(genjet->eta()) < jetEtaCut_)
       genjetVec.push_back(genjet);
   }
 
   fillTemplates(jetMultTemp_, (float)genjetVec.size());
 
   if (!genjetVec.empty()) {
     std::sort(genjetVec.begin(), genjetVec.end(), HepMCValidationHelper::sortByPtRef<reco::GenJetRef>);
 
     auto leadJet = genjetVec.at(0);
     fillTemplates(leadJetPtTemp_, leadJet->pt());
     fillTemplates(leadJetEtaTemp_, leadJet->eta());
   }
 }  //analyze

◆ bookHistograms()

void GenWeightValidation::bookHistograms	(	DQMStore::IBooker &	iBook,
		edm::Run const &	,
		edm::EventSetup const &
	)

overridevirtual

Setting the DQM top directories

Implements DQMEDAnalyzer.

Definition at line 33 of file GenWeightValidation.cc.

References DQMHelper::book1dHisto(), bookTemplates(), ALCARECODTCalibSynchCosmicsDQM_cff::folderName, jetMultTemp_, jetPtNbin_, jetPtRange_, leadJetEtaTemp_, leadJetPtTemp_, leadLepEtaTemp_, leadLepPtNbin_, leadLepPtRange_, leadLepPtTemp_, nEvt_, nJetsNbin_, nlogWgt_, rapidityNbin_, rapidityRange_, dqm::implementation::NavigatorBase::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, and wgtVal_.

                                                                                                     {
   std::string folderName = "Generator/GenWeight";
   DQMHelper aDqmHelper(&iBook);
   iBook.setCurrentFolder(folderName);
 
   // Number of analyzed events
   nEvt_ = aDqmHelper.book1dHisto("nEvt", "n analyzed Events", 1, 0., 1., "bin", "Number of Events");
   nlogWgt_ = aDqmHelper.book1dHisto("nlogWgt", "Log10(n weights)", 100, 0., 3., "log_{10}(nWgts)", "Number of Events");
   wgtVal_ = aDqmHelper.book1dHisto("wgtVal", "weights", 100, -1.5, 3., "weight", "Number of Weights");
   bookTemplates(aDqmHelper,
                 leadLepPtTemp_,
                 "leadLepPt",
                 "leading lepton Pt",
                 leadLepPtNbin_,
                 0.,
                 leadLepPtRange_,
                 "Pt_{l} (GeV)",
                 "Number of Events");
   bookTemplates(aDqmHelper,
                 leadLepEtaTemp_,
                 "leadLepEta",
                 "leading lepton #eta",
                 rapidityNbin_,
                 -rapidityRange_,
                 rapidityRange_,
                 "#eta_{l}",
                 "Number of Events");
   bookTemplates(aDqmHelper,
                 jetMultTemp_,
                 "JetMultiplicity",
                 "Gen jet multiplicity",
                 nJetsNbin_,
                 0,
                 nJetsNbin_,
                 "n",
                 "Number of Events");
   bookTemplates(aDqmHelper,
                 leadJetPtTemp_,
                 "leadJetPt",
                 "leading Gen jet Pt",
                 jetPtNbin_,
                 0.,
                 jetPtRange_,
                 "Pt_{j} (GeV)",
                 "Number of Events");
   bookTemplates(aDqmHelper,
                 leadJetEtaTemp_,
                 "leadJetEta",
                 "leading Gen jet #eta",
                 rapidityNbin_,
                 -rapidityRange_,
                 rapidityRange_,
                 "#eta_{j}",
                 "Number of Events");
 
   return;
 }

◆ bookTemplates()

void GenWeightValidation::bookTemplates	(	DQMHelper &	aDqmHelper,
		std::vector< MonitorElement *> &	tmps,
		const std::string &	name,
		const std::string &	title,
		int	nbin,
		float	low,
		float	high,
		const std::string &	xtitle,
		const std::string &	ytitle
	)

private

Definition at line 92 of file GenWeightValidation.cc.

References DQMHelper::book1dHisto(), LaserClient_cfi::high, LaserClient_cfi::low, Skims_PA_cff::name, runGCPTkAlMap::title, compareTotals::xtitle, and compareTotals::ytitle.

Referenced by bookHistograms().

                                                                  {
   tmps.push_back(aDqmHelper.book1dHisto(name, title, nbin, low, high, xtitle, ytitle));
   tmps.push_back(aDqmHelper.book1dHisto(name + "FSRup", title + " FSR up", nbin, low, high, xtitle, ytitle));
   tmps.push_back(aDqmHelper.book1dHisto(name + "FSRdn", title + " FSR down", nbin, low, high, xtitle, ytitle));
   tmps.push_back(aDqmHelper.book1dHisto(
       name + "FSRup_ratio", "Ratio of " + title + " FSR up / Nominal", nbin, low, high, xtitle, ytitle));
   tmps.at(3)->setEfficiencyFlag();
   tmps.push_back(aDqmHelper.book1dHisto(
       name + "FSRdn_ratio", "Ratio of " + title + " FSR down / Nominal", nbin, low, high, xtitle, ytitle));
   tmps.at(4)->setEfficiencyFlag();
   tmps.push_back(aDqmHelper.book1dHisto(name + "ISRup", title + " ISR up", nbin, low, high, xtitle, ytitle));
   tmps.push_back(aDqmHelper.book1dHisto(name + "ISRdn", title + " ISR down", nbin, low, high, xtitle, ytitle));
   tmps.push_back(aDqmHelper.book1dHisto(
       name + "ISRup_ratio", "Ratio of " + title + " ISR up / Nominal", nbin, low, high, xtitle, ytitle));
   tmps.at(7)->setEfficiencyFlag();
   tmps.push_back(aDqmHelper.book1dHisto(
       name + "ISRdn_ratio", "Ratio of " + title + " ISR down / Nominal", nbin, low, high, xtitle, ytitle));
   tmps.at(8)->setEfficiencyFlag();
 }  // to get ratio plots correctly - need to modify PostProcessor_cff.py as well!

◆ dqmBeginRun()

void GenWeightValidation::dqmBeginRun	(	const edm::Run &	,
		const edm::EventSetup &
	)

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 120 of file GenWeightValidation.cc.

120 {}

◆ fillTemplates()

void GenWeightValidation::fillTemplates	(	std::vector< MonitorElement *> &	tmps,
		float	obs
	)

private

Definition at line 182 of file GenWeightValidation.cc.

References idxFSRdown_, idxFSRup_, idxGenEvtInfo_, idxISRdown_, idxISRup_, weight_, and weights_.

Referenced by analyze().

                                                                                    {
   tmps.at(0)->Fill(obs, weight_);
   tmps.at(1)->Fill(obs, weights_.at(idxGenEvtInfo_)[idxFSRup_] / weights_.at(idxGenEvtInfo_)[1]);
   tmps.at(2)->Fill(obs, weights_.at(idxGenEvtInfo_)[idxFSRdown_] / weights_.at(idxGenEvtInfo_)[1]);
   tmps.at(5)->Fill(obs, weights_.at(idxGenEvtInfo_)[idxISRup_] / weights_.at(idxGenEvtInfo_)[1]);
   tmps.at(6)->Fill(obs, weights_.at(idxGenEvtInfo_)[idxISRdown_] / weights_.at(idxGenEvtInfo_)[1]);
 }

Member Data Documentation

◆ genJetToken_

const edm::EDGetTokenT<reco::GenJetCollection> GenWeightValidation::genJetToken_

private

Definition at line 60 of file GenWeightValidation.h.

Referenced by analyze().

◆ genParticleToken_

const edm::EDGetTokenT<reco::GenParticleCollection> GenWeightValidation::genParticleToken_

private

Definition at line 59 of file GenWeightValidation.h.

Referenced by analyze().

◆ idxFSRdown_

const int GenWeightValidation::idxFSRdown_

private

Definition at line 62 of file GenWeightValidation.h.

Referenced by fillTemplates(), and GenWeightValidation().

◆ idxFSRup_

const int GenWeightValidation::idxFSRup_

private

Definition at line 62 of file GenWeightValidation.h.

Referenced by fillTemplates(), and GenWeightValidation().

◆ idxGenEvtInfo_

const int GenWeightValidation::idxGenEvtInfo_

private

Definition at line 62 of file GenWeightValidation.h.

Referenced by analyze(), and fillTemplates().

◆ idxISRdown_

const int GenWeightValidation::idxISRdown_

private

Definition at line 62 of file GenWeightValidation.h.

Referenced by fillTemplates(), and GenWeightValidation().

◆ idxISRup_

const int GenWeightValidation::idxISRup_

private

Definition at line 62 of file GenWeightValidation.h.

Referenced by fillTemplates(), and GenWeightValidation().

◆ idxMax_

int GenWeightValidation::idxMax_

private

Definition at line 66 of file GenWeightValidation.h.

Referenced by analyze(), and GenWeightValidation().

◆ jetEtaCut_

const double GenWeightValidation::jetEtaCut_

private

Definition at line 65 of file GenWeightValidation.h.

Referenced by analyze().

◆ jetMultTemp_

std::vector<MonitorElement *> GenWeightValidation::jetMultTemp_

private

Definition at line 55 of file GenWeightValidation.h.

Referenced by analyze(), and bookHistograms().

◆ jetPtCut_

const double GenWeightValidation::jetPtCut_

private

Definition at line 65 of file GenWeightValidation.h.

Referenced by analyze().

◆ jetPtNbin_

const int GenWeightValidation::jetPtNbin_

private

Definition at line 64 of file GenWeightValidation.h.

Referenced by bookHistograms().

◆ jetPtRange_

const double GenWeightValidation::jetPtRange_

private

Definition at line 65 of file GenWeightValidation.h.

Referenced by bookHistograms().

◆ leadJetEtaTemp_

std::vector<MonitorElement *> GenWeightValidation::leadJetEtaTemp_

private

Definition at line 57 of file GenWeightValidation.h.

Referenced by analyze(), and bookHistograms().

◆ leadJetPtTemp_

std::vector<MonitorElement *> GenWeightValidation::leadJetPtTemp_

private

Definition at line 56 of file GenWeightValidation.h.

Referenced by analyze(), and bookHistograms().

◆ leadLepEtaTemp_

std::vector<MonitorElement *> GenWeightValidation::leadLepEtaTemp_

private

Definition at line 54 of file GenWeightValidation.h.

Referenced by analyze(), and bookHistograms().

◆ leadLepPtCut_

const double GenWeightValidation::leadLepPtCut_

private

Definition at line 63 of file GenWeightValidation.h.

Referenced by analyze().

◆ leadLepPtNbin_

const int GenWeightValidation::leadLepPtNbin_

private

Definition at line 62 of file GenWeightValidation.h.

Referenced by bookHistograms().

◆ leadLepPtRange_

const double GenWeightValidation::leadLepPtRange_

private

Definition at line 63 of file GenWeightValidation.h.

Referenced by bookHistograms().

◆ leadLepPtTemp_

std::vector<MonitorElement *> GenWeightValidation::leadLepPtTemp_

private

Definition at line 53 of file GenWeightValidation.h.

Referenced by analyze(), and bookHistograms().

◆ lepEtaCut_

const double GenWeightValidation::lepEtaCut_

private

Definition at line 63 of file GenWeightValidation.h.

Referenced by analyze().

◆ nEvt_

MonitorElement* GenWeightValidation::nEvt_

private

Definition at line 50 of file GenWeightValidation.h.

Referenced by analyze(), and bookHistograms().

◆ nJetsNbin_

const int GenWeightValidation::nJetsNbin_

private

Definition at line 64 of file GenWeightValidation.h.

Referenced by bookHistograms().

◆ nlogWgt_

MonitorElement* GenWeightValidation::nlogWgt_

private

Definition at line 51 of file GenWeightValidation.h.

Referenced by analyze(), and bookHistograms().

◆ rapidityNbin_

const int GenWeightValidation::rapidityNbin_

private

Definition at line 62 of file GenWeightValidation.h.

Referenced by bookHistograms().

◆ rapidityRange_

const double GenWeightValidation::rapidityRange_

private

Definition at line 63 of file GenWeightValidation.h.

Referenced by bookHistograms().

◆ weight_

double GenWeightValidation::weight_

private

Definition at line 47 of file GenWeightValidation.h.

Referenced by analyze(), and fillTemplates().

◆ weights_

std::vector<std::vector<double> > GenWeightValidation::weights_

private

Definition at line 48 of file GenWeightValidation.h.

Referenced by analyze(), and fillTemplates().

◆ wgtVal_

MonitorElement* GenWeightValidation::wgtVal_

private

Definition at line 52 of file GenWeightValidation.h.

Referenced by analyze(), and bookHistograms().

◆ wmanager_

WeightManager GenWeightValidation::wmanager_

private

Definition at line 45 of file GenWeightValidation.h.

Referenced by analyze().

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

◆ GenWeightValidation()

◆ ~GenWeightValidation()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ bookTemplates()

◆ dqmBeginRun()

◆ fillTemplates()

Member Data Documentation

◆ genJetToken_

◆ genParticleToken_

◆ idxFSRdown_

◆ idxFSRup_

◆ idxGenEvtInfo_

◆ idxISRdown_

◆ idxISRup_

◆ idxMax_

◆ jetEtaCut_

◆ jetMultTemp_

◆ jetPtCut_

◆ jetPtNbin_

◆ jetPtRange_

◆ leadJetEtaTemp_

◆ leadJetPtTemp_

◆ leadLepEtaTemp_

◆ leadLepPtCut_

◆ leadLepPtNbin_

◆ leadLepPtRange_

◆ leadLepPtTemp_

◆ lepEtaCut_

◆ nEvt_

◆ nJetsNbin_

◆ nlogWgt_

◆ rapidityNbin_

◆ rapidityRange_

◆ weight_

◆ weights_

◆ wgtVal_

◆ wmanager_