CMS 3D CMS Logo
CTPPSCommonDQMSource.cc
Go to the documentation of this file.
1 /****************************************************************************
2  *
3  * This is a part of TotemDQM and TOTEM offline software.
4  * Authors:
5  * Jan Kašpar (jan.kaspar@gmail.com)
6  *
7  ****************************************************************************/
8 
9 #include "FWCore/Framework/interface/ESHandle.h"
10 #include "FWCore/Framework/interface/MakerMacros.h"
11 #include "FWCore/Framework/interface/Event.h"
12 #include "FWCore/Framework/interface/EventSetup.h"
13 #include "FWCore/Framework/interface/LuminosityBlock.h"
14 #include "FWCore/ParameterSet/interface/ParameterSet.h"
15 #include "FWCore/Utilities/interface/InputTag.h"
16 
17 #include "DQMServices/Core/interface/DQMOneEDAnalyzer.h"
18 #include "DQMServices/Core/interface/DQMStore.h"
19 
20 #include "DataFormats/CTPPSDetId/interface/CTPPSDetId.h"
21 #include "DataFormats/CTPPSReco/interface/CTPPSLocalTrackLite.h"
22 #include "DataFormats/ProtonReco/interface/ForwardProton.h"
23 #include "DataFormats/OnlineMetaData/interface/CTPPSRecord.h"
24 
25 #include <string>
26 
27 //----------------------------------------------------------------------------------------------------
28 
29 class CTPPSCommonDQMSource : public DQMOneEDAnalyzer<edm::LuminosityBlockCache<std::vector<int>>> {
30 public:
31  CTPPSCommonDQMSource(const edm::ParameterSet &ps);
32  ~CTPPSCommonDQMSource() override;
33 
34 protected:
35  void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override;
36  void analyze(edm::Event const &e, edm::EventSetup const &eSetup) override;
37  std::shared_ptr<std::vector<int>> globalBeginLuminosityBlock(const edm::LuminosityBlock &iLumi,
38  const edm::EventSetup &c) const override;
39  void globalEndLuminosityBlock(const edm::LuminosityBlock &iLumi, const edm::EventSetup &c) override;
40 
41  void analyzeCTPPSRecord(edm::Event const &event, edm::EventSetup const &eventSetup);
42  void analyzeTracks(edm::Event const &event, edm::EventSetup const &eventSetup);
43  void analyzeProtons(edm::Event const &event, edm::EventSetup const &eventSetup);
44 
45 private:
46  const unsigned int verbosity;
47  constexpr static int MAX_LUMIS = 6000;
48  constexpr static int MAX_VBINS = 18;
49 
50  const edm::EDGetTokenT<CTPPSRecord> ctppsRecordToken;
51  const edm::EDGetTokenT<std::vector<CTPPSLocalTrackLite>> tokenLocalTrackLite;
52  const edm::EDGetTokenT<std::vector<reco::ForwardProton>> tokenRecoProtons;
53 
54  bool makeProtonRecoPlots_;
55 
56  int currentLS;
57  int endLS;
58 
59  std::vector<int> rpstate;
60 
62  struct GlobalPlots {
63  MonitorElement *RPState = nullptr;
64  MonitorElement *events_per_bx = nullptr, *events_per_bx_short = nullptr;
65  MonitorElement *h_trackCorr_hor = nullptr, *h_trackCorr_vert = nullptr;
66 
67  void Init(DQMStore::IBooker &ibooker);
68  };
69 
70  GlobalPlots globalPlots;
71 
73  struct ArmPlots {
74  int id;
75 
76  MonitorElement *h_numRPWithTrack_top = nullptr, *h_numRPWithTrack_hor = nullptr, *h_numRPWithTrack_bot = nullptr;
77  MonitorElement *h_trackCorr = nullptr, *h_trackCorr_overlap = nullptr;
78 
79  MonitorElement *h_proton_xi = nullptr, *h_proton_t = nullptr, *h_proton_time = nullptr;
80 
81  struct TrackingRPPlots {
82  MonitorElement *h_x, *h_y;
83  };
84 
85  std::map<unsigned int, TrackingRPPlots> trackingRPPlots;
86 
87  struct TimingRPPlots {
88  MonitorElement *h_x, *h_time;
89  };
90 
91  std::map<unsigned int, TimingRPPlots> timingRPPlots;
92 
93  ArmPlots() {}
94 
95  ArmPlots(DQMStore::IBooker &ibooker, int _id, bool makeProtonRecoPlots);
96  };
97 
98  std::map<unsigned int, ArmPlots> armPlots;
99 };
100 
101 //----------------------------------------------------------------------------------------------------
102 //----------------------------------------------------------------------------------------------------
103 
104 using namespace std;
105 using namespace edm;
106 
107 const int CTPPSCommonDQMSource::MAX_LUMIS;
108 const int CTPPSCommonDQMSource::MAX_VBINS;
109 
110 //----------------------------------------------------------------------------------------------------
111 
112 void CTPPSCommonDQMSource::GlobalPlots::Init(DQMStore::IBooker &ibooker) {
113  ibooker.setCurrentFolder("CTPPS/common");
114 
115  events_per_bx = ibooker.book1D("events per BX", "rp;Event.BX", 4002, -1.5, 4000. + 0.5);
116  events_per_bx_short = ibooker.book1D("events per BX (short)", "rp;Event.BX", 102, -1.5, 100. + 0.5);
117 
118  /*
119  RP State (HV & LV & Insertion):
120  0 -> not used
121  1 -> bad
122  2 -> warning
123  3 -> ok
124  */
125  RPState = ibooker.book2D("rpstate per LS",
126  "RP State per Lumisection;Luminosity Section;",
127  MAX_LUMIS,
128  0,
129  MAX_LUMIS,
130  MAX_VBINS,
131  0.,
132  MAX_VBINS);
133  {
134  TH2F *hist = RPState->getTH2F();
135  hist->SetCanExtend(TH1::kAllAxes);
136  TAxis *ya = hist->GetYaxis();
137  ya->SetBinLabel(1, "45, 210, FR-BT");
138  ya->SetBinLabel(2, "45, 210, FR-HR");
139  ya->SetBinLabel(3, "45, 210, FR-TP");
140  ya->SetBinLabel(4, "45, 220, C1");
141  ya->SetBinLabel(5, "45, 220, FR-BT");
142  ya->SetBinLabel(6, "45, 220, FR-HR");
143  ya->SetBinLabel(7, "45, 220, FR-TP");
144  ya->SetBinLabel(8, "45, 220, NR-BP");
145  ya->SetBinLabel(9, "45, 220, NR-TP");
146  ya->SetBinLabel(10, "56, 210, FR-BT");
147  ya->SetBinLabel(11, "56, 210, FR-HR");
148  ya->SetBinLabel(12, "56, 210, FR-TP");
149  ya->SetBinLabel(13, "56, 220, C1");
150  ya->SetBinLabel(14, "56, 220, FR-BT");
151  ya->SetBinLabel(15, "56, 220, FR-HR");
152  ya->SetBinLabel(16, "56, 220, FR-TP");
153  ya->SetBinLabel(17, "56, 220, NR-BP");
154  ya->SetBinLabel(18, "56, 220, NR-TP");
155  }
156 
157  h_trackCorr_hor = ibooker.book2D("track correlation hor", "ctpps_common_rp_hor", 6, -0.5, 5.5, 6, -0.5, 5.5);
158  {
159  TH2F *hist = h_trackCorr_hor->getTH2F();
160  TAxis *xa = hist->GetXaxis(), *ya = hist->GetYaxis();
161  xa->SetBinLabel(1, "45, 210, far");
162  ya->SetBinLabel(1, "45, 210, far");
163  xa->SetBinLabel(2, "45, 220, far");
164  ya->SetBinLabel(2, "45, 220, far");
165  xa->SetBinLabel(3, "45, 220, cyl");
166  ya->SetBinLabel(3, "45, 220, cyl");
167  xa->SetBinLabel(4, "56, 210, far");
168  ya->SetBinLabel(4, "56, 210, far");
169  xa->SetBinLabel(5, "56, 220, far");
170  ya->SetBinLabel(5, "56, 220, far");
171  xa->SetBinLabel(6, "56, 220, cyl");
172  ya->SetBinLabel(6, "56, 220, cyl");
173  }
174 
175  h_trackCorr_vert = ibooker.book2D("track correlation vert", "ctpps_common_rp_vert", 8, -0.5, 7.5, 8, -0.5, 7.5);
176  {
177  TH2F *hist = h_trackCorr_vert->getTH2F();
178  TAxis *xa = hist->GetXaxis(), *ya = hist->GetYaxis();
179  xa->SetBinLabel(1, "45, 210, far, top");
180  ya->SetBinLabel(1, "45, 210, far, top");
181  xa->SetBinLabel(2, "45, 210, far, bot");
182  ya->SetBinLabel(2, "45, 210, far, bot");
183  xa->SetBinLabel(3, "45, 220, far, top");
184  ya->SetBinLabel(3, "45, 220, far, top");
185  xa->SetBinLabel(4, "45, 220, far, bot");
186  ya->SetBinLabel(4, "45, 220, far, bot");
187  xa->SetBinLabel(5, "56, 210, far, top");
188  ya->SetBinLabel(5, "56, 210, far, top");
189  xa->SetBinLabel(6, "56, 210, far, bot");
190  ya->SetBinLabel(6, "56, 210, far, bot");
191  xa->SetBinLabel(7, "56, 220, far, top");
192  ya->SetBinLabel(7, "56, 220, far, top");
193  xa->SetBinLabel(8, "56, 220, far, bot");
194  ya->SetBinLabel(8, "56, 220, far, bot");
195  }
196 }
197 
198 //----------------------------------------------------------------------------------------------------
199 
200 CTPPSCommonDQMSource::ArmPlots::ArmPlots(DQMStore::IBooker &ibooker, int _id, bool makeProtonRecoPlots) : id(_id) {
201  string name;
202  CTPPSDetId(CTPPSDetId::sdTrackingStrip, id, 0).armName(name, CTPPSDetId::nShort);
203 
204  ibooker.setCurrentFolder("CTPPS/common/sector " + name);
205 
206  string title = "ctpps_common_sector_" + name;
207 
208  h_numRPWithTrack_top =
209  ibooker.book1D("number of top RPs with tracks", title + ";number of top RPs with tracks", 5, -0.5, 4.5);
210  h_numRPWithTrack_hor =
211  ibooker.book1D("number of hor RPs with tracks", title + ";number of hor RPs with tracks", 5, -0.5, 4.5);
212  h_numRPWithTrack_bot =
213  ibooker.book1D("number of bot RPs with tracks", title + ";number of bot RPs with tracks", 5, -0.5, 4.5);
214 
215  h_trackCorr = ibooker.book2D("track correlation", title, 7, -0.5, 6.5, 7, -0.5, 6.5);
216  TH2F *h_trackCorr_h = h_trackCorr->getTH2F();
217  TAxis *xa = h_trackCorr_h->GetXaxis(), *ya = h_trackCorr_h->GetYaxis();
218  xa->SetBinLabel(1, "210, far, hor");
219  ya->SetBinLabel(1, "210, far, hor");
220  xa->SetBinLabel(2, "210, far, top");
221  ya->SetBinLabel(2, "210, far, top");
222  xa->SetBinLabel(3, "210, far, bot");
223  ya->SetBinLabel(3, "210, far, bot");
224  xa->SetBinLabel(4, "220, cyl");
225  ya->SetBinLabel(4, "220, cyl");
226  xa->SetBinLabel(5, "220, far, hor");
227  ya->SetBinLabel(5, "220, far, hor");
228  xa->SetBinLabel(6, "220, far, top");
229  ya->SetBinLabel(6, "220, far, top");
230  xa->SetBinLabel(7, "220, far, bot");
231  ya->SetBinLabel(7, "220, far, bot");
232 
233  h_trackCorr_overlap = ibooker.book2D("track correlation hor-vert overlaps", title, 7, -0.5, 6.5, 7, -0.5, 6.5);
234  h_trackCorr_h = h_trackCorr_overlap->getTH2F();
235  xa = h_trackCorr_h->GetXaxis();
236  ya = h_trackCorr_h->GetYaxis();
237  xa->SetBinLabel(1, "210, far, hor");
238  ya->SetBinLabel(1, "210, far, hor");
239  xa->SetBinLabel(2, "210, far, top");
240  ya->SetBinLabel(2, "210, far, top");
241  xa->SetBinLabel(3, "210, far, bot");
242  ya->SetBinLabel(3, "210, far, bot");
243  xa->SetBinLabel(4, "220, cyl");
244  ya->SetBinLabel(4, "220, cyl");
245  xa->SetBinLabel(5, "220, far, hor");
246  ya->SetBinLabel(5, "220, far, hor");
247  xa->SetBinLabel(6, "220, far, top");
248  ya->SetBinLabel(6, "220, far, top");
249  xa->SetBinLabel(7, "220, far, bot");
250  ya->SetBinLabel(7, "220, far, bot");
251 
252  if (makeProtonRecoPlots) {
253  h_proton_xi = ibooker.book1D("proton xi", title + ";xi", 100, 0., 0.3);
254  h_proton_t = ibooker.book1D("proton t", title + ";|t| GeV^{2}", 100, 0., 5.);
255  h_proton_time = ibooker.book1D("proton time", title + ";time (ns)", 100, -1., 1.);
256  }
257 
258  for (const unsigned int rpDecId : {2, 3, 16, 23}) {
259  unsigned int st = rpDecId / 10, rp = rpDecId % 10, rpFullDecId = id * 100 + rpDecId;
260  CTPPSDetId rpId(CTPPSDetId::sdTrackingStrip, id, st, rp);
261  string stName, rpName;
262  rpId.stationName(stName, CTPPSDetId::nShort);
263  rpId.rpName(rpName, CTPPSDetId::nShort);
264  rpName = stName + "_" + rpName;
265 
266  if (rp == 6) {
267  timingRPPlots[rpFullDecId] = {
268  ibooker.book1D(rpName + " - track x histogram", title + "/" + rpName + ";track x (mm)", 200, 0., 40.),
269  ibooker.book1D(
270  rpName + " - track time histogram", title + "/" + rpName + ";track time (ns)", 100, -25., +50.)};
271  } else {
272  trackingRPPlots[rpFullDecId] = {
273  ibooker.book1D(rpName + " - track x histogram", title + "/" + rpName + ";track x (mm)", 200, 0., 40.),
274  ibooker.book1D(rpName + " - track y histogram", title + "/" + rpName + ";track y (mm)", 200, -20., +20.)};
275  }
276  }
277 }
278 
279 //----------------------------------------------------------------------------------------------------
280 //----------------------------------------------------------------------------------------------------
281 
282 CTPPSCommonDQMSource::CTPPSCommonDQMSource(const edm::ParameterSet &ps)
283  : verbosity(ps.getUntrackedParameter<unsigned int>("verbosity", 0)),
284  ctppsRecordToken(consumes<CTPPSRecord>(ps.getUntrackedParameter<edm::InputTag>("ctppsmetadata"))),
285  tokenLocalTrackLite(consumes<vector<CTPPSLocalTrackLite>>(ps.getParameter<edm::InputTag>("tagLocalTrackLite"))),
286  tokenRecoProtons(consumes<std::vector<reco::ForwardProton>>(ps.getParameter<InputTag>("tagRecoProtons"))),
287  makeProtonRecoPlots_(ps.getParameter<bool>("makeProtonRecoPlots")) {
288  currentLS = 0;
289  endLS = 0;
290  rpstate.clear();
291 }
292 
293 //----------------------------------------------------------------------------------------------------
294 
295 CTPPSCommonDQMSource::~CTPPSCommonDQMSource() {}
296 
297 //----------------------------------------------------------------------------------------------------
298 
299 void CTPPSCommonDQMSource::bookHistograms(DQMStore::IBooker &ibooker, edm::Run const &, edm::EventSetup const &) {
300  // global plots
301  globalPlots.Init(ibooker);
302 
303  // loop over arms
304  for (unsigned int arm = 0; arm < 2; arm++) {
305  armPlots[arm] = ArmPlots(ibooker, arm, makeProtonRecoPlots_);
306  }
307 }
308 
309 //----------------------------------------------------------------------------------------------------
310 
311 void CTPPSCommonDQMSource::analyze(edm::Event const &event, edm::EventSetup const &eventSetup) {
312  analyzeCTPPSRecord(event, eventSetup);
313  analyzeTracks(event, eventSetup);
314 
315  if (makeProtonRecoPlots_)
316  analyzeProtons(event, eventSetup);
317 }
318 
319 //----------------------------------------------------------------------------------------------------
320 
321 void CTPPSCommonDQMSource::analyzeCTPPSRecord(edm::Event const &event, edm::EventSetup const &eventSetup) {
322  Handle<CTPPSRecord> hCTPPSRecord;
323  event.getByToken(ctppsRecordToken, hCTPPSRecord);
324 
325  if (!hCTPPSRecord.isValid()) {
326  if (verbosity)
327  LogProblem("CTPPSCommonDQMSource") << "ERROR in CTPPSCommonDQMSource::analyzeCTPPSRecord > input not available.";
328 
329  return;
330  }
331 
332  auto &rpstate = *luminosityBlockCache(event.getLuminosityBlock().index());
333  if (rpstate.empty()) {
334  rpstate.reserve(CTPPSRecord::RomanPot::Last);
335  for (uint8_t i = 0; i < CTPPSRecord::RomanPot::Last; ++i)
336  rpstate.push_back(hCTPPSRecord->status(i));
337  }
338 }
339 
340 //----------------------------------------------------------------------------------------------------
341 
342 void CTPPSCommonDQMSource::analyzeTracks(edm::Event const &event, edm::EventSetup const &eventSetup) {
343  // get event data
344  Handle<vector<CTPPSLocalTrackLite>> hTracks;
345  event.getByToken(tokenLocalTrackLite, hTracks);
346 
347  // check validity
348  if (!hTracks.isValid()) {
349  if (verbosity)
350  LogProblem("CTPPSCommonDQMSource") << "ERROR in CTPPSCommonDQMSource::analyzeTracks > input not available.";
351 
352  return;
353  }
354 
355  //------------------------------
356  // collect indeces of RP with tracks, for each correlation plot
357  set<signed int> s_rp_idx_global_hor, s_rp_idx_global_vert;
358  map<unsigned int, set<signed int>> ms_rp_idx_arm;
359 
360  for (auto &tr : *hTracks) {
361  const CTPPSDetId rpId(tr.rpId());
362  const unsigned int arm = rpId.arm();
363  const unsigned int stNum = rpId.station();
364  const unsigned int rpNum = rpId.rp();
365  const unsigned int stRPNum = stNum * 10 + rpNum;
366 
367  {
368  signed int idx = -1;
369  if (stRPNum == 3)
370  idx = 0;
371  if (stRPNum == 23)
372  idx = 1;
373  if (stRPNum == 16)
374  idx = 2;
375 
376  if (idx >= 0)
377  s_rp_idx_global_hor.insert(3 * arm + idx);
378  }
379 
380  {
381  signed int idx = -1;
382  if (stRPNum == 4)
383  idx = 0;
384  if (stRPNum == 5)
385  idx = 1;
386  if (stRPNum == 24)
387  idx = 2;
388  if (stRPNum == 25)
389  idx = 3;
390 
391  if (idx >= 0)
392  s_rp_idx_global_vert.insert(4 * arm + idx);
393  }
394 
395  {
396  signed int idx = -1;
397  if (stRPNum == 3)
398  idx = 0;
399  if (stRPNum == 4)
400  idx = 1;
401  if (stRPNum == 5)
402  idx = 2;
403  if (stRPNum == 16)
404  idx = 3;
405  if (stRPNum == 23)
406  idx = 4;
407  if (stRPNum == 24)
408  idx = 5;
409  if (stRPNum == 25)
410  idx = 6;
411 
412  const signed int hor = ((rpNum == 2) || (rpNum == 3) || (rpNum == 6)) ? 1 : 0;
413 
414  if (idx >= 0)
415  ms_rp_idx_arm[arm].insert(idx * 10 + hor);
416  }
417  }
418 
419  //------------------------------
420  // Global Plots
421 
422  globalPlots.events_per_bx->Fill(event.bunchCrossing());
423  globalPlots.events_per_bx_short->Fill(event.bunchCrossing());
424 
425  for (const auto &idx1 : s_rp_idx_global_hor)
426  for (const auto &idx2 : s_rp_idx_global_hor)
427  globalPlots.h_trackCorr_hor->Fill(idx1, idx2);
428 
429  for (const auto &idx1 : s_rp_idx_global_vert)
430  for (const auto &idx2 : s_rp_idx_global_vert)
431  globalPlots.h_trackCorr_vert->Fill(idx1, idx2);
432 
433  //------------------------------
434  // Arm Plots
435 
436  map<unsigned int, set<unsigned int>> mTop, mHor, mBot;
437 
438  for (auto &tr : *hTracks) {
439  CTPPSDetId rpId(tr.rpId());
440  const unsigned int rpNum = rpId.rp();
441  const unsigned int armIdx = rpId.arm();
442 
443  if (rpNum == 0 || rpNum == 4)
444  mTop[armIdx].insert(rpId);
445  if (rpNum == 2 || rpNum == 3 || rpNum == 6)
446  mHor[armIdx].insert(rpId);
447  if (rpNum == 1 || rpNum == 5)
448  mBot[armIdx].insert(rpId);
449 
450  auto &ap = armPlots[rpId.arm()];
451  unsigned int rpDecId = rpId.arm() * 100 + rpId.station() * 10 + rpId.rp();
452 
453  // fill in reference tracking-RP plots
454  {
455  auto it = ap.trackingRPPlots.find(rpDecId);
456  if (it != ap.trackingRPPlots.end()) {
457  it->second.h_x->Fill(tr.x());
458  it->second.h_y->Fill(tr.y());
459  }
460  }
461 
462  // fill in reference timing-RP plots
463  {
464  auto it = ap.timingRPPlots.find(rpDecId);
465  if (it != ap.timingRPPlots.end()) {
466  it->second.h_x->Fill(tr.x());
467  it->second.h_time->Fill(tr.time());
468  }
469  }
470  }
471 
472  for (auto &p : armPlots) {
473  p.second.h_numRPWithTrack_top->Fill(mTop[p.first].size());
474  p.second.h_numRPWithTrack_hor->Fill(mHor[p.first].size());
475  p.second.h_numRPWithTrack_bot->Fill(mBot[p.first].size());
476  }
477 
478  //------------------------------
479  // Correlation plots
480 
481  for (const auto &ap : ms_rp_idx_arm) {
482  auto &plots = armPlots[ap.first];
483 
484  for (const auto &idx1 : ap.second) {
485  for (const auto &idx2 : ap.second) {
486  plots.h_trackCorr->Fill(idx1 / 10, idx2 / 10);
487 
488  if ((idx1 % 10) != (idx2 % 10))
489  plots.h_trackCorr_overlap->Fill(idx1 / 10, idx2 / 10);
490  }
491  }
492  }
493 }
494 
495 //----------------------------------------------------------------------------------------------------
496 
497 void CTPPSCommonDQMSource::analyzeProtons(edm::Event const &event, edm::EventSetup const &eventSetup) {
498  // get event data
499  Handle<vector<reco::ForwardProton>> hRecoProtons;
500  event.getByToken(tokenRecoProtons, hRecoProtons);
501 
502  // check validity
503  if (!hRecoProtons.isValid()) {
504  if (verbosity)
505  LogProblem("CTPPSCommonDQMSource") << "ERROR in CTPPSCommonDQMSource::analyzeProtons > input not available.";
506 
507  return;
508  }
509 
510  // loop over protons
511  for (auto &p : *hRecoProtons) {
512  if (!p.validFit())
513  continue;
514 
515  signed int armIndex = -1;
516  if (p.lhcSector() == reco::ForwardProton::LHCSector::sector45)
517  armIndex = 0;
518  if (p.lhcSector() == reco::ForwardProton::LHCSector::sector56)
519  armIndex = 1;
520  if (armIndex < 0)
521  continue;
522 
523  auto &plots = armPlots[armIndex];
524 
525  plots.h_proton_xi->Fill(p.xi());
526  plots.h_proton_t->Fill(fabs(p.t()));
527  plots.h_proton_time->Fill(p.time());
528  }
529 }
530 
531 //----------------------------------------------------------------------------------------------------
532 
533 std::shared_ptr<std::vector<int>> CTPPSCommonDQMSource::globalBeginLuminosityBlock(const edm::LuminosityBlock &,
534  const edm::EventSetup &) const {
535  return std::make_shared<std::vector<int>>();
536 }
537 
538 //----------------------------------------------------------------------------------------------------
539 
540 void CTPPSCommonDQMSource::globalEndLuminosityBlock(const edm::LuminosityBlock &iLumi, const edm::EventSetup &c) {
541  auto const &rpstate = *luminosityBlockCache(iLumi.index());
542  auto currentLS = iLumi.id().luminosityBlock();
543  for (std::vector<int>::size_type i = 0; i < rpstate.size(); i++)
544  globalPlots.RPState->setBinContent(currentLS, i + 1, rpstate[i]);
545 }
546 
547 //----------------------------------------------------------------------------------------------------
548 
549 DEFINE_FWK_MODULE(CTPPSCommonDQMSource);
CTPPSCommonDQMSource::analyzeTracks
void analyzeTracks(edm::Event const &event, edm::EventSetup const &eventSetup)
Definition: CTPPSCommonDQMSource.cc:341
ctppsCommonDQMSource_cfi.makeProtonRecoPlots
makeProtonRecoPlots
Definition: ctppsCommonDQMSource_cfi.py:10
CTPPSLocalTrackLite
Local (=single RP) track with essential information only.
Definition: CTPPSLocalTrackLite.h:19
dqm::impl::MonitorElement
Definition: MonitorElement.h:98
CTPPSCommonDQMSource::tokenRecoProtons
const edm::EDGetTokenT< std::vector< reco::ForwardProton > > tokenRecoProtons
Definition: CTPPSCommonDQMSource.cc:56
CTPPSRecord::status
Status status(const uint8_t rp) const
Return the status of the given roman pot.
Definition: CTPPSRecord.h:61
CTPPSCommonDQMSource::ArmPlots::h_proton_xi
MonitorElement * h_proton_xi
Definition: CTPPSCommonDQMSource.cc:83
electrons_cff.bool
bool
Definition: electrons_cff.py:372
mps_fire.i
i
Definition: mps_fire.py:355
CTPPSCommonDQMSource::ArmPlots::h_proton_time
MonitorElement * h_proton_time
Definition: CTPPSCommonDQMSource.cc:83
CTPPSCommonDQMSource::ArmPlots::TimingRPPlots
Definition: CTPPSCommonDQMSource.cc:91
ESHandle.h
CTPPSCommonDQMSource::GlobalPlots::Init
void Init(DQMStore::IBooker &ibooker)
Definition: CTPPSCommonDQMSource.cc:111
CTPPSCommonDQMSource::ArmPlots::TrackingRPPlots::h_y
MonitorElement * h_y
Definition: CTPPSCommonDQMSource.cc:86
edm::LuminosityBlock
Definition: LuminosityBlock.h:50
CTPPSCommonDQMSource::analyze
void analyze(edm::Event const &e, edm::EventSetup const &eSetup) override
Definition: CTPPSCommonDQMSource.cc:310
edm::Run
Definition: Run.h:45
edm::EDGetTokenT< CTPPSRecord >
LuminosityBlock.h
edm::LuminosityBlock::index
LuminosityBlockIndex index() const
Definition: LuminosityBlock.cc:20
CTPPSCommonDQMSource::ArmPlots::trackingRPPlots
std::map< unsigned int, TrackingRPPlots > trackingRPPlots
Definition: CTPPSCommonDQMSource.cc:89
edm
HLT enums.
Definition: AlignableModifier.h:19
CTPPSCommonDQMSource::CTPPSCommonDQMSource
CTPPSCommonDQMSource(const edm::ParameterSet &ps)
Definition: CTPPSCommonDQMSource.cc:281
AlCaHLTBitMon_ParallelJobs.p
p
Definition: AlCaHLTBitMon_ParallelJobs.py:153
CTPPSCommonDQMSource::rpstate
std::vector< int > rpstate
Definition: CTPPSCommonDQMSource.cc:63
DQMOneEDAnalyzer
Definition: DQMOneEDAnalyzer.h:20
CTPPSCommonDQMSource::ArmPlots::TimingRPPlots::h_time
MonitorElement * h_time
Definition: CTPPSCommonDQMSource.cc:92
DQMStore.h
CTPPSDetId::nShort
Definition: CTPPSDetId.h:87
CTPPSLocalTrackLite.h
CTPPSCommonDQMSource::GlobalPlots::h_trackCorr_hor
MonitorElement * h_trackCorr_hor
Definition: CTPPSCommonDQMSource.cc:69
reco
fixed size matrix
Definition: AlignmentAlgorithmBase.h:45
year_2016_postTS2_cff.rpId
rpId
Definition: year_2016_postTS2_cff.py:23
CTPPSCommonDQMSource::GlobalPlots::events_per_bx_short
MonitorElement * events_per_bx_short
Definition: CTPPSCommonDQMSource.cc:68
edm::Handle
Definition: AssociativeIterator.h:50
training_settings.idx
idx
Definition: training_settings.py:16
DQMOneEDAnalyzer.h
CTPPSRecord
Class to contain CTPPS information from soft FED 1022.
Definition: CTPPSRecord.h:20
CTPPSDetId::sdTrackingStrip
Definition: CTPPSDetId.h:44
dqm::impl::MonitorElement::getTH2F
virtual TH2F * getTH2F()
Definition: MonitorElement.cc:992
dqm::legacy::DQMStore::IBooker
dqm::implementation::IBooker IBooker
Definition: DQMStore.h:729
MakerMacros.h
trigger::size_type
uint16_t size_type
Definition: TriggerTypeDefs.h:18
CTPPSCommonDQMSource::ArmPlots::ArmPlots
ArmPlots()
Definition: CTPPSCommonDQMSource.cc:97
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:16
CTPPSCommonDQMSource::globalBeginLuminosityBlock
std::shared_ptr< std::vector< int > > globalBeginLuminosityBlock(const edm::LuminosityBlock &iLumi, const edm::EventSetup &c) const override
Definition: CTPPSCommonDQMSource.cc:532
compare.hist
hist
Definition: compare.py:376
CTPPSCommonDQMSource::analyzeCTPPSRecord
void analyzeCTPPSRecord(edm::Event const &event, edm::EventSetup const &eventSetup)
Definition: CTPPSCommonDQMSource.cc:320
CTPPSCommonDQMSource::tokenLocalTrackLite
const edm::EDGetTokenT< std::vector< CTPPSLocalTrackLite > > tokenLocalTrackLite
Definition: CTPPSCommonDQMSource.cc:55
TtSemiLepHitFitProducer_Electrons_cfi.mTop
mTop
Definition: TtSemiLepHitFitProducer_Electrons_cfi.py:36
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:290
CTPPSCommonDQMSource::MAX_LUMIS
constexpr static int MAX_LUMIS
Definition: CTPPSCommonDQMSource.cc:51
CTPPSCommonDQMSource::ctppsRecordToken
const edm::EDGetTokenT< CTPPSRecord > ctppsRecordToken
Definition: CTPPSCommonDQMSource.cc:54
CTPPSCommonDQMSource::GlobalPlots::h_trackCorr_vert
MonitorElement * h_trackCorr_vert
Definition: CTPPSCommonDQMSource.cc:69
CTPPSCommonDQMSource::ArmPlots::h_trackCorr_overlap
MonitorElement * h_trackCorr_overlap
Definition: CTPPSCommonDQMSource.cc:81
CTPPSCommonDQMSource::ArmPlots::TimingRPPlots::h_x
MonitorElement * h_x
Definition: CTPPSCommonDQMSource.cc:92
CTPPSCommonDQMSource::analyzeProtons
void analyzeProtons(edm::Event const &event, edm::EventSetup const &eventSetup)
Definition: CTPPSCommonDQMSource.cc:496
CTPPSCommonDQMSource::ArmPlots::h_numRPWithTrack_bot
MonitorElement * h_numRPWithTrack_bot
Definition: CTPPSCommonDQMSource.cc:80
CTPPSRecord.h
CTPPSCommonDQMSource::currentLS
int currentLS
Definition: CTPPSCommonDQMSource.cc:60
edm::ParameterSet
Definition: ParameterSet.h:36
HLTObjectsMonitor_cfi.plots
plots
Definition: HLTObjectsMonitor_cfi.py:17
Event.h
CTPPSCommonDQMSource::MAX_VBINS
constexpr static int MAX_VBINS
Definition: CTPPSCommonDQMSource.cc:52
CTPPSCommonDQMSource::GlobalPlots
plots related to the whole system
Definition: CTPPSCommonDQMSource.cc:66
CTPPSCommonDQMSource::ArmPlots::h_proton_t
MonitorElement * h_proton_t
Definition: CTPPSCommonDQMSource.cc:83
edm::LogProblem
Definition: MessageLogger.h:387
ForwardProton.h
reco::ForwardProton::LHCSector::sector56
CTPPSDetId
Base class for CTPPS detector IDs.
Definition: CTPPSDetId.h:31
createfilelist.int
int
Definition: createfilelist.py:10
edm::LuminosityBlockID::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: LuminosityBlockID.h:42
edm::LuminosityBlockBase::id
LuminosityBlockID id() const
Definition: LuminosityBlockBase.h:44
CTPPSCommonDQMSource::globalPlots
GlobalPlots globalPlots
Definition: CTPPSCommonDQMSource.cc:74
edm::EventSetup
Definition: EventSetup.h:57
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
CTPPSCommonDQMSource::ArmPlots::h_trackCorr
MonitorElement * h_trackCorr
Definition: CTPPSCommonDQMSource.cc:81
InputTag.h
CTPPSCommonDQMSource::ArmPlots::h_numRPWithTrack_hor
MonitorElement * h_numRPWithTrack_hor
Definition: CTPPSCommonDQMSource.cc:80
overlapproblemtsosanalyzer_cfi.title
title
Definition: overlapproblemtsosanalyzer_cfi.py:7
dqm::impl::MonitorElement::setBinContent
virtual void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:691
std
Definition: JetResolutionObject.h:76
triggerObjects_cff.id
id
Definition: triggerObjects_cff.py:31
reco::ForwardProton::LHCSector::sector45
CTPPSCommonDQMSource
Definition: CTPPSCommonDQMSource.cc:28
CTPPSCommonDQMSource::globalEndLuminosityBlock
void globalEndLuminosityBlock(const edm::LuminosityBlock &iLumi, const edm::EventSetup &c) override
Definition: CTPPSCommonDQMSource.cc:539
CTPPSCommonDQMSource::~CTPPSCommonDQMSource
~CTPPSCommonDQMSource() override
Definition: CTPPSCommonDQMSource.cc:294
CTPPSCommonDQMSource::endLS
int endLS
Definition: CTPPSCommonDQMSource.cc:61
CTPPSCommonDQMSource::GlobalPlots::RPState
MonitorElement * RPState
Definition: CTPPSCommonDQMSource.cc:67
CTPPSCommonDQMSource::ArmPlots::TrackingRPPlots
Definition: CTPPSCommonDQMSource.cc:85
CTPPSCommonDQMSource::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: CTPPSCommonDQMSource.cc:298
Skims_PA_cff.name
name
Definition: Skims_PA_cff.py:17
EventSetup.h
CTPPSDetId.h
CTPPSCommonDQMSource::makeProtonRecoPlots_
bool makeProtonRecoPlots_
Definition: CTPPSCommonDQMSource.cc:58
ParameterSet.h
edm::EDConsumerBase::consumes
EDGetTokenT< ProductType > consumes(edm::InputTag const &tag)
Definition: EDConsumerBase.h:126
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:70
event
Definition: event.py:1
edm::Event
Definition: Event.h:73
CTPPSCommonDQMSource::GlobalPlots::events_per_bx
MonitorElement * events_per_bx
Definition: CTPPSCommonDQMSource.cc:68
edm::InputTag
Definition: InputTag.h:15
CTPPSCommonDQMSource::verbosity
const unsigned int verbosity
Definition: CTPPSCommonDQMSource.cc:50
CTPPSCommonDQMSource::ArmPlots::id
int id
Definition: CTPPSCommonDQMSource.cc:78
CTPPSCommonDQMSource::ArmPlots::TrackingRPPlots::h_x
MonitorElement * h_x
Definition: CTPPSCommonDQMSource.cc:86
CTPPSCommonDQMSource::armPlots
std::map< unsigned int, ArmPlots > armPlots
Definition: CTPPSCommonDQMSource.cc:102
CTPPSCommonDQMSource::ArmPlots::timingRPPlots
std::map< unsigned int, TimingRPPlots > timingRPPlots
Definition: CTPPSCommonDQMSource.cc:95
CTPPSCommonDQMSource::ArmPlots
plots related to one arm
Definition: CTPPSCommonDQMSource.cc:77
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
CTPPSDetId::armName
void armName(std::string &name, NameFlag flag=nFull) const
Definition: CTPPSDetId.h:96
CTPPSCommonDQMSource::ArmPlots::h_numRPWithTrack_top
MonitorElement * h_numRPWithTrack_top
Definition: CTPPSCommonDQMSource.cc:80
Generated for CMSSW Reference Manual by doxygen 1.8.16