CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalZDCMonitor.cc
Go to the documentation of this file.
1 #include "DQM/HcalMonitorTasks/interface/HcalZDCMonitor.h"
2 
3 HcalZDCMonitor::HcalZDCMonitor() {
4 }
5 HcalZDCMonitor::~HcalZDCMonitor() {
6 }
7 void HcalZDCMonitor::reset() {
8 }
9 
10 void HcalZDCMonitor::setup(const edm::ParameterSet & ps, DQMStore * dbe) {
11  HcalBaseMonitor::setup(ps, dbe);
12 
13  baseFolder_ = rootFolder_ + "ZDCMonitor_Hcal";
14 
15  if (showTiming) {
16  cpu_timer.reset();
17  cpu_timer.start();
18  }
19 
20  if (fVerbosity > 0)
21  std::cout << "<HcalZDCMonitor::setup> Setting up histograms" << std::endl;
22 
23  if (fVerbosity > 1)
24  std::cout << "<HcalZDCMonitor::setup> Getting variable values from cfg files" << std::endl;
25 
26  // Set initial event # to 0
27  ievt_ = 0;
28 
29 //Histograms
30  if (m_dbe) {
31  if (fVerbosity > 1)
32  std::cout << "<HcalZDCMonitor::setup> Setting up Histograms" << std::endl;
33 
34  m_dbe->setCurrentFolder(baseFolder_);
35  meEVT_ = m_dbe->bookInt("ZDC Event Number");
36  meEVT_->Fill(ievt_);
37  char name[128];
38  char title[128];
39 
40  h_2D_charge = m_dbe->book2D("2D_DigiCharge", "Digi Charge (fC)", 2, 0, 2, 9, 0, 9);
41  h_2D_charge->setBinLabel(1,"ZDC+",1);
42  h_2D_charge->setBinLabel(2,"ZDC-",1);
43  h_2D_charge->setBinLabel(1,"EM1",2);
44  h_2D_charge->setBinLabel(2,"EM2",2);
45  h_2D_charge->setBinLabel(3,"EM3",2);
46  h_2D_charge->setBinLabel(4,"EM4",2);
47  h_2D_charge->setBinLabel(5,"EM5",2);
48  h_2D_charge->setBinLabel(6,"HAD1",2);
49  h_2D_charge->setBinLabel(7,"HAD2",2);
50  h_2D_charge->setBinLabel(8,"HAD3",2);
51  h_2D_charge->setBinLabel(9,"HAD4",2);
52 
53  h_2D_TSMean = m_dbe->book2D("2D_DigiTiming", "Digi Timing", 2, 0, 2, 9, 0, 9);
54  h_2D_TSMean->setBinLabel(1,"ZDC+",1);
55  h_2D_TSMean->setBinLabel(2,"ZDC-",1);
56  h_2D_TSMean->setBinLabel(1,"EM1",2);
57  h_2D_TSMean->setBinLabel(2,"EM2",2);
58  h_2D_TSMean->setBinLabel(3,"EM3",2);
59  h_2D_TSMean->setBinLabel(4,"EM4",2);
60  h_2D_TSMean->setBinLabel(5,"EM5",2);
61  h_2D_TSMean->setBinLabel(6,"HAD1",2);
62  h_2D_TSMean->setBinLabel(7,"HAD2",2);
63  h_2D_TSMean->setBinLabel(8,"HAD3",2);
64  h_2D_TSMean->setBinLabel(9,"HAD4",2);
65 
66  h_2D_RecHitEnergy = m_dbe->book2D("2D_RecHitEnergy", "Rechit Energy", 2, 0, 2, 9, 0, 9);
67  h_2D_RecHitEnergy->setBinLabel(1,"ZDC+",1);
68  h_2D_RecHitEnergy->setBinLabel(2,"ZDC-",1);
69  h_2D_RecHitEnergy->setBinLabel(1,"EM1",2);
70  h_2D_RecHitEnergy->setBinLabel(2,"EM2",2);
71  h_2D_RecHitEnergy->setBinLabel(3,"EM3",2);
72  h_2D_RecHitEnergy->setBinLabel(4,"EM4",2);
73  h_2D_RecHitEnergy->setBinLabel(5,"EM5",2);
74  h_2D_RecHitEnergy->setBinLabel(6,"HAD1",2);
75  h_2D_RecHitEnergy->setBinLabel(7,"HAD2",2);
76  h_2D_RecHitEnergy->setBinLabel(8,"HAD3",2);
77  h_2D_RecHitEnergy->setBinLabel(9,"HAD4",2);
78 
79  h_2D_RecHitTime = m_dbe->book2D("2D_RecHitTime", "Rechit Timing", 2, 0, 2, 9, 0, 9);
80  h_2D_RecHitTime->setBinLabel(1,"ZDC+",1);
81  h_2D_RecHitTime->setBinLabel(2,"ZDC-",1);
82  h_2D_RecHitTime->setBinLabel(1,"EM1",2);
83  h_2D_RecHitTime->setBinLabel(2,"EM2",2);
84  h_2D_RecHitTime->setBinLabel(3,"EM3",2);
85  h_2D_RecHitTime->setBinLabel(4,"EM4",2);
86  h_2D_RecHitTime->setBinLabel(5,"EM5",2);
87  h_2D_RecHitTime->setBinLabel(6,"HAD1",2);
88  h_2D_RecHitTime->setBinLabel(7,"HAD2",2);
89  h_2D_RecHitTime->setBinLabel(8,"HAD3",2);
90  h_2D_RecHitTime->setBinLabel(9,"HAD4",2);
91 
92  h_2D_saturation = m_dbe->book2D("h_2D_QIE", "Saturation Check", 2, 0, 2, 9, 0, 9);
93  h_2D_saturation->setBinLabel(1,"ZDC+",1);
94  h_2D_saturation->setBinLabel(2,"ZDC-",1);
95  h_2D_saturation->setBinLabel(1,"EM1",2);
96  h_2D_saturation->setBinLabel(2,"EM2",2);
97  h_2D_saturation->setBinLabel(3,"EM3",2);
98  h_2D_saturation->setBinLabel(4,"EM4",2);
99  h_2D_saturation->setBinLabel(5,"EM5",2);
100  h_2D_saturation->setBinLabel(6,"HAD1",2);
101  h_2D_saturation->setBinLabel(7,"HAD2",2);
102  h_2D_saturation->setBinLabel(8,"HAD3",2);
103  h_2D_saturation->setBinLabel(9,"HAD4",2);
104 
105  m_dbe->setCurrentFolder(baseFolder_ + "/Digis");
106 
107  for (int i = 0; i < 5; ++i) {
108  // pulse Plus Side
109  sprintf(title, "h_ZDCP_EMChan_%i_Pulse", i + 1);
110  sprintf(name, "ZDC Plus EM Section Pulse for channel %i", i + 1);
111  h_ZDCP_EM_Pulse[i] = m_dbe->book1D(title, name, 10, -0.5, 9.5);
112  h_ZDCP_EM_Pulse[i]->setAxisTitle("Time Slice id",1);
113  h_ZDCP_EM_Pulse[i]->setAxisTitle("Pulse Height",2);
114  // pulse Minus Side
115  sprintf(title, "h_ZDCM_EMChan_%i_Pulse", i + 1);
116  sprintf(name, "ZDC Minus EM Section Pulse for channel %i", i + 1);
117  h_ZDCM_EM_Pulse[i] = m_dbe->book1D(title, name, 10, -0.5, 9.5);
118  h_ZDCM_EM_Pulse[i]->setAxisTitle("Time Slice id",1);
119  h_ZDCM_EM_Pulse[i]->setAxisTitle("Pulse Height",2);
120  // integrated charge over 10 time samples
121  sprintf(title, "h_ZDCP_EMChan_%i_Charge", i + 1);
122  sprintf(name, "ZDC Plus EM Section Charge for channel %i", i + 1);
123  h_ZDCP_EM_Charge[i] = m_dbe->book1D(title, name, 1000, 0., 30000.);
124  h_ZDCP_EM_Charge[i]->setAxisTitle("Charge (fC)",1);
125  h_ZDCP_EM_Charge[i]->setAxisTitle("Events",2);
126  // integrated charge over 10 time samples
127  sprintf(title, "h_ZDCM_EMChan_%i_Charge", i + 1);
128  sprintf(name, "ZDC Minus EM Section Charge for channel %i", i + 1);
129  h_ZDCM_EM_Charge[i] = m_dbe->book1D(title, name, 1000, 0., 30000.);
130  h_ZDCM_EM_Charge[i]->setAxisTitle("Charge (fC)",1);
131  h_ZDCM_EM_Charge[i]->setAxisTitle("Events",2);
132  // charge weighted time slice
133  sprintf(title, "h_ZDCP_EMChan_%i_TSMean", i + 1);
134  sprintf(name, "ZDC Plus EM Section TSMean for channel %i", i + 1);
135  h_ZDCP_EM_TSMean[i] = m_dbe->book1D(title, name, 100, -0.5, 9.5);
136  h_ZDCP_EM_TSMean[i]->setAxisTitle("Timing",1);
137  h_ZDCP_EM_TSMean[i]->setAxisTitle("Events",2);
138  // charge weighted time slice
139  sprintf(title, "h_ZDCM_EMChan_%i_TSMean", i + 1);
140  sprintf(name, "ZDC Minus EM Section TSMean for channel %i", i + 1);
141  h_ZDCM_EM_TSMean[i] = m_dbe->book1D(title, name, 100, -0.5, 9.5);
142  h_ZDCM_EM_TSMean[i]->setAxisTitle("Timing",1);
143  h_ZDCM_EM_TSMean[i]->setAxisTitle("Events",2);
144  }
145 
146  for (int i = 0; i < 4; ++i) {
147  // pulse Plus Side
148  sprintf(title, "h_ZDCP_HADChan_%i_Pulse", i + 1);
149  sprintf(name, "ZDC Plus HAD Section Pulse for channel %i", i + 1);
150  h_ZDCP_HAD_Pulse[i] = m_dbe->book1D(title, name, 10, -0.5, 9.5);
151  h_ZDCP_HAD_Pulse[i]->setAxisTitle("Time Slice id",1);
152  h_ZDCP_HAD_Pulse[i]->setAxisTitle("Pulse Height",2);
153  // pulse Minus Side
154  sprintf(title, "h_ZDCM_HADChan_%i_Pulse", i + 1);
155  sprintf(name, "ZDC Minus HAD Section Pulse for channel %i", i + 1);
156  h_ZDCM_HAD_Pulse[i] = m_dbe->book1D(title, name, 10, -0.5, 9.5);
157  h_ZDCP_HAD_Pulse[i]->setAxisTitle("Time Slice id",1);
158  h_ZDCP_HAD_Pulse[i]->setAxisTitle("Pulse Height",2);
159  // integrated charge over 10 time samples
160  sprintf(title, "h_ZDCP_HADChan_%i_Charge", i + 1);
161  sprintf(name, "ZDC Plus HAD Section Charge for channel %i", i + 1);
162  h_ZDCP_HAD_Charge[i] = m_dbe->book1D(title, name, 1000, 0., 30000.);
163  h_ZDCP_HAD_Charge[i]->setAxisTitle("Charge (fC)",1);
164  h_ZDCP_HAD_Charge[i]->setAxisTitle("Events",2);
165  // integrated charge over 10 time samples
166  sprintf(title, "h_ZDCM_HADChan_%i_Charge", i + 1);
167  sprintf(name, "ZDC Minus HAD Section Charge for channel %i", i + 1);
168  h_ZDCM_HAD_Charge[i] = m_dbe->book1D(title, name, 1000, 0., 30000.);
169  h_ZDCM_HAD_Charge[i]->setAxisTitle("Charge (fC)",1);
170  h_ZDCM_HAD_Charge[i]->setAxisTitle("Events",2);
171  // charge weighted time slice
172  sprintf(title, "h_ZDCP_HADChan_%i_TSMean", i + 1);
173  sprintf(name, "ZDC Plus HAD Section TSMean for channel %i", i + 1);
174  h_ZDCP_HAD_TSMean[i] = m_dbe->book1D(title, name, 100, -0.5, 9.5);
175  h_ZDCP_HAD_TSMean[i]->setAxisTitle("Timing",1);
176  h_ZDCP_HAD_TSMean[i]->setAxisTitle("Events",2);
177  // charge weighted time slice
178  sprintf(title, "h_ZDCM_HADChan_%i_TSMean", i + 1);
179  sprintf(name, "ZDC Minus HAD Section TSMean for channel %i", i + 1);
180  h_ZDCM_HAD_TSMean[i] = m_dbe->book1D(title, name, 100, -0.5, 9.5);
181  h_ZDCM_HAD_TSMean[i]->setAxisTitle("Timing",1);
182  h_ZDCM_HAD_TSMean[i]->setAxisTitle("Events",2);
183  }
184 
185  m_dbe->setCurrentFolder(baseFolder_ + "/RecHits");
186 
187  for (int i = 0; i < 5; ++i) {
188  //RecHitEnergy Plus Side
189  sprintf(title,"h_ZDCP_EMChan_%i_RecHit_Energy",i+1);
190  sprintf(name,"ZDC EM Section Rechit Energy for channel %i",i+1);
191  h_ZDCP_EM_RecHitEnergy[i] = m_dbe->book1D(title, name, 1010, -100., 10000.);
192  h_ZDCP_EM_RecHitEnergy[i]->setAxisTitle("Energy (GeV)",1);
193  h_ZDCP_EM_RecHitEnergy[i]->setAxisTitle("Events",2);
194  //RecHitEnergy Minus Side
195  sprintf(title,"h_ZDCM_EMChan_%i_RecHit_Energy",i+1);
196  sprintf(name,"ZDC EM Section Rechit Energy for channel %i",i+1);
197  h_ZDCM_EM_RecHitEnergy[i] = m_dbe->book1D(title, name, 1010, -100., 10000.);
198  h_ZDCM_EM_RecHitEnergy[i]->setAxisTitle("Energy (GeV)",1);
199  h_ZDCM_EM_RecHitEnergy[i]->setAxisTitle("Events",2);
200  //RecHit Timing Plus Side
201  sprintf(title,"h_ZDCP_EMChan_%i_RecHit_Timing",i+1);
202  sprintf(name,"ZDC EM Section Rechit Timing for channel %i",i+1);
203  h_ZDCP_EM_RecHitTiming[i] = m_dbe->book1D(title, name, 100, -100., 100.);
204  h_ZDCP_EM_RecHitTiming[i]->setAxisTitle("RecHit Time",1);
205  h_ZDCP_EM_RecHitTiming[i]->setAxisTitle("Events",2);
206  //RecHit Timing Minus Side
207  sprintf(title,"h_ZDCM_EMChan_%i_RecHit_Timing",i+1);
208  sprintf(name,"ZDC EM Section Rechit Timing for channel %i",i+1);
209  h_ZDCM_EM_RecHitTiming[i] = m_dbe->book1D(title, name, 100, -100., 100.);
210  h_ZDCM_EM_RecHitTiming[i]->setAxisTitle("RecHit Time",1);
211  h_ZDCM_EM_RecHitTiming[i]->setAxisTitle("Events",2);
212  }
213 
214  for (int i = 0; i < 4; ++i) {
215  //RecHitEnergy Plus Side
216  sprintf(title,"h_ZDCP_HADChan_%i_RecHit_Energy",i+1);
217  sprintf(name,"ZDC HAD Section Rechit Energy for channel %i",i+1);
218  h_ZDCP_HAD_RecHitEnergy[i] = m_dbe->book1D(title, name, 1010, -100., 10000.);
219  h_ZDCP_HAD_RecHitEnergy[i]->setAxisTitle("Energy (GeV)",1);
220  h_ZDCP_HAD_RecHitEnergy[i]->setAxisTitle("Events",2);
221  //RecHitEnergy Minus Side
222  sprintf(title,"h_ZDCM_HADChan_%i_RecHit_Energy",i+1);
223  sprintf(name,"ZDC HAD Section Rechit Energy for channel %i",i+1);
224  h_ZDCM_HAD_RecHitEnergy[i] = m_dbe->book1D(title, name, 1010, -100., 10000.);
225  h_ZDCM_HAD_RecHitEnergy[i]->setAxisTitle("Energy (GeV)",1);
226  h_ZDCM_HAD_RecHitEnergy[i]->setAxisTitle("Events",2);
227  //RecHit Timing Plus Side
228  sprintf(title,"h_ZDCP_HADChan_%i_RecHit_Timing",i+1);
229  sprintf(name,"ZDC HAD Section Rechit Timing for channel %i",i+1);
230  h_ZDCP_HAD_RecHitTiming[i] = m_dbe->book1D(title, name, 100, -100., 100.);
231  h_ZDCP_HAD_RecHitTiming[i]->setAxisTitle("RecHit Time",1);
232  h_ZDCP_HAD_RecHitTiming[i]->setAxisTitle("Events",2);
233  //RecHit Timing Minus Side
234  sprintf(title,"h_ZDCM_HADChan_%i_RecHit_Timing",i+1);
235  sprintf(name,"ZDC HAD Section Rechit Timing for channel %i",i+1);
236  h_ZDCM_HAD_RecHitTiming[i] = m_dbe->book1D(title, name, 100, -100., 100.);
237  h_ZDCM_HAD_RecHitTiming[i]->setAxisTitle("RecHit Time",1);
238  h_ZDCM_HAD_RecHitTiming[i]->setAxisTitle("Events",2);
239  }
240 
241  }
242  return;
243 }
244 
245 void HcalZDCMonitor::processEvent(const ZDCDigiCollection& digi, const ZDCRecHitCollection& rechit) {
246  if (fVerbosity > 0)
247  std::cout << "<HcalZDCMonitor::processEvent> Processing Event..." << std::endl;
248  if (showTiming)
249  {
250  cpu_timer.reset();
251  cpu_timer.start();
252  }
253  ++ievt_;
254  meEVT_->Fill(ievt_);
255 
256  //--------------------------------------
257  // ZDC Digi part
258  //--------------------------------------
259  double fSum = 0.;
260  std::vector<double> fData;
261  double digiThresh = 99.5; //corresponds to 40 ADC counts
262  //int digiThreshADC = 40;
263  int digiSaturation = 127;
264  //double ZDCQIEConst = 2.6;
265 
266  for (ZDCDigiCollection::const_iterator digi_iter = digi.begin();
267  digi_iter != digi.end(); ++digi_iter)
268  {
269  const ZDCDataFrame digi = (const ZDCDataFrame) (*digi_iter);
270  //HcalZDCDetId id(digi_iter->id());
271  int iSide = digi_iter->id().zside();
272  int iSection = digi_iter->id().section();
273  int iChannel = digi_iter->id().channel();
274 
275  unsigned int fTS = digi_iter->size();
276  while (fData.size()<fTS)
277  fData.push_back(-999);
278  while (fData.size()>fTS)
279  fData.pop_back(); // delete last elements
280 
281  fSum = 0.;
282  bool saturated = false;
283  for (unsigned int i = 0; i < fTS; ++i)
284  {
285  //fData[i]=digi[i].nominal_fC() * ZDCQIEConst;
286  fData[i]=digi[i].nominal_fC();
287  if (digi[i].adc()==digiSaturation){
288  saturated=true;
289  }
290  }
291 
292  double fTSMean = 0;
293  if (fData.size()>6)
294  fTSMean = getTime(fData, 4, 6, fSum); // tsmin = 4, tsmax = 6.
295  //std::cout << "Side= " << iSide << " Section= " << iSection << " Channel= " << iChannel << "\tCharge\t" << fSum <<std::endl;
296  if (saturated==true){
297  h_2D_saturation->Fill(iSide==1?0:1,iSection==1?iChannel-1:iChannel+4,1);
298  }
299 
300  if (iSection == 1)
301  { // EM
302  if (iSide == 1) { // Plus
303  for (unsigned int i = 0; i < fTS; ++i) {
304  if (fData[i] > digiThresh) h_ZDCP_EM_Pulse[iChannel - 1]->Fill(i, fData[i]);
305  }
306  if (fSum > digiThresh) {
307  h_ZDCP_EM_Charge[iChannel - 1]->Fill(fSum);
308  h_ZDCP_EM_TSMean[iChannel - 1]->Fill(fTSMean);
309  //std::cout<< "fSum " << fSum << " fTSMean " << fTSMean <<std::endl;
310  }
311  } // Plus
312  if (iSide == -1) { // Minus
313  for (unsigned int i = 0; i < fTS; ++i) {
314  if (fData[i] > digiThresh) h_ZDCM_EM_Pulse[iChannel - 1]->Fill(i, fData[i]);
315  }
316  if (fSum > digiThresh) {
317  h_ZDCM_EM_Charge[iChannel - 1]->Fill(fSum);
318  h_ZDCM_EM_TSMean[iChannel - 1]->Fill(fTSMean);
319  }
320  } // Minus
321  }// EM
322 
323  else if (iSection == 2)
324  { // HAD
325  if (iSide == 1) { // Plus
326  for (unsigned int i = 0; i < fTS; ++i) {
327  if (fData[i] > digiThresh) h_ZDCP_HAD_Pulse[iChannel - 1]->Fill(i, fData[i]);
328  }
329  if (fSum > digiThresh) {
330  h_ZDCP_HAD_Charge[iChannel - 1]->Fill(fSum);
331  h_ZDCP_HAD_TSMean[iChannel - 1]->Fill(fTSMean);
332  }
333  } // Plus
334  if (iSide == -1) { // Minus
335  for (unsigned int i = 0; i < fTS; ++i) {
336  if (fData[i] > digiThresh) h_ZDCM_HAD_Pulse[iChannel - 1]->Fill(i, fData[i]);
337  }
338  if (fSum > digiThresh) {
339  h_ZDCM_HAD_Charge[iChannel - 1]->Fill(fSum);
340  h_ZDCM_HAD_TSMean[iChannel - 1]->Fill(fTSMean);
341  }
342  }// Minus
343  } // HAD
344  } // loop on zdc digi collection
345 
346 
347  //--------------------------------------
348  // ZDC RecHit part
349  //--------------------------------------
350  for (ZDCRecHitCollection::const_iterator rechit_iter = rechit.begin();
351  rechit_iter != rechit.end(); ++rechit_iter)
352  {
353  HcalZDCDetId id(rechit_iter->id());
354  int Side = (rechit_iter->id()).zside();
355  int Section = (rechit_iter->id()).section();
356  int Channel = (rechit_iter->id()).channel();
357  //std::cout << "RecHitEnergy " << zhit->energy() << " RecHitTime " << zhit->time() << std::endl;
358 
359  if(Section==1)
360  { //EM
361  if (Side ==1 ){ // Plus
362  h_ZDCP_EM_RecHitEnergy[Channel-1]->Fill(rechit_iter->energy());
363  h_ZDCP_EM_RecHitTiming[Channel-1]->Fill(rechit_iter->time());
364  }
365  if (Side == -1 ){ //Minus
366  h_ZDCM_EM_RecHitEnergy[Channel-1]->Fill(rechit_iter->energy());
367  h_ZDCM_EM_RecHitTiming[Channel-1]->Fill(rechit_iter->time());
368  }
369  } //EM
370  else if(Section==2)
371  { //HAD
372  if (Side ==1 ){ //Plus
373  h_ZDCP_HAD_RecHitEnergy[Channel-1]->Fill(rechit_iter->energy());
374  h_ZDCP_HAD_RecHitTiming[Channel-1]->Fill(rechit_iter->time());
375  }
376  if (Side == -1 ){ //Minus
377  h_ZDCM_HAD_RecHitEnergy[Channel-1]->Fill(rechit_iter->energy());
378  h_ZDCM_HAD_RecHitTiming[Channel-1]->Fill(rechit_iter->time());
379  }
380  } // HAD
381  } // loop on rechits
382 
383 } // end of event processing
384 /*
385 ------------------------------------------------------------------------------------
386 // This is what we did to find the good signal. After we've started to use only time slice 4,5,6.
387 bool HcalZDCMonitor::isGood(std::vector<double>fData, double fCut, double fPercentage) {
388  bool dec = false;
389  int ts_max = -1;
390 
391  ts_max = getTSMax(fData);
392  if (ts_max == 0 || ts_max == (int)(fData.size() - 1))
393  return false;
394  float sum = fData[ts_max - 1] + fData[ts_max + 1];
395 
396  // cout << "tsMax " << ts_max << " data[tsmax] " << mData[ts_max] << " sum " << sum << endl;
397  if (fData[ts_max] > fCut && sum > (fData[ts_max] * fPercentage))
398  dec = true;
399  return dec;
400 } // bool HcalZDCMonitor::isGood
401 
402 int HcalZDCMonitor::getTSMax(std::vector<double>fData)
403 {
404  int ts_max = -100;
405  double max = -999.;
406 
407  for (unsigned int j = 0; j < fData.size(); ++j) {
408  if (max < fData[j]) {
409  max = fData[j];
410  ts_max = j;
411  }
412  }
413  return ts_max;
414 } // int HcalZDCMonitor::getTSMax()
415 ------------------------------------------------------------------------------------
416 */
417 double HcalZDCMonitor::getTime(const std::vector<double>& fData, unsigned int ts_min, unsigned int ts_max, double &fSum) {
418  double weightedTime = 0.;
419  double SumT = 0.;
420  double Time = -999.;
421  double digiThreshf = 99.5;
422 
423  for (unsigned int ts=ts_min; ts<=ts_max; ++ts) {
424  if (fData[ts] > digiThreshf){
425  weightedTime += ts * fData[ts];
426  SumT += fData[ts];
427  }
428  }
429 
430  if (SumT > 0.) {
431  Time = weightedTime / SumT;
432  }
433 
434  fSum = SumT;
435 
436  return Time;
437 
438 } //double HcalZDCMonitor::getTime()
439 
440 
441 void HcalZDCMonitor::endLuminosityBlock()
442 {
443 
444  for (int i = 0; i < 5; ++i) { // EM Channels
445  // ZDC Plus
446  h_2D_charge->setBinContent(1, i + 1, h_ZDCP_EM_Charge[i]->getMean());
447  h_2D_TSMean->setBinContent(1, i + 1, h_ZDCP_EM_TSMean[i]->getMean());
448  h_2D_RecHitEnergy->setBinContent(1, i + 1, h_ZDCP_EM_RecHitEnergy[i]->getMean());
449  h_2D_RecHitTime->setBinContent(1, i + 1, h_ZDCP_EM_RecHitTiming[i]->getMean());
450  // ZDC Minus
451  h_2D_charge->setBinContent(2, i + 1, h_ZDCM_EM_Charge[i]->getMean());
452  h_2D_TSMean->setBinContent(2, i + 1, h_ZDCM_EM_TSMean[i]->getMean());
453  h_2D_RecHitEnergy->setBinContent(2, i + 1, h_ZDCM_EM_RecHitEnergy[i]->getMean());
454  h_2D_RecHitTime->setBinContent(2, i + 1, h_ZDCM_EM_RecHitTiming[i]->getMean());
455  }
456 
457  for (int i = 0; i < 4; ++i) { // HAD channels
458  // ZDC Plus
459  h_2D_charge->setBinContent(1, i + 6, h_ZDCP_HAD_Charge[i]->getMean());
460  h_2D_TSMean->setBinContent(1, i + 6, h_ZDCP_HAD_TSMean[i]->getMean());
461  h_2D_RecHitEnergy->setBinContent(1, i + 6, h_ZDCP_HAD_RecHitEnergy[i]->getMean());
462  h_2D_RecHitTime->setBinContent(1, i + 6, h_ZDCP_HAD_RecHitTiming[i]->getMean());
463  // ZDC Minus
464  //h_ZDCM_HAD_Pulse[i]->Scale(10. / h_ZDCM_HAD_Pulse[i]->getEntries());
465  h_2D_charge->setBinContent(2, i + 6, h_ZDCM_HAD_Charge[i]->getMean());
466  h_2D_TSMean->setBinContent(2, i + 6, h_ZDCM_HAD_TSMean[i]->getMean());
467  h_2D_RecHitEnergy->setBinContent(2, i + 6, h_ZDCM_HAD_RecHitEnergy[i]->getMean());
468  h_2D_RecHitTime->setBinContent(2, i + 6, h_ZDCM_HAD_RecHitTiming[i]->getMean());
469  }
470 } // void HcalZDCMonitor::endLuminosityBlock()
471 
ecalMGPA::adc
int adc(sample_type sample)
get the ADC sample (12 bits)
Definition: EcalMGPASample.h:11
HcalZDCMonitor::h_ZDCP_HAD_TSMean
MonitorElement * h_ZDCP_HAD_TSMean[4]
Definition: HcalZDCMonitor.h:63
i
int i
Definition: DBlmapReader.cc:9
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:825
HcalZDCMonitor::h_ZDCP_HAD_RecHitTiming
MonitorElement * h_ZDCP_HAD_RecHitTiming[4]
Definition: HcalZDCMonitor.h:71
HcalZDCMonitor::h_ZDCM_EM_Pulse
MonitorElement * h_ZDCM_EM_Pulse[5]
Definition: HcalZDCMonitor.h:54
HcalZDCMonitor::h_ZDCM_HAD_TSMean
MonitorElement * h_ZDCM_HAD_TSMean[4]
Definition: HcalZDCMonitor.h:64
edm::CPUTimer::start
void start()
Definition: CPUTimer.cc:74
HcalZDCMonitor::h_ZDCM_EM_TSMean
MonitorElement * h_ZDCM_EM_TSMean[5]
Definition: HcalZDCMonitor.h:58
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:872
HcalZDCMonitor::h_ZDCM_EM_RecHitTiming
MonitorElement * h_ZDCM_EM_RecHitTiming[5]
Definition: HcalZDCMonitor.h:68
HcalZDCMonitor::h_2D_TSMean
MonitorElement * h_2D_TSMean
Definition: HcalZDCMonitor.h:50
HcalZDCMonitor::h_ZDCM_HAD_RecHitTiming
MonitorElement * h_ZDCM_HAD_RecHitTiming[4]
Definition: HcalZDCMonitor.h:72
HcalZDCMonitor::h_ZDCP_EM_RecHitEnergy
MonitorElement * h_ZDCP_EM_RecHitEnergy[5]
Definition: HcalZDCMonitor.h:65
edm::SortedCollection< ZDCDataFrame >::const_iterator
std::vector< ZDCDataFrame >::const_iterator const_iterator
Definition: SortedCollection.h:80
HcalZDCMonitor::h_ZDCP_EM_Pulse
MonitorElement * h_ZDCP_EM_Pulse[5]
Definition: HcalZDCMonitor.h:53
MonitorElement::setBinLabel
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:901
edm::CPUTimer::reset
void reset()
Definition: CPUTimer.cc:107
HcalZDCMonitor::h_ZDCP_EM_TSMean
MonitorElement * h_ZDCP_EM_TSMean[5]
Definition: HcalZDCMonitor.h:57
HcalZDCMonitor::h_ZDCM_EM_Charge
MonitorElement * h_ZDCM_EM_Charge[5]
Definition: HcalZDCMonitor.h:56
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:161
HcalZDCMonitor::h_ZDCP_HAD_Charge
MonitorElement * h_ZDCP_HAD_Charge[4]
Definition: HcalZDCMonitor.h:61
HcalZDCMonitor::setup
void setup(const edm::ParameterSet &ps, DQMStore *dbe)
Definition: HcalZDCMonitor.cc:10
HcalBaseMonitor::rootFolder_
std::string rootFolder_
Definition: HcalBaseMonitor.h:115
HcalZDCMonitor::h_2D_saturation
MonitorElement * h_2D_saturation
Definition: HcalZDCMonitor.h:48
HcalZDCMonitor::h_2D_charge
MonitorElement * h_2D_charge
Definition: HcalZDCMonitor.h:49
HcalBaseMonitor::cpu_timer
edm::CPUTimer cpu_timer
Definition: HcalBaseMonitor.h:108
HcalZDCMonitor::h_2D_RecHitTime
MonitorElement * h_2D_RecHitTime
Definition: HcalZDCMonitor.h:52
HcalZDCMonitor::h_ZDCM_HAD_Charge
MonitorElement * h_ZDCM_HAD_Charge[4]
Definition: HcalZDCMonitor.h:62
HcalBaseMonitor::baseFolder_
std::string baseFolder_
Definition: HcalBaseMonitor.h:116
HcalZDCMonitor::h_ZDCM_HAD_RecHitEnergy
MonitorElement * h_ZDCM_HAD_RecHitEnergy[4]
Definition: HcalZDCMonitor.h:70
edm::SortedCollection::end
const_iterator end() const
Definition: SortedCollection.h:301
HcalZDCMonitor::endLuminosityBlock
void endLuminosityBlock(void)
Definition: HcalZDCMonitor.cc:441
HcalZDCMonitor::getTime
double getTime(const std::vector< double > &fData, unsigned int ts_min, unsigned int ts_max, double &fSum)
Definition: HcalZDCMonitor.cc:417
HcalZDCMonitor::h_2D_RecHitEnergy
MonitorElement * h_2D_RecHitEnergy
Definition: HcalZDCMonitor.h:51
HcalZDCMonitor::h_ZDCP_HAD_Pulse
MonitorElement * h_ZDCP_HAD_Pulse[4]
Definition: HcalZDCMonitor.h:59
HcalZDCMonitor::meEVT_
MonitorElement * meEVT_
Definition: HcalZDCMonitor.h:47
HcalZDCMonitor::h_ZDCM_EM_RecHitEnergy
MonitorElement * h_ZDCM_EM_RecHitEnergy[5]
Definition: HcalZDCMonitor.h:66
edm::SortedCollection::size
size_type size() const
Definition: SortedCollection.h:227
HcalZDCMonitor::h_ZDCP_HAD_RecHitEnergy
MonitorElement * h_ZDCP_HAD_RecHitEnergy[4]
Definition: HcalZDCMonitor.h:69
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:121
dqm::me_util::Channel
DQMChannel Channel
Definition: DQMDefinitions.h:24
DQMStore::bookInt
MonitorElement * bookInt(const char *name)
Book int.
Definition: DQMStore.cc:779
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:1000
HcalZDCMonitor::h_ZDCP_EM_Charge
MonitorElement * h_ZDCP_EM_Charge[5]
Definition: HcalZDCMonitor.h:55
MonitorElement::setAxisTitle
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:928
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:584
edm::SortedCollection::begin
const_iterator begin() const
Definition: SortedCollection.h:294
HcalBaseMonitor::setup
virtual void setup(const edm::ParameterSet &ps, DQMStore *dbe)
Definition: HcalBaseMonitor.cc:21
HcalZDCMonitor::h_ZDCP_EM_RecHitTiming
MonitorElement * h_ZDCP_EM_RecHitTiming[5]
Definition: HcalZDCMonitor.h:67
HcalZDCMonitor::h_ZDCM_HAD_Pulse
MonitorElement * h_ZDCM_HAD_Pulse[4]
Definition: HcalZDCMonitor.h:60
Generated for CMSSW Reference Manual by doxygen 1.8.5