#include <QcdHighPtDQM.h>

Inheritance diagram for QcdHighPtDQM:

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

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

	QcdHighPtDQM (const edm::ParameterSet &)

	~QcdHighPtDQM () override

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

virtual void	analyze (edm::Event const &, edm::EventSetup const &)

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

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

virtual void	dqmBeginLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

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

	DQMEDAnalyzer ()

virtual void	dqmEndLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &)

virtual void	dqmEndRun (edm::Run const &, edm::EventSetup const &)

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

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

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

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

virtual bool	getCanSaveByLumi ()

Public Member Functions inherited from edm::one::EDProducer< edm::EndRunProducer, edm::one::WatchRuns, edm::EndLuminosityBlockProducer, edm::one::WatchLuminosityBlocks, edm::Accumulator >
	EDProducer ()=default

SerialTaskQueue *	globalLuminosityBlocksQueue () final

SerialTaskQueue *	globalRunsQueue () final

bool	hasAbilityToProduceInBeginLumis () const final

bool	hasAbilityToProduceInBeginRuns () const final

bool	hasAbilityToProduceInEndLumis () const final

bool	hasAbilityToProduceInEndRuns () const final

bool	wantsGlobalLuminosityBlocks () const final

bool	wantsGlobalRuns () const final

Public Member Functions inherited from edm::one::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

bool	wantsStreamLuminosityBlocks () const

bool	wantsStreamRuns () const

	~EDProducerBase () override

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)

	~ProducerBase () noexcept(false) 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

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 const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

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
float	moverl (const reco::CaloMETCollection &metcollection, float &ljpt)

float	movers (const reco::CaloMETCollection &metcollection)

Private Attributes
edm::EDGetTokenT< reco::CaloJetCollection >	jetToken_

std::map< std::string, MonitorElement * >	MEcontainer_

edm::EDGetTokenT< reco::CaloMETCollection >	metToken1_

edm::EDGetTokenT< reco::CaloMETCollection >	metToken2_

edm::EDGetTokenT< reco::CaloMETCollection >	metToken3_

edm::EDGetTokenT< reco::CaloMETCollection >	metToken4_

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

typedef dqm::reco::MonitorElement	MonitorElement

Public Types inherited from edm::one::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const , const char , edm::ProductResolverIndex >>

typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::one::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::ProducerBase
ProducesCollector	producesCollector ()

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 ESProduct , typename ESRecord , Transition Tr = Transition::Event>
auto	esConsumes ()

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

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)

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

DQM Physics Module for High Pt QCD group

Based on DQM/SiPixel and DQM/Physics code Version 1.0, 7/7/09 By Keith Rose

Definition at line 21 of file QcdHighPtDQM.h.

Constructor & Destructor Documentation

QcdHighPtDQM::QcdHighPtDQM ( const edm::ParameterSet & iConfig )

Definition at line 31 of file QcdHighPtDQM.cc.

     : jetToken_(consumes<CaloJetCollection>(iConfig.getUntrackedParameter<edm::InputTag>("jetTag"))),
       metToken1_(consumes<CaloMETCollection>(iConfig.getUntrackedParameter<edm::InputTag>("metTag1"))),
       metToken2_(consumes<CaloMETCollection>(iConfig.getUntrackedParameter<edm::InputTag>("metTag2"))),
       metToken3_(consumes<CaloMETCollection>(iConfig.getUntrackedParameter<edm::InputTag>("metTag3"))),
       metToken4_(consumes<CaloMETCollection>(iConfig.getUntrackedParameter<edm::InputTag>("metTag4"))) {}

QcdHighPtDQM::~QcdHighPtDQM ( )

override

Definition at line 38 of file QcdHighPtDQM.cc.

38 {}

Member Function Documentation

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

override

Definition at line 121 of file QcdHighPtDQM.cc.

References reco::LeafCandidate::energy(), defaults_cfi::etabins, edm::Event::getByToken(), singleTopDQM_cfi::jets, jetToken_, M_PI, MEcontainer_, BTaggingMonitor_cfi::met, metToken1_, metToken2_, metToken3_, metToken4_, movers(), BTaggingMonitoring_cff::njets, particlelevel_cff::njets30, and reco::LeafCandidate::p4().

                                                                         {
   // Get Jets
   edm::Handle<CaloJetCollection> jetHandle;
   iEvent.getByToken(jetToken_, jetHandle);
   const CaloJetCollection &jets = *jetHandle;
   CaloJetCollection::const_iterator jet_iter;
 
   // Get MET collections
   edm::Handle<CaloMETCollection> metHandle;
   iEvent.getByToken(metToken1_, metHandle);
   const CaloMETCollection &met = *metHandle;
 
   edm::Handle<CaloMETCollection> metHOHandle;
   iEvent.getByToken(metToken2_, metHOHandle);
   const CaloMETCollection &metHO = *metHOHandle;
 
   edm::Handle<CaloMETCollection> metNoHFHandle;
   iEvent.getByToken(metToken3_, metNoHFHandle);
   const CaloMETCollection &metNoHF = *metNoHFHandle;
 
   edm::Handle<CaloMETCollection> metNoHFHOHandle;
   iEvent.getByToken(metToken4_, metNoHFHOHandle);
   const CaloMETCollection &metNoHFHO = *metNoHFHOHandle;
 
   // initialize leading jet value and jet multiplicity counter
   int njets = 0;
   int njets30 = 0;
   float leading_jetpt = 0;
   float leading_jeteta = 0;
 
   // initialize variables for picking out leading 2 barrel jets
   reco::CaloJet leadingbarreljet;
   reco::CaloJet secondbarreljet;
   int nbarreljets = 0;
 
   // get bins in eta.
   // Bins correspond to calotower regions.
 
   const float etabins[83] = {
       -5.191, -4.889, -4.716, -4.538, -4.363, -4.191, -4.013, -3.839, -3.664, -3.489, -3.314, -3.139, -2.964, -2.853,
       -2.650, -2.500, -2.322, -2.172, -2.043, -1.930, -1.830, -1.740, -1.653, -1.566, -1.479, -1.392, -1.305, -1.218,
       -1.131, -1.044, -.957,  -.879,  -.783,  -.696,  -.609,  -.522,  -.435,  -.348,  -.261,  -.174,  -.087,  0,
       .087,   .174,   .261,   .348,   .435,   .522,   .609,   .696,   .783,   .879,   .957,   1.044,  1.131,  1.218,
       1.305,  1.392,  1.479,  1.566,  1.653,  1.740,  1.830,  1.930,  2.043,  2.172,  2.322,  2.500,  2.650,  2.853,
       2.964,  3.139,  3.314,  3.489,  3.664,  3.839,  4.013,  4.191,  4.363,  4.538,  4.889,  5.191};
 
   for (jet_iter = jets.begin(); jet_iter != jets.end(); ++jet_iter) {
     njets++;
 
     // get Jet stats
     float jet_pt = jet_iter->pt();
     float jet_eta = jet_iter->eta();
     float jet_phi = jet_iter->phi();
 
     // fill jet Pt and jet EMF
     MEcontainer_["inclusive_jet_pt"]->Fill(jet_pt);
     MEcontainer_["inclusive_jet_EMF"]->Fill(jet_iter->emEnergyFraction());
 
     // pick out up to the first 2 leading barrel jets
     // for use in calculating dijet mass in barrel region
     // also fill jet Pt histogram for barrel
 
     if (jet_eta <= 1.3) {
       MEcontainer_["inclusive_jet_pt_barrel"]->Fill(jet_pt);
       if (nbarreljets == 0) {
         leadingbarreljet = jets[(njets - 1)];
         nbarreljets++;
       } else if (nbarreljets == 1) {
         secondbarreljet = jets[(njets - 1)];
         nbarreljets++;
       }
 
     }
 
     // fill jet Pt for endcap and forward regions
     else if (jet_eta <= 3.0 && jet_eta > 1.3) {
       MEcontainer_["inclusive_jet_pt_endcap"]->Fill(jet_pt);
     } else if (jet_eta <= 5.0 && jet_eta > 3.0) {
       MEcontainer_["inclusive_jet_pt_forward"]->Fill(jet_pt);
     }
 
     // count jet multiplicity for jets with Pt > 30
     if ((jet_pt) > 30)
       njets30++;
 
     // check leading jet quantities
     if (jet_pt > leading_jetpt) {
       leading_jetpt = jet_pt;
       leading_jeteta = jet_eta;
     }
 
     // fill eta-phi plot
     for (int eit = 0; eit < 81; eit++) {
       for (int pit = 0; pit < 72; pit++) {
         float low_eta = etabins[eit];
         float high_eta = etabins[eit + 1];
         float low_phi = (-M_PI) + pit * (M_PI / 36);
         float high_phi = low_phi + (M_PI / 36);
         if (jet_eta > low_eta && jet_eta < high_eta && jet_phi > low_phi && jet_phi < high_phi) {
           MEcontainer_["etaphi"]->Fill((eit - 41), jet_phi);
         }
       }
     }
   }
 
   // after iterating over all jets, fill leading jet quantity histograms
   // and jet multiplicity histograms
 
   MEcontainer_["leading_jet_pt"]->Fill(leading_jetpt);
 
   if (leading_jeteta <= 1.3) {
     MEcontainer_["leading_jet_pt_barrel"]->Fill(leading_jetpt);
   } else if (leading_jeteta <= 3.0 && leading_jeteta > 1.3) {
     MEcontainer_["leading_jet_pt_endcap"]->Fill(leading_jetpt);
   } else if (leading_jeteta <= 5.0 && leading_jeteta > 3.0) {
     MEcontainer_["leading_jet_pt_forward"]->Fill(leading_jetpt);
   }
 
   MEcontainer_["njets"]->Fill(njets);
 
   MEcontainer_["njets30"]->Fill(njets30);
 
   // fill MET over Sum ET and Leading jet PT for all MET flavors
   MEcontainer_["movers_met"]->Fill(movers(met));
   MEcontainer_["moverl_met"]->Fill(movers(met), leading_jetpt);
   MEcontainer_["movers_metho"]->Fill(movers(metHO));
   MEcontainer_["moverl_metho"]->Fill(movers(metHO), leading_jetpt);
   MEcontainer_["movers_metnohf"]->Fill(movers(metNoHF));
   MEcontainer_["moverl_metnohf"]->Fill(movers(metNoHF), leading_jetpt);
   MEcontainer_["movers_metnohfho"]->Fill(movers(metNoHFHO));
   MEcontainer_["moverl_metnohfho"]->Fill(movers(metNoHFHO), leading_jetpt);
 
   // fetch first 3 jet EMF
 
   if (!jets.empty()) {
     MEcontainer_["leading_jet_EMF"]->Fill(jets[0].emEnergyFraction());
     if (jets.size() >= 2) {
       MEcontainer_["second_jet_EMF"]->Fill(jets[1].emEnergyFraction());
       if (jets.size() >= 3) {
         MEcontainer_["third_jet_EMF"]->Fill(jets[2].emEnergyFraction());
       }
     }
   }
 
   // if 2 nontrivial barrel jets, reconstruct dijet mass
 
   if (nbarreljets == 2) {
     if (leadingbarreljet.energy() > 0 && secondbarreljet.energy() > 0) {
       math::XYZTLorentzVector DiJet = leadingbarreljet.p4() + secondbarreljet.p4();
       float dijet_mass = DiJet.mass();
       MEcontainer_["dijet_mass"]->Fill(dijet_mass);
     }
   }
 }

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

overridevirtual

Implements DQMEDAnalyzer.

Definition at line 40 of file QcdHighPtDQM.cc.

References dqm::dqmstoreimpl::DQMStore::IBooker::book1D(), dqm::dqmstoreimpl::DQMStore::IBooker::book2D(), M_PI, MEcontainer_, and dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder().

                                                                                                  {
   iBooker.setCurrentFolder("Physics/QcdHighPt");
 
   MEcontainer_["dijet_mass"] =
       iBooker.book1D("dijet_mass", "dijet resonance invariant mass, barrel region", 100, 0, 1000);
   MEcontainer_["njets"] = iBooker.book1D("njets", "jet multiplicity", 10, 0, 10);
   MEcontainer_["etaphi"] = iBooker.book2D("etaphi", "eta/phi distribution", 83, -42, 42, 72, -M_PI, M_PI);
   MEcontainer_["njets30"] = iBooker.book1D("njets30", "jet multiplicity, pt > 30 GeV", 10, 0, 10);
 
   // book histograms for inclusive jet quantities
   MEcontainer_["inclusive_jet_pt"] = iBooker.book1D("inclusive_jet_pt", "inclusive jet Pt spectrum", 100, 0, 1000);
   MEcontainer_["inclusive_jet_pt_barrel"] =
       iBooker.book1D("inclusive_jet_pt_barrel", "inclusive jet Pt, eta < 1.3", 100, 0, 1000);
   MEcontainer_["inclusive_jet_pt_forward"] =
       iBooker.book1D("inclusive_jet_pt_forward", "inclusive jet Pt, 3.0 < eta < 5.0", 100, 0, 1000);
   MEcontainer_["inclusive_jet_pt_endcap"] =
       iBooker.book1D("inclusive_jet_pt_endcap", "inclusive jet Pt, 1.3 < eta < 3.0", 100, 0, 1000);
 
   // book histograms for leading jet quantities
   MEcontainer_["leading_jet_pt"] = iBooker.book1D("leading_jet_pt", "leading jet Pt", 100, 0, 1000);
   MEcontainer_["leading_jet_pt_barrel"] =
       iBooker.book1D("leading_jet_pt_barrel", "leading jet Pt, eta < 1.3", 100, 0, 1000);
   MEcontainer_["leading_jet_pt_forward"] =
       iBooker.book1D("leading_jet_pt_forward", "leading jet Pt, 3.0 < eta < 5.0", 100, 0, 1000);
   MEcontainer_["leading_jet_pt_endcap"] =
       iBooker.book1D("leading_jet_pt_endcap", "leading jet Pt, 1.3 < eta < 3.0", 100, 0, 1000);
 
   // book histograms for met over sum et and met over leading jet pt for various
   // flavors of MET
   MEcontainer_["movers_met"] = iBooker.book1D("movers_met", "MET over Sum ET for basic MET collection", 50, 0, 1);
   MEcontainer_["moverl_met"] =
       iBooker.book1D("moverl_met", "MET over leading jet Pt for basic MET collection", 50, 0, 2);
 
   MEcontainer_["movers_metho"] = iBooker.book1D("movers_metho", "MET over Sum ET for MET HO collection", 50, 0, 1);
   MEcontainer_["moverl_metho"] =
       iBooker.book1D("moverl_metho", "MET over leading jet Pt for MET HO collection", 50, 0, 2);
 
   MEcontainer_["movers_metnohf"] =
       iBooker.book1D("movers_metnohf", "MET over Sum ET for MET no HF collection", 50, 0, 1);
   MEcontainer_["moverl_metnohf"] =
       iBooker.book1D("moverl_metnohf", "MET over leading jet Pt for MET no HF collection", 50, 0, 2);
 
   MEcontainer_["movers_metnohfho"] =
       iBooker.book1D("movers_metnohfho", "MET over Sum ET for MET no HF HO collection", 50, 0, 1);
   MEcontainer_["moverl_metnohfho"] =
       iBooker.book1D("moverl_metnohfho", "MET over leading jet Pt for MET no HF HO collection", 50, 0, 2);
 
   // book histograms for EMF fraction for all jets and first 3 jets
   MEcontainer_["inclusive_jet_EMF"] = iBooker.book1D("inclusive_jet_EMF", "inclusive jet EMF", 50, -1, 1);
   MEcontainer_["leading_jet_EMF"] = iBooker.book1D("leading_jet_EMF", "leading jet EMF", 50, -1, 1);
   MEcontainer_["second_jet_EMF"] = iBooker.book1D("second_jet_EMF", "second jet EMF", 50, -1, 1);
   MEcontainer_["third_jet_EMF"] = iBooker.book1D("third_jet_EMF", "third jet EMF", 50, -1, 1);
 }

float QcdHighPtDQM::moverl	(	const reco::CaloMETCollection &	metcollection,
		float &	ljpt
	)

private

Definition at line 109 of file QcdHighPtDQM.cc.

References BTaggingMonitor_cfi::met, and mathSSE::sqrt().

                                                                               {
   float metoverl = 0;
   CaloMETCollection::const_iterator met_iter;
   for (met_iter = metcollection.begin(); met_iter != metcollection.end(); ++met_iter) {
     float mex = met_iter->momentum().x();
     float mey = met_iter->momentum().y();
     float met = sqrt(mex * mex + mey * mey);
     metoverl = (met / ljpt);
   }
   return metoverl;
 }

float QcdHighPtDQM::movers ( const reco::CaloMETCollection & metcollection )

private

Definition at line 95 of file QcdHighPtDQM.cc.

References BTaggingMonitor_cfi::met, and mathSSE::sqrt().

Referenced by analyze().

                                                                  {
   float metovers = 0;
   CaloMETCollection::const_iterator met_iter;
   for (met_iter = metcollection.begin(); met_iter != metcollection.end(); ++met_iter) {
     float mex = met_iter->momentum().x();
     float mey = met_iter->momentum().y();
     float met = sqrt(mex * mex + mey * mey);
     float sumet = met_iter->sumEt();
     metovers = (met / sumet);
   }
   return metovers;
 }