CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
DQMHcalPhiSymAlCaReco.cc
Go to the documentation of this file.
1 /*
2  * \file DQMHcalPhiSymAlCaReco.cc
3  *
4  * \author Olga Kodolova
5  *
6  *
7  *
8  * Description: Monitoring of Phi Symmetry Calibration Stream
9 */
10 
15 
17 
18 // DQM include files
19 
21 
22 // work on collections
23 
28 
34 
37 
38 using namespace std;
39 using namespace edm;
40 
41 // ******************************************
42 // constructors
43 // *****************************************
44 
46 eventCounter_(0)
47 {
48  //
49  // Input from configurator file
50  //
51  folderName_ = ps.getUntrackedParameter<string>("FolderName","ALCAStreamHcalPhiSym");
52 
53  hbherecoMB = consumes<HBHERecHitCollection>(ps.getParameter<edm::InputTag>("hbheInputMB"));
54  horecoMB = ps.getParameter<edm::InputTag>("hoInputMB");
55  hfrecoMB = consumes<HFRecHitCollection>(ps.getParameter<edm::InputTag>("hfInputMB"));
56 
57  hbherecoNoise = consumes<HBHERecHitCollection>(ps.getParameter<edm::InputTag>("hbheInputNoise"));
58  horecoNoise = ps.getParameter<edm::InputTag>("hoInputNoise");
59  hfrecoNoise = consumes<HFRecHitCollection>(ps.getParameter<edm::InputTag>("hfInputNoise"));
60 
61  rawInLabel_ = consumes<FEDRawDataCollection>(ps.getParameter<edm::InputTag>("rawInputLabel"));
62 
63  period_ = ps.getParameter<unsigned int>("period") ;
64 
65  saveToFile_ = ps.getUntrackedParameter<bool>("SaveToFile",false);
66  fileName_ = ps.getUntrackedParameter<string>("FileName","MonitorAlCaHcalPhiSym.root");
67 
68  // histogram parameters
69 
70  // Distribution of rechits in iPhi, iEta
71  hiDistr_y_nbin_ = ps.getUntrackedParameter<int>("hiDistr_y_nbin",72);
72  hiDistr_y_min_ = ps.getUntrackedParameter<double>("hiDistr_y_min",0.5);
73  hiDistr_y_max_ = ps.getUntrackedParameter<double>("hiDistr_y_max",72.5);
74  hiDistr_x_nbin_ = ps.getUntrackedParameter<int>("hiDistr_x_nbin",41);
75  hiDistr_x_min_ = ps.getUntrackedParameter<double>("hiDistr_x_min",0.5);
76  hiDistr_x_max_ = ps.getUntrackedParameter<double>("hiDistr_x_max",41.5);
77  // Check for NZS
78  hiDistr_r_nbin_ = ps.getUntrackedParameter<int>("hiDistr_r_nbin",100);
79  ihbhe_size_ = ps.getUntrackedParameter<double>("ihbhe_size_",5184.);
80  ihf_size_ = ps.getUntrackedParameter<double>("ihf_size_",1728.);
81 
82 }
83 
85 {}
86 
87 //--------------------------------------------------------
89 
90  // create and cd into new folder
92 
93  eventCounter_ = 0;
94 
95  hFEDsize = ibooker.book1D("hFEDsize","HCAL FED size (kB)",200,-0.5,20.5);
96  hFEDsize->setAxisTitle("kB",1);
97 
98  hHcalIsZS = ibooker.book1D("hHcalIsZS", "Hcal Is ZS", 4, -1.5, 2.5);
99  hHcalIsZS->setBinLabel(2, "NZS");
100  hHcalIsZS->setBinLabel(3, "ZS");
101 
102  char hname[50];
103  sprintf(hname, "L1 Event Number %% %i", period_);
104  hL1Id = ibooker.book1D("hL1Id", hname,4200,-99.5,4099.5);
105  hL1Id->setAxisTitle(hname);
106 
107 
108  // book some histograms 1D
109  double xmin = 0.1;
110  double xmax = 1.1;
112  ibooker.book1D("DistrHBHEsize","Size of HBHE Collection",
114  xmin,
115  xmax
116  );
118  ibooker.book1D("DistrHFsize","Size of HF Collection",
120  xmin,
121  xmax
122  );
123 
124 
125  // First moment
126  hiDistrMBPl2D_ =
127  ibooker.book2D("MBdepthPl1", "iphi- +ieta signal distribution at depth1",
134  );
135 
136  hiDistrMBPl2D_->setAxisTitle("i#phi ", 2);
137  hiDistrMBPl2D_->setAxisTitle("i#eta ", 1);
138 
139 
141  ibooker.book2D("NoisedepthPl1", "iphi-ieta noise distribution at depth1",
142  hiDistr_x_nbin_+1,
143  hiDistr_x_min_-1.,
145  hiDistr_y_nbin_+1,
146  hiDistr_y_min_-1.,
148  );
149 
150  hiDistrNoisePl2D_->setAxisTitle("i#phi ", 2);
151  hiDistrNoisePl2D_->setAxisTitle("i#eta ", 1);
152 // Second moment
154  ibooker.book2D("MB2depthPl1", "iphi- +ieta signal distribution at depth1",
161  );
162 
163  hiDistrMB2Pl2D_->setAxisTitle("i#phi ", 2);
164  hiDistrMB2Pl2D_->setAxisTitle("i#eta ", 1);
165 
166 
168  ibooker.book2D("Noise2depthPl1", "iphi-ieta noise distribution at depth1",
175  );
176 
177  hiDistrNoise2Pl2D_->setAxisTitle("i#phi ", 2);
178  hiDistrNoise2Pl2D_->setAxisTitle("i#eta ", 1);
179 
180 // Variance
182  ibooker.book2D("VarMBdepthPl1", "iphi- +ieta signal distribution at depth1",
189  );
190 
191  hiDistrVarMBPl2D_->setAxisTitle("i#phi ", 2);
192  hiDistrVarMBPl2D_->setAxisTitle("i#eta ", 1);
193 
194 
196  ibooker.book2D("VarNoisedepthPl1", "iphi-ieta noise distribution at depth1",
203  );
204 
205  hiDistrVarNoisePl2D_->setAxisTitle("i#phi ", 2);
206  hiDistrVarNoisePl2D_->setAxisTitle("i#eta ", 1);
207 
208 //==================================================================================
209 // First moment
210  hiDistrMBMin2D_ =
211  ibooker.book2D("MBdepthMin1", "iphi- +ieta signal distribution at depth1",
218  );
219 
220  hiDistrMBMin2D_->setAxisTitle("i#phi ", 2);
221  hiDistrMBMin2D_->setAxisTitle("i#eta ", 1);
222 
223 
225  ibooker.book2D("NoisedepthMin1", "iphi-ieta noise distribution at depth1",
232  );
233 
234  hiDistrNoiseMin2D_->setAxisTitle("i#phi ", 2);
235  hiDistrNoiseMin2D_->setAxisTitle("i#eta ", 1);
236 // Second moment
238  ibooker.book2D("MB2depthMin1", "iphi- +ieta signal distribution at depth1",
245  );
246 
247  hiDistrMB2Min2D_->setAxisTitle("i#phi ", 2);
248  hiDistrMB2Min2D_->setAxisTitle("i#eta ", 1);
249 
250 
252  ibooker.book2D("Noise2depthMin1", "iphi-ieta noise distribution at depth1",
259  );
260 
261  hiDistrNoise2Min2D_->setAxisTitle("i#phi ", 2);
262  hiDistrNoise2Min2D_->setAxisTitle("i#eta ", 1);
263 
264 // Variance
266  ibooker.book2D("VarMBdepthMin1", "iphi- +ieta signal distribution at depth1",
273  );
274 
275  hiDistrVarMBMin2D_->setAxisTitle("i#phi ", 2);
276  hiDistrVarMBMin2D_->setAxisTitle("i#eta ", 1);
277 
278 
280  ibooker.book2D("VarNoisedepthMin1", "iphi-ieta noise distribution at depth1",
287  );
288 
289  hiDistrVarNoiseMin2D_->setAxisTitle("i#phi ", 2);
290  hiDistrVarNoiseMin2D_->setAxisTitle("i#eta ", 1);
291 
292 
293 }
294 
295 //--------------------------------------------------------
296 //void DQMHcalPhiSymAlCaReco::beginRun(const edm::Run& r, const EventSetup& context) {
298 //}
299 
300 //--------------------------------------------------------
302  const EventSetup& context) {
303 
304 }
305 
306 //-------------------------------------------------------------
307 
309  const EventSetup& iSetup ){
310 
311 
312  eventCounter_++;
313 
315  iEvent.getByToken(rawInLabel_,rawIn);
316 
317  if(!rawIn.isValid()){
318  LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
319  return ;
320  }
321 
322  //get HCAL FEDs:
323  std::vector<int> selFEDs;
325  {
326  selFEDs.push_back(i);
327  }
328 
329  // std::cout<<" Size of FED "<<selFEDs.size()<<std::endl;
330 
331  const FEDRawDataCollection *rdc=rawIn.product();
332 
333  bool hcalIsZS = false;
334  int lvl1ID = 0;
335  bool lvl1IDFound = false;
336  for (unsigned int k=0; k<selFEDs.size(); k++) {
337  const FEDRawData & fedData = rdc->FEDData(selFEDs[k]);
338  // std::cout<<fedData.size()*std::pow(1024.,-1)<<std::endl;
339  hFEDsize->Fill(fedData.size()*std::pow(1024.,-1),1);
340 
341  // get HCAL DCC Header for each FEDRawData
342  const HcalDCCHeader* dccHeader=(const HcalDCCHeader*)(fedData.data());
343  if (dccHeader) {
344  // walk through the HTR data...
345  HcalHTRData htr;
346 
347  int nspigot =0;
348  for (int spigot=0; spigot<HcalDCCHeader::SPIGOT_COUNT; spigot++) {
349  nspigot++;
350 
351  if (!dccHeader->getSpigotPresent(spigot)) continue;
352 
353  // Load the given decoder with the pointer and length from this spigot.
354  dccHeader->getSpigotData(spigot,htr, fedData.size());
355 
356  if(k != 20 && nspigot !=14 ) {
357  if ( !htr.isUnsuppressed() ) { hcalIsZS = true; }
358  }
359  }
360  }
361 
362  // try to get the lvl1ID from the HCAL fed
363  if (!lvl1IDFound && (fedData.size() > 0)) {
364  // get FED Header for FEDRawData
365  FEDHeader fedHeader(fedData.data());
366  lvl1ID = fedHeader.lvl1ID();
367  lvl1IDFound = true;
368  }
369  } // loop over HcalFEDs
370 
371  hHcalIsZS->Fill( hcalIsZS );
372  hL1Id->Fill( lvl1ID%period_ );
373 
375  iEvent.getByToken(hbherecoNoise, hbheNS);
376 
377  if(!hbheNS.isValid()){
378  LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
379  return ;
380  }
381 
383  iEvent.getByToken(hbherecoMB, hbheMB);
384 
385  if(!hbheMB.isValid()){
386  LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
387  return ;
388  }
389 
391  iEvent.getByToken(hfrecoNoise, hfNS);
392 
393  if(!hfNS.isValid()){
394  LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
395  return ;
396  }
397 
399  iEvent.getByToken(hfrecoMB, hfMB);
400 
401  if(!hfMB.isValid()){
402  LogDebug("") << "HcalCalibAlgos: Error! can't get hbhe product!" << std::endl;
403  return ;
404  }
405 
406  const HBHERecHitCollection HithbheNS = *(hbheNS.product());
407 
408  hiDistrHBHEsize1D_->Fill(HithbheNS.size()/ihbhe_size_);
409 
410 
411  for(HBHERecHitCollection::const_iterator hbheItr=HithbheNS.begin(); hbheItr!=HithbheNS.end(); hbheItr++)
412  {
413  DetId id = (*hbheItr).detid();
414  HcalDetId hid=HcalDetId(id);
415 
416  if(hid.depth() == 1) {
417  if( hid.ieta() > 0 ) {
418  hiDistrNoisePl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy());
419  hiDistrNoise2Pl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy()*hbheItr->energy());
420  } else {
421  hiDistrNoiseMin2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy());
422  hiDistrNoise2Min2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy()*hbheItr->energy());
423  }
424  }
425  }
426 
427  const HBHERecHitCollection HithbheMB = *(hbheMB.product());
428 
429  for(HBHERecHitCollection::const_iterator hbheItr=HithbheMB.begin(); hbheItr!=HithbheMB.end(); hbheItr++)
430  {
431  DetId id = (*hbheItr).detid();
432  HcalDetId hid=HcalDetId(id);
433 
434  if(hid.depth() == 1) {
435  if( hid.ieta() > 0 ) {
436  hiDistrMBPl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy());
437  hiDistrMB2Pl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy()*hbheItr->energy());
438  } else {
439  hiDistrMBMin2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy());
440  hiDistrMB2Min2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy()*hbheItr->energy());
441  }
442  }
443 
444  }
445 
446  const HFRecHitCollection HithfNS = *(hfNS.product());
447 
448  hiDistrHFsize1D_->Fill(HithfNS.size()/ihf_size_);
449 
450  for(HFRecHitCollection::const_iterator hbheItr=HithfNS.begin(); hbheItr!=HithfNS.end(); hbheItr++)
451  {
452 
453  DetId id = (*hbheItr).detid();
454  HcalDetId hid=HcalDetId(id);
455 
456  if(hid.depth() == 1) {
457  if( hid.ieta() > 0 ) {
458  hiDistrNoisePl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy());
459  hiDistrNoise2Pl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy()*hbheItr->energy());
460  } else {
461  hiDistrNoiseMin2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy());
462  hiDistrNoise2Min2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy()*hbheItr->energy());
463  }
464  }
465 
466  }
467 
468  const HFRecHitCollection HithfMB = *(hfMB.product());
469 
470  for(HFRecHitCollection::const_iterator hbheItr=HithfMB.begin(); hbheItr!=HithfMB.end(); hbheItr++)
471  {
472  DetId id = (*hbheItr).detid();
473  HcalDetId hid=HcalDetId(id);
474 
475  if(hid.depth() == 1) {
476  if( hid.ieta() > 0 ) {
477  hiDistrMBPl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy());
478  hiDistrMB2Pl2D_->Fill(hid.ieta(),hid.iphi(),hbheItr->energy()*hbheItr->energy());
479  } else {
480  hiDistrMBMin2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy());
481  hiDistrMB2Min2D_->Fill(fabs(hid.ieta()),hid.iphi(),hbheItr->energy()*hbheItr->energy());
482  }
483  }
484  }
485 
486 
487 } //analyze
488 
489 
490 
491 
492 //--------------------------------------------------------
494  const EventSetup& context) {
495 }
496 //--------------------------------------------------------
497 void DQMHcalPhiSymAlCaReco::endRun(const Run& r, const EventSetup& context){
498 // Keep Variances
499  if(eventCounter_ > 0) {
500  for(int k=0; k<=hiDistr_x_nbin_;k++)
501  {
502  for(int j=0; j<=hiDistr_y_nbin_;j++)
503  {
504 // First moment
505  float cc1=hiDistrMBPl2D_->getBinContent(k,j);
506  cc1 = cc1 * 1./eventCounter_;
507  float cc2=hiDistrNoisePl2D_->getBinContent(k,j);
508  cc2 = cc2 * 1./eventCounter_;
509  float cc3=hiDistrMBMin2D_->getBinContent(k,j);
510  cc3 = cc3 * 1./eventCounter_;
511  float cc4=hiDistrNoiseMin2D_->getBinContent(k,j);
512  cc4 = cc4 * 1./eventCounter_;
513 // Second moment
514  float cc11=hiDistrMB2Pl2D_->getBinContent(k,j);
515  cc11 = cc11 * 1./eventCounter_;
516  hiDistrVarMBPl2D_->setBinContent(k,j,cc11-cc1*cc1);
517  float cc22=hiDistrNoise2Pl2D_->getBinContent(k,j);
518  cc22 = cc22 * 1./eventCounter_;
519  hiDistrVarNoisePl2D_->setBinContent(k,j,cc22-cc2*cc2);
520  float cc33=hiDistrMB2Min2D_->getBinContent(k,j);
521  cc33 = cc33 * 1./eventCounter_;
522  hiDistrVarMBMin2D_->setBinContent(k,j,cc33-cc3*cc3);
523  float cc44=hiDistrNoise2Min2D_->getBinContent(k,j);
524  cc44 = cc44 * 1./eventCounter_;
525  hiDistrVarNoiseMin2D_->setBinContent(k,j,cc44-cc4*cc4);
526  }
527  }
528  }
529 }
530 
531 
#define LogDebug(id)
T getParameter(std::string const &) const
MonitorElement * hiDistrNoisePl2D_
T getUntrackedParameter(std::string const &, T const &) const
bool saveToFile_
Write to file.
int i
Definition: DBlmapReader.cc:9
void setBinContent(int binx, double content)
set content of bin (1-D)
MonitorElement * hiDistrVarMBPl2D_
MonitorElement * hiDistrVarNoiseMin2D_
DQMHcalPhiSymAlCaReco(const edm::ParameterSet &)
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
edm::EDGetTokenT< HBHERecHitCollection > hbherecoNoise
std::vector< HBHERecHit >::const_iterator const_iterator
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)
edm::EDGetTokenT< FEDRawDataCollection > rawInLabel_
void endRun(const edm::Run &r, const edm::EventSetup &c)
int getSpigotData(int nspigot, HcalHTRData &decodeTool, int validSize) const
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:47
return((rh^lh)&mask)
void Fill(long long x)
bool isUnsuppressed() const
Is this event an unsuppresed event?
Definition: HcalHTRData.cc:353
MonitorElement * hiDistrMBPl2D_
MonitorElement * hiDistrHFsize1D_
int depth() const
get the tower depth
Definition: HcalDetId.h:55
edm::EDGetTokenT< HFRecHitCollection > hfrecoNoise
int iEvent
Definition: GenABIO.cc:230
const FEDRawData & FEDData(int fedid) const
retrieve data for fed
int ieta() const
get the cell ieta
Definition: HcalDetId.h:51
MonitorElement * hiDistrVarMBMin2D_
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
int j
Definition: DBlmapReader.cc:9
MonitorElement * hiDistrVarNoisePl2D_
bool getSpigotPresent(unsigned int nspigot) const
Read the &quot;PRESENT&quot; bit for this spigot.
std::string folderName_
DQM folder name.
MonitorElement * hiDistrMB2Pl2D_
bool isValid() const
Definition: HandleBase.h:75
MonitorElement * hiDistrNoise2Min2D_
void analyze(const edm::Event &e, const edm::EventSetup &c)
const_iterator end() const
int iphi() const
get the cell iphi
Definition: HcalDetId.h:53
Definition: DetId.h:18
void beginLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context)
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
T const * product() const
Definition: Handle.h:81
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
MonitorElement * hiDistrMBMin2D_
edm::EDGetTokenT< HFRecHitCollection > hfrecoMB
static const int SPIGOT_COUNT
Definition: HcalDCCHeader.h:19
double getBinContent(int binx) const
get content of bin (1-D)
void endLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
size_type size() const
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:28
MonitorElement * hiDistrNoise2Pl2D_
MonitorElement * hiDistrNoiseMin2D_
std::string fileName_
Output file name if required.
int lvl1ID()
Level-1 event number generated by the TTC system.
Definition: FEDHeader.cc:20
edm::EDGetTokenT< HBHERecHitCollection > hbherecoMB
object to monitor
void setAxisTitle(const std::string &title, int axis=1)
set x-, y- or z-axis title (axis=1, 2, 3 respectively)
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Power< A, B >::type pow(const A &a, const B &b)
Definition: Power.h:40
const_iterator begin() const
Definition: Run.h:41
MonitorElement * hiDistrHBHEsize1D_
MonitorElement * hiDistrMB2Min2D_