CMS 3D CMS Logo
List of all members | Public Types | Public Member Functions | Static Public Member Functions | Protected Member Functions | Private Attributes
DiJetMonitor Class Reference
Inheritance diagram for DiJetMonitor:
DQMEDAnalyzer TriggerDQMBase edm::stream::EDProducer< edm::GlobalCache< DQMEDAnalyzerGlobalCache >, edm::EndRunProducer, edm::EndLuminosityBlockProducer, edm::Accumulator >

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 std::vector< double > &binningX, 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 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 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 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 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_
 
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 106 of file DiJetMonitor.cc.

107  : folderName_(iConfig.getParameter<std::string>("FolderName")),
108  requireValidHLTPaths_(iConfig.getParameter<bool>("requireValidHLTPaths")),
109  hltPathsAreValid_(false),
110  dijetSrc_(mayConsume<reco::PFJetCollection>(iConfig.getParameter<edm::InputTag>("dijetSrc"))),
111  dijetpt_binning_(getHistoPSet(
112  iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("dijetPSet"))),
113  dijetptThr_binning_(getHistoPSet(
114  iConfig.getParameter<edm::ParameterSet>("histoPSet").getParameter<edm::ParameterSet>("dijetPtThrPSet"))),
115  num_genTriggerEventFlag_(new GenericTriggerEventFlag(
116  iConfig.getParameter<edm::ParameterSet>("numGenericTriggerEventPSet"), consumesCollector(), *this)),
117  den_genTriggerEventFlag_(new GenericTriggerEventFlag(
118  iConfig.getParameter<edm::ParameterSet>("denGenericTriggerEventPSet"), consumesCollector(), *this)),
119  ptcut_(iConfig.getParameter<double>("ptcut")) {}

◆ ~DiJetMonitor()

DiJetMonitor::~DiJetMonitor ( )
throw (
)
override

Definition at line 121 of file DiJetMonitor.cc.

121  {
122  if (num_genTriggerEventFlag_) {
123  num_genTriggerEventFlag_.reset();
124  }
125  if (den_genTriggerEventFlag_) {
126  den_genTriggerEventFlag_.reset();
127  }
128 }

References den_genTriggerEventFlag_, and num_genTriggerEventFlag_.

Member Function Documentation

◆ analyze()

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

Reimplemented from DQMEDAnalyzer.

Definition at line 266 of file DiJetMonitor.cc.

266  {
267  // if valid HLT paths are required,
268  // analyze event only if all paths are valid
269  if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
270  return;
271  }
272 
273  // Filter out events if Trigger Filtering is requested
274  if (den_genTriggerEventFlag_->on() && !den_genTriggerEventFlag_->accept(iEvent, iSetup))
275  return;
276 
277  std::vector<double> v_jetpt;
278  std::vector<double> v_jeteta;
279  std::vector<double> v_jetphi;
280 
281  v_jetpt.clear();
282  v_jeteta.clear();
283  v_jetphi.clear();
284 
285  edm::Handle<reco::PFJetCollection> offjets;
286  iEvent.getByToken(dijetSrc_, offjets);
287  if (!offjets.isValid()) {
288  edm::LogWarning("DiJetMonitor") << "DiJet handle not valid \n";
289  return;
290  }
291  for (reco::PFJetCollection::const_iterator ibegin = offjets->begin(), iend = offjets->end(), ijet = ibegin;
292  ijet != iend;
293  ++ijet) {
294  if (ijet->pt() < ptcut_) {
295  continue;
296  }
297  v_jetpt.push_back(ijet->pt());
298  v_jeteta.push_back(ijet->eta());
299  v_jetphi.push_back(ijet->phi());
300  }
301  if (v_jetpt.size() < 2) {
302  return;
303  }
304  double pt_1 = v_jetpt[0];
305  double eta_1 = v_jeteta[0];
306  double phi_1 = v_jetphi[0];
307  double pt_2 = v_jetpt[1];
308  double eta_2 = v_jeteta[1];
309  double phi_2 = v_jetphi[1];
310  double pt_avg_b = (pt_1 + pt_2) * 0.5;
311  int tag_id = -999, probe_id = -999;
312 
313  jetpt1ME_.denominator->Fill(pt_1);
314  jetpt2ME_.denominator->Fill(pt_2);
315  jetptAvgbME_.denominator->Fill(pt_avg_b);
316 
317  if (dijet_selection(eta_1, phi_1, eta_2, phi_2, pt_1, pt_2, tag_id, probe_id, iEvent.id().event())) {
318  if (tag_id == 0 && probe_id == 1) {
319  double pt_asy = (pt_2 - pt_1) / (pt_1 + pt_2);
320  double pt_avg = (pt_1 + pt_2) * 0.5;
321  jetptAvgaME_.denominator->Fill(pt_avg);
322  jetptAvgaThrME_.denominator->Fill(pt_avg);
323  jetptTagME_.denominator->Fill(pt_1);
324  jetptPrbME_.denominator->Fill(pt_2);
325  jetetaPrbME_.denominator->Fill(eta_2);
326  jetetaTagME_.denominator->Fill(eta_1);
327  jetptAsyME_.denominator->Fill(pt_asy);
328  jetphiPrbME_.denominator->Fill(phi_2);
329  jetAsyEtaME_.denominator->Fill(pt_asy, eta_2);
330  jetEtaPhiME_.denominator->Fill(eta_2, phi_2);
331  }
332  if (tag_id == 1 && probe_id == 0) {
333  double pt_asy = (pt_1 - pt_2) / (pt_2 + pt_1);
334  double pt_avg = (pt_2 + pt_1) * 0.5;
335  jetptAvgaME_.denominator->Fill(pt_avg);
336  jetptAvgaThrME_.denominator->Fill(pt_avg);
337  jetptTagME_.denominator->Fill(pt_2);
338  jetptPrbME_.denominator->Fill(pt_1);
339  jetetaPrbME_.denominator->Fill(eta_1);
340  jetetaTagME_.denominator->Fill(eta_2);
341  jetptAsyME_.denominator->Fill(pt_asy);
342  jetphiPrbME_.denominator->Fill(phi_1);
343  jetAsyEtaME_.denominator->Fill(pt_asy, eta_1);
344  jetEtaPhiME_.denominator->Fill(eta_1, phi_1);
345  }
346 
347  if (num_genTriggerEventFlag_->on() && !num_genTriggerEventFlag_->accept(iEvent, iSetup))
348  return;
349 
350  jetpt1ME_.numerator->Fill(pt_1);
351  jetpt2ME_.numerator->Fill(pt_2);
352  jetptAvgbME_.numerator->Fill(pt_avg_b);
353 
354  if (tag_id == 0 && probe_id == 1) {
355  double pt_asy = (pt_2 - pt_1) / (pt_1 + pt_2);
356  double pt_avg = (pt_1 + pt_2) * 0.5;
357  jetptAvgaME_.numerator->Fill(pt_avg);
358  jetptAvgaThrME_.numerator->Fill(pt_avg);
359  jetptTagME_.numerator->Fill(pt_1);
360  jetptPrbME_.numerator->Fill(pt_2);
361  jetetaPrbME_.numerator->Fill(eta_2);
362  jetetaTagME_.numerator->Fill(eta_1);
363  jetptAsyME_.numerator->Fill(pt_asy);
364  jetphiPrbME_.numerator->Fill(phi_2);
365  jetAsyEtaME_.numerator->Fill(pt_asy, eta_2);
366  jetEtaPhiME_.numerator->Fill(eta_2, phi_2);
367  }
368  if (tag_id == 1 && probe_id == 0) {
369  double pt_asy = (pt_1 - pt_2) / (pt_2 + pt_1);
370  double pt_avg = (pt_2 + pt_1) * 0.5;
371  jetptAvgaME_.numerator->Fill(pt_avg);
372  jetptAvgaThrME_.numerator->Fill(pt_avg);
373  jetptTagME_.numerator->Fill(pt_2);
374  jetptPrbME_.numerator->Fill(pt_1);
375  jetetaPrbME_.numerator->Fill(eta_1);
376  jetetaTagME_.numerator->Fill(eta_2);
377  jetptAsyME_.numerator->Fill(pt_asy);
378  jetphiPrbME_.numerator->Fill(phi_1);
379  jetAsyEtaME_.numerator->Fill(pt_asy, eta_1);
380  jetEtaPhiME_.numerator->Fill(eta_1, phi_1);
381  }
382  }
383 }

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

◆ bookHistograms()

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

Implements DQMEDAnalyzer.

Definition at line 130 of file DiJetMonitor.cc.

130  {
131  // Initialize the GenericTriggerEventFlag
132  if (num_genTriggerEventFlag_ && num_genTriggerEventFlag_->on()) {
133  num_genTriggerEventFlag_->initRun(iRun, iSetup);
134  }
135  if (den_genTriggerEventFlag_ && den_genTriggerEventFlag_->on()) {
136  den_genTriggerEventFlag_->initRun(iRun, iSetup);
137  }
138 
139  // check if every HLT path specified in numerator and denominator has a valid match in the HLT Menu
140  hltPathsAreValid_ = (num_genTriggerEventFlag_ && den_genTriggerEventFlag_ && num_genTriggerEventFlag_->on() &&
141  den_genTriggerEventFlag_->on() && num_genTriggerEventFlag_->allHLTPathsAreValid() &&
142  den_genTriggerEventFlag_->allHLTPathsAreValid());
143 
144  // if valid HLT paths are required,
145  // create DQM outputs only if all paths are valid
146  if (requireValidHLTPaths_ and (not hltPathsAreValid_)) {
147  return;
148  }
149 
150  std::string histname, histtitle;
151  std::string currentFolder = folderName_;
152  ibooker.setCurrentFolder(currentFolder);
153 
154  histname = "jetpt1";
155  histtitle = "leading Jet Pt";
156  bookME(ibooker, jetpt1ME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
157  setMETitle(jetpt1ME_, "Pt_1 [GeV]", "events");
158 
159  histname = "jetpt2";
160  histtitle = "second leading Jet Pt";
161  bookME(ibooker, jetpt2ME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
162  setMETitle(jetpt2ME_, "Pt_2 [GeV]", "events");
163 
164  histname = "jetptAvgB";
165  histtitle = "Pt average before offline selection";
166  bookME(
167  ibooker, jetptAvgbME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
168  setMETitle(jetptAvgbME_, "(pt_1 + pt_2)*0.5 [GeV]", "events");
169 
170  histname = "jetptAvgA";
171  histtitle = "Pt average after offline selection";
172  bookME(
173  ibooker, jetptAvgaME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
174  setMETitle(jetptAvgaME_, "(pt_1 + pt_2)*0.5 [GeV]", "events");
175 
176  histname = "jetptAvgAThr";
177  histtitle = "Pt average after offline selection";
178  bookME(ibooker,
179  jetptAvgaThrME_,
180  histname,
181  histtitle,
182  dijetptThr_binning_.nbins,
183  dijetptThr_binning_.xmin,
184  dijetptThr_binning_.xmax);
185  setMETitle(jetptAvgaThrME_, "(pt_1 + pt_2)*0.5 [GeV]", "events");
186 
187  histname = "jetptTag";
188  histtitle = "Tag Jet Pt";
189  bookME(
190  ibooker, jetptTagME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
191  setMETitle(jetptTagME_, "Pt_tag [GeV]", "events ");
192 
193  histname = "jetptPrb";
194  histtitle = "Probe Jet Pt";
195  bookME(
196  ibooker, jetptPrbME_, histname, histtitle, dijetpt_binning_.nbins, dijetpt_binning_.xmin, dijetpt_binning_.xmax);
197  setMETitle(jetptPrbME_, "Pt_prb [GeV]", "events");
198 
199  histname = "jetptAsym";
200  histtitle = "Jet Pt Asymetry";
201  bookME(ibooker, jetptAsyME_, histname, histtitle, asy_binning.nbins, asy_binning.xmin, asy_binning.xmax);
202  setMETitle(jetptAsyME_, "(pt_prb - pt_tag)/(pt_prb + pt_tag)", "events");
203 
204  histname = "jetetaPrb";
205  histtitle = "Probe Jet eta";
206  bookME(ibooker,
207  jetetaPrbME_,
208  histname,
209  histtitle,
210  dijet_eta_binning.nbins,
211  dijet_eta_binning.xmin,
212  dijet_eta_binning.xmax);
213  setMETitle(jetetaPrbME_, "Eta_probe #eta", "events");
214 
215  histname = "jetetaTag";
216  histtitle = "Tag Jet eta";
217  bookME(ibooker,
218  jetetaTagME_,
219  histname,
220  histtitle,
221  dijet_eta_binning.nbins,
222  dijet_eta_binning.xmin,
223  dijet_eta_binning.xmax);
224  setMETitle(jetetaTagME_, "Eta_tag #eta", "events");
225 
226  histname = "ptAsymVSetaPrb";
227  histtitle = "Pt_Asym vs eta_prb";
228  bookME(ibooker,
229  jetAsyEtaME_,
230  histname,
231  histtitle,
232  asy_binning.nbins,
233  asy_binning.xmin,
234  asy_binning.xmax,
235  dijet_eta_binning.nbins,
236  dijet_eta_binning.xmin,
237  dijet_eta_binning.xmax);
238  setMETitle(jetAsyEtaME_, "(pt_prb - pt_tag)/(pt_prb + pt_tag)", "Eta_probe #eta");
239 
240  histname = "etaPrbVSphiPrb";
241  histtitle = "eta_prb vs phi_prb";
242  bookME(ibooker,
243  jetEtaPhiME_,
244  histname,
245  histtitle,
246  dijet_eta_binning.nbins,
247  dijet_eta_binning.xmin,
248  dijet_eta_binning.xmax,
249  dijet_phi_binning.nbins,
250  dijet_phi_binning.xmin,
251  dijet_phi_binning.xmax);
252  setMETitle(jetEtaPhiME_, "Eta_probe #eta", "Phi_probe #phi");
253 
254  histname = "jetphiPrb";
255  histtitle = "Probe Jet phi";
256  bookME(ibooker,
257  jetphiPrbME_,
258  histname,
259  histtitle,
260  dijet_phi_binning.nbins,
261  dijet_phi_binning.xmin,
262  dijet_phi_binning.xmax);
263  setMETitle(jetphiPrbME_, "Phi_probe #phi", "events");
264 }

References asy_binning, TriggerDQMBase::bookME(), den_genTriggerEventFlag_, dijet_eta_binning, dijet_phi_binning, dijetpt_binning_, dijetptThr_binning_, folderName_, hltPathsAreValid_, jetAsyEtaME_, jetEtaPhiME_, jetetaPrbME_, jetetaTagME_, jetphiPrbME_, 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.

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

394  {
395  double etacut = 1.7;
396  double phicut = 2.7;
397 
398  bool passeta = (std::abs(eta_1) < etacut || std::abs(eta_2) < etacut); //check that one of the jets in the barrel
399 
400  float delta_phi_1_2 = (phi_1 - phi_2);
401  bool other_cuts = (std::abs(delta_phi_1_2) >= phicut); //check that jets are back to back
402 
403  if (std::abs(eta_1) < etacut && std::abs(eta_2) > etacut) {
404  tag_id = 0;
405  probe_id = 1;
406  } else if (std::abs(eta_2) < etacut && std::abs(eta_1) > etacut) {
407  tag_id = 1;
408  probe_id = 0;
409  } else if (std::abs(eta_2) < etacut && std::abs(eta_1) < etacut) {
410  int numb = Event % 2;
411  if (numb == 0) {
412  tag_id = 0;
413  probe_id = 1;
414  }
415  if (numb == 1) {
416  tag_id = 1;
417  probe_id = 0;
418  }
419  }
420 
421  return (passeta && other_cuts);
422 }

References funct::abs().

Referenced by analyze().

◆ fillDescriptions()

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

Definition at line 424 of file DiJetMonitor.cc.

424  {
425  edm::ParameterSetDescription desc;
426  desc.add<std::string>("FolderName", "HLT/JME/Jets/AK4/PF");
427  desc.add<bool>("requireValidHLTPaths", true);
428 
429  desc.add<edm::InputTag>("met", edm::InputTag("pfMet"));
430  desc.add<edm::InputTag>("dijetSrc", edm::InputTag("ak4PFJets"));
431  desc.add<edm::InputTag>("electrons", edm::InputTag("gedGsfElectrons"));
432  desc.add<edm::InputTag>("muons", edm::InputTag("muons"));
433  desc.add<int>("njets", 0);
434  desc.add<int>("nelectrons", 0);
435  desc.add<double>("ptcut", 20);
436 
437  edm::ParameterSetDescription genericTriggerEventPSet;
438  genericTriggerEventPSet.add<bool>("andOr");
439  genericTriggerEventPSet.add<edm::InputTag>("dcsInputTag", edm::InputTag("scalersRawToDigi"));
440  genericTriggerEventPSet.add<std::vector<int> >("dcsPartitions", {});
441  genericTriggerEventPSet.add<bool>("andOrDcs", false);
442  genericTriggerEventPSet.add<bool>("errorReplyDcs", true);
443  genericTriggerEventPSet.add<std::string>("dbLabel", "");
444  genericTriggerEventPSet.add<bool>("andOrHlt", true);
445  genericTriggerEventPSet.add<edm::InputTag>("hltInputTag", edm::InputTag("TriggerResults::HLT"));
446  genericTriggerEventPSet.add<std::vector<std::string> >("hltPaths", {});
447  genericTriggerEventPSet.add<bool>("errorReplyHlt", false);
448  genericTriggerEventPSet.add<unsigned int>("verbosityLevel", 1);
449 
450  desc.add<edm::ParameterSetDescription>("numGenericTriggerEventPSet", genericTriggerEventPSet);
451  desc.add<edm::ParameterSetDescription>("denGenericTriggerEventPSet", genericTriggerEventPSet);
452 
453  edm::ParameterSetDescription histoPSet;
454  edm::ParameterSetDescription dijetPSet;
455  edm::ParameterSetDescription dijetPtThrPSet;
456  fillHistoPSetDescription(dijetPSet);
457  fillHistoPSetDescription(dijetPtThrPSet);
458  histoPSet.add<edm::ParameterSetDescription>("dijetPSet", dijetPSet);
459  histoPSet.add<edm::ParameterSetDescription>("dijetPtThrPSet", dijetPtThrPSet);
460  std::vector<double> bins = {
461  0., 20., 40., 60., 80., 90., 100., 110., 120., 130., 140., 150., 160.,
462  170., 180., 190., 200., 220., 240., 260., 280., 300., 350., 400., 450., 1000.}; // DiJet pT Binning
463  histoPSet.add<std::vector<double> >("jetptBinning", bins);
464 
465  edm::ParameterSetDescription lsPSet;
466  fillHistoLSPSetDescription(lsPSet);
467  histoPSet.add<edm::ParameterSetDescription>("lsPSet", lsPSet);
468  desc.add<edm::ParameterSetDescription>("histoPSet", histoPSet);
469 
470  descriptions.add("dijetMonitoring", desc);
471 }

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

Member Data Documentation

◆ asy_binning

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

Definition at line 103 of file DiJetMonitor.cc.

Referenced by bookHistograms().

◆ den_genTriggerEventFlag_

std::unique_ptr<GenericTriggerEventFlag> DiJetMonitor::den_genTriggerEventFlag_
private

Definition at line 83 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 97 of file DiJetMonitor.cc.

Referenced by bookHistograms().

◆ DiJet_MAX_ETA

const double DiJetMonitor::DiJet_MAX_ETA = 5
private

Definition at line 94 of file DiJetMonitor.cc.

◆ DiJet_MAX_PHI

const double DiJetMonitor::DiJet_MAX_PHI = 3.2
private

Definition at line 89 of file DiJetMonitor.cc.

◆ DiJet_N_ETA

unsigned int DiJetMonitor::DiJet_N_ETA = 20
private

Definition at line 96 of file DiJetMonitor.cc.

◆ DiJet_N_PHI

unsigned int DiJetMonitor::DiJet_N_PHI = 32
private

Definition at line 91 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 92 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 79 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetEtaPhiME_

ObjME DiJetMonitor::jetEtaPhiME_
private

Definition at line 80 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetetaPrbME_

ObjME DiJetMonitor::jetetaPrbME_
private

Definition at line 76 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetetaTagME_

ObjME DiJetMonitor::jetetaTagME_
private

Definition at line 77 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetphiPrbME_

ObjME DiJetMonitor::jetphiPrbME_
private

Definition at line 78 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 75 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptAvgaME_

ObjME DiJetMonitor::jetptAvgaME_
private

Definition at line 70 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptAvgaThrME_

ObjME DiJetMonitor::jetptAvgaThrME_
private

Definition at line 71 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptAvgbME_

ObjME DiJetMonitor::jetptAvgbME_
private

Definition at line 72 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptPrbME_

ObjME DiJetMonitor::jetptPrbME_
private

Definition at line 74 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ jetptTagME_

ObjME DiJetMonitor::jetptTagME_
private

Definition at line 73 of file DiJetMonitor.cc.

Referenced by analyze(), and bookHistograms().

◆ MAX_asy

const double DiJetMonitor::MAX_asy = 1
private

Definition at line 99 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 100 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 102 of file DiJetMonitor.cc.

◆ nmuons_

int DiJetMonitor::nmuons_
private

Definition at line 85 of file DiJetMonitor.cc.

◆ num_genTriggerEventFlag_

std::unique_ptr<GenericTriggerEventFlag> DiJetMonitor::num_genTriggerEventFlag_
private

Definition at line 82 of file DiJetMonitor.cc.

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

◆ ptcut_

double DiJetMonitor::ptcut_
private

Definition at line 86 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().

DiJetMonitor::jetptAvgbME_
ObjME jetptAvgbME_
Definition: DiJetMonitor.cc:72
DiJetMonitor::jetEtaPhiME_
ObjME jetEtaPhiME_
Definition: DiJetMonitor.cc:80
HLT_FULL_cff.InputTag
InputTag
Definition: HLT_FULL_cff.py:89301
edm::ParameterSetDescription
Definition: ParameterSetDescription.h:52
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:32
DiJetMonitor::folderName_
const std::string folderName_
Definition: DiJetMonitor.cc:55
SiStripSourceConfigTier0_cff.genericTriggerEventPSet
genericTriggerEventPSet
Definition: SiStripSourceConfigTier0_cff.py:96
edm::Handle< reco::PFJetCollection >
DiJetMonitor::jetetaPrbME_
ObjME jetetaPrbME_
Definition: DiJetMonitor.cc:76
edm::LogWarning
Log< level::Warning, false > LogWarning
Definition: MessageLogger.h:122
DiJetMonitor::jetpt2ME_
ObjME jetpt2ME_
Definition: DiJetMonitor.cc:69
TriggerDQMBase::ObjME::numerator
MonitorElement * numerator
Definition: TriggerDQMBase.h:27
DiJetMonitor::jetptAvgaME_
ObjME jetptAvgaME_
Definition: DiJetMonitor.cc:70
LumiMonitor_cff.lsPSet
lsPSet
Definition: LumiMonitor_cff.py:33
edm::ConfigurationDescriptions::add
void add(std::string const &label, ParameterSetDescription const &psetDescription)
Definition: ConfigurationDescriptions.cc:57
TriggerDQMBase::setMETitle
void setMETitle(ObjME &me, const std::string &titleX, const std::string &titleY)
Definition: TriggerDQMBase.cc:3
DiJetMonitor::jetpt1ME_
ObjME jetpt1ME_
Definition: DiJetMonitor.cc:68
DiJetMonitor::jetptAvgaThrME_
ObjME jetptAvgaThrME_
Definition: DiJetMonitor.cc:71
TriggerDQMBase::MEbinning::xmax
double xmax
Definition: TriggerDQMBase.h:19
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
Event
TriggerDQMBase::bookME
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)
Definition: TriggerDQMBase.cc:15
DiJetMonitor::jetptTagME_
ObjME jetptTagME_
Definition: DiJetMonitor.cc:73
DiJetMonitor::dijetSrc_
edm::EDGetTokenT< reco::PFJetCollection > dijetSrc_
Definition: DiJetMonitor.cc:63
DiJetMonitor::jetetaTagME_
ObjME jetetaTagME_
Definition: DiJetMonitor.cc:77
TriggerDQMBase::fillHistoPSetDescription
static void fillHistoPSetDescription(edm::ParameterSetDescription &pset)
Definition: TriggerDQMBase.cc:108
DiJetMonitor::dijet_phi_binning
MEbinning dijet_phi_binning
Definition: DiJetMonitor.cc:92
DiJetMonitor::jetphiPrbME_
ObjME jetphiPrbME_
Definition: DiJetMonitor.cc:78
DiJetMonitor_cfi.dijetPSet
dijetPSet
Definition: DiJetMonitor_cfi.py:6
edm::ParameterSet
Definition: ParameterSet.h:47
GenericTriggerEventFlag
Provides a code based selection for trigger and DCS information in order to have no failing filters i...
Definition: GenericTriggerEventFlag.h:42
DiJetMonitor::dijet_selection
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)
Definition: DiJetMonitor.cc:386
DiJetMonitor::dijet_eta_binning
MEbinning dijet_eta_binning
Definition: DiJetMonitor.cc:97
iEvent
int iEvent
Definition: GenABIO.cc:224
DiJetMonitor_cff.dijetPtThrPSet
dijetPtThrPSet
Definition: DiJetMonitor_cff.py:8
LumiMonitor_cff.histoPSet
histoPSet
Definition: LumiMonitor_cff.py:12
DiJetMonitor::num_genTriggerEventFlag_
std::unique_ptr< GenericTriggerEventFlag > num_genTriggerEventFlag_
Definition: DiJetMonitor.cc:82
DiJetMonitor::jetptPrbME_
ObjME jetptPrbME_
Definition: DiJetMonitor.cc:74
DiJetMonitor::ptcut_
double ptcut_
Definition: DiJetMonitor.cc:86
DiJetMonitor::jetAsyEtaME_
ObjME jetAsyEtaME_
Definition: DiJetMonitor.cc:79
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
TriggerDQMBase::getHistoPSet
static MEbinning getHistoPSet(const edm::ParameterSet &pset)
Definition: TriggerDQMBase.cc:120
DiJetMonitor::den_genTriggerEventFlag_
std::unique_ptr< GenericTriggerEventFlag > den_genTriggerEventFlag_
Definition: DiJetMonitor.cc:83
DiJetMonitor::dijetpt_binning_
MEbinning dijetpt_binning_
Definition: DiJetMonitor.cc:65
submitPVResolutionJobs.desc
string desc
Definition: submitPVResolutionJobs.py:251
TriggerDQMBase::ObjME::denominator
MonitorElement * denominator
Definition: TriggerDQMBase.h:28
TriggerDQMBase::MEbinning::nbins
uint nbins
Definition: TriggerDQMBase.h:17
TriggerDQMBase::MEbinning::xmin
double xmin
Definition: TriggerDQMBase.h:18
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
Definition: ParameterSet.h:303
DiJetMonitor::hltPathsAreValid_
bool hltPathsAreValid_
Definition: DiJetMonitor.cc:58
DiJetMonitor::dijetptThr_binning_
MEbinning dijetptThr_binning_
Definition: DiJetMonitor.cc:66
DiJetMonitor::jetptAsyME_
ObjME jetptAsyME_
Definition: DiJetMonitor.cc:75
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
trigObjTnPSource_cfi.bins
bins
Definition: trigObjTnPSource_cfi.py:20
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
TriggerDQMBase::fillHistoLSPSetDescription
static void fillHistoLSPSetDescription(edm::ParameterSetDescription &pset)
Definition: TriggerDQMBase.cc:114
edm::InputTag
Definition: InputTag.h:15
DiJetMonitor::asy_binning
MEbinning asy_binning
Definition: DiJetMonitor.cc:103
DiJetMonitor::requireValidHLTPaths_
const bool requireValidHLTPaths_
Definition: DiJetMonitor.cc:57
Generated for CMSSW Reference Manual by doxygen 1.8.16