test
CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
DTLocalTriggerBaseTask.cc
Go to the documentation of this file.
1 /*
2  * \file DTLocalTriggerBaseTask.cc
3  *
4  * \author C. Battilana - CIEMAT
5  *
6 */
7 
8 #include "DQM/DTMonitorModule/src/DTLocalTriggerBaseTask.h"
9 
10 // Framework
11 #include "FWCore/Framework/interface/EventSetup.h"
12 
13 // DT DQM
14 #include "DQM/DTMonitorModule/interface/DTTimeEvolutionHisto.h"
15 
16 // DT trigger
17 #include "DQM/DTMonitorModule/interface/DTTrigGeomUtils.h"
18 
19 // Geometry
20 #include "DataFormats/GeometryVector/interface/Pi.h"
21 #include "Geometry/Records/interface/MuonGeometryRecord.h"
22 #include "Geometry/DTGeometry/interface/DTGeometry.h"
23 #include "Geometry/DTGeometry/interface/DTLayer.h"
24 #include "Geometry/DTGeometry/interface/DTTopology.h"
25 
26 //Root
27 #include"TH1.h"
28 #include"TAxis.h"
29 
30 #include <sstream>
31 #include <iostream>
32 #include <fstream>
33 
34 using namespace edm;
35 using namespace std;
36 
37 class DTTPGCompareUnit {
38 
39 public:
40 
41  DTTPGCompareUnit() { theQual[0]=-1 ; theQual[1]=-1; }
42  ~DTTPGCompareUnit() { };
43 
44  void setDDU(int qual, int bx) { theQual[0] = qual; theBX[0] = bx; }
45  void setTM(int qual, int bx) { theQual[1] = qual; theBX[1] = bx; }
46 
47  bool hasOne() const { return theQual[0]!=-1 || theQual[1]!=-1; };
48  bool hasBoth() const { return theQual[0]!=-1 && theQual[1]!=-1; };
49  bool hasSameQual() const { return hasBoth() && theQual[0]==theQual[1]; };
50  int deltaBX() const { return theBX[0] - theBX[1]; }
51  int qualDDU() const { return theQual[0]; }
52  int qualTM() const { return theQual[1]; }
53 
54 private:
55 
56  int theQual[2]; // 0=DDU 1=TM
57  int theBX[2]; // 0=DDU 1=TM
58 
59 };
60 
61 
62 
63 DTLocalTriggerBaseTask::DTLocalTriggerBaseTask(const edm::ParameterSet& ps) :
64  nEvents(0), nLumis(0), theTrigGeomUtils(0) {
65 
66  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
67  << "[DTLocalTriggerBaseTask]: Constructor"<<endl;
68 
69  tpMode = ps.getUntrackedParameter<bool>("testPulseMode");
70  detailedAnalysis = ps.getUntrackedParameter<bool>("detailedAnalysis");
71 
72  targetBXTM = ps.getUntrackedParameter<int>("targetBXTM");
73  targetBXDDU = ps.getUntrackedParameter<int>("targetBXDDU");
74  bestAccRange = ps.getUntrackedParameter<int>("bestTrigAccRange");
75 
76  processTM = ps.getUntrackedParameter<bool>("processTM");
77  processDDU = ps.getUntrackedParameter<bool>("processDDU");
78 
79  tm_phiIn_Token_ = consumes<L1MuDTChambPhContainer>(
80  ps.getUntrackedParameter<InputTag>("inputTagTMphIn"));
81  tm_phiOut_Token_ = consumes<L1MuDTChambPhContainer>(
82  ps.getUntrackedParameter<InputTag>("inputTagTMphOut"));
83  tm_theta_Token_ = consumes<L1MuDTChambThContainer>(
84  ps.getUntrackedParameter<InputTag>("inputTagTMth"));
85  trig_Token_ = consumes<DTLocalTriggerCollection>(
86  ps.getUntrackedParameter<InputTag>("inputTagDDU"));
87 
88  if (processTM) theTypes.push_back("TM");
89  if (processDDU) theTypes.push_back("DDU");
90 
91  if (tpMode) {
92  topFolder("TM") = "DT/11-LocalTriggerTP-TM/";
93  topFolder("DDU") = "DT/12-LocalTriggerTP-DDU/";
94  } else {
95  topFolder("TM") = "DT/03-LocalTrigger-TM/";
96  topFolder("DDU") = "DT/04-LocalTrigger-DDU/";
97  }
98 
99  theParams = ps;
100 
101 }
102 
103 
104 DTLocalTriggerBaseTask::~DTLocalTriggerBaseTask() {
105 
106  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
107  << "[DTLocalTriggerBaseTask]: analyzed " << nEvents << " events" << endl;
108  if (theTrigGeomUtils) { delete theTrigGeomUtils; }
109 
110 }
111 
112 void DTLocalTriggerBaseTask::bookHistograms(DQMStore::IBooker & ibooker, edm::Run const & iRun, edm::EventSetup const & context) {
113 
114  ibooker.setCurrentFolder("DT/EventInfo/Counters");
115  nEventMonitor = ibooker.bookFloat("nProcessedEventsTrigger");
116  for (int wh=-2;wh<3;++wh){
117  for (int stat=1;stat<5;++stat){
118  for (int sect=1;sect<13;++sect){
119  bookHistos(ibooker, DTChamberId(wh,stat,sect));
120  }
121  }
122  bookHistos(ibooker, wh);
123  }
124 }
125 
126 
127 void DTLocalTriggerBaseTask::beginLuminosityBlock(const LuminosityBlock& lumiSeg, const EventSetup& context) {
128 
129  nEventsInLS=0;
130  nLumis++;
131  int resetCycle = theParams.getUntrackedParameter<int>("ResetCycle");
132 
133  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
134  << "[DTLocalTriggerBaseTask]: Begin of LS transition" << endl;
135 
136  if( nLumis%resetCycle == 0 ) {
137  map<uint32_t,map<string,MonitorElement*> >::const_iterator chambIt = chamberHistos.begin();
138  map<uint32_t,map<string,MonitorElement*> >::const_iterator chambEnd = chamberHistos.end();
139  for(;chambIt!=chambEnd;++chambIt) {
140  map<string,MonitorElement*>::const_iterator histoIt = chambIt->second.begin();
141  map<string,MonitorElement*>::const_iterator histoEnd = chambIt->second.end();
142  for(;histoIt!=histoEnd;++histoIt) {
143  histoIt->second->Reset();
144  }
145  }
146  }
147 
148 }
149 
150 void DTLocalTriggerBaseTask::endLuminosityBlock(const LuminosityBlock& lumiSeg, const EventSetup& context) {
151 
152  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
153  << "[DTLocalTriggerBaseTask]: End of LS transition" << endl;
154 
155  map<uint32_t,DTTimeEvolutionHisto* >::const_iterator chambIt = trendHistos.begin();
156  map<uint32_t,DTTimeEvolutionHisto* >::const_iterator chambEnd = trendHistos.end();
157  for(;chambIt!=chambEnd;++chambIt) {
158  chambIt->second->updateTimeSlot(lumiSeg.luminosityBlock(), nEventsInLS);
159  }
160 
161 }
162 
163 
164 void DTLocalTriggerBaseTask::dqmBeginRun(const edm::Run& run, const edm::EventSetup& context) {
165 
166  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
167  << "[DTLocalTriggerBaseTask]: BeginRun" << endl;
168 
169  ESHandle<DTGeometry> theGeom;
170  context.get<MuonGeometryRecord>().get(theGeom);
171  theTrigGeomUtils = new DTTrigGeomUtils(theGeom);
172 
173 }
174 
175 
176 void DTLocalTriggerBaseTask::analyze(const edm::Event& e, const edm::EventSetup& c){
177 
178  nEvents++;
179  nEventsInLS++;
180  nEventMonitor->Fill(nEvents);
181 
182  theCompMapIn.clear();
183  theCompMapOut.clear();
184 
185  Handle<L1MuDTChambPhContainer> phiInTrigsTM;
186  Handle<L1MuDTChambPhContainer> phiOutTrigsTM;
187  Handle<L1MuDTChambThContainer> thetaTrigsTM;
188  Handle<DTLocalTriggerCollection> trigsDDU;
189 
190  if (processTM) {
191  InputTag inputTagTM = theParams.getUntrackedParameter<InputTag>("inputTagTM");
192 
193  e.getByToken(tm_phiIn_Token_, phiInTrigsTM);
194  e.getByToken(tm_phiOut_Token_, phiOutTrigsTM);
195  e.getByToken(tm_theta_Token_, thetaTrigsTM);
196 
197  if (phiInTrigsTM.isValid() && phiOutTrigsTM.isValid() && thetaTrigsTM.isValid()) {
198  runTMAnalysis(phiInTrigsTM->getContainer(),phiOutTrigsTM->getContainer(),thetaTrigsTM->getContainer());
199  } else {
200  LogVerbatim("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
201  << "[DTLocalTriggerBaseTask]: one or more TM handles for Input Tag "
202  << inputTagTM <<" not found!" << endl;
203  return;
204  }
205  }
206 
207  if (processDDU) {
208  InputTag inputTagDDU = theParams.getUntrackedParameter<InputTag>("inputTagDDU");
209  e.getByToken(trig_Token_, trigsDDU);
210 
211  if (trigsDDU.isValid()) {
212  runDDUAnalysis(trigsDDU);
213  } else {
214  LogVerbatim("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
215  << "[DTLocalTriggerBaseTask]: one or more DDU handles for Input Tag "
216  << inputTagDDU <<" not found!" << endl;
217  return;
218  }
219  }
220 
221  if (processTM && processDDU)
222  runDDUvsTMAnalysis();
223 
224 }
225 
226 
227 void DTLocalTriggerBaseTask::bookHistos(DQMStore::IBooker & ibooker, const DTChamberId& dtCh) {
228 
229  uint32_t rawId = dtCh.rawId();
230 
231  stringstream wheel; wheel << dtCh.wheel();
232  stringstream station; station << dtCh.station();
233  stringstream sector; sector << dtCh.sector();
234 
235  map<string,int> minBX;
236  map<string,int> maxBX;
237 
238  minBX["TM"] = theParams.getUntrackedParameter<int>("minBXTM");
239  maxBX["TM"] = theParams.getUntrackedParameter<int>("maxBXTM");
240  minBX["DDU"] = theParams.getUntrackedParameter<int>("minBXDDU");
241  maxBX["DDU"] = theParams.getUntrackedParameter<int>("maxBXDDU");
242 
243  int nTimeBins = theParams.getUntrackedParameter<int>("nTimeBins");
244  int nLSTimeBin = theParams.getUntrackedParameter<int>("nLSTimeBin");
245 
246  string chTag = "_W" + wheel.str() + "_Sec" + sector.str() + "_St" + station.str();
247  string labelInOut = "";
248 
249  vector<string>::const_iterator typeIt = theTypes.begin();
250  vector<string>::const_iterator typeEnd = theTypes.end();
251 
252  for (; typeIt!=typeEnd; ++typeIt) {
253 
254  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
255  << "[DTLocalTriggerBaseTask]: booking histos for " << topFolder((*typeIt)) << "Wheel"
256  << wheel.str() << "/Sector" << sector.str() << "/Station"<< station.str() << endl;
257 
258  for (int InOut=0; InOut<2;InOut++){
259  // Book Phi View Related Plots
260 
261  if(InOut==0)
262  {ibooker.setCurrentFolder(topFolder(*typeIt) + "Wheel" + wheel.str() + "/Sector"
263  + sector.str() + "/Station" + station.str() + "/LocalTriggerPhiIn");
264  labelInOut = "_In";}
265  else if (InOut==1){
266  ibooker.setCurrentFolder(topFolder(*typeIt) + "Wheel" + wheel.str() + "/Sector"
267  + sector.str() + "/Station" + station.str() + "/LocalTriggerPhiOut");
268  labelInOut = "_Out";}
269 
270 
271  string histoTag = (*typeIt) + "_BXvsQual" + labelInOut;
272  chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,"BX vs trigger quality",
273  7,-0.5,6.5,(int)(maxBX[(*typeIt)]-minBX[*typeIt]+1),minBX[*typeIt]-.5,maxBX[*typeIt]+.5);
274  setQLabels((chamberHistos[rawId])[histoTag],1);
275 
276  if (!tpMode) {
277  histoTag = (*typeIt) + "_BestQual" + labelInOut;
278  chamberHistos[rawId][histoTag] = ibooker.book1D(histoTag+chTag,
279  "Trigger quality of best primitives",7,-0.5,6.5);
280  setQLabels(chamberHistos[rawId][histoTag],1);
281 
282  histoTag = (*typeIt) + "_Flag1stvsQual" + labelInOut;
283  chamberHistos[dtCh.rawId()][histoTag] = ibooker.book2D(histoTag+chTag,
284  "1st/2nd trig flag vs quality",7,-0.5,6.5,2,-0.5,1.5);
285  setQLabels(chamberHistos[rawId][histoTag],1);
286  }
287 
288  if (*typeIt=="TM") {
289  float minPh, maxPh; int nBinsPh;
290  theTrigGeomUtils->phiRange(dtCh,minPh,maxPh,nBinsPh);
291 
292  histoTag = (*typeIt) + "_QualvsPhirad" + labelInOut;
293  chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,
294  "Trigger quality vs local position",nBinsPh,minPh,maxPh,7,-0.5,6.5);
295  setQLabels(chamberHistos[rawId][histoTag],2);
296 
297  if (detailedAnalysis && !tpMode) {
298  histoTag = (*typeIt) + "_QualvsPhibend" + labelInOut;
299  chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,
300  "Trigger quality vs local direction",200,-40.,40.,7,-0.5,6.5);
301  setQLabels((chamberHistos[dtCh.rawId()])[histoTag],2);
302  }
303  }
304  } //InOut loop
305 
306  // Book Theta View Related Plots
307  ibooker.setCurrentFolder(topFolder(*typeIt) + "Wheel" + wheel.str() + "/Sector"
308  + sector.str() + "/Station" + station.str() + "/LocalTriggerTheta");
309 
310  string histoTag = "";
311  if((*typeIt)=="TM" && dtCh.station()!=4) {
312  histoTag = (*typeIt) + "_PositionvsBX";
313  chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,"Theta trigger position vs BX",
314  (int)(maxBX[(*typeIt)]-minBX[*typeIt]+1),minBX[*typeIt]-.5,maxBX[*typeIt]+.5,7,-0.5,6.5);
315  histoTag = (*typeIt) + "_PositionvsQual";
316  chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,"Theta trigger position vs quality",
317  2,0.5,2.5,7,-0.5,6.5);
318  setQLabelsTheta(chamberHistos[rawId][histoTag],1);
319  histoTag = (*typeIt) + "_ThetaBXvsQual";
320  chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,"BX vs trigger quality",
321  2,0.5,2.5,(int)(maxBX[(*typeIt)]-minBX[*typeIt]+1),minBX[*typeIt]-.5,maxBX[*typeIt]+.5);
322  setQLabelsTheta(chamberHistos[rawId][histoTag],1);
323 // histoTag = (*typeIt) + "_ThetaBestQual";
324 // chamberHistos[rawId][histoTag] = ibooker.book1D(histoTag+chTag,
325 // "Trigger quality of best primitives (theta)",2,0.5,2.5); // 0 = not fired, 1 = L, 2 = H
326 // setQLabelsTheta(chamberHistos[rawId][histoTag],1);
327  } else {
328  if(dtCh.station()!=4){
329  histoTag = (*typeIt) + "_ThetaBXvsQual";
330  chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,"BX vs trigger quality",7,-0.5,6.5,
331  (int)(maxBX[(*typeIt)]-minBX[*typeIt]+1),minBX[*typeIt]-.5,maxBX[*typeIt]+.5);
332  setQLabels((chamberHistos[dtCh.rawId()])[histoTag],1);
333 
334  histoTag = (*typeIt) + "_ThetaBestQual";
335  chamberHistos[rawId][histoTag] = ibooker.book1D(histoTag+chTag,
336  "Trigger quality of best primitives (theta)",7,-0.5,6.5);
337  setQLabels((chamberHistos[dtCh.rawId()])[histoTag],1);
338  }
339  }
340 
341  }
342 
343  if (processTM && processDDU) {
344  // Book TM/DDU Comparison Plots
345  for (int InOut=0; InOut<2;InOut++){
346  if(InOut==0)
347  ibooker.setCurrentFolder(topFolder("DDU") + "Wheel" + wheel.str() + "/Sector"
348  + sector.str() + "/Station" + station.str() + "/LocalTriggerPhiIn");
349  else if(InOut==1)
350  ibooker.setCurrentFolder(topFolder("DDU") + "Wheel" + wheel.str() + "/Sector"
351  + sector.str() + "/Station" + station.str() + "/LocalTriggerPhiOut");
352 
353 
354  string histoTag = "COM_QualDDUvsQualTM";
355  chamberHistos[rawId][histoTag] = ibooker.book2D(histoTag+chTag,
356  "DDU quality vs TM quality",8,-1.5,6.5,8,-1.5,6.5);
357  setQLabels((chamberHistos[rawId])[histoTag],1);
358  setQLabels((chamberHistos[rawId])[histoTag],2);
359 
360  histoTag = "COM_MatchingTrend";
361  trendHistos[rawId] = new DTTimeEvolutionHisto(ibooker,histoTag+chTag,
362  "Fraction of DDU-TM matches w.r.t. proc evts",
363  nTimeBins,nLSTimeBin,true,0);
364  } //InOut loop
365  }
366 
367 }
368 
369 void DTLocalTriggerBaseTask::bookHistos(DQMStore::IBooker & ibooker, int wh) {
370 
371  stringstream wheel; wheel << wh;
372  ibooker.setCurrentFolder(topFolder("DDU") + "Wheel" + wheel.str() + "/");
373  string whTag = "_W" + wheel.str();
374 
375  LogTrace("DTDQM|DTMonitorModule|DTLocalTriggerBaseTask")
376  << "[DTLocalTriggerBaseTask]: booking wheel histos for "
377  << topFolder("DDU") << "Wheel" << wh << endl;
378 
379  string histoTag = "COM_BXDiff";
380  MonitorElement *me = ibooker.bookProfile2D(histoTag+whTag,
381  "DDU-TM BX Difference",12,1,13,4,1,5,0.,20.);
382  me->setAxisTitle("Sector",1);
383  me->setAxisTitle("station",2);
384  wheelHistos[wh][histoTag] = me;
385 
386 }
387 
388 
389 void DTLocalTriggerBaseTask::runTMAnalysis( std::vector<L1MuDTChambPhDigi> const* phInTrigs,
390  std::vector<L1MuDTChambPhDigi> const* phOutTrigs,
391  std::vector<L1MuDTChambThDigi> const* thTrigs )
392 {
393  vector<L1MuDTChambPhDigi>::const_iterator iph = phInTrigs->begin();
394  vector<L1MuDTChambPhDigi>::const_iterator iphe = phInTrigs->end();
395 
396  for(; iph !=iphe ; ++iph) {
397 
398  int wh = iph->whNum();
399  int sec = iph->scNum() + 1; // DTTF->DT Convention
400  int st = iph->stNum();
401  int qual = iph->code();
402  int is1st = iph->Ts2Tag() ? 1 : 0;
403  int bx = iph->bxNum() - is1st;
404 
405  if (qual <0 || qual>6) continue; // Check that quality is in a valid range
406 
407  DTChamberId dtChId(wh,st,sec);
408  uint32_t rawId = dtChId.rawId();
409 
410  float pos = theTrigGeomUtils->trigPos(&(*iph));
411  float dir = theTrigGeomUtils->trigDir(&(*iph));
412 
413  if (abs(bx-targetBXTM)<= bestAccRange &&
414  theCompMapIn[rawId].qualTM() <= qual)
415  theCompMapIn[rawId].setTM(qual,bx);
416 
417  map<string, MonitorElement*> &innerME = chamberHistos[rawId];
418  if (tpMode) {
419  innerME["TM_BXvsQual_In"]->Fill(qual,bx); // SM BX vs Qual Phi view (1st tracks)
420  innerME["TM_QualvsPhirad_In"]->Fill(pos,qual); // SM Qual vs radial angle Phi view
421  } else {
422  innerME["TM_BXvsQual_In"]->Fill(qual,bx); // SM BX vs Qual Phi view (1st tracks)
423  innerME["TM_Flag1stvsQual_In"]->Fill(qual,is1st); // SM Qual 1st/2nd track flag Phi view
424  if (!is1st) innerME["TM_QualvsPhirad_In"]->Fill(pos,qual); // SM Qual vs radial angle Phi view ONLY for 1st tracks
425  if (detailedAnalysis) {
426  innerME["TM_QualvsPhibend_In"]->Fill(dir,qual); // SM Qual vs bending Phi view
427  }
428  }
429 
430  }
431 
432  iph = phOutTrigs->begin();
433  iphe = phOutTrigs->end();
434 
435  for(; iph !=iphe ; ++iph) {
436 
437  int wh = iph->whNum();
438  int sec = iph->scNum() + 1; // DTTF->DT Convention
439  int st = iph->stNum();
440  int qual = iph->code();
441  int is1st = iph->Ts2Tag() ? 1 : 0;
442  int bx = iph->bxNum() - is1st;
443 
444  if (qual <0 || qual>6) continue; // Check that quality is in a valid range
445 
446  DTChamberId dtChId(wh,st,sec);
447  uint32_t rawId = dtChId.rawId();
448 
449  float pos = theTrigGeomUtils->trigPos(&(*iph));
450  float dir = theTrigGeomUtils->trigDir(&(*iph));
451 
452  if (abs(bx-targetBXTM)<= bestAccRange &&
453  theCompMapOut[rawId].qualTM() <= qual)
454  theCompMapOut[rawId].setTM(qual,bx);
455 
456  map<string, MonitorElement*> &innerME = chamberHistos[rawId];
457  if (tpMode) {
458  innerME["TM_BXvsQual_Out"]->Fill(qual,bx); // SM BX vs Qual Phi view (1st tracks)
459  innerME["TM_QualvsPhirad_Out"]->Fill(pos,qual); // SM Qual vs radial angle Phi view
460  } else {
461  innerME["TM_BXvsQual_Out"]->Fill(qual,bx); // SM BX vs Qual Phi view (1st tracks)
462  innerME["TM_Flag1stvsQual_Out"]->Fill(qual,is1st); // SM Qual 1st/2nd track flag Phi view
463  if (!is1st) innerME["TM_QualvsPhirad_Out"]->Fill(pos,qual); // SM Qual vs radial angle Phi view ONLY for 1st tracks
464  if (detailedAnalysis) {
465  innerME["TM_QualvsPhibend_Out"]->Fill(dir,qual); // SM Qual vs bending Phi view
466  }
467  }
468 
469  }
470 
471  vector<L1MuDTChambThDigi>::const_iterator ith = thTrigs->begin();
472  vector<L1MuDTChambThDigi>::const_iterator ithe = thTrigs->end();
473 
474  for(; ith != ithe; ++ith) {
475  int wh = ith->whNum();
476  int sec = ith->scNum() + 1; // DTTF -> DT Convention
477  int st = ith->stNum();
478  int bx = ith->bxNum();
479 
480  int thcode[7];
481 
482  for (int pos=0; pos<7; pos++)
483  thcode[pos] = ith->code(pos);
484 
485  DTChamberId dtChId(wh,st,sec);
486  uint32_t rawId = dtChId.rawId();
487 
488  map<string, MonitorElement*> &innerME = chamberHistos[rawId];
489 
490  for (int pos=0; pos<7; pos++)
491  if (thcode[pos]>0) {//Fired
492  innerME["TM_PositionvsBX"]->Fill(bx,pos); // SM BX vs Position Theta view
493  innerME["TM_PositionvsQual"]->Fill(thcode[pos],pos); //code = pos + qual; so 0, 1, 2 for 0, L, H resp.
494  innerME["TM_ThetaBXvsQual"]->Fill(thcode[pos],bx); //code = pos + qual; so 0, 1, 2 for 0, L, H resp.
495  }
496  }
497  // Fill Quality plots with best TM triggers (phi view In)
498  if (!tpMode) {
499  map<uint32_t,DTTPGCompareUnit>::const_iterator compIt = theCompMapIn.begin();
500  map<uint32_t,DTTPGCompareUnit>::const_iterator compEnd = theCompMapIn.end();
501  for (; compIt!=compEnd; ++compIt) {
502  int bestQual = compIt->second.qualTM();
503  if (bestQual > -1)
504  chamberHistos[compIt->first]["TM_BestQual_In"]->Fill(bestQual); // SM Best Qual Trigger Phi view
505  }
506  }
507 
508  // Fill Quality plots with best TM triggers (phi view Out)
509  if (!tpMode) {
510  map<uint32_t,DTTPGCompareUnit>::const_iterator compIt = theCompMapOut.begin();
511  map<uint32_t,DTTPGCompareUnit>::const_iterator compEnd = theCompMapOut.end();
512  for (; compIt!=compEnd; ++compIt) {
513  int bestQual = compIt->second.qualTM();
514  if (bestQual > -1)
515  chamberHistos[compIt->first]["TM_BestQual_Out"]->Fill(bestQual); // SM Best Qual Trigger Phi view
516  }
517  }
518 
519 }
520 
521 
522 void DTLocalTriggerBaseTask::runDDUAnalysis(Handle<DTLocalTriggerCollection>& trigsDDU){
523 
524  DTLocalTriggerCollection::DigiRangeIterator detUnitIt = trigsDDU->begin();
525  DTLocalTriggerCollection::DigiRangeIterator detUnitEnd = trigsDDU->end();
526 
527  for (; detUnitIt!=detUnitEnd; ++detUnitIt){
528 
529  const DTChamberId& chId = (*detUnitIt).first;
530  uint32_t rawId = chId.rawId();
531 
532  const DTLocalTriggerCollection::Range& range = (*detUnitIt).second;
533  DTLocalTriggerCollection::const_iterator trigIt = range.first;
534  map<string, MonitorElement*> &innerME = chamberHistos[rawId];
535 
536  int bestQualTheta = -1;
537 
538  for (; trigIt!=range.second; ++trigIt){
539 
540  int qualPhi = trigIt->quality();
541  int qualTheta = trigIt->trTheta();
542  int flag1st = trigIt->secondTrack() ? 1 : 0;
543  int bx = trigIt->bx();
544  int bxPhi = bx - flag1st; // phi BX assign is different for 1st & 2nd tracks
545 
546  if( qualPhi>-1 && qualPhi<7 ) { // it is a phi trigger
547  if (abs(bx-targetBXDDU) <= bestAccRange &&
548  theCompMapIn[rawId].qualDDU()<= qualPhi)
549  theCompMapIn[rawId].setDDU(qualPhi,bxPhi);
550  if(tpMode) {
551  innerME["DDU_BXvsQual"]->Fill(qualPhi,bxPhi); // SM BX vs Qual Phi view
552  } else {
553  innerME["DDU_BXvsQual"]->Fill(qualPhi,bxPhi); // SM BX vs Qual Phi view
554  innerME["DDU_Flag1stvsQual"]->Fill(qualPhi,flag1st); // SM Quality vs 1st/2nd track flag Phi view
555  }
556  }
557 
558  if( qualTheta>0 && !tpMode ){// it is a theta trigger & is not TP
559  if (qualTheta > bestQualTheta){
560  bestQualTheta = qualTheta;
561  }
562  innerME["DDU_ThetaBXvsQual"]->Fill(qualTheta,bx); // SM BX vs Qual Theta view
563  }
564  }
565 
566  // Fill Quality plots with best ddu triggers
567  if (!tpMode && theCompMapIn.find(rawId)!= theCompMapIn.end()) {
568  int bestQualPhi = theCompMapIn[rawId].qualDDU();
569  if (bestQualPhi>-1)
570  innerME["DDU_BestQual"]->Fill(bestQualPhi); // SM Best Qual Trigger Phi view
571  if(bestQualTheta>0) {
572  innerME["DDU_ThetaBestQual"]->Fill(bestQualTheta); // SM Best Qual Trigger Theta view
573  }
574  }
575  }
576 
577 }
578 
579 
580 void DTLocalTriggerBaseTask::runDDUvsTMAnalysis(){
581 
582  map<uint32_t,DTTPGCompareUnit>::const_iterator compIt = theCompMapIn.begin();
583  map<uint32_t,DTTPGCompareUnit>::const_iterator compEnd = theCompMapIn.end();
584 
585  for (; compIt!=compEnd; ++compIt) {
586 
587  uint32_t rawId = compIt->first;
588  DTChamberId chId(rawId);
589  map<string, MonitorElement*> &innerME = chamberHistos[rawId];
590 
591  const DTTPGCompareUnit & compUnit = compIt->second;
592  if ( compUnit.hasOne() ){
593  innerME["COM_QualDDUvsQualTM"]->Fill(compUnit.qualTM(),compUnit.qualDDU());
594  }
595  if ( compUnit.hasBoth() ){
596  wheelHistos[chId.wheel()]["COM_BXDiff"]->Fill(chId.sector(),chId.station(),compUnit.deltaBX());
597  if ( compUnit.hasSameQual() ) {
598  trendHistos[rawId]->accumulateValueTimeSlot(1);
599  }
600  }
601  }
602 
603 }
604 
605 
606 void DTLocalTriggerBaseTask::setQLabels(MonitorElement* me, short int iaxis){
607 
608  TH1* histo = me->getTH1();
609  if (!histo) return;
610 
611  TAxis* axis=0;
612  if (iaxis==1) {
613  axis=histo->GetXaxis();
614  }
615  else if(iaxis==2) {
616  axis=histo->GetYaxis();
617  }
618  if (!axis) return;
619 
620  string labels[7] = {"LI","LO","HI","HO","LL","HL","HH"};
621  int istart = axis->GetXmin()<-1 ? 2 : 1;
622  for (int i=0;i<7;i++) {
623  axis->SetBinLabel(i+istart,labels[i].c_str());
624  }
625 
626 }
627 
628 void DTLocalTriggerBaseTask::setQLabelsTheta(MonitorElement* me, short int iaxis){
629 
630  TH1* histo = me->getTH1();
631  if (!histo) return;
632 
633  TAxis* axis=0;
634  if (iaxis==1) {
635  axis=histo->GetXaxis();
636  }
637  else if(iaxis==2) {
638  axis=histo->GetYaxis();
639  }
640  if (!axis) return;
641 
642  string labels[2] = {"L","H"};
643  int istart = axis->GetXmin()<-1 ? 2 : 1;
644  for (int i=0;i<2;i++) {
645  axis->SetBinLabel(i+istart,labels[i].c_str());
646  }
647 
648 }
649 // Local Variables:
650 // show-trailing-whitespace: t
651 // truncate-lines: t
652 // End:
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
i
int i
Definition: DBlmapReader.cc:9
DTLocalTriggerBaseTask::runDDUvsTMAnalysis
void runDDUvsTMAnalysis()
Run analysis on ROS data.
Definition: DTLocalTriggerBaseTask.cc:580
DTChamberId
DTChamberId
Definition: RecoLocalMuon_RECO.doi:23
DTLocalTriggerBaseTask::dqmBeginRun
void dqmBeginRun(const edm::Run &, const edm::EventSetup &) override
Beginrun.
Definition: DTLocalTriggerBaseTask.cc:164
DTLocalTriggerBaseTask::bookHistos
void bookHistos(DQMStore::IBooker &, const DTChamberId &chamb)
Book the histograms.
Definition: DTLocalTriggerBaseTask.cc:227
edm::Event::getByToken
bool getByToken(EDGetToken token, Handle< PROD > &result) const
Definition: Event.h:457
DTTrigGeomUtils::trigDir
float trigDir(const L1MuDTChambPhDigi *trig)
Return local direction (trigger RF) for a given trigger primitive.
Definition: DTTrigGeomUtils.cc:151
DTTPGCompareUnit::setDDU
void setDDU(int qual, int bx)
Definition: DTLocalTriggerBaseTask.cc:44
DTLocalTriggerBaseTask::wheelHistos
std::map< int, std::map< std::string, MonitorElement * > > wheelHistos
Definition: DTLocalTriggerBaseTask.h:117
DTTPGCompareUnit::setTM
void setTM(int qual, int bx)
Definition: DTLocalTriggerBaseTask.cc:45
hlt_dqm_clientPB-live_cfg.me
tuple me
Definition: hlt_dqm_clientPB-live_cfg.py:43
DTLocalTriggerBaseTask::~DTLocalTriggerBaseTask
virtual ~DTLocalTriggerBaseTask()
Destructor.
Definition: DTLocalTriggerBaseTask.cc:104
DTLocalTriggerBaseTask::theTypes
std::vector< std::string > theTypes
Definition: DTLocalTriggerBaseTask.h:113
MonitorElement::Fill
void Fill(long long x)
Definition: MonitorElement.h:170
edm::LuminosityBlockBase::luminosityBlock
LuminosityBlockNumber_t luminosityBlock() const
Definition: LuminosityBlockBase.h:40
DetId::rawId
uint32_t rawId() const
get the raw id
Definition: DetId.h:43
DTLocalTriggerBaseTask::theTrigGeomUtils
DTTrigGeomUtils * theTrigGeomUtils
Definition: DTLocalTriggerBaseTask.h:112
DTLocalTriggerBaseTask::endLuminosityBlock
void endLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context) override
Perform trend plot operations.
Definition: DTLocalTriggerBaseTask.cc:150
DTLocalTriggerBaseTask::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c) override
Analyze.
Definition: DTLocalTriggerBaseTask.cc:176
DQMStore::IBooker::bookProfile2D
MonitorElement * bookProfile2D(Args &&...args)
Definition: DQMStore.h:163
DTLocalTriggerBaseTask::tm_theta_Token_
edm::EDGetTokenT< L1MuDTChambThContainer > tm_theta_Token_
Definition: DTLocalTriggerBaseTask.h:124
DTLocalTriggerBaseTask::bookHistograms
void bookHistograms(DQMStore::IBooker &, edm::Run const &, edm::EventSetup const &) override
Definition: DTLocalTriggerBaseTask.cc:112
DQMStore::IBooker::book1D
MonitorElement * book1D(Args &&...args)
Definition: DQMStore.h:115
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
MonitorElement::getTH1
TH1 * getTH1(void) const
Definition: MonitorElement.cc:1450
DTLocalTriggerBaseTask::theCompMapOut
std::map< uint32_t, DTTPGCompareUnit > theCompMapOut
Definition: DTLocalTriggerBaseTask.h:116
DTLocalTriggerBaseTask::trig_Token_
edm::EDGetTokenT< DTLocalTriggerCollection > trig_Token_
Definition: DTLocalTriggerBaseTask.h:125
DTLocalTriggerBaseTask::runTMAnalysis
void runTMAnalysis(std::vector< L1MuDTChambPhDigi > const *phInTrigs, std::vector< L1MuDTChambPhDigi > const *phOutTrigs, std::vector< L1MuDTChambThDigi > const *thTrigs)
Run analysis on TM data.
Definition: DTLocalTriggerBaseTask.cc:389
edm::HandleBase::isValid
bool isValid() const
Definition: HandleBase.h:75
DTLocalTriggerBaseTask::beginLuminosityBlock
void beginLuminosityBlock(const edm::LuminosityBlock &lumiSeg, const edm::EventSetup &context) override
To reset the MEs.
Definition: DTLocalTriggerBaseTask.cc:127
LogTrace
#define LogTrace(id)
Definition: PFTrackAlgoTools.cc:502
DTTrigGeomUtils::phiRange
void phiRange(const DTChamberId &id, float &min, float &max, int &nbins, float step=15)
Compute phi range in local chamber coordinates.
Definition: DTTrigGeomUtils.cc:75
DTLocalTriggerBaseTask::theCompMapIn
std::map< uint32_t, DTTPGCompareUnit > theCompMapIn
Definition: DTLocalTriggerBaseTask.h:115
DTLocalTriggerBaseTask::topFolder
std::string & topFolder(std::string const &type)
Get the Top folder (different between Physics and TP and TM/DDU)
Definition: DTLocalTriggerBaseTask.h:83
DTLocalTriggerBaseTask::chamberHistos
std::map< uint32_t, std::map< std::string, MonitorElement * > > chamberHistos
Definition: DTLocalTriggerBaseTask.h:118
DTLocalTriggerBaseTask::setQLabelsTheta
void setQLabelsTheta(MonitorElement *me, short int iaxis)
Definition: DTLocalTriggerBaseTask.cc:628
DQMStore::IBooker::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:277
DQMStore::IBooker::book2D
MonitorElement * book2D(Args &&...args)
Definition: DQMStore.h:133
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:56
MuonDigiCollection< DTChamberId, DTLocalTrigger >::const_iterator
std::vector< DTLocalTrigger >::const_iterator const_iterator
Definition: MuonDigiCollection.h:110
DTChamberId::sector
int sector() const
Definition: DTChamberId.h:61
DTLocalTriggerBaseTask::setQLabels
void setQLabels(MonitorElement *me, short int iaxis)
Set Quality labels.
Definition: DTLocalTriggerBaseTask.cc:606
DTLocalTriggerBaseTask::tm_phiIn_Token_
edm::EDGetTokenT< L1MuDTChambPhContainer > tm_phiIn_Token_
Definition: DTLocalTriggerBaseTask.h:122
DTLocalTriggerBaseTask::trendHistos
std::map< uint32_t, DTTimeEvolutionHisto * > trendHistos
Definition: DTLocalTriggerBaseTask.h:119
MuonDigiCollection< DTChamberId, DTLocalTrigger >::Range
std::pair< const_iterator, const_iterator > Range
Definition: MuonDigiCollection.h:111
DQMStore::IBooker::bookFloat
MonitorElement * bookFloat(Args &&...args)
Definition: DQMStore.h:109
dir
dbl *** dir
Definition: mlp_gen.cc:35
DTLocalTriggerBaseTask::runDDUAnalysis
void runDDUAnalysis(edm::Handle< DTLocalTriggerCollection > &trigsDDU)
Run analysis on ROS data.
Definition: DTLocalTriggerBaseTask.cc:522
DTTrigGeomUtils::trigPos
float trigPos(const L1MuDTChambPhDigi *trig)
Return local position (trigger RF) for a given trigger primitive.
Definition: DTTrigGeomUtils.cc:112
DTChamberId::station
int station() const
Return the station number.
Definition: DTChamberId.h:51
nEvents
UInt_t nEvents
Definition: hcalCalib.cc:42
DTLocalTriggerBaseTask::tm_phiOut_Token_
edm::EDGetTokenT< L1MuDTChambPhContainer > tm_phiOut_Token_
Definition: DTLocalTriggerBaseTask.h:123
DTLocalTriggerBaseTask::DTLocalTriggerBaseTask
DTLocalTriggerBaseTask(const edm::ParameterSet &ps)
Constructor.
Definition: DTLocalTriggerBaseTask.cc:63
DTChamberId::wheel
int wheel() const
Return the wheel number.
Definition: DTChamberId.h:45
HLT_25ns10e33_v2_cff.labels
tuple labels
Definition: HLT_25ns10e33_v2_cff.py:5989
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
edm::Run
Definition: Run.h:42
Generated for CMSSW Reference Manual by doxygen 1.8.5