CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalTimingMonitorModule.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: HcalTimingMonitorModule
4 // Class: HcalTimingMonitorModule
5 //
13 //
14 // Original Author: Dmitry Vishnevskiy
15 // Created: Thu Mar 27 08:12:02 CET 2008
16 //
17 //
18 
19 static const int MAXRPC =20;
20 static const int MAXDTBX=20;
21 static const int MAXCSC =20;
22 static const int TRIG_DT =1;
23 static const int TRIG_RPC=2;
24 static const int TRIG_GCT=4;
25 static const int TRIG_CSC=8;
26 
27 // system include files
28 #include <iostream>
29 #include <fstream>
30 #include <cmath>
31 #include <iosfwd>
32 #include <bitset>
33 #include <memory>
34 
35 // user include files
36 #include "FWCore/Framework/interface/Frameworkfwd.h"
37 #include "FWCore/Framework/interface/EDAnalyzer.h"
38 #include "FWCore/Framework/interface/Event.h"
39 #include "FWCore/Framework/interface/ESHandle.h"
40 #include "FWCore/Framework/interface/MakerMacros.h"
41 #include "FWCore/ParameterSet/interface/ParameterSet.h"
42 #include "FWCore/Utilities/interface/InputTag.h"
43 
44 // this is to retrieve HCAL digi's
45 #include "DataFormats/HcalDigi/interface/HcalDigiCollections.h"
46 #include "CalibFormats/HcalObjects/interface/HcalDbService.h"
47 #include "CalibFormats/HcalObjects/interface/HcalDbRecord.h"
48 #include "CalibFormats/HcalObjects/interface/HcalCalibrations.h"
49 #include "CalibFormats/HcalObjects/interface/HcalCoderDb.h"
50 
51 // this is to retrieve GT digi's
52 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuRegionalCand.h"
53 #include "DataFormats/L1GlobalMuonTrigger/interface/L1MuGMTReadoutCollection.h"
54 #include "DataFormats/L1GlobalTrigger/interface/L1GlobalTriggerReadoutRecord.h"
55 #include "DataFormats/L1GlobalTrigger/interface/L1GtPsbWord.h"
56 #include "DataFormats/L1GlobalTrigger/interface/L1GtFdlWord.h"
57 
58 // DQM stuff
59 #include "DQMServices/Core/interface/DQMStore.h"
60 #include "DQMServices/Core/interface/MonitorElement.h"
61 #include "FWCore/ServiceRegistry/interface/Service.h"
62 
63 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,
64  13.5,15.,17.,19.,21.,23.,25.,27.,29.5,32.5,35.5,38.5,42.,46.,50.,54.5,59.5,
65  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,
66  124.5,129.5,137.,147.,157.,167.,177.,187.,197.,209.5,224.5,239.5,254.5,272.,
67  292.,312.,334.5,359.5,384.5,359.5,384.5,409.5,434.5,459.5,484.5,509.5,534.5,
68  559.5,584.5,609.5,634.5,659.5,684.5,709.5,747.,797.,847.,897.,947.,997.,
69  1047.,1109.5,1184.5,1259.5,1334.5,1422.,1522.,1622.,1734.5,1859.5,1984.5,
70  1859.5,1984.5,2109.5,2234.5,2359.5,2484.5,2609.5,2734.5,2859.5,2984.5,
71  3109.5,3234.5,3359.5,3484.5,3609.5,3797.,4047.,4297.,4547.,4797.,5047.,
72  5297.,5609.5,5984.5,6359.5,6734.5,7172.,7672.,8172.,8734.5,9359.5,9984.5};
73 
74 
75 class HcalTimingMonitorModule : public edm::EDAnalyzer {
76  public:
77  explicit HcalTimingMonitorModule(const edm::ParameterSet&);
78  ~HcalTimingMonitorModule();
79  void initialize();
80 
81  private:
82  virtual void beginJob() override ;
83  virtual void analyze(const edm::Event&, const edm::EventSetup&) override;
84  virtual void endJob() override ;
85 
86  double GetTime(double *data,int n){
87  int MaxI=-100; double Time=0,SumT=0,MaxT=-10;
88  for(int j=0;j<n;++j) if(MaxT<data[j]){ MaxT=data[j]; MaxI=j; }
89  if (MaxI>=0)
90  {
91  Time=MaxI*data[MaxI];
92  SumT=data[MaxI];
93  if(MaxI>0){ Time+=(MaxI-1)*data[MaxI-1]; SumT+=data[MaxI-1]; }
94  if(MaxI<(n-1)){ Time+=(MaxI+1)*data[MaxI+1]; SumT+=data[MaxI+1]; }
95  Time=Time/SumT;
96  }
97  return Time;
98  }
99  bool isSignal(double *data,int n){
100  int Imax=-1; double max=-100;
101  for(int i=0;i<n;i++) if(data[i]>max){max=data[i]; Imax=i;}
102  if(Imax==0 && Imax==(n-1)) return false;
103  float sum=data[Imax-1]+data[Imax+1];
104  if(data[Imax]>5.5 && sum>(data[Imax]*0.20)) return true;
105  return false;
106  }
109  // to be undependent from th DB we want to calculate pedestals using first 100 events
110  void set_hbhe(int eta,int phi,int depth,int cap,float val){
111  HBHE[eta+50][phi][depth][cap]+=val;
112  nHBHE[eta+50][phi][depth][cap]+=1.0;
113  }
114  void set_ho(int eta,int phi,int depth,int cap,float val){
115  HO[eta+50][phi][depth][cap]+=val;
116  nHO[eta+50][phi][depth][cap]+=1.0;
117  }
118  void set_hf(int eta,int phi,int depth,int cap,float val){
119  HF[eta+50][phi][depth][cap]+=val;
120  nHF[eta+50][phi][depth][cap]+=1.0;
121  }
122  double get_ped_hbhe(int eta,int phi,int depth,int cup){
123  if(nHBHE[eta+50][phi][depth][cup]<10) return 2.5;
124  if(nHBHE[eta+50][phi][depth][cup]!=0){
125  double ped=HBHE[eta+50][phi][depth][cup]/nHBHE[eta+50][phi][depth][cup];
126  if(ped>1.5 && ped<4.5) return ped;
127  }
128  return 99999;
129  }
130  double get_ped_ho(int eta,int phi,int depth,int cup){
131  if(nHO[eta+50][phi][depth][cup]<10) return 2.5;
132  if(nHO[eta+50][phi][depth][cup]!=0){
133  double ped=HO[eta+50][phi][depth][cup]/nHO[eta+50][phi][depth][cup];
134  if(ped>1.5 && ped<4.5) return ped;
135  }
136  return 99999;
137  }
138  double get_ped_hf(int eta,int phi,int depth,int cup){
139  if(nHF[eta+50][phi][depth][cup]<10) return 2.5;
140  if(nHF[eta+50][phi][depth][cup]!=0){
141  double ped=HF[eta+50][phi][depth][cup]/nHF[eta+50][phi][depth][cup];
142  if(ped>1.5 && ped<4.5) return ped;
143  }
144  return 99999;
145  }
146  double HBHE[100][73][5][4];
147  double nHBHE[100][73][5][4];
148  double HO[100][73][5][4];
149  double nHO[100][73][5][4];
150  double HF[100][73][5][4];
151  double nHF[100][73][5][4];
154  int counterEvt_;
155  int run_number;
156 
157  int TrigCSC,TrigDT,TrigRPC,TrigGCT;
158 
159  edm::ParameterSet parameters_;
160  DQMStore *dbe_;
161  std::string monitorName_;
162  int prescaleLS_,prescaleEvt_;
163  int GCTTriggerBit1_;
164  int GCTTriggerBit2_;
165  int GCTTriggerBit3_;
166  int GCTTriggerBit4_;
167  int GCTTriggerBit5_;
168  bool CosmicsCorr_;
169  bool Debug_;
170 
171  MonitorElement *HBEnergy,*HEEnergy,*HOEnergy,*HFEnergy;
172  MonitorElement *DTcand,*RPCbcand,*RPCfcand,*CSCcand,*OR;
173 
174  MonitorElement *HBShapeDT;
175  MonitorElement *HBShapeRPC;
176  MonitorElement *HBShapeGCT;
177  MonitorElement *HOShapeDT;
178  MonitorElement *HOShapeRPC;
179  MonitorElement *HOShapeGCT;
180  MonitorElement *HEShapeCSCp;
181  MonitorElement *HEShapeCSCm;
182  MonitorElement *HFShapeCSCp;
183  MonitorElement *HFShapeCSCm;
184 
185  MonitorElement *HBTimeDT;
186  MonitorElement *HBTimeRPC;
187  MonitorElement *HBTimeGCT;
188  MonitorElement *HOTimeDT;
189  MonitorElement *HOTimeRPC;
190  MonitorElement *HOTimeGCT;
191  MonitorElement *HETimeCSCp;
192  MonitorElement *HETimeCSCm;
193  MonitorElement *HFTimeCSCp;
194  MonitorElement *HFTimeCSCm;
195 
196  edm::EDGetTokenT<L1GlobalTriggerReadoutRecord> tok_gtro_;
197  edm::EDGetTokenT<L1MuGMTReadoutCollection> tok_L1mu_;
198 
199  edm::EDGetTokenT<HBHEDigiCollection> tok_hbhe_;
200  edm::EDGetTokenT<HODigiCollection> tok_ho_;
201  edm::EDGetTokenT<HFDigiCollection> tok_hf_;
202 };
203 
204 HcalTimingMonitorModule::HcalTimingMonitorModule(const edm::ParameterSet& iConfig) {
205 
206  tok_hbhe_ = consumes<HBHEDigiCollection>(iConfig.getParameter<edm::InputTag>("hbheDigiCollectionTag"));
207  tok_ho_ = consumes<HODigiCollection>(iConfig.getParameter<edm::InputTag>("hoDigiCollectionTag"));
208  tok_hf_ = consumes<HFDigiCollection>(iConfig.getParameter<edm::InputTag>("hfDigiCollectionTag"));
209 
210  std::string str;
211  parameters_ = iConfig;
212  dbe_ = edm::Service<DQMStore>().operator->();
213  // Base folder for the contents of this job
214  std::string subsystemname = parameters_.getUntrackedParameter<std::string>("subSystemFolder", "HcalTiming") ;
215 
216  monitorName_ = parameters_.getUntrackedParameter<std::string>("monitorName","HcalTiming");
217  if (monitorName_ != "" ) monitorName_ =subsystemname+"/"+monitorName_+"/" ;
218  counterEvt_=0;
219 
220  // some currently dummy things for compartability with GUI
221  dbe_->setCurrentFolder(subsystemname+"/EventInfo/");
222  str="reportSummary";
223  dbe_->bookFloat(str)->Fill(1); // Unknown status by default
224  str="reportSummaryMap";
225  MonitorElement* me=dbe_->book2D(str,str,5,0,5,1,0,1); // Unknown status by default
226  TH2F* myhist=me->getTH2F();
227  myhist->GetXaxis()->SetBinLabel(1,"HB");
228  myhist->GetXaxis()->SetBinLabel(2,"HE");
229  myhist->GetXaxis()->SetBinLabel(3,"HO");
230  myhist->GetXaxis()->SetBinLabel(4,"HF");
231  myhist->GetYaxis()->SetBinLabel(1,"Status");
232  // Unknown status by default
233  myhist->SetBinContent(1,1,-1);
234  myhist->SetBinContent(2,1,-1);
235  myhist->SetBinContent(3,1,-1);
236  myhist->SetBinContent(4,1,-1);
237  // Add ZDC at some point
238  myhist->GetXaxis()->SetBinLabel(5,"ZDC");
239  myhist->SetBinContent(5,1,-1); // no ZDC info known
240  myhist->SetOption("textcolz");
241 
242  run_number=0;
243  TrigCSC=TrigDT=TrigRPC=TrigGCT=0;
244  std::string sLabel = iConfig.getUntrackedParameter<std::string>("L1ADataLabel" , "l1GtUnpack");
245 
246  tok_gtro_ = consumes<L1GlobalTriggerReadoutRecord>(edm::InputTag(sLabel));
247  tok_L1mu_ = consumes<L1MuGMTReadoutCollection>(edm::InputTag(sLabel));
248 
249  prescaleLS_ = parameters_.getUntrackedParameter<int>("prescaleLS", 1);
250  prescaleEvt_ = parameters_.getUntrackedParameter<int>("prescaleEvt", 1);
251  GCTTriggerBit1_= parameters_.getUntrackedParameter<int>("GCTTriggerBit1", -1);
252  GCTTriggerBit2_= parameters_.getUntrackedParameter<int>("GCTTriggerBit2", -1);
253  GCTTriggerBit3_= parameters_.getUntrackedParameter<int>("GCTTriggerBit3", -1);
254  GCTTriggerBit4_= parameters_.getUntrackedParameter<int>("GCTTriggerBit4", -1);
255  GCTTriggerBit5_= parameters_.getUntrackedParameter<int>("GCTTriggerBit5", -1);
256  CosmicsCorr_ = parameters_.getUntrackedParameter<bool>("CosmicsCorr", true);
257  Debug_ = parameters_.getUntrackedParameter<bool>("Debug", true);
258  initialize();
259 }
260 
261 HcalTimingMonitorModule::~HcalTimingMonitorModule(){}
262 
263 // ------------ method called once each job just before starting event loop ------------
264 void HcalTimingMonitorModule::beginJob(){}
265 // ------------ method called once each job just after ending the event loop ------------
266 void HcalTimingMonitorModule::endJob(){}
267 
268 void HcalTimingMonitorModule::initialize(){
269  std::string str;
270  dbe_->setCurrentFolder(monitorName_+"DebugPlots");
271  str="L1MuGMTReadoutRecord_getDTBXCands"; DTcand =dbe_->book1D(str,str,5,-0.5,4.5);
272  str="L1MuGMTReadoutRecord_getBrlRPCCands"; RPCbcand=dbe_->book1D(str,str,5,-0.5,4.5);
273  str="L1MuGMTReadoutRecord_getFwdRPCCands"; RPCfcand=dbe_->book1D(str,str,5,-0.5,4.5);
274  str="L1MuGMTReadoutRecord_getCSCCands"; CSCcand =dbe_->book1D(str,str,5,-0.5,4.5);
275  str="DT_OR_RPCb_OR_RPCf_OR_CSC"; OR =dbe_->book1D(str,str,5,-0.5,4.5);
276 
277  str="HB Tower Energy (LinADC-PED)"; HBEnergy=dbe_->book1D(str,str,1000,-10,90);
278  str="HE Tower Energy (LinADC-PED)"; HEEnergy=dbe_->book1D(str,str,1000,-10,90);
279  str="HO Tower Energy (LinADC-PED)"; HOEnergy=dbe_->book1D(str,str,1000,-10,90);
280  str="HF Tower Energy (LinADC-PED)"; HFEnergy=dbe_->book1D(str,str,1000,-10,90);
281 
282  dbe_->setCurrentFolder(monitorName_+"ShapePlots");
283  str="HB Shape (DT Trigger)"; HBShapeDT =dbe_->book1D(str,str,10,-0.5,9.5);
284  str="HB Shape (RPC Trigger)"; HBShapeRPC =dbe_->book1D(str,str,10,-0.5,9.5);
285  str="HB Shape (GCT Trigger)"; HBShapeGCT =dbe_->book1D(str,str,10,-0.5,9.5);
286  str="HO Shape (DT Trigger)"; HOShapeDT =dbe_->book1D(str,str,10,-0.5,9.5);
287  str="HO Shape (RPC Trigger)"; HOShapeRPC =dbe_->book1D(str,str,10,-0.5,9.5);
288  str="HO Shape (GCT Trigger)"; HOShapeGCT =dbe_->book1D(str,str,10,-0.5,9.5);
289  str="HE+ Shape (CSC Trigger)"; HEShapeCSCp=dbe_->book1D(str,str,10,-0.5,9.5);
290  str="HE- Shape (CSC Trigger)"; HEShapeCSCm=dbe_->book1D(str,str,10,-0.5,9.5);
291  str="HF+ Shape (CSC Trigger)"; HFShapeCSCp=dbe_->book1D(str,str,10,-0.5,9.5);
292  str="HF- Shape (CSC Trigger)"; HFShapeCSCm=dbe_->book1D(str,str,10,-0.5,9.5);
293 
294  dbe_->setCurrentFolder(monitorName_+"TimingPlots");
295  str="HB Timing (DT Trigger)"; HBTimeDT =dbe_->book1D(str,str,100,0,10);
296  str="HB Timing (RPC Trigger)"; HBTimeRPC =dbe_->book1D(str,str,100,0,10);
297  str="HB Timing (GCT Trigger)"; HBTimeGCT =dbe_->book1D(str,str,100,0,10);
298  str="HO Timing (DT Trigger)"; HOTimeDT =dbe_->book1D(str,str,100,0,10);
299  str="HO Timing (RPC Trigger)"; HOTimeRPC =dbe_->book1D(str,str,100,0,10);
300  str="HO Timing (GCT Trigger)"; HOTimeGCT =dbe_->book1D(str,str,100,0,10);
301  str="HE+ Timing (CSC Trigger)"; HETimeCSCp =dbe_->book1D(str,str,100,0,10);
302  str="HE- Timing (CSC Trigger)"; HETimeCSCm =dbe_->book1D(str,str,100,0,10);
303  str="HF+ Timing (CSC Trigger)"; HFTimeCSCp =dbe_->book1D(str,str,100,0,10);
304  str="HF- Timing (CSC Trigger)"; HFTimeCSCm =dbe_->book1D(str,str,100,0,10);
305 }
306 
307 void HcalTimingMonitorModule::analyze(const edm::Event& iEvent, const edm::EventSetup& iSetup){
308 int HBcnt=0,HEcnt=0,HOcnt=0,HFcnt=0,eta,phi,depth,nTS;
309 int TRIGGER=0;
310  counterEvt_++;
311  if (prescaleEvt_<1) return;
312  if (counterEvt_%prescaleEvt_!=0) return;
313 
314  run_number=iEvent.id().run();
315  // Check GCT trigger bits
316  edm::Handle< L1GlobalTriggerReadoutRecord > gtRecord;
317 
318  if (!iEvent.getByToken( tok_gtro_, gtRecord))
319  return;
320  const TechnicalTriggerWord tWord = gtRecord->technicalTriggerWord();
321  const DecisionWord dWord = gtRecord->decisionWord();
322  //bool HFselfTrigger = tWord.at(9);
323  //bool HOselfTrigger = tWord.at(11);
324  bool GCTTrigger1 = dWord.at(GCTTriggerBit1_);
325  bool GCTTrigger2 = dWord.at(GCTTriggerBit2_);
326  bool GCTTrigger3 = dWord.at(GCTTriggerBit3_);
327  bool GCTTrigger4 = dWord.at(GCTTriggerBit4_);
328  bool GCTTrigger5 = dWord.at(GCTTriggerBit5_);
329  if(GCTTrigger1 || GCTTrigger2 || GCTTrigger3 || GCTTrigger4 || GCTTrigger5){ TrigGCT++; TRIGGER=+TRIG_GCT; }
330 
333  // define trigger trigger source (example from GMT group)
334  edm::Handle<L1MuGMTReadoutCollection> gmtrc_handle;
335  if (!iEvent.getByToken(tok_L1mu_,gmtrc_handle)) return;
336  L1MuGMTReadoutCollection const* gmtrc = gmtrc_handle.product();
337 
338  int idt =0;
339  int icsc =0;
340  int irpcb =0;
341  int irpcf =0;
342  int ndt[5] = {0,0,0,0,0};
343  int ncsc[5] = {0,0,0,0,0};
344  int nrpcb[5] = {0,0,0,0,0};
345  int nrpcf[5] = {0,0,0,0,0};
346  int N;
347 
348  std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
349  std::vector<L1MuGMTReadoutRecord>::const_iterator igmtrr;
350  N=0;
351  for(igmtrr=gmt_records.begin(); igmtrr!=gmt_records.end(); igmtrr++) {
352  std::vector<L1MuRegionalCand>::const_iterator iter1;
353  std::vector<L1MuRegionalCand> rmc;
354  // DTBX Trigger
355  rmc = igmtrr->getDTBXCands();
356  for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
357  if ( idt < MAXDTBX && !(*iter1).empty() ) {
358  idt++;
359  if(N<5) ndt[N]++;
360 
361  }
362  }
363  // CSC Trigger
364  rmc = igmtrr->getCSCCands();
365  for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
366  if ( icsc < MAXCSC && !(*iter1).empty() ) {
367  icsc++;
368  if(N<5) ncsc[N]++;
369  }
370  }
371  // RPCb Trigger
372  rmc = igmtrr->getBrlRPCCands();
373  for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
374  if ( irpcb < MAXRPC && !(*iter1).empty() ) {
375  irpcb++;
376  if(N<5) nrpcb[N]++;
377 
378  }
379  }
380  // RPCfwd Trigger
381  rmc = igmtrr->getFwdRPCCands();
382  for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
383  if ( irpcf < MAXRPC && !(*iter1).empty() ) {
384  irpcf++;
385  if(N<5) nrpcf[N]++;
386 
387  }
388  }
389 
390  N++;
391  }
392  if(ndt[0]) DTcand->Fill(0);
393  if(ndt[1]) DTcand->Fill(1);
394  if(ndt[2]) DTcand->Fill(2);
395  if(ndt[3]) DTcand->Fill(3);
396  if(ndt[4]) DTcand->Fill(4);
397  if(ncsc[0]) CSCcand->Fill(0);
398  if(ncsc[1]) CSCcand->Fill(1);
399  if(ncsc[2]) CSCcand->Fill(2);
400  if(ncsc[3]) CSCcand->Fill(3);
401  if(ncsc[4]) CSCcand->Fill(4);
402  if(nrpcb[0]) RPCbcand->Fill(0);
403  if(nrpcb[1]) RPCbcand->Fill(1);
404  if(nrpcb[2]) RPCbcand->Fill(2);
405  if(nrpcb[3]) RPCbcand->Fill(3);
406  if(nrpcb[4]) RPCbcand->Fill(4);
407  if(nrpcf[0]) RPCfcand->Fill(0);
408  if(nrpcf[1]) RPCfcand->Fill(1);
409  if(nrpcf[2]) RPCfcand->Fill(2);
410  if(nrpcf[3]) RPCfcand->Fill(3);
411  if(nrpcf[4]) RPCfcand->Fill(4);
412  if(ndt[0]||nrpcb[0]||nrpcf[0]||ncsc[0]) OR->Fill(0);
413  if(ndt[1]||nrpcb[1]||nrpcf[1]||ncsc[1]) OR->Fill(1);
414  if(ndt[2]||nrpcb[2]||nrpcf[2]||ncsc[2]) OR->Fill(2);
415  if(ndt[3]||nrpcb[3]||nrpcf[3]||ncsc[3]) OR->Fill(3);
416  if(ndt[4]||nrpcb[4]||nrpcf[4]||ncsc[4]) OR->Fill(4);
417 
418  if(ncsc[1]>0 ) { TrigCSC++; TRIGGER=+TRIG_CSC; }
419  if(ndt[1]>0 ) { TrigDT++; TRIGGER=+TRIG_DT; }
420  if(nrpcb[1]>0) { TrigRPC++; TRIGGER=+TRIG_RPC; }
421 
424  if(counterEvt_<100){
425  edm::Handle<HBHEDigiCollection> hbhe;
426  iEvent.getByToken(tok_hbhe_, hbhe);
427  if (hbhe.isValid())
428  {
429  for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
430  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
431  if(digi->id().subdet()==HcalBarrel) HBcnt++;
432  if(digi->id().subdet()==HcalEndcap) HEcnt++;
433  for(int i=0;i<nTS;i++)
434  if(digi->sample(i).adc()<20) set_hbhe(eta,phi,depth,digi->sample(i).capid(),adc2fC[digi->sample(i).adc()]);
435  }
436  }
437  edm::Handle<HODigiCollection> ho;
438  iEvent.getByToken(tok_ho_, ho);
439  if (ho.isValid())
440  {
441  for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
442  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
443  HOcnt++;
444  for(int i=0;i<nTS;i++)
445  if(digi->sample(i).adc()<20) set_ho(eta,phi,depth,digi->sample(i).capid(),adc2fC[digi->sample(i).adc()]);
446  }
447  } // if
448 
449  edm::Handle<HFDigiCollection> hf;
450  iEvent.getByToken(tok_hf_, hf);
451  if (hf.isValid())
452  {
453  for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
454  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
455  HFcnt++;
456  for(int i=0;i<nTS;i++)
457  if(digi->sample(i).adc()<20) set_hf(eta,phi,depth,digi->sample(i).capid(),adc2fC[digi->sample(i).adc()]);
458  }
459  }
460  } // if (counterEvt<100)
461  else{
463  double data[10];
464 
465  edm::Handle<HBHEDigiCollection> hbhe;
466  iEvent.getByToken(tok_hbhe_, hbhe);
467  if (hbhe.isValid())
468  {
469  for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
470  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
471  if(nTS>10) nTS=10;
472  if(digi->id().subdet()==HcalBarrel) HBcnt++;
473  if(digi->id().subdet()==HcalEndcap) HEcnt++;
474  double energy=0;
475  for(int i=0;i<nTS;i++){
476  data[i]=adc2fC[digi->sample(i).adc()]-get_ped_hbhe(eta,phi,depth,digi->sample(i).capid());
477  energy+=data[i];
478  }
479  if(digi->id().subdet()==HcalBarrel) HBEnergy->Fill(energy);
480  if(digi->id().subdet()==HcalEndcap) HEEnergy->Fill(energy);
481  if(!isSignal(data,nTS)) continue;
482  for(int i=0;i<nTS;i++){
483  if(data[i]>-1.0){
484  if(digi->id().subdet()==HcalBarrel && (TRIGGER|TRIG_DT)==TRIG_DT) HBShapeDT->Fill(i,data[i]);
485  if(digi->id().subdet()==HcalBarrel && (TRIGGER|TRIG_RPC)==TRIG_RPC) HBShapeRPC->Fill(i,data[i]);
486  if(digi->id().subdet()==HcalBarrel && (TRIGGER|TRIG_GCT)==TRIG_GCT) HBShapeGCT->Fill(i,data[i]);
487  if(digi->id().subdet()==HcalEndcap && (TRIGGER|TRIG_CSC)==TRIG_CSC && eta>0) HEShapeCSCp->Fill(i,data[i]);
488  if(digi->id().subdet()==HcalEndcap && (TRIGGER|TRIG_CSC)==TRIG_CSC && eta<0) HEShapeCSCm->Fill(i,data[i]);
489  }
490  }
491  double Time=GetTime(data,nTS);
492  if(digi->id().subdet()==HcalBarrel){
493  if(CosmicsCorr_) Time+=(7.5*sin((phi*5.0)/180.0*3.14159))/25.0;
494  if((TRIGGER&TRIG_DT)==TRIG_DT) HBTimeDT ->Fill(GetTime(data,nTS));
495  if((TRIGGER&TRIG_RPC)==TRIG_RPC) HBTimeRPC->Fill(GetTime(data,nTS));
496  if((TRIGGER&TRIG_GCT)==TRIG_GCT) HBTimeGCT->Fill(GetTime(data,nTS));
497  }else{
498  if(CosmicsCorr_) Time+=(3.5*sin((phi*5.0)/180.0*3.14159))/25.0;
499  if(digi->id().subdet()==HcalEndcap && (TRIGGER&TRIG_CSC)==TRIG_CSC && eta>0) HETimeCSCp->Fill(Time);
500  if(digi->id().subdet()==HcalEndcap && (TRIGGER&TRIG_CSC)==TRIG_CSC && eta<0) HETimeCSCm->Fill(Time);
501  }
502  }
503  } // if (...)
504 
505  edm::Handle<HODigiCollection> ho;
506  iEvent.getByToken(tok_ho_, ho);
507  if (ho.isValid())
508  {
509  for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
510  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
511  if(nTS>10) nTS=10;
512  HOcnt++;
513  double energy=0;
514  for(int i=0;i<nTS;i++){
515  data[i]=adc2fC[digi->sample(i).adc()]-get_ped_ho(eta,phi,depth,digi->sample(i).capid());
516  energy+=data[i];
517  }
518  HOEnergy->Fill(energy);
519  if(!isSignal(data,nTS)) continue;
520  for(int i=0;i<nTS;i++){
521  if(data[i]>-1.0){
522  if((TRIGGER&TRIG_DT)==TRIG_DT) HOShapeDT->Fill(i,data[i]);
523  if((TRIGGER&TRIG_RPC)==TRIG_RPC) HOShapeRPC->Fill(i,data[i]);
524  if((TRIGGER&TRIG_GCT)==TRIG_GCT) HOShapeGCT->Fill(i,data[i]);
525  }
526  }
527  double Time=GetTime(data,nTS);
528  if(CosmicsCorr_) Time+=(12.0*sin((phi*5.0)/180.0*3.14159))/25.0;
529  if((TRIGGER&TRIG_DT)==TRIG_DT) HOTimeDT->Fill(Time);
530  if((TRIGGER&TRIG_RPC)==TRIG_RPC) HOTimeRPC->Fill(Time);
531  if((TRIGGER&TRIG_GCT)==TRIG_GCT) HOTimeGCT->Fill(Time);
532  }
533  }// if (ho)
534 
535  edm::Handle<HFDigiCollection> hf;
536  iEvent.getByToken(tok_hf_, hf);
537  if (hf.isValid())
538  {
539  for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
540  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
541  if(nTS>10) nTS=10;
542  HFcnt++;
543  double energy=0;
544  for(int i=0;i<nTS;i++){
545  data[i]=adc2fC[digi->sample(i).adc()]-get_ped_hf(eta,phi,depth,digi->sample(i).capid());
546  energy+=data[i];
547  }
548  HFEnergy->Fill(energy);
549  if(energy<15.0) continue;
550  for(int i=0;i<nTS;i++){
551  if(data[i]>-1.0){
552  if((TRIGGER&TRIG_CSC)==TRIG_CSC && eta>0) HFShapeCSCp->Fill(i,data[i]);
553  if((TRIGGER&TRIG_CSC)==TRIG_CSC && eta<0) HFShapeCSCm->Fill(i,data[i]);
554  }
555  }
556  if((TRIGGER&TRIG_CSC)==TRIG_CSC && eta>0) HFTimeCSCp->Fill(GetTime(data,nTS));
557  if((TRIGGER&TRIG_CSC)==TRIG_CSC && eta<0) HFTimeCSCm->Fill(GetTime(data,nTS));
558  }
559  } // if (hf)
562  }
563  if(Debug_) if((counterEvt_%100)==0) printf("Run: %i,Events processed: %i (HB: %i towers,HE: %i towers,HO: %i towers,HF: %i towers)"
564  " CSC: %i DT: %i RPC: %i GCT: %i\n",
565  run_number,counterEvt_,HBcnt,HEcnt,HOcnt,HFcnt,TrigCSC,TrigDT,TrigRPC,TrigGCT);
568 }
569 //define this as a plug-in
570 DEFINE_FWK_MODULE(HcalTimingMonitorModule);
571 
572 
edm::EventID::run
RunNumber_t run() const
Definition: EventID.h:42
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
MAXCSC
static const int MAXCSC
Definition: HcalTimingMonitorModule.cc:21
TRIG_RPC
static const int TRIG_RPC
Definition: HcalTimingMonitorModule.cc:23
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:954
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:434
HcalTimingMonitorModule::HcalTimingMonitorModule
HcalTimingMonitorModule(const edm::ParameterSet &)
Definition: HcalTimingMonitorModule.cc:204
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
HcalTimingMonitorModule::tok_ho_
edm::EDGetTokenT< HODigiCollection > tok_ho_
Definition: HcalTimingMonitorModule.cc:200
HcalTimingMonitorModule::beginJob
virtual void beginJob() override
Definition: HcalTimingMonitorModule.cc:264
funct::sin
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
HcalTimingMonitorModule::GetTime
double GetTime(double *data, int n)
Definition: HcalTimingMonitorModule.cc:86
HcalTimingMonitorModule::get_ped_ho
double get_ped_ho(int eta, int phi, int depth, int cup)
Definition: HcalTimingMonitorModule.cc:130
edm::SortedCollection< HBHEDataFrame >::const_iterator
std::vector< HBHEDataFrame >::const_iterator const_iterator
Definition: SortedCollection.h:82
TRIG_GCT
static const int TRIG_GCT
Definition: HcalTimingMonitorModule.cc:24
MAXRPC
static const int MAXRPC
Definition: HcalTimingMonitorModule.cc:19
HcalTimingMonitorModule::set_hf
void set_hf(int eta, int phi, int depth, int cap, float val)
Definition: HcalTimingMonitorModule.cc:118
HcalTimingMonitorModule::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: HcalTimingMonitorModule.cc:307
eta
T eta() const
Definition: Basic3DVectorLD.h:177
adc2fC
static const float adc2fC[128]
Definition: HcalTimingMonitorModule.cc:63
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
DQMStore::bookFloat
MonitorElement * bookFloat(const char *name)
Book float.
Definition: DQMStore.cc:891
HcalTimingMonitorModule::isSignal
bool isSignal(double *data, int n)
Definition: HcalTimingMonitorModule.cc:99
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:166
TRIG_CSC
static const int TRIG_CSC
Definition: HcalTimingMonitorModule.cc:25
iEvent
int iEvent
Definition: GenABIO.cc:230
MAXDTBX
static const int MAXDTBX
Definition: HcalTimingMonitorModule.cc:20
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
DecisionWord
std::vector< bool > DecisionWord
typedefs
Definition: L1GlobalTriggerReadoutSetupFwd.h:34
HcalTimingMonitorModule::tok_hf_
edm::EDGetTokenT< HFDigiCollection > tok_hf_
Definition: HcalTimingMonitorModule.cc:201
j
int j
Definition: DBlmapReader.cc:9
TechnicalTriggerWord
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
Definition: L1GlobalTriggerReadoutSetupFwd.h:41
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:76
HcalTimingMonitorModule::get_ped_hbhe
double get_ped_hbhe(int eta, int phi, int depth, int cup)
Definition: HcalTimingMonitorModule.cc:122
HcalTimingMonitorModule::set_hbhe
void set_hbhe(int eta, int phi, int depth, int cap, float val)
Definition: HcalTimingMonitorModule.cc:110
N
#define N
Definition: blowfish.cc:9
HcalTimingMonitorModule::set_ho
void set_ho(int eta, int phi, int depth, int cap, float val)
Definition: HcalTimingMonitorModule.cc:114
HcalTimingMonitorModule::tok_gtro_
edm::EDGetTokenT< L1GlobalTriggerReadoutRecord > tok_gtro_
Definition: HcalTimingMonitorModule.cc:196
edm::Handle::product
T const * product() const
Definition: Handle.h:81
edm::EventBase::id
edm::EventID id() const
Definition: EventBase.h:56
data
char data[epos_bytes_allocation]
Definition: EPOS_Wrapper.h:82
HcalTimingMonitorModule::tok_L1mu_
edm::EDGetTokenT< L1MuGMTReadoutCollection > tok_L1mu_
Definition: HcalTimingMonitorModule.cc:197
L1MuGMTReadoutCollection::getRecords
std::vector< L1MuGMTReadoutRecord > const & getRecords() const
Definition: L1MuGMTReadoutCollection.h:51
relval_parameters_module.energy
string energy
Definition: relval_parameters_module.py:29
HcalTimingMonitorModule::endJob
virtual void endJob() override
Definition: HcalTimingMonitorModule.cc:266
HcalTimingMonitorModule::tok_hbhe_
edm::EDGetTokenT< HBHEDigiCollection > tok_hbhe_
Definition: HcalTimingMonitorModule.cc:199
MonitorElement::getTH2F
TH2F * getTH2F(void) const
Definition: MonitorElement.cc:1563
TRIG_DT
static const int TRIG_DT
Definition: HcalTimingMonitorModule.cc:22
DQMStore::book2D
MonitorElement * book2D(const char *name, const char *title, int nchX, double lowX, double highX, int nchY, double lowY, double highY)
Book 2D histogram.
Definition: DQMStore.cc:1082
HcalTimingMonitorModule::get_ped_hf
double get_ped_hf(int eta, int phi, int depth, int cup)
Definition: HcalTimingMonitorModule.cc:138
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:667
phi
Definition: DDAxes.h:10
Generated for CMSSW Reference Manual by doxygen 1.8.5