CMS 3D CMS Logo
L1TdeCSCTPGShower.cc
Go to the documentation of this file.
1 #include <string>
2 
3 #include "DQM/L1TMonitor/interface/L1TdeCSCTPGShower.h"
4 #include "DataFormats/MuonDetId/interface/CSCTriggerNumbering.h"
5 
6 L1TdeCSCTPGShower::L1TdeCSCTPGShower(const edm::ParameterSet& ps)
7  : dataALCTShower_token_(consumes<CSCShowerDigiCollection>(ps.getParameter<edm::InputTag>("dataALCTShower"))),
8  emulALCTShower_token_(consumes<CSCShowerDigiCollection>(ps.getParameter<edm::InputTag>("emulALCTShower"))),
9  dataCLCTShower_token_(consumes<CSCShowerDigiCollection>(ps.getParameter<edm::InputTag>("dataCLCTShower"))),
10  emulCLCTShower_token_(consumes<CSCShowerDigiCollection>(ps.getParameter<edm::InputTag>("emulCLCTShower"))),
11  dataLCTShower_token_(consumes<CSCShowerDigiCollection>(ps.getParameter<edm::InputTag>("dataLCTShower"))),
12  emulLCTShower_token_(consumes<CSCShowerDigiCollection>(ps.getParameter<edm::InputTag>("emulLCTShower"))),
13  monitorDir_(ps.getUntrackedParameter<std::string>("monitorDir")) {}
14 
15 L1TdeCSCTPGShower::~L1TdeCSCTPGShower() {}
16 
17 void L1TdeCSCTPGShower::bookHistograms(DQMStore::IBooker& iBooker, const edm::Run&, const edm::EventSetup&) {
18  iBooker.setCurrentFolder(monitorDir_);
19 
20  // 2D summary plots
21  lctShowerDataNomSummary_denom_ =
22  iBooker.book2D("lct_cscshower_data_nom_summary_denom", "Data LCT Nominal Shower All", 36, 1, 37, 18, 0, 18);
23  lctShowerDataNomSummary_num_ = iBooker.book2D(
24  "lct_cscshower_data_nom_summary_num", "Data LCT Nominal Shower Emul Matched", 36, 1, 37, 18, 0, 18);
25  alctShowerDataNomSummary_denom_ =
26  iBooker.book2D("alct_cscshower_data_nom_summary_denom", "Data ALCT Nominal Shower All", 36, 1, 37, 18, 0, 18);
27  alctShowerDataNomSummary_num_ = iBooker.book2D(
28  "alct_cscshower_data_nom_summary_num", "Data ALCT Nominal Shower Emul Matched", 36, 1, 37, 18, 0, 18);
29  clctShowerDataNomSummary_denom_ =
30  iBooker.book2D("clct_cscshower_data_nom_summary_denom", "Data CLCT Nominal Shower All", 36, 1, 37, 18, 0, 18);
31  clctShowerDataNomSummary_num_ = iBooker.book2D(
32  "clct_cscshower_data_nom_summary_num", "Data CLCT Nominal Shower Emul Matched", 36, 1, 37, 18, 0, 18);
33 
34  lctShowerEmulNomSummary_denom_ =
35  iBooker.book2D("lct_cscshower_emul_nom_summary_denom", "Emul LCT Nominal Shower All", 36, 1, 37, 18, 0, 18);
36  lctShowerEmulNomSummary_num_ = iBooker.book2D(
37  "lct_cscshower_emul_nom_summary_num", "Emul LCT Nominal Shower Not Matched to Data", 36, 1, 37, 18, 0, 18);
38  alctShowerEmulNomSummary_denom_ =
39  iBooker.book2D("alct_cscshower_emul_nom_summary_denom", "Emul ALCT Nominal Shower All", 36, 1, 37, 18, 0, 18);
40  alctShowerEmulNomSummary_num_ = iBooker.book2D(
41  "alct_cscshower_emul_nom_summary_num", "Emul ALCT Nominal Shower Not Matched to Data", 36, 1, 37, 18, 0, 18);
42  clctShowerEmulNomSummary_denom_ =
43  iBooker.book2D("clct_cscshower_emul_nom_summary_denom", "Emul CLCT Nominal Shower All", 36, 1, 37, 18, 0, 18);
44  clctShowerEmulNomSummary_num_ = iBooker.book2D(
45  "clct_cscshower_emul_nom_summary_num", "Emul CLCT Nominal Shower Not Matched to Data", 36, 1, 37, 18, 0, 18);
46 
47  lctShowerDataTightSummary_denom_ =
48  iBooker.book2D("lct_cscshower_data_tight_summary_denom", "Data LCT Tight Shower All", 36, 1, 37, 18, 0, 18);
49  lctShowerDataTightSummary_num_ = iBooker.book2D(
50  "lct_cscshower_data_tight_summary_num", "Data LCT Tight Shower Emul Matched", 36, 1, 37, 18, 0, 18);
51  alctShowerDataTightSummary_denom_ =
52  iBooker.book2D("alct_cscshower_data_tight_summary_denom", "Data ALCT Tight Shower All", 36, 1, 37, 18, 0, 18);
53  alctShowerDataTightSummary_num_ = iBooker.book2D(
54  "alct_cscshower_data_tight_summary_num", "Data ALCT Tight Shower Emul Matched", 36, 1, 37, 18, 0, 18);
55  clctShowerDataTightSummary_denom_ =
56  iBooker.book2D("clct_cscshower_data_tight_summary_denom", "Data CLCT Tight Shower All", 36, 1, 37, 18, 0, 18);
57  clctShowerDataTightSummary_num_ = iBooker.book2D(
58  "clct_cscshower_data_tight_summary_num", "Data CLCT Tight Shower Emul Matched", 36, 1, 37, 18, 0, 18);
59 
60  lctShowerEmulTightSummary_denom_ =
61  iBooker.book2D("lct_cscshower_emul_tight_summary_denom", "Emul LCT Tight Shower All", 36, 1, 37, 18, 0, 18);
62  lctShowerEmulTightSummary_num_ = iBooker.book2D(
63  "lct_cscshower_emul_tight_summary_num", "Emul LCT Tight Shower Not Matched to Data", 36, 1, 37, 18, 0, 18);
64  alctShowerEmulTightSummary_denom_ =
65  iBooker.book2D("alct_cscshower_emul_tight_summary_denom", "Emul ALCT Tight Shower All", 36, 1, 37, 18, 0, 18);
66  alctShowerEmulTightSummary_num_ = iBooker.book2D(
67  "alct_cscshower_emul_tight_summary_num", "Emul ALCT Tight Shower Not Matched to Data", 36, 1, 37, 18, 0, 18);
68  clctShowerEmulTightSummary_denom_ =
69  iBooker.book2D("clct_cscshower_emul_tight_summary_denom", "Emul CLCT Tight Shower All", 36, 1, 37, 18, 0, 18);
70  clctShowerEmulTightSummary_num_ = iBooker.book2D(
71  "clct_cscshower_emul_tight_summary_num", "Emul CLCT Tight Shower Not Matched to Data", 36, 1, 37, 18, 0, 18);
72 
73  // x labels
74  lctShowerDataNomSummary_denom_->setAxisTitle("Chamber", 1);
75  lctShowerDataNomSummary_num_->setAxisTitle("Chamber", 1);
76  alctShowerDataNomSummary_denom_->setAxisTitle("Chamber", 1);
77  alctShowerDataNomSummary_num_->setAxisTitle("Chamber", 1);
78  clctShowerDataNomSummary_denom_->setAxisTitle("Chamber", 1);
79  clctShowerDataNomSummary_num_->setAxisTitle("Chamber", 1);
80 
81  lctShowerEmulNomSummary_denom_->setAxisTitle("Chamber", 1);
82  lctShowerEmulNomSummary_num_->setAxisTitle("Chamber", 1);
83  alctShowerEmulNomSummary_denom_->setAxisTitle("Chamber", 1);
84  alctShowerEmulNomSummary_num_->setAxisTitle("Chamber", 1);
85  clctShowerEmulNomSummary_denom_->setAxisTitle("Chamber", 1);
86  clctShowerEmulNomSummary_num_->setAxisTitle("Chamber", 1);
87 
88  lctShowerDataTightSummary_denom_->setAxisTitle("Chamber", 1);
89  lctShowerDataTightSummary_num_->setAxisTitle("Chamber", 1);
90  alctShowerDataTightSummary_denom_->setAxisTitle("Chamber", 1);
91  alctShowerDataTightSummary_num_->setAxisTitle("Chamber", 1);
92  clctShowerDataTightSummary_denom_->setAxisTitle("Chamber", 1);
93  clctShowerDataTightSummary_num_->setAxisTitle("Chamber", 1);
94 
95  lctShowerEmulTightSummary_denom_->setAxisTitle("Chamber", 1);
96  lctShowerEmulTightSummary_num_->setAxisTitle("Chamber", 1);
97  alctShowerEmulTightSummary_denom_->setAxisTitle("Chamber", 1);
98  alctShowerEmulTightSummary_num_->setAxisTitle("Chamber", 1);
99  clctShowerEmulTightSummary_denom_->setAxisTitle("Chamber", 1);
100  clctShowerEmulTightSummary_num_->setAxisTitle("Chamber", 1);
101 
102  // plotting option
103  lctShowerDataNomSummary_denom_->setOption("colz");
104  lctShowerDataNomSummary_num_->setOption("colz");
105  alctShowerDataNomSummary_denom_->setOption("colz");
106  alctShowerDataNomSummary_num_->setOption("colz");
107  clctShowerDataNomSummary_denom_->setOption("colz");
108  clctShowerDataNomSummary_num_->setOption("colz");
109 
110  lctShowerEmulNomSummary_denom_->setOption("colz");
111  lctShowerEmulNomSummary_num_->setOption("colz");
112  alctShowerEmulNomSummary_denom_->setOption("colz");
113  alctShowerEmulNomSummary_num_->setOption("colz");
114  clctShowerEmulNomSummary_denom_->setOption("colz");
115  clctShowerEmulNomSummary_num_->setOption("colz");
116 
117  lctShowerDataTightSummary_denom_->setOption("colz");
118  lctShowerDataTightSummary_num_->setOption("colz");
119  alctShowerDataTightSummary_denom_->setOption("colz");
120  alctShowerDataTightSummary_num_->setOption("colz");
121  clctShowerDataTightSummary_denom_->setOption("colz");
122  clctShowerDataTightSummary_num_->setOption("colz");
123 
124  lctShowerEmulTightSummary_denom_->setOption("colz");
125  lctShowerEmulTightSummary_num_->setOption("colz");
126  alctShowerEmulTightSummary_denom_->setOption("colz");
127  alctShowerEmulTightSummary_num_->setOption("colz");
128  clctShowerEmulTightSummary_denom_->setOption("colz");
129  clctShowerEmulTightSummary_num_->setOption("colz");
130 
131  const std::array<std::string, 9> suffix_label{{"4/2", "4/1", "3/2", "3/1", " 2/2", "2/1", "1/3", "1/2", "1/1"}};
132 
133  // y labels
134  for (int ybin = 1; ybin <= 9; ++ybin) {
135  lctShowerDataNomSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
136  lctShowerDataNomSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
137  alctShowerDataNomSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
138  alctShowerDataNomSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
139  clctShowerDataNomSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
140  clctShowerDataNomSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
141 
142  lctShowerEmulNomSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
143  lctShowerEmulNomSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
144  alctShowerEmulNomSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
145  alctShowerEmulNomSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
146  clctShowerEmulNomSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
147  clctShowerEmulNomSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
148 
149  lctShowerDataNomSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
150  lctShowerDataNomSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
151  alctShowerDataNomSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
152  alctShowerDataNomSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
153  clctShowerDataNomSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
154  clctShowerDataNomSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
155 
156  lctShowerEmulNomSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
157  lctShowerEmulNomSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
158  alctShowerEmulNomSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
159  alctShowerEmulNomSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
160  clctShowerEmulNomSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
161  clctShowerEmulNomSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
162 
163  lctShowerDataTightSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
164  lctShowerDataTightSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
165  alctShowerDataTightSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
166  alctShowerDataTightSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
167  clctShowerDataTightSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
168  clctShowerDataTightSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
169 
170  lctShowerEmulTightSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
171  lctShowerEmulTightSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
172  alctShowerEmulTightSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
173  alctShowerEmulTightSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
174  clctShowerEmulTightSummary_denom_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
175  clctShowerEmulTightSummary_num_->setBinLabel(ybin, "ME-" + suffix_label[ybin - 1], 2);
176 
177  lctShowerDataTightSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
178  lctShowerDataTightSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
179  alctShowerDataTightSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
180  alctShowerDataTightSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
181  clctShowerDataTightSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
182  clctShowerDataTightSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
183 
184  lctShowerEmulTightSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
185  lctShowerEmulTightSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
186  alctShowerEmulTightSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
187  alctShowerEmulTightSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
188  clctShowerEmulTightSummary_denom_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
189  clctShowerEmulTightSummary_num_->setBinLabel(19 - ybin, "ME+" + suffix_label[ybin - 1], 2);
190  }
191 }
192 
193 void L1TdeCSCTPGShower::analyze(const edm::Event& e, const edm::EventSetup& c) {
194  // handles
195  edm::Handle<CSCShowerDigiCollection> dataALCTshs;
196  edm::Handle<CSCShowerDigiCollection> emulALCTshs;
197  edm::Handle<CSCShowerDigiCollection> dataCLCTshs;
198  edm::Handle<CSCShowerDigiCollection> emulCLCTshs;
199  edm::Handle<CSCShowerDigiCollection> dataLCTshs;
200  edm::Handle<CSCShowerDigiCollection> emulLCTshs;
201 
202  e.getByToken(dataALCTShower_token_, dataALCTshs);
203  e.getByToken(emulALCTShower_token_, emulALCTshs);
204  e.getByToken(dataCLCTShower_token_, dataCLCTshs);
205  e.getByToken(emulCLCTShower_token_, emulCLCTshs);
206  e.getByToken(dataLCTShower_token_, dataLCTshs);
207  e.getByToken(emulLCTShower_token_, emulLCTshs);
208 
209  const std::map<std::pair<int, int>, int> histIndexCSC = {{{1, 1}, 8},
210  {{1, 2}, 7},
211  {{1, 3}, 6},
212  {{2, 1}, 5},
213  {{2, 2}, 4},
214  {{3, 1}, 3},
215  {{3, 2}, 2},
216  {{4, 1}, 1},
217  {{4, 2}, 0}};
218 
219  const int min_endcap = CSCDetId::minEndcapId();
220  const int max_endcap = CSCDetId::maxEndcapId();
221  const int min_station = CSCDetId::minStationId();
222  const int max_station = CSCDetId::maxStationId();
223  const int min_sector = CSCTriggerNumbering::minTriggerSectorId();
224  const int max_sector = CSCTriggerNumbering::maxTriggerSectorId();
225  const int min_subsector = CSCTriggerNumbering::minTriggerSubSectorId();
226  const int max_subsector = CSCTriggerNumbering::maxTriggerSubSectorId();
227  const int min_chamber = CSCTriggerNumbering::minTriggerCscId();
228  const int max_chamber = CSCTriggerNumbering::maxTriggerCscId();
229 
230  for (int endc = min_endcap; endc <= max_endcap; endc++) {
231  // loop on all stations
232  for (int stat = min_station; stat <= max_station; stat++) {
233  int numsubs = ((stat == 1) ? max_subsector : 1);
234  // loop on sectors and subsectors
235  for (int sect = min_sector; sect <= max_sector; sect++) {
236  for (int subs = min_subsector; subs <= numsubs; subs++) {
237  // loop on all chambers
238  for (int cham = min_chamber; cham <= max_chamber; cham++) {
239  // extract the ring number
240  int ring = CSCTriggerNumbering::ringFromTriggerLabels(stat, cham);
241 
242  // actual chamber number =/= trigger chamber number
243  int chid = CSCTriggerNumbering::chamberFromTriggerLabels(sect, subs, stat, cham);
244 
245  // 0th layer means whole chamber.
246  CSCDetId detid(endc, stat, ring, chid, 0);
247 
248  int chamber = detid.chamber();
249 
250  int sr = histIndexCSC.at({stat, ring});
251  if (endc == 1)
252  sr = 17 - sr;
253  bool chamber20 = (sr == 1 or sr == 3 or sr == 5 or sr == 12 or sr == 14 or sr == 16);
254 
255  // ALCT analysis
256  auto range_dataALCT = dataALCTshs->get(detid);
257  auto range_emulALCT = emulALCTshs->get(detid);
258 
259  for (auto dalct = range_dataALCT.first; dalct != range_dataALCT.second; dalct++) {
260  if (dalct->isValid() and dalct->isNominalInTime()) {
261  if (dalct->isTightInTime()) {
262  if (chamber20) {
263  alctShowerDataTightSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
264  alctShowerDataTightSummary_denom_->Fill(chamber * 2, sr, 0.5);
265  } else
266  alctShowerDataTightSummary_denom_->Fill(chamber, sr);
267  }
268  if (chamber20) {
269  alctShowerDataNomSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
270  alctShowerDataNomSummary_denom_->Fill(chamber * 2, sr, 0.5);
271  } else
272  alctShowerDataNomSummary_denom_->Fill(chamber, sr);
273  // check for least one matching ALCT
274  for (auto ealct = range_emulALCT.first; ealct != range_emulALCT.second; ealct++) {
275  if (ealct->isValid() and areSameShowers(*dalct, *ealct)) {
276  if (dalct->isTightInTime()) {
277  if (chamber20) {
278  alctShowerDataTightSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
279  alctShowerDataTightSummary_num_->Fill(chamber * 2, sr, 0.5);
280  } else
281  alctShowerDataTightSummary_num_->Fill(chamber, sr);
282  }
283  if (chamber20) {
284  alctShowerDataNomSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
285  alctShowerDataNomSummary_num_->Fill(chamber * 2, sr, 0.5);
286  } else
287  alctShowerDataNomSummary_num_->Fill(chamber, sr);
288  }
289  }
290  }
291  } // End of for (auto dalct = range_dataALCT.first; dalct != range_dataALCT.second; dalct++)
292 
293  for (auto ealct = range_emulALCT.first; ealct != range_emulALCT.second; ealct++) {
294  bool isMatched = false;
295  if (ealct->isValid() and ealct->isNominalInTime()) {
296  if (ealct->isTightInTime()) {
297  if (chamber20) {
298  alctShowerEmulTightSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
299  alctShowerEmulTightSummary_denom_->Fill(chamber * 2, sr, 0.5);
300  } else
301  alctShowerEmulTightSummary_denom_->Fill(chamber, sr);
302  }
303  if (chamber20) {
304  alctShowerEmulNomSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
305  alctShowerEmulNomSummary_denom_->Fill(chamber * 2, sr, 0.5);
306  } else
307  alctShowerEmulNomSummary_denom_->Fill(chamber, sr);
308  // check for least one matching ALCT
309  for (auto dalct = range_dataALCT.first; dalct != range_dataALCT.second; dalct++) {
310  if (areSameShowers(*dalct, *ealct))
311  isMatched = true;
312  }
313  // only fill when it is not matched to an ALCT
314  // to understand if the emulator is producing too many ALCTs
315  if (!isMatched) {
316  if (ealct->isTightInTime()) {
317  if (chamber20) {
318  alctShowerEmulTightSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
319  alctShowerEmulTightSummary_num_->Fill(chamber * 2, sr, 0.5);
320  } else
321  alctShowerEmulTightSummary_num_->Fill(chamber, sr);
322  }
323  if (chamber20) {
324  alctShowerEmulNomSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
325  alctShowerEmulNomSummary_num_->Fill(chamber * 2, sr, 0.5);
326  } else
327  alctShowerEmulNomSummary_num_->Fill(chamber, sr);
328  }
329  }
330  } // End of for (auto ealct = range_emulALCT.first; ealct != range_emulALCT.second; ealct++)
331 
332  // CLCT analysis
333  auto range_dataCLCT = dataCLCTshs->get(detid);
334  auto range_emulCLCT = emulCLCTshs->get(detid);
335 
336  for (auto dclct = range_dataCLCT.first; dclct != range_dataCLCT.second; dclct++) {
337  if (dclct->isValid() and dclct->isNominalInTime()) {
338  if (dclct->isTightInTime()) {
339  if (chamber20) {
340  clctShowerDataTightSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
341  clctShowerDataTightSummary_denom_->Fill(chamber * 2, sr, 0.5);
342  } else
343  clctShowerDataTightSummary_denom_->Fill(chamber, sr);
344  }
345  if (chamber20) {
346  clctShowerDataNomSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
347  clctShowerDataNomSummary_denom_->Fill(chamber * 2, sr, 0.5);
348  } else
349  clctShowerDataNomSummary_denom_->Fill(chamber, sr);
350  // check for least one matching CLCT
351  for (auto eclct = range_emulCLCT.first; eclct != range_emulCLCT.second; eclct++) {
352  if (eclct->isValid() and areSameShowers(*dclct, *eclct)) {
353  if (dclct->isTightInTime()) {
354  if (chamber20) {
355  clctShowerDataTightSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
356  clctShowerDataTightSummary_num_->Fill(chamber * 2, sr, 0.5);
357  } else
358  clctShowerDataTightSummary_num_->Fill(chamber, sr);
359  }
360  if (chamber20) {
361  clctShowerDataNomSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
362  clctShowerDataNomSummary_num_->Fill(chamber * 2, sr, 0.5);
363  } else
364  clctShowerDataNomSummary_num_->Fill(chamber, sr);
365  }
366  }
367  }
368  } // End of for (auto dclct = range_dataCLCT.first; dclct != range_dataCLCT.second; dclct++)
369 
370  for (auto eclct = range_emulCLCT.first; eclct != range_emulCLCT.second; eclct++) {
371  bool isMatched = false;
372  if (eclct->isValid() and eclct->isNominalInTime()) {
373  if (eclct->isTightInTime()) {
374  if (chamber20) {
375  clctShowerEmulTightSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
376  clctShowerEmulTightSummary_denom_->Fill(chamber * 2, sr, 0.5);
377  } else
378  clctShowerEmulTightSummary_denom_->Fill(chamber, sr);
379  }
380  if (chamber20) {
381  clctShowerEmulNomSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
382  clctShowerEmulNomSummary_denom_->Fill(chamber * 2, sr, 0.5);
383  } else
384  clctShowerEmulNomSummary_denom_->Fill(chamber, sr);
385  // check for least one matching CLCT
386  for (auto dclct = range_dataCLCT.first; dclct != range_dataCLCT.second; dclct++) {
387  if (areSameShowers(*dclct, *eclct))
388  isMatched = true;
389  }
390  // only fill when it is not matched to an CLCT
391  // to understand if the emulator is producing too many CLCTs
392  if (!isMatched) {
393  if (eclct->isTightInTime()) {
394  if (chamber20) {
395  clctShowerEmulTightSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
396  clctShowerEmulTightSummary_num_->Fill(chamber * 2, sr, 0.5);
397  } else
398  clctShowerEmulTightSummary_num_->Fill(chamber, sr);
399  }
400  if (chamber20) {
401  clctShowerEmulNomSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
402  clctShowerEmulNomSummary_num_->Fill(chamber * 2, sr, 0.5);
403  } else
404  clctShowerEmulNomSummary_num_->Fill(chamber, sr);
405  }
406  }
407  } // End of for (auto eclct = range_emulCLCT.first; eclct != range_emulCLCT.second; eclct++)
408 
409  // LCT analysis
410  auto range_dataLCT = dataLCTshs->get(detid);
411  auto range_emulLCT = emulLCTshs->get(detid);
412 
413  for (auto dlct = range_dataLCT.first; dlct != range_dataLCT.second; dlct++) {
414  if (dlct->isValid() and dlct->isNominalInTime()) {
415  if (dlct->isTightInTime()) {
416  if (chamber20) {
417  lctShowerDataTightSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
418  lctShowerDataTightSummary_denom_->Fill(chamber * 2, sr, 0.5);
419  } else
420  lctShowerDataTightSummary_denom_->Fill(chamber, sr);
421  }
422  if (chamber20) {
423  lctShowerDataNomSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
424  lctShowerDataNomSummary_denom_->Fill(chamber * 2, sr, 0.5);
425  } else
426  lctShowerDataNomSummary_denom_->Fill(chamber, sr);
427  // check for least one matching LCT
428  for (auto elct = range_emulLCT.first; elct != range_emulLCT.second; elct++) {
429  if (elct->isValid() and areSameShowers(*dlct, *elct)) {
430  if (dlct->isTightInTime()) {
431  if (chamber20) {
432  lctShowerDataTightSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
433  lctShowerDataTightSummary_num_->Fill(chamber * 2, sr, 0.5);
434  } else
435  lctShowerDataTightSummary_num_->Fill(chamber, sr);
436  }
437  if (chamber20) {
438  lctShowerDataNomSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
439  lctShowerDataNomSummary_num_->Fill(chamber * 2, sr, 0.5);
440  } else
441  lctShowerDataNomSummary_num_->Fill(chamber, sr);
442  }
443  }
444  }
445  } // End of for (auto dlct = range_dataLCT.first; dlct != range_dataLCT.second; dlct++)
446 
447  for (auto elct = range_emulLCT.first; elct != range_emulLCT.second; elct++) {
448  bool isMatched = false;
449  if (elct->isValid() and elct->isNominalInTime()) {
450  if (elct->isTightInTime()) {
451  if (chamber20) {
452  lctShowerEmulTightSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
453  lctShowerEmulTightSummary_denom_->Fill(chamber * 2, sr, 0.5);
454  } else
455  lctShowerEmulTightSummary_denom_->Fill(chamber, sr);
456  }
457  if (chamber20) {
458  lctShowerEmulNomSummary_denom_->Fill(chamber * 2 - 1, sr, 0.5);
459  lctShowerEmulNomSummary_denom_->Fill(chamber * 2, sr, 0.5);
460  } else
461  lctShowerEmulNomSummary_denom_->Fill(chamber, sr);
462  // check for least one matching LCT
463  for (auto dlct = range_dataLCT.first; dlct != range_dataLCT.second; dlct++) {
464  if (areSameShowers(*dlct, *elct))
465  isMatched = true;
466  }
467  // only fill when it is not matched to an LCT
468  // to understand if the emulator is producing too many LCTs
469  if (!isMatched) {
470  if (elct->isTightInTime()) {
471  if (chamber20) {
472  lctShowerEmulTightSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
473  lctShowerEmulTightSummary_num_->Fill(chamber * 2, sr, 0.5);
474  } else
475  lctShowerEmulTightSummary_num_->Fill(chamber, sr);
476  }
477  if (chamber20) {
478  lctShowerEmulNomSummary_num_->Fill(chamber * 2 - 1, sr, 0.5);
479  lctShowerEmulNomSummary_num_->Fill(chamber * 2, sr, 0.5);
480  } else
481  lctShowerEmulNomSummary_num_->Fill(chamber, sr);
482  }
483  }
484  } // End of for (auto elct = range_emulLCT.first; elct != range_emulLCT.second; elct++) {
485  }
486  }
487  }
488  }
489  }
490 }
491 
492 bool L1TdeCSCTPGShower::areSameShowers(const CSCShowerDigi& lhs, const CSCShowerDigi& rhs) const {
493  bool returnValue = false;
494  if (lhs.isValid() == rhs.isValid() && lhs.getCSCID() == rhs.getCSCID() && lhs.bitsInTime() == rhs.bitsInTime() &&
495  lhs.bitsOutOfTime() == rhs.bitsOutOfTime()) {
496  returnValue = true;
497  }
498  return returnValue;
499 }
L1TdeCSCTPGShower::lctShowerEmulNomSummary_num_
MonitorElement * lctShowerEmulNomSummary_num_
Definition: L1TdeCSCTPGShower.h:42
L1TdeCSCTPGShower::lctShowerEmulNomSummary_denom_
MonitorElement * lctShowerEmulNomSummary_denom_
Definition: L1TdeCSCTPGShower.h:41
L1TdeCSCTPGShower::alctShowerDataNomSummary_num_
MonitorElement * alctShowerDataNomSummary_num_
Definition: L1TdeCSCTPGShower.h:37
dqm::implementation::NavigatorBase::setCurrentFolder
virtual void setCurrentFolder(std::string const &fullpath)
Definition: DQMStore.cc:36
CSCDetId::minEndcapId
static int minEndcapId()
Definition: CSCDetId.h:234
L1TdeCSCTPGShower::clctShowerDataNomSummary_num_
MonitorElement * clctShowerDataNomSummary_num_
Definition: L1TdeCSCTPGShower.h:39
L1TdeCSCTPGShower::clctShowerDataTightSummary_denom_
MonitorElement * clctShowerDataTightSummary_denom_
Definition: L1TdeCSCTPGShower.h:52
L1TdeCSCTPGShower::dataLCTShower_token_
edm::EDGetTokenT< CSCShowerDigiCollection > dataLCTShower_token_
Definition: L1TdeCSCTPGShower.h:29
dqm::impl::MonitorElement::setOption
virtual void setOption(const char *option)
Definition: MonitorElement.cc:888
CSCTriggerNumbering::maxTriggerSectorId
static int maxTriggerSectorId()
Definition: CSCTriggerNumbering.h:112
CSCShowerDigi::isValid
bool isValid() const
data
Definition: CSCShowerDigi.cc:38
L1TdeCSCTPGShower::dataCLCTShower_token_
edm::EDGetTokenT< CSCShowerDigiCollection > dataCLCTShower_token_
Definition: L1TdeCSCTPGShower.h:27
HLT_2022v15_cff.InputTag
InputTag
Definition: HLT_2022v15_cff.py:60425
L1TdeCSCTPGShower::emulCLCTShower_token_
edm::EDGetTokenT< CSCShowerDigiCollection > emulCLCTShower_token_
Definition: L1TdeCSCTPGShower.h:28
MillePedeFileConverter_cfg.e
e
Definition: MillePedeFileConverter_cfg.py:37
CSCShowerDigi::bitsOutOfTime
uint16_t bitsOutOfTime() const
Definition: CSCShowerDigi.h:49
L1TdeCSCTPGShower::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: L1TdeCSCTPGShower.cc:193
CSCDetId::maxStationId
static int maxStationId()
Definition: CSCDetId.h:237
CSCTriggerNumbering::minTriggerSectorId
static int minTriggerSectorId()
Definition: CSCTriggerNumbering.h:113
HltBtagPostValidation_cff.c
c
Definition: HltBtagPostValidation_cff.py:31
CSCTriggerNumbering::ringFromTriggerLabels
static int ringFromTriggerLabels(int station, int triggerCSCID)
Definition: CSCTriggerNumbering.cc:5
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:256
dqm::impl::MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:292
L1TdeCSCTPGShower::lctShowerDataTightSummary_num_
MonitorElement * lctShowerDataTightSummary_num_
Definition: L1TdeCSCTPGShower.h:49
L1TdeCSCTPGShower::lctShowerEmulTightSummary_num_
MonitorElement * lctShowerEmulTightSummary_num_
Definition: L1TdeCSCTPGShower.h:56
L1TdeCSCTPGShower::L1TdeCSCTPGShower
L1TdeCSCTPGShower(const edm::ParameterSet &ps)
Definition: L1TdeCSCTPGShower.cc:6
CSCDetId::minStationId
static int minStationId()
Definition: CSCDetId.h:236
L1TdeCSCTPGShower::bookHistograms
void bookHistograms(DQMStore::IBooker &, const edm::Run &, const edm::EventSetup &) override
Definition: L1TdeCSCTPGShower.cc:17
or
The Signals That Services Can Subscribe To This is based on ActivityRegistry and is current per Services can connect to the signals distributed by the ActivityRegistry in order to monitor the activity of the application Each possible callback has some defined which we here list in angle e< void, edm::EventID const &, edm::Timestamp const & > We also list in braces which AR_WATCH_USING_METHOD_ is used for those or
Definition: Activities.doc:12
CSCDetId::maxEndcapId
static int maxEndcapId()
Definition: CSCDetId.h:235
L1TdeCSCTPGShower::alctShowerEmulNomSummary_denom_
MonitorElement * alctShowerEmulNomSummary_denom_
Definition: L1TdeCSCTPGShower.h:43
L1TdeCSCTPGShower::alctShowerDataTightSummary_denom_
MonitorElement * alctShowerDataTightSummary_denom_
Definition: L1TdeCSCTPGShower.h:50
edm_modernize_messagelogger.stat
stat
Definition: edm_modernize_messagelogger.py:27
MuonDigiCollection::get
Range get(const IndexType &index) const
return the digis for a given DetUnit
Definition: MuonDigiCollection.h:121
CSCDetId::chamber
int chamber() const
Definition: CSCDetId.h:62
trackerHitRTTI::isMatched
bool isMatched(TrackingRecHit const &hit)
Definition: trackerHitRTTI.h:33
L1TdeCSCTPGShower::lctShowerEmulTightSummary_denom_
MonitorElement * lctShowerEmulTightSummary_denom_
Definition: L1TdeCSCTPGShower.h:55
L1TdeCSCTPGShower::clctShowerDataNomSummary_denom_
MonitorElement * clctShowerDataNomSummary_denom_
Definition: L1TdeCSCTPGShower.h:38
dqm::impl::MonitorElement::setBinLabel
virtual void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:772
L1TdeCSCTPGShower::emulALCTShower_token_
edm::EDGetTokenT< CSCShowerDigiCollection > emulALCTShower_token_
Definition: L1TdeCSCTPGShower.h:26
L1TdeCSCTPGShower::alctShowerEmulNomSummary_num_
MonitorElement * alctShowerEmulNomSummary_num_
Definition: L1TdeCSCTPGShower.h:44
CSCShowerDigi::getCSCID
uint16_t getCSCID() const
Definition: CSCShowerDigi.h:52
L1TdeCSCTPGShower::alctShowerDataTightSummary_num_
MonitorElement * alctShowerDataTightSummary_num_
Definition: L1TdeCSCTPGShower.h:51
L1TdeCSCTPGShower::clctShowerEmulTightSummary_denom_
MonitorElement * clctShowerEmulTightSummary_denom_
Definition: L1TdeCSCTPGShower.h:59
L1TdeCSCTPGShower::alctShowerEmulTightSummary_num_
MonitorElement * alctShowerEmulTightSummary_num_
Definition: L1TdeCSCTPGShower.h:58
CSCTriggerNumbering::minTriggerSubSectorId
static int minTriggerSubSectorId()
Definition: CSCTriggerNumbering.h:115
L1TdeCSCTPGShower::monitorDir_
std::string monitorDir_
Definition: L1TdeCSCTPGShower.h:32
dqm::implementation::IBooker::book2D
MonitorElement * book2D(TString const &name, TString const &title, int nchX, double lowX, double highX, int nchY, double lowY, double highY, FUNC onbooking=NOOP())
Definition: DQMStore.h:212
L1TdeCSCTPGShower::areSameShowers
bool areSameShowers(const CSCShowerDigi &lhs, const CSCShowerDigi &rhs) const
Definition: L1TdeCSCTPGShower.cc:492
L1TdeCSCTPGShower::clctShowerEmulNomSummary_num_
MonitorElement * clctShowerEmulNomSummary_num_
Definition: L1TdeCSCTPGShower.h:46
CSCShowerDigi::bitsInTime
uint16_t bitsInTime() const
Definition: CSCShowerDigi.h:48
edm
HLT enums.
Definition: AlignableModifier.h:19
L1TdeCSCTPGShower::clctShowerEmulNomSummary_denom_
MonitorElement * clctShowerEmulNomSummary_denom_
Definition: L1TdeCSCTPGShower.h:45
L1TdeCSCTPGShower::lctShowerDataTightSummary_denom_
MonitorElement * lctShowerDataTightSummary_denom_
Definition: L1TdeCSCTPGShower.h:48
L1TdeCSCTPGShower::dataALCTShower_token_
edm::EDGetTokenT< CSCShowerDigiCollection > dataALCTShower_token_
Definition: L1TdeCSCTPGShower.h:25
L1TdeCSCTPGShower::alctShowerDataNomSummary_denom_
MonitorElement * alctShowerDataNomSummary_denom_
Definition: L1TdeCSCTPGShower.h:36
L1TdeCSCTPGShower::lctShowerDataNomSummary_num_
MonitorElement * lctShowerDataNomSummary_num_
Definition: L1TdeCSCTPGShower.h:35
L1TdeCSCTPGShower::lctShowerDataNomSummary_denom_
MonitorElement * lctShowerDataNomSummary_denom_
Definition: L1TdeCSCTPGShower.h:34
L1TdeCSCTPGShower::emulLCTShower_token_
edm::EDGetTokenT< CSCShowerDigiCollection > emulLCTShower_token_
Definition: L1TdeCSCTPGShower.h:30
L1TdeCSCTPGShower::clctShowerDataTightSummary_num_
MonitorElement * clctShowerDataTightSummary_num_
Definition: L1TdeCSCTPGShower.h:53
L1TdeCSCTPGShower::clctShowerEmulTightSummary_num_
MonitorElement * clctShowerEmulTightSummary_num_
Definition: L1TdeCSCTPGShower.h:60
L1TdeCSCTPGShower::~L1TdeCSCTPGShower
~L1TdeCSCTPGShower() override
Definition: L1TdeCSCTPGShower.cc:15
L1TdeCSCTPGShower::alctShowerEmulTightSummary_denom_
MonitorElement * alctShowerEmulTightSummary_denom_
Definition: L1TdeCSCTPGShower.h:57
CSCTriggerNumbering::chamberFromTriggerLabels
static int chamberFromTriggerLabels(int TriggerSector, int TriggerSubSector, int station, int TriggerCSCID)
Definition: CSCTriggerNumbering.cc:30
CSCTriggerNumbering::maxTriggerSubSectorId
static int maxTriggerSubSectorId()
Definition: CSCTriggerNumbering.h:114
edm::Run
Definition: Run.h:45
dqm::impl::MonitorElement::setAxisTitle
virtual void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:801
MuonDigiCollection
A container for a generic type of digis indexed by some index, implemented with a map<IndexType...
Definition: MuonDigiCollection.h:87
Generated for CMSSW Reference Manual by doxygen 1.8.14