CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalDetDiagLaserMonitor.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: HcalDetDiagLaserMonitor
4 // Class: HcalDetDiagLaserMonitor
5 //
13 //
14 // Original Author: Dmitry Vishnevskiy,591 R-013,+41227674265,
15 // Created: Wed Mar 3 12:14:16 CET 2010
16 //
17 //
18 
19 #include "DQM/HcalMonitorTasks/interface/HcalBaseDQMonitor.h"
20 #include "FWCore/Framework/interface/Frameworkfwd.h"
21 #include "FWCore/Framework/interface/Event.h"
22 #include "FWCore/Framework/interface/MakerMacros.h"
23 #include "FWCore/ParameterSet/interface/ParameterSet.h"
24 #include "FWCore/Utilities/interface/InputTag.h"
25 
26 #include "DQMServices/Core/interface/DQMStore.h"
27 #include "DQMServices/Core/interface/MonitorElement.h"
28 #include "DQM/HcalMonitorTasks/interface/HcalBaseMonitor.h"
29 
30 #include "CalibFormats/HcalObjects/interface/HcalDbService.h"
31 #include "CalibFormats/HcalObjects/interface/HcalDbRecord.h"
32 
33 #include <math.h>
34 
35 // this is to retrieve HCAL digi's
36 #include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
37 // to retrive trigger information (local runs only)
38 #include "TBDataFormats/HcalTBObjects/interface/HcalTBTriggerData.h"
39 
40 #include "CalibCalorimetry/HcalAlgos/interface/HcalLogicalMapGenerator.h"
41 #include "CondFormats/HcalObjects/interface/HcalLogicalMap.h"
42 
43 #include "DataFormats/HcalDigi/interface/HcalCalibrationEventTypes.h"
44 #include "DataFormats/FEDRawData/interface/FEDRawData.h"
45 #include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
46 #include "DataFormats/FEDRawData/interface/FEDNumbering.h"
47 #include "DataFormats/HcalDigi/interface/HcalCalibrationEventTypes.h"
48 #include "EventFilter/HcalRawToDigi/interface/HcalDCCHeader.h"
49 
50 #include "TFile.h"
51 #include "TTree.h"
52 #include "TSystem.h"
53 #include "TF1.h"
54 
55 #include "DataFormats/HcalDetId/interface/HcalDetId.h"
56 #include "DQM/HcalMonitorTasks/interface/HcalEtaPhiHists.h"
57 
58 #include <iostream>
59 #include <fstream>
60 
62 static const float adc2fC[128]={-0.5,0.5,1.5,2.5,3.5,4.5,5.5,6.5,7.5,8.5,9.5, 10.5,11.5,12.5,
63  13.5,15.,17.,19.,21.,23.,25.,27.,29.5,32.5,35.5,38.5,42.,46.,50.,54.5,59.5,
64  64.5,59.5,64.5,69.5,74.5,79.5,84.5,89.5,94.5,99.5,104.5,109.5,114.5,119.5,
65  124.5,129.5,137.,147.,157.,167.,177.,187.,197.,209.5,224.5,239.5,254.5,272.,
66  292.,312.,334.5,359.5,384.5,359.5,384.5,409.5,434.5,459.5,484.5,509.5,534.5,
67  559.5,584.5,609.5,634.5,659.5,684.5,709.5,747.,797.,847.,897.,947.,997.,
68  1047.,1109.5,1184.5,1259.5,1334.5,1422.,1522.,1622.,1734.5,1859.5,1984.5,
69  1859.5,1984.5,2109.5,2234.5,2359.5,2484.5,2609.5,2734.5,2859.5,2984.5,
70  3109.5,3234.5,3359.5,3484.5,3609.5,3797.,4047.,4297.,4547.,4797.,5047.,
71  5297.,5609.5,5984.5,6359.5,6734.5,7172.,7672.,8172.,8734.5,9359.5,9984.5};
73 typedef struct{
74 int eta;
75 int phi;
76 int depth;
77 }Raddam_ch;
78 
79 static const Raddam_ch RADDAM_CH[56]={{-30,35,1},{-30,71,1},{-32,15,1},{-32,51,1},{-34,35,1},{-34,71,1},{-36,15,1},
80  {-36,51,1},{-38,35,1},{-38,71,1},{-40,15,1},{-40,51,1},{-41,35,1},{-41,71,1},
81  {30,21,1}, {30,57,1}, {32,1,1}, {32,37,1}, {34,21,1}, {34,57,1}, {36,1,1 },
82  {36,37,1}, {38,21,1}, {38,57,1}, {40,35,1}, {40,71,1}, {41,19,1}, {41,55,1 },
83  {-30,15,2},{-30,51,2},{-32,35,2},{-32,71,2},{-34,15,2},{-34,51,2},{-36,35,2},
84  {-36,71,2},{-38,15,2},{-38,51,2},{-40,35,2},{-40,71,2},{-41,15,2},{-41,51,2},
85  {30,1,2}, {30,37,2}, {32,21,2}, {32,57,2}, {34,1,2}, {34,37,2}, {36,21,2 },
86  {36,57,2}, {38,1,2}, {38,37,2}, {40,19,2}, {40,55,2}, {41,35,2}, {41,71,2}};
87 class HcalRaddamData{
88  public:
89  HcalRaddamData(){
90  for(int i=0;i<128;i++) s1_adc[i]=s2_adc[i]=0;
91  TOTEVNT=CUT1EVNT=CUT2EVNT=S1MEAN=S2MEAN=S1RMS=S2RMS=0;
92  S1FITMEAN=S2FITMEAN=S1FITMEANER=S2FITMEANER=S1FITSIGMA=S2FITSIGMA=0;
93  S1CHI2=S2CHI2=S1NDF=S2NDF=S1BINWIDTH=S2BINWIDTH=0;
94  }
95  ~HcalRaddamData(){};
96  int TOTEVNT;
97  int CUT1EVNT;
98  int CUT2EVNT;
99  float S1MEAN,S2MEAN;
100  float S1RMS,S2RMS;
101  float S1FITMEAN,S2FITMEAN;
102  float S1FITMEANER,S2FITMEANER;
103  float S1FITSIGMA,S2FITSIGMA;
104  float S1CHI2,S2CHI2;
105  float S1NDF,S2NDF;
106  float S1BINWIDTH,S2BINWIDTH;
107  int s1_adc[128];
108  int s2_adc[128];
109 };
110 
111 HcalRaddamData Raddam_data[56];
112 
113 class HcalDetDiagLaserData{
114 public:
115  HcalDetDiagLaserData(){
116  IsRefetence=false;
117  ds_amp=ds_rms=ds_time=ds_time_rms=-100;
118  nChecks=0;nBadTime=0;nBadEnergy=0;
119  status=0;
120  reset();
121  reset1();
122  }
123  void reset(){
124  Xe=XXe=Xt=XXt=n=0;
125  overflow=0;
126  undeflow=0;
127  }
128  void reset1(){
129  Xe1=0;
130  XXe1=0;
131  Xt1=0;
132  XXt1=0;
133  n1=0;
134  overflow1=0;
135  undeflow1=0;
136  }
137  void add_statistics(double *data,int nTS){
138  double e=GetEnergy(data,nTS);
139  double t=GetTime(data,nTS);
140  if(e<20){ undeflow++;undeflow1++; }else if(e>10000){ overflow++;overflow1++; }else{
141  n++; Xe+=e; XXe+=e*e; Xt+=t; XXt+=t*t;
142  n1++; Xe1+=e; XXe1+=e*e; Xt1+=t; XXt1+=t*t;
143  }
144  }
145  void set_reference(float val,float rms,float time,float time_rms){
146  ref_amp=val; ref_rms=rms;
147  ref_time=time; ref_time_rms=time_rms;
148  IsRefetence=true;
149  }
150  void set_data(float val,float rms,float time,float time_rms){
151  ds_amp=val; ds_rms=rms;
152  ds_time=time; ds_time_rms=time_rms;
153  }
154  void change_status(int val){
155  status|=val;
156  }
157  int get_status(){
158  return status;
159  }
160  bool get_reference(double *val,double *rms,double *time,double *time_rms){
161  *val=ref_amp; *rms=ref_rms;
162  *time=ref_time; *time_rms=ref_time_rms;
163  return IsRefetence;
164  }
165  bool get_average_amp(double *ave,double *rms){
166  if(ds_time>-10){ *ave=ds_amp; *rms=ds_rms; return true; }
167  if(n>0){ *ave=Xe/n; *rms=sqrt(XXe/n-(Xe*Xe)/(n*n));} else return false;
168  return true;
169  }
170  bool get_average_time(double *ave,double *rms){
171  if(ds_time>-10){ *ave=ds_time; *rms=ds_time_rms; return true; }
172  if(n>0){ *ave=Xt/n; *rms=sqrt(XXt/n-(Xt*Xt)/(n*n));} else return false;
173  return true;
174  }
175  int get_statistics(){
176  return (int)n;
177  }
178  void set_statistics(int stat){
179  n=stat;
180  }
181  void set_statistics1(int stat){
182  n1=stat;
183  }
184  int get_overflow(){
185  return overflow;
186  }
187  int get_undeflow(){
188  return undeflow;
189  }
190  bool get_average_amp1(double *ave,double *rms){
191  if(ds_time>-10){ *ave=ds_amp; *rms=ds_rms; return true; }
192  if(n1>0){ *ave=Xe1/n1; *rms=sqrt(XXe1/n1-(Xe1*Xe1)/(n1*n1));} else return false;
193  return true;
194  }
195  bool get_average_time1(double *ave,double *rms){
196  if(ds_time>-10){ *ave=ds_time; *rms=ds_time_rms; return true; }
197  if(n1>0){ *ave=Xt1/n1; *rms=sqrt(XXt1/n1-(Xt1*Xt1)/(n1*n1));} else return false;
198  return true;
199  }
200  int get_statistics1(){
201  return (int)n1;
202  }
203  int get_overflow1(){
204  return overflow1;
205  }
206  int get_undeflow1(){
207  return undeflow1;
208  }
209  double GetEnergy(double *data,int n){
210  int MaxI=0; double Energy,MaxE=0;
211  for(int j=0;j<n;++j) if(MaxE<data[j]){ MaxE=data[j]; MaxI=j; }
212  Energy=data[MaxI];
213  if(MaxI>0) Energy+=data[MaxI-1];
214  if(MaxI>1) Energy+=data[MaxI-2];
215  if(MaxI<(n-1)) Energy+=data[MaxI+1];
216  if(MaxI<(n-2)) Energy+=data[MaxI+2];
217  return Energy;
218  }
219  double GetTime(double *data,int n=10){
220  int MaxI=-100; double Time=-9999,SumT=0,MaxT=-10;
221  for(int j=0;j<n;++j) if(MaxT<data[j]){ MaxT=data[j]; MaxI=j; }
222  if (MaxI>=0) // dummy protection so that compiler doesn't think MaxI=-100
223  {
224  Time=MaxI*data[MaxI];
225  SumT=data[MaxI];
226  if(MaxI>0){ Time+=(MaxI-1)*data[MaxI-1]; SumT+=data[MaxI-1]; }
227  if(MaxI<(n-1)){ Time+=(MaxI+1)*data[MaxI+1]; SumT+=data[MaxI+1]; }
228  Time=Time/SumT;
229  }
230  return Time;
231  }
232  int overflow;
233  int undeflow;
234  int overflow1;
235  int undeflow1;
236  double Xe,XXe,Xt,XXt,n;
237  double Xe1,XXe1,Xt1,XXt1,n1;
238  bool IsRefetence;
239  float ref_amp;
240  float ref_rms;
241  float ref_time;
242  float ref_time_rms;
243  float ds_amp;
244  float ds_rms;
245  float ds_time;
246  float ds_time_rms;
247  int status;
248  float nChecks,nBadTime,nBadEnergy;
249 };
250 
251 class HcalDetDiagLaserMonitor : public HcalBaseDQMonitor {
252  public:
253  explicit HcalDetDiagLaserMonitor(const edm::ParameterSet&);
254  ~HcalDetDiagLaserMonitor();
255 
256  private:
257  HcalDetDiagLaserData* GetCalib(std::string sd,int eta,int phi){
258  int SD=0,ETA=0,PHI=0;
259  if(sd.compare("HB")==0) SD=1;
260  if(sd.compare("HE")==0) SD=2;
261  if(sd.compare("HO")==0) SD=3;
262  if(sd.compare("HF")==0) SD=4;
263  if(SD==1 || SD==2){
264  if(eta>0) ETA=1; else ETA=-1;
265  if(phi==71 ||phi==72 || phi==1 || phi==2) PHI=71; else PHI=((phi-3)/4)*4+3;
266  }else if(SD==3){
267  if(abs(eta)<=4){
268  ETA=0;
269  if(phi==71 ||phi==72 || phi==1 || phi==2 || phi==3 || phi==4) PHI=71; else PHI=((phi-5)/6)*6+5;
270  }else{
271  if(abs(eta)>4 && abs(eta)<=10) ETA=1;
272  if(abs(eta)>10 && abs(eta)<=15) ETA=2;
273  if(eta<0) ETA=-ETA;
274  if(phi==71 ||phi==72 || (phi>=1 && phi<=10)) PHI=71; else PHI=((phi-11)/12)*12+11;
275  }
276  }else if(SD==4){
277  if(eta>0) ETA=1; else ETA=-1;
278  if(phi>=1 && phi<=18) PHI=1;
279  if(phi>=19 && phi<=36) PHI=19;
280  if(phi>=37 && phi<=54) PHI=37;
281  if(phi>=55 && phi<=72) PHI=55;
282  }
283  return &calib_data[SD][ETA+2][PHI-1];
284  };
285  void bookHistograms(DQMStore::IBooker &ib, const edm::Run& run, const edm::EventSetup& c) override;
286  void endRun(const edm::Run& run, const edm::EventSetup& c) override;
287  void beginLuminosityBlock(const edm::LuminosityBlock& lumiSeg,const edm::EventSetup& c) override ;
288  void endLuminosityBlock(const edm::LuminosityBlock& lumiSeg,const edm::EventSetup& c) override;
289  void analyze(const edm::Event&, const edm::EventSetup&) override;
290 
291  const HcalElectronicsMap *emap;
292  edm::InputTag inputLabelDigi_;
293 
294  edm::EDGetTokenT<FEDRawDataCollection> tok_raw_;
295  edm::EDGetTokenT<HcalCalibDigiCollection> tok_calib_;
296  edm::EDGetTokenT<HcalTBTriggerData> tok_tb_;
297  edm::EDGetTokenT<HBHEDigiCollection> tok_hbhe_;
298  edm::EDGetTokenT<HODigiCollection> tok_ho_;
299  edm::EDGetTokenT<HFDigiCollection> tok_hf_;
300 
301  void SaveReference();
302  void SaveRaddamData();
303  void LoadReference();
304  void LoadDataset();
305 
306  bool get_ave_rbx(int sd,int side,int rbx,float *ave,float *rms);
307  bool get_ave_subdet(int sd,float *ave_t,float *ave_e,float *ave_t_r,float *ave_e_r);
308  void fillHistos(int sd);
309  void fillProblems(int sd);
310  int nHBHEchecks,nHOchecks,nHFchecks;
311  double LaserTimingThreshold,LaserEnergyThreshold,RaddamThreshold1,RaddamThreshold2;
312 
313  int ievt_;
314  int run_number;
315  int dataset_seq_number;
316  bool IsReference;
317  bool LocalRun,RaddamRun;
318  int nHB,nHE,nHO,nHF;
319  std::string ReferenceData;
320  std::string ReferenceRun;
321  std::string OutputFilePath;
322  std::string XmlFilePath;
323  std::string baseFolder_;
324  std::string prefixME_;
325  bool Online_;
326  bool Overwrite;
327 
328  // to create html from processed dataset
329  std::string DatasetName;
330  std::string htmlOutputPath;
331  bool createHTMLonly;
332 
333  MonitorElement *meEVT_,*meRUN_,*htmlFolder;
334  MonitorElement *RefRun_;
335  MonitorElement *hbheEnergy;
336  MonitorElement *hbheTime;
337  MonitorElement *hbheEnergyRMS;
338  MonitorElement *hbheTimeRMS;
339  MonitorElement *hoEnergy;
340  MonitorElement *hoTime;
341  MonitorElement *hoEnergyRMS;
342  MonitorElement *hoTimeRMS;
343  MonitorElement *hfEnergy;
344  MonitorElement *hfTime;
345  MonitorElement *hfEnergyRMS;
346  MonitorElement *hfTimeRMS;
347 
348  MonitorElement *Time2Dhbhehf;
349  MonitorElement *Time2Dho;
350  MonitorElement *Energy2Dhbhehf;
351  MonitorElement *Energy2Dho;
352  MonitorElement *refTime2Dhbhehf;
353  MonitorElement *refTime2Dho;
354  MonitorElement *refEnergy2Dhbhehf;
355  MonitorElement *refEnergy2Dho;
356 
357  MonitorElement *hb_time_rbx;
358  MonitorElement *he_time_rbx;
359  MonitorElement *ho_time_rbx;
360  MonitorElement *hf_time_rbx;
361 
362  MonitorElement *Raddam[56];
363  TH1F *S1[56],*S2[56];
364 
365  EtaPhiHists* ProblemCellsByDepth_timing;
366  EtaPhiHists* ProblemCellsByDepth_energy;
367  std::vector<std::string> problemnames_;
368 
369  EtaPhiHists* ProblemCellsByDepth_timing_val;
370  EtaPhiHists* ProblemCellsByDepth_energy_val;
371 
372  HcalDetDiagLaserData hb_data[85][72][4];
373  HcalDetDiagLaserData he_data[85][72][4];
374  HcalDetDiagLaserData ho_data[85][72][4];
375  HcalDetDiagLaserData hf_data[85][72][4];
376  HcalDetDiagLaserData calib_data[5][5][72];
377 
378  std::map<unsigned int, int> KnownBadCells_;
379 };
380 
381 HcalDetDiagLaserMonitor::HcalDetDiagLaserMonitor(const edm::ParameterSet& iConfig):
382  HcalBaseDQMonitor(iConfig)
383  {
384 
385 
386  ievt_=-1;
387  emap=0;
388  dataset_seq_number=1;
389  run_number=-1;
390  IsReference=false;
391  LocalRun=RaddamRun=false;
392  createHTMLonly=false;
393  nHB=nHE=nHO=nHF=0;
394  nHBHEchecks=nHOchecks=nHFchecks=0;
395 
396  inputLabelDigi_ = iConfig.getUntrackedParameter<edm::InputTag>("digiLabel",edm::InputTag("hcalDigis"));
397 
398  ReferenceData = iConfig.getUntrackedParameter<std::string>("LaserReferenceData" ,"");
399  OutputFilePath = iConfig.getUntrackedParameter<std::string>("OutputFilePath", "");
400  DatasetName = iConfig.getUntrackedParameter<std::string>("LaserDatasetName", "");
401  htmlOutputPath = iConfig.getUntrackedParameter<std::string>("htmlOutputPath", "");
402  XmlFilePath = iConfig.getUntrackedParameter<std::string>("XmlFilePath", "");
403  Online_ = iConfig.getUntrackedParameter<bool> ("online",false);
404  Overwrite = iConfig.getUntrackedParameter<bool> ("Overwrite",true);
405  prefixME_ = iConfig.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
406  if (prefixME_.size()>0 && prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
407  prefixME_.append("/");
408  subdir_ = iConfig.getUntrackedParameter<std::string>("TaskFolder","DetDiagPedestalMonitor_Hcal/");
409  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
410  subdir_.append("/");
411  subdir_=prefixME_+subdir_;
412  debug_ = iConfig.getUntrackedParameter<int>("debug",0);
413 
414  LaserTimingThreshold = iConfig.getUntrackedParameter<double>("LaserTimingThreshold",0.2);
415  LaserEnergyThreshold = iConfig.getUntrackedParameter<double>("LaserEnergyThreshold",0.1);
416  RaddamThreshold1 = iConfig.getUntrackedParameter<double>("RaddamThreshold1",10.0);
417  RaddamThreshold2 = iConfig.getUntrackedParameter<double>("RaddamThreshold2",0.95);
418 
419  // register for data access
420  tok_tb_ = consumes<HcalTBTriggerData>(iConfig.getParameter<edm::InputTag>("hcalTBTriggerDataTag"));
421  tok_raw_ = consumes<FEDRawDataCollection>(iConfig.getUntrackedParameter<edm::InputTag>("RawDataLabel",edm::InputTag("source")));
422  tok_calib_ = consumes<HcalCalibDigiCollection>(iConfig.getUntrackedParameter<edm::InputTag>("calibDigiLabel",edm::InputTag("hcalDigis")));
423  tok_hbhe_ = consumes<HBHEDigiCollection>(inputLabelDigi_);
424  tok_ho_ = consumes<HODigiCollection>(inputLabelDigi_);
425  tok_hf_ = consumes<HFDigiCollection>(inputLabelDigi_);
426 
427 
428  ProblemCellsByDepth_timing=0;
429  ProblemCellsByDepth_energy=0;
430  ProblemCellsByDepth_timing_val=0;
431  ProblemCellsByDepth_energy_val=0;
432 
433 }
434 
435 
436 void HcalDetDiagLaserMonitor::bookHistograms(DQMStore::IBooker &ib, const edm::Run& run, const edm::EventSetup& c){
437  edm::ESHandle<HcalChannelQuality> p;
438  c.get<HcalChannelQualityRcd>().get("withTopo",p);
439  const HcalChannelQuality* chanquality= p.product();
440  std::vector<DetId> mydetids = chanquality->getAllChannels();
441  KnownBadCells_.clear();
442 
443  for (std::vector<DetId>::const_iterator i = mydetids.begin();i!=mydetids.end();++i){
444  if (i->det()!=DetId::Hcal) continue; // not an hcal cell
445  HcalDetId id=HcalDetId(*i);
446  int status=(chanquality->getValues(id))->getValue();
447  if((status & HcalChannelStatus::HcalCellOff) || (status & HcalChannelStatus::HcalCellMask)){
448  KnownBadCells_[id.rawId()]=status;
449  }
450  }
451 
452  edm::ESHandle<HcalDbService> conditions_;
453  c.get<HcalDbRecord>().get(conditions_);
454  emap=conditions_->getHcalMapping();
455 
456  HcalBaseDQMonitor::setup(ib);
457  std::string name;
458 
459  ib.setCurrentFolder(subdir_);
460  meEVT_ = ib.bookInt("HcalDetDiagLaserMonitor Event Number");
461  meRUN_ = ib.bookInt("HcalDetDiagLaserMonitor Run Number");
462 
463  ReferenceRun="UNKNOWN";
464  LoadReference();
465  LoadDataset();
466  if(DatasetName.size()>0 && createHTMLonly){
467  char str[200]; sprintf(str,"%sHcalDetDiagLaserData_run%i_%i/",htmlOutputPath.c_str(),run_number,dataset_seq_number);
468  htmlFolder=ib.bookString("HcalDetDiagLaserMonitor HTML folder",str);
469  MonitorElement *me;
470  ib.setCurrentFolder(prefixME_+"HcalInfo");
471  me=ib.bookInt("HBpresent");
472  if(nHB>0) me->Fill(1);
473  me=ib.bookInt("HEpresent");
474  if(nHE>0) me->Fill(1);
475  me=ib.bookInt("HOpresent");
476  if(nHO>0) me->Fill(1);
477  me=ib.bookInt("HFpresent");
478  if(nHF>0) me->Fill(1);
479  }
480  ProblemCellsByDepth_timing = new EtaPhiHists();
481  ProblemCellsByDepth_timing->setup(ib," Problem Bad Laser Timing");
482  for(unsigned int i=0;i<ProblemCellsByDepth_timing->depth.size();i++)
483  problemnames_.push_back(ProblemCellsByDepth_timing->depth[i]->getName());
484  ProblemCellsByDepth_energy = new EtaPhiHists();
485  ProblemCellsByDepth_energy->setup(ib," Problem Bad Laser Energy");
486  for(unsigned int i=0;i<ProblemCellsByDepth_energy->depth.size();i++)
487  problemnames_.push_back(ProblemCellsByDepth_energy->depth[i]->getName());
488 
489  ib.setCurrentFolder(subdir_+"Summary Plots");
490 
491  name="HBHE Laser Energy Distribution"; hbheEnergy = ib.book1D(name,name,200,0,3000);
492  name="HBHE Laser Timing Distribution"; hbheTime = ib.book1D(name,name,200,0,10);
493  name="HBHE Laser Energy RMS_div_Energy Distribution"; hbheEnergyRMS = ib.book1D(name,name,200,0,0.5);
494  name="HBHE Laser Timing RMS Distribution"; hbheTimeRMS = ib.book1D(name,name,200,0,1);
495  name="HO Laser Energy Distribution"; hoEnergy = ib.book1D(name,name,200,0,3000);
496  name="HO Laser Timing Distribution"; hoTime = ib.book1D(name,name,200,0,10);
497  name="HO Laser Energy RMS_div_Energy Distribution"; hoEnergyRMS = ib.book1D(name,name,200,0,0.5);
498  name="HO Laser Timing RMS Distribution"; hoTimeRMS = ib.book1D(name,name,200,0,1);
499  name="HF Laser Energy Distribution"; hfEnergy = ib.book1D(name,name,200,0,3000);
500  name="HF Laser Timing Distribution"; hfTime = ib.book1D(name,name,200,0,10);
501  name="HF Laser Energy RMS_div_Energy Distribution"; hfEnergyRMS = ib.book1D(name,name,200,0,0.7);
502  name="HF Laser Timing RMS Distribution"; hfTimeRMS = ib.book1D(name,name,200,0,1);
503 
504  name="Laser Timing HBHEHF"; Time2Dhbhehf = ib.book2D(name,name,87,-43,43,74,0,73);
505  name="Laser Timing HO"; Time2Dho = ib.book2D(name,name,33,-16,16,74,0,73);
506  name="Laser Energy HBHEHF"; Energy2Dhbhehf = ib.book2D(name,name,87,-43,43,74,0,73);
507  name="Laser Energy HO"; Energy2Dho = ib.book2D(name,name,33,-16,16,74,0,73);
508  name="HBHEHF Laser (Timing-Ref)+1"; refTime2Dhbhehf = ib.book2D(name,name,87,-43,43,74,0,73);
509  name="HO Laser (Timing-Ref)+1"; refTime2Dho = ib.book2D(name,name,33,-16,16,74,0,73);
510  name="HBHEHF Laser Energy_div_Ref"; refEnergy2Dhbhehf = ib.book2D(name,name,87,-43,43,74,0,73);
511  name="HO Laser Energy_div_Ref"; refEnergy2Dho = ib.book2D(name,name,33,-16,16,74,0,73);
512 
513  name="HB RBX average Time-Ref"; hb_time_rbx = ib.book1D(name,name,36,0.5,36.5);
514  name="HE RBX average Time-Ref"; he_time_rbx = ib.book1D(name,name,36,0.5,36.5);
515  name="HO RBX average Time-Ref"; ho_time_rbx = ib.book1D(name,name,36,0.5,36.5);
516  name="HF RoBox average Time-Ref"; hf_time_rbx = ib.book1D(name,name,24,0.5,24.5);
517 
518  char str[200];
519  for(int i=1;i<=18;i++){ sprintf(str,"HBM%02i",i); hb_time_rbx->setBinLabel(i,str); }
520  for(int i=1;i<=18;i++){ sprintf(str,"HBP%02i",i); hb_time_rbx->setBinLabel(i+18,str); }
521  for(int i=1;i<=18;i++){ sprintf(str,"HEM%02i",i); he_time_rbx->setBinLabel(i,str); }
522  for(int i=1;i<=18;i++){ sprintf(str,"HEP%02i",i); he_time_rbx->setBinLabel(i+18,str); }
523  for(int i=1;i<=12;i++){ sprintf(str,"HFM%02i",i); hf_time_rbx->setBinLabel(i,str); }
524  for(int i=1;i<=12;i++){ sprintf(str,"HFP%02i",i); hf_time_rbx->setBinLabel(i+12,str); }
525  for(int i=1;i<=6;i++){ sprintf(str,"HO2M%02i",i*2); ho_time_rbx->setBinLabel(i,str); }
526  for(int i=1;i<=6;i++){ sprintf(str,"HO1M%02i",i*2); ho_time_rbx->setBinLabel(i+6,str); }
527  for(int i=1;i<=12;i++){ sprintf(str,"HO0%02i",i); ho_time_rbx->setBinLabel(i+12,str); }
528  for(int i=1;i<=6;i++){ sprintf(str,"HO1P%02i",i*2); ho_time_rbx->setBinLabel(i+24,str); }
529  for(int i=1;i<=6;i++){ sprintf(str,"HO2P%02i",i*2); ho_time_rbx->setBinLabel(i+30,str); }
530 
531  Time2Dhbhehf->setAxisTitle("i#eta",1);
532  Time2Dhbhehf->setAxisTitle("i#phi",2);
533  Time2Dho->setAxisTitle("i#eta",1);
534  Time2Dho->setAxisTitle("i#phi",2);
535  Energy2Dhbhehf->setAxisTitle("i#eta",1);
536  Energy2Dhbhehf->setAxisTitle("i#phi",2);
537  Energy2Dho->setAxisTitle("i#eta",1);
538  Energy2Dho->setAxisTitle("i#phi",2);
539  refTime2Dhbhehf->setAxisTitle("i#eta",1);
540  refTime2Dhbhehf->setAxisTitle("i#phi",2);
541  refTime2Dho->setAxisTitle("i#eta",1);
542  refTime2Dho->setAxisTitle("i#phi",2);
543  refEnergy2Dhbhehf->setAxisTitle("i#eta",1);
544  refEnergy2Dhbhehf->setAxisTitle("i#phi",2);
545  refEnergy2Dho->setAxisTitle("i#eta",1);
546  refEnergy2Dho->setAxisTitle("i#phi",2);
547 
548  refTime2Dhbhehf->setAxisRange(0,2,3);
549  refTime2Dho->setAxisRange(0,2,3);
550  refEnergy2Dhbhehf->setAxisRange(0.5,1.5,3);
551  refEnergy2Dho->setAxisRange(0.5,1.5,3);
552 
553  ib.setCurrentFolder(subdir_);
554  RefRun_= ib.bookString("HcalDetDiagLaserMonitor Reference Run",ReferenceRun);
555 
556  ib.setCurrentFolder(subdir_+"Raddam Plots");
557  for(int i=0;i<56;i++){
558  sprintf(str,"RADDAM (%i %i)",RADDAM_CH[i].eta,RADDAM_CH[i].phi);
559  Raddam[i] = ib.book1D(str,str,10,-0.5,9.5);
560  }
561 
562  ib.setCurrentFolder(subdir_+"Plots for client");
563  ProblemCellsByDepth_timing_val = new EtaPhiHists();
564  ProblemCellsByDepth_timing_val->setup(ib," Laser Timing difference");
565  ProblemCellsByDepth_energy_val = new EtaPhiHists();
566  ProblemCellsByDepth_energy_val->setup(ib," Laser Energy difference");
567 }
568 
569 
570 HcalDetDiagLaserMonitor::~HcalDetDiagLaserMonitor(){
571 
572  if ( ProblemCellsByDepth_timing ) delete ProblemCellsByDepth_timing;
573  if ( ProblemCellsByDepth_energy ) delete ProblemCellsByDepth_energy;
574  if ( ProblemCellsByDepth_timing_val ) delete ProblemCellsByDepth_timing_val;
575  if ( ProblemCellsByDepth_energy_val ) delete ProblemCellsByDepth_energy_val;
576 
577 }
578 
579 // ------------ method called to for each event ------------
580 void HcalDetDiagLaserMonitor::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup){
581  if(createHTMLonly) return;
582  HcalBaseDQMonitor::analyze(iEvent,iSetup); // base class increments ievt_, etc. counters
583 
584 int eta,phi,depth,nTS;
585 static bool HBHEseq,HOseq,HFseq;
586 static int lastHBHEorbit,lastHOorbit,lastHForbit,nChecksHBHE,nChecksHO,nChecksHF,ievt_hbhe,ievt_ho,ievt_hf;
587  if(ievt_==-1){
588  ievt_=0;HBHEseq=HOseq=HFseq=false; lastHBHEorbit=lastHOorbit=lastHForbit=-1;nChecksHBHE=nChecksHO=nChecksHF=0;
589  ievt_hbhe=0,ievt_ho=0,ievt_hf=0;
590  }
591 
592  bool LaserEvent=false;
593  bool LaserRaddam=false;
594  int orbit=iEvent.orbitNumber();
595  meRUN_->Fill(iEvent.id().run());
596  // for local runs
597  edm::Handle<HcalTBTriggerData> trigger_data;
598  iEvent.getByToken(tok_tb_, trigger_data);
599  if(trigger_data.isValid()){
600  if(trigger_data->wasLaserTrigger()) LaserEvent=true;
601  LocalRun=true;
602  }
603  if(!LocalRun && Online_){
604  if(HBHEseq && (orbit-lastHBHEorbit)>(11223*10) && ievt_hbhe>40){
605  HBHEseq=false;
606  fillHistos(HcalBarrel);
607  fillProblems(HcalBarrel);
608  fillProblems(HcalEndcap);
609  nChecksHBHE++;
610  ievt_hbhe=0;
611  for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) hb_data[i][j][k].reset();
612  for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) he_data[i][j][k].reset();
613  }
614  if(HOseq && (orbit-lastHOorbit)>(11223*10) && ievt_ho>40){
615  HOseq=false;
616  fillHistos(HcalOuter);
617  fillProblems(HcalOuter);
618  nChecksHO++;
619  ievt_ho=0;
620  for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) ho_data[i][j][k].reset();
621  }
622  if(HFseq && (orbit-lastHForbit)>(11223*10) && ievt_hf>40){
623  HFseq=false;
624  fillHistos(HcalForward);
625  fillProblems(HcalForward);
626  nChecksHF++;
627  ievt_hf=0;
628  if(nChecksHF==1 || (nChecksHF>1 && ((nChecksHF-1)%12)==0)){
629  SaveReference();
630  }
631  for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) hf_data[i][j][k].reset();
632  }
633  }
634 
635  // Abort Gap laser
636  if(LocalRun==false || LaserEvent==false){
637  edm::Handle<FEDRawDataCollection> rawdata;
638  iEvent.getByToken(tok_raw_ ,rawdata);
639  // edm::Handle<FEDRawDataCollection> rawdata;
640  // iEvent.getByType(rawdata);
641  //checking FEDs for calibration information
642  for (int i=FEDNumbering::MINHCALFEDID;i<=FEDNumbering::MAXHCALFEDID; i++) {
643  const FEDRawData& fedData = rawdata->FEDData(i) ;
644  if ( fedData.size() < 24 ) continue ;
645  int value = ((const HcalDCCHeader*)(fedData.data()))->getCalibType() ;
646  if(value==hc_HBHEHPD){ HBHEseq=true; HOseq=HFseq=false; lastHBHEorbit=orbit; ievt_hbhe++; }
647  if(value==hc_HOHPD){ HOseq=true; HBHEseq=HFseq=false; lastHOorbit=orbit; ievt_ho++; }
648  if(value==hc_HFPMT){ HFseq=true; HBHEseq=HOseq=false; lastHForbit=orbit; ievt_hf++; }
649 
650  if(value==hc_HBHEHPD || value==hc_HOHPD || value==hc_HFPMT){ LaserEvent=true; break;}
651  if(value==hc_RADDAM){ LaserEvent=true; LaserRaddam=true; break;}
652  }
653  }
654  if(!LaserEvent) return;
655  edm::Handle<HBHEDigiCollection> hbhe;
656  iEvent.getByToken(tok_hbhe_,hbhe);
657  edm::Handle<HODigiCollection> ho;
658  iEvent.getByToken(tok_ho_,ho);
659  edm::Handle<HFDigiCollection> hf;
660  iEvent.getByToken(tok_hf_,hf);
661  edm::Handle<HcalCalibDigiCollection> calib;
662  iEvent.getByToken(tok_calib_, calib);
663 
664  if(LocalRun && LaserEvent){
665  int N=0;
666  if(hf.isValid()){
667  for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
668  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth();
669  float e=0;
670  for(int i=0;i<digi->size();i++) e+=adc2fC[digi->sample(i).adc()&0xff]-2.5;
671  if(e>40){ N++;}
672  }
673  }
674  if(N>50 && N<57){ RaddamRun=true; /*LaserRaddam=true;*/}
675  }
676  if(RaddamRun){
677  if(hf.isValid()){
678  for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
679  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
680  int N;
681  for(N=0;N<56;N++)if(eta==RADDAM_CH[N].eta && phi==RADDAM_CH[N].phi && depth==RADDAM_CH[N].depth) break;
682  if(N==56) continue;
683  float max1=0,max2=0;
684  int nmax1=0,nmax2=0;
685  for(int i=0;i<nTS;i++){
686  if(max1<adc2fC[digi->sample(i).adc()&0xff]){ nmax1=i; max1=adc2fC[digi->sample(i).adc()&0xff]; }
687  }
688  Raddam_data[N].TOTEVNT++;
689  for(int i=0;i<nTS;i++){
690  if(i==nmax1) continue;
691  if(max2<adc2fC[digi->sample(i).adc()&0xff]){ nmax2=i; max2=adc2fC[digi->sample(i).adc()&0xff]; }
692  }
693  if(nmax1>nmax2){
694  int tmp1=nmax2;
695  nmax2=nmax1;nmax1=tmp1;
696  }
697  if(nmax1==0 || nmax2==(nTS-1)) continue;
698  if(nmax2!=(nmax1+1)) continue;
699 
700  if(max1<RaddamThreshold1 || max2<RaddamThreshold1) continue;
701  Raddam_data[N].CUT1EVNT++;
702  max1-=2.5; max2-=2.5;
703  float S2=max1+max2;
704  float S4=S2+adc2fC[digi->sample(nmax1-1).adc()&0xff]+adc2fC[digi->sample(nmax2+1).adc()&0xff]-5.0;
705  if((S2/S4)<RaddamThreshold2) continue;
706  Raddam_data[N].CUT2EVNT++;
707  Raddam_data[N].s1_adc[digi->sample(nmax1).adc()&0xff]++;
708  Raddam_data[N].s2_adc[digi->sample(nmax2).adc()&0xff]++;
709  }
710  }
711  }
712 
713  meEVT_->Fill(++ievt_);
714  run_number=iEvent.id().run();
715  double data[20];
716  if(!LaserRaddam){
717  if(hbhe.isValid()){
718  for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
719  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
720  double ped=(adc2fC[digi->sample(0).adc()&0xff]+adc2fC[digi->sample(1).adc()&0xff])/2.0;
721  if(digi->id().subdet()==HcalBarrel){
722  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
723  hb_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
724  }
725  if(digi->id().subdet()==HcalEndcap){
726  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
727  he_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
728  }
729  }
730  }
731  if(ho.isValid()){
732  for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
733  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
734  double ped=(adc2fC[digi->sample(0).adc()&0xff]+adc2fC[digi->sample(1).adc()&0xff])/2.0;
735  if((eta>=11 && eta<=15 && phi>=59 && phi<=70) || (eta>=5 && eta<=10 && phi>=47 && phi<=58)){
736  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
737  }else{
738  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
739  }
740  ho_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
741  }
742  }
743  if(hf.isValid()){
744  for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
745  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
746  double ped=adc2fC[digi->sample(0).adc()&0xff];
747  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()&0xff]-ped;
748  hf_data[eta+42][phi-1][depth-1].add_statistics(data,nTS);
749  }
750  }
751  if(calib.isValid())for(HcalCalibDigiCollection::const_iterator digi=calib->begin();digi!=calib->end();digi++){
752  if(digi->id().cboxChannel()!=0 || digi->id().hcalSubdet()==0) continue;
753  nTS=digi->size();
754  float e=0;
755  for(int i=0;i<nTS;i++){ data[i]=adc2fC[digi->sample(i).adc()&0xff]; e+=data[i];}
756  if(e<15000) calib_data[digi->id().hcalSubdet()][digi->id().ieta()+2][digi->id().iphi()-1].add_statistics(data,nTS);
757  }
758  }else{ //Raddam
759  if(hf.isValid()){
760  for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
761  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
762  int N;
763  for(N=0;N<56;N++)if(eta==RADDAM_CH[N].eta && phi==RADDAM_CH[N].phi && depth==RADDAM_CH[N].depth) break;
764  if(N==56) continue;
765  for(int i=0;i<nTS;i++) Raddam[N]->Fill(i,adc2fC[digi->sample(i).adc()&0xff]-2.5);
766 
767  }
768  }
769  }
770 }
771 bool HcalDetDiagLaserMonitor::get_ave_subdet(int sd,float *ave_t,float *ave_e,float *ave_t_r,float *ave_e_r){
772 double T=0,nT=0,E=0,nE=0,Tr=0,nTr=0,Er=0,nEr=0;
773  if(sd==HcalBarrel) for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++){
774  for(int depth=1;depth<=2;depth++){
775  if(hb_data[eta+42][phi-1][depth-1].get_statistics()>10){
776  double ave=0,rms=0,time=0,time_rms=0;
777  hb_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
778  hb_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
779  T+=time; nT++; E+=ave; nE++;
780  if(hb_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms,&time,&time_rms)){
781  Tr+=time; nTr++; Er+=ave; nEr++;}
782  }
783  }
784  }
785  // HE histograms
786  if(sd==HcalEndcap) for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++){
787  for(int depth=1;depth<=3;depth++){
788  if(he_data[eta+42][phi-1][depth-1].get_statistics()>10){
789  double ave=0; double rms=0; double time=0; double time_rms=0;
790  he_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
791  he_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
792  T+=time; nT++; E+=ave; nE++;
793  if(he_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms,&time,&time_rms)){
794  Tr+=time; nTr++; Er+=ave; nEr++;}
795  }
796  }
797  }
798  // HF histograms
799  if(sd==HcalForward) for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++){
800  for(int depth=1;depth<=2;depth++){
801  if(hf_data[eta+42][phi-1][depth-1].get_statistics()>10){
802  double ave=0; double rms=0; double time=0; double time_rms=0;
803  hf_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
804  hf_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
805  T+=time; nT++; E+=ave; nE++;
806  if(hf_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms,&time,&time_rms)){
807  Tr+=time; nTr++; Er+=ave; nEr++;}
808  }
809  }
810  }
811  // HO histograms
812  if(sd==HcalOuter) for(int eta=-15;eta<=15;eta++) for(int phi=1;phi<=72;phi++){
813  for(int depth=4;depth<=4;depth++){
814  if(ho_data[eta+42][phi-1][depth-1].get_statistics()>10){
815  double ave=0; double rms=0; double time=0; double time_rms=0;
816  ho_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
817  ho_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
818  T+=time; nT++; E+=ave; nE++;
819  if(ho_data[eta+42][phi-1][depth-1].get_reference(&ave,&rms,&time,&time_rms)){
820  Tr+=time; nTr++; Er+=ave; nEr++;}
821  }
822  }
823  }
824  if(nT<200 || nE<200 || nTr<200 || nEr<200) return false;
825  *ave_t=T/nT;
826  *ave_e=E/nE;
827  *ave_t_r=Tr/nTr;
828  *ave_e_r=Er/nEr;
829  return true;
830 }
831 
832 void HcalDetDiagLaserMonitor::fillProblems(int sd){
833 float ave_t,ave_e,ave_t_r,ave_e_r;
834  if(!get_ave_subdet(sd,&ave_t,&ave_e,&ave_t_r,&ave_e_r)) return;
835 
836  for(int i=0;i<4;i++){
837  ProblemCellsByDepth_energy->depth[i]->Reset();
838  ProblemCellsByDepth_timing->depth[i]->Reset();
839  }
840 
841  std::vector <HcalElectronicsId> AllElIds = emap->allElectronicsIdPrecision();
842  for (std::vector <HcalElectronicsId>::iterator eid = AllElIds.begin(); eid != AllElIds.end(); eid++){
843  DetId detid=emap->lookup(*eid);
844  if (detid.det()!=DetId::Hcal) continue;
845  HcalGenericDetId gid(emap->lookup(*eid));
846  if(!(!(gid.null()) &&
847  (gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel ||
848  gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap ||
849  gid.genericSubdet()==HcalGenericDetId::HcalGenForward ||
850  gid.genericSubdet()==HcalGenericDetId::HcalGenOuter))) continue;
851  int eta=0,phi=0,depth=0;
852  HcalDetId hid(detid);
853  if(KnownBadCells_.find(hid.rawId())==KnownBadCells_.end()) continue;
854  eta=hid.ieta();
855  phi=hid.iphi();
856  depth=hid.depth();
857 
858  int e=CalcEtaBin(sd,eta,depth)+1;
859  if(detid.subdetId()==HcalBarrel && sd==HcalBarrel){
860  double val=0,rms=0,time=0,time_rms=0,VAL=0,RMS=0,TIME=0,TIME_RMS=0;
861  if(!hb_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
862  if(!hb_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
863  if(!hb_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
864  hb_data[eta+42][phi-1][depth-1].nChecks++;
865  float diff_t=(TIME-ave_t)-(time-ave_t_r); if(diff_t<0) diff_t=-diff_t;
866  if(diff_t>LaserTimingThreshold){
867  hb_data[eta+42][phi-1][depth-1].nBadTime++;
868  ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,(TIME-ave_t)-(time-ave_t_r));
869  }else ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,0);
870  if(VAL!=0 && val!=0 && ave_e!=0 && ave_e_r!=0){
871  float diff_e=((VAL/ave_e))/(val/ave_e_r);
872  if(diff_e>(1+LaserEnergyThreshold) ||diff_e<(1-LaserEnergyThreshold) ){
873  hb_data[eta+42][phi-1][depth-1].nBadEnergy++;
874  ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,((VAL/ave_e))/(val/ave_e_r));
875  }else ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,0);
876  }
877  }
878  if(detid.subdetId()==HcalEndcap && sd==HcalEndcap){
879  double val=0,rms=0,time=0,time_rms=0,VAL=0,RMS=0,TIME=0,TIME_RMS=0;
880  if(!he_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
881  if(!he_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
882  if(!he_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
883  he_data[eta+42][phi-1][depth-1].nChecks++;
884  float diff_t=(TIME-ave_t)-(time-ave_t_r); if(diff_t<0) diff_t=-diff_t;
885  if(diff_t>LaserTimingThreshold){
886  he_data[eta+42][phi-1][depth-1].nBadTime++;
887  ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,(TIME-ave_t)-(time-ave_t_r));
888  }else ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,0);
889  if(VAL!=0 && val!=0 && ave_e!=0 && ave_e_r!=0){
890  float diff_e=((VAL/ave_e))/(val/ave_e_r);
891  if(diff_e>(1+LaserEnergyThreshold) ||diff_e<(1-LaserEnergyThreshold) ){
892  he_data[eta+42][phi-1][depth-1].nBadEnergy++;
893  ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,((VAL/ave_e))/(val/ave_e_r));
894  }else ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,0);
895  }
896  }
897  if(detid.subdetId()==HcalOuter && sd==HcalOuter){
898  double val=0,rms=0,time=0,time_rms=0,VAL=0,RMS=0,TIME=0,TIME_RMS=0;
899  if(!ho_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
900  if(!ho_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
901  if(!ho_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
902  ho_data[eta+42][phi-1][depth-1].nChecks++;
903  float diff_t=(TIME-ave_t)-(time-ave_t_r); if(diff_t<0) diff_t=-diff_t;
904  if(diff_t>LaserTimingThreshold){
905  ho_data[eta+42][phi-1][depth-1].nBadTime++;
906  ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,(TIME-ave_t)-(time-ave_t_r));
907  }else ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,0);
908  if(VAL!=0 && val!=0 && ave_e!=0 && ave_e_r!=0){
909  float diff_e=((VAL/ave_e))/(val/ave_e_r);
910  if(diff_e>(1+LaserEnergyThreshold) ||diff_e<(1-LaserEnergyThreshold) ){
911  ho_data[eta+42][phi-1][depth-1].nBadEnergy++;
912  ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,((VAL/ave_e))/(val/ave_e_r));
913  }else ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,0);
914  }
915  }
916  if(detid.subdetId()==HcalForward && sd==HcalForward){
917  double val=0,rms=0,time=0,time_rms=0,VAL=0,RMS=0,TIME=0,TIME_RMS=0;
918  if(!hf_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
919  if(!hf_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
920  if(!hf_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
921  hf_data[eta+42][phi-1][depth-1].nChecks++;
922  float diff_t=(TIME-ave_t)-(time-ave_t_r); if(diff_t<0) diff_t=-diff_t;
923  if(diff_t>LaserTimingThreshold){
924  hf_data[eta+42][phi-1][depth-1].nBadTime++;
925  ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,(TIME-ave_t)-(time-ave_t_r));
926  }else ProblemCellsByDepth_timing_val->depth[depth-1]->setBinContent(e,phi,0);
927  if(VAL!=0 && val!=0 && ave_e!=0 && ave_e_r!=0){
928  float diff_e=((VAL/ave_e))/(val/ave_e_r);
929  if(diff_e>(1+LaserEnergyThreshold) ||diff_e<(1-LaserEnergyThreshold) ){
930  hf_data[eta+42][phi-1][depth-1].nBadEnergy++;
931  ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,((VAL/ave_e))/(val/ave_e_r));
932  }else ProblemCellsByDepth_energy_val->depth[depth-1]->setBinContent(e,phi,0);
933  }
934  }
935  }
936  for (std::vector <HcalElectronicsId>::iterator eid = AllElIds.begin(); eid != AllElIds.end(); eid++){
937  DetId detid=emap->lookup(*eid);
938  if (detid.det()!=DetId::Hcal) continue;
939  HcalGenericDetId gid(emap->lookup(*eid));
940  if(!(!(gid.null()) &&
941  (gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel ||
942  gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap ||
943  gid.genericSubdet()==HcalGenericDetId::HcalGenForward ||
944  gid.genericSubdet()==HcalGenericDetId::HcalGenOuter))) continue;
945  int eta=0,phi=0,depth=0;
946  HcalDetId hid(detid);
947  eta=hid.ieta();
948  phi=hid.iphi();
949  depth=hid.depth();
950 
951  if(detid.subdetId()==HcalBarrel){
952  if(hb_data[eta+42][phi-1][depth-1].nBadTime>0){
953  int e=CalcEtaBin(HcalBarrel,eta,depth)+1;
954  double val=hb_data[eta+42][phi-1][depth-1].nBadTime/hb_data[eta+42][phi-1][depth-1].nChecks;
955  ProblemCellsByDepth_timing->depth[depth-1]->setBinContent(e,phi,val);
956  }
957  if(hb_data[eta+42][phi-1][depth-1].nBadEnergy>0){
958  int e=CalcEtaBin(HcalBarrel,eta,depth)+1;
959  double val=hb_data[eta+42][phi-1][depth-1].nBadEnergy/hb_data[eta+42][phi-1][depth-1].nChecks;
960  ProblemCellsByDepth_energy->depth[depth-1]->setBinContent(e,phi,val);
961  }
962  }
963  if(detid.subdetId()==HcalEndcap){
964  if(he_data[eta+42][phi-1][depth-1].nBadTime>0){
965  int e=CalcEtaBin(HcalEndcap,eta,depth)+1;
966  double val=he_data[eta+42][phi-1][depth-1].nBadTime/he_data[eta+42][phi-1][depth-1].nChecks;
967  ProblemCellsByDepth_timing->depth[depth-1]->setBinContent(e,phi,val);
968  }
969  if(he_data[eta+42][phi-1][depth-1].nBadEnergy>0){
970  int e=CalcEtaBin(HcalEndcap,eta,depth)+1;
971  double val=he_data[eta+42][phi-1][depth-1].nBadEnergy/he_data[eta+42][phi-1][depth-1].nChecks;
972  ProblemCellsByDepth_energy->depth[depth-1]->setBinContent(e,phi,val);
973  }
974  }
975  if(detid.subdetId()==HcalOuter){
976  if(ho_data[eta+42][phi-1][depth-1].nBadTime>0){
977  int e=CalcEtaBin(HcalOuter,eta,depth)+1;
978  double val=ho_data[eta+42][phi-1][depth-1].nBadTime/ho_data[eta+42][phi-1][depth-1].nChecks;
979  ProblemCellsByDepth_timing->depth[depth-1]->setBinContent(e,phi,val);
980  }
981  if(ho_data[eta+42][phi-1][depth-1].nBadEnergy>0){
982  int e=CalcEtaBin(HcalOuter,eta,depth)+1;
983  double val=ho_data[eta+42][phi-1][depth-1].nBadEnergy/ho_data[eta+42][phi-1][depth-1].nChecks;
984  ProblemCellsByDepth_energy->depth[depth-1]->setBinContent(e,phi,val);
985  }
986  }
987  if(detid.subdetId()==HcalForward){
988  if(hf_data[eta+42][phi-1][depth-1].nBadTime>0){
989  int e=CalcEtaBin(HcalForward,eta,depth)+1;
990  double val=hf_data[eta+42][phi-1][depth-1].nBadTime/hf_data[eta+42][phi-1][depth-1].nChecks;
991  ProblemCellsByDepth_timing->depth[depth-1]->setBinContent(e,phi,val);
992  }
993  if(hf_data[eta+42][phi-1][depth-1].nBadEnergy>0){
994  int e=CalcEtaBin(HcalForward,eta,depth)+1;
995  double val=hf_data[eta+42][phi-1][depth-1].nBadEnergy/hf_data[eta+42][phi-1][depth-1].nChecks;
996  ProblemCellsByDepth_energy->depth[depth-1]->setBinContent(e,phi,val);
997  }
998  }
999  }
1000 }
1001 
1002 bool HcalDetDiagLaserMonitor::get_ave_rbx(int sd,int side,int rbx,float *ave,float *rms){
1003  double xt=0,xxt=0;
1004  int eta_min=0,eta_max=0,n=0;
1005  if(sd==HcalBarrel){
1006  if(side>0){eta_min=1; eta_max=29;}
1007  if(side<0){eta_min=-29; eta_max=-1;}
1008  if(rbx==1){
1009  for(int i=eta_min;i<=eta_max;i++) for(int j=71;j<=72;j++)for(int k=1;k<=3;k++){
1010  double val,rms,time,time_rms;
1011  double TIME,TIME_RMS;
1012  if(!hb_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1013  if(!hb_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
1014  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1015  }
1016  for(int i=eta_min;i<=eta_max;i++) for(int j=1;j<=2;j++)for(int k=1;k<=3;k++){
1017  double val,rms,time,time_rms;
1018  double TIME,TIME_RMS;
1019  if(!hb_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1020  if(!hb_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
1021  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1022  }
1023 
1024  }else{
1025  for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*4-1);j<=((rbx-1)*4+2);j++)for(int k=1;k<=3;k++){
1026  double val,rms,time,time_rms;
1027  double TIME,TIME_RMS;
1028  if(!hb_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1029  if(!hb_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
1030  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1031  }
1032  }
1033  }
1034  if(sd==HcalEndcap){
1035  if(side>0){eta_min=1; eta_max=29;}
1036  if(side<0){eta_min=-29; eta_max=-1;}
1037  if(rbx==1){
1038  for(int i=eta_min;i<=eta_max;i++) for(int j=71;j<=72;j++)for(int k=1;k<=3;k++){
1039  double val,rms,time,time_rms;
1040  double TIME,TIME_RMS;
1041  if(!he_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1042  if(!he_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
1043  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1044  }
1045  for(int i=eta_min;i<=eta_max;i++) for(int j=1;j<=2;j++)for(int k=1;k<=3;k++){
1046  double val,rms,time,time_rms;
1047  double TIME,TIME_RMS;
1048  if(!he_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1049  if(!he_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
1050  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1051  }
1052 
1053  }else{
1054  for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*4-1);j<=((rbx-1)*4+2);j++)for(int k=1;k<=3;k++){
1055  double val,rms,time,time_rms;
1056  double TIME,TIME_RMS;
1057  if(!he_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1058  if(!he_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
1059  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1060  }
1061  }
1062  }
1063  if(sd==HcalForward){
1064  if(side>0){eta_min=29; eta_max=40;}
1065  if(side<0){eta_min=-40; eta_max=-29;}
1066  for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*6+1);j<=((rbx-1)*6+6);j++)for(int k=1;k<=2;k++){
1067  double val,rms,time,time_rms;
1068  double TIME,TIME_RMS;
1069  if(!hf_data[i+42][j-1][k-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1070  if(!hf_data[i+42][j-1][k-1].get_average_time(&TIME,&TIME_RMS)) continue;
1071  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1072  }
1073  }
1074  if(sd==HcalOuter){
1075  if(side==0){
1076  eta_min=-4,eta_max=4;
1077  if(rbx==1){
1078  for(int i=eta_min;i<=eta_max;i++) for(int j=71;j<=72;j++){
1079  double val,rms,time,time_rms;
1080  double TIME,TIME_RMS;
1081  if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1082  if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
1083  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1084  }
1085  for(int i=eta_min;i<=eta_max;i++) for(int j=1;j<=4;j++){
1086  double val,rms,time,time_rms;
1087  double TIME,TIME_RMS;
1088  if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1089  if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
1090  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1091  }
1092 
1093  }else{
1094  for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*6-1);j<=((rbx-1)*6+4);j++){
1095  double val,rms,time,time_rms;
1096  double TIME,TIME_RMS;
1097  if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1098  if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
1099  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1100  }
1101  }
1102  }
1103  if(side==-1){ eta_min=-10,eta_max=-5;}
1104  if(side==-2){ eta_min=-15,eta_max=-11;}
1105  if(side==1) { eta_min=5, eta_max=10;}
1106  if(side==2) { eta_min=11, eta_max=15;}
1107  if(side!=0){
1108  if(rbx==1){
1109  for(int i=eta_min;i<=eta_max;i++) for(int j=71;j<=72;j++){
1110  double val,rms,time,time_rms;
1111  double TIME,TIME_RMS;
1112  if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1113  if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
1114  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1115  }
1116  for(int i=eta_min;i<=eta_max;i++) for(int j=1;j<=10;j++){
1117  double val,rms,time,time_rms;
1118  double TIME,TIME_RMS;
1119  if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1120  if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
1121  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1122  }
1123 
1124  }else{
1125  for(int i=eta_min;i<=eta_max;i++) for(int j=((rbx-1)*12-1);j<=((rbx-1)*12+10);j++){
1126  double val,rms,time,time_rms;
1127  double TIME,TIME_RMS;
1128  if(!ho_data[i+42][j-1][4-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1129  if(!ho_data[i+42][j-1][4-1].get_average_time(&TIME,&TIME_RMS)) continue;
1130  xt+=TIME-time; xxt+=(TIME-time)*(TIME-time); n++;
1131  }
1132  }
1133  }
1134  }
1135  if(n<10) return false;
1136  *ave=xt/n;
1137  *rms=sqrt(xxt/n-(xt*xt)/(n*n));
1138  return true;
1139 }
1140 
1141 
1142 void HcalDetDiagLaserMonitor::fillHistos(int sd){
1143  if(sd==HcalBarrel || sd==HcalEndcap){
1144  hbheEnergy->Reset();
1145  hbheTime->Reset();
1146  hbheEnergyRMS->Reset();
1147  hbheTimeRMS->Reset();
1148  hb_time_rbx->Reset();
1149  he_time_rbx->Reset();
1150  }
1151  if(sd==HcalOuter){
1152  hoEnergy->Reset();
1153  hoTime->Reset();
1154  hoEnergyRMS->Reset();
1155  hoTimeRMS->Reset();
1156  Time2Dho->Reset();
1157  Energy2Dho->Reset();
1158  refTime2Dho->Reset();
1159  refEnergy2Dho->Reset();
1160  ho_time_rbx->Reset();
1161  }
1162  if(sd==HcalForward){
1163  hfEnergy->Reset();
1164  hfTime->Reset();
1165  hfEnergyRMS->Reset();
1166  hfTimeRMS->Reset();
1167  hf_time_rbx->Reset();
1168  }
1169  if(sd==HcalBarrel || sd==HcalEndcap){
1170  // HB histograms
1171  for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++){
1172  double T=0,nT=0,E=0,nE=0;
1173  for(int depth=1;depth<=2;depth++){
1174  if(hb_data[eta+42][phi-1][depth-1].get_statistics()>10){
1175  double ave=0;
1176  double rms=0;
1177  double time=0;
1178  double time_rms=0;
1179  hb_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
1180  hb_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
1181  hbheEnergy->Fill(ave);
1182  if(ave>0)hbheEnergyRMS->Fill(rms/ave);
1183  hbheTime->Fill(time);
1184  hbheTimeRMS->Fill(time_rms);
1185  T+=time; nT++; E+=ave; nE++;
1186  }
1187  }
1188  if(nT>0){Time2Dhbhehf->setBinContent(eta+44,phi+1,T/nT);Energy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE); }
1189  }
1190  // HE histograms
1191  for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++){
1192  double T=0,nT=0,E=0,nE=0;
1193  for(int depth=1;depth<=3;depth++){
1194  if(he_data[eta+42][phi-1][depth-1].get_statistics()>10){
1195  double ave=0; double rms=0; double time=0; double time_rms=0;
1196  he_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
1197  he_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
1198  hbheEnergy->Fill(ave);
1199  if(ave>0)hbheEnergyRMS->Fill(rms/ave);
1200  hbheTime->Fill(time);
1201  hbheTimeRMS->Fill(time_rms);
1202  T+=time; nT++; E+=ave; nE++;
1203  }
1204  }
1205  if(nT>0 && abs(eta)>16 ){Time2Dhbhehf->setBinContent(eta+44,phi+1,T/nT); Energy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE); }
1206  if(nT>0 && abs(eta)>20 ){Time2Dhbhehf->setBinContent(eta+44,phi+2,T/nT); Energy2Dhbhehf->setBinContent(eta+44,phi+2,E/nE);}
1207  }
1208  }
1209  if(sd==HcalForward){
1210  // HF histograms
1211  for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++){
1212  double T=0,nT=0,E=0,nE=0;
1213  for(int depth=1;depth<=2;depth++){
1214  if(hf_data[eta+42][phi-1][depth-1].get_statistics()>10){
1215  double ave=0; double rms=0; double time=0; double time_rms=0;
1216  hf_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
1217  hf_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
1218  hfEnergy->Fill(ave);
1219  if(ave>0)hfEnergyRMS->Fill(rms/ave);
1220  hfTime->Fill(time);
1221  T+=time; nT++; E+=ave; nE++;
1222  hfTimeRMS->Fill(time_rms);
1223  }
1224  }
1225  if(nT>0 && abs(eta)>29 ){ Time2Dhbhehf->setBinContent(eta+44,phi+1,T/nT); Time2Dhbhehf->setBinContent(eta+44,phi+2,T/nT);}
1226  if(nT>0 && abs(eta)>29 ){ Energy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE); Energy2Dhbhehf->setBinContent(eta+44,phi+2,E/nE);}
1227  }
1228  }
1229  if(sd==HcalOuter){
1230  // HO histograms
1231  for(int eta=-10;eta<=15;eta++) for(int phi=1;phi<=72;phi++){
1232  if(eta>10 && !isSiPM(eta,phi,4)) continue;
1233  double T=0,nT=0,E=0,nE=0;
1234  for(int depth=4;depth<=4;depth++){
1235  if(ho_data[eta+42][phi-1][depth-1].get_statistics()>10){
1236  double ave=0; double rms=0; double time=0; double time_rms=0;
1237  ho_data[eta+42][phi-1][depth-1].get_average_amp(&ave,&rms);
1238  ho_data[eta+42][phi-1][depth-1].get_average_time(&time,&time_rms);
1239  hoEnergy->Fill(ave);
1240  if(ave>0)hoEnergyRMS->Fill(rms/ave);
1241  hoTime->Fill(time);
1242  T+=time; nT++; E+=ave; nE++;
1243  hoTimeRMS->Fill(time_rms);
1244  }
1245  }
1246  if(nT>0){ Time2Dho->Fill(eta,phi,T/nT); Energy2Dho->Fill(eta,phi+1,E/nE) ;}
1247  }
1248  }
1249 
1250  // compare with reference...
1251  if(sd==HcalBarrel || sd==HcalEndcap){
1252  for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++){
1253  double T=0,nT=0,E=0,nE=0;
1254  for(int depth=1;depth<=2;depth++){
1255  if(hb_data[eta+42][phi-1][depth-1].get_statistics()>10){
1256  double val=0,rms=0,time=0,time_rms=0;
1257  double VAL=0,RMS=0,TIME=0,TIME_RMS=0;
1258  if(!hb_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1259  if(!hb_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
1260  if(!hb_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
1261  E+=VAL/val; nE++;
1262  T+=TIME-time; nT++;
1263  }
1264  }
1265  if(nE>0) refEnergy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE);
1266  if(nT>0){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01; refTime2Dhbhehf->setBinContent(eta+44,phi+1,TTT); }
1267  }
1268  for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++){
1269  double T=0,nT=0,E=0,nE=0;
1270  for(int depth=1;depth<=3;depth++){
1271  if(he_data[eta+42][phi-1][depth-1].get_statistics()>10){
1272  double val=0,rms=0,time=0,time_rms=0;
1273  double VAL=0,RMS=0,TIME=0,TIME_RMS=0;
1274  if(!he_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1275  if(!he_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
1276  if(!he_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
1277  E+=VAL/val; nE++;
1278  T+=TIME-time; nT++;
1279  }
1280  }
1281  if(nE>0 && abs(eta)>16) refEnergy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE);
1282  if(nT>0 && abs(eta)>16){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01; refTime2Dhbhehf->setBinContent(eta+44,phi+1,TTT); }
1283  if(nE>0 && abs(eta)>20) refEnergy2Dhbhehf->setBinContent(eta+44,phi+2,E/nE);
1284  if(nT>0 && abs(eta)>20){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01; refTime2Dhbhehf->setBinContent(eta+44,phi+2,TTT); }
1285  }
1286  }
1287  if(sd==HcalForward){
1288  for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++){
1289  double T=0,nT=0,E=0,nE=0;
1290  for(int depth=1;depth<=2;depth++){
1291  if(hf_data[eta+42][phi-1][depth-1].get_statistics()>10){
1292  double val=0,rms=0,time=0,time_rms=0;
1293  double VAL=0,RMS=0,TIME=0,TIME_RMS=0;
1294  if(!hf_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1295  if(!hf_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
1296  if(!hf_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
1297  E+=VAL/val; nE++;
1298  T+=TIME-time; nT++;
1299  }
1300  }
1301  if(nE>0 && abs(eta)>29) refEnergy2Dhbhehf->setBinContent(eta+44,phi+1,E/nE);
1302  if(nT>0 && abs(eta)>29){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01; refTime2Dhbhehf->setBinContent(eta+44,phi+1,TTT); }
1303  if(nE>0 && abs(eta)>29) refEnergy2Dhbhehf->setBinContent(eta+44,phi+2,E/nE);
1304  if(nT>0 && abs(eta)>29){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01; refTime2Dhbhehf->setBinContent(eta+44,phi+2,TTT); }
1305  }
1306  }
1307  if(sd==HcalOuter){
1308  for(int eta=-15;eta<=15;eta++) for(int phi=1;phi<=72;phi++){
1309  if(eta>10 && !isSiPM(eta,phi,4)) continue;
1310  double T=0,nT=0,E=0,nE=0;
1311  for(int depth=4;depth<=4;depth++){
1312  if(ho_data[eta+42][phi-1][depth-1].get_statistics()>10){
1313  double val=0,rms=0,time=0,time_rms=0;
1314  double VAL=0,RMS=0,TIME=0,TIME_RMS=0;
1315  if(!ho_data[eta+42][phi-1][depth-1].get_reference(&val,&rms,&time,&time_rms)) continue;
1316  if(!ho_data[eta+42][phi-1][depth-1].get_average_amp(&VAL,&RMS)) continue;
1317  if(!ho_data[eta+42][phi-1][depth-1].get_average_time(&TIME,&TIME_RMS)) continue;
1318  E+=VAL/val; nE++;
1319  T+=TIME-time; nT++;
1320  }
1321  }
1322  if(nE>0) refEnergy2Dho->Fill(eta,phi,E/nE);
1323  if(nT>0){ double TTT=T/nT+1; if(TTT<0.01) TTT=0.01; refTime2Dho->Fill(eta,phi,TTT);}
1324  }
1325  }
1327  float min,max;
1328  if(sd==HcalBarrel || sd==HcalEndcap){
1329  min=100;max=-100;
1330  for(int i=1;i<=18;i++){
1331  float ave=-10,rms=-10;
1332  if(get_ave_rbx(HcalBarrel,-1,i,&ave,&rms)){
1333  hb_time_rbx->setBinContent(i,ave);
1334  hb_time_rbx->setBinError(i,rms);
1335  if(ave<min) min=ave;
1336  if(ave>max) max=ave;
1337  }
1338  }
1339  for(int i=1;i<=18;i++){
1340  float ave=-10,rms=-10;
1341  if(get_ave_rbx(HcalBarrel,1,i,&ave,&rms)){
1342  hb_time_rbx->setBinContent(i+18,ave);
1343  hb_time_rbx->setBinError(i+18,rms);
1344  if(ave<min) min=ave;
1345  if(ave>max) max=ave;
1346  }
1347  }
1348  if(max>-100)hb_time_rbx->setAxisRange(min-1,max+1,2);
1349  min=100;max=-100;
1350  for(int i=1;i<=18;i++){
1351  float ave=-10,rms=-10;
1352  if(get_ave_rbx(HcalEndcap,-1,i,&ave,&rms)){
1353  he_time_rbx->setBinContent(i,ave);
1354  he_time_rbx->setBinError(i,rms);
1355  if(ave<min) min=ave;
1356  if(ave>max) max=ave;
1357  }
1358  }
1359  for(int i=1;i<=18;i++){
1360  float ave=-10,rms=-10;
1361  if(get_ave_rbx(HcalEndcap,1,i,&ave,&rms)){
1362  he_time_rbx->setBinContent(i+18,ave);
1363  he_time_rbx->setBinError(i+18,rms);
1364  if(ave<min) min=ave;
1365  if(ave>max) max=ave;
1366  }
1367  }
1368  if(max>-100)he_time_rbx->setAxisRange(min-1,max+1,2);
1369  }
1371  if(sd==HcalOuter){
1372  min=100;max=-100;
1373  for(int i=1;i<=6;i++){
1374  float ave=-10,rms=-10;
1375  if(get_ave_rbx(HcalOuter,-2,i,&ave,&rms)){
1376  ho_time_rbx->setBinContent(i,ave);
1377  ho_time_rbx->setBinError(i,rms);
1378  if(ave<min) min=ave;
1379  if(ave>max) max=ave;
1380  }
1381  }
1382  for(int i=1;i<=6;i++){
1383  float ave=-10,rms=-10;
1384  if(get_ave_rbx(HcalOuter,-1,i,&ave,&rms)){
1385  ho_time_rbx->setBinContent(i+6,ave);
1386  ho_time_rbx->setBinError(i+6,rms);
1387  if(ave<min) min=ave;
1388  if(ave>max) max=ave;
1389  }
1390  }
1391  for(int i=1;i<=12;i++){
1392  float ave=-10,rms=-10;
1393  if(get_ave_rbx(HcalOuter,0,i,&ave,&rms)){
1394  ho_time_rbx->setBinContent(i+12,ave);
1395  ho_time_rbx->setBinError(i+12,rms);
1396  if(ave<min) min=ave;
1397  if(ave>max) max=ave;
1398  }
1399  }
1400  for(int i=1;i<=6;i++){
1401  float ave=-10,rms=-10;
1402  if(get_ave_rbx(HcalOuter,1,i,&ave,&rms)){
1403  ho_time_rbx->setBinContent(i+24,ave);
1404  ho_time_rbx->setBinError(i+24,rms);
1405  if(ave<min) min=ave;
1406  if(ave>max) max=ave;
1407  }
1408  }
1409  for(int i=1;i<=6;i++){
1410  float ave=-10,rms=-10;
1411  if(get_ave_rbx(HcalOuter,2,i,&ave,&rms)){
1412  ho_time_rbx->setBinContent(i+30,ave);
1413  ho_time_rbx->setBinError(i+30,rms);
1414  if(ave<min) min=ave;
1415  if(ave>max) max=ave;
1416  }
1417  }
1418  if(max>-100)ho_time_rbx->setAxisRange(min-1,max+1,2);
1419  }
1421  if(sd==HcalForward){
1422  min=100;max=-100;
1423  for(int i=1;i<=12;i++){
1424  float ave=-10,rms=-10;
1425  if(get_ave_rbx(HcalForward,-1,i,&ave,&rms)){
1426  hf_time_rbx->setBinContent(i,ave);
1427  hf_time_rbx->setBinError(i,rms);
1428  if(ave<min) min=ave;
1429  if(ave>max) max=ave;
1430  }
1431  }
1432  for(int i=1;i<=12;i++){
1433  float ave=-10,rms=-10;
1434  if(get_ave_rbx(HcalForward,1,i,&ave,&rms)){
1435  hf_time_rbx->setBinContent(i+12,ave);
1436  hf_time_rbx->setBinError(i+12,rms);
1437  if(ave<min) min=ave;
1438  if(ave>max) max=ave;
1439  }
1440  }
1441  if(max>-100)hf_time_rbx->setAxisRange(min-1,max+1,2);
1442  }
1443 }
1444 
1445 void HcalDetDiagLaserMonitor::SaveReference(){
1446 double amp,rms,Time,time_rms;
1447 int Eta,Phi,Depth,Statistic,Status=0;
1448 char Subdet[10],str[500];
1449  if(OutputFilePath.size()>0){
1450  if(!Overwrite){
1451  sprintf(str,"%sHcalDetDiagLaserData_run%06i_%i.root",OutputFilePath.c_str(),run_number,dataset_seq_number);
1452  }else{
1453  sprintf(str,"%sHcalDetDiagLaserData.root",OutputFilePath.c_str());
1454  }
1455  TFile *theFile = new TFile(str, "RECREATE");
1456  if(!theFile->IsOpen()) return;
1457  theFile->cd();
1458  sprintf(str,"%d",run_number); TObjString run(str); run.Write("run number");
1459  sprintf(str,"%d",ievt_); TObjString events(str); events.Write("Total events processed");
1460  sprintf(str,"%d",dataset_seq_number); TObjString dsnum(str); dsnum.Write("Dataset number");
1461  Long_t t; t=time(0); strftime(str,30,"%F %T",localtime(&t)); TObjString tm(str); tm.Write("Dataset creation time");
1462 
1463  TTree *tree =new TTree("HCAL Laser data","HCAL Laser data");
1464  if(tree==0) return;
1465  tree->Branch("Subdet", &Subdet, "Subdet/C");
1466  tree->Branch("eta", &Eta, "Eta/I");
1467  tree->Branch("phi", &Phi, "Phi/I");
1468  tree->Branch("depth", &Depth, "Depth/I");
1469  tree->Branch("statistic",&Statistic, "Statistic/I");
1470  tree->Branch("status", &Status, "Status/I");
1471  tree->Branch("amp", &amp, "amp/D");
1472  tree->Branch("rms", &rms, "rms/D");
1473  tree->Branch("time", &Time, "time/D");
1474  tree->Branch("time_rms", &time_rms, "time_rms/D");
1475  sprintf(Subdet,"HB");
1476  for(int eta=-16;eta<=16;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=2;depth++){
1477  if((Statistic=hb_data[eta+42][phi-1][depth-1].get_statistics1())>10){
1478  Eta=eta; Phi=phi; Depth=depth;
1479  Status=hb_data[eta+42][phi-1][depth-1].get_status();
1480  hb_data[eta+42][phi-1][depth-1].get_average_amp1(&amp,&rms);
1481  hb_data[eta+42][phi-1][depth-1].get_average_time1(&Time,&time_rms);
1482  tree->Fill();
1483  }
1484  }
1485  sprintf(Subdet,"HE");
1486  for(int eta=-29;eta<=29;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=3;depth++){
1487  if((Statistic=he_data[eta+42][phi-1][depth-1].get_statistics1())>10){
1488  Eta=eta; Phi=phi; Depth=depth;
1489  Status=he_data[eta+42][phi-1][depth-1].get_status();
1490  he_data[eta+42][phi-1][depth-1].get_average_amp1(&amp,&rms);
1491  he_data[eta+42][phi-1][depth-1].get_average_time1(&Time,&time_rms);
1492  tree->Fill();
1493  }
1494  }
1495  sprintf(Subdet,"HO");
1496  for(int eta=-15;eta<=15;eta++) for(int phi=1;phi<=72;phi++) for(int depth=4;depth<=4;depth++){
1497  if((Statistic=ho_data[eta+42][phi-1][depth-1].get_statistics1())>10){
1498  Eta=eta; Phi=phi; Depth=depth;
1499  Status=ho_data[eta+42][phi-1][depth-1].get_status();
1500  ho_data[eta+42][phi-1][depth-1].get_average_amp1(&amp,&rms);
1501  ho_data[eta+42][phi-1][depth-1].get_average_time1(&Time,&time_rms);
1502  tree->Fill();
1503  }
1504  }
1505  sprintf(Subdet,"HF");
1506  for(int eta=-42;eta<=42;eta++) for(int phi=1;phi<=72;phi++) for(int depth=1;depth<=2;depth++){
1507  if((Statistic=hf_data[eta+42][phi-1][depth-1].get_statistics1())>10){
1508  Eta=eta; Phi=phi; Depth=depth;
1509  Status=hf_data[eta+42][phi-1][depth-1].get_status();
1510  hf_data[eta+42][phi-1][depth-1].get_average_amp1(&amp,&rms);
1511  hf_data[eta+42][phi-1][depth-1].get_average_time1(&Time,&time_rms);
1512  tree->Fill();
1513  }
1514  }
1515  sprintf(Subdet,"CALIB_HB");
1516  for(int eta=-1;eta<=1;eta++) for(int phi=1;phi<=72;phi++){
1517  if((Statistic=calib_data[1][eta+2][phi-1].get_statistics1())>10){
1518  Eta=eta; Phi=phi; Depth=0;
1519  Status=calib_data[1][eta+2][phi-1].get_status();
1520  calib_data[1][eta+2][phi-1].get_average_amp1(&amp,&rms);
1521  calib_data[1][eta+2][phi-1].get_average_time1(&Time,&time_rms);
1522  tree->Fill();
1523  }
1524  }
1525  sprintf(Subdet,"CALIB_HE");
1526  for(int eta=-1;eta<=1;eta++) for(int phi=1;phi<=72;phi++){
1527  if((Statistic=calib_data[2][eta+2][phi-1].get_statistics1())>10){
1528  Eta=eta; Phi=phi; Depth=0;
1529  Status=calib_data[2][eta+2][phi-1].get_status();
1530  calib_data[2][eta+2][phi-1].get_average_amp1(&amp,&rms);
1531  calib_data[2][eta+2][phi-1].get_average_time1(&Time,&time_rms);
1532  tree->Fill();
1533  }
1534  }
1535  sprintf(Subdet,"CALIB_HO");
1536  for(int eta=-2;eta<=2;eta++) for(int phi=1;phi<=72;phi++){
1537  if((Statistic=calib_data[3][eta+2][phi-1].get_statistics1())>10){
1538  Eta=eta; Phi=phi; Depth=0;
1539  Status=calib_data[3][eta+2][phi-1].get_status();
1540  calib_data[3][eta+2][phi-1].get_average_amp1(&amp,&rms);
1541  calib_data[3][eta+2][phi-1].get_average_time1(&Time,&time_rms);
1542  tree->Fill();
1543  }
1544  }
1545  sprintf(Subdet,"CALIB_HF");
1546  for(int eta=-2;eta<=2;eta++) for(int phi=1;phi<=72;phi++){
1547  if((Statistic=calib_data[4][eta+2][phi-1].get_statistics1())>10){
1548  Eta=eta; Phi=phi; Depth=0;
1549  Status=calib_data[4][eta+2][phi-1].get_status();
1550  calib_data[4][eta+2][phi-1].get_average_amp1(&amp,&rms);
1551  calib_data[4][eta+2][phi-1].get_average_time1(&Time,&time_rms);
1552  tree->Fill();
1553  }
1554  }
1555  theFile->Write();
1556  theFile->Close();
1557  theFile->Delete();
1558  }
1559  if(XmlFilePath.size()>0){
1560  char TIME[40];
1561  Long_t t; t=time(0); strftime(TIME,30,"%F %T",localtime(&t));
1562  //create XML file
1563  if(!Overwrite){
1564  sprintf(str,"HcalDetDiagLaser_%i_%i.xml",run_number,dataset_seq_number);
1565  }else{
1566  sprintf(str,"HcalDetDiagLaser.xml");
1567  }
1568  std::string xmlName=str;
1569  std::ofstream xmlFile;
1570  xmlFile.open(xmlName.c_str());
1571 
1572  xmlFile<<"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"yes\"?>\n";
1573  xmlFile<<"<ROOT>\n";
1574  xmlFile<<" <HEADER>\n";
1575  xmlFile<<" <HINTS mode='only-det-root'/>\n";
1576  xmlFile<<" <TYPE>\n";
1577  xmlFile<<" <EXTENSION_TABLE_NAME>HCAL_DETMON_LED_LASER_V1</EXTENSION_TABLE_NAME>\n";
1578  xmlFile<<" <NAME>HCAL Laser HBHE HPD [abort gap global]</NAME>\n";
1579  xmlFile<<" </TYPE>\n";
1580  xmlFile<<" <!-- run details -->\n";
1581  xmlFile<<" <RUN>\n";
1582  xmlFile<<" <RUN_TYPE>GLOBAL-RUN</RUN_TYPE>\n";
1583  xmlFile<<" <RUN_NUMBER>"<<run_number<<"</RUN_NUMBER>\n";
1584  xmlFile<<" <RUN_BEGIN_TIMESTAMP>2009-01-01 00:00:00</RUN_BEGIN_TIMESTAMP>\n";
1585  xmlFile<<" <COMMENT_DESCRIPTION>hcal laser data</COMMENT_DESCRIPTION>\n";
1586  xmlFile<<" <LOCATION>P5</LOCATION>\n";
1587  xmlFile<<" <INITIATED_BY_USER>dma</INITIATED_BY_USER>\n";
1588  xmlFile<<" </RUN>\n";
1589  xmlFile<<" </HEADER>\n";
1590  xmlFile<<" <DATA_SET>\n";
1591  xmlFile<<" <!-- optional dataset metadata -->\n\n";
1592  xmlFile<<" <SET_NUMBER>"<<dataset_seq_number<<"</SET_NUMBER>\n";
1593  xmlFile<<" <SET_BEGIN_TIMESTAMP>2009-01-01 00:00:00</SET_BEGIN_TIMESTAMP>\n";
1594  xmlFile<<" <SET_END_TIMESTAMP>2009-01-01 00:00:00</SET_END_TIMESTAMP>\n";
1595  xmlFile<<" <NUMBER_OF_EVENTS_IN_SET>"<<ievt_<<"</NUMBER_OF_EVENTS_IN_SET>\n";
1596  xmlFile<<" <COMMENT_DESCRIPTION>Automatic DQM output</COMMENT_DESCRIPTION>\n";
1597  xmlFile<<" <DATA_FILE_NAME>"<< xmlName <<"</DATA_FILE_NAME>\n";
1598  xmlFile<<" <IMAGE_FILE_NAME>data plot url or file path</IMAGE_FILE_NAME>\n";
1599  xmlFile<<" <!-- who and when created this dataset-->\n\n";
1600  xmlFile<<" <CREATE_TIMESTAMP>"<<TIME<<"</CREATE_TIMESTAMP>\n";
1601  xmlFile<<" <CREATED_BY_USER>dma</CREATED_BY_USER>\n";
1602  xmlFile<<" <!-- version (string) and subversion (number) -->\n";
1603  xmlFile<<" <!-- fields are used to read data back from the database -->\n\n";
1604  xmlFile<<" <VERSION>"<<run_number<<dataset_seq_number<<"</VERSION>\n";
1605  xmlFile<<" <SUBVERSION>1</SUBVERSION>\n";
1606  xmlFile<<" <!-- Assign predefined dataset attributes -->\n\n";
1607  xmlFile<<" <PREDEFINED_ATTRIBUTES>\n";
1608  xmlFile<<" <ATTRIBUTE>\n";
1609  xmlFile<<" <NAME>HCAL Dataset Status</NAME>\n";
1610  xmlFile<<" <VALUE>VALID</VALUE>\n";
1611  xmlFile<<" </ATTRIBUTE>\n";
1612  xmlFile<<" </PREDEFINED_ATTRIBUTES>\n";
1613  xmlFile<<" <!-- multiple data block records -->\n\n";
1614 
1615  std::vector <HcalElectronicsId> AllElIds = emap->allElectronicsIdPrecision();
1616  for(std::vector <HcalElectronicsId>::iterator eid = AllElIds.begin(); eid != AllElIds.end(); eid++){
1617  DetId detid=emap->lookup(*eid);
1618  if (detid.det()!=DetId::Hcal) continue;
1619  HcalGenericDetId gid(emap->lookup(*eid));
1620  if(!(!(gid.null()) &&
1621  (gid.genericSubdet()==HcalGenericDetId::HcalGenBarrel ||
1622  gid.genericSubdet()==HcalGenericDetId::HcalGenEndcap ||
1623  gid.genericSubdet()==HcalGenericDetId::HcalGenForward ||
1624  gid.genericSubdet()==HcalGenericDetId::HcalGenOuter))) continue;
1625  int eta,phi,depth;
1626  std::string subdet="";
1627  HcalDetId hid(detid);
1628  eta=hid.ieta();
1629  phi=hid.iphi();
1630  depth=hid.depth();
1631 
1632  double e=0,e_rms=0,t=0,t_rms=0;
1633  if(detid.subdetId()==HcalBarrel){
1634  subdet="HB";
1635  Statistic=hb_data[eta+42][phi-1][depth-1].get_statistics1();
1636  Status =hb_data[eta+42][phi-1][depth-1].get_status();
1637  hb_data[eta+42][phi-1][depth-1].get_average_amp1(&e,&e_rms);
1638  hb_data[eta+42][phi-1][depth-1].get_average_time1(&t,&t_rms);
1639  }else if(detid.subdetId()==HcalEndcap){
1640  subdet="HE";
1641  Statistic=he_data[eta+42][phi-1][depth-1].get_statistics1();
1642  Status =he_data[eta+42][phi-1][depth-1].get_status();
1643  he_data[eta+42][phi-1][depth-1].get_average_amp1(&e,&e_rms);
1644  he_data[eta+42][phi-1][depth-1].get_average_time1(&t,&t_rms);
1645  }else if(detid.subdetId()==HcalForward){
1646  subdet="HF";
1647  Statistic=hf_data[eta+42][phi-1][depth-1].get_statistics1();
1648  Status =hf_data[eta+42][phi-1][depth-1].get_status();
1649  hf_data[eta+42][phi-1][depth-1].get_average_amp1(&e,&e_rms);
1650  hf_data[eta+42][phi-1][depth-1].get_average_time1(&t,&t_rms);
1651  }else if(detid.subdetId()==HcalOuter){
1652  subdet="HO";
1653  Statistic=ho_data[eta+42][phi-1][depth-1].get_statistics1();
1654  Status =ho_data[eta+42][phi-1][depth-1].get_status();
1655  ho_data[eta+42][phi-1][depth-1].get_average_amp1(&e,&e_rms);
1656  ho_data[eta+42][phi-1][depth-1].get_average_time1(&t,&t_rms);
1657  }else continue;
1658  xmlFile<<" <DATA>\n";
1659  xmlFile<<" <NUMBER_OF_EVENTS_USED>"<<Statistic<<"</NUMBER_OF_EVENTS_USED>\n";
1660  xmlFile<<" <SIGNAL_MEAN>"<<e<<"</SIGNAL_MEAN>\n";
1661  xmlFile<<" <SIGNAL_RMS>"<<e_rms<<"</SIGNAL_RMS>\n";
1662  xmlFile<<" <TIME_MEAN>"<<t<<"</TIME_MEAN>\n";
1663  xmlFile<<" <TIME_RMS>"<<t_rms<<"</TIME_RMS>\n";
1664  xmlFile<<" <CHANNEL_STATUS_WORD>"<<Status<<"</CHANNEL_STATUS_WORD>\n";
1665  xmlFile<<" <CHANNEL_OBJECTNAME>HcalDetId</CHANNEL_OBJECTNAME>\n";
1666  xmlFile<<" <SUBDET>"<<subdet<<"</SUBDET>\n";
1667  xmlFile<<" <IETA>"<<eta<<"</IETA>\n";
1668  xmlFile<<" <IPHI>"<<phi<<"</IPHI>\n";
1669  xmlFile<<" <DEPTH>"<<depth<<"</DEPTH>\n";
1670  xmlFile<<" <TYPE>0</TYPE>\n";
1671  xmlFile<<" </DATA>\n";
1672  }
1674  xmlFile<<" </DATA_SET>\n";
1675  xmlFile<<"</ROOT>\n";
1676  xmlFile.close();
1677 
1678 
1679  //create CALIB XML file
1680  sprintf(str,"HcalDetDiagLaserCalib_%i_%i.xml",run_number,dataset_seq_number);
1681  std::string xmlNameCalib=str;
1682  std::ofstream xmlFileCalib;
1683  xmlFileCalib.open(xmlNameCalib.c_str());
1684 
1685  xmlFileCalib<<"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"yes\"?>\n";
1686  xmlFileCalib<<"<ROOT>\n";
1687  xmlFileCalib<<" <HEADER>\n";
1688  xmlFileCalib<<" <HINTS mode='only-det-root'/>\n";
1689  xmlFileCalib<<" <TYPE>\n";
1690  xmlFileCalib<<" <EXTENSION_TABLE_NAME>HCAL_DETMON_LED_LASER_V1</EXTENSION_TABLE_NAME>\n";
1691  xmlFileCalib<<" <NAME>HCAL Laser CALIB [abort gap global]</NAME>\n";
1692  xmlFileCalib<<" </TYPE>\n";
1693  xmlFileCalib<<" <!-- run details -->\n";
1694  xmlFileCalib<<" <RUN>\n";
1695  xmlFileCalib<<" <RUN_TYPE>Global-RUN</RUN_TYPE>\n";
1696  xmlFileCalib<<" <RUN_NUMBER>"<<run_number<<"</RUN_NUMBER>\n";
1697  xmlFileCalib<<" <RUN_BEGIN_TIMESTAMP>2009-01-01 00:00:00</RUN_BEGIN_TIMESTAMP>\n";
1698  xmlFileCalib<<" <COMMENT_DESCRIPTION>hcal Laser CALIB data</COMMENT_DESCRIPTION>\n";
1699  xmlFileCalib<<" <LOCATION>P5</LOCATION>\n";
1700  xmlFileCalib<<" <INITIATED_BY_USER>dma</INITIATED_BY_USER>\n";
1701  xmlFileCalib<<" </RUN>\n";
1702  xmlFileCalib<<" </HEADER>\n";
1703  xmlFileCalib<<" <DATA_SET>\n";
1704  xmlFileCalib<<" <!-- optional dataset metadata -->\n\n";
1705  xmlFileCalib<<" <SET_NUMBER>"<<dataset_seq_number<<"</SET_NUMBER>\n";
1706  xmlFileCalib<<" <SET_BEGIN_TIMESTAMP>2009-01-01 00:00:00</SET_BEGIN_TIMESTAMP>\n";
1707  xmlFileCalib<<" <SET_END_TIMESTAMP>2009-01-01 00:00:00</SET_END_TIMESTAMP>\n";
1708  xmlFileCalib<<" <NUMBER_OF_EVENTS_IN_SET>"<<ievt_<<"</NUMBER_OF_EVENTS_IN_SET>\n";
1709  xmlFileCalib<<" <COMMENT_DESCRIPTION>Automatic DQM output</COMMENT_DESCRIPTION>\n";
1710  xmlFileCalib<<" <DATA_FILE_NAME>"<< xmlNameCalib <<"</DATA_FILE_NAME>\n";
1711  xmlFileCalib<<" <IMAGE_FILE_NAME>data plot url or file path</IMAGE_FILE_NAME>\n";
1712  xmlFileCalib<<" <!-- who and when created this dataset-->\n\n";
1713  xmlFileCalib<<" <CREATE_TIMESTAMP>"<<TIME<<"</CREATE_TIMESTAMP>\n";
1714  xmlFileCalib<<" <CREATED_BY_USER>dma</CREATED_BY_USER>\n";
1715  xmlFileCalib<<" <!-- version (string) and subversion (number) -->\n";
1716  xmlFileCalib<<" <!-- fields are used to read data back from the database -->\n\n";
1717  xmlFileCalib<<" <VERSION>"<<run_number<<dataset_seq_number<<"</VERSION>\n";
1718  xmlFileCalib<<" <SUBVERSION>1</SUBVERSION>\n";
1719  xmlFileCalib<<" <!-- Assign predefined dataset attributes -->\n\n";
1720  xmlFileCalib<<" <PREDEFINED_ATTRIBUTES>\n";
1721  xmlFileCalib<<" <ATTRIBUTE>\n";
1722  xmlFileCalib<<" <NAME>HCAL Dataset Status</NAME>\n";
1723  xmlFileCalib<<" <VALUE>VALID</VALUE>\n";
1724  xmlFileCalib<<" </ATTRIBUTE>\n";
1725  xmlFileCalib<<" </PREDEFINED_ATTRIBUTES>\n";
1726  xmlFileCalib<<" <!-- multiple data block records -->\n\n";
1727 
1728  for(int sd=1;sd<=4;sd++) for(int eta=-2;eta<=2;eta++) for(int phi=1;phi<=72;phi++){
1729  std::string subdet="";
1730  if(sd==1) subdet="HB";
1731  if(sd==2) subdet="HE";
1732  if(sd==3) subdet="HO";
1733  if(sd==4) subdet="HF";
1734  if((calib_data[sd][eta+2][phi-1].get_statistics())>100){
1735  double e=0,e_rms=0,t=0,t_rms=0;
1736  Status=calib_data[sd][eta+2][phi-1].get_status();
1737  Statistic=calib_data[sd][eta+2][phi-1].get_statistics1();
1738  calib_data[sd][eta+2][phi-1].get_average_amp1(&e,&e_rms);
1739  calib_data[sd][eta+2][phi-1].get_average_time1(&t,&t_rms);
1740  xmlFileCalib<<" <DATA>\n";
1741  xmlFileCalib<<" <NUMBER_OF_EVENTS_USED>"<<Statistic<<"</NUMBER_OF_EVENTS_USED>\n";
1742  xmlFileCalib<<" <SIGNAL_MEAN>"<<e<<"</SIGNAL_MEAN>\n";
1743  xmlFileCalib<<" <SIGNAL_RMS>"<<e_rms<<"</SIGNAL_RMS>\n";
1744  xmlFileCalib<<" <TIME_MEAN>"<<t<<"</TIME_MEAN>\n";
1745  xmlFileCalib<<" <TIME_RMS>"<<t_rms<<"</TIME_RMS>\n";
1746  xmlFileCalib<<" <CHANNEL_STATUS_WORD>"<<Status<<"</CHANNEL_STATUS_WORD>\n";
1747  xmlFileCalib<<" <CHANNEL_OBJECTNAME>HcalDetId</CHANNEL_OBJECTNAME>\n";
1748  xmlFileCalib<<" <SUBDET>"<<subdet<<"</SUBDET>\n";
1749  xmlFileCalib<<" <IETA>"<<eta<<"</IETA>\n";
1750  xmlFileCalib<<" <IPHI>"<<phi<<"</IPHI>\n";
1751  xmlFileCalib<<" <DEPTH>"<<0<<"</DEPTH>\n";
1752  xmlFileCalib<<" <TYPE>0</TYPE>\n";
1753  xmlFileCalib<<" </DATA>\n";
1754  }
1755  }
1757  xmlFileCalib<<" </DATA_SET>\n";
1758  xmlFileCalib<<"</ROOT>\n";
1759  xmlFileCalib.close();
1760 
1761  sprintf(str,"zip %s.zip %s %s",xmlName.c_str(),xmlName.c_str(),xmlNameCalib.c_str());
1762  system(str);
1763  sprintf(str,"rm -f %s %s",xmlName.c_str(),xmlNameCalib.c_str());
1764  system(str);
1765  sprintf(str,"mv -f %s.zip %s",xmlName.c_str(),XmlFilePath.c_str());
1766  system(str);
1767 
1768  }
1769  for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) hb_data[i][j][k].reset1();
1770  for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) he_data[i][j][k].reset1();
1771  for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) ho_data[i][j][k].reset1();
1772  for(int i=0;i<85;i++)for(int j=0;j<72;j++)for(int k=0;k<4;k++) hf_data[i][j][k].reset1();
1773  for(int i=0;i<5;i++)for(int j=0;j<5;j++)for(int k=0;k<72;k++)calib_data[i][j][k].reset1();
1774  ievt_=0;
1775  dataset_seq_number++;
1776 }
1777 void HcalDetDiagLaserMonitor::LoadReference(){
1778 double amp,rms,time,time_rms;
1779 int Eta,Phi,Depth;
1780 char subdet[10];
1781 TFile *f;
1782  if(gSystem->AccessPathName(ReferenceData.c_str())) return;
1783  f = new TFile(ReferenceData.c_str(),"READ");
1784 
1785  if(!f->IsOpen()) return ;
1786  TObjString *STR=(TObjString *)f->Get("run number");
1787 
1788  if(STR){ std::string Ref(STR->String()); ReferenceRun=Ref;}
1789 
1790  TTree* t=(TTree*)f->Get("HCAL Laser data");
1791  if(!t) return;
1792  t->SetBranchAddress("Subdet", subdet);
1793  t->SetBranchAddress("eta", &Eta);
1794  t->SetBranchAddress("phi", &Phi);
1795  t->SetBranchAddress("depth", &Depth);
1796  t->SetBranchAddress("amp", &amp);
1797  t->SetBranchAddress("rms", &rms);
1798  t->SetBranchAddress("time", &time);
1799  t->SetBranchAddress("time_rms", &time_rms);
1800  for(int ievt=0;ievt<t->GetEntries();ievt++){
1801  t->GetEntry(ievt);
1802  if(strcmp(subdet,"HB")==0) hb_data[Eta+42][Phi-1][Depth-1].set_reference(amp,rms,time,time_rms);
1803  if(strcmp(subdet,"HE")==0) he_data[Eta+42][Phi-1][Depth-1].set_reference(amp,rms,time,time_rms);
1804  if(strcmp(subdet,"HO")==0) ho_data[Eta+42][Phi-1][Depth-1].set_reference(amp,rms,time,time_rms);
1805  if(strcmp(subdet,"HF")==0) hf_data[Eta+42][Phi-1][Depth-1].set_reference(amp,rms,time,time_rms);
1806  if(strcmp(subdet,"CALIB_HB")==0) calib_data[1][Eta+2][Phi-1].set_reference(amp,rms,time,time_rms);
1807  if(strcmp(subdet,"CALIB_HE")==0) calib_data[2][Eta+2][Phi-1].set_reference(amp,rms,time,time_rms);
1808  if(strcmp(subdet,"CALIB_HO")==0) calib_data[3][Eta+2][Phi-1].set_reference(amp,rms,time,time_rms);
1809  if(strcmp(subdet,"CALIB_HF")==0) calib_data[4][Eta+2][Phi-1].set_reference(amp,rms,time,time_rms);
1810  }
1811  f->Close();
1812  f->Delete();
1813  IsReference=true;
1814 }
1815 void HcalDetDiagLaserMonitor::LoadDataset(){
1816 double amp,rms,time,time_rms;
1817 int Eta,Phi,Depth,Statistic;
1818 char subdet[10];
1819 TFile *f;
1820  if(DatasetName.size()==0) return;
1821  createHTMLonly=true;
1822  if(gSystem->AccessPathName(DatasetName.c_str())) return;
1823  f = new TFile(DatasetName.c_str(),"READ");
1824 
1825  if(!f->IsOpen()) return ;
1826 
1827  TTree* t=0;
1828  t=(TTree*)f->Get("HCAL Laser data");
1829  if(!t) return;
1830  t->SetBranchAddress("Subdet", subdet);
1831  t->SetBranchAddress("eta", &Eta);
1832  t->SetBranchAddress("phi", &Phi);
1833  t->SetBranchAddress("depth", &Depth);
1834  t->SetBranchAddress("amp", &amp);
1835  t->SetBranchAddress("rms", &rms);
1836  t->SetBranchAddress("time", &time);
1837  t->SetBranchAddress("time_rms", &time_rms);
1838  t->SetBranchAddress("statistic",&Statistic);
1839  for(int ievt=0;ievt<t->GetEntries();ievt++){
1840  t->GetEntry(ievt);
1841  if(strcmp(subdet,"HB")==0){ nHB++;
1842  hb_data[Eta+42][Phi-1][Depth-1].set_data(amp,rms,time,time_rms);
1843  hb_data[Eta+42][Phi-1][Depth-1].set_statistics(Statistic);
1844  hb_data[Eta+42][Phi-1][Depth-1].set_statistics1(Statistic);
1845  }
1846  if(strcmp(subdet,"HE")==0){ nHE++;
1847  he_data[Eta+42][Phi-1][Depth-1].set_data(amp,rms,time,time_rms);
1848  he_data[Eta+42][Phi-1][Depth-1].set_statistics(Statistic);
1849  he_data[Eta+42][Phi-1][Depth-1].set_statistics1(Statistic);
1850  }
1851  if(strcmp(subdet,"HO")==0){ nHO++;
1852  ho_data[Eta+42][Phi-1][Depth-1].set_data(amp,rms,time,time_rms);
1853  ho_data[Eta+42][Phi-1][Depth-1].set_statistics(Statistic);
1854  ho_data[Eta+42][Phi-1][Depth-1].set_statistics1(Statistic);
1855  }
1856  if(strcmp(subdet,"HF")==0){ nHF++;
1857  hf_data[Eta+42][Phi-1][Depth-1].set_data(amp,rms,time,time_rms);
1858  hf_data[Eta+42][Phi-1][Depth-1].set_statistics(Statistic);
1859  hf_data[Eta+42][Phi-1][Depth-1].set_statistics1(Statistic);
1860  }
1861  if(strcmp(subdet,"CALIB_HB")==0){
1862  calib_data[1][Eta+2][Phi-1].set_data(amp,rms,time,time_rms);
1863  calib_data[1][Eta+2][Phi-1].set_statistics(Statistic);
1864  calib_data[1][Eta+2][Phi-1].set_statistics1(Statistic);
1865  }
1866  if(strcmp(subdet,"CALIB_HE")==0){
1867  calib_data[2][Eta+2][Phi-1].set_data(amp,rms,time,time_rms);
1868  calib_data[2][Eta+2][Phi-1].set_statistics(Statistic);
1869  calib_data[2][Eta+2][Phi-1].set_statistics1(Statistic);
1870  }
1871  if(strcmp(subdet,"CALIB_HO")==0){
1872  calib_data[3][Eta+2][Phi-1].set_data(amp,rms,time,time_rms);
1873  calib_data[3][Eta+2][Phi-1].set_statistics(Statistic);
1874  calib_data[3][Eta+2][Phi-1].set_statistics1(Statistic);
1875  }
1876  if(strcmp(subdet,"CALIB_HF")==0){
1877  calib_data[4][Eta+2][Phi-1].set_data(amp,rms,time,time_rms);
1878  calib_data[4][Eta+2][Phi-1].set_statistics(Statistic);
1879  calib_data[4][Eta+2][Phi-1].set_statistics1(Statistic);
1880  }
1881  }
1882  TObjString *STR1=(TObjString *)f->Get("run number");
1883  if(STR1){ int run; sscanf(STR1->String(),"%i",&run); meRUN_->Fill(run); run_number=run;}
1884 
1885  TObjString *STR2=(TObjString *)f->Get("Total events processed");
1886  if(STR2){ int events; sscanf(STR2->String(),"%i",&events); meEVT_->Fill(events); ievt_=events;}
1887 
1888  TObjString *STR3=(TObjString *)f->Get("Dataset number");
1889  if(STR3){ int ds; sscanf(STR3->String(),"%i",&ds); dataset_seq_number=ds;}
1890  f->Close();
1891  f->Delete();
1892 }
1893 
1894 void HcalDetDiagLaserMonitor::SaveRaddamData(){
1895 float adc_range[20]={14,28,40,52,67,132,202,262,322,397,722,1072,1372,1672,2047,3672,5422,6922,8422,10297};
1896 int adc_bins[20]={1,2,3,4,5,5,10,15,20,25,25,50,75,100,125,125,250,375,500,625};
1897 char str[100];
1898  TF1 *fitFunc = new TF1("fitFunc","gaus");
1899  if(fitFunc==0) return;
1900  for(int i=0;i<56;i++){
1901  float sum1=0,sum2=0,n=0;
1902  S1[i]=S2[i]=0;
1903  for(int j=0;j<128;j++){
1904  sum1+=(adc2fC[j]-2.5)*Raddam_data[i].s1_adc[j];
1905  sum2+=(adc2fC[j]-2.5)*Raddam_data[i].s2_adc[j];
1906  n+=Raddam_data[i].s1_adc[j];
1907  }
1908  if(n<100) continue;
1909  sum1=sum1/n;
1910  sum2=sum2/n;
1911  int N=0;
1912  int Ws1=1,Ws2=1;
1913  for(N=1;N<19;N++) if(sum1>adc_range[N-1] && sum1<adc_range[N]) break;
1914  Ws1=adc_bins[N+1];
1915  for(N=1;N<19;N++) if(sum2>adc_range[N-1] && sum2<adc_range[N]) break;
1916  Ws2=adc_bins[N+1];
1917  sprintf(str,"Raddam(%i,%i,%i) S1",RADDAM_CH[i].eta,RADDAM_CH[i].phi,RADDAM_CH[i].depth);
1918  S1[i]=new TH1F(str,str,10000/Ws1,0,10000);
1919  sprintf(str,"Raddam(%i,%i,%i) S2",RADDAM_CH[i].eta,RADDAM_CH[i].phi,RADDAM_CH[i].depth);
1920  S2[i]=new TH1F(str,str,10000/Ws1,0,10000);
1921  for(int j=0;j<128;j++){
1922  S1[i]->Fill(adc2fC[j]-2.5,Raddam_data[i].s1_adc[j]);
1923  S2[i]->Fill(adc2fC[j]-2.5,Raddam_data[i].s2_adc[j]);
1924  }
1925  double parm[3];
1926  S1[i]->Fit(fitFunc);
1927  S1[i]->GetFunction("fitFunc")->GetParameters(parm);
1928  Raddam_data[i].S1MEAN=S1[i]->GetMean();
1929  Raddam_data[i].S1RMS=S1[i]->GetRMS();
1930  Raddam_data[i].S1FITMEAN=parm[1];
1931  Raddam_data[i].S1FITMEANER=S1[i]->GetFunction("fitFunc")->GetParError(1);
1932  Raddam_data[i].S1FITSIGMA=parm[2];
1933  Raddam_data[i].S1CHI2=S1[i]->GetFunction("fitFunc")->GetChisquare();
1934  Raddam_data[i].S1NDF=S1[i]->GetFunction("fitFunc")->GetNDF();
1935  Raddam_data[i].S1BINWIDTH=Ws1;
1936  S2[i]->Fit(fitFunc);
1937  S2[i]->GetFunction("fitFunc")->GetParameters(parm);
1938  Raddam_data[i].S2MEAN=S2[i]->GetMean();
1939  Raddam_data[i].S2RMS=S2[i]->GetRMS();
1940  Raddam_data[i].S2FITMEAN=parm[1];
1941  Raddam_data[i].S2FITMEANER=S2[i]->GetFunction("fitFunc")->GetParError(1);
1942  Raddam_data[i].S2FITSIGMA=parm[2];
1943  Raddam_data[i].S2CHI2=S2[i]->GetFunction("fitFunc")->GetChisquare();
1944  Raddam_data[i].S2NDF=S2[i]->GetFunction("fitFunc")->GetNDF();
1945  Raddam_data[i].S2BINWIDTH=Ws2;
1946  }
1947  if(XmlFilePath.size()>0){
1948  char TIME[40];
1949  Long_t t; t=time(0); strftime(TIME,30,"%F %T",localtime(&t));
1950  //create XML file
1951  if(!Overwrite){
1952  sprintf(str,"HcalDetDiagRaddam_%i_%i.xml",run_number,dataset_seq_number);
1953  }else{
1954  sprintf(str,"HcalDetDiagRaddam.xml");
1955  }
1956  std::string xmlName=str;
1957  std::ofstream xmlFile;
1958  xmlFile.open(xmlName.c_str());
1959  xmlFile<<"<?xml version=\"1.0\" encoding=\"UTF-8\" standalone=\"no\" ?>\n";
1960  xmlFile<<"<ROOT>\n";
1961  xmlFile<<" <HEADER>\n";
1962  xmlFile<<" <TYPE>\n";
1963  xmlFile<<" <EXTENSION_TABLE_NAME>HCAL_RADDAM</EXTENSION_TABLE_NAME>\n";
1964  xmlFile<<" <NAME>HCAL Raddam</NAME>\n";
1965  xmlFile<<" </TYPE>\n";
1966  xmlFile<<" <!-- run details -->\n";
1967  xmlFile<<" <RUN>\n";
1968  xmlFile<<" <RUN_TYPE>TEST LOCAL-RUN</RUN_TYPE>\n";
1969  xmlFile<<" <RUN_NUMBER>"<<run_number<<"</RUN_NUMBER>\n";
1970  xmlFile<<" <RUN_BEGIN_TIMESTAMP>"<<TIME<<"</RUN_BEGIN_TIMESTAMP>\n";
1971  xmlFile<<" <COMMENT_DESCRIPTION>hcal raddam data</COMMENT_DESCRIPTION>\n";
1972  xmlFile<<" <LOCATION>P5</LOCATION>\n";
1973  xmlFile<<" <INITIATED_BY_USER>dma</INITIATED_BY_USER>\n";
1974  xmlFile<<" </RUN>\n";
1975  xmlFile<<" </HEADER>\n";
1976  xmlFile<<" <DATA_SET>\n";
1977  xmlFile<<" <COMMENT_DESCRIPTION>Test Raddam data</COMMENT_DESCRIPTION>\n";
1978  xmlFile<<" <CREATE_TIMESTAMP>"<<TIME<<"</CREATE_TIMESTAMP>\n";
1979  xmlFile<<" <CREATED_BY_USER>dma</CREATED_BY_USER>\n";
1980  xmlFile<<" <VERSION>Test_Version_1</VERSION>\n";
1981 
1982  for(int i=0;i<56;i++){
1983  xmlFile<<" <DATA>\n";
1984  xmlFile<<" <SUBDET>HF</SUBDET>\n";
1985  xmlFile<<" <IETA>"<<RADDAM_CH[i].eta<<"</IETA>\n";
1986  xmlFile<<" <IPHI>"<<RADDAM_CH[i].phi<<"</IPHI>\n";
1987  xmlFile<<" <DEPTH>"<<RADDAM_CH[i].depth<<"</DEPTH>\n";
1988 
1989  xmlFile<<" <TOTEVNT>"<<Raddam_data[i].TOTEVNT<<"</TOTEVNT>\n";
1990  xmlFile<<" <CUT1EVNT>"<<Raddam_data[i].CUT1EVNT<<"</CUT1EVNT>\n";
1991  xmlFile<<" <CUT2EVNT>"<<Raddam_data[i].CUT2EVNT<<"</CUT2EVNT>\n";
1992 
1993  xmlFile<<" <S1MEAN>"<<Raddam_data[i].S1MEAN <<"</S1MEAN>\n";
1994  xmlFile<<" <S1RMS>"<<Raddam_data[i].S1RMS <<"</S1RMS>\n";
1995  xmlFile<<" <S1FITMEAN>"<<Raddam_data[i].S1FITMEAN <<"</S1FITMEAN>\n";
1996  xmlFile<<" <S1FITMEANER>"<<Raddam_data[i].S1FITMEANER <<"</S1FITMEANER>\n";
1997  xmlFile<<" <S1FITSIGMA>"<<Raddam_data[i].S1FITSIGMA <<"</S1FITSIGMA>\n";
1998  xmlFile<<" <S1CHI2>"<<Raddam_data[i].S1CHI2 <<"</S1CHI2>\n";
1999  xmlFile<<" <S1NDF>"<<Raddam_data[i].S1NDF <<"</S1NDF>\n";
2000  xmlFile<<" <S1BINWIDTH>"<<Raddam_data[i].S1BINWIDTH <<"</S1BINWIDTH>\n";
2001 
2002  xmlFile<<" <S2MEAN>"<<Raddam_data[i].S2MEAN <<"</S2MEAN>\n";
2003  xmlFile<<" <S2RMS>"<<Raddam_data[i].S2RMS <<"</S2RMS>\n";
2004  xmlFile<<" <S2FITMEAN>"<<Raddam_data[i].S2FITMEAN <<"</S2FITMEAN>\n";
2005  xmlFile<<" <S2FITMEANER>"<<Raddam_data[i].S2FITMEANER <<"</S2FITMEANER>\n";
2006  xmlFile<<" <S2FITSIGMA>"<<Raddam_data[i].S2FITSIGMA <<"</S2FITSIGMA>\n";
2007  xmlFile<<" <S2CHI2>"<<Raddam_data[i].S2CHI2 <<"</S2CHI2>\n";
2008  xmlFile<<" <S2NDF>"<<Raddam_data[i].S2NDF <<"</S2NDF>\n";
2009  xmlFile<<" <S2BINWIDTH>"<<Raddam_data[i].S2BINWIDTH <<"</S2BINWIDTH>\n";
2010  xmlFile<<" </DATA>\n";
2011  }
2012  xmlFile<<" </DATA_SET>\n";
2013  xmlFile<<"</ROOT>\n";
2014 
2015  xmlFile.close();
2016  sprintf(str,"zip %s.zip %s",xmlName.c_str(),xmlName.c_str());
2017  system(str);
2018  sprintf(str,"rm -f %s",xmlName.c_str());
2019  system(str);
2020  sprintf(str,"mv -f %s.zip %s",xmlName.c_str(),XmlFilePath.c_str());
2021  system(str);
2022  }
2023  if(OutputFilePath.size()>0){
2024  if(!Overwrite){
2025  sprintf(str,"%sHcalDetDiagRaddamData_run%06i_%i.root",OutputFilePath.c_str(),run_number,dataset_seq_number);
2026  }else{
2027  sprintf(str,"%sHcalDetDiagRaddamData.root",OutputFilePath.c_str());
2028  }
2029  TFile *theFile = new TFile(str, "RECREATE");
2030  if(!theFile->IsOpen()) return;
2031  theFile->cd();
2032  for(int i=0;i<56;i++){
2033  if(S1[i]!=0)S1[i]->Write();
2034  if(S2[i]!=0)S2[i]->Write();
2035  }
2036  theFile->Write();
2037  theFile->Close();
2038  theFile->Delete();
2039  }
2040 }
2041 
2042 void HcalDetDiagLaserMonitor::endRun(const edm::Run& run, const edm::EventSetup& c) {
2043  if(RaddamRun){
2044  SaveRaddamData();
2045  }
2046  if((LocalRun || !Online_ || createHTMLonly) && ievt_>10){
2047  fillHistos(HcalBarrel);
2048  fillHistos(HcalOuter);
2049  fillHistos(HcalForward);
2050  fillProblems(HcalBarrel);
2051  fillProblems(HcalEndcap);
2052  fillProblems(HcalOuter);
2053  fillProblems(HcalForward);
2054  if(!RaddamRun)SaveReference();
2055  }
2056 }
2057 
2058 void HcalDetDiagLaserMonitor::beginLuminosityBlock(const edm::LuminosityBlock& lumiSeg,const edm::EventSetup& c){}
2059 void HcalDetDiagLaserMonitor::endLuminosityBlock(const edm::LuminosityBlock& lumiSeg,const edm::EventSetup& c){}
2060 
2061 //define this as a plug-in
2062 DEFINE_FWK_MODULE(HcalDetDiagLaserMonitor);
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
HcalDetDiagLaserMonitor::analyze
void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: HcalDetDiagLaserMonitor.cc:580
edm::ParameterSet::getParameter
T getParameter(std::string const &) const
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
HcalDetDiagLaserMonitor::tok_hbhe_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
Definition: HcalDetDiagLaserMonitor.cc:297
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:854
HcalDetDiagLaserData::add_statistics
void add_statistics(double *data, int nTS)
Definition: HcalDetDiagLaserMonitor.cc:137
HcalDetDiagLaserMonitor::endRun
void endRun(const edm::Run &run, const edm::EventSetup &c) override
Definition: HcalDetDiagLaserMonitor.cc:2042
cuy.ib
int ib
Definition: cuy.py:660
HcalElectronicsMap::allElectronicsIdPrecision
std::vector< HcalElectronicsId > allElectronicsIdPrecision() const
Definition: HcalElectronicsMap.cc:153
HcalDetDiagLaserData::GetEnergy
double GetEnergy(double *data, int n)
Definition: HcalDetDiagLaserMonitor.cc:209
HcalDetDiagLaserMonitor::GetCalib
HcalDetDiagLaserData * GetCalib(std::string sd, int eta, int phi)
Definition: HcalDetDiagLaserMonitor.cc:257
HcalBaseDQMonitor::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:244
HcalDetDiagLaserData::set_reference
void set_reference(float val, float rms, float time, float time_rms)
Definition: HcalDetDiagLaserMonitor.cc:145
MonitorElement::setAxisRange
void setAxisRange(double xmin, double xmax, int axis=1)
set x-, y- or z-axis range (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:948
HcalDetDiagLaserMonitor::tok_calib_
edm::EDGetTokenT< HcalCalibDigiCollection > tok_calib_
Definition: HcalDetDiagLaserMonitor.cc:295
HcalDetDiagLaserMonitor::tok_raw_
edm::EDGetTokenT< FEDRawDataCollection > tok_raw_
Definition: HcalDetDiagLaserMonitor.cc:294
HcalDetDiagLaserData::get_average_amp
bool get_average_amp(double *ave, double *rms)
Definition: HcalDetDiagLaserMonitor.cc:165
HcalDetDiagLaserData::get_average_amp1
bool get_average_amp1(double *ave, double *rms)
Definition: HcalDetDiagLaserMonitor.cc:190
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:446
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
Raddam_ch::phi
int phi
Definition: HcalDetDiagLaserClient.cc:21
HcalDetDiagLaserMonitor::emap
const HcalElectronicsMap * emap
Definition: HcalDetDiagLaserMonitor.cc:291
DQMStore::IBooker::bookInt
MonitorElement * bookInt(Args &&...args)
Definition: DQMStore.h:103
edm::SortedCollection< HFDataFrame >::const_iterator
std::vector< HFDataFrame >::const_iterator const_iterator
Definition: SortedCollection.h:82
DQMStore::IBooker::bookString
MonitorElement * bookString(Args &&...args)
Definition: DQMStore.h:97
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:930
HcalCondObjectContainer::getValues
const Item * getValues(DetId fId, bool throwOnFail=true) const
Definition: HcalCondObjectContainer.h:132
HcalDetDiagLaserMonitor::get_ave_subdet
bool get_ave_subdet(int sd, float *ave_t, float *ave_e, float *ave_t_r, float *ave_e_r)
Definition: HcalDetDiagLaserMonitor.cc:771
HcalDetDiagLaserMonitor::tok_hf_
edm::EDGetTokenT< HFDigiCollection > tok_hf_
Definition: HcalDetDiagLaserMonitor.cc:299
Raddam_data
HcalRaddamData Raddam_data[56]
Definition: HcalDetDiagLaserMonitor.cc:111
eta
T eta() const
Definition: Basic3DVectorLD.h:177
RADDAM_CH
static const Raddam_ch RADDAM_CH[56]
Definition: HcalDetDiagLaserMonitor.cc:79
adc2fC
static const float adc2fC[128]
Definition: HcalTimingMonitorModule.cc:63
HcalDetDiagLaserMonitor::get_ave_rbx
bool get_ave_rbx(int sd, int side, int rbx, float *ave, float *rms)
Definition: HcalDetDiagLaserMonitor.cc:1002
Raddam_ch::eta
int eta
Definition: HcalDetDiagLaserClient.cc:20
FEDRawData::size
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:47
HcalDetDiagLaserData::get_average_time
bool get_average_time(double *ave, double *rms)
Definition: HcalDetDiagLaserMonitor.cc:170
HcalDetDiagLaserMonitor::KnownBadCells_
std::map< unsigned int, int > KnownBadCells_
Definition: HcalDetDiagLaserMonitor.cc:378
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:164
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
SD
static const float SD[]
Definition: sicif.h:67
isSiPM
bool isSiPM(int ieta, int iphi, int depth)
Definition: HcalEtaPhiHists.h:489
HcalDetId::depth
int depth() const
get the tower depth
Definition: HcalDetId.h:40
ETA
#define ETA
Definition: GenericBenchmark.cc:32
roll_playback.k
tuple k
Definition: roll_playback.py:143
calib
MVATrainerComputer * calib
Definition: MVATrainer.cc:64
iEvent
int iEvent
Definition: GenABIO.cc:230
EtaPhiHists::depth
std::vector< MonitorElement * > depth
Definition: HcalEtaPhiHists.h:183
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
HcalCondObjectContainer::getAllChannels
std::vector< DetId > getAllChannels() const
Definition: HcalCondObjectContainer.h:242
mathSSE::sqrt
T sqrt(T t)
Definition: SSEVec.h:48
HcalDetDiagLaserMonitor::beginLuminosityBlock
void beginLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c) override
Definition: HcalDetDiagLaserMonitor.cc:2058
HcalDetId::ieta
int ieta() const
get the cell ieta
Definition: HcalDetId.h:36
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
j
int j
Definition: DBlmapReader.cc:9
f
double f[11][100]
Definition: MuScleFitUtils.cc:78
EtaPhiHists::setup
void setup(DQMStore::IBooker &m_dbe, std::string Name, std::string Units="")
Definition: HcalEtaPhiHists.h:24
MonitorElement::setBinError
void setBinError(int binx, double error)
set uncertainty on content of bin (1-D)
Definition: MonitorElement.cc:879
edm::EventBase::orbitNumber
int orbitNumber() const
Definition: EventBase.h:63
min
T min(T a, T b)
Definition: MathUtil.h:58
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
DetId::subdetId
int subdetId() const
get the contents of the subdetector field (not cast into any detector&#39;s numbering enum) ...
Definition: DetId.h:37
HcalDetDiagLaserMonitor::bookHistograms
void bookHistograms(DQMStore::IBooker &ib, const edm::Run &run, const edm::EventSetup &c) override
Definition: HcalDetDiagLaserMonitor.cc:436
HcalDetDiagLaserMonitor::tok_ho_
edm::EDGetTokenT< HODigiCollection > tok_ho_
Definition: HcalDetDiagLaserMonitor.cc:298
step1_ZMM_7Tev.Status
tuple Status
Definition: step1_ZMM_7Tev.py:112
HcalDetId::iphi
int iphi() const
get the cell iphi
Definition: HcalDetId.h:38
DetId
Definition: DetId.h:18
HcalDetDiagLaserMonitor::hb_data
HcalDetDiagLaserData hb_data[85][72][4]
Definition: HcalDetDiagLaserMonitor.cc:372
HcalDetDiagLaserMonitor::HcalDetDiagLaserMonitor
HcalDetDiagLaserMonitor(const edm::ParameterSet &)
Definition: HcalDetDiagLaserMonitor.cc:381
HcalDetDiagLaserData::get_average_time1
bool get_average_time1(double *ave, double *rms)
Definition: HcalDetDiagLaserMonitor.cc:195
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:274
N
#define N
Definition: blowfish.cc:9
HcalDetDiagLaserMonitor::endLuminosityBlock
void endLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c) override
Definition: HcalDetDiagLaserMonitor.cc:2059
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
sd
double sd
Definition: CascadeWrapper.h:137
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
PHI
#define PHI
Definition: GenericBenchmark.cc:33
Statistic
Statistic
Definition: trackSplitPlot.h:33
edm::ESHandle::product
T const * product() const
Definition: ESHandle.h:86
HcalBaseDQMonitor::reset
virtual void reset(void)
Definition: HcalBaseDQMonitor.cc:127
xmlFile
bool xmlFile(const std::string fParam)
Definition: hcalCalibrationsCopy.cc:214
CalcEtaBin
int CalcEtaBin(int subdet, int ieta, int depth)
Definition: HcalEtaPhiHists.h:202
HcalDetDiagLaserData::GetTime
double GetTime(double *data, int n=10)
Definition: HcalDetDiagLaserMonitor.cc:219
AlCaHLTBitMon_ParallelJobs.p
tuple p
Definition: AlCaHLTBitMon_ParallelJobs.py:152
Raddam_ch::depth
int depth
Definition: HcalDetDiagLaserMonitor.cc:76
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
HcalDetDiagLaserMonitor::hf_data
HcalDetDiagLaserData hf_data[85][72][4]
Definition: HcalDetDiagLaserMonitor.cc:375
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
patZpeak.events
tuple events
Definition: patZpeak.py:19
FEDRawData::data
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:28
reco::JetExtendedAssociation::getValue
const JetExtendedData & getValue(const Container &, const reco::JetBaseRef &)
get value for the association. Throw exception if no association found
Definition: JetExtendedAssociation.cc:51
HcalDetDiagLaserData::get_reference
bool get_reference(double *val, double *rms, double *time, double *time_rms)
Definition: HcalDetDiagLaserMonitor.cc:160
HcalDetDiagLaserData::set_data
void set_data(float val, float rms, float time, float time_rms)
Definition: HcalDetDiagLaserMonitor.cc:150
lumiPlot.rawdata
dictionary rawdata
Definition: lumiPlot.py:393
ntuplemaker.status
tuple status
Definition: ntuplemaker.py:245
HcalDetDiagLaserMonitor::he_data
HcalDetDiagLaserData he_data[85][72][4]
Definition: HcalDetDiagLaserMonitor.cc:373
DetId::det
Detector det() const
get the detector field from this detid
Definition: DetId.h:35
Raddam_ch
Definition: HcalDetDiagLaserClient.cc:19
T
long double T
Definition: Basic3DVectorLD.h:57
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:957
HcalDetDiagLaserMonitor::problemnames_
std::vector< std::string > problemnames_
Definition: HcalDetDiagLaserMonitor.cc:367
MonitorElement::Reset
void Reset(void)
reset ME (ie. contents, errors, etc)
Definition: MonitorElement.cc:463
HcalDetDiagLaserMonitor::ho_data
HcalDetDiagLaserData ho_data[85][72][4]
Definition: HcalDetDiagLaserMonitor.cc:374
HcalBaseDQMonitor::setup
virtual void setup(DQMStore::IBooker &)
Definition: HcalBaseDQMonitor.cc:151
HcalDetDiagLaserMonitor::tok_tb_
edm::EDGetTokenT< HcalTBTriggerData > tok_tb_
Definition: HcalDetDiagLaserMonitor.cc:296
HcalDetDiagLaserMonitor::calib_data
HcalDetDiagLaserData calib_data[5][5][72]
Definition: HcalDetDiagLaserMonitor.cc:376
HcalElectronicsMap::lookup
const DetId lookup(HcalElectronicsId fId) const
lookup the logical detid associated with the given electronics id
Definition: HcalElectronicsMap.cc:104
TTT
TrackTransientTrack TTT
Definition: TransientTrack.cc:9
edm::Run
Definition: Run.h:41
phi
Definition: DDAxes.h:10
Generated for CMSSW Reference Manual by doxygen 1.8.5