CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
DTResolutionAnalysisTest.cc
Go to the documentation of this file.
1 
2 
3 /*
4  * See header file for a description of this class.
5  *
6  * \author G. Mila - INFN Torino
7  */
8 
9 
10 #include <DQM/DTMonitorClient/src/DTResolutionAnalysisTest.h>
11 
12 // Framework
13 #include <FWCore/Framework/interface/Event.h>
14 #include <FWCore/Framework/interface/EventSetup.h>
15 #include <FWCore/ParameterSet/interface/ParameterSet.h>
16 
17 
18 // Geometry
19 #include "Geometry/Records/interface/MuonGeometryRecord.h"
20 #include "Geometry/DTGeometry/interface/DTGeometry.h"
21 
22 #include "FWCore/ServiceRegistry/interface/Service.h"
23 #include "DQMServices/Core/interface/DQMStore.h"
24 #include "DQMServices/Core/interface/MonitorElement.h"
25 #include "FWCore/MessageLogger/interface/MessageLogger.h"
26 
27 #include <string>
28 #include <sstream>
29 #include <math.h>
30 
31 
32 using namespace edm;
33 using namespace std;
34 
35 
36 DTResolutionAnalysisTest::DTResolutionAnalysisTest(const ParameterSet& ps){
37 
38  LogTrace ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") << "[DTResolutionAnalysisTest]: Constructor";
39 
40  dbe = Service<DQMStore>().operator->();
41 
42  prescaleFactor = ps.getUntrackedParameter<int>("diagnosticPrescale", 1);
43  // permitted test range
44  maxGoodMeanValue = ps.getUntrackedParameter<double>("maxGoodMeanValue",0.02);
45  minBadMeanValue = ps.getUntrackedParameter<double>("minBadMeanValue",0.04);
46  maxGoodSigmaValue = ps.getUntrackedParameter<double>("maxGoodSigmaValue",0.08);
47  minBadSigmaValue = ps.getUntrackedParameter<double>("minBadSigmaValue",0.16);
48  // top folder for the histograms in DQMStore
49  topHistoFolder = ps.getUntrackedParameter<string>("topHistoFolder","DT/02-Segments");
50 
51  doCalibAnalysis = ps.getUntrackedParameter<bool>("doCalibAnalysis",false);
52 }
53 
54 
55 DTResolutionAnalysisTest::~DTResolutionAnalysisTest(){
56 
57  LogTrace ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") << "DTResolutionAnalysisTest: analyzed " << nevents << " events";
58 
59 }
60 
61 
62 void DTResolutionAnalysisTest::beginJob(){
63 
64  LogTrace ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") <<"[DTResolutionAnalysisTest]: BeginJob";
65 
66  nevents = 0;
67 
68 }
69 
70 void DTResolutionAnalysisTest::beginRun(const Run& run, const EventSetup& context){
71 
72  LogTrace ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") <<"[DTResolutionAnalysisTest]: BeginRun";
73 
74  // Get the geometry
75  context.get<MuonGeometryRecord>().get(muonGeom);
76 
77 }
78 
79 
80 void DTResolutionAnalysisTest::beginLuminosityBlock(LuminosityBlock const& lumiSeg, EventSetup const& context) {
81 
82  LogTrace ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest") <<"[DTResolutionAnalysisTest]: Begin of LS transition";
83 
84  // Get the run number
85  run = lumiSeg.run();
86 
87 }
88 
89 
90 
91 void DTResolutionAnalysisTest::analyze(const Event& e, const EventSetup& context){
92 
93  nevents++;
94 
95 }
96 
97 
98 
99 void DTResolutionAnalysisTest::endLuminosityBlock(LuminosityBlock const& lumiSeg, EventSetup const& context) {
100 
101 }
102 
103 
104 void DTResolutionAnalysisTest::bookHistos() {
105 
106  // global residual summary
107  dbe->setCurrentFolder(topHistoFolder);
108  globalResSummary = dbe->book2D("ResidualsGlbSummary", "# of SLs with good mean and good sigma of residuals",12,1,13,5,-2,3);
109 
110  // book summaries for mean and sigma
111  dbe->setCurrentFolder(topHistoFolder + "/00-MeanRes");
112  meanDistr[-2] = dbe->book1D("MeanDistr","Mean value of the residuals all (cm)",
113  100,-0.1,0.1);
114  meanDistr[-1] = dbe->book1D("MeanDistr_Phi","Mean value of the residuals #phi SL (cm)",
115  100,-0.1,0.1);
116  meanDistr[0] = dbe->book1D("MeanDistr_ThetaWh0","Mean values of the residuals #theta SL Wh 0 (cm)",
117  100,-0.1,0.1);
118  meanDistr[1] = dbe->book1D("MeanDistr_ThetaWh1","Mean value of the residuals #theta SL Wh +/-1 (cm)",
119  100,-0.1,0.1);
120  meanDistr[2] = dbe->book1D("MeanDistr_ThetaWh2","Mean value of the residuals #theta SL Wh +/-2 (cm)",
121  100,-0.1,0.1);
122 
123 
124  string histoTitle = "# of SLs with good mean of residuals";
125  wheelMeanHistos[3] = dbe->book2D("MeanResGlbSummary",histoTitle.c_str(),12,1,13,5,-2,3);
126  wheelMeanHistos[3]->setAxisTitle("Sector",1);
127  wheelMeanHistos[3]->setAxisTitle("Wheel",2);
128 
129 
130  dbe->setCurrentFolder(topHistoFolder + "/01-SigmaRes");
131  sigmaDistr[-2] = dbe->book1D("SigmaDistr","Sigma value of the residuals all (cm)",
132  50,0.0,0.2);
133  sigmaDistr[-1] = dbe->book1D("SigmaDistr_Phi","Sigma value of the residuals #phi SL (cm)",
134  50,0.0,0.2);
135  sigmaDistr[0] = dbe->book1D("SigmaDistr_ThetaWh0","Sigma value of the residuals #theta SL Wh 0 (cm)",
136  50,0.0,0.2);
137  sigmaDistr[1] = dbe->book1D("SigmaDistr_ThetaWh1","Sigma value of the residuals #theta SL Wh +/-1 (cm)",
138  50,0.0,0.2);
139  sigmaDistr[2] = dbe->book1D("SigmaDistr_ThetaWh2","Sigma value of the residuals #theta SL Wh +/-2 (cm)",
140  50,0.0,0.2);
141 
142  histoTitle = "# of SLs with good sigma of residuals";
143  wheelSigmaHistos[3] = dbe->book2D("SigmaResGlbSummary",histoTitle.c_str(),12,1,13,5,-2,3);
144  wheelSigmaHistos[3]->setAxisTitle("Sector",1);
145  wheelSigmaHistos[3]->setAxisTitle("Wheel",2);
146 
147 
148  // loop over all the CMS wheels, sectors & book the summary histos
149  for (int wheel=-2; wheel<=2; wheel++){
150  bookHistos(wheel);
151  for (int sector=1; sector<=12; sector++){
152  bookHistos(wheel, sector);
153  }
154  }
155 
156 }
157 
158 void DTResolutionAnalysisTest::endRun(Run const& run, EventSetup const& context) {
159 
160  if (!dbe->dirExists(topHistoFolder)) {
161  LogTrace ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest")
162  <<"[DTResolutionAnalysisTest]: Base folder " << topHistoFolder
163  << " does not exist. Skipping client operation." << endl;
164  return;
165  }
166 
167  bookHistos(); // histos booked only if top histo folder exist
168  // as Standard/AlcaReco Harvest is performed in the same step
169 
170  LogTrace ("DTDQM|DTMonitorClient|DTResolutionAnalysisTest")
171  << "[DTResolutionAnalysisTest]: End of Run transition, performing the DQM client operation" << endl;
172 
173  // reset the ME with fixed scale
174  resetMEs();
175 
176  for (vector<DTChamber*>::const_iterator ch_it = muonGeom->chambers().begin();
177  ch_it != muonGeom->chambers().end(); ++ch_it) { // loop over the chambers
178 
179  DTChamberId chID = (*ch_it)->id();
180 
181  // Fill the test histos
182  for(vector<const DTSuperLayer*>::const_iterator sl_it = (*ch_it)->superLayers().begin();
183  sl_it != (*ch_it)->superLayers().end(); ++sl_it) { // loop over SLs
184 
185 
186  DTSuperLayerId slID = (*sl_it)->id();
187  MonitorElement * res_histo = dbe->get(getMEName(slID));
188 
189  if(res_histo) { // Gaussian Fit
190  float statMean = res_histo->getMean(1);
191  float statSigma = res_histo->getRMS(1);
192  double mean = -1;
193  double sigma = -1;
194  TH1F * histo_root = res_histo->getTH1F();
195 
196  // fill the summaries
197  int entry= (chID.station() - 1) * 3;
198  int binSect = slID.sector();
199  if(slID.sector() == 13) binSect = 4;
200  else if(slID.sector() == 14) binSect = 10;
201  int binSL = entry+slID.superLayer();
202  if(chID.station() == 4 && slID.superLayer() == 3) binSL--;
203  if((slID.sector()==13 || slID.sector()==14) && slID.superLayer()==1) binSL=12;
204  if((slID.sector()==13 || slID.sector()==14) && slID.superLayer()==3) binSL=13;
205 
206  if(histo_root->GetEntries()>20) {
207  TF1 *gfit = new TF1("Gaussian","gaus",(statMean-(2*statSigma)),(statMean+(2*statSigma)));
208  try {
209  histo_root->Fit(gfit, "Q0", "", -0.1, 0.1);
210  } catch (cms::Exception& iException) {
211  LogWarning ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask")
212  << "[DTResolutionAnalysisTask]: Exception when fitting SL : " << slID;
213  // FIXME: the SL is set as OK in the summary
214  double weight = 1/11.;
215  if((binSect == 4 || binSect == 10) && slID.station() == 4) weight = 1/22.;
216  globalResSummary->Fill(binSect, slID.wheel(), weight);
217  continue;
218  }
219 
220  if(gfit){
221  // get the mean and the sigma of the distribution
222  mean = gfit->GetParameter(1);
223  sigma = gfit->GetParameter(2);
224 
225  // fill the distributions
226  meanDistr[-2]->Fill(mean);
227  sigmaDistr[-2]->Fill(sigma);
228  if(slID.superlayer() == 2) {
229  meanDistr[abs(slID.wheel())]->Fill(mean);
230  sigmaDistr[abs(slID.wheel())]->Fill(sigma);
231  } else {
232  meanDistr[-1]->Fill(mean);
233  sigmaDistr[-1]->Fill(sigma);
234  }
235 
236  // sector summaries
237  MeanHistos[make_pair(slID.wheel(),binSect)]->setBinContent(binSL, mean);
238  SigmaHistos[make_pair(slID.wheel(),binSect)]->setBinContent(binSL, sigma);
239 
240  if((slID.sector() == 13 || slID.sector() == 14) && binSL == 12) binSL=10;
241  if((slID.sector() == 13 || slID.sector() == 14) && binSL == 13) binSL=11;
242 
243 
244  if((slID.sector() == 13 || slID.sector() == 14) ) {
245 
246  double MeanVal = wheelMeanHistos[slID.wheel()]->getBinContent(binSect,binSL);
247  double MeanBinVal = (MeanVal > 0. && MeanVal < meanInRange(mean)) ? MeanVal : meanInRange(mean);
248  wheelMeanHistos[slID.wheel()]->setBinContent(binSect,binSL,MeanBinVal);
249 
250  double SigmaVal = wheelSigmaHistos[slID.wheel()]->getBinContent(binSect,binSL);
251  double SigmaBinVal = (SigmaVal > 0. && SigmaVal < sigmaInRange(sigma)) ? SigmaVal : sigmaInRange(sigma);
252  wheelSigmaHistos[slID.wheel()]->setBinContent(binSect,binSL,SigmaBinVal);
253 
254  } else {
255  wheelMeanHistos[slID.wheel()]->setBinContent(binSect,binSL,meanInRange(mean));
256  wheelSigmaHistos[slID.wheel()]->setBinContent(binSect,binSL,sigmaInRange(sigma));
257  }
258 
259  // set the weight
260  double weight = 1/11.;
261  if((binSect == 4 || binSect == 10) && slID.station() == 4) weight = 1/22.;
262 
263  // test the values of mean and sigma
264  if( (meanInRange(mean) > 0.85) && (sigmaInRange(sigma) > 0.85) ) { // sigma and mean ok
265  globalResSummary->Fill(binSect, slID.wheel(), weight);
266  wheelMeanHistos[3]->Fill(binSect,slID.wheel(),weight);
267  wheelSigmaHistos[3]->Fill(binSect,slID.wheel(),weight);
268  } else {
269  if( (meanInRange(mean) < 0.85) && (sigmaInRange(sigma) > 0.85) ) { // only sigma ok
270  wheelSigmaHistos[3]->Fill(binSect,slID.wheel(),weight);
271  }
272  if((meanInRange(mean) > 0.85) && (sigmaInRange(sigma) < 0.85) ) { // only mean ok
273  wheelMeanHistos[3]->Fill(binSect,slID.wheel(),weight);
274  }
275  }
276  }
277  delete gfit;
278  }
279  else{
280  LogVerbatim ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask")
281  << "[DTResolutionAnalysisTask] Fit of " << slID
282  << " not performed because # entries < 20 ";
283  // FIXME: the SL is set as OK in the summary
284  double weight = 1/11.;
285  if((binSect == 4 || binSect == 10) && slID.station() == 4) weight = 1/22.;
286  globalResSummary->Fill(binSect, slID.wheel(), weight);
287  wheelMeanHistos[3]->Fill(binSect,slID.wheel(),weight);
288  wheelSigmaHistos[3]->Fill(binSect,slID.wheel(),weight);
289  wheelMeanHistos[slID.wheel()]->setBinContent(binSect,binSL,1.);
290  wheelSigmaHistos[slID.wheel()]->setBinContent(binSect,binSL,1.);
291  }
292  } else {
293  LogWarning ("DTDQM|DTMonitorModule|DTResolutionAnalysisTask")
294  << "[DTResolutionAnalysisTask] Histo: " << getMEName(slID) << " not found" << endl;
295  }
296  } // loop on SLs
297  } // Loop on Stations
298 
299 }
300 
301 
302 
303 void DTResolutionAnalysisTest::bookHistos(int wh) {
304 
305  stringstream wheel; wheel <<wh;
306 
307  dbe->setCurrentFolder(topHistoFolder + "/00-MeanRes");
308  string histoName = "MeanSummaryRes_W" + wheel.str();
309  string histoTitle = "# of SLs with wrong mean of residuals (Wheel " + wheel.str() + ")";
310  wheelMeanHistos[wh] = dbe->book2D(histoName.c_str(),histoTitle.c_str(),12,1,13,11,1,12);
311  wheelMeanHistos[wh]->setAxisTitle("Sector",1);
312  wheelMeanHistos[wh]->setBinLabel(1,"MB1_SL1",2);
313  wheelMeanHistos[wh]->setBinLabel(2,"MB1_SL2",2);
314  wheelMeanHistos[wh]->setBinLabel(3,"MB1_SL3",2);
315  wheelMeanHistos[wh]->setBinLabel(4,"MB2_SL1",2);
316  wheelMeanHistos[wh]->setBinLabel(5,"MB2_SL2",2);
317  wheelMeanHistos[wh]->setBinLabel(6,"MB2_SL3",2);
318  wheelMeanHistos[wh]->setBinLabel(7,"MB3_SL1",2);
319  wheelMeanHistos[wh]->setBinLabel(8,"MB3_SL2",2);
320  wheelMeanHistos[wh]->setBinLabel(9,"MB3_SL3",2);
321  wheelMeanHistos[wh]->setBinLabel(10,"MB4_SL1",2);
322  wheelMeanHistos[wh]->setBinLabel(11,"MB4_SL3",2);
323  // wheelMeanHistos[wh]->setBinLabel(12,"MB4_SL1",2);
324  // wheelMeanHistos[wh]->setBinLabel(13,"MB4_SL3",2);
325 
326 
327 
328  dbe->setCurrentFolder(topHistoFolder + "/01-SigmaRes");
329  histoName = "SigmaSummaryRes_W" + wheel.str();
330  histoTitle = "# of SLs with wrong sigma of residuals (Wheel " + wheel.str() + ")";
331  wheelSigmaHistos[wh] = dbe->book2D(histoName.c_str(),histoTitle.c_str(),12,1,13,11,1,12);
332  wheelSigmaHistos[wh]->setAxisTitle("Sector",1);
333  wheelSigmaHistos[wh]->setBinLabel(1,"MB1_SL1",2);
334  wheelSigmaHistos[wh]->setBinLabel(2,"MB1_SL2",2);
335  wheelSigmaHistos[wh]->setBinLabel(3,"MB1_SL3",2);
336  wheelSigmaHistos[wh]->setBinLabel(4,"MB2_SL1",2);
337  wheelSigmaHistos[wh]->setBinLabel(5,"MB2_SL2",2);
338  wheelSigmaHistos[wh]->setBinLabel(6,"MB2_SL3",2);
339  wheelSigmaHistos[wh]->setBinLabel(7,"MB3_SL1",2);
340  wheelSigmaHistos[wh]->setBinLabel(8,"MB3_SL2",2);
341  wheelSigmaHistos[wh]->setBinLabel(9,"MB3_SL3",2);
342  wheelSigmaHistos[wh]->setBinLabel(10,"MB4_SL1",2);
343  wheelSigmaHistos[wh]->setBinLabel(11,"MB4_SL3",2);
344  // wheelSigmaHistos[wh]->setBinLabel(12,"MB4_SL1",2);
345  // wheelSigmaHistos[wh]->setBinLabel(13,"MB4_SL3",2);
346 
347 }
348 
349 
350 void DTResolutionAnalysisTest::bookHistos(int wh, int sect) {
351 
352  stringstream wheel; wheel << wh;
353  stringstream sector; sector << sect;
354 
355 
356  string MeanHistoName = "MeanTest_W" + wheel.str() + "_Sec" + sector.str();
357  string SigmaHistoName = "SigmaTest_W" + wheel.str() + "_Sec" + sector.str();
358 
359  string folder = topHistoFolder + "/Wheel" + wheel.str() + "/Sector" + sector.str();
360  dbe->setCurrentFolder(folder);
361 
362  if(sect!=4 && sect!=10) {
363  MeanHistos[make_pair(wh,sect)] =
364  dbe->book1D(MeanHistoName.c_str(),"Mean (from gaussian fit) of the residuals distribution",11,1,12);
365  } else {
366  MeanHistos[make_pair(wh,sect)] =
367  dbe->book1D(MeanHistoName.c_str(),"Mean (from gaussian fit) of the residuals distribution",13,1,14);
368  }
369  (MeanHistos[make_pair(wh,sect)])->setBinLabel(1,"MB1_SL1",1);
370  (MeanHistos[make_pair(wh,sect)])->setBinLabel(2,"MB1_SL2",1);
371  (MeanHistos[make_pair(wh,sect)])->setBinLabel(3,"MB1_SL3",1);
372  (MeanHistos[make_pair(wh,sect)])->setBinLabel(4,"MB2_SL1",1);
373  (MeanHistos[make_pair(wh,sect)])->setBinLabel(5,"MB2_SL2",1);
374  (MeanHistos[make_pair(wh,sect)])->setBinLabel(6,"MB2_SL3",1);
375  (MeanHistos[make_pair(wh,sect)])->setBinLabel(7,"MB3_SL1",1);
376  (MeanHistos[make_pair(wh,sect)])->setBinLabel(8,"MB3_SL2",1);
377  (MeanHistos[make_pair(wh,sect)])->setBinLabel(9,"MB3_SL3",1);
378  (MeanHistos[make_pair(wh,sect)])->setBinLabel(10,"MB4_SL1",1);
379  (MeanHistos[make_pair(wh,sect)])->setBinLabel(11,"MB4_SL3",1);
380  if(sect==4){
381  (MeanHistos[make_pair(wh,sect)])->setBinLabel(12,"MB4S13_SL1",1);
382  (MeanHistos[make_pair(wh,sect)])->setBinLabel(13,"MB4S13_SL3",1);
383  }
384  if(sect==10){
385  (MeanHistos[make_pair(wh,sect)])->setBinLabel(12,"MB4S14_SL1",1);
386  (MeanHistos[make_pair(wh,sect)])->setBinLabel(13,"MB4S14_SL3",1);
387  }
388 
389  if(sect!=4 && sect!=10) {
390  SigmaHistos[make_pair(wh,sect)] =
391  dbe->book1D(SigmaHistoName.c_str(),"Sigma (from gaussian fit) of the residuals distribution",11,1,12);
392  } else {
393  SigmaHistos[make_pair(wh,sect)] =
394  dbe->book1D(SigmaHistoName.c_str(),"Sigma (from gaussian fit) of the residuals distribution",13,1,14);
395  }
396  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(1,"MB1_SL1",1);
397  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(2,"MB1_SL2",1);
398  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(3,"MB1_SL3",1);
399  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(4,"MB2_SL1",1);
400  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(5,"MB2_SL2",1);
401  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(6,"MB2_SL3",1);
402  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(7,"MB3_SL1",1);
403  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(8,"MB3_SL2",1);
404  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(9,"MB3_SL3",1);
405  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(10,"MB4_SL1",1);
406  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(11,"MB4_SL3",1);
407  if(sect==4){
408  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(12,"MB4S13_SL1",1);
409  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(13,"MB4S13_SL3",1);
410  }
411  if(sect==10){
412  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(12,"MB4S14_SL1",1);
413  (SigmaHistos[make_pair(wh,sect)])->setBinLabel(13,"MB4S14_SL3",1);
414  }
415 
416 
417 }
418 
419 
420 string DTResolutionAnalysisTest::getMEName(const DTSuperLayerId & slID) {
421 
422  stringstream wheel; wheel << slID.wheel();
423  stringstream station; station << slID.station();
424  stringstream sector; sector << slID.sector();
425  stringstream superLayer; superLayer << slID.superlayer();
426 
427  string folderName =
428  topHistoFolder + "/Wheel" + wheel.str() +
429  "/Sector" + sector.str() +
430  "/Station" + station.str() + "/";
431 
432  if(doCalibAnalysis) folderName =
433  "DT/DTCalibValidation/Wheel" + wheel.str() +
434  "/Station" + station.str() + "/Sector" + sector.str() + "/";
435 
436  string histoname = folderName + "hResDist"
437  + "_W" + wheel.str()
438  + "_St" + station.str()
439  + "_Sec" + sector.str()
440  + "_SL" + superLayer.str();
441 
442  if(doCalibAnalysis) histoname = folderName + "hResDist_STEP3"
443  + "_W" + wheel.str()
444  + "_St" + station.str()
445  + "_Sec" + sector.str()
446  + "_SL" + superLayer.str();
447 
448  return histoname;
449 
450 }
451 
452 
453 
454 int DTResolutionAnalysisTest::stationFromBin(int bin) const {
455  return (int) (bin /3.1)+1;
456 }
457 
458 
459 int DTResolutionAnalysisTest::slFromBin(int bin) const {
460  int ret = bin%3;
461  if(ret == 0 || bin == 11) ret = 3;
462 
463  return ret;
464 }
465 
466 double DTResolutionAnalysisTest::meanInRange(double mean) const {
467  double value(0.);
468  if( fabs(mean) <= maxGoodMeanValue ) {value = 1.;}
469  else if( fabs(mean) > maxGoodMeanValue && fabs(mean) < minBadMeanValue ) {value = 0.9;}
470  else if( fabs(mean) >= minBadMeanValue ) {value = 0.1;}
471  return value;
472 }
473 
474 double DTResolutionAnalysisTest::sigmaInRange(double sigma) const {
475  double value(0.);
476  if( sigma <= maxGoodSigmaValue ) {value = 1.;}
477  else if( sigma > maxGoodSigmaValue && sigma < minBadSigmaValue ) {value = 0.9;}
478  else if( sigma >= minBadSigmaValue ) {value = 0.1;}
479  return value;
480 }
481 
482 void DTResolutionAnalysisTest::resetMEs() {
483  globalResSummary->Reset();
484  // Reset the summary histo
485  for(map<int, MonitorElement*> ::const_iterator histo = wheelMeanHistos.begin();
486  histo != wheelMeanHistos.end();
487  histo++) {
488  (*histo).second->Reset();
489  }
490  for(map<int, MonitorElement*> ::const_iterator histo = wheelSigmaHistos.begin();
491  histo != wheelSigmaHistos.end();
492  histo++) {
493  (*histo).second->Reset();
494  }
495 
496  for(int indx = -2; indx != 3; ++indx) {
497  meanDistr[indx]->Reset();
498  sigmaDistr[indx]->Reset();
499  }
500 }
501 
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
DTResolutionAnalysisTest::meanInRange
double meanInRange(double mean) const
Definition: DTResolutionAnalysisTest.cc:466
DTResolutionAnalysisTest::DTResolutionAnalysisTest
DTResolutionAnalysisTest(const edm::ParameterSet &ps)
Constructor.
Definition: DTResolutionAnalysisTest.cc:36
run_regression.ret
int ret
Definition: run_regression.py:388
DTResolutionAnalysisTest::endLuminosityBlock
void endLuminosityBlock(edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &c)
DQM Client Diagnostic.
Definition: DTResolutionAnalysisTest.cc:99
DTResolutionAnalysisTest::beginLuminosityBlock
void beginLuminosityBlock(edm::LuminosityBlock const &lumiSeg, edm::EventSetup const &context)
Definition: DTResolutionAnalysisTest.cc:80
DTResolutionAnalysisTest::getMEName
std::string getMEName(const DTSuperLayerId &slID)
Get the ME name.
Definition: DTResolutionAnalysisTest.cc:420
DTResolutionAnalysisTest::~DTResolutionAnalysisTest
virtual ~DTResolutionAnalysisTest()
Destructor.
Definition: DTResolutionAnalysisTest.cc:55
MonitorElement::getMean
double getMean(int axis=1) const
get mean value of histogram along x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:681
bookHistos
void bookHistos()
Definition: Histogram.h:33
HcalObjRepresent::Fill
void Fill(HcalDetId &id, double val, std::vector< TH2F > &depth)
Definition: HcalObjRepresent.h:611
funct::abs
Abs< T >::type abs(const T &t)
Definition: Abs.h:22
DTSuperLayerId::superLayer
int superLayer() const
Return the superlayer number.
Definition: DTSuperLayerId.h:49
edm::LuminosityBlockBase::run
RunNumber_t run() const
Definition: LuminosityBlockBase.h:49
nevents
int nevents
Definition: DTDataIntegrityTask.cc:31
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
dtDQMClient_cfg.prescaleFactor
tuple prescaleFactor
Definition: dtDQMClient_cfg.py:51
DTSuperLayerId::superlayer
int superlayer() const
Return the superlayer number (deprecated method name)
Definition: DTSuperLayerId.h:55
DTResolutionAnalysisTest::bookHistos
void bookHistos()
book the summary histograms
Definition: DTResolutionAnalysisTest.cc:104
edm::EventSetup::get
const T & get() const
Definition: EventSetup.h:55
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1498
MonitorElement::getRMS
double getRMS(int axis=1) const
get RMS of histogram along x, y or z axis (axis=1, 2, 3 respectively)
Definition: MonitorElement.cc:694
DTResolutionAnalysisTest::sigmaInRange
double sigmaInRange(double sigma) const
Definition: DTResolutionAnalysisTest.cc:474
DTChamberId::sector
int sector() const
Definition: DTChamberId.h:61
histoStyle.weight
int weight
Definition: histoStyle.py:50
DTChamberId::station
int station() const
Return the station number.
Definition: DTChamberId.h:51
DTChamberId::wheel
int wheel() const
Return the wheel number.
Definition: DTChamberId.h:45
DTResolutionAnalysisTest::beginRun
void beginRun(const edm::Run &r, const edm::EventSetup &c)
BeginRun.
Definition: DTResolutionAnalysisTest.cc:70
edm::Run
Definition: Run.h:41
DTResolutionAnalysisTest::analyze
void analyze(const edm::Event &e, const edm::EventSetup &c)
Analyze.
Definition: DTResolutionAnalysisTest.cc:91
DTResolutionAnalysisTest::endRun
void endRun(edm::Run const &run, edm::EventSetup const &c)
Definition: DTResolutionAnalysisTest.cc:158
Generated for CMSSW Reference Manual by doxygen 1.8.5