CMS 3D CMS Logo

 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
HcalDetDiagTimingMonitor.cc
Go to the documentation of this file.
2 
8 
9 #include <math.h>
10 
11 // this is to retrieve HCAL digi's
13 // this is to retrieve GT digi's
19 
20 
21 
23 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,
24  13.5,15.,17.,19.,21.,23.,25.,27.,29.5,32.5,35.5,38.5,42.,46.,50.,54.5,59.5,
25  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,
26  124.5,129.5,137.,147.,157.,167.,177.,187.,197.,209.5,224.5,239.5,254.5,272.,
27  292.,312.,334.5,359.5,384.5,359.5,384.5,409.5,434.5,459.5,484.5,509.5,534.5,
28  559.5,584.5,609.5,634.5,659.5,684.5,709.5,747.,797.,847.,897.,947.,997.,
29  1047.,1109.5,1184.5,1259.5,1334.5,1422.,1522.,1622.,1734.5,1859.5,1984.5,
30  1859.5,1984.5,2109.5,2234.5,2359.5,2484.5,2609.5,2734.5,2859.5,2984.5,
31  3109.5,3234.5,3359.5,3484.5,3609.5,3797.,4047.,4297.,4547.,4797.,5047.,
32  5297.,5609.5,5984.5,6359.5,6734.5,7172.,7672.,8172.,8734.5,9359.5,9984.5};
34 static const int MAXGEN =10;
35 static const int MAXRPC =20;
36 static const int MAXDTBX=20;
37 static const int MAXCSC =20;
38 static const int MAXGMT =20;
39 static const int TRIG_DT =1;
40 static const int TRIG_RPC=2;
41 static const int TRIG_GCT=4;
42 static const int TRIG_CSC=8;
43 static const int TRIG_RPCF=16;
44 
46 {
47  Online_ = ps.getUntrackedParameter<bool>("online",false);
48  mergeRuns_ = ps.getUntrackedParameter<bool>("mergeRuns",false);
49  enableCleanup_ = ps.getUntrackedParameter<bool>("enableCleanup",false);
50  debug_ = ps.getUntrackedParameter<int>("debug",0);
51  prefixME_ = ps.getUntrackedParameter<std::string>("subSystemFolder","Hcal/");
52  if (prefixME_.substr(prefixME_.size()-1,prefixME_.size())!="/")
53  prefixME_.append("/");
54  subdir_ = ps.getUntrackedParameter<std::string>("TaskFolder","DetDiagTimingMonitor_Hcal");
55  if (subdir_.size()>0 && subdir_.substr(subdir_.size()-1,subdir_.size())!="/")
56  subdir_.append("/");
57  subdir_=prefixME_+subdir_;
58  AllowedCalibTypes_ = ps.getUntrackedParameter<std::vector<int> > ("AllowedCalibTypes");
59  skipOutOfOrderLS_ = ps.getUntrackedParameter<bool>("skipOutOfOrderLS","false");
60  NLumiBlocks_ = ps.getUntrackedParameter<int>("NLumiBlocks",4000);
61  makeDiagnostics_ = ps.getUntrackedParameter<bool>("makeDiagnostics",false);
62 
63  GCTTriggerBit1_= ps.getUntrackedParameter<int>("GCTTriggerBit1", 15);
64  GCTTriggerBit2_= ps.getUntrackedParameter<int>("GCTTriggerBit2", 16);
65  GCTTriggerBit3_= ps.getUntrackedParameter<int>("GCTTriggerBit3", 17);
66  GCTTriggerBit4_= ps.getUntrackedParameter<int>("GCTTriggerBit4", 18);
67  GCTTriggerBit5_= ps.getUntrackedParameter<int>("GCTTriggerBit5", 16);
68  CosmicsCorr_ = ps.getUntrackedParameter<bool>("CosmicsCorr", true);
69 
72  FEDRawDataCollection_ = ps.getUntrackedParameter<edm::InputTag>("FEDRawDataCollection",edm::InputTag("source",""));
73 }
74 
76 
78  if(dbe_){
81  dbe_ = 0;
82  }
83 }
85 
87 {
88  if (debug_>1) std::cout <<"HcalDetDiagTimingMonitor::beginRun"<<std::endl;
90 
91  if (tevt_==0) this->setup(); // set up histograms if they have not been created before
92  if (mergeRuns_==false)
93  this->reset();
94 
95  return;
96 
97 } // void HcalNDetDiagTimingMonitor::beginRun(...)
98 
99 
100 
102 {
103 
105 
106  std::string str;
107  if(dbe_!=NULL){
109  str="Hcal Timing summary"; Summary = dbe_->book2D(str,str,6,0,6,6,0,6);
110  Summary->setBinLabel(1,"DT",1);
111  Summary->setBinLabel(2,"RPC",1);
112  Summary->setBinLabel(3,"GCT",1);
113  Summary->setBinLabel(4,"CSC",1);
114  Summary->setBinLabel(5,"RPCf",1);
115  Summary->setBinLabel(6,"bit11",1);
116  Summary->setBinLabel(1,"HB",2);
117  Summary->setBinLabel(2,"HO",2);
118  Summary->setBinLabel(3,"HEM",2);
119  Summary->setBinLabel(4,"HEP",2);
120  Summary->setBinLabel(5,"HFM",2);
121  Summary->setBinLabel(6,"HFP",2);
122  for(int i=1;i<=6;i++) for(int j=1;j<=6;j++) Summary->setBinContent(i,j,-1);
123 
124  dbe_->setCurrentFolder(subdir_+"Timing Plots");
125  str="HB Timing (DT Trigger)"; HBTimeDT = dbe_->book1D(str,str,100,0,10);
126  str="HO Timing (DT Trigger)"; HOTimeDT = dbe_->book1D(str,str,100,0,10);
127  str="HB Timing (RPC Trigger)"; HBTimeRPC = dbe_->book1D(str,str,100,0,10);
128  str="HO Timing (RPC Trigger)"; HOTimeRPC = dbe_->book1D(str,str,100,0,10);
129  str="HB Timing (HO SelfTrigger tech bit 11)"; HBTimeHO = dbe_->book1D(str,str,100,0,10);
130  str="HO Timing (HO SelfTrigger tech bit 11)"; HOTimeHO = dbe_->book1D(str,str,100,0,10);
131 
132  str="HB Timing (GCT Trigger alg bit 15 16 17 18)"; HBTimeGCT =dbe_->book1D(str,str,100,0,10);
133  str="HB Timing (GCT Trigger alg bit 15 16 17 18)"; HOTimeGCT =dbe_->book1D(str,str,100,0,10);
134 
135  str="HEP Timing (CSC Trigger)"; HETimeCSCp =dbe_->book1D(str,str,100,0,10);
136  str="HEM Timing (CSC Trigger)"; HETimeCSCm =dbe_->book1D(str,str,100,0,10);
137  str="HEP Timing (RPCf Trigger)"; HETimeRPCp =dbe_->book1D(str,str,100,0,10);
138  str="HEM Timing (RPCf Trigger)"; HETimeRPCm =dbe_->book1D(str,str,100,0,10);
139  str="HFP Timing (CSC Trigger)"; HFTimeCSCp =dbe_->book1D(str,str,100,0,10);
140  str="HFM Timing (CSC Trigger)"; HFTimeCSCm =dbe_->book1D(str,str,100,0,10);
141  }
142 }
143 
145 {
146  if (!IsAllowedCalibType()) return;
147  if (LumiInOrder(iEvent.luminosityBlock())==false) return;
148  HcalBaseDQMonitor::analyze(iEvent, iSetup);
149 
150 
151  int eta,phi,depth,nTS,BXinEVENT=1,TRIGGER=0;
152 
153  if(!dbe_) return;
154  // We do not want to look at Abort Gap events
156  iEvent.getByLabel(FEDRawDataCollection_,rawdata);
157  //checking FEDs for calibration information
158  if(!rawdata.isValid()) return;
160  const FEDRawData& fedData = rawdata->FEDData(i) ;
161  if ( fedData.size() < 24 ) continue ;
162  if(((const HcalDCCHeader*)(fedData.data()))->getCalibType()!=hc_Null) return;
163  }
164 
167  bool GCTTrigger1=false,GCTTrigger2=false,GCTTrigger3=false,GCTTrigger4=false,GCTTrigger5=false,HOselfTrigger=false;
168  // Check GCT trigger bits
170  iEvent.getByLabel(L1ADataLabel_, gtRecord);
171  if(gtRecord.isValid()){
172  const TechnicalTriggerWord tWord = gtRecord->technicalTriggerWord();
173  const DecisionWord dWord = gtRecord->decisionWord();
174  //bool HFselfTrigger = tWord.at(9);
175  HOselfTrigger = tWord.at(11);
176 
177  GCTTrigger1 = dWord.at(GCTTriggerBit1_);
178  GCTTrigger2 = dWord.at(GCTTriggerBit2_);
179  GCTTrigger3 = dWord.at(GCTTriggerBit3_);
180  GCTTrigger4 = dWord.at(GCTTriggerBit4_);
181  GCTTrigger5 = dWord.at(GCTTriggerBit5_);
182 
183  // define trigger trigger source (example from GMT group)
185  iEvent.getByLabel(L1ADataLabel_,gmtrc_handle);
186  if(!gmtrc_handle.isValid()) return;
187  L1MuGMTReadoutCollection const* gmtrc = gmtrc_handle.product();
188 
189  int idt =0;
190  int icsc =0;
191  int irpcb =0;
192  int irpcf =0;
193  int ndt[5] = {0,0,0,0,0};
194  int ncsc[5] = {0,0,0,0,0};
195  int nrpcb[5] = {0,0,0,0,0};
196  int nrpcf[5] = {0,0,0,0,0};
197  int N;
198  std::vector<L1MuGMTReadoutRecord> gmt_records = gmtrc->getRecords();
199  std::vector<L1MuGMTReadoutRecord>::const_iterator igmtrr;
200  N=0;
201  int NN=0;
202  for(igmtrr=gmt_records.begin(); igmtrr!=gmt_records.end(); igmtrr++) {
203  if(igmtrr->getBxInEvent()==0) BXinEVENT=NN;
204  NN++;
205  std::vector<L1MuRegionalCand>::const_iterator iter1;
206  std::vector<L1MuRegionalCand> rmc;
207  // DTBX Trigger
208  rmc = igmtrr->getDTBXCands();
209  for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
210  if ( idt < MAXDTBX && !(*iter1).empty() ) {
211  idt++;
212  if(N<5) ndt[N]++;
213 
214  }
215  }
216  // CSC Trigger
217  rmc = igmtrr->getCSCCands();
218  for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
219  if ( icsc < MAXCSC && !(*iter1).empty() ) {
220  icsc++;
221  if(N<5) ncsc[N]++;
222  }
223  }
224  // RPCb Trigger
225  rmc = igmtrr->getBrlRPCCands();
226  for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
227  if ( irpcb < MAXRPC && !(*iter1).empty() ) {
228  irpcb++;
229  if(N<5) nrpcb[N]++;
230 
231  }
232  }
233  // RPCfwd Trigger
234  rmc = igmtrr->getFwdRPCCands();
235  for(iter1=rmc.begin(); iter1!=rmc.end(); iter1++) {
236  if ( irpcf < MAXRPC && !(*iter1).empty() ) {
237  irpcf++;
238  if(N<5) nrpcf[N]++;
239 
240  }
241  }
242 
243  N++;
244  }
245  if(ncsc[BXinEVENT]>0 ) { TRIGGER=+TRIG_CSC; }
246  if(ndt[BXinEVENT]>0 ) { TRIGGER=+TRIG_DT; }
247  if(nrpcb[BXinEVENT]>0) { TRIGGER=+TRIG_RPC; }
248  if(nrpcf[BXinEVENT]>0) { TRIGGER=+TRIG_RPCF; }
249  if(GCTTrigger1 || GCTTrigger2 || GCTTrigger3 || GCTTrigger4 || GCTTrigger5) { TRIGGER=+TRIG_GCT; }
250  }
253  if(ievt_<100){
255  iEvent.getByLabel(inputLabelDigi_,hbhe);
256  if(hbhe.isValid()){
257  for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
258  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
259  for(int i=0;i<nTS;i++) if(digi->sample(i).adc()<20) set_hbhe(eta,phi,depth,digi->sample(i).capid(),adc2fC[digi->sample(i).adc()]);
260  }
261  }
263  iEvent.getByLabel(inputLabelDigi_,ho);
264  if(ho.isValid()){
265  for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
266  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
267  for(int i=0;i<nTS;i++) if(digi->sample(i).adc()<20) set_ho(eta,phi,depth,digi->sample(i).capid(),adc2fC[digi->sample(i).adc()]);
268  }
269  }
271  iEvent.getByLabel(inputLabelDigi_,hf);
272  if(hf.isValid()){
273  for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
274  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
275  for(int i=0;i<nTS;i++) if(digi->sample(i).adc()<20) set_hf(eta,phi,depth,digi->sample(i).capid(),adc2fC[digi->sample(i).adc()]);
276  }
277  }
278  return;
279  }
282  double data[20];
284  iEvent.getByLabel(inputLabelDigi_,hbhe);
285  if(hbhe.isValid()){
286  for(HBHEDigiCollection::const_iterator digi=hbhe->begin();digi!=hbhe->end();digi++){
287  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
288  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()]-get_ped_hbhe(eta,phi,depth,digi->sample(i).capid());
289  if(!isSignal(data,nTS)) continue;
290 
291  occHBHE[eta+50][phi][depth]+=1.0; occSum+=1.0;
292  if((occHBHE[eta+50][phi][depth]/(double)(ievt_))>0.001) continue;
293 
294  double Time=GetTime(data,nTS);
295  if(digi->id().subdet()==HcalBarrel){
296  if(TRIGGER==TRIG_GCT) HBTimeGCT->Fill(Time);
297  if(CosmicsCorr_) Time+=(7.5*sin((phi*5.0)/180.0*3.14159))/25.0;
298  if(TRIGGER==TRIG_DT) HBTimeDT->Fill(Time);
299  if(HOselfTrigger) HBTimeHO->Fill(Time);
300  if(TRIGGER==TRIG_RPC) HBTimeRPC->Fill(Time);
301  }
302  if(digi->id().subdet()==HcalEndcap){
303  if(CosmicsCorr_) Time+=(3.5*sin((phi*5.0)/180.0*3.14159))/25.0;
304  if(TRIGGER==TRIG_CSC && eta>0) HETimeCSCp->Fill(Time);
305  if(TRIGGER==TRIG_CSC && eta<0) HETimeCSCm->Fill(Time);
306  if(TRIGGER==TRIG_RPCF && eta>0) HETimeRPCp->Fill(Time);
307  if(TRIGGER==TRIG_RPCF && eta<0) HETimeRPCm->Fill(Time);
308  }
309 
310  }
311  }
313  iEvent.getByLabel(inputLabelDigi_,ho);
314  if(ho.isValid()){
315  for(HODigiCollection::const_iterator digi=ho->begin();digi!=ho->end();digi++){
316  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
317  for(int i=0;i<nTS;i++) data[i]=adc2fC[digi->sample(i).adc()]-get_ped_ho(eta,phi,depth,digi->sample(i).capid());
318  if(!isSignal(data,nTS)) continue;
319  occHO[eta+50][phi][depth]+=1.0;
320  occSum+=1.0;
321  if((occHO[eta+50][phi][depth]/(double)(ievt_))>0.001) continue;
322 
323  double Time=GetTime(data,nTS);
324  if(CosmicsCorr_) Time+=(12.0*sin((phi*5.0)/180.0*3.14159))/25.0;
325  if(TRIGGER==TRIG_DT) HOTimeDT ->Fill(Time);
326  if(HOselfTrigger) HOTimeHO ->Fill(Time);
327  if(TRIGGER==TRIG_RPC) HOTimeRPC->Fill(Time);
328  if(TRIGGER==TRIG_GCT) HOTimeGCT->Fill(Time);
329  }
330  }
332  iEvent.getByLabel(inputLabelDigi_,hf);
333  if(hf.isValid()){
334  for(HFDigiCollection::const_iterator digi=hf->begin();digi!=hf->end();digi++){
335  eta=digi->id().ieta(); phi=digi->id().iphi(); depth=digi->id().depth(); nTS=digi->size();
336  double energy=0;
337  for(int i=0;i<nTS;i++){
338  data[i]=adc2fC[digi->sample(i).adc()]-get_ped_hf(eta,phi,depth,digi->sample(i).capid());
339  energy+=data[i];
340  }
341 
342  if(energy<25.0) continue;
343  occHF[eta+50][phi][depth]+=1.0;
344  occSum+=1.0;
345 
346  double Time=GetTime(data,nTS);
347  if((occHF[eta+50][phi][depth]/(double)(ievt_))>0.01) continue;
348 
349  if(TRIGGER==TRIG_CSC && eta>0) HFTimeCSCp->Fill(Time);
350  if(TRIGGER==TRIG_CSC && eta<0) HFTimeCSCm->Fill(Time);
351  }
352  }
353  if((ievt_%500)==0){
354  CheckTiming();
355  //printf("Run: %i, Evants processed: %i\n",iEvent.run(),ievt_);
356  }
357 }
358 
360  if(HBTimeDT->getEntries()>10){
361  Summary->setBinContent(1,1,1);
362  }
363  if(HBTimeRPC->getEntries()>10){
364  Summary->setBinContent(2,1,1);
365  }
366  if(HBTimeGCT->getEntries()>10){
367  Summary->setBinContent(3,1,1);
368  }
369  if(HBTimeHO->getEntries()>10){
370  Summary->setBinContent(6,1,1);
371  }
372  if(HOTimeDT->getEntries()>10){
373  Summary->setBinContent(1,2,1);
374  }
375  if(HOTimeRPC->getEntries()>10){
376  Summary->setBinContent(2,2,1);
377  }
378  if(HOTimeGCT->getEntries()>10){
379  Summary->setBinContent(3,2,1);
380  }
381  if(HOTimeHO->getEntries()>10){
382  Summary->setBinContent(6,2,1);
383  }
384  if(HETimeCSCp->getEntries()>10){
385  Summary->setBinContent(4,4,1);
386  }
387  if(HETimeCSCm->getEntries()>10){
388  Summary->setBinContent(4,3,1);
389  }
390  if(HETimeRPCp->getEntries()>10){
391  Summary->setBinContent(5,4,1);
392  }
393  if(HETimeRPCm->getEntries()>10){
394  Summary->setBinContent(5,3,1);
395  }
396  if(HFTimeCSCp->getEntries()>10){
397  Summary->setBinContent(4,6,1);
398  }
399  if(HFTimeCSCm->getEntries()>10){
400  Summary->setBinContent(4,4,1);
401  }
402 }
403 
T getUntrackedParameter(std::string const &, T const &) const
int i
Definition: DBlmapReader.cc:9
bool LumiInOrder(int lumisec)
static const int MAXCSC
void setBinContent(int binx, double content)
set content of bin (1-D)
virtual void analyze(const edm::Event &e, const edm::EventSetup &c)
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:519
void set_hf(int eta, int phi, int depth, int cap, float val)
static const int TRIG_CSC
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
std::vector< int > AllowedCalibTypes_
Sin< T >::type sin(const T &t)
Definition: Sin.h:22
std::vector< T >::const_iterator const_iterator
void setBinLabel(int bin, const std::string &label, int axis=1)
set bin label for x, y or z axis (axis=1, 2, 3 respectively)
edm::LuminosityBlockNumber_t luminosityBlock() const
Definition: EventBase.h:59
#define NULL
Definition: scimark2.h:8
static const int MAXRPC
void analyze(const edm::Event &iEvent, const edm::EventSetup &iSetup)
double getEntries(void) const
get # of entries
static const int MAXGEN
T eta() const
size_t size() const
Lenght of the data buffer in bytes.
Definition: FEDRawData.h:47
void Fill(long long x)
static const int TRIG_RPC
int iEvent
Definition: GenABIO.cc:243
virtual void beginRun(const edm::Run &run, const edm::EventSetup &c)
void removeContents(void)
erase all monitoring elements in current directory (not including subfolders);
Definition: DQMStore.cc:2330
static const int MAXDTBX
bool isSignal(double *data, int n)
double get_ped_ho(int eta, int phi, int depth, int cup)
std::vector< bool > DecisionWord
typedefs
static const int TRIG_GCT
int j
Definition: DBlmapReader.cc:9
std::vector< bool > TechnicalTriggerWord
technical trigger bits (64 bits)
bool isValid() const
Definition: HandleBase.h:76
bool getByLabel(InputTag const &tag, Handle< PROD > &result) const
Definition: Event.h:359
void set_ho(int eta, int phi, int depth, int cap, float val)
std::vector< L1MuGMTReadoutRecord > getRecords() const
static const int TRIG_RPCF
static const float adc2fC[128]
HcalDetDiagTimingMonitor(const edm::ParameterSet &ps)
double get_ped_hbhe(int eta, int phi, int depth, int cup)
double get_ped_hf(int eta, int phi, int depth, int cup)
T const * product() const
Definition: Handle.h:74
void beginRun(const edm::Run &run, const edm::EventSetup &c)
void set_hbhe(int eta, int phi, int depth, int cap, float val)
const unsigned char * data() const
Return a const pointer to the beginning of the data buffer.
Definition: FEDRawData.cc:29
static const int MAXGMT
tuple cout
Definition: gather_cfg.py:41
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:647
virtual void setup(void)
double GetTime(double *data, int n)
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:237
Definition: Run.h:31
static const int TRIG_DT
Definition: DDAxes.h:10