#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 ( const edm::ParameterSet & iPSet )

explicit

Definition at line 8 of file GenWeightValidation.cc.

References idxFSRdown_, idxFSRup_, idxISRdown_, idxISRup_, and idxMax_.

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

overridedefault

Member Function Documentation

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_, edm::Event::getByToken(), idxGenEvtInfo_, idxMax_, jetEtaCut_, jetMultTemp_, jetPtCut_, leadJetEtaTemp_, leadJetPtTemp_, leadLepEtaTemp_, leadLepPtCut_, leadLepPtTemp_, lepEtaCut_, HLT_FULL_cff::leptons, nEvt_, nlogWgt_, 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

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

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

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!

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

overridevirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 120 of file GenWeightValidation.cc.

120 {}

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