#include <DiJetMonitor.h>

Inheritance diagram for DiJetMonitor:

Public Types
typedef dqm::reco::DQMStore	DQMStore

typedef dqm::reco::MonitorElement	MonitorElement

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

Public Types inherited from TriggerDQMBase
typedef dqm::legacy::DQMStore	DQMStore

typedef dqm::legacy::MonitorElement	MonitorElement

Public Member Functions
	DiJetMonitor (const edm::ParameterSet &)

	~DiJetMonitor () override throw ()

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

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)

Public Member Functions inherited from TriggerDQMBase
void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, unsigned nbins, double xmin, double xmax)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const std::vector< double > &binningX)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, unsigned nbinsX, double xmin, double xmax, double ymin, double ymax)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, unsigned nbinsX, double xmin, double xmax, unsigned nbinsY, double ymin, double ymax)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const std::vector< double > &binningX, const std::vector< double > &binningY)

void	setMETitle (ObjME &me, const std::string &titleX, const std::string &titleY)

	TriggerDQMBase ()=default

virtual	~TriggerDQMBase ()=default

Static Public Member Functions
static void	fillDescriptions (edm::ConfigurationDescriptions &descriptions)

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

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Static Public Member Functions inherited from TriggerDQMBase
static void	fillHistoLSPSetDescription (edm::ParameterSetDescription &pset)

static void	fillHistoPSetDescription (edm::ParameterSetDescription &pset)

static MEbinning	getHistoLSPSet (const edm::ParameterSet &pset)

static MEbinning	getHistoPSet (const edm::ParameterSet &pset)

Protected Member Functions
void	analyze (edm::Event const &iEvent, edm::EventSetup const &iSetup) override

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

bool	dijet_selection (double eta_1, double phi_1, double eta_2, double phi_2, double pt_1, double pt_2, int &tag_id, int &probe_id, int Event)

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)

Private Attributes
MEbinning	asy_binning {N_asy, MIN_asy, MAX_asy}

std::unique_ptr< GenericTriggerEventFlag >	den_genTriggerEventFlag_

MEbinning	dijet_eta_binning {DiJet_N_ETA, -DiJet_MAX_ETA, DiJet_MAX_ETA}

const double	DiJet_MAX_ETA = 5

const double	DiJet_MAX_PHI = 3.2

unsigned int	DiJet_N_ETA = 20

unsigned int	DiJet_N_PHI = 32

MEbinning	dijet_phi_binning {DiJet_N_PHI, -DiJet_MAX_PHI, DiJet_MAX_PHI}

MEbinning	dijetpt_binning_

MEbinning	dijetptThr_binning_

edm::EDGetTokenT< reco::PFJetCollection >	dijetSrc_

edm::EDGetTokenT< reco::GsfElectronCollection >	eleToken_

std::string	folderName_

std::string	histoSuffix_

ObjME	jetAsyEtaME_

ObjME	jetEtaPhiME_

ObjME	jetetaPrbME_

ObjME	jetetaTagME_

ObjME	jetphiPrbME_

ObjME	jetpt1ME_

ObjME	jetpt2ME_

ObjME	jetptAsyME_

ObjME	jetptAvgaME_

ObjME	jetptAvgaThrME_

ObjME	jetptAvgbME_

ObjME	jetptPrbME_

ObjME	jetptTagME_

const double	MAX_asy = 1

edm::EDGetTokenT< reco::PFMETCollection >	metToken_

const double	MIN_asy = -1

edm::EDGetTokenT< reco::MuonCollection >	muoToken_

unsigned int	N_asy = 50

int	nmuons_

std::unique_ptr< GenericTriggerEventFlag >	num_genTriggerEventFlag_

double	ptcut_

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

edm::EDPutTokenT< DQMToken >	runToken_

Detailed Description

Definition at line 49 of file DiJetMonitor.h.

Member Typedef Documentation

typedef dqm::reco::DQMStore DiJetMonitor::DQMStore

Definition at line 52 of file DiJetMonitor.h.

typedef dqm::reco::MonitorElement DiJetMonitor::MonitorElement

Definition at line 51 of file DiJetMonitor.h.

Constructor & Destructor Documentation

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

Definition at line 13 of file DiJetMonitor.cc.

References dijetpt_binning_, dijetptThr_binning_, dijetSrc_, folderName_, TriggerDQMBase::getHistoPSet(), edm::ParameterSet::getParameter(), ptcut_, and AlCaHLTBitMon_QueryRunRegistry::string.

     : num_genTriggerEventFlag_(new GenericTriggerEventFlag(
           iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"), consumesCollector(), *this)),
       den_genTriggerEventFlag_(new GenericTriggerEventFlag(
           iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet"), consumesCollector(), *this)) {
   folderName_ = iConfig.getParameter<std::string>("FolderName");
   dijetSrc_ = mayConsume<reco::PFJetCollection>(iConfig.getParameter<edm::InputTag>("dijetSrc"));  //jet
 
   dijetpt_binning_ =
       getHistoPSet(iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("dijetPSet"));
   dijetptThr_binning_ = getHistoPSet(
       iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("dijetPtThrPSet"));
 
   ptcut_ = iConfig.getParameter<double>("ptcut");
 }

DiJetMonitor::~DiJetMonitor	(	)
throw	(
	)

inlineoverride

Definition at line 55 of file DiJetMonitor.h.

References analyze(), bookHistograms(), dijet_selection(), fillDescriptions(), and iEvent.

55 {};

Member Function Documentation

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 158 of file DiJetMonitor.cc.

References den_genTriggerEventFlag_, TriggerDQMBase::ObjME::denominator, dijet_selection(), dijetSrc_, edm::EventID::event(), dqm::impl::MonitorElement::Fill(), edm::Event::getByToken(), edm::EventBase::id(), edm::HandleBase::isValid(), jetAsyEtaME_, jetEtaPhiME_, jetetaPrbME_, jetetaTagME_, jetphiPrbME_, jetpt1ME_, jetpt2ME_, jetptAsyME_, jetptAvgaME_, jetptAvgaThrME_, jetptAvgbME_, jetptPrbME_, jetptTagME_, num_genTriggerEventFlag_, TriggerDQMBase::ObjME::numerator, and ptcut_.

Referenced by ~DiJetMonitor().

                                                                               {
   int Event = -999;
   Event = iEvent.id().event();
   // Filter out events if Trigger Filtering is requested
   if (den_genTriggerEventFlag_->on() && !den_genTriggerEventFlag_->accept(iEvent, iSetup))
     return;
 
   std::vector<double> v_jetpt;
   std::vector<double> v_jeteta;
   std::vector<double> v_jetphi;
 
   v_jetpt.clear();
   v_jeteta.clear();
   v_jetphi.clear();
 
   edm::Handle<reco::PFJetCollection> offjets;
   iEvent.getByToken(dijetSrc_, offjets);
   if (!offjets.isValid()) {
     edm::LogWarning("DiJetMonitor") << "DiJet handle not valid \n";
     return;
   }
   for (reco::PFJetCollection::const_iterator ibegin = offjets->begin(), iend = offjets->end(), ijet = ibegin;
        ijet != iend;
        ++ijet) {
     if (ijet->pt() < ptcut_) {
       continue;
     }
     v_jetpt.push_back(ijet->pt());
     v_jeteta.push_back(ijet->eta());
     v_jetphi.push_back(ijet->phi());
   }
   if (v_jetpt.size() < 2) {
     return;
   }
   double pt_1 = v_jetpt[0];
   double eta_1 = v_jeteta[0];
   double phi_1 = v_jetphi[0];
   double pt_2 = v_jetpt[1];
   double eta_2 = v_jeteta[1];
   double phi_2 = v_jetphi[1];
   double pt_avg_b = (pt_1 + pt_2) * 0.5;
   int tag_id = -999, probe_id = -999;
 
   jetpt1ME_.denominator->Fill(pt_1);
   jetpt2ME_.denominator->Fill(pt_2);
   jetptAvgbME_.denominator->Fill(pt_avg_b);
 
   if (dijet_selection(eta_1, phi_1, eta_2, phi_2, pt_1, pt_2, tag_id, probe_id, Event)) {
     if (tag_id == 0 && probe_id == 1) {
       double pt_asy = (pt_2 - pt_1) / (pt_1 + pt_2);
       double pt_avg = (pt_1 + pt_2) * 0.5;
       jetptAvgaME_.denominator->Fill(pt_avg);
       jetptAvgaThrME_.denominator->Fill(pt_avg);
       jetptTagME_.denominator->Fill(pt_1);
       jetptPrbME_.denominator->Fill(pt_2);
       jetetaPrbME_.denominator->Fill(eta_2);
       jetetaTagME_.denominator->Fill(eta_1);
       jetptAsyME_.denominator->Fill(pt_asy);
       jetphiPrbME_.denominator->Fill(phi_2);
       jetAsyEtaME_.denominator->Fill(pt_asy, eta_2);
       jetEtaPhiME_.denominator->Fill(eta_2, phi_2);
     }
     if (tag_id == 1 && probe_id == 0) {
       double pt_asy = (pt_1 - pt_2) / (pt_2 + pt_1);
       double pt_avg = (pt_2 + pt_1) * 0.5;
       jetptAvgaME_.denominator->Fill(pt_avg);
       jetptAvgaThrME_.denominator->Fill(pt_avg);
       jetptTagME_.denominator->Fill(pt_2);
       jetptPrbME_.denominator->Fill(pt_1);
       jetetaPrbME_.denominator->Fill(eta_1);
       jetetaTagME_.denominator->Fill(eta_2);
       jetptAsyME_.denominator->Fill(pt_asy);
       jetphiPrbME_.denominator->Fill(phi_1);
       jetAsyEtaME_.denominator->Fill(pt_asy, eta_1);
       jetEtaPhiME_.denominator->Fill(eta_1, phi_1);
     }
 
     if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup))
       return;
 
     jetpt1ME_.numerator->Fill(pt_1);
     jetpt2ME_.numerator->Fill(pt_2);
     jetptAvgbME_.numerator->Fill(pt_avg_b);
 
     if (tag_id == 0 && probe_id == 1) {
       double pt_asy = (pt_2 - pt_1) / (pt_1 + pt_2);
       double pt_avg = (pt_1 + pt_2) * 0.5;
       jetptAvgaME_.numerator->Fill(pt_avg);
       jetptAvgaThrME_.numerator->Fill(pt_avg);
       jetptTagME_.numerator->Fill(pt_1);
       jetptPrbME_.numerator->Fill(pt_2);
       jetetaPrbME_.numerator->Fill(eta_2);
       jetetaTagME_.numerator->Fill(eta_1);
       jetptAsyME_.numerator->Fill(pt_asy);
       jetphiPrbME_.numerator->Fill(phi_2);
       jetAsyEtaME_.numerator->Fill(pt_asy, eta_2);
       jetEtaPhiME_.numerator->Fill(eta_2, phi_2);
     }
     if (tag_id == 1 && probe_id == 0) {
       double pt_asy = (pt_1 - pt_2) / (pt_2 + pt_1);
       double pt_avg = (pt_2 + pt_1) * 0.5;
       jetptAvgaME_.numerator->Fill(pt_avg);
       jetptAvgaThrME_.numerator->Fill(pt_avg);
       jetptTagME_.numerator->Fill(pt_2);
       jetptPrbME_.numerator->Fill(pt_1);
       jetetaPrbME_.numerator->Fill(eta_1);
       jetetaTagME_.numerator->Fill(eta_2);
       jetptAsyME_.numerator->Fill(pt_asy);
       jetphiPrbME_.numerator->Fill(phi_1);
       jetAsyEtaME_.numerator->Fill(pt_asy, eta_1);
       jetEtaPhiME_.numerator->Fill(eta_1, phi_1);
     }
   }
 }

void DiJetMonitor::bookHistograms	(	DQMStore::IBooker &	ibooker,
		edm::Run const &	iRun,
		edm::EventSetup const &	iSetup
	)

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 29 of file DiJetMonitor.cc.

References asy_binning, TriggerDQMBase::bookME(), den_genTriggerEventFlag_, dijet_eta_binning, dijet_phi_binning, dijetpt_binning_, dijetptThr_binning_, folderName_, jetAsyEtaME_, jetEtaPhiME_, jetetaPrbME_, jetetaTagME_, jetphiPrbME_, jetpt1ME_, jetpt2ME_, jetptAsyME_, jetptAvgaME_, jetptAvgaThrME_, jetptAvgbME_, jetptPrbME_, jetptTagME_, TriggerDQMBase::MEbinning::nbins, num_genTriggerEventFlag_, dqm::dqmstoreimpl::DQMStore::IBooker::setCurrentFolder(), TriggerDQMBase::setMETitle(), AlCaHLTBitMon_QueryRunRegistry::string, TriggerDQMBase::MEbinning::xmax, and TriggerDQMBase::MEbinning::xmin.

Referenced by ~DiJetMonitor().

                                                                                                            {
   std::string histname, histtitle;
   std::string currentFolder = folderName_;
   ibooker.setCurrentFolder(currentFolder);
 
   histname = "jetpt1";
   histtitle = "leading Jet Pt";
   bookME(ibooker, jetpt1ME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetpt1ME_, "Pt_1 [GeV]", "events");
 
   histname = "jetpt2";
   histtitle = "second leading Jet Pt";
   bookME(ibooker, jetpt2ME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetpt2ME_, "Pt_2 [GeV]", "events");
 
   histname = "jetptAvgB";
   histtitle = "Pt average before offline selection";
   bookME(
       ibooker, jetptAvgbME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetptAvgbME_, "(pt_1 + pt_2)*0.5 [GeV]", "events");
 
   histname = "jetptAvgA";
   histtitle = "Pt average after offline selection";
   bookME(
       ibooker, jetptAvgaME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetptAvgaME_, "(pt_1 + pt_2)*0.5 [GeV]", "events");
 
   histname = "jetptAvgAThr";
   histtitle = "Pt average after offline selection";
   bookME(ibooker,
          jetptAvgaThrME_,
          histname,
          histtitle,
          dijetptThr_binning_.nbins,
          dijetptThr_binning_.xmin,
          dijetptThr_binning_.xmax);
   setMETitle(jetptAvgaThrME_, "(pt_1 + pt_2)*0.5 [GeV]", "events");
 
   histname = "jetptTag";
   histtitle = "Tag Jet Pt";
   bookME(
       ibooker, jetptTagME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetptTagME_, "Pt_tag [GeV]", "events ");
 
   histname = "jetptPrb";
   histtitle = "Probe Jet Pt";
   bookME(
       ibooker, jetptPrbME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
   setMETitle(jetptPrbME_, "Pt_prb [GeV]", "events");
 
   histname = "jetptAsym";
   histtitle = "Jet Pt Asymetry";
   bookME(ibooker, jetptAsyME_, histname, histtitle, asy_binning.nbins, asy_binning.xmin, asy_binning.xmax);
   setMETitle(jetptAsyME_, "(pt_prb - pt_tag)/(pt_prb + pt_tag)", "events");
 
   histname = "jetetaPrb";
   histtitle = "Probe Jet eta";
   bookME(ibooker,
          jetetaPrbME_,
          histname,
          histtitle,
          dijet_eta_binning.nbins,
          dijet_eta_binning.xmin,
          dijet_eta_binning.xmax);
   setMETitle(jetetaPrbME_, "Eta_probe #eta", "events");
 
   histname = "jetetaTag";
   histtitle = "Tag Jet eta";
   bookME(ibooker,
          jetetaTagME_,
          histname,
          histtitle,
          dijet_eta_binning.nbins,
          dijet_eta_binning.xmin,
          dijet_eta_binning.xmax);
   setMETitle(jetetaTagME_, "Eta_tag #eta", "events");
 
   histname = "ptAsymVSetaPrb";
   histtitle = "Pt_Asym vs eta_prb";
   bookME(ibooker,
          jetAsyEtaME_,
          histname,
          histtitle,
          asy_binning.nbins,
          asy_binning.xmin,
          asy_binning.xmax,
          dijet_eta_binning.nbins,
          dijet_eta_binning.xmin,
          dijet_eta_binning.xmax);
   setMETitle(jetAsyEtaME_, "(pt_prb - pt_tag)/(pt_prb + pt_tag)", "Eta_probe #eta");
 
   histname = "etaPrbVSphiPrb";
   histtitle = "eta_prb vs phi_prb";
   bookME(ibooker,
          jetEtaPhiME_,
          histname,
          histtitle,
          dijet_eta_binning.nbins,
          dijet_eta_binning.xmin,
          dijet_eta_binning.xmax,
          dijet_phi_binning.nbins,
          dijet_phi_binning.xmin,
          dijet_phi_binning.xmax);
   setMETitle(jetEtaPhiME_, "Eta_probe #eta", "Phi_probe #phi");
 
   histname = "jetphiPrb";
   histtitle = "Probe Jet phi";
   bookME(ibooker,
          jetphiPrbME_,
          histname,
          histtitle,
          dijet_phi_binning.nbins,
          dijet_phi_binning.xmin,
          dijet_phi_binning.xmax);
   setMETitle(jetphiPrbME_, "Phi_probe #phi", "events");
 
   // Initialize the GenericTriggerEventFlag
   // Initialize the GenericTriggerEventFlag
   if (num_genTriggerEventFlag_ && num_genTriggerEventFlag_->on())
     num_genTriggerEventFlag_->initRun(iRun, iSetup);
   if (den_genTriggerEventFlag_ && den_genTriggerEventFlag_->on())
     den_genTriggerEventFlag_->initRun(iRun, iSetup);
 }

bool DiJetMonitor::dijet_selection	(	double	eta_1,
		double	phi_1,
		double	eta_2,
		double	phi_2,
		double	pt_1,
		double	pt_2,
		int &	tag_id,
		int &	probe_id,
		int	Event
	)

protected

Definition at line 322 of file DiJetMonitor.cc.

References funct::abs(), and DEFINE_FWK_MODULE.

Referenced by analyze(), and ~DiJetMonitor().

                                               {
   double etacut = 1.7;
   double phicut = 2.7;
 
   bool passeta = (std::abs(eta_1) < etacut || std::abs(eta_2) < etacut);  //check that one of the jets in the barrel
 
   float delta_phi_1_2 = (phi_1 - phi_2);
   bool other_cuts = (std::abs(delta_phi_1_2) >= phicut);  //check that jets are back to back
 
   if (std::abs(eta_1) < etacut && std::abs(eta_2) > etacut) {
     tag_id = 0;
     probe_id = 1;
   } else if (std::abs(eta_2) < etacut && std::abs(eta_1) > etacut) {
     tag_id = 1;
     probe_id = 0;
   } else if (std::abs(eta_2) < etacut && std::abs(eta_1) < etacut) {
     int numb = Event % 2;
     if (numb == 0) {
       tag_id = 0;
       probe_id = 1;
     }
     if (numb == 1) {
       tag_id = 1;
       probe_id = 0;
     }
   }
   if (passeta && other_cuts)
     return true;
   else
     return false;
 }

void DiJetMonitor::fillDescriptions ( edm::ConfigurationDescriptions & descriptions )

static

Definition at line 273 of file DiJetMonitor.cc.

References edm::ConfigurationDescriptions::add(), edm::ParameterSetDescription::add(), trigObjTnPSource_cfi::bins, DiJetMonitor_cfi::dijetPSet, DiJetMonitor_cff::dijetPtThrPSet, TriggerDQMBase::fillHistoLSPSetDescription(), TriggerDQMBase::fillHistoPSetDescription(), SiStripSourceConfigTier0_cff::genericTriggerEventPSet, LumiMonitor_cff::histoPSet, HLT_2018_cff::InputTag, LumiMonitor_cff::lsPSet, and AlCaHLTBitMon_QueryRunRegistry::string.

Referenced by ~DiJetMonitor().

                                                                               {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("FolderName", "HLT/JME/Jets/AK4/PF");
 
   desc.add<edm::InputTag>("met", edm::InputTag("pfMet"));
   desc.add<edm::InputTag>("dijetSrc", edm::InputTag("ak4PFJets"));
   desc.add<edm::InputTag>("electrons", edm::InputTag("gedGsfElectrons"));
   desc.add<edm::InputTag>("muons", edm::InputTag("muons"));
   desc.add<int>("njets", 0);
   desc.add<int>("nelectrons", 0);
   desc.add<double>("ptcut", 20);
 
   edm::ParameterSetDescription genericTriggerEventPSet;
   genericTriggerEventPSet.add<bool>("andOr");
   genericTriggerEventPSet.add<edm::InputTag>("dcsInputTag", edm::InputTag("scalersRawToDigi"));
   genericTriggerEventPSet.add<std::vector<int> >("dcsPartitions", {});
   genericTriggerEventPSet.add<bool>("andOrDcs", false);
   genericTriggerEventPSet.add<bool>("errorReplyDcs", true);
   genericTriggerEventPSet.add<std::string>("dbLabel", "");
   genericTriggerEventPSet.add<bool>("andOrHlt", true);
   genericTriggerEventPSet.add<edm::InputTag>("hltInputTag", edm::InputTag("TriggerResults::HLT"));
   genericTriggerEventPSet.add<std::vector<std::string> >("hltPaths", {});
   genericTriggerEventPSet.add<bool>("errorReplyHlt", false);
   genericTriggerEventPSet.add<unsigned int>("verbosityLevel", 1);
 
   desc.add<edm::ParameterSetDescription>("numGenericTriggerEventPSet", genericTriggerEventPSet);
   desc.add<edm::ParameterSetDescription>("denGenericTriggerEventPSet", genericTriggerEventPSet);
 
   edm::ParameterSetDescription histoPSet;
   edm::ParameterSetDescription dijetPSet;
   edm::ParameterSetDescription dijetPtThrPSet;
   fillHistoPSetDescription(dijetPSet);
   fillHistoPSetDescription(dijetPtThrPSet);
   histoPSet.add<edm::ParameterSetDescription>("dijetPSet", dijetPSet);
   histoPSet.add<edm::ParameterSetDescription>("dijetPtThrPSet", dijetPtThrPSet);
   std::vector<double> bins = {
       0.,   20.,  40.,  60.,  80.,  90.,  100., 110., 120., 130., 140., 150., 160.,
       170., 180., 190., 200., 220., 240., 260., 280., 300., 350., 400., 450., 1000.};  // DiJet pT Binning
   histoPSet.add<std::vector<double> >("jetptBinning", bins);
 
   edm::ParameterSetDescription lsPSet;
   fillHistoLSPSetDescription(lsPSet);
   histoPSet.add<edm::ParameterSetDescription>("lsPSet", lsPSet);
   desc.add<edm::ParameterSetDescription>("histoPSet", histoPSet);
 
   descriptions.add("dijetMonitoring", desc);
 }