CMS 3D CMS Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Pages
EcalZmassClient.cc
Go to the documentation of this file.
1 // -*- C++ -*-
2 //
3 // Package: ZfitterAnalyzer
4 // Class: ZfitterAnalyzer
5 //
13 //
14 // Original Author: Vieri Candelise
15 // Created: Mon Jun 13 09:49:08 CEST 2011
16 //
17 //
18 
19 
20 // system include files
21 #include <memory>
22 
23 // user include files
24 #include "FWCore/Framework/interface/Frameworkfwd.h"
25 #include "FWCore/Framework/interface/EDAnalyzer.h"
26 #include "FWCore/ServiceRegistry/interface/Service.h"
27 #include "FWCore/Framework/interface/Event.h"
28 #include "FWCore/Framework/interface/MakerMacros.h"
29 #include "DQM/Physics/src/EwkDQM.h"
30 #include "FWCore/ParameterSet/interface/ParameterSet.h"
31 #include "TMath.h"
32 #include <string>
33 #include <cmath>
34 #include "TH1.h"
35 #include "DQMServices/Core/interface/DQMStore.h"
36 #include "DQMServices/Core/interface/MonitorElement.h"
37 #include "FWCore/Framework/interface/MakerMacros.h"
38 #include "FWCore/Framework/interface/EventSetup.h"
39 #include <iostream>
40 #include "FWCore/ParameterSet/interface/ParameterSet.h"
41 #include "FWCore/Utilities/interface/Exception.h"
42 
43 class DQMStore;
44 class MonitorElement;
45 
46 //
47 // class declaration
48 //
49 
50 class EcalZmassClient:public
51  edm::EDAnalyzer
52 {
53 public:
54  explicit
55  EcalZmassClient (const edm::ParameterSet &);
56  ~
57  EcalZmassClient ();
58 
59  static void
60  fillDescriptions (edm::ConfigurationDescriptions & descriptions);
61 
62 
63 private:
64  virtual void
65  beginJob () override;
66  virtual void
67  analyze (const edm::Event &, const edm::EventSetup &) override;
68  virtual void
69  endJob () override;
70 
71  virtual void
72  beginRun (edm::Run const &, edm::EventSetup const &) override;
73  virtual void
74  endRun (edm::Run const &, edm::EventSetup const &) override;
75  virtual void
76  beginLuminosityBlock (edm::LuminosityBlock const &,
77  edm::EventSetup const &) override;
78  virtual void
79  endLuminosityBlock (edm::LuminosityBlock const &, edm::EventSetup const &) override;
80 
81  // ----------member data ---------------------------
82 
83  std::string prefixME_;
84 
85  MonitorElement *
86  h_ee_invMass_EB;
87  MonitorElement *
88  h_ee_invMass_EE;
89  MonitorElement *
90  h_ee_invMass_BB;
91  MonitorElement *
92  h_95_ee_invMass_EB;
93  MonitorElement *
94  h_95_ee_invMass_EE;
95  MonitorElement *
96  h_95_ee_invMass_BB;
97  MonitorElement *
98  h_ee_invMass;
99  MonitorElement *
100  h_e1_et;
101  MonitorElement *
102  h_e2_et;
103  MonitorElement *
104  h_e1_eta;
105  MonitorElement *
106  h_e2_eta;
107  MonitorElement *
108  h_e1_phi;
109  MonitorElement *
110  h_e2_phi;
111  MonitorElement *
112  h_fitres1;
113  MonitorElement *
114  h_fitres1bis;
115  MonitorElement *
116  h_fitres1Chi2;
117  MonitorElement *
118  h_fitres2;
119  MonitorElement *
120  h_fitres2bis;
121  MonitorElement *
122  h_fitres2Chi2;
123  MonitorElement *
124  h_fitres3;
125  MonitorElement *
126  h_fitres3bis;
127  MonitorElement *
128  h_fitres3Chi2;
129 
130 
131 };
132 
133 EcalZmassClient::EcalZmassClient (const edm::ParameterSet & iConfig)
134 {
135  prefixME_ = iConfig.getUntrackedParameter < std::string > ("prefixME", "");
136 
137 }
138 
139 
140 EcalZmassClient::~EcalZmassClient ()
141 {
142 
143  // do anything here that needs to be done at desctruction time
144  // (e.g. close files, deallocate resources etc.)
145 
146 }
147 
148 
149 //
150 // member functions
151 //
152 
153 // ------------ method called for each event ------------
154 void
155 EcalZmassClient::analyze (const edm::Event & iEvent,
156  const edm::EventSetup & iSetup)
157 {
158 
159 }
160 
161 
162 // ------------ method called once each job just before starting event loop ------------
163 void
164 EcalZmassClient::beginJob ()
165 {
166 }
167 
168 // ------------ method called once each job just after ending the event loop ------------
169 void
170 EcalZmassClient::endJob ()
171 {
172 }
173 
174 // ------------ method called when starting to processes a run ------------
175 void
176 EcalZmassClient::beginRun (edm::Run const &, edm::EventSetup const &)
177 {
178  DQMStore *theDbe = edm::Service < DQMStore > ().operator-> ();
179 
180  theDbe->setCurrentFolder (prefixME_ + "/Zmass");
181  h_fitres1 =
182  theDbe->book1D ("Gaussian mean WP80 EB-EB",
183  "Gaussian mean WP80 EB-EB", 1, 0, 1);
184  h_fitres1bis =
185  theDbe->book1D ("Gaussian sigma WP80 EB-EB",
186  "Gaussian sigma WP80 EB-EB", 1, 0, 1);
187  h_fitres1Chi2 =
188  theDbe->book1D ("Gaussian Chi2 result over NDF WP80 EB-EB",
189  "Gaussian Chi2 result over NDF WP80 EB-EB", 1, 0, 1);
190 
191  h_fitres3 =
192  theDbe->book1D ("Gaussian mean WP80 EB-EE",
193  "Gaussian mean result WP80 EB-EE", 1, 0, 1);
194  h_fitres3bis =
195  theDbe->book1D ("Gaussian sigma WP80 EB-EE",
196  "Gaussian sigma WP80 EB-EE", 1, 0, 1);
197  h_fitres3Chi2 =
198  theDbe->book1D ("Gaussian Chi2 result over NDF WP80 EB-EE",
199  "Gaussian Chi2 result over NDF WP80 EB-EE", 1, 0, 1);
200 
201  h_fitres2 =
202  theDbe->book1D ("Gaussian mean WP80 EE-EE",
203  "Gaussian mean WP80 EE-EE", 1, 0, 1);
204  h_fitres2bis =
205  theDbe->book1D ("Gaussian sigma WP80 EE-EE",
206  "Gaussian sigma WP80 EE-EE", 1, 0, 1);
207  h_fitres2Chi2 =
208  theDbe->book1D ("Gaussian Chi2 result over NDF WP80 EE-EE",
209  "Gaussian Chi2 result over NDF WP80 EE-EE", 1, 0, 1);
210 }
211 
212 
213 //Breit-Wigner function
214 Double_t
215 mybw (Double_t * x, Double_t * par)
216 {
217  Double_t arg1 = 14.0 / 22.0; // 2 over pi
218  Double_t arg2 = par[1] * par[1] * par[2] * par[2]; //Gamma=par[2] M=par[1]
219  Double_t arg3 =
220  ((x[0] * x[0]) - (par[1] * par[1])) * ((x[0] * x[0]) - (par[1] * par[1]));
221  Double_t arg4 =
222  x[0] * x[0] * x[0] * x[0] * ((par[2] * par[2]) / (par[1] * par[1]));
223  return par[0] * arg1 * arg2 / (arg3 + arg4);
224 }
225 
226 //Gaussian
227 Double_t
228 mygauss (Double_t * x, Double_t * par)
229 {
230  Double_t arg = 0;
231  if (par[2] < 0)
232  par[2] = -par[2]; // par[2]: sigma
233  if (par[2] != 0)
234  arg = (x[0] - par[1]) / par[2]; // par[1]: mean
235  return par[0] * TMath::Exp (-0.5 * arg * arg) / (TMath::Sqrt (2 * TMath::Pi ()) * par[2]); // par[0] is constant
236 
237 }
238 
239 
240 
241 // ------------ method called when ending the processing of a run ------------
242 void
243 EcalZmassClient::endRun (edm::Run const &, edm::EventSetup const &)
244 {
245  DQMStore *theDbe = edm::Service < DQMStore > ().operator-> ();
246 
247  LogTrace ("EwkAnalyzer") << "Parameters initialization";
248 
249  MonitorElement *me1 = theDbe->get (prefixME_ + "/Zmass/Z peak - WP80 EB-EB");
250  MonitorElement *me2 = theDbe->get (prefixME_ + "/Zmass/Z peak - WP80 EE-EE");
251  MonitorElement *me3 = theDbe->get (prefixME_ + "/Zmass/Z peak - WP80 EB-EE");
252 
253 
254  if (me1 != 0)
255  {
256  TH1F *B = me1->getTH1F ();
257  TH1F *R1 = h_fitres1->getTH1F ();
258  int division = B->GetNbinsX ();
259  float massMIN = B->GetBinLowEdge (1);
260  float massMAX = B->GetBinLowEdge (division + 1);
261  //float BIN_SIZE = B->GetBinWidth(1);
262 
263  TF1 *func = new TF1 ("mygauss", mygauss, massMIN, massMAX, 3);
264  func->SetParameter (0, 1.0);
265  func->SetParName (0, "const");
266  func->SetParameter (1, 95.0);
267  func->SetParName (1, "mean");
268  func->SetParameter (2, 5.0);
269  func->SetParName (2, "sigma");
270 
271  double stats[4];
272  R1->GetStats (stats);
273  float N = 0;
274  float mean = 0;
275  float sigma = 0;
276  N = B->GetEntries ();
277 
278  try
279  {
280  if (N != 0)
281  {
282 
283  B->Fit ("mygauss", "QR");
284  mean = fabs (func->GetParameter (1));
285  sigma = fabs (func->GetParError (1));
286  }
287 
288  if (N == 0 || mean < 50 || mean > 100 || sigma <= 0 || sigma > 20)
289  {
290  N = 1;
291  mean = 0;
292  sigma = 0;
293  }
294 
295  }
296  catch (cms::Exception& e)
297  {
298  edm::LogError ("ZFitter") << "[Zfitter]: Exception when fitting..." << e.what();
299  N = 1;
300  mean = 40;
301  sigma = 0;
302  }
303 
304  stats[0] = N;
305  stats[1] = N;
306  stats[2] = mean * N;
307  stats[3] = sigma * sigma * N + mean * mean * N;
308 
309  R1->SetEntries (N);
310  R1->PutStats (stats);
311  }
312 /*******************************************************/
313  if (me1 != 0)
314  {
315  TH1F *Bbis = me1->getTH1F ();
316  TH1F *R1bis = h_fitres1bis->getTH1F ();
317  int division = Bbis->GetNbinsX ();
318  float massMIN = Bbis->GetBinLowEdge (1);
319  float massMAX = Bbis->GetBinLowEdge (division + 1);
320  //float BIN_SIZE = B->GetBinWidth(1);
321 
322  TF1 *func = new TF1 ("mygauss", mygauss, massMIN, massMAX, 3);
323  func->SetParameter (0, 1.0);
324  func->SetParName (0, "const");
325  func->SetParameter (1, 95.0);
326  func->SetParName (1, "mean");
327  func->SetParameter (2, 5.0);
328  func->SetParName (2, "sigma");
329 
330  double stats[4];
331  R1bis->GetStats (stats);
332  float N = 0;
333  float rms = 0;
334  float rmsErr = 0;
335  N = Bbis->GetEntries ();
336 
337  try
338  {
339  if (N != 0)
340  {
341 
342  Bbis->Fit ("mygauss", "QR");
343  rms = fabs (func->GetParameter (2));
344  rmsErr = fabs (func->GetParError (2));
345  }
346 
347  if (N == 0 || rms < 0 || rms > 50 || rmsErr <= 0 || rmsErr > 50)
348  {
349  N = 1;
350  rms = 0;
351  rmsErr = 0;
352  }
353 
354  }
355  catch (cms::Exception& e)
356  {
357  edm::LogError ("ZFitter") << "[Zfitter]: Exception when fitting..." << e.what();
358  N = 1;
359  rms = 40;
360  rmsErr = 0;
361  }
362 
363  stats[0] = N;
364  stats[1] = N;
365  stats[2] = rms * N;
366  stats[3] = rmsErr * rmsErr * N + rms * rms * N;
367 
368  R1bis->SetEntries (N);
369  R1bis->PutStats (stats);
370  }
371 /****************************************/
372 
373  if (me2 != 0)
374  {
375  TH1F *E = me2->getTH1F ();
376  TH1F *R2 = h_fitres2->getTH1F ();
377  int division = E->GetNbinsX ();
378  float massMIN = E->GetBinLowEdge (1);
379  float massMAX = E->GetBinLowEdge (division + 1);
380  //float BIN_SIZE = E->GetBinWidth(1);
381 
382  TF1 *func = new TF1 ("mygauss", mygauss, massMIN, massMAX, 3);
383  func->SetParameter (0, 1.0);
384  func->SetParName (0, "const");
385  func->SetParameter (1, 95.0);
386  func->SetParName (1, "mean");
387  func->SetParameter (2, 5.0);
388  func->SetParName (2, "sigma");
389 
390  double stats[4];
391  R2->GetStats (stats);
392  float N = 0;
393  float mean = 0;
394  float sigma = 0;
395  N = E->GetEntries ();
396 
397  try
398  {
399  if (N != 0)
400  {
401  E->Fit ("mygauss", "QR");
402  mean = fabs (func->GetParameter (1));
403  sigma = fabs (func->GetParError (1));
404  }
405 
406  if (N == 0 || mean < 50 || mean > 100 || sigma <= 0 || sigma > 20)
407  {
408  N = 1;
409  mean = 0;
410  sigma = 0;
411 
412  }
413 
414  }
415  catch (cms::Exception& e)
416  {
417  edm::LogError ("ZFitter") << "[Zfitter]: Exception when fitting..." << e.what();
418  N = 1;
419  mean = 40;
420  sigma = 0;
421  }
422 
423  stats[0] = N;
424  stats[1] = N;
425  stats[2] = mean * N;
426  stats[3] = sigma * sigma * N + mean * mean * N;
427 
428  R2->SetEntries (N);
429  R2->PutStats (stats);
430  }
431 /**************************************************************************/
432 
433  if (me2 != 0)
434  {
435  TH1F *Ebis = me2->getTH1F ();
436  TH1F *R2bis = h_fitres2bis->getTH1F ();
437  int division = Ebis->GetNbinsX ();
438  float massMIN = Ebis->GetBinLowEdge (1);
439  float massMAX = Ebis->GetBinLowEdge (division + 1);
440  //float BIN_SIZE = B->GetBinWidth(1);
441 
442  TF1 *func = new TF1 ("mygauss", mygauss, massMIN, massMAX, 3);
443  func->SetParameter (0, 1.0);
444  func->SetParName (0, "const");
445  func->SetParameter (1, 95.0);
446  func->SetParName (1, "mean");
447  func->SetParameter (2, 5.0);
448  func->SetParName (2, "sigma");
449 
450  double stats[4];
451  R2bis->GetStats (stats);
452  float N = 0;
453  float rms = 0;
454  float rmsErr = 0;
455  N = Ebis->GetEntries ();
456 
457  try
458  {
459  if (N != 0)
460  {
461 
462  Ebis->Fit ("mygauss", "QR");
463  rms = fabs (func->GetParameter (2));
464  rmsErr = fabs (func->GetParError (2));
465  }
466 
467  if (N == 0 || rms < 0 || rms > 50 || rmsErr <= 0 || rmsErr > 50)
468  {
469  N = 1;
470  rms = 0;
471  rmsErr = 0;
472  }
473 
474  }
475  catch (cms::Exception& e)
476  {
477  edm::LogError ("ZFitter") << "[Zfitter]: Exception when fitting..." << e.what();
478  N = 1;
479  rms = 40;
480  rmsErr = 0;
481  }
482 
483  stats[0] = N;
484  stats[1] = N;
485  stats[2] = rms * N;
486  stats[3] = rmsErr * rmsErr * N + rms * rms * N;
487 
488  R2bis->SetEntries (N);
489  R2bis->PutStats (stats);
490  }
491 /*********************************************************************************************/
492 
493  if (me3 != 0)
494  {
495  TH1F *R3 = h_fitres3->getTH1F ();
496  TH1F *M = me3->getTH1F ();
497  int division = M->GetNbinsX ();
498  float massMIN = M->GetBinLowEdge (1);
499  float massMAX = M->GetBinLowEdge (division + 1);
500  //float BIN_SIZE = M->GetBinWidth(1);
501 
502  TF1 *func = new TF1 ("mygauss", mygauss, massMIN, massMAX, 3);
503  func->SetParameter (0, 1.0);
504  func->SetParName (0, "const");
505  func->SetParameter (1, 95.0);
506  func->SetParName (1, "mean");
507  func->SetParameter (2, 5.0);
508  func->SetParName (2, "sigma");
509 
510  double stats[4];
511  R3->GetStats (stats);
512  float N = 0;
513  float mean = 0;
514  float sigma = 0;
515  N = M->GetEntries ();
516 
517 
518  try
519  {
520  if (N != 0)
521  {
522 
523  M->Fit ("mygauss", "QR");
524  mean = fabs (func->GetParameter (1));
525  sigma = fabs (func->GetParError (1));
526  }
527  if (N == 0 || mean < 50 || mean > 100 || sigma <= 0 || sigma > 20)
528  {
529  N = 1;
530  mean = 0;
531  sigma = 0;
532  }
533 
534  }
535  catch (cms::Exception& e)
536  {
537  edm::LogError ("ZFitter") << "[Zfitter]: Exception when fitting..." << e.what();
538  N = 1;
539  mean = 40;
540  sigma = 0;
541  }
542 
543  stats[0] = N;
544  stats[1] = N;
545  stats[2] = mean * N;
546  stats[3] = sigma * sigma * N + mean * mean * N;
547 
548  R3->SetEntries (N);
549  R3->PutStats (stats);
550  }
551  /********************************************************************************/
552 
553  if (me3 != 0)
554  {
555  TH1F *Mbis = me3->getTH1F ();
556  TH1F *R3bis = h_fitres3bis->getTH1F ();
557  int division = Mbis->GetNbinsX ();
558  float massMIN = Mbis->GetBinLowEdge (1);
559  float massMAX = Mbis->GetBinLowEdge (division + 1);
560  //float BIN_SIZE = B->GetBinWidth(1);
561 
562  TF1 *func = new TF1 ("mygauss", mygauss, massMIN, massMAX, 3);
563  func->SetParameter (0, 1.0);
564  func->SetParName (0, "const");
565  func->SetParameter (1, 95.0);
566  func->SetParName (1, "mean");
567  func->SetParameter (2, 5.0);
568  func->SetParName (2, "sigma");
569 
570  double stats[4];
571  R3bis->GetStats (stats);
572  float N = 0;
573  float rms = 0;
574  float rmsErr = 0;
575  N = Mbis->GetEntries ();
576 
577  try
578  {
579  if (N != 0)
580  {
581 
582  Mbis->Fit ("mygauss", "QR");
583  rms = fabs (func->GetParameter (2));
584  rmsErr = fabs (func->GetParError (2));
585  }
586 
587  if (N == 0 || rms < 0 || rms > 50 || rmsErr <= 0 || rmsErr > 50)
588  {
589  N = 1;
590  rms = 0;
591  rmsErr = 0;
592  }
593 
594  }
595  catch (cms::Exception& e)
596  {
597  edm::LogError ("ZFitter") << "[Zfitter]: Exception when fitting..." << e.what();
598  N = 1;
599  rms = 40;
600  rmsErr = 0;
601  }
602 
603  stats[0] = N;
604  stats[1] = N;
605  stats[2] = rms * N;
606  stats[3] = rmsErr * rmsErr * N + rms * rms * N;
607 
608  R3bis->SetEntries (N);
609  R3bis->PutStats (stats);
610  }
611 
612  /*Chi2 */
613 
614  if (me1 != 0)
615  {
616  TH1F *C1 = me1->getTH1F ();
617  TH1F *S1 = h_fitres1Chi2->getTH1F ();
618  int division = C1->GetNbinsX ();
619  float massMIN = C1->GetBinLowEdge (1);
620  float massMAX = C1->GetBinLowEdge (division + 1);
621  //float BIN_SIZE = B->GetBinWidth(1);
622 
623  TF1 *func = new TF1 ("mygauss", mygauss, massMIN, massMAX, 3);
624  func->SetParameter (0, 1.0);
625  func->SetParName (0, "const");
626  func->SetParameter (1, 95.0);
627  func->SetParName (1, "mean");
628  func->SetParameter (2, 5.0);
629  func->SetParName (2, "sigma");
630 
631  double stats[4];
632  S1->GetStats (stats);
633  float N = 0;
634  float Chi2 = 0;
635  float NDF = 0;
636  N = C1->GetEntries ();
637 
638  try
639  {
640  if (N != 0)
641  {
642 
643  C1->Fit ("mygauss", "QR");
644  if ((func->GetNDF () != 0))
645  {
646  Chi2 = fabs (func->GetChisquare ()) / fabs (func->GetNDF ());
647  NDF = 0.1;
648  }
649  }
650 
651  if (N == 0 || Chi2 < 0 || NDF < 0)
652  {
653  N = 1;
654  Chi2 = 0;
655  NDF = 0;
656  }
657 
658  }
659  catch (cms::Exception& e)
660  {
661  edm::LogError ("ZFitter") << "[Zfitter]: Exception when fitting..." << e.what();
662  N = 1;
663  Chi2 = 40;
664  NDF = 0;
665  }
666 
667  stats[0] = N;
668  stats[1] = N;
669  stats[2] = Chi2 * N;
670  stats[3] = NDF * NDF * N + Chi2 * Chi2 * N;
671 
672  S1->SetEntries (N);
673  S1->PutStats (stats);
674  }
675  /**********************************************/
676 
677  if (me2 != 0)
678  {
679  TH1F *C2 = me2->getTH1F ();
680  TH1F *S2 = h_fitres2Chi2->getTH1F ();
681  int division = C2->GetNbinsX ();
682  float massMIN = C2->GetBinLowEdge (1);
683  float massMAX = C2->GetBinLowEdge (division + 1);
684  //float BIN_SIZE = B->GetBinWidth(1);
685 
686  TF1 *func = new TF1 ("mygauss", mygauss, massMIN, massMAX, 3);
687  func->SetParameter (0, 1.0);
688  func->SetParName (0, "const");
689  func->SetParameter (1, 95.0);
690  func->SetParName (1, "mean");
691  func->SetParameter (2, 5.0);
692  func->SetParName (2, "sigma");
693 
694  double stats[4];
695  S2->GetStats (stats);
696  float N = 0;
697  float Chi2 = 0;
698  float NDF = 0;
699  N = C2->GetEntries ();
700 
701  try
702  {
703  if (N != 0)
704  {
705  C2->Fit ("mygauss", "QR");
706  if (func->GetNDF () != 0)
707  {
708  Chi2 = fabs (func->GetChisquare ()) / fabs (func->GetNDF ());
709  NDF = 0.1;
710  }
711  }
712 
713  if (N == 0 || Chi2 < 0 || NDF < 0)
714  {
715  N = 1;
716  Chi2 = 0;
717  NDF = 0;
718  }
719 
720  }
721  catch (cms::Exception& e)
722  {
723  edm::LogError ("ZFitter") << "[Zfitter]: Exception when fitting..." << e.what();
724  N = 1;
725  Chi2 = 40;
726  NDF = 0;
727  }
728 
729  stats[0] = N;
730  stats[1] = N;
731  stats[2] = Chi2 * N;
732  stats[3] = NDF * NDF * N + Chi2 * Chi2 * N;
733 
734  S2->SetEntries (N);
735  S2->PutStats (stats);
736  }
737  /**************************************************************************/
738  if (me3 != 0)
739  {
740  TH1F *C3 = me3->getTH1F ();
741  TH1F *S3 = h_fitres3Chi2->getTH1F ();
742  int division = C3->GetNbinsX ();
743  float massMIN = C3->GetBinLowEdge (1);
744  float massMAX = C3->GetBinLowEdge (division + 1);
745  //float BIN_SIZE = B->GetBinWidth(1);
746 
747  TF1 *func = new TF1 ("mygauss", mygauss, massMIN, massMAX, 3);
748  func->SetParameter (0, 1.0);
749  func->SetParName (0, "const");
750  func->SetParameter (1, 95.0);
751  func->SetParName (1, "mean");
752  func->SetParameter (2, 5.0);
753  func->SetParName (2, "sigma");
754 
755  double stats[4];
756  S3->GetStats (stats);
757  float N = 0;
758  float Chi2 = 0;
759  float NDF = 0;
760  N = C3->GetEntries ();
761 
762  try
763  {
764  if (N != 0)
765  {
766  C3->Fit ("mygauss", "QR");
767  if ((func->GetNDF () != 0))
768  {
769  Chi2 = fabs (func->GetChisquare ()) / fabs (func->GetNDF ());
770  NDF = 0.1;
771  }
772  }
773 
774 
775  if (N == 0 || Chi2 < 0 || NDF < 0)
776  {
777  N = 1;
778  Chi2 = 0;
779  NDF = 0;
780  }
781 
782  }
783  catch (cms::Exception& e)
784  {
785  edm::LogError ("ZFitter") << "[Zfitter]: Exception when fitting..." << e.what();
786  N = 1;
787  Chi2 = 40;
788  NDF = 0;
789  }
790 
791  stats[0] = N;
792  stats[1] = N;
793  stats[2] = Chi2 * N;
794  stats[3] = NDF * NDF * N + Chi2 * Chi2 * N;
795 
796  S3->SetEntries (N);
797  S3->PutStats (stats);
798  }
799 }
800 
801  // ------------ method called when starting to processes a luminosity block ------------
802 void
803 EcalZmassClient::beginLuminosityBlock (edm::LuminosityBlock const &,
804  edm::EventSetup const &)
805 {
806 }
807 
808  // ------------ method called when ending the processing of a luminosity block ------------
809 void
810 EcalZmassClient::endLuminosityBlock (edm::LuminosityBlock const &,
811  edm::EventSetup const &)
812 {
813 }
814 
815  // ------------ method fills 'descriptions' with the allowed parameters for the module ------------
816 void
817 EcalZmassClient::fillDescriptions (edm::
818  ConfigurationDescriptions & descriptions)
819 {
820  //The following says we do not know what parameters are allowed so do no validation
821  // Please change this to state exactly what you do use, even if it is no parameters
822  edm::ParameterSetDescription desc;
823  desc.setUnknown ();
824  descriptions.addDefault (desc);
825 }
826 
827  //define this as a plug-in
828 DEFINE_FWK_MODULE (EcalZmassClient);
EcalZmassClient::h_e1_phi
MonitorElement * h_e1_phi
Definition: EcalZmassClient.cc:108
cms::Exception::what
virtual char const * what() const
Definition: Exception.cc:141
Pi
const double Pi
Definition: CosmicMuonParameters.h:18
EcalZmassClient::h_e2_et
MonitorElement * h_e2_et
Definition: EcalZmassClient.cc:102
edm::ParameterSet::getUntrackedParameter
T getUntrackedParameter(std::string const &, T const &) const
EcalZmassClient::endRun
virtual void endRun(edm::Run const &, edm::EventSetup const &) override
Definition: EcalZmassClient.cc:243
EcalZmassClient::h_fitres1Chi2
MonitorElement * h_fitres1Chi2
Definition: EcalZmassClient.cc:116
EcalZmassClient::h_fitres2
MonitorElement * h_fitres2
Definition: EcalZmassClient.cc:118
EcalZmassClient::h_ee_invMass_EE
MonitorElement * h_ee_invMass_EE
Definition: EcalZmassClient.cc:88
EcalZmassClient::beginLuminosityBlock
virtual void beginLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
Definition: EcalZmassClient.cc:803
EcalZmassClient::h_ee_invMass
MonitorElement * h_ee_invMass
Definition: EcalZmassClient.cc:98
EcalZmassClient::h_ee_invMass_EB
MonitorElement * h_ee_invMass_EB
Definition: EcalZmassClient.cc:86
DQMStore::book1D
MonitorElement * book1D(const char *name, const char *title, int nchX, double lowX, double highX)
Book 1D histogram.
Definition: DQMStore.cc:873
EcalZmassClient::h_95_ee_invMass_EE
MonitorElement * h_95_ee_invMass_EE
Definition: EcalZmassClient.cc:94
DEFINE_FWK_MODULE
#define DEFINE_FWK_MODULE(type)
Definition: MakerMacros.h:17
EcalZmassClient::h_e2_phi
MonitorElement * h_e2_phi
Definition: EcalZmassClient.cc:110
funct::arg
A arg
Definition: Factorize.h:36
mybw
Double_t mybw(Double_t *x, Double_t *par)
Definition: EcalZmassClient.cc:215
EcalZmassClient::h_e1_eta
MonitorElement * h_e1_eta
Definition: EcalZmassClient.cc:104
AlCaHLTBitMon_QueryRunRegistry.string
string string
Definition: AlCaHLTBitMon_QueryRunRegistry.py:255
EcalZmassClient::endJob
virtual void endJob() override
Definition: EcalZmassClient.cc:170
EcalZmassClient::h_fitres1
MonitorElement * h_fitres1
Definition: EcalZmassClient.cc:112
EcalZmassClient::beginJob
virtual void beginJob() override
Definition: EcalZmassClient.cc:164
double_binary
Definition: DDStreamer.cc:227
EcalZmassClient::h_e1_et
MonitorElement * h_e1_et
Definition: EcalZmassClient.cc:100
mygauss
Double_t mygauss(Double_t *x, Double_t *par)
Definition: EcalZmassClient.cc:228
iEvent
int iEvent
Definition: GenABIO.cc:243
EcalZmassClient::h_ee_invMass_BB
MonitorElement * h_ee_invMass_BB
Definition: EcalZmassClient.cc:90
EcalZmassClient::h_95_ee_invMass_BB
MonitorElement * h_95_ee_invMass_BB
Definition: EcalZmassClient.cc:96
EcalZmassClient::h_fitres2Chi2
MonitorElement * h_fitres2Chi2
Definition: EcalZmassClient.cc:122
DQMStore::get
MonitorElement * get(const std::string &path) const
get ME from full pathname (e.g. &quot;my/long/dir/my_histo&quot;)
Definition: DQMStore.cc:1624
LogTrace
#define LogTrace(id)
Definition: MessageLogger.h:502
EcalZmassClient::h_95_ee_invMass_EB
MonitorElement * h_95_ee_invMass_EB
Definition: EcalZmassClient.cc:92
EcalZmassClient::h_fitres3Chi2
MonitorElement * h_fitres3Chi2
Definition: EcalZmassClient.cc:128
EcalZmassClient::h_fitres2bis
MonitorElement * h_fitres2bis
Definition: EcalZmassClient.cc:120
N
#define N
Definition: blowfish.cc:9
EcalZmassClient::h_e2_eta
MonitorElement * h_e2_eta
Definition: EcalZmassClient.cc:106
EcalZmassClient::h_fitres3
MonitorElement * h_fitres3
Definition: EcalZmassClient.cc:124
EcalZmassClient::h_fitres3bis
MonitorElement * h_fitres3bis
Definition: EcalZmassClient.cc:126
MonitorElement::getTH1F
TH1F * getTH1F(void) const
Definition: MonitorElement.cc:1498
EcalZmassClient::analyze
virtual void analyze(const edm::Event &, const edm::EventSetup &) override
Definition: EcalZmassClient.cc:155
EcalZmassClient::beginRun
virtual void beginRun(edm::Run const &, edm::EventSetup const &) override
Definition: EcalZmassClient.cc:176
statics.func
string func
Definition: statics.py:48
EcalZmassClient::fillDescriptions
static void fillDescriptions(edm::ConfigurationDescriptions &descriptions)
Definition: EcalZmassClient.cc:817
x
Definition: DDAxes.h:10
EcalZmassClient::endLuminosityBlock
virtual void endLuminosityBlock(edm::LuminosityBlock const &, edm::EventSetup const &) override
Definition: EcalZmassClient.cc:810
Chi2
Definition: Chi2.h:17
EcalZmassClient::EcalZmassClient
EcalZmassClient(const edm::ParameterSet &)
Definition: EcalZmassClient.cc:133
DQMStore::setCurrentFolder
void setCurrentFolder(const std::string &fullpath)
Definition: DQMStore.cc:585
EcalZmassClient::h_fitres1bis
MonitorElement * h_fitres1bis
Definition: EcalZmassClient.cc:114
edm::Run
Definition: Run.h:41
EcalZmassClient::prefixME_
std::string prefixME_
Definition: EcalZmassClient.cc:83
Generated for CMSSW Reference Manual by doxygen 1.8.5