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

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

void	beginStream (edm::StreamID id) final

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

	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

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

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

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const uint nbinsX, const double xmin, const double xmax, const double ymin, const double ymax, const bool bookDen=true)

void	bookME (DQMStore::IBooker &, ObjME &me, const std::string &histname, const std::string &histtitle, const uint nbinsX, const double xmin, const double xmax, const uint nbinsY, const double ymin, const double ymax, const bool bookDen=true)

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

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

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 DQMEDAnalyzer
uint64_t	meId () const

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_

const std::string	folderName_

bool	hltPathsAreValid_

ObjME	jetAsyEtaME_

ObjME	jetEta1ME_

ObjME	jetEta2ME_

ObjME	jetEtaPhiME_

ObjME	jetetaPrbME_

ObjME	jetetaTagME_

ObjME	jetPhi1ME_

ObjME	jetPhi2ME_

ObjME	jetphiPrbME_

ObjME	jetphiTagME_

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_

const bool	requireValidHLTPaths_

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

edm::EDPutTokenT< DQMToken >	runToken_

unsigned int	streamId_

Detailed Description

Definition at line 32 of file DiJetMonitor.cc.

Member Typedef Documentation

◆ DQMStore

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

Definition at line 35 of file DiJetMonitor.cc.

◆ MonitorElement

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

Definition at line 34 of file DiJetMonitor.cc.

Constructor & Destructor Documentation

◆ DiJetMonitor()

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

Definition at line 111 of file DiJetMonitor.cc.

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

◆ ~DiJetMonitor()

DiJetMonitor::~DiJetMonitor	(	)
throw	(
	)

override

Definition at line 126 of file DiJetMonitor.cc.

References den_genTriggerEventFlag_, and num_genTriggerEventFlag_.

                                     {
   if (num_genTriggerEventFlag_) {
     num_genTriggerEventFlag_.reset();
   }
   if (den_genTriggerEventFlag_) {
     den_genTriggerEventFlag_.reset();
   }
 }

Member Function Documentation

◆ analyze()

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

overrideprotectedvirtual

Reimplemented from DQMEDAnalyzer.

Definition at line 306 of file DiJetMonitor.cc.

References den_genTriggerEventFlag_, TriggerDQMBase::ObjME::denominator, dijet_selection(), dijetSrc_, gmt_cfi::eta_1, gmt_cfi::eta_2, dqm::impl::MonitorElement::Fill(), hltPathsAreValid_, iEvent, edm::HandleBase::isValid(), jetAsyEtaME_, jetEta1ME_, jetEta2ME_, jetEtaPhiME_, jetetaPrbME_, jetetaTagME_, jetPhi1ME_, jetPhi2ME_, jetphiPrbME_, jetphiTagME_, jetpt1ME_, jetpt2ME_, jetptAsyME_, jetptAvgaME_, jetptAvgaThrME_, jetptAvgbME_, jetptPrbME_, jetptTagME_, num_genTriggerEventFlag_, TriggerDQMBase::ObjME::numerator, ptcut_, and requireValidHLTPaths_.

                                                                               {
   // if valid HLT paths are required,
   // analyze event only if all paths are valid
   if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
     return;
   }
 
   // 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);
   jetPhi1ME_.denominator->Fill(phi_1);
   jetPhi2ME_.denominator->Fill(phi_2);
   jetEta1ME_.denominator->Fill(eta_1);
   jetEta2ME_.denominator->Fill(eta_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, iEvent.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);
       jetphiTagME_.denominator->Fill(phi_1);
       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);
       jetphiTagME_.denominator->Fill(phi_2);
       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);
     jetPhi1ME_.numerator->Fill(phi_1);
     jetPhi2ME_.numerator->Fill(phi_2);
     jetEta1ME_.numerator->Fill(eta_1);
     jetEta2ME_.numerator->Fill(eta_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);
       jetphiTagME_.numerator->Fill(phi_1);
       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);
       jetphiTagME_.numerator->Fill(phi_2);
       jetAsyEtaME_.numerator->Fill(pt_asy, eta_1);
       jetEtaPhiME_.numerator->Fill(eta_1, phi_1);
     }
   }
 }

◆ bookHistograms()

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

overrideprotectedvirtual

Implements DQMEDAnalyzer.

Definition at line 135 of file DiJetMonitor.cc.

References asy_binning, TriggerDQMBase::bookME(), den_genTriggerEventFlag_, dijet_eta_binning, dijet_phi_binning, dijetpt_binning_, dijetptThr_binning_, folderName_, hltPathsAreValid_, jetAsyEtaME_, jetEta1ME_, jetEta2ME_, jetEtaPhiME_, jetetaPrbME_, jetetaTagME_, jetPhi1ME_, jetPhi2ME_, jetphiPrbME_, jetphiTagME_, jetpt1ME_, jetpt2ME_, jetptAsyME_, jetptAvgaME_, jetptAvgaThrME_, jetptAvgbME_, jetptPrbME_, jetptTagME_, TriggerDQMBase::MEbinning::nbins, num_genTriggerEventFlag_, requireValidHLTPaths_, dqm::implementation::NavigatorBase::setCurrentFolder(), TriggerDQMBase::setMETitle(), AlCaHLTBitMon_QueryRunRegistry::string, TriggerDQMBase::MEbinning::xmax, and TriggerDQMBase::MEbinning::xmin.

                                                                                                            {
   // 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);
   }
 
   // check if every HLT path specified in numerator and denominator has a valid match in the HLT Menu
   hltPathsAreValid_ = (num_genTriggerEventFlag_ && den_genTriggerEventFlag_ && num_genTriggerEventFlag_->on() &&
                        den_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->allHLTPathsAreValid() &&
                        den_genTriggerEventFlag_->allHLTPathsAreValid());
 
   // if valid HLT paths are required,
   // create DQM outputs only if all paths are valid
   if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
     return;
   }
 
   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 = "jetphi1";
   histtitle = "leading Jet Phi";
   bookME(
       ibooker, jetPhi1ME_, histname, histtitle, dijet_phi_binning.nbins, dijet_phi_binning.xmin, dijet_phi_binning.xmax);
   setMETitle(jetPhi1ME_, "Jet_Phi_1", "events");
 
   histname = "jetphi2";
   histtitle = "second leading Jet Phi";
   bookME(
       ibooker, jetPhi2ME_, histname, histtitle, dijet_phi_binning.nbins, dijet_phi_binning.xmin, dijet_phi_binning.xmax);
   setMETitle(jetPhi2ME_, "Jet_Phi_2", "events");
 
   histname = "jeteta1";
   histtitle = "leading Jet Eta";
   bookME(
       ibooker, jetEta1ME_, histname, histtitle, dijet_eta_binning.nbins, dijet_eta_binning.xmin, dijet_eta_binning.xmax);
   setMETitle(jetEta1ME_, "Jet_Eta_1", "events");
 
   histname = "jeteta2";
   histtitle = "second leading Jet Eta";
   bookME(
       ibooker, jetEta2ME_, histname, histtitle, dijet_eta_binning.nbins, dijet_eta_binning.xmin, dijet_eta_binning.xmax);
   setMETitle(jetEta2ME_, "Jet_Eta_2", "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");
 
   histname = "jetphiTag";
   histtitle = "Tag Jet phi";
   bookME(ibooker,
          jetphiTagME_,
          histname,
          histtitle,
          dijet_phi_binning.nbins,
          dijet_phi_binning.xmin,
          dijet_phi_binning.xmax);
   setMETitle(jetphiTagME_, "Phi_tag #phi", "events");
 }

◆ dijet_selection()

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 438 of file DiJetMonitor.cc.

References funct::abs(), gmt_cfi::eta_1, and gmt_cfi::eta_2.

Referenced by analyze().

                                               {
   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;
     }
   }
 
   return (passeta && other_cuts);
 }

◆ fillDescriptions()

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

static

Definition at line 476 of file DiJetMonitor.cc.

References edm::ConfigurationDescriptions::add(), trigObjTnPSource_cfi::bins, submitPVResolutionJobs::desc, DiJetMonitor_cfi::dijetPSet, DiJetMonitor_cff::dijetPtThrPSet, TriggerDQMBase::fillHistoLSPSetDescription(), TriggerDQMBase::fillHistoPSetDescription(), GenericTriggerEventFlag::fillPSetDescription(), SiStripSourceConfigTier0_cff::genericTriggerEventPSet, V0Monitor_cff::histoPSet, ProducerED_cfi::InputTag, V0Monitor_cfi::lsPSet, and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                               {
   edm::ParameterSetDescription desc;
   desc.add<std::string>("FolderName", "HLT/JME/Jets/AK4/PF");
   desc.add<bool>("requireValidHLTPaths", true);
 
   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;
   GenericTriggerEventFlag::fillPSetDescription(genericTriggerEventPSet);
 
   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);
 }

Member Data Documentation

◆ asy_binning

MEbinning DiJetMonitor::asy_binning {N_asy, MIN_asy, MAX_asy}

private

Definition at line 108 of file DiJetMonitor.cc.

Referenced by bookHistograms().

◆ den_genTriggerEventFlag_

std::unique_ptr<GenericTriggerEventFlag> DiJetMonitor::den_genTriggerEventFlag_

private

Definition at line 88 of file DiJetMonitor.cc.

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

◆ dijet_eta_binning

MEbinning DiJetMonitor::dijet_eta_binning {DiJet_N_ETA, -DiJet_MAX_ETA, DiJet_MAX_ETA}

private

Definition at line 102 of file DiJetMonitor.cc.

Referenced by bookHistograms().

◆ DiJet_MAX_ETA

const double DiJetMonitor::DiJet_MAX_ETA = 5

private

Definition at line 99 of file DiJetMonitor.cc.

◆ DiJet_MAX_PHI

const double DiJetMonitor::DiJet_MAX_PHI = 3.2

private

Definition at line 94 of file DiJetMonitor.cc.

◆ DiJet_N_ETA

unsigned int DiJetMonitor::DiJet_N_ETA = 20

private

Definition at line 101 of file DiJetMonitor.cc.

◆ DiJet_N_PHI

unsigned int DiJetMonitor::DiJet_N_PHI = 32

private

Definition at line 96 of file DiJetMonitor.cc.

◆ dijet_phi_binning

MEbinning DiJetMonitor::dijet_phi_binning {DiJet_N_PHI, -DiJet_MAX_PHI, DiJet_MAX_PHI}

private

Definition at line 97 of file DiJetMonitor.cc.

Referenced by bookHistograms().

◆ dijetpt_binning_

MEbinning DiJetMonitor::dijetpt_binning_

private

Definition at line 65 of file DiJetMonitor.cc.

Referenced by bookHistograms().

◆ dijetptThr_binning_

MEbinning DiJetMonitor::dijetptThr_binning_

private

Definition at line 66 of file DiJetMonitor.cc.

Referenced by bookHistograms().

◆ dijetSrc_

edm::EDGetTokenT<reco::PFJetCollection> DiJetMonitor::dijetSrc_

private

Definition at line 63 of file DiJetMonitor.cc.

Referenced by analyze().

◆ eleToken_

edm::EDGetTokenT<reco::GsfElectronCollection> DiJetMonitor::eleToken_

private

Definition at line 61 of file DiJetMonitor.cc.

◆ folderName_

const std::string DiJetMonitor::folderName_

private

Definition at line 55 of file DiJetMonitor.cc.

Referenced by bookHistograms().

◆ hltPathsAreValid_

bool DiJetMonitor::hltPathsAreValid_

private

Definition at line 58 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetAsyEtaME_

ObjME DiJetMonitor::jetAsyEtaME_

private

Definition at line 84 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetEta1ME_

ObjME DiJetMonitor::jetEta1ME_

private

Definition at line 72 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetEta2ME_

ObjME DiJetMonitor::jetEta2ME_

private

Definition at line 73 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetEtaPhiME_

ObjME DiJetMonitor::jetEtaPhiME_

private

Definition at line 85 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetetaPrbME_

ObjME DiJetMonitor::jetetaPrbME_

private

Definition at line 81 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetetaTagME_

ObjME DiJetMonitor::jetetaTagME_

private

Definition at line 82 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetPhi1ME_

ObjME DiJetMonitor::jetPhi1ME_

private

Definition at line 70 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetPhi2ME_

ObjME DiJetMonitor::jetPhi2ME_

private

Definition at line 71 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetphiPrbME_

ObjME DiJetMonitor::jetphiPrbME_

private

Definition at line 83 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetphiTagME_

ObjME DiJetMonitor::jetphiTagME_

private

Definition at line 74 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetpt1ME_

ObjME DiJetMonitor::jetpt1ME_

private

Definition at line 68 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetpt2ME_

ObjME DiJetMonitor::jetpt2ME_

private

Definition at line 69 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptAsyME_

ObjME DiJetMonitor::jetptAsyME_

private

Definition at line 80 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptAvgaME_

ObjME DiJetMonitor::jetptAvgaME_

private

Definition at line 75 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptAvgaThrME_

ObjME DiJetMonitor::jetptAvgaThrME_

private

Definition at line 76 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptAvgbME_

ObjME DiJetMonitor::jetptAvgbME_

private

Definition at line 77 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptPrbME_

ObjME DiJetMonitor::jetptPrbME_

private

Definition at line 79 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptTagME_

ObjME DiJetMonitor::jetptTagME_

private

Definition at line 78 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ MAX_asy

const double DiJetMonitor::MAX_asy = 1

private

Definition at line 104 of file DiJetMonitor.cc.

◆ metToken_

edm::EDGetTokenT<reco::PFMETCollection> DiJetMonitor::metToken_

private

Definition at line 60 of file DiJetMonitor.cc.

◆ MIN_asy

const double DiJetMonitor::MIN_asy = -1

private

Definition at line 105 of file DiJetMonitor.cc.

◆ muoToken_

edm::EDGetTokenT<reco::MuonCollection> DiJetMonitor::muoToken_

private

Definition at line 62 of file DiJetMonitor.cc.

◆ N_asy

unsigned int DiJetMonitor::N_asy = 50

private

Definition at line 107 of file DiJetMonitor.cc.

◆ nmuons_

int DiJetMonitor::nmuons_

private

Definition at line 90 of file DiJetMonitor.cc.

◆ num_genTriggerEventFlag_

std::unique_ptr<GenericTriggerEventFlag> DiJetMonitor::num_genTriggerEventFlag_

private

Definition at line 87 of file DiJetMonitor.cc.

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

◆ ptcut_

double DiJetMonitor::ptcut_

private

Definition at line 91 of file DiJetMonitor.cc.

Referenced by analyze().

◆ requireValidHLTPaths_

const bool DiJetMonitor::requireValidHLTPaths_

private

Definition at line 57 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

Public Types

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Member Typedef Documentation

◆ DQMStore

◆ MonitorElement

Constructor & Destructor Documentation

◆ DiJetMonitor()

◆ ~DiJetMonitor()

Member Function Documentation

◆ analyze()

◆ bookHistograms()

◆ dijet_selection()

◆ fillDescriptions()

Member Data Documentation

◆ asy_binning

◆ den_genTriggerEventFlag_

◆ dijet_eta_binning

◆ DiJet_MAX_ETA

◆ DiJet_MAX_PHI

◆ DiJet_N_ETA

◆ DiJet_N_PHI

◆ dijet_phi_binning

◆ dijetpt_binning_

◆ dijetptThr_binning_

◆ dijetSrc_

◆ eleToken_

◆ folderName_

◆ hltPathsAreValid_

◆ jetAsyEtaME_

◆ jetEta1ME_

◆ jetEta2ME_

◆ jetEtaPhiME_

◆ jetetaPrbME_

◆ jetetaTagME_

◆ jetPhi1ME_

◆ jetPhi2ME_

◆ jetphiPrbME_

◆ jetphiTagME_

◆ jetpt1ME_

◆ jetpt2ME_

◆ jetptAsyME_

◆ jetptAvgaME_

◆ jetptAvgaThrME_

◆ jetptAvgbME_

◆ jetptPrbME_

◆ jetptTagME_

◆ MAX_asy

◆ metToken_

◆ MIN_asy

◆ muoToken_

◆ N_asy

◆ nmuons_

◆ num_genTriggerEventFlag_

◆ ptcut_

◆ requireValidHLTPaths_