CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalDetDiagNoiseMonitor.cc
Go to the documentation of this file.
1 #include "DQM/HcalMonitorTasks/interface/HcalDetDiagNoiseMonitor.h"
2 #include "DataFormats/FEDRawData/interface/FEDRawData.h"
3 #include "DataFormats/FEDRawData/interface/FEDNumbering.h"
4 #include "DataFormats/HcalDigi/interface/HcalCalibrationEventTypes.h"
5 #include "EventFilter/HcalRawToDigi/interface/HcalDCCHeader.h"
6 #include "FWCore/Common/interface/TriggerNames.h"
7 
8 #include "TFile.h"
9 #include "TTree.h"
10 #include <TVector2.h>
11 #include <TVector3.h>
12 
13 #include "FWCore/MessageLogger/interface/MessageLogger.h"
14 
15 #include "DataFormats/Math/interface/deltaR.h"
16 #include "DataFormats/Common/interface/TriggerResults.h"
17 
18 #include "DataFormats/L1Trigger/interface/L1EmParticle.h"
19 #include "DataFormats/L1Trigger/interface/L1EmParticleFwd.h"
20 #include "DataFormats/L1Trigger/interface/L1MuonParticle.h"
21 #include "DataFormats/L1Trigger/interface/L1MuonParticleFwd.h"
22 #include "DataFormats/L1Trigger/interface/L1JetParticle.h"
23 #include "DataFormats/L1Trigger/interface/L1JetParticleFwd.h"
24 #include "DataFormats/L1Trigger/interface/L1EtMissParticle.h"
25 #include "DataFormats/L1Trigger/interface/L1EtMissParticleFwd.h"
26 #include "DataFormats/L1GlobalCaloTrigger/interface/L1GctHFRingEtSums.h"
27 #include "DataFormats/L1GlobalCaloTrigger/interface/L1GctHFBitCounts.h"
28 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTExtendedCand.h"
29 #include "DataFormats/Candidate/interface/Candidate.h"
30 #include "L1Trigger/RegionalCaloTrigger/interface/L1RCTProducer.h"
31 #include "DataFormats/L1CaloTrigger/interface/L1CaloCollections.h"
32 #include "DataFormats/EcalDigi/interface/EcalDigiCollections.h"
33 #include "HLTrigger/HLTanalyzers/interface/JetUtil.h"
34 #include "DataFormats/METReco/interface/CaloMET.h"
35 #include "DataFormats/METReco/interface/CaloMETCollection.h"
36 // #include "DataFormats/L1Trigger/interface/L1ParticleMap.h"
37 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutSetupFwd.h"
38 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapRecord.h"
39 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMapFwd.h"
40 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerObjectMap.h"
41 //#include "DataFormats/L1GlobalTrigger/interface/L1GtLogicParser.h"
42 #include "DataFormats/TrackReco/interface/Track.h"
43 #include "DataFormats/TrackReco/interface/TrackFwd.h"
44 #include "DataFormats/TrackReco/interface/TrackBase.h"
45 #include "RecoMET/METAlgorithms/interface/HcalNoiseRBXArray.h"
46 #include "DataFormats/VertexReco/interface/Vertex.h"
47 #include "DataFormats/VertexReco/interface/VertexFwd.h"
48 #include "DataFormats/HcalRecHit/interface/HcalRecHitCollections.h"
49 #include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
50 #include "FWCore/ParameterSet/interface/ParameterSet.h"
51 #include "DQMServices/Core/interface/DQMStore.h"
52 #include "DQMServices/Core/interface/MonitorElement.h"
53 #include "CalibFormats/HcalObjects/interface/HcalDbService.h"
54 #include "CalibFormats/HcalObjects/interface/HcalDbRecord.h"
55 
56 // this is to retrieve HCAL LogicalMap
57 #include "CondFormats/HcalObjects/interface/HcalLogicalMap.h"
58 
59 
60 #include <math.h>
61 
62 using namespace reco;
63 
65 constexpr 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,
66  13.5,15.,17.,19.,21.,23.,25.,27.,29.5,32.5,35.5,38.5,42.,46.,50.,54.5,59.5,
67  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,
68  124.5,129.5,137.,147.,157.,167.,177.,187.,197.,209.5,224.5,239.5,254.5,272.,
69  292.,312.,334.5,359.5,384.5,359.5,384.5,409.5,434.5,459.5,484.5,509.5,534.5,
70  559.5,584.5,609.5,634.5,659.5,684.5,709.5,747.,797.,847.,897.,947.,997.,
71  1047.,1109.5,1184.5,1259.5,1334.5,1422.,1522.,1622.,1734.5,1859.5,1984.5,
72  1859.5,1984.5,2109.5,2234.5,2359.5,2484.5,2609.5,2734.5,2859.5,2984.5,
73  3109.5,3234.5,3359.5,3484.5,3609.5,3797.,4047.,4297.,4547.,4797.,5047.,
74  5297.,5609.5,5984.5,6359.5,6734.5,7172.,7672.,8172.,8734.5,9359.5,9984.5};
76 constexpr char const *subdets[11]={"HBM","HBP","HEM","HEP","HO1M","HO0","HO1P","HO2M","HO2P","HFM","HFP"};
77 constexpr char const *HB_RBX[36]={
78 "HBM01","HBM02","HBM03","HBM04","HBM05","HBM06","HBM07","HBM08","HBM09","HBM10","HBM11","HBM12","HBM13","HBM14","HBM15","HBM16","HBM17","HBM18",
79 "HBP01","HBP02","HBP03","HBP04","HBP05","HBP06","HBP07","HBP08","HBP09","HBP10","HBP11","HBP12","HBP13","HBP14","HBP15","HBP16","HBP17","HBP18"};
80 constexpr char const *HE_RBX[36]={
81 "HEM01","HEM02","HEM03","HEM04","HEM05","HEM06","HEM07","HEM08","HEM09","HEM10","HEM11","HEM12","HEM13","HEM14","HEM15","HEM16","HEM17","HEM18",
82 "HEP01","HEP02","HEP03","HEP04","HEP05","HEP06","HEP07","HEP08","HEP09","HEP10","HEP11","HEP12","HEP13","HEP14","HEP15","HEP16","HEP17","HEP18"};
83 constexpr char const *HO_RBX[36]={
84 "HO2M02","HO2M04","HO2M06","HO2M08","HO2M10","HO2M12","HO1M02","HO1M04","HO1M06","HO1M08","HO1M10","HO1M12",
85 "HO001","HO002","HO003","HO004","HO005","HO006","HO007","HO008","HO009","HO010","HO011","HO012",
86 "HO1P02","HO1P04","HO1P06","HO1P08","HO1P10","HO1P12","HO2P02","HO2P04","HO2P06","HO2P08","HO2P10","HO2P12",
87 };
88 
89 
90 class HcalDetDiagNoiseRMData{
91 public:
92  HcalDetDiagNoiseRMData(){
93  reset();
94  reset_LS();
95  };
96  void reset(){
97  n_th_hi=n_th_300=n_pix_1=n_pix_8=pix=n_pix=0;
98  }
99  void reset_LS(){
100  n_th_hi_LS=n_th_300_LS=0;
101  }
102  int n_th_hi;
103  int n_th_300;
104  int n_pix_1;
105  int n_pix_8;
106  int pix;
107  int n_pix;
108  int n_th_hi_LS;
109  int n_th_300_LS;
110 };
111 class HcalDetDiagNoiseRMEvent{
112 public:
113  HcalDetDiagNoiseRMEvent(){
114  reset();
115  }
116  void reset(){
117  n_pix_hi=n_pix_lo=energy=n_zero=0;
118  }
119  int n_pix_hi;
120  int n_pix_lo;
121  int n_zero;
122  float energy;
123 };
124 
125 class HcalDetDiagNoiseRMSummary{
126 public:
127  HcalDetDiagNoiseRMSummary(){
128  reset();
129  }
130  void reset(){
131  for(int i=0;i<HcalFrontEndId::maxRmIndex;i++) rm[i].reset();
132  }
133  void reset_LS(){
134  for(int i=0;i<HcalFrontEndId::maxRmIndex;i++) rm[i].reset_LS();
135  }
136  int GetRMindex(const std::string& rbx,int rm){
137  if(rbx.substr(0,3)=="HO0"){
138  int sect=atoi(rbx.substr(3,2).c_str());
139  if(sect>12) return -1;
140  if(rm==1 && (sect==2 || sect==3 || sect==6 || sect==7 || sect==10 || sect==11)) return -1;
141  if(rm==4 && (sect==12 || sect==1 || sect==4 || sect==5 || sect==8 || sect==9 )) return -1;
142  }
143  if(rbx.substr(0,3)=="HO1" || rbx.substr(0,3)=="HO2"){
144  int sect=atoi(rbx.substr(4,2).c_str());
145  if(sect>12) return -1;
146  if(sect==1 || sect==3 || sect==5 || sect==7 || sect==9 || sect==11) return -1;
147  }
148  HcalFrontEndId id(rbx,rm,1,1,1,1,1);
149  if(id.rawId()==0) return -1;
150  return id.rmIndex();
151  }
152  HcalDetDiagNoiseRMData rm[HcalFrontEndId::maxRmIndex];
153 };
154 
155 HcalDetDiagNoiseMonitor::HcalDetDiagNoiseMonitor(const edm::ParameterSet& ps):
156  HcalBaseDQMonitor(ps)
157  {
158 
159  tok_tb_ = consumes<HcalTBTriggerData>(ps.getParameter<edm::InputTag>("hcalTBTriggerDataTag"));
160 
161  ievt_=0;
162  run_number=-1;
163  NoisyEvents=0;
164  LocalRun=false;
165  dataset_seq_number=1;
166  FirstOrbit=FirstOrbitLS=0xFFFFFFFF;
167  LastOrbit=LastOrbitLS=0;
168 
169  Online_ = ps.getUntrackedParameter<bool>("online",false);
170  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns",false);
171  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup",false);
172  debug_ = ps.getUntrackedParameter<int>("debug",0);
173  prefixME_ = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
174  if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
175  prefixME_.append("/");
176  subdir_ = ps.getUntrackedParameter<std::string>("TaskFolder","DetDiagNoiseMonitor_Hcal");
177  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
178  subdir_.append("/");
179  subdir_=prefixME_+subdir_;
180  AllowedCalibTypes_ = ps.getUntrackedParameter<std::vector<int> > ("AllowedCalibTypes");
181  skipOutOfOrderLS_ = ps.getUntrackedParameter<bool>("skipOutOfOrderLS",false);
182  NLumiBlocks_ = ps.getUntrackedParameter<int>("NLumiBlocks",4000);
183  makeDiagnostics_ = ps.getUntrackedParameter<bool>("makeDiagnostics",false);
184 
185  UseDB = ps.getUntrackedParameter<bool> ("UseDB" , false);
186  OutputFilePath = ps.getUntrackedParameter<std::string>("OutputFilePath", "");
187  HPDthresholdHi = ps.getUntrackedParameter<double>("NoiseThresholdHPDhi",49.0);
188  HPDthresholdLo = ps.getUntrackedParameter<double>("NoiseThresholdHPDlo",10.0);
189  SpikeThreshold = ps.getUntrackedParameter<double>("NoiseSpikeThreshold",0.5);
190  Overwrite = ps.getUntrackedParameter<bool> ("Overwrite",true);
191 
192  tok_raw_ = consumes<FEDRawDataCollection>(ps.getUntrackedParameter<edm::InputTag>("RawDataLabel",edm::InputTag("source","")));
193  digiLabel_ = ps.getUntrackedParameter<edm::InputTag>("digiLabel",edm::InputTag("hcalDigis"));
194  tok_l1_ = consumes<L1GlobalTriggerReadoutRecord>(ps.getUntrackedParameter<edm::InputTag>("gtLabel"));
195 
196 
197  tok_hbhe_ = consumes<HBHEDigiCollection>(digiLabel_);
198  tok_ho_ = consumes<HODigiCollection>(digiLabel_);
199 
200  RMSummary = 0;
201  needLogicalMap_=true;
202  setupDone_ = false;
203 }
204 
205 void HcalDetDiagNoiseMonitor::reset(){}
206 
207 
208 void HcalDetDiagNoiseMonitor::bookHistograms(DQMStore::IBooker &ib, const edm::Run& run, const edm::EventSetup& c)
209 {
210  if (debug_>1) std::cout <<"HcalDetDiagNoiseMonitor::bookHistograms"<<std::endl;
211  HcalBaseDQMonitor::bookHistograms(ib,run,c);
212 
213  if (tevt_==0) this->setup(ib); // set up histograms if they have not been created before
214  if (mergeRuns_==false)
215  this->reset();
216 
217  return;
218 
219 }
220 
221 void HcalDetDiagNoiseMonitor::setup(DQMStore::IBooker &ib){
222  if (setupDone_)
223  return;
224  setupDone_ = true;
225  // Call base class setup
226  HcalBaseDQMonitor::setup(ib);
227  RMSummary = new HcalDetDiagNoiseRMSummary();
228 
229  std::string name;
230  ib.setCurrentFolder(subdir_);
231  meEVT_ = ib.bookInt("HcalNoiseMonitor Event Number");
232  ib.setCurrentFolder(subdir_+"Common Plots");
233 
234  name="RBX Pixel multiplicity"; PixelMult = ib.book1D(name,name,73,0,73);
235  name="HPD energy"; HPDEnergy = ib.book1D(name,name,200,0,2500);
236  name="RBX energy"; RBXEnergy = ib.book1D(name,name,200,0,3500);
237  name="Number of zero TS per RBX"; NZeroes = ib.book1D(name,name,100,0,100);
238  name="Trigger BX Tbit11"; TriggerBx11 = ib.book1D(name,name,4000,0,4000);
239  name="Trigger BX Tbit12"; TriggerBx12 = ib.book1D(name,name,4000,0,4000);
240 
241  ib.setCurrentFolder(subdir_+"HBHE Plots");
242  name="HBP HPD Noise Rate Pixel above 50fC"; HBP_Rate50 = ib.book1D(name,name,73,0,73);
243  name="HBM HPD Noise Rate Pixel above 50fC"; HBM_Rate50 = ib.book1D(name,name,73,0,73);
244  name="HEP HPD Noise Rate Pixel above 50fC"; HEP_Rate50 = ib.book1D(name,name,73,0,73);
245  name="HEM HPD Noise Rate Pixel above 50fC"; HEM_Rate50 = ib.book1D(name,name,73,0,73);
246  name="HBP HPD Noise Rate HPD above 300fC"; HBP_Rate300 = ib.book1D(name,name,73,0,73);
247  name="HBM HPD Noise Rate HPD above 300fC"; HBM_Rate300 = ib.book1D(name,name,73,0,73);
248  name="HEP HPD Noise Rate HPD above 300fC"; HEP_Rate300 = ib.book1D(name,name,73,0,73);
249  name="HEM HPD Noise Rate HPD above 300fC"; HEM_Rate300 = ib.book1D(name,name,73,0,73);
250 
251  ib.setCurrentFolder(subdir_+"HO Plots");
252  name="HO0 HPD Noise Rate Pixel above 50fC"; HO0_Rate50 = ib.book1D(name,name,49,0,49);
253  name="HO1P HPD Noise Rate Pixel above 50fC"; HO1P_Rate50 = ib.book1D(name,name,48,0,48);
254  name="HO1M HPD Noise Rate Pixel above 50fC"; HO1M_Rate50 = ib.book1D(name,name,48,0,48);
255  name="HO0 HPD Noise Rate HPD above 300fC"; HO0_Rate300 = ib.book1D(name,name,48,0,48);
256  name="HO1P HPD Noise Rate HPD abGetRMindexove 300fC"; HO1P_Rate300 = ib.book1D(name,name,48,0,48);
257  name="HO1M HPD Noise Rate HPD above 300fC"; HO1M_Rate300 = ib.book1D(name,name,48,0,48);
258 
259 
260  ib.setCurrentFolder(subdir_+"Noise Spike Plots");
261 
262  name="HB RM Spike Map"; HB_RBXmapSpikeCnt= ib.book2D(name,name,4,0.5,4.5,36,0.5,36.5);
263  name="HE RM Spike Map"; HE_RBXmapSpikeCnt= ib.book2D(name,name,4,0.5,4.5,36,0.5,36.5);
264  name="HO RM Spike Map"; HO_RBXmapSpikeCnt= ib.book2D(name,name,4,0.5,4.5,36,0.5,36.5);
265 
266  std::string title="RM";
267  HB_RBXmapSpikeCnt->setAxisTitle(title);
268  HE_RBXmapSpikeCnt->setAxisTitle(title);
269  HO_RBXmapSpikeCnt->setAxisTitle(title);
270 
271  for(int i=0;i<36;i++){
272  HB_RBXmapSpikeCnt->setBinLabel(i+1,HB_RBX[i],2);
273  HE_RBXmapSpikeCnt->setBinLabel(i+1,HE_RBX[i],2);
274  HO_RBXmapSpikeCnt->setBinLabel(i+1,HO_RBX[i],2);
275  }
276 
277 
278  return;
279 }
280 
281 void HcalDetDiagNoiseMonitor::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup){
282  getLogicalMap(iSetup);
283  HcalBaseDQMonitor::analyze(iEvent, iSetup);
284  if (!IsAllowedCalibType()) return;
285  if (LumiInOrder(iEvent.luminosityBlock())==false) return;
286 // HcalBaseDQMonitor::analyze(iEvent, iSetup);
287  bool isNoiseEvent=false;
288  int orbit=-1111;
289  int bx=-1111;
290 
291  // for local runs
292  edm::Handle<HcalTBTriggerData> trigger_data;
293  iEvent.getByToken(tok_tb_, trigger_data);
294  if(trigger_data.isValid()){
295  if(trigger_data->triggerWord()>1000) isNoiseEvent=true;
296  LocalRun=true;
297  }
298 
299  // We do not want to look at Abort Gap events
300  edm::Handle<FEDRawDataCollection> rawdata;
301  iEvent.getByToken(tok_raw_,rawdata);
302  //checking FEDs for calibration information
303  for(int i=FEDNumbering::MINHCALFEDID;
304  i<=FEDNumbering::MAXHCALuTCAFEDID; i++)
305  {
306  if (i>FEDNumbering::MAXHCALFEDID && i<FEDNumbering::MINHCALuTCAFEDID)
307  continue;
308 
309  const FEDRawData& fedData = rawdata->FEDData(i) ;
310  if ( fedData.size() < 24 ) continue ;
311  orbit= ((const HcalDCCHeader*)(fedData.data()))->getOrbitNumber();
312  bx=((const HcalDCCHeader*)(fedData.data()))->getBunchId();
313  if(((const HcalDCCHeader*)(fedData.data()))->getCalibType()!=hc_Null) return;
314  }
315 
316  // Check GCT trigger bits
317  edm::Handle< L1GlobalTriggerReadoutRecord > gtRecord;
318  iEvent.getByToken(tok_l1_, gtRecord);
319  if(gtRecord.isValid()){
320  const TechnicalTriggerWord tWord = gtRecord->technicalTriggerWord();
321  if(tWord.at(11) || tWord.at(12)) isNoiseEvent=true;
322  if(tWord.at(11)){ TriggerBx11->Fill(bx);}
323  if(tWord.at(12)){ TriggerBx12->Fill(bx);}
324  }
325 
326  if(!isNoiseEvent) return;
327  if(ievt_==0){ FirstOrbit=orbit; FirstOrbitLS=orbit; newLS=true;}
328  if(LastOrbit <orbit) LastOrbit=orbit;
329  if(FirstOrbit>orbit) FirstOrbit=orbit;
330  if(LastOrbitLS <orbit) LastOrbitLS=orbit;
331  if(FirstOrbitLS>orbit) FirstOrbitLS=orbit;
332  if(newLS){
333  FirstOrbitLS=orbit;
334  newLS=false;
335  }
336 
337  if(!LocalRun){
338  double TIME=(double)(LastOrbit-FirstOrbit)/11223.0;
339  if(TIME>1800.0){
340  UpdateHistos();
341  SaveRates();
342  RMSummary->reset();
343  FirstOrbit=orbit;
344  }
345  }
346 
347  meEVT_->Fill(++ievt_);
348 
349  run_number=iEvent.id().run();
350 
351  HcalDetDiagNoiseRMEvent RMs[HcalFrontEndId::maxRmIndex];
352 
353  edm::Handle<HBHEDigiCollection> hbhe;
354  iEvent.getByToken(tok_hbhe_,hbhe);
355  for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
356  double max=-100/*,sum*/,energy=0; int n_zero=0;
357  for(int i=0;i<digi->size();i++){
358  // sum=adc2fC[digi->sample(i).adc()&0xff];
359  if(max<adc2fC[digi->sample(i).adc()&0xff]) max=adc2fC[digi->sample(i).adc()&0xff];
360  if(adc2fC[digi->sample(i).adc()&0xff]==0) n_zero++;
361  }
362  HcalFrontEndId lmap_entry=logicalMap_->getHcalFrontEndId(digi->id());
363  int index=lmap_entry.rmIndex(); if(index>=HcalFrontEndId::maxRmIndex) continue;
364  RMs[index].n_zero++;
365  if(max>HPDthresholdLo){
366  for(int i=0;i<digi->size();i++) energy+=adc2fC[digi->sample(i).adc()&0xff]-2.5;
367  RMs[index].n_pix_lo++;
368  if(max>HPDthresholdHi){ RMs[index].n_pix_hi++; isNoiseEvent=true;}
369  RMs[index].energy+=energy;
370  }
371  }
372 
373  edm::Handle<HODigiCollection> ho;
374  iEvent.getByToken(tok_ho_,ho);
375  for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
376  double max=-100,energy=0; int Eta=digi->id().ieta(); int Phi=digi->id().iphi(); int n_zero=0;
377  for(int i=0;i<digi->size()-1;i++){
378  if(max<adc2fC[digi->sample(i).adc()&0xff]) max=adc2fC[digi->sample(i).adc()&0xff];
379  if(adc2fC[digi->sample(i).adc()&0xff]==0) n_zero++;
380  }
381  if((Eta>=11 && Eta<=15 && Phi>=59 && Phi<=70) || (Eta>=5 && Eta<=10 && Phi>=47 && Phi<=58)){
382  continue; // ignory SiPMs
383  }else{
384  HcalFrontEndId lmap_entry=logicalMap_->getHcalFrontEndId(digi->id());
385  int index=lmap_entry.rmIndex(); if(index>=HcalFrontEndId::maxRmIndex) continue;
386  RMs[index].n_zero++;
387  if(max>HPDthresholdLo){
388  for(int i=0;i<digi->size();i++) energy+=adc2fC[digi->sample(i).adc()&0xff]-2.5;
389  RMs[index].n_pix_lo++;
390  if(max>HPDthresholdHi){ RMs[index].n_pix_hi++; isNoiseEvent=true;}
391  RMs[index].energy+=energy;
392  }
393  }
394  }
395 
396  NoisyEvents++;
397 
398  // RMs loop
399  for(int i=0;i<HcalFrontEndId::maxRmIndex;i++){
400  if(RMs[i].n_pix_hi>0){
401  HPDEnergy->Fill(RMs[i].energy);
402  RMSummary->rm[i].n_th_hi++;
403  RMSummary->rm[i].n_th_hi_LS++;
404  if(RMs[i].energy>300) RMSummary->rm[i].n_th_300++;
405  if(RMs[i].energy>300) RMSummary->rm[i].n_th_300_LS++;
406  if(RMs[i].n_pix_lo>1) RMSummary->rm[i].n_pix_1++;
407  if(RMs[i].n_pix_lo>8) RMSummary->rm[i].n_pix_8++;
408  RMSummary->rm[i].pix+=RMs[i].n_pix_lo;
409  RMSummary->rm[i].n_pix++;
410  }
411  }
412 
413  // RBX loop
414  for(int sd=0;sd<7;sd++) for(int sect=1;sect<=18;sect++){
415  std::stringstream tempss;
416  tempss << std::setw(2) << std::setfill('0') << sect;
417  std::string rbx= subdets[sd]+tempss.str();
418 
419  double rbx_energy=0;int pix_mult=0; int n_zero=0; bool isValidRBX=false;
420  for(int rm=1;rm<=4;rm++){
421  int index=RMSummary->GetRMindex(rbx,rm);
422  if(index>0 && index<HcalFrontEndId::maxRmIndex){
423  rbx_energy+=RMs[index].energy;
424  pix_mult+=RMs[index].n_pix_lo;
425  n_zero+=RMs[index].n_zero;
426  isValidRBX=true;
427  }
428  }
429  if(isValidRBX){
430  PixelMult->Fill(pix_mult);
431  RBXEnergy->Fill(rbx_energy);
432  NZeroes->Fill(n_zero);
433  }
434  }
435 
436  if((ievt_%100)==0 && debug_>0)
437  std::cout <<ievt_<<"\t"<<NoisyEvents<<std::endl;
438  return;
439 }
440 
441 void HcalDetDiagNoiseMonitor::endLuminosityBlock(const edm::LuminosityBlock& lumiSeg,const edm::EventSetup& c){
442 int first_rbx=0,last_rbx=0;
443  //double TIME=(double)(LastOrbitLS-FirstOrbitLS)/11223.0;
444  double TIME=23.0;
445  newLS=true;
446  if(TIME==0) return;
447 
448  for(int sd=0;sd<9;sd++){
449  if(sd==0){ first_rbx=0; last_rbx=18;} //HBM
450  if(sd==1){ first_rbx=18; last_rbx=36;} //HBP
451  if(sd==0 || sd==1){ // update HB plots
452  for(int rbx=first_rbx;rbx<last_rbx;rbx++)for(int rm=1;rm<=4;rm++){
453  int index=RMSummary->GetRMindex(HB_RBX[rbx],rm);
454  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
455  double val=RMSummary->rm[index].n_th_hi_LS/TIME;
456  if(val>SpikeThreshold){
457  HB_RBXmapSpikeCnt->Fill(rm,rbx+1,1);
458  //printf("%s %i %f (%f)\n",HO_RBX[rbx].c_str(),rm,RMSummary->rm[index].n_th_hi_LS/TIME,TIME);
459  }
460  RMSummary->rm[index].reset_LS();
461  }
462  }
463  if(sd==2){ first_rbx=0; last_rbx=18;} //HEM
464  if(sd==3){ first_rbx=18; last_rbx=36;} //HEP
465  if(sd==2 || sd==3){ // update HB plots
466  for(int rbx=first_rbx;rbx<last_rbx;rbx++)for(int rm=1;rm<=4;rm++){
467  int index=RMSummary->GetRMindex(HE_RBX[rbx],rm);
468  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
469  double val=RMSummary->rm[index].n_th_hi_LS/TIME;
470  if(val>SpikeThreshold){
471  HE_RBXmapSpikeCnt->Fill(rm,rbx+1,1);
472  //printf("%s %i %f (%f)\n",HO_RBX[rbx].c_str(),rm,RMSummary->rm[index].n_th_hi_LS/TIME,TIME);
473  }
474  RMSummary->rm[index].reset_LS();
475  }
476  }
477  if(sd==4){ first_rbx=6; last_rbx=12;} //HO1M
478  if(sd==5){ first_rbx=12; last_rbx=24;} //HO0
479  if(sd==6){ first_rbx=24; last_rbx=30;} //HO1P
480  if(sd>3){
481  for(int rbx=first_rbx;rbx<last_rbx;rbx++)for(int rm=1;rm<=4;rm++){
482  int index=RMSummary->GetRMindex(HO_RBX[rbx],rm);
483  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
484  double val=RMSummary->rm[index].n_th_hi_LS/TIME;
485  if(val>SpikeThreshold){
486  HO_RBXmapSpikeCnt->Fill(rm,rbx+1,1);
487  //printf("%s %i %f (%f)\n",HO_RBX[rbx].c_str(),rm,RMSummary->rm[index].n_th_hi_LS/TIME,TIME);
488  }
489  RMSummary->rm[index].reset_LS();
490  }
491  }
492  } //sd=0;sd<9
493 }
494 
495 void HcalDetDiagNoiseMonitor::UpdateHistos(){
496  int first_rbx=0;
497  double TIME=(double)(LastOrbit-FirstOrbit)/11223.0;
498  if(TIME==0) return;
499  for(int sd=0;sd<9;sd++){
500  if(sd==0){ first_rbx=0; } //HBM
501  if(sd==1){ first_rbx=18;} //HBP
502  if(sd==0 || sd==1){ // update HB plots
503  for(int rbx=0;rbx<18;rbx++)for(int rm=1;rm<=4;rm++){
504  int index=RMSummary->GetRMindex(HB_RBX[rbx+first_rbx],rm);
505  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
506  if(sd==0){
507  HBM_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
508  HBM_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
509  }
510  if(sd==1){
511  HBP_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
512  HBP_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
513  }
514  }
515  }
516  if(sd==2){ first_rbx=0;} //HEM
517  if(sd==3){ first_rbx=18;} //HEP
518  if(sd==2 || sd==3){ // update HB plots
519  for(int rbx=0;rbx<18;rbx++)for(int rm=1;rm<=4;rm++){
520  int index=RMSummary->GetRMindex(HE_RBX[rbx+first_rbx],rm);
521  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
522  if(sd==2){
523  HEM_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
524  HEM_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
525  }
526  if(sd==3){
527  HEP_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
528  HEP_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
529  }
530  }
531  }
532  int n=0;
533  if(sd==4){ first_rbx=6; n=6;} //HO1M
534  if(sd==5){ first_rbx=12;n=12;} //HO0
535  if(sd==6){ first_rbx=24;n=6;} //HO1P
536  if(sd>3){ // update HO plots
537  for(int rbx=0;rbx<n;rbx++)for(int rm=1;rm<=4;rm++){
538  int index=RMSummary->GetRMindex(HO_RBX[rbx+first_rbx],rm);
539  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
540  if(sd==4){
541  HO1M_Rate50->setBinContent(rbx*4*2+rm,RMSummary->rm[index].n_th_hi/TIME);
542  HO1M_Rate300->setBinContent(rbx*4*2+rm,RMSummary->rm[index].n_th_300/TIME);
543  }
544  if(sd==5){
545  HO0_Rate50->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_hi/TIME);
546  HO0_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
547  }
548  if(sd==5){
549  HO1P_Rate50->setBinContent(rbx*4*2+rm,RMSummary->rm[index].n_th_hi/TIME);
550  HO1P_Rate300->setBinContent(rbx*4+rm,RMSummary->rm[index].n_th_300/TIME);
551  }
552  }
553  }
554  } //sd=0;sd<9
555 }
556 
557 void HcalDetDiagNoiseMonitor::SaveRates(){
558 char RBX[20];
559 int RM;
560 float VAL1,VAL2,VAL3,VAL4,VAL5;
561 char str[500];
562  double TIME=(double)(LastOrbit-FirstOrbit)/11223.0;
563  if(TIME==0) return;
564  if(OutputFilePath.size()>0){
565  if(!Overwrite){
566  sprintf(str,"%sHcalDetDiagNoiseData_run%06i_%i.root",OutputFilePath.c_str(),run_number,dataset_seq_number);
567  }else{
568  sprintf(str,"%sHcalDetDiagNoiseData.root",OutputFilePath.c_str());
569  }
570  TFile *theFile = new TFile(str, "RECREATE");
571  if(!theFile->IsOpen()) return;
572  theFile->cd();
573  sprintf(str,"%d",run_number); TObjString run(str); run.Write("run number");
574  sprintf(str,"%d",ievt_); TObjString events(str); events.Write("Total events processed");
575  sprintf(str,"%d",dataset_seq_number); TObjString dsnum(str); dsnum.Write("Dataset number");
576  Long_t t; t=time(0); strftime(str,30,"%F %T",localtime(&t)); TObjString tm(str); tm.Write("Dataset creation time");
577 
578  TTree *tree =new TTree("HCAL Noise data","HCAL Noise data");
579  if(tree==0) return;
580  tree->Branch("RBX", &RBX, "RBX/C");
581  tree->Branch("rm", &RM, "rm/I");
582  tree->Branch("RATE_50", &VAL1, "RATE_50");
583  tree->Branch("RATE_300", &VAL2, "RATE_300");
584  tree->Branch("RATE_PIX1", &VAL3, "RATE_PIX1");
585  tree->Branch("RATE_PIX8", &VAL4, "RATE_PIX8");
586  tree->Branch("RATE_PIXMEAN", &VAL5, "RATE_PIXMEAN");
587  for(int rbx=0;rbx<36;rbx++) for(int rm=1;rm<=4;rm++){
588  int index=RMSummary->GetRMindex(HB_RBX[rbx],rm);
589  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
590  sprintf(RBX,"%s",HB_RBX[rbx]);
591  RM=rm;
592  VAL1=RMSummary->rm[index].n_th_hi/TIME;
593  VAL2=RMSummary->rm[index].n_th_300/TIME;
594  VAL3=RMSummary->rm[index].n_pix_1/TIME;
595  VAL4=RMSummary->rm[index].n_pix_8/TIME;
596  if(RMSummary->rm[index].n_pix>0)VAL5=RMSummary->rm[index].pix/RMSummary->rm[index].n_pix; else VAL5=0;
597  tree->Fill();
598  }
599  for(int rbx=0;rbx<36;rbx++) for(int rm=1;rm<=4;rm++){
600  int index=RMSummary->GetRMindex(HE_RBX[rbx],rm);
601  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
602  sprintf(RBX,"%s",HE_RBX[rbx]);
603  RM=rm;
604  VAL1=RMSummary->rm[index].n_th_hi/TIME;
605  VAL2=RMSummary->rm[index].n_th_300/TIME;
606  VAL3=RMSummary->rm[index].n_pix_1/TIME;
607  VAL4=RMSummary->rm[index].n_pix_8/TIME;
608  if(RMSummary->rm[index].n_pix>0)VAL5=RMSummary->rm[index].pix/RMSummary->rm[index].n_pix; else VAL5=0;
609  tree->Fill();
610  }
611  for(int rbx=0;rbx<36;rbx++) for(int rm=1;rm<=4;rm++){
612  int index=RMSummary->GetRMindex(HO_RBX[rbx],rm);
613  if(index<0 || index>=HcalFrontEndId::maxRmIndex) continue;
614  sprintf(RBX,"%s",HO_RBX[rbx]);
615  RM=rm;
616  VAL1=RMSummary->rm[index].n_th_hi/TIME;
617  VAL2=RMSummary->rm[index].n_th_300/TIME;
618  VAL3=RMSummary->rm[index].n_pix_1/TIME;
619  VAL4=RMSummary->rm[index].n_pix_8/TIME;
620  if(RMSummary->rm[index].n_pix>0)VAL5=RMSummary->rm[index].pix/RMSummary->rm[index].n_pix; else VAL5=0;
621  tree->Fill();
622  }
623  theFile->Write();
624  theFile->Close();
625  theFile->Delete();
626  dataset_seq_number++;
627 
628  }
629 }
630 
631 
632 void HcalDetDiagNoiseMonitor::done(){}
633 
634 HcalDetDiagNoiseMonitor::~HcalDetDiagNoiseMonitor()
635 {
636  if(LocalRun) UpdateHistos(); SaveRates();
637 
638  if ( RMSummary ) delete RMSummary;
639 }
640 
641 DEFINE_FWK_MODULE (HcalDetDiagNoiseMonitor);
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:39
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
HcalBaseDQMonitor::LumiInOrder
bool LumiInOrder(int lumisec)
Definition: HcalBaseDQMonitor.cc:203
MonitorElement::setBinContent
void setBinContent(int binx, double content)
set content of bin (1-D)
Definition: MonitorElement.cc:844
cuy.ib
int ib
Definition: cuy.py:660
HcalBaseDQMonitor::analyze
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:244
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:462
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
HcalBaseDQMonitor::AllowedCalibTypes_
std::vector< int > AllowedCalibTypes_
Definition: HcalBaseDQMonitor.h:84
HcalDetDiagNoiseRMSummary::rm
HcalDetDiagNoiseRMData rm[HcalFrontEndId::maxRmIndex]
Definition: HcalDetDiagNoiseMonitor.cc:152
HcalBaseDQMonitor::bookHistograms
virtual void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &)
Definition: HcalBaseDQMonitor.cc:107
DQMStore::IBooker::bookInt
MonitorElement * bookInt(Args &&...args)
Definition: DQMStore.h:103
edm::SortedCollection< HBHEDataFrame >::const_iterator
std::vector< HBHEDataFrame >::const_iterator const_iterator
Definition: SortedCollection.h:82
HO_RBX
char const * HO_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:83
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:920
edm::EventBase::luminosityBlock
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:62
HcalBaseDQMonitor::getLogicalMap
void getLogicalMap(const edm::EventSetup &c)
Definition: HcalBaseDQMonitor.cc:235
dt_dqm_sourceclient_common_cff.reco
tuple reco
Definition: dt_dqm_sourceclient_common_cff.py:107
HcalDetDiagNoiseMonitor::RMSummary
HcalDetDiagNoiseRMSummary * RMSummary
Definition: HcalDetDiagNoiseMonitor.h:113
constexpr
#define constexpr
FEDRawData::size
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:47
HB_RBX
char const * HB_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:77
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
HLT_25ns10e33_v2_cff.InputTag
tuple InputTag
Definition: HLT_25ns10e33_v2_cff.py:58467
iEvent
int iEvent
Definition: GenABIO.cc:230
HcalDetDiagNoiseMonitor::analyze
void analyze(edm::Event const &e, edm::EventSetup const &s)
Definition: HcalDetDiagNoiseMonitor.cc:281
HcalDetDiagNoiseRMSummary::GetRMindex
int GetRMindex(const std::string &rbx, int rm)
Definition: HcalDetDiagNoiseMonitor.cc:136
RBX_struct
Definition: myFilter.cc:50
HcalDetDiagNoiseMonitor::tok_l1_
edm::EDGetTokenT< L1GlobalTriggerReadoutRecord > tok_l1_
Definition: HcalDetDiagNoiseMonitor.h:55
HcalFrontEndId::rmIndex
int rmIndex() const
Definition: HcalFrontEndId.h:20
HcalBaseDQMonitor::logicalMap_
HcalLogicalMap * logicalMap_
Definition: HcalBaseDQMonitor.h:120
HcalDetDiagNoiseMonitor::tok_hbhe_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
Definition: HcalDetDiagNoiseMonitor.h:56
submit.rm
string rm
Definition: submit.py:76
HcalFrontEndId::maxRmIndex
static const int maxRmIndex
Definition: HcalFrontEndId.h:21
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
TechnicalTriggerWord
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
Definition: L1GlobalTriggerReadoutSetupFwd.h:41
HcalDetDiagNoiseMonitor::tok_ho_
edm::EDGetTokenT< HODigiCollection > tok_ho_
Definition: HcalDetDiagNoiseMonitor.h:57
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
subdets
char const * subdets[11]
Definition: HcalDetDiagNoiseMonitor.cc:76
HcalLogicalMap::getHcalFrontEndId
const HcalFrontEndId getHcalFrontEndId(const DetId &)
Definition: HcalLogicalMap.cc:263
HcalDetDiagNoiseMonitor::endLuminosityBlock
void endLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &c)
Definition: HcalDetDiagNoiseMonitor.cc:441
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:276
HE_RBX
char const * HE_RBX[36]
Definition: HcalDetDiagNoiseMonitor.cc:80
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
sd
double sd
Definition: CascadeWrapper.h:137
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:59
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
relval_parameters_module.energy
string energy
Definition: relval_parameters_module.py:29
gather_cfg.cout
tuple cout
Definition: gather_cfg.py:145
HcalDetDiagNoiseMonitor::tok_raw_
edm::EDGetTokenT< FEDRawDataCollection > tok_raw_
Definition: HcalDetDiagNoiseMonitor.h:53
HcalDetDiagNoiseMonitor::setup
void setup(DQMStore::IBooker &)
Definition: HcalDetDiagNoiseMonitor.cc:221
lumiPlot.rawdata
dictionary rawdata
Definition: lumiPlot.py:393
HcalDetDiagNoiseMonitor::HcalDetDiagNoiseMonitor
HcalDetDiagNoiseMonitor(const edm::ParameterSet &ps)
Definition: HcalDetDiagNoiseMonitor.cc:155
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:947
reset
void reset(double vett[256])
Definition: TPedValues.cc:11
HcalDetDiagNoiseMonitor::bookHistograms
void bookHistograms(DQMStore::IBooker &ib, const edm::Run &run, const edm::EventSetup &c)
Definition: HcalDetDiagNoiseMonitor.cc:208
HcalBaseDQMonitor::setup
virtual void setup(DQMStore::IBooker &)
Definition: HcalBaseDQMonitor.cc:151
edm::Run
Definition: Run.h:43
adc2fC
float adc2fC[128]
Definition: HcalDetDiagNoiseMonitor.cc:65
HcalDetDiagNoiseMonitor::tok_tb_
edm::EDGetTokenT< HcalTBTriggerData > tok_tb_
Definition: HcalDetDiagNoiseMonitor.h:54
Generated for CMSSW Reference Manual by doxygen 1.8.5