ATLAS Offline Software
T0MTHistos.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 #include "MdtCalibT0/T0MTHistos.h"
6 
7 #include <cmath>
8 
9 #include "AthenaKernel/getMessageSvc.h"
10 #include "CxxUtils/checker_macros.h"
11 #include "GaudiKernel/MsgStream.h"
12 #include "MdtCalibT0/MTT0PatternRecognition.h"
13 #include "MdtCalibT0/MTTmaxPatternRecognition.h"
14 #include "TDirectory.h"
15 #include "TLine.h"
16 #include "TRandom3.h"
17 #include "list"
18 #include "boost/thread/tss.hpp"
19 
20 namespace MuonCalib {
21 
22  TRandom3* getTLSRandomGen()
23  {
24  static boost::thread_specific_ptr<TRandom3> rnd ATLAS_THREAD_SAFE;
25  TRandom3* random = rnd.get();
26  if (!random) {
27  random = new TRandom3();
28  rnd.reset(random);
29  }
30  return random;
31  }
32 
34  inline Double_t mt_t0_fermi(Double_t *x, Double_t *par) {
35  // more convenient parameters
36  const Double_t &t(x[0]);
37  const Double_t &t_0(par[T0MTHistos::T0_PAR_NR_T0]), &T(par[T0MTHistos::T0_PAR_NR_T]), &back(par[T0MTHistos::T0_PAR_NR_BACK]),
38  &A(par[T0MTHistos::T0_PAR_NR_A]);
39  // the formula
40  return (back + A / (1 + std::exp(-(t - t_0) / T)));
41  }
42 
44  inline Double_t mt_tmax_fermi(Double_t *x, Double_t *par) {
45  // more convenient parameters
46  Double_t &t(x[0]);
47  Double_t &t_max(par[T0MTHistos::TMAX_PAR_NR_TMAX]), &T(par[T0MTHistos::TMAX_PAR_NR_T]), &back(par[T0MTHistos::TMAX_PAR_NR_BACK]),
48  &a(par[T0MTHistos::TMAX_PAR_NR_A]), &b(par[T0MTHistos::TMAX_PAR_NR_B]), &t_0(par[T0MTHistos::TMAX_PAR_NR_T0]);
49  // the formula
50  return (back + (std::exp(a + b * (t - t_0))) / (1 + std::exp((t - t_max) / T)));
51  }
52 
54  // Fill.. //
56  void T0MTHistos::FillT(double t) { m_time->Fill(static_cast<Axis_t>(t)); }
57 
59  // Initialize(int id, const T0MTSettings & settings) //
61  void T0MTHistos::Initialize(int id, const T0MTSettings *settings, const char *hname) {
62  m_settings = settings;
63 #ifndef NDEBUG
64  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
65  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "T0MTHistos::Initialize: called" << endmsg;
66 #endif
67  char buf[100];
68  if (hname == nullptr)
69  snprintf(buf, 100, "t_spec_%d", id);
70  else
71  snprintf(buf, 100, "t_spec_%s", hname);
72 #ifndef NDEBUG
73  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "directory=" << gDirectory->GetName() << endmsg;
74 #endif
75  m_time = std::make_unique<TH1F>(buf, "", settings->NBinsTime(), settings->TimeMin(), settings->TimeMax());
76  m_id = id;
77  if (settings->DrawDebugGraphs()) {
78 #ifndef NDEBUG
79  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "T0MTHistos::Initialize: debug directory created" << endmsg;
80 #endif
81  TDirectory *cwd = gDirectory;
82  snprintf(buf, 100, "t0_tmax_dir_%d", id);
83  m_dir = gDirectory->mkdir(buf, buf);
84  cwd->cd();
85  } else {
86  m_dir = nullptr;
87  }
88  m_t0_ok = false;
89  m_tmax_ok = false;
90  } // end T0MTHistos::Initialize
91 
93  // SetTSpec(int id, TH1F *spec) //
95  void T0MTHistos::SetTSpec(int id, TH1F *spec, const T0MTSettings *settings, bool copy_spec) {
96  m_settings = settings;
97  TDirectory *cwd = gDirectory;
98  if (copy_spec)
99  m_time = std::make_unique<TH1F>(*spec);
100  else
101  m_time.reset(spec);
102  m_id = id;
103  if (settings->DrawDebugGraphs()) {
104  char buffer[100];
105  snprintf(buffer, 100, "t0_tmax_dir_%d", id);
106  m_dir = gDirectory->mkdir(buffer, buffer);
107  cwd->cd();
108  }
109  } // end T0MTHistos::SetTSpec
110 
112  // FitT0() //
114  bool T0MTHistos::FitT0() {
115  if (m_time->GetEntries() < 1000) {
116  m_status_code = 2;
117  return false;
118  }
119  if (!NormalFit()) return false;
120  if (m_time->GetEntries() < 10000) {
121  m_status_code = 0;
122  return true;
123  }
124  if (m_settings->T0Settings()->ScrambleThreshold() > 0 && m_chi2 > m_settings->T0Settings()->ScrambleThreshold()) {
125  if (!T0Scramble()) {
126  m_status_code = 3;
127  return false;
128  }
129  }
130  if (m_settings->T0Settings()->SlicingThreshold() > 0 && m_chi2 > m_settings->T0Settings()->SlicingThreshold()) { TopSlicing(); }
131 #ifndef NDEBUG
132  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
133  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << m_time->GetName() << " " << m_chi2 << endmsg;
134 #endif
135  m_status_code = 0;
136  return true;
137  } // end T0MTHistos::FitT0
138 
140  // FitTmax() //
142  bool T0MTHistos::FitTmax() {
143  TDirectory *cwd = gDirectory;
144  if (m_dir != nullptr) m_dir->cd();
145  if (!m_time) {
146  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
147  log << MSG::WARNING << "T0MTHistos::FitTmax: Class is not initialized!" << endmsg;
148  m_tmax_ok = false;
149  cwd->cd();
150  return false;
151  }
152  // check if t0-fit was successfull t0 is needed for tmax-pattern recognition
153  if (!m_t0_ok || m_t0_fermi == nullptr) {
154  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
155  log << MSG::WARNING << "T0MTHistos::FitTmax for tube " << m_id << ": No valid t0-value!" << endmsg;
156  m_tmax_ok = false;
157  cwd->cd();
158  return false;
159  }
160  // Create pattern Recognition Class
161  MTTmaxPatternRecognition rec;
162  // perform pattern recognition
163  if (!rec.Initialize(m_time.get(), m_t0_fermi->GetParameter(T0_PAR_NR_T0), m_settings)) {
164  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
165  log << MSG::WARNING << "T0MTHistos::FitTmax for tube " << m_id << ": Pattern recognition failed!" << endmsg;
166  m_tmax_ok = false;
167  cwd->cd();
168  return false;
169  }
170  // create function object
171  char buffer[100];
172  snprintf(buffer, 100, "mt_tmax_fermi");
173  if (!m_tmax_fermi) {
174  m_tmax_fermi = std::make_unique<TF1>(buffer, mt_tmax_fermi, rec.GetFitRangeMin(), rec.GetFitRangeMax(), N_TMAX_FIT_PAR);
175  // set parameter names
176  m_tmax_fermi->SetParName(TMAX_PAR_NR_TMAX, "t_{max}");
177  m_tmax_fermi->SetParName(TMAX_PAR_NR_T, "T");
178  m_tmax_fermi->SetParName(TMAX_PAR_NR_BACK, "r_{b}");
179  m_tmax_fermi->SetParName(TMAX_PAR_NR_A, "a");
180  m_tmax_fermi->SetParName(TMAX_PAR_NR_B, "b");
181  // set fixed values
182  m_tmax_fermi->FixParameter(TMAX_PAR_NR_BACK, rec.GetBackground());
183  m_tmax_fermi->FixParameter(TMAX_PAR_NR_A, rec.GetA());
184  m_tmax_fermi->FixParameter(TMAX_PAR_NR_B, rec.GetB());
185  m_tmax_fermi->FixParameter(TMAX_PAR_NR_T0, m_t0_fermi->GetParameter(T0_PAR_NR_T0));
186  // set start values
187  m_tmax_fermi->SetParameter(TMAX_PAR_NR_TMAX, rec.GetEstimatedTMax());
188  m_tmax_fermi->SetParameter(TMAX_PAR_NR_T, 3.0);
189  }
190  // perform fit
191  if (m_dir != nullptr) {
192  m_tmax_fermi->SetLineColor(3);
193  m_tmax_fermi->Write();
194  }
195  m_tmax_fermi->SetLineColor(4);
196  std::string fitopt("LR");
197  if (m_settings->VerboseLevel() == 0) fitopt += "Q";
198  if (m_settings->AddFitfun()) {
199  fitopt += "+";
200  } else {
201  fitopt += "0";
202  }
203  m_time->Fit(m_tmax_fermi.get(), fitopt.c_str());
204  m_tmax_ok = true;
205  cwd->cd();
206  return true;
207  } // end T0MTHistos::FitTmax
208 
210  // FitT0() //
212  bool T0MTHistos::NormalFit() {
213 #ifndef NDEBUG
214  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
215  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "T0MTHistos::FitT0(): called" << endmsg;
216 #endif
217  TDirectory *cwd = gDirectory;
218  if (m_dir) m_dir->cd();
219  // check if class is initialized
220  if (!m_time) {
221 #ifdef NDEBUG
222  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
223 #endif
224  log << MSG::WARNING << "T0MTHistos::FitT0: Class is not initialized!" << endmsg;
225  m_t0_ok = false;
226  cwd->cd();
227  m_status_code = 3;
228  return false;
229  }
230  // create pattern recognition class
231  MTT0PatternRecognition rec;
232  // perform pattern recognition
233  if (!rec.Initialize(m_time.get(), m_settings)) {
234  m_t0_ok = false;
235 #ifdef NDEBUG
236  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
237 #endif
238  log << MSG::WARNING << "T0MTHistos::FitT0 for tube " << m_id << ": Pattern recognition failed!" << endmsg;
239  cwd->cd();
240  m_status_code = 3;
241  return false;
242  }
243  m_t0_ok = true;
244  // create function class
245  char buffer[100];
246  // sprintf(buffer,"mt_t0_fermi_%d", m_id);
247  snprintf(buffer, 100, "mt_t0_fermi");
248  m_t0_fermi = std::make_unique<TF1>(buffer, mt_t0_fermi, rec.GetFitRangeMin(), rec.GetFitRangeMax(), N_T0_FIT_PAR);
249  // set parameter names
250  m_t0_fermi->SetParName(T0_PAR_NR_T0, "t_{0}");
251  m_t0_fermi->SetParName(T0_PAR_NR_T, "T");
252  m_t0_fermi->SetParName(T0_PAR_NR_BACK, "r_{b}");
253  m_t0_fermi->SetParName(T0_PAR_NR_A, "A");
254  // set fixed values
255  m_t0_fermi->FixParameter(T0_PAR_NR_BACK, rec.GetBackground());
256  m_t0_fermi->SetParameter(T0_PAR_NR_A, rec.GetHeight() - rec.GetBackground());
257  // set estimates as start values
258  m_t0_fermi->SetParameter(T0_PAR_NR_T0, rec.GetEstimatedT0());
259  // set resonable start value for T
260  m_t0_fermi->SetParameter(T0_PAR_NR_T, 3.0);
261  // perform fit - NOTE: The return value of the Fit function is not documented!
262  if (m_dir != nullptr) {
263  std::unique_ptr<TLine> ln = std::make_unique<TLine>(rec.GetFitRangeMin(), 0, rec.GetFitRangeMin(), m_time->GetMaximum());
264  ln->Write("t0_range_min");
265  ln = std::make_unique<TLine>(rec.GetFitRangeMax(), 0, rec.GetFitRangeMax(), m_time->GetMaximum());
266  ln->Write("t0_range_max");
267  m_t0_fermi->SetLineColor(3);
268  m_t0_fermi->Write();
269  }
270  m_t0_fermi->SetLineColor(2);
271  std::string fitopt("BLR");
272  if (m_settings->VerboseLevel() == 0) fitopt += "Q";
273  if (m_settings->AddFitfun()) {
274  fitopt += "+";
275  } else {
276  fitopt += "0";
277  }
278  m_time->Fit(m_t0_fermi.get(), fitopt.c_str(), "", rec.GetFitRangeMin(), rec.GetFitRangeMax());
279  if (m_settings->T0Settings()->UseTopChi2())
280  TopChi2();
281  else
282  m_chi2 = m_time->GetFunction("mt_t0_fermi")->GetChisquare() / m_time->GetFunction("mt_t0_fermi")->GetNDF();
283  m_t0_ok = true;
284  cwd->cd();
285  m_status_code = 0;
286  return true;
287  } // end T0MTHistos::NormalFit
288 
290  // T0Scramble() //
292  bool T0MTHistos::T0Scramble() {
293 #ifndef NDEBUG
294  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
295  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "Scrambling for " << m_time->GetName() << endmsg;
296 #endif
297  std::string fitopt("BLR");
298  if (m_settings->VerboseLevel() == 0) fitopt += "Q";
299  if (m_settings->AddFitfun()) {
300  fitopt += "+";
301  } else {
302  fitopt += "0";
303  }
304  // create scrambled histogram
305  char scramhistname[100];
306  snprintf(scramhistname, 100, "%s_scram", m_time->GetName());
307  std::unique_ptr<TH1F> scramhist = std::make_unique<TH1F>(scramhistname, "scrambled histogram", m_time->GetSize() - 2,
308  m_time->GetXaxis()->GetXmin(), m_time->GetXaxis()->GetXmax());
309 
310  for (int binnr = 0; binnr < m_time->GetSize(); binnr++) {
311  scramhist->SetBinContent(binnr, m_time->GetBinContent(binnr) + getTLSRandomGen()->Gaus(0, m_time->GetBinError(binnr)));
312  scramhist->SetBinError(binnr, m_time->GetBinError(binnr) * 1.41421356);
313  if (scramhist->GetBinContent(binnr) < 0) scramhist->SetBinContent(binnr, 0);
314  }
315  TDirectory *cwd = gDirectory;
316  if (m_dir) m_dir->cd();
317  MTT0PatternRecognition scramrec;
318  // perform pattern recognition
319  if (!scramrec.Initialize(scramhist.get(), m_settings)) {
320 #ifdef NDEBUG
321  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
322 #endif
323  log << MSG::WARNING << "T0MTHistos::FitT0 for tube " << m_id << ": Scrambed pattern recognition failed!" << endmsg;
324  cwd->cd();
325  return false;
326  }
327  char scrambuffer[100];
328  snprintf(scrambuffer, 100, "scrammt_t0_fermi");
329  std::unique_ptr<TF1> scramm_t0_fermi = std::make_unique<TF1>();
330  m_t0_fermi->Copy(*scramm_t0_fermi);
331  scramm_t0_fermi->SetName(scrambuffer);
332  scramm_t0_fermi->SetRange(scramrec.GetFitRangeMin(), scramrec.GetFitRangeMax());
333  // set parameter names
334  scramm_t0_fermi->SetParName(T0_PAR_NR_T0, "t_{0}");
335  scramm_t0_fermi->SetParName(T0_PAR_NR_T, "T");
336  scramm_t0_fermi->SetParName(T0_PAR_NR_BACK, "r_{b}");
337  scramm_t0_fermi->SetParName(T0_PAR_NR_A, "A");
338  // set fixed values
339  scramm_t0_fermi->FixParameter(T0_PAR_NR_BACK, scramrec.GetBackground());
340  scramm_t0_fermi->SetParameter(T0_PAR_NR_A, scramrec.GetHeight() - scramrec.GetBackground());
341  // set estimates as start values
342  scramm_t0_fermi->SetParameter(T0_PAR_NR_T0, scramrec.GetEstimatedT0());
343  // set resonable start value for T
344  scramm_t0_fermi->SetParameter(T0_PAR_NR_T, 3.0);
345  // perform fit - NOTE: The return value of the Fit function is not documented!
346  scramhist->Fit(scrambuffer, fitopt.c_str(), "", scramrec.GetFitRangeMin(), scramrec.GetFitRangeMax());
347  // set parameter for the new fit of the original histogram
348  m_time->GetListOfFunctions()->Clear();
349  // set fixed values
350  m_t0_fermi->FixParameter(T0_PAR_NR_BACK, scramm_t0_fermi->GetParameter(T0_PAR_NR_BACK));
351  m_t0_fermi->SetParameter(T0_PAR_NR_A, scramm_t0_fermi->GetParameter(T0_PAR_NR_A));
352  // set estimates as start values
353  m_t0_fermi->SetParameter(T0_PAR_NR_T0, scramm_t0_fermi->GetParameter(T0_PAR_NR_T0));
354  // set resonable start value for T
355  m_t0_fermi->SetParameter(T0_PAR_NR_T, scramm_t0_fermi->GetParameter(T0_PAR_NR_T));
356  m_time->GetListOfFunctions()->Clear();
357  m_time->Fit("mt_t0_fermi", fitopt.c_str(), "", scramrec.GetFitRangeMin(), scramrec.GetFitRangeMax());
358  if (m_settings->T0Settings()->UseTopChi2())
359  TopChi2();
360  else
361  m_chi2 = m_time->GetFunction("mt_t0_fermi")->GetChisquare() / m_time->GetFunction("mt_t0_fermi")->GetNDF();
362  cwd->cd();
363  return true;
364  } // end T0MTHistos::T0Scramble
365 
366  void T0MTHistos::TopChi2() {
367  // calculate topchi2
368  m_chi2 = 0;
369  int topndf = 0;
370  TF1 *t0_fermi = m_time->GetFunction("mt_t0_fermi");
371  Double_t min, max;
372  t0_fermi->GetRange(min, max);
373  int startbin = m_time->FindBin(min) + 1;
374  int endbin = m_time->FindBin(max) - 1;
375  for (int bin = startbin; bin < endbin; bin++) {
376  float measval = m_time->GetBinContent(bin);
377  float funcval = t0_fermi->Eval(m_time->GetBinCenter(bin));
378  float errval = m_time->GetBinError(bin);
379  // take only chi2 from top part or if the bin content is min 10 or if the function >10
380  if (measval < 10 && funcval < 10 &&
381  m_time->GetBinCenter(bin) < t0_fermi->GetParameter(T0_PAR_NR_T0) + 2 * t0_fermi->GetParameter(T0_PAR_NR_T))
382  continue;
383  if (errval == 0)
384  m_chi2 += (measval - funcval) * (measval - funcval);
385  else
386  m_chi2 += (measval - funcval) * (measval - funcval) / (errval * errval);
387  topndf++;
388  }
389  if (topndf != 0) m_chi2 = m_chi2 / topndf;
390  } // end T0MTHistos::TopChi2
391 
392  void T0MTHistos::TopSlicing() {
393 #ifndef NDEBUG
394  MsgStream log(Athena::getMessageSvc(), "T0MTHistos");
395  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "Slicing for " << m_time->GetName() << endmsg;
396 #endif
397  // vector with slice chi2
398  std::list<Slice> slice_chi2;
399  Slice current;
400  current.chi_2 = 0.0;
401  current.n_bins = 0;
402  TF1 *t0_fermi = m_time->GetFunction("mt_t0_fermi");
403  current.min_bin = m_time->FindBin(t0_fermi->GetParameter(T0_PAR_NR_T0) + 2 * t0_fermi->GetParameter(T0_PAR_NR_T));
404  Double_t min, max;
405  t0_fermi->GetRange(min, max);
406 #ifndef NDEBUG
407  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << current.min_bin << " " << max << endmsg;
408 #endif
409  for (int bin = m_time->FindBin(t0_fermi->GetParameter(T0_PAR_NR_T0) + 2 * t0_fermi->GetParameter(T0_PAR_NR_T));
410  bin < m_time->FindBin(max) - 1; bin++) {
411  if (current.n_bins == 10) {
412 #ifndef NDEBUG
413  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << current.chi_2 / current.n_bins << endmsg;
414 #endif
415  current.max_bin = bin;
416  slice_chi2.push_back(current);
417  current.chi_2 = 0.0;
418  current.min_bin = bin;
419  current.n_bins = 0;
420  }
421  double measval = m_time->GetBinContent(bin);
422  double funcval = t0_fermi->Eval(m_time->GetBinCenter(bin));
423  double errval = m_time->GetBinError(bin);
424  if (errval == 0)
425  current.chi_2 += std::pow(measval - funcval, 2.0);
426  else
427  current.chi_2 += std::pow(measval - funcval, 2.0) / std::pow(errval, 2);
428  current.n_bins++;
429  }
430 #ifndef NDEBUG
431  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "number of slices: " << slice_chi2.size() << endmsg;
432 #endif
433  std::list<Slice>::iterator it = slice_chi2.end();
434  do {
435  --it;
436  if (it == slice_chi2.begin()) {
437 #ifndef NDEBUG
438  if (log.level() <= MSG::VERBOSE) log << MSG::VERBOSE << "No gain in slicing!" << endmsg;
439 #endif
440  return;
441  }
442  } while (it->chi_2 / static_cast<double>(it->n_bins) > 3);
443  max = m_time->GetBinCenter(it->min_bin);
444  m_time->GetListOfFunctions()->Clear();
445  std::string fitopt("BLR");
446  if (m_settings->VerboseLevel() == 0) fitopt += "Q";
447  if (m_settings->AddFitfun()) {
448  fitopt += "+";
449  } else {
450  fitopt += "0";
451  }
452  m_time->Fit("mt_t0_fermi", fitopt.c_str(), "", min, max);
453  if (m_settings->T0Settings()->UseTopChi2())
454  TopChi2();
455  else
456  m_chi2 = m_time->GetFunction("mt_t0_fermi")->GetChisquare() / m_time->GetFunction("mt_t0_fermi")->GetNDF();
457  } // end T0MTHistos::TopSlicing
458 
459 } // namespace MuonCalib
xAOD::iterator
JetConstituentVector::iterator iterator
Definition: JetConstituentVector.cxx:68
fillPileUpNoiseLumi.current
current
Definition: fillPileUpNoiseLumi.py:52
MuonCalib::T0MTSettingsT0::ScrambleThreshold
const double & ScrambleThreshold() const
the chi2 threshold at which the scrambling method is used
Definition: T0MTSettingsT0.h:63
MuonCalib::T0MTHistos::FillT
void FillT(double t)
fill drift time spectrum
Definition: T0MTHistos.cxx:56
MuonCalib::MTTmaxPatternRecognition::GetFitRangeMin
double GetFitRangeMin() const
get fit range
Definition: MTTmaxPatternRecognition.h:72
MuonCalib::T0MTHistos::m_t0_fermi
std::unique_ptr< TF1 > m_t0_fermi
function fitted to the riding edghe of the spectrum
Definition: T0MTHistos.h:123
TH1F::GetBinContent
double GetBinContent(int) const
Definition: rootspy.cxx:326
max
#define max(a, b)
Definition: cfImp.cxx:41
MuonCalib::T0MTHistos::m_dir
TDirectory * m_dir
TDirectory where debug and result histograms are stored.
Definition: T0MTHistos.h:133
T0MTHistos.h
getMessageSvc.h
singleton-like access to IMessageSvc via open function and helper
MuonCalib::T0MTHistos::TopChi2
void TopChi2()
top chi2 calculation
Definition: T0MTHistos.cxx:366
MuonCalib::MTTmaxPatternRecognition::GetBackground
double GetBackground() const
get the background level
Definition: MTTmaxPatternRecognition.h:60
MuonCalib::T0MTHistos::N_T0_FIT_PAR
static constexpr int N_T0_FIT_PAR
number of parameters in t0 fit
Definition: T0MTHistos.h:55
MuonCalib::mt_t0_fermi
Double_t mt_t0_fermi(Double_t *x, Double_t *par)
The fermi function to be fitted at the rising edge of the spectrum.
Definition: T0MTHistos.cxx:34
dqt_zlumi_pandas.hname
string hname
Definition: dqt_zlumi_pandas.py:272
CaloCondBlobAlgs_fillNoiseFromASCII.spec
spec
Definition: CaloCondBlobAlgs_fillNoiseFromASCII.py:47
MuonCalib::T0MTSettingsT0::UseTopChi2
const bool & UseTopChi2() const
If true use only the top part of the function for the chi2 calculation.
Definition: T0MTSettingsT0.h:60
RootHelpers::FindBin
Int_t FindBin(const TAxis *axis, const double x)
Definition: RootHelpers.cxx:14
MuonCalib::getTLSRandomGen
TRandom3 * getTLSRandomGen()
Definition: T0MTHistos.cxx:22
conifer::pow
constexpr int pow(int x)
Definition: conifer.h:20
MuonCalib::T0MTSettings::NBinsTime
const int & NBinsTime() const
Number of bins for time histogram and range.
Definition: T0MTSettings.h:45
skel.it
it
Definition: skel.GENtoEVGEN.py:423
bin
Definition: BinsDiffFromStripMedian.h:43
MuonCalib::MTT0PatternRecognition::GetEstimatedT0
double GetEstimatedT0() const
get estimated t0
Definition: MTT0PatternRecognition.h:70
MuonCalib::T0MTHistos::m_chi2
double m_chi2
chi2/NDF value
Definition: T0MTHistos.h:137
MuonCalib::T0MTSettings::T0Settings
const T0MTSettingsT0 * T0Settings() const
get settings for the t0-fit
Definition: T0MTSettings.h:78
MuonCalib::MTT0PatternRecognition::GetHeight
double GetHeight() const
get height
Definition: MTT0PatternRecognition.h:67
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
MuonCalib::T0MTSettingsT0::SlicingThreshold
const double & SlicingThreshold() const
the chi2 threshold at which the slicing method is used
Definition: T0MTSettingsT0.h:70
drawFromPickle.exp
exp
Definition: drawFromPickle.py:36
x
#define x
Athena::getMessageSvc
IMessageSvc * getMessageSvc(bool quiet=false)
Definition: getMessageSvc.cxx:20
MuonCalib::T0MTHistos::T0_PAR_NR_T
static constexpr int T0_PAR_NR_T
Definition: T0MTHistos.h:57
MuonCalib::T0MTHistos::m_tmax_fermi
std::unique_ptr< TF1 > m_tmax_fermi
function fitted to the falling edge of the spectrum
Definition: T0MTHistos.h:129
MuonCalib::T0MTHistos::FitTmax
bool FitTmax()
Performs tmax-fit Returns true if fit is successfull.
Definition: T0MTHistos.cxx:142
MTT0PatternRecognition.h
createCoolChannelIdFile.buffer
buffer
Definition: createCoolChannelIdFile.py:12
MuonCalib::MTT0PatternRecognition::GetBackground
double GetBackground() const
get the background level
Definition: MTT0PatternRecognition.h:64
MuonCalib::T0MTHistos::T0_PAR_NR_BACK
static constexpr int T0_PAR_NR_BACK
Definition: T0MTHistos.h:57
MuonCalib::T0MTHistos::m_status_code
int m_status_code
status code for t0 fit (0 ok, 1 not fitted, 2 low statistics, 3 failed)
Definition: T0MTHistos.h:127
ParseInputs.gDirectory
gDirectory
Definition: Final2012/ParseInputs.py:133
MuonCalib::T0MTHistos::SetTSpec
void SetTSpec(int id, TH1F *spec, const T0MTSettings *settings, bool copy_spec=true)
set the pointer of the drift-time spectrum to an existing spectrum.
Definition: T0MTHistos.cxx:95
MTTmaxPatternRecognition.h
MuonCalib::T0MTSettings::VerboseLevel
const int & VerboseLevel() const
verbose level 0: no output 1: Fitter output
Definition: T0MTSettings.h:73
MuonCalib::MTT0PatternRecognition::GetFitRangeMin
double GetFitRangeMin() const
get fit range
Definition: MTT0PatternRecognition.h:73
endmsg
#define endmsg
Definition: AnalysisConfig_Ntuple.cxx:63
MuonCalib::T0MTSettings::TimeMin
const double & TimeMin() const
Definition: T0MTSettings.h:47
MuonCalib::T0MTHistos::m_t0_ok
bool m_t0_ok
is true if t0 fit was successful
Definition: T0MTHistos.h:125
MuonCalib::T0MTHistos::TMAX_PAR_NR_A
static constexpr int TMAX_PAR_NR_A
Definition: T0MTHistos.h:61
MuonCalib
CscCalcPed - algorithm that finds the Cathode Strip Chamber pedestals from an RDO.
Definition: CscCalcPed.cxx:22
MuonCalib::T0MTHistos::Initialize
void Initialize(int id, const T0MTSettings *settings, const char *hname=nullptr)
Initialize class.
Definition: T0MTHistos.cxx:61
TH1F::SetBinContent
void SetBinContent(int, double)
Definition: rootspy.cxx:327
MuonCalib::T0MTHistos::TMAX_PAR_NR_B
static constexpr int TMAX_PAR_NR_B
Definition: T0MTHistos.h:61
MuonCalib::T0MTHistos::TMAX_PAR_NR_T0
static constexpr int TMAX_PAR_NR_T0
Definition: T0MTHistos.h:62
MuonCalib::MTT0PatternRecognition::GetFitRangeMax
double GetFitRangeMax() const
get fit range
Definition: MTT0PatternRecognition.h:76
MuonCalib::T0MTHistos::T0_PAR_NR_T0
static constexpr int T0_PAR_NR_T0
parameter numbers in t0 fit
Definition: T0MTHistos.h:57
min
#define min(a, b)
Definition: cfImp.cxx:40
MuonCalib::MTTmaxPatternRecognition::GetA
double GetA() const
get parameter a in exp-function representing the end of the spectrum
Definition: MTTmaxPatternRecognition.h:63
MuonCalib::T0MTHistos::m_settings
const T0MTSettings * m_settings
Pointer to settings class.
Definition: T0MTHistos.h:135
id
SG::auxid_t id
Definition: Control/AthContainers/Root/debug.cxx:194
createCoolChannelIdFile.par
par
Definition: createCoolChannelIdFile.py:29
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:77
MuonCalib::MTT0PatternRecognition::Initialize
bool Initialize(TH1F *hist, const T0MTSettings *settings)
Initialize class - returns true if pattern recognition was successfull.
Definition: MTT0PatternRecognition.h:53
plotBeamSpotVxVal.bin
int bin
Definition: plotBeamSpotVxVal.py:83
MuonCalib::T0MTHistos::T0Scramble
bool T0Scramble()
try to get better start values from a scrambled histogram
Definition: T0MTHistos.cxx:292
MuonCalib::T0MTHistos::FitT0
bool FitT0()
Perform t0-fit Returns true if fit is successfull.
Definition: T0MTHistos.cxx:114
MuonCalib::T0MTHistos::m_tmax_ok
double m_tmax_ok
is true if tmax fit was successful
Definition: T0MTHistos.h:131
MuonCalib::T0MTHistos::m_time
std::unique_ptr< TH1F > m_time
time spectrum
Definition: T0MTHistos.h:119
cwd
std::string cwd
Definition: listroot.cxx:38
MuonCalib::T0MTSettings::AddFitfun
const bool & AddFitfun() const
If set to true the fitted functions are added to the histograms.
Definition: T0MTSettings.h:63
MuonCalib::T0MTHistos::T0_PAR_NR_A
static constexpr int T0_PAR_NR_A
Definition: T0MTHistos.h:57
MuonCalib::T0MTHistos::TMAX_PAR_NR_T
static constexpr int TMAX_PAR_NR_T
Definition: T0MTHistos.h:61
MuonCalib::T0MTHistos::m_id
int m_id
tube id;
Definition: T0MTHistos.h:121
a
TList * a
Definition: liststreamerinfos.cxx:10
MuonCalib::T0MTSettings::TimeMax
const double & TimeMax() const
Definition: T0MTSettings.h:49
TH1F
Definition: rootspy.cxx:320
MuonCalib::T0MTHistos::N_TMAX_FIT_PAR
static constexpr int N_TMAX_FIT_PAR
number of parameters for tmax fit
Definition: T0MTHistos.h:59
MuonCalib::mt_tmax_fermi
Double_t mt_tmax_fermi(Double_t *x, Double_t *par)
The fermi function to be fitted at the trailing slope of the spectrum.
Definition: T0MTHistos.cxx:44
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
MuonCalib::MTTmaxPatternRecognition::Initialize
bool Initialize(TH1F *hist, double t0, const T0MTSettings *settings)
Initialize class.
Definition: MTTmaxPatternRecognition.cxx:20
MuonCalib::MTT0PatternRecognition
Definition: MTT0PatternRecognition.h:35
MuonCalib::T0MTHistos::TMAX_PAR_NR_TMAX
static constexpr int TMAX_PAR_NR_TMAX
parameters numbers for tmax fit
Definition: T0MTHistos.h:61
MuonCalib::T0MTHistos::TopSlicing
void TopSlicing()
top slicing metyhod
Definition: T0MTHistos.cxx:392
MuonCalib::T0MTSettings::DrawDebugGraphs
const bool & DrawDebugGraphs() const
If set to true for every tube a TDirectory will be created.
Definition: T0MTSettings.h:60
MuonCalib::MTTmaxPatternRecognition::GetEstimatedTMax
double GetEstimatedTMax() const
get estimated t0
Definition: MTTmaxPatternRecognition.h:69
MuonCalib::T0MTHistos::TMAX_PAR_NR_BACK
static constexpr int TMAX_PAR_NR_BACK
Definition: T0MTHistos.h:61
ATLAS_THREAD_SAFE
#define ATLAS_THREAD_SAFE
Definition: checker_macros.h:211
python.Constants.VERBOSE
int VERBOSE
Definition: Control/AthenaCommon/python/Constants.py:14
checker_macros.h
Define macros for attributes used to control the static checker.
python.compareTCTs.ln
ln
Definition: compareTCTs.py:296
MuonCalib::MTTmaxPatternRecognition
Definition: MTTmaxPatternRecognition.h:36
MuonCalib::MTTmaxPatternRecognition::GetFitRangeMax
double GetFitRangeMax() const
get fit range
Definition: MTTmaxPatternRecognition.h:75
MuonCalib::MTTmaxPatternRecognition::GetB
double GetB() const
get parameter a in exp-function representing the end of the spectrum
Definition: MTTmaxPatternRecognition.h:66
MuonCalib::T0MTSettings
Definition: T0MTSettings.h:28
MuonCalib::T0MTHistos::Slice
Definition: T0MTHistos.h:146
MuonCalib::T0MTHistos::NormalFit
bool NormalFit()
normal t0 fit
Definition: T0MTHistos.cxx:212
Generated on Tue Jul 2 2024 21:27:35 for ATLAS Offline Software by   doxygen 1.8.18