DiJetMonitor Class Reference

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_, gmtStubs_cfi::eta_1, gmtStubs_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(), gmtStubs_cfi::eta_1, and gmtStubs_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().