CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
TPTask.cc
Go to the documentation of this file.
1 
2 #include "DQM/HcalTasks/interface/TPTask.h"
3 
4 using namespace hcaldqm;
5 TPTask::TPTask(edm::ParameterSet const& ps):
6  DQTask(ps)
7 {
8  // Et
9  _cEtData_SubDet.initialize(_name+"/Et/SubDet_Data", "Et",
10  mapper::fTPSubDet,
11  new axis::ValueAxis(axis::fXaxis, axis::fEt_256),
12  new axis::ValueAxis(axis::fYaxis, axis::fEntries, true));
13  _cEtEmul_SubDet.initialize(_name+"/Et/SubDet_Data", "Et",
14  mapper::fTPSubDet,
15  new axis::ValueAxis(axis::fXaxis, axis::fEt_256),
16  new axis::ValueAxis(axis::fYaxis, axis::fEntries, true));
17  _cEtCorr_TPSubDet.initialize(_name+"/EtCorrelation/TPSubDet",
18  "EtCorrelation",
19  mapper::fTPSubDet,
20  new axis::ValueAxis(axis::fXaxis, axis::fEt_128),
21  new axis::ValueAxis(axis::fYaxis, axis::fEt_128));
22  _cEtCorrRatiovsLS_TPSubDet.initialize(_name+"/EtCorrRatio/vsLS_TPSubDet",
23  "EtCorrelationRatio", mapper::fTPSubDet,
24  new axis::ValueAxis(axis::fXaxis, axis::fLS),
25  new axis::ValueAxis(axis::fYaxis, axis::fRatio));
26  _cEtCorrRatio_TPSubDet.initialize(_name+"/EtCorrRatio/TPSubDet",
27  "EtCorrelationRatio", mapper::fTPSubDet,
28  new axis::ValueAxis(axis::fXaxis, axis::fRatio),
29  new axis::ValueAxis(axis::fYaxis, axis::fEntries, true));
30  _cEtCorrRatiovsiphi_TPSubDetPM.initialize(
31  _name+"/EtCorrRatio/vsiphi_TPSubDetPM",
32  "EtCorrelationRatio", mapper::fTPSubDetPM,
33  new axis::CoordinateAxis(axis::fXaxis, axis::fiphi),
34  new axis::ValueAxis(axis::fXaxis, axis::fRatio));
35  _cEtMsm.initialize(_name+"/EtMisMatched", "EtMismatched",
36  new axis::CoordinateAxis(axis::fXaxis, axis::fTPieta),
37  new axis::CoordinateAxis(axis::fYaxis, axis::fiphi),
38  new axis::ValueAxis(axis::fZaxis, axis::fEntries));
39  _cNumEtMsmvsLS_TPSubDet.initialize(_name+"/EtMismatched/vsLS_TPSubDet",
40  "EtMismatched", mapper::fTPSubDet,
41  new axis::ValueAxis(axis::fXaxis, axis::fLS),
42  new axis::ValueAxis(axis::fYaxis, axis::fEntries100));
43  _cNumEtMsm_TPSubDet.initialize(_name+"/EtMismatched/TPSubDet",
44  "EtMismatched", mapper::fTPSubDet,
45  new axis::ValueAxis(axis::fXaxis, axis::fEntries100),
46  new axis::ValueAxis(axis::fYaxis, axis::fEntries, true));
47  _cNumEtMsmvsiphi_TPSubDetPM.initialize(
48  _name+"/EtMismatched/vsiphi_TPSubDetPM",
49  "EtMismatched", mapper::fTPSubDetPM,
50  new axis::CoordinateAxis(axis::fXaxis, axis::fiphi),
51  new axis::ValueAxis(axis::fYaxis, axis::fEntries));
52  _cSumdEtvsLS_TPSubDet.initialize(_name+"/SumdEt/vsLS_TPSubDet",
53  "SumdEt", mapper::fTPSubDet,
54  new axis::ValueAxis(axis::fXaxis, axis::fLS),
55  new axis::ValueAxis(axis::fYaxis, axis::fSumdEt));
56  _cSumdEt_TPSubDet.initialize(_name+"/SumdEt/TPSubDet",
57  "SumdEt", mapper::fTPSubDet,
58  new axis::ValueAxis(axis::fXaxis, axis::fSumdEt),
59  new axis::ValueAxis(axis::fYaxis, axis::fEntries, true));
60  _cEtData_SubDetPM_iphi.initialize(_name+"/Et/SubDetPM_iphi_Data", "Et",
61  mapper::fTPSubDetPM_iphi,
62  new axis::ValueAxis(axis::fXaxis, axis::fEt_256),
63  new axis::ValueAxis(axis::fYaxis, axis::fEntries, true));
64  _cEtData_SubDet_ieta.initialize(_name+"/Et/SubDet_ieta_Data", "Et",
65  mapper::fTPSubDet_ieta,
66  new axis::ValueAxis(axis::fXaxis, axis::fEt_256),
67  new axis::ValueAxis(axis::fYaxis, axis::fEntries, true));
68 
69  // FG
70  _cFGCorr_SubDet.initialize(_name+"/FGCorrelation/SubDet", "FG",
71  mapper::fTPSubDet,
72  new axis::ValueAxis(axis::fXaxis, axis::fFG),
73  new axis::ValueAxis(axis::fYaxis, axis::fFG));
74  _cFGMsm.initialize(_name+"/FGMisMatched", "FGMismatched",
75  new axis::CoordinateAxis(axis::fXaxis, axis::fTPieta),
76  new axis::CoordinateAxis(axis::fYaxis, axis::fiphi),
77  new axis::ValueAxis(axis::fZaxis, axis::fEntries));
78 
79  // Occupancy
80  _cOccupancyData.initialize(_name+"/Occupancy/Data", "Occupancy",
81  new axis::CoordinateAxis(axis::fXaxis, axis::fTPieta),
82  new axis::CoordinateAxis(axis::fYaxis, axis::fiphi),
83  new axis::ValueAxis(axis::fZaxis, axis::fEntries));
84  _cOccupancyEmul.initialize(_name+"/Occupancy/Emul", "Occupancy",
85  new axis::CoordinateAxis(axis::fXaxis, axis::fTPieta),
86  new axis::CoordinateAxis(axis::fYaxis, axis::fiphi),
87  new axis::ValueAxis(axis::fZaxis, axis::fEntries));
88  _cOccupancyDatavsiphi_TPSubDetPM.initialize(
89  _name+"/Occupancy/vsiphi_TPSubDetPM_Data",
90  "Occupancy", mapper::fTPSubDetPM,
91  new axis::CoordinateAxis(axis::fXaxis, axis::fiphi));
92  _cOccupancyEmulvsiphi_TPSubDetPM.initialize(
93  _name+"/Occupancy/vsiphi_TPSubDetPM_Emul",
94  "Occupancy", mapper::fTPSubDetPM,
95  new axis::CoordinateAxis(axis::fXaxis, axis::fiphi));
96  _cOccupancyDatavsLS_TPSubDet.initialize(
97  _name+"/Occupancy/vsLS_TPSubDet_Data",
98  "Occupancy", mapper::fTPSubDet,
99  new axis::ValueAxis(axis::fXaxis, axis::fLS),
100  new axis::ValueAxis(axis::fYaxis, axis::fEntries));
101  _cOccupancyEmulvsLS_TPSubDet.initialize(
102  _name+"/Occupancy/vsLS_TPSubDet_Emul",
103  "Occupancy", mapper::fTPSubDet,
104  new axis::ValueAxis(axis::fXaxis, axis::fLS),
105  new axis::ValueAxis(axis::fYaxis, axis::fEntries));
106  _cMsData.initialize(_name+"/Missing/Data", "Missing",
107  new axis::CoordinateAxis(axis::fXaxis, axis::fTPieta),
108  new axis::CoordinateAxis(axis::fYaxis, axis::fiphi),
109  new axis::ValueAxis(axis::fZaxis, axis::fEntries));
110  _cMsDatavsiphi_TPSubDetPM.initialize(
111  _name+"/Missing/vsiphi_TPSubDetPM_Data",
112  "Missing", mapper::fTPSubDetPM,
113  new axis::CoordinateAxis(axis::fXaxis, axis::fiphi),
114  new axis::ValueAxis(axis::fYaxis, axis::fEntries));
115  _cMsEmul.initialize(_name+"/Missing/Emul", "Missing",
116  new axis::CoordinateAxis(axis::fXaxis, axis::fTPieta),
117  new axis::CoordinateAxis(axis::fYaxis, axis::fiphi),
118  new axis::ValueAxis(axis::fZaxis, axis::fEntries));
119 
120  // Digi Sizes
121  _cDigiSizeDatavsLS_TPSubDet.initialize(_name+"/DigiSize/vsLS_TPSubDet_Data",
122  "DigiSize", mapper::fTPSubDet,
123  new axis::ValueAxis(axis::fXaxis, axis::fLS),
124  new axis::ValueAxis(axis::fYaxis, axis::fDigiSize));
125  _cDigiSizeEmulvsLS_TPSubDet.initialize(_name+"/DigiSize/vsLS_TPSubDet_Emul",
126  "DigiSize", mapper::fTPSubDet,
127  new axis::ValueAxis(axis::fXaxis, axis::fLS),
128  new axis::ValueAxis(axis::fYaxis, axis::fDigiSize));
129 
130  // Summary
131  _cSummary.initialize(_name+"/Summary", "Summary",
132  new axis::CoordinateAxis(axis::fXaxis, axis::fTPSubDet),
133  new axis::FlagAxis(axis::fYaxis, "Flag", int(nTPFlag)));
134  _cSummaryvsLS_TPSubDet.initialize(_name+"/Summary/vsLS_TPSubDet", "Summary",
135  mapper::fTPSubDet,
136  new axis::ValueAxis(axis::fXaxis, axis::fLS),
137  new axis::FlagAxis(axis::fYaxis, "Flag", int(nTPFlag)));
138 
139  // tags
140  _tagData = ps.getUntrackedParameter<edm::InputTag>("tagData",
141  edm::InputTag("hcalDigis"));
142  _tagEmul = ps.getUntrackedParameter<edm::InputTag>("tagEmul",
143  edm::InputTag("emulDigis"));
144  _tokData = consumes<HcalTrigPrimDigiCollection>(_tagData);
145  _tokEmul = consumes<HcalTrigPrimDigiCollection>(_tagEmul);
146 
147  // Special
148  _skip1x1 = ps.getUntrackedParameter<bool>("skip1x1", true);
149 
150  // flags
151  _fNames.push_back("Occ. iphi Uniform. D");
152  _fNames.push_back("Occ. iphi Uniform. E");
153  _fNames.push_back("Low Occupancy E");
154  _fNames.push_back("Corr. Ratio");
155  _fNames.push_back("Corr. Ratio iphi Uniform");
156  _fNames.push_back("Msm. Et iphi Uniform. ");
157  _fNames.push_back("Msm. Et Number");
158  _fNames.push_back("Missing iphi Uniform. D");
159  _cSummary.loadLabels(_fNames);
160  _cSummaryvsLS_TPSubDet.loadLabels(_fNames);
161 }
162 
163 /* virtual */ void TPTask::bookHistograms(DQMStore::IBooker &ib,
164  edm::Run const& r, edm::EventSetup const& es)
165 {
166  DQTask::bookHistograms(ib, r, es);
167  _cEtData_SubDet.book(ib);
168  _cEtEmul_SubDet.book(ib);
169  _cEtCorr_TPSubDet.book(ib);
170  _cEtCorrRatiovsLS_TPSubDet.book(ib);
171  _cEtCorrRatio_TPSubDet.book(ib);
172  _cEtCorrRatiovsiphi_TPSubDetPM.book(ib);
173  _cEtMsm.book(ib);
174  _cSumdEtvsLS_TPSubDet.book(ib);
175  _cSumdEt_TPSubDet.book(ib);
176  _cNumEtMsmvsLS_TPSubDet.book(ib);
177  _cNumEtMsm_TPSubDet.book(ib);
178  _cNumEtMsmvsiphi_TPSubDetPM.book(ib);
179  _cEtData_SubDetPM_iphi.book(ib);
180  _cEtData_SubDet_ieta.book(ib);
181  _cFGCorr_SubDet.book(ib);
182  _cFGMsm.book(ib);
183  _cOccupancyData.book(ib);
184  _cOccupancyEmul.book(ib);
185  _cOccupancyDatavsiphi_TPSubDetPM.book(ib);
186  _cOccupancyEmulvsiphi_TPSubDetPM.book(ib);
187  _cOccupancyDatavsLS_TPSubDet.book(ib);
188  _cOccupancyEmulvsLS_TPSubDet.book(ib);
189  _cMsData.book(ib);
190  _cMsEmul.book(ib);
191  _cMsDatavsiphi_TPSubDetPM.book(ib);
192 
193  _cDigiSizeDatavsLS_TPSubDet.book(ib);
194  _cDigiSizeEmulvsLS_TPSubDet.book(ib);
195 
196  _cSummary.book(ib);
197  _cSummaryvsLS_TPSubDet.book(ib);
198 }
199 
200 /* virtual */ void TPTask::_process(edm::Event const& e,
201  edm::EventSetup const &es)
202 {
203  edm::Handle<HcalTrigPrimDigiCollection> ctpd;
204  edm::Handle<HcalTrigPrimDigiCollection> ctpe;
205 
206  if (!e.getByToken(_tokData, ctpd))
207  _logger.dqmthrow("Collection HcalTrigPrimDigiCollection isn't available"
208  + _tagData.label() + " " + _tagData.instance());
209  if (!e.getByToken(_tokEmul, ctpe))
210  _logger.dqmthrow("Collection HcalTrigPrimDigiCollection isn't available"
211  + _tagEmul.label() + " " + _tagEmul.instance());
212 
213  // tmp
214  bool useD1 = false;
215  // tmp
216 
217 
218  // iterate thru the Data Digis
219  double sumdEt[2] = {0, 0};
220  for (HcalTrigPrimDigiCollection::const_iterator hddigi=ctpd->begin();
221  hddigi!=ctpd->end(); ++hddigi)
222  {
223  // tmpt
224  if (_skip1x1)
225  if (hddigi->id().depth()==10) // 10 is the depth for 1x1 TPs
226  continue;
227  //\tmp
228 
229  // get quantities
230  int soiEt_d = hddigi->SOI_compressedEt();
231  int soiFG_d = hddigi->SOI_fineGrain() ? 1 : 0;
232  HcalTrigTowerDetId tid = hddigi->id();
233  // NOTE: skip1x1 should've done that - in the emap it's depth10
234  // in data in 80X it's 0 for some reason!!!
235  // skip using version number
236  if (tid.version()>0)
237  continue;
238 
239  // tmp
240  if (hddigi->id().depth()==1)
241  useD1 = true;
242  //\tmp
243 
244  // fill individual
245  _cEtData_SubDet.fill(tid, soiEt_d);
246  _cEtData_SubDetPM_iphi.fill(tid, soiEt_d);
247  _cEtData_SubDet_ieta.fill(tid, soiEt_d);
248  _cOccupancyData.fill(tid);
249  _cOccupancyDatavsiphi_TPSubDetPM.fill(tid);
250  _nTPs_Data[utilities::getTPSubDet(tid)]++;
251  _cDigiSizeDatavsLS_TPSubDet.fill(tid, _currentLS, hddigi->size());
252 
253  // now, find the emulator digi and compare
254  HcalTrigPrimDigiCollection::const_iterator hedigi =
255  ctpe->find(HcalTrigTowerDetId(tid.ieta(),
256  tid.iphi(), 0));
257 
258  if (hedigi!=ctpe->end())
259  {
260  int soiEt_e = hedigi->SOI_compressedEt();
261  int soiFG_e = hedigi->SOI_fineGrain() ? 1 : 0;
262  double rEt = double(std::min(soiEt_d, soiEt_e))/
263  double(std::max(soiEt_e, soiEt_d));
264  sumdEt[utilities::getTPSubDet(tid)]+=std::abs(soiEt_e-soiEt_d);
265 
266  // fill correlations
267  _cEtEmul_SubDet.fill(tid, soiEt_e);
268  _cEtCorr_TPSubDet.fill(tid, soiEt_e, soiEt_d);
269  _cFGCorr_SubDet.fill(tid, soiFG_e, soiFG_d);
270  _cOccupancyEmul.fill(tid);
271  _cEtCorrRatiovsLS_TPSubDet.fill(tid, _currentLS, rEt);
272  _cEtCorrRatio_TPSubDet.fill(tid, rEt);
273  _cOccupancyEmulvsiphi_TPSubDetPM.fill(tid);
274  _nTPs_Emul[utilities::getTPSubDet(tid)]++;
275  _cEtCorrRatiovsiphi_TPSubDetPM.fill(tid, rEt);
276  _cDigiSizeEmulvsLS_TPSubDet.fill(tid, _currentLS, hedigi->size());
277 
278  // fill those that mismatched
279  if (soiEt_d!=soiEt_e)
280  {
281  _cEtMsm.fill(tid);
282  _nMsmEt[utilities::getTPSubDet(tid)]++;
283  _cNumEtMsmvsiphi_TPSubDetPM.fill(tid);
284  }
285  if (soiFG_d!=soiEt_e)
286  {
287  _cFGMsm.fill(hddigi->id());
288  _nMsmFG[utilities::getTPSubDet(tid)]++;
289  }
290  }
291  else
292  {
293  _cMsEmul.fill(hddigi->id());
294  }
295  }
296 
297  // iterate thru the Emulator Digis
298  for (HcalTrigPrimDigiCollection::const_iterator hedigi=ctpe->begin();
299  hedigi!=ctpe->end(); ++hedigi)
300  {
301  // tmp
302  if (_skip1x1)
303  if (hedigi->id().depth()==10) // 10 for 1x1
304  continue;
305  // \tmp
306 
307  // now, find the emulator digi and compare
308  // tmp
309  HcalTrigPrimDigiCollection::const_iterator hddigi =
310  ctpd->find(HcalTrigTowerDetId(hedigi->id().ieta(),
311  hedigi->id().iphi(), useD1 ? 1: 0));
312  // \tmp
313 
314  // we have to only check those that are missing
315  if (hddigi==ctpd->end())
316  {
317  _cMsData.fill(hedigi->id());
318  _cMsDatavsiphi_TPSubDetPM.fill(hedigi->id());
319  _cOccupancyEmulvsiphi_TPSubDetPM.fill(hedigi->id());
320  _nTPs_Emul[utilities::getTPSubDet(hedigi->id())]++;
321  _cDigiSizeEmulvsLS_TPSubDet.fill(hedigi->id(), _currentLS,
322  hedigi->size());
323  }
324  }
325 
326  // fill per event quantities
327  _cOccupancyDatavsLS_TPSubDet.fill(HcalTrigTowerDetId(5, 5),
328  _currentLS, _nTPs_Data[0]);
329  _cOccupancyDatavsLS_TPSubDet.fill(HcalTrigTowerDetId(30, 5),
330  _currentLS, _nTPs_Data[1]);
331  _cOccupancyEmulvsLS_TPSubDet.fill(HcalTrigTowerDetId(5, 5),
332  _currentLS, _nTPs_Emul[0]);
333  _cOccupancyEmulvsLS_TPSubDet.fill(HcalTrigTowerDetId(30, 5),
334  _currentLS, _nTPs_Emul[1]);
335  _cNumEtMsmvsLS_TPSubDet.fill(HcalTrigTowerDetId(5, 5),
336  _currentLS, _nMsmEt[0]);
337  _cNumEtMsmvsLS_TPSubDet.fill(HcalTrigTowerDetId(30, 5),
338  _currentLS, _nMsmEt[1]);
339  _cNumEtMsm_TPSubDet.fill(HcalTrigTowerDetId(5, 5), _nMsmEt[0]);
340  _cNumEtMsm_TPSubDet.fill(HcalTrigTowerDetId(30, 5), _nMsmEt[1]);
341  _cSumdEtvsLS_TPSubDet.fill(HcalTrigTowerDetId(5, 5), _currentLS,
342  sumdEt[0]);
343  _cSumdEtvsLS_TPSubDet.fill(HcalTrigTowerDetId(30, 5), _currentLS,
344  sumdEt[1]);
345  _cSumdEt_TPSubDet.fill(HcalTrigTowerDetId(5, 5),
346  sumdEt[0]);
347  _cSumdEt_TPSubDet.fill(HcalTrigTowerDetId(30, 5),
348  sumdEt[1]);
349 }
350 
351 /* virtual */ void TPTask::endLuminosityBlock(edm::LuminosityBlock const& l,
352  edm::EventSetup const& es)
353 {
354  // statuses
355  double status[constants::TPSUBDET_NUM][nTPFlag];
356  for (unsigned int i=0; i<constants::TPSUBDET_NUM; i++)
357  for (int j=fOccUniphi_Data; j<nTPFlag; j++)
358  status[i][j] = constants::NOT_APPLICABLE;
359 
360  /*
361  * Do the checks here
362  * -> iphi Uniformity of the Data HF Occupancy
363  * -> iphi Uniformity of the Emul HF Occupancy
364  * -> Low Occupancy for HF Emulator (Data is ZSs)
365  * -> Correlation Ratio vs LS check
366  * -> Correlation Ratio Uniformity
367  * -> iphi Uniformity for the number of mismatches
368  * -> iphi Uniformity for Missing Data
369  * -> iphi Uniformity for Missing Emul
370  */
371 
372  // Low Occp for HF Emul
373  MonitorElement *meocpHF = _cOccupancyEmulvsLS_TPSubDet.at(1);
374  double numChs = meocpHF->getBinContent(_currentLS);
375  if (constants::TPCHS_NUM[1] - numChs>=4)
376  status[1][fLowOcp_Emul] = constants::VERY_LOW;
377  else if (constants::TPCHS_NUM[1] - numChs>=1)
378  status[1][fLowOcp_Emul] = constants::LOW;
379  else if (constants::TPCHS_NUM[1]-numChs<0)
380  status[1][fLowOcp_Emul] = constants::LOW;
381  else if (constants::TPCHS_NUM[1]==numChs)
382  status[1][fLowOcp_Emul] = constants::GOOD;
383 
384  // iphi Occupancy Uniformity for Data/Emul HF
385  bool nonUni_d = false;
386  bool nonUni_e = false;
387  bool nonUniMsm = false;
388  bool nonUniMsn = false;
389  for (int i=1; i<=69; i+=8)
390  {
391  int i1= i;
392  int i2 = (i+4)%IPHI_NUM;
393 
394  double occ1_m_d = _cOccupancyDatavsiphi_TPSubDetPM.getBinContent(2, i1);
395  double occ2_m_d = _cOccupancyDatavsiphi_TPSubDetPM.getBinContent(2, i2);
396  double ratio_m_d = std::min(occ1_m_d, occ2_m_d)/
397  std::max(occ1_m_d, occ2_m_d);
398  double occ1_p_d = _cOccupancyDatavsiphi_TPSubDetPM.getBinContent(3, i1);
399  double occ2_p_d = _cOccupancyDatavsiphi_TPSubDetPM.getBinContent(3, i2);
400  double ratio_p_d = std::min(occ1_p_d, occ2_p_d)/
401  std::max(occ1_p_d, occ2_p_d);
402 
403  double occ1_m_e = _cOccupancyEmulvsiphi_TPSubDetPM.getBinContent(2, i1);
404  double occ2_m_e = _cOccupancyEmulvsiphi_TPSubDetPM.getBinContent(2, i2);
405  double ratio_m_e = std::min(occ1_m_e, occ2_m_e)/
406  std::max(occ1_m_e, occ2_m_e);
407  double occ1_p_e = _cOccupancyEmulvsiphi_TPSubDetPM.getBinContent(3, i1);
408  double occ2_p_e = _cOccupancyEmulvsiphi_TPSubDetPM.getBinContent(3, i2);
409  double ratio_p_e = std::min(occ1_p_e, occ2_p_e)/
410  std::max(occ1_p_e, occ2_p_e);
411 
412  double msm1_m = _cNumEtMsmvsiphi_TPSubDetPM.getBinContent(2, i1);
413  double msm2_m = _cNumEtMsmvsiphi_TPSubDetPM.getBinContent(2, i2);
414  double msmratio_m = std::min(msm1_m, msm2_m)/
415  std::max(msm1_m, msm2_m);
416  double msm1_p = _cNumEtMsmvsiphi_TPSubDetPM.getBinContent(3, i1);
417  double msm2_p = _cNumEtMsmvsiphi_TPSubDetPM.getBinContent(3, i2);
418  double msmratio_p = std::min(msm1_p, msm2_p)/
419  std::max(msm1_p, msm2_p);
420 
421  double msn1_m = _cMsDatavsiphi_TPSubDetPM.getBinContent(2, i1);
422  double msn2_m = _cMsDatavsiphi_TPSubDetPM.getBinContent(2, i2);
423  double msnratio_m = std::min(msn1_m, msn2_m)/
424  std::max(msn1_m, msn2_m);
425  double msn1_p = _cMsDatavsiphi_TPSubDetPM.getBinContent(3, i1);
426  double msn2_p = _cMsDatavsiphi_TPSubDetPM.getBinContent(3, i2);
427  double msnratio_p = std::min(msn1_p, msn2_p)/
428  std::max(msn1_p, msn2_p);
429 
430  if (ratio_m_d<0.8 || ratio_p_d<0.8)
431  nonUni_d = true;
432  else
433  status[1][fOccUniphi_Data] = constants::GOOD;
434  if (ratio_m_e<0.8 || ratio_p_e<0.8)
435  nonUni_e = true;
436  else
437  status[1][fOccUniphi_Emul] = constants::GOOD;
438  if (msmratio_m<0.8 || msmratio_p<0.8)
439  nonUniMsm = true;
440  else
441  status[1][fMsmEtUniphi] = constants::GOOD;
442  if (msnratio_m<0.8 || msnratio_p<0.8)
443  nonUniMsn = true;
444  else
445  status[1][fMsnUniphi_Data] = constants::GOOD;
446  }
447 
448  // set status if any of the flags was set
449  if (nonUni_d)
450  status[1][fOccUniphi_Data] = constants::VERY_LOW;
451  if (nonUni_e)
452  status[1][fOccUniphi_Emul] = constants::VERY_LOW;
453  if (nonUniMsm)
454  status[1][fMsmEtUniphi] = constants::LOW;
455  if (nonUniMsn)
456  status[1][fMsnUniphi_Data] = constants::LOW;
457 
458 
459  // Correlation Ratio
460  double ratio_HBHE = _cEtCorrRatiovsLS_TPSubDet.at(0)
461  ->getBinContent(_currentLS);
462  double ratio_HF = _cEtCorrRatiovsLS_TPSubDet.at(1)
463  ->getBinContent(_currentLS);
464  if (ratio_HBHE<0.7)
465  status[0][fCorrRatio] = constants::LOW;
466  else if (ratio_HBHE<0.8)
467  status[0][fCorrRatio] = constants::PROBLEMATIC;
468  else
469  status[0][fCorrRatio] = constants::GOOD;
470  if (ratio_HF<0.8)
471  status[1][fCorrRatio] = constants::LOW;
472  else if (ratio_HF<0.9)
473  status[1][fCorrRatio] = constants::PROBLEMATIC;
474  else
475  status[1][fCorrRatio] = constants::GOOD;
476 
477  // Number of Mismatches
478  double ratio_nmsm_HBHE = _cNumEtMsmvsLS_TPSubDet.at(0)
479  ->getBinContent(_currentLS)/constants::TPCHS_NUM[0];
480  double ratio_nmsm_HF = _cNumEtMsmvsLS_TPSubDet.at(1)
481  ->getBinContent(_currentLS)/constants::TPCHS_NUM[1];
482  if (ratio_nmsm_HBHE>=0.2)
483  status[0][fMsmEtNum] = constants::VERY_LOW;
484  else if (ratio_nmsm_HBHE>=0.1)
485  status[0][fMsmEtNum] = constants::LOW;
486  else if (ratio_nmsm_HBHE>=0.05)
487  status[0][fMsmEtNum] = constants::PROBLEMATIC;
488  else
489  status[0][fMsmEtNum] = constants::GOOD;
490  if (ratio_nmsm_HF>=0.1)
491  status[1][fMsmEtNum] = constants::LOW;
492  else if (ratio_nmsm_HF>=0.05)
493  status[1][fMsmEtNum] = constants::PROBLEMATIC;
494  else
495  status[1][fMsmEtNum] = constants::GOOD;
496 
497  // finally set all the statuses
498  for (unsigned int i=0; i<constants::TPSUBDET_NUM; i++)
499  for (int j=fOccUniphi_Data; j<nTPFlag; j++)
500  {
501  _cSummary.setBinContent(i, j, status[i][j]);
502  _cSummaryvsLS_TPSubDet.setBinContent(i, _currentLS, j,
503  status[i][j]);
504  }
505 
506  DQTask::endLuminosityBlock(l, es);
507 }
508 
509 /* virtual */ void TPTask::_resetMonitors(UpdateFreq uf)
510 {
511  switch (uf)
512  {
513  case fEvent:
514  for (unsigned int idet=0; idet<constants::TPSUBDET_NUM; idet++)
515  {
516  _nMsmEt[idet] = 0;
517  _nMsmFG[idet] = 0;
518  _nTPs_Data[idet] = 0;
519  _nTPs_Emul[idet] = 0;
520  }
521  break;
522  case hcaldqm::fLS:
523  break;
524  case hcaldqm::f10LS:
525  break;
526  default:
527  break;
528  }
529 
530  DQTask::_resetMonitors(uf);
531 }
532 
533 DEFINE_FWK_MODULE(TPTask);
534 
535 
hcaldqm::constants::LOW
double const LOW
Definition: Constants.h:17
hcaldqm::ContainerProf1D::initialize
virtual void initialize(std::string const &folder, std::string const &nametitle, mapper::MapperType mt, axis::Axis *xaxis, axis::Axis *yaxis=new axis::ValueAxis(axis::fYaxis, axis::fEntries), int debug=0)
Definition: ContainerProf1D.cc:23
TPTask::_cSumdEt_TPSubDet
Container1D _cSumdEt_TPSubDet
Definition: TPTask.h:79
TPTask::_cNumEtMsmvsiphi_TPSubDetPM
Container1D _cNumEtMsmvsiphi_TPSubDetPM
Definition: TPTask.h:77
hcaldqm::ContainerSingle2D::initialize
virtual void initialize(std::string const &folder, std::string const &nametitle, axis::Axis *xaxis, axis::Axis *yaxis=new CoordinateAxis(fYaxis, axis::fiphi), axis::Axis *zaxis=new ValueAxis(fZaxis, fEntries), int debug=0)
Definition: ContainerSingle2D.cc:30
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
TPTask::_skip1x1
bool _skip1x1
Definition: TPTask.h:62
hcaldqm::Container2D::book
virtual void book(DQMStore::IBooker &, std::string subsystem="Hcal", std::string aux="")
Definition: Container2D.cc:233
hcaldqm::constants::IPHI_NUM
int const IPHI_NUM
Definition: Constants.h:100
TPTask::_nMsmEt
int _nMsmEt[constants::TPSUBDET_NUM]
Definition: TPTask.h:56
TPTask::_tagData
edm::InputTag _tagData
Definition: TPTask.h:50
hcaldqm::Container2D::initialize
virtual void initialize(std::string const &folder, std::string nametitle, mapper::MapperType mt, axis::Axis *xaxis, axis::Axis *yaxis=new axis::CoordinateAxis(axis::fYaxis, axis::fiphi), axis::Axis *zaxis=new axis::ValueAxis(axis::fZaxis, axis::fEntries), int debug=0)
Definition: Container2D.cc:30
cuy.ib
int ib
Definition: cuy.py:660
TPTask::_cEtData_SubDetPM_iphi
Container1D _cEtData_SubDetPM_iphi
Definition: TPTask.h:81
hcaldqm::Container1D::book
virtual void book(DQMStore::IBooker &, std::string subsystem="Hcal", std::string aux="")
Definition: Container1D.cc:198
TPTask::_cMsData
ContainerSingle2D _cMsData
Definition: TPTask.h:95
hcaldqm::constants::VERY_LOW
double const VERY_LOW
Definition: Constants.h:18
hcaldqm::constants::TPCHS_NUM
int const TPCHS_NUM[TPSUBDET_NUM]
Definition: Constants.h:137
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
TPTask::_tokData
edm::EDGetTokenT< HcalTrigPrimDigiCollection > _tokData
Definition: TPTask.h:52
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
TPTask::_cOccupancyEmulvsiphi_TPSubDetPM
Container1D _cOccupancyEmulvsiphi_TPSubDetPM
Definition: TPTask.h:92
TPTask::_fNames
std::vector< std::string > _fNames
Definition: TPTask.h:65
TPTask::_cNumEtMsmvsLS_TPSubDet
ContainerProf1D _cNumEtMsmvsLS_TPSubDet
Definition: TPTask.h:75
hcaldqm::UpdateFreq
UpdateFreq
Definition: DQTask.h:16
edm::SortedCollection< HcalTriggerPrimitiveDigi >::const_iterator
std::vector< HcalTriggerPrimitiveDigi >::const_iterator const_iterator
Definition: SortedCollection.h:82
TPTask::_cEtEmul_SubDet
Container1D _cEtEmul_SubDet
Definition: TPTask.h:69
HcalTrigTowerDetId::ieta
int ieta() const
get the tower ieta
Definition: HcalTrigTowerDetId.h:53
hcaldqm::ContainerSingle2D::setBinContent
virtual void setBinContent(int, int, double)
Definition: ContainerSingle2D.cc:104
hcaldqm::ContainerSingle2D::loadLabels
virtual void loadLabels(std::vector< std::string > const &)
Definition: ContainerSingle2D.cc:291
bookHistograms
void bookHistograms(fwlite::EventContainer &eventCont)
Definition: secvtxMassTemplates.cc:293
TPTask::endLuminosityBlock
virtual void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &)
Definition: TPTask.cc:351
TPTask::_cOccupancyEmul
ContainerSingle2D _cOccupancyEmul
Definition: TPTask.h:90
hcaldqm::ContainerSingle2D::book
virtual void book(DQMStore::IBooker &, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerSingle2D.cc:40
TPTask::_cFGMsm
ContainerSingle2D _cFGMsm
Definition: TPTask.h:86
hcaldqm::Container2D::fill
virtual void fill(HcalDetId const &)
Definition: Container2D.cc:39
TPTask::_cMsDatavsiphi_TPSubDetPM
Container1D _cMsDatavsiphi_TPSubDetPM
Definition: TPTask.h:96
TPTask::_cEtData_SubDet
Container1D _cEtData_SubDet
Definition: TPTask.h:68
hcaldqm::ContainerProf1D::book
virtual void book(DQMStore::IBooker &, std::string subsystem="Hcal", std::string aux="")
Definition: ContainerProf1D.cc:30
Logger::dqmthrow
void dqmthrow(std::string const &msg) const
Definition: Logger.h:15
TPTask::_cOccupancyDatavsLS_TPSubDet
ContainerProf1D _cOccupancyDatavsLS_TPSubDet
Definition: TPTask.h:93
TPTask::_cEtMsm
ContainerSingle2D _cEtMsm
Definition: TPTask.h:74
TPTask::_cOccupancyEmulvsLS_TPSubDet
ContainerProf1D _cOccupancyEmulvsLS_TPSubDet
Definition: TPTask.h:94
hcaldqm::Container1D::fill
virtual void fill(double)
Definition: Container1D.h:53
hcaldqm::Container2D::loadLabels
virtual void loadLabels(std::vector< std::string > const &)
Definition: Container2D.cc:294
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
j
int j
Definition: DBlmapReader.cc:9
TPTask::_cNumEtMsm_TPSubDet
Container1D _cNumEtMsm_TPSubDet
Definition: TPTask.h:76
hcaldqm::DQModule::_logger
Logger _logger
Definition: DQModule.h:62
TPTask::_cEtCorrRatio_TPSubDet
Container1D _cEtCorrRatio_TPSubDet
Definition: TPTask.h:72
TPTask::_cSummaryvsLS_TPSubDet
Container2D _cSummaryvsLS_TPSubDet
Definition: TPTask.h:105
hcaldqm::constants::GOOD
double const GOOD
Definition: Constants.h:15
min
T min(T a, T b)
Definition: MathUtil.h:58
hcaldqm::Container1D::initialize
virtual void initialize(std::string const &folder, std::string const &nametitle, mapper::MapperType mt, axis::Axis *xaxis, axis::Axis *yaxis=new ValueAxis(fYaxis, fEntries), int debug=0)
Definition: Container1D.cc:28
TPTask::_nMsmFG
int _nMsmFG[constants::TPSUBDET_NUM]
Definition: TPTask.h:57
hcaldqm::DQModule::_name
std::string _name
Definition: DQModule.h:50
TPTask::_cEtCorrRatiovsiphi_TPSubDetPM
ContainerProf1D _cEtCorrRatiovsiphi_TPSubDetPM
Definition: TPTask.h:73
TPTask::TPTask
TPTask(edm::ParameterSet const &ps)
Definition: TPTask.cc:5
hcaldqm::utilities::getTPSubDet
int getTPSubDet(HcalTrigTowerDetId const &)
Definition: Utilities.cc:58
TPTask::_cDigiSizeDatavsLS_TPSubDet
ContainerProf1D _cDigiSizeDatavsLS_TPSubDet
Definition: TPTask.h:100
TPTask::_cDigiSizeEmulvsLS_TPSubDet
ContainerProf1D _cDigiSizeEmulvsLS_TPSubDet
Definition: TPTask.h:101
hcaldqm::Container1D::at
virtual MonitorElement * at(unsigned int i)
Definition: Container1D.cc:236
HcalTrigTowerDetId::version
int version() const
get the version code for the trigger tower
Definition: HcalTrigTowerDetId.h:59
TPTask::_resetMonitors
virtual void _resetMonitors(UpdateFreq)
Definition: TPTask.cc:509
TPTask::_cOccupancyDatavsiphi_TPSubDetPM
Container1D _cOccupancyDatavsiphi_TPSubDetPM
Definition: TPTask.h:91
edm::InputTag::label
std::string const & label() const
Definition: InputTag.h:36
hcaldqm::Container2D::setBinContent
virtual void setBinContent(int, int, int, double)
Definition: Container2D.cc:263
TPTask::_cFGCorr_SubDet
Container2D _cFGCorr_SubDet
Definition: TPTask.h:85
HLT_FULL_cff.InputTag
tuple InputTag
Definition: HLT_FULL_cff.py:69759
hcaldqm::constants::NOT_APPLICABLE
double const NOT_APPLICABLE
Definition: Constants.h:19
TPTask::_nTPs_Emul
int _nTPs_Emul[constants::TPSUBDET_NUM]
Definition: TPTask.h:59
TPTask::_cSumdEtvsLS_TPSubDet
ContainerProf1D _cSumdEtvsLS_TPSubDet
Definition: TPTask.h:78
MonitorElement::getBinContent
double getBinContent(int binx) const
get content of bin (1-D)
Definition: MonitorElement.cc:743
TPTask::_cEtCorrRatiovsLS_TPSubDet
ContainerProf1D _cEtCorrRatiovsLS_TPSubDet
Definition: TPTask.h:71
TPTask::_tokEmul
edm::EDGetTokenT< HcalTrigPrimDigiCollection > _tokEmul
Definition: TPTask.h:53
hcaldqm::constants::TPSUBDET_NUM
int const TPSUBDET_NUM
Definition: Constants.h:88
TPTask::_cEtData_SubDet_ieta
Container1D _cEtData_SubDet_ieta
Definition: TPTask.h:82
TPTask
Definition: TPTask.h:18
TPTask::_cEtCorr_TPSubDet
Container2D _cEtCorr_TPSubDet
Definition: TPTask.h:70
TPTask::_nTPs_Data
int _nTPs_Data[constants::TPSUBDET_NUM]
Definition: TPTask.h:58
TPTask::_cOccupancyData
ContainerSingle2D _cOccupancyData
Definition: TPTask.h:89
hcaldqm::ContainerSingle2D::fill
virtual void fill(int, int)
Definition: ContainerSingle2D.cc:64
ntuplemaker.status
tuple status
Definition: ntuplemaker.py:245
TPTask::_tagEmul
edm::InputTag _tagEmul
Definition: TPTask.h:51
hcaldqm::Container1D::getBinContent
virtual double getBinContent(unsigned int, int)
Definition: Container1D.cc:225
hcaldqm::constants::PROBLEMATIC
double const PROBLEMATIC
Definition: Constants.h:16
edm::InputTag::instance
std::string const & instance() const
Definition: InputTag.h:37
TPTask::_cSummary
ContainerSingle2D _cSummary
Definition: TPTask.h:104
HcalTrigTowerDetId::iphi
int iphi() const
get the tower iphi
Definition: HcalTrigTowerDetId.h:55
TPTask::_process
virtual void _process(edm::Event const &, edm::EventSetup const &)
Definition: TPTask.cc:200
TPTask::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
Definition: TPTask.cc:163
edm::Run
Definition: Run.h:43
TPTask::_cMsEmul
ContainerSingle2D _cMsEmul
Definition: TPTask.h:97
Generated for CMSSW Reference Manual by doxygen 1.8.5