ATLAS Offline Software
DigitMonitor.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 #include "LArSamplesMon/DigitMonitor.h"
6 
7 #include "LArSamplesMon/Interface.h"
8 #include "LArSamplesMon/History.h"
9 #include "LArSamplesMon/Data.h"
10 #include "LArSamplesMon/OFC.h"
11 #include "LArSamplesMon/UniformShapeErrorGetter.h"
12 #include "LArSamplesMon/CombinedShapeErrorGetter.h"
13 #include "LArSamplesMon/TreeShapeErrorGetter.h"
14 #include "LArSamplesMon/ScaledErrorData.h"
15 #include "LArCafJobs/SimpleShape.h"
16 #include "LArSamplesMon/Residual.h"
17 #include "LArCafJobs/Geometry.h"
18 
19 #include "TFile.h"
20 #include "TCanvas.h"
21 #include "TF1.h"
22 #include "TH1I.h"
23 #include "TH2D.h"
24 #include "TLatex.h"
25 #include "TMath.h"
26 #include "TPaveText.h"
27 #include "TROOT.h"
28 #include "TVectorD.h"
29 #include <map>
30 #include <utility>
31 
32 #include <vector>
33 
34 #include <iostream>
35 using std::cout;
36 using std::endl;
37 
38 using namespace LArSamples;
39 
40 
41 TH1D* DigitMonitor::energyDist(const TString& name, int nBins, double eMax) const
42 {return dist(&Data::_energy, DataFuncArgs(), name, nBins, 0, eMax, "Energy distribution", "Energy (MeV)"); }
43 
44 TH1D* DigitMonitor::timeDist(const TString& name, int nBins, double tMin, double tMax) const
45 { return dist(&Data::_ofcTime, DataFuncArgs(), name, nBins, tMin, tMax, "Time distribution", "Time (ns)"); }
46 
47 TH1D* DigitMonitor::adcMaxDist(const TString& name, int nBins, double aMin, double aMax) const
48 { return dist(&Data::_adcMax, DataFuncArgs(), name, nBins, aMin, aMax, "ADCMax distribution", "ADCMax (counts)"); }
49 
50 TH1D* DigitMonitor::noiseDist(const TString& name, int nBins, double max) const
51 { return dist(&Data::_noise, DataFuncArgs(), name, nBins, 0, max, "Noise distribution", "Noise (ADC counts)"); }
52 
53 TH2D* DigitMonitor::maxValueMap(TString name, PartitionId partition) const
54 { return partitionMap(&Data::_maxValue, DataFuncArgs(), std::move(name), partition, "Max sample value (ADC counts)", MaxValue); }
55 
56 TH2D* DigitMonitor::minValueMap(TString name, PartitionId partition) const
57 { return partitionMap(&Data::_minValue, DataFuncArgs(), std::move(name), partition, "Min sample value (ADC counts)", MinValue); }
58 
59 
60 TH1D* DigitMonitor::chi2Dist(const TString& name, int nBins, double max, double kFactor,
61  double fitMax, int lwb, int upb, unsigned int chi2Pars,
62  ShapeErrorType shapeErrorType, unsigned int nDof) const
63 {
64  TH1D* h = new TH1D(name, "#chi^{2} distribution", nBins, 0, max);
65  unsigned int nDofEff = 0, nDofLast = 0;
66  std::unique_ptr<UniformShapeErrorGetter> kFactorGetter;
67  std::unique_ptr<CombinedShapeErrorGetter> shapeErrorGetter;
68  if (kFactor > 0) {
69  kFactorGetter = std::make_unique<UniformShapeErrorGetter>(kFactor);
70  shapeErrorGetter = std::make_unique<CombinedShapeErrorGetter>();
71  if (interface().shapeErrorGetter()) shapeErrorGetter->add(*interface().shapeErrorGetter());
72  shapeErrorGetter->add(*kFactorGetter);
73  }
74  for (unsigned int i = 0; i < nChannels(); i++) {
75  if ((i+1) % 10000 == 0) cout << "Processing entry # " << i+1 << endl;
76  const History* history = cellHistory(i);
77  if (!history) continue;
78  if (shapeErrorGetter) history->setShapeErrorGetter(shapeErrorGetter.get()); // Do it at history level so it doesn't "stick" afterwards,,,
79  for (unsigned int j = 0; j < history->nData(); j++) {
80  h->Fill(history->chi2(j, lwb, upb, chi2Pars, shapeErrorType, &nDofEff));
81  if (nDofLast > 0 && nDofLast != nDofEff) cout << "WARNING @ hash = " << i << ", index = " << j << " : nDof varied from " << nDofLast << " to " << nDofEff << endl;
82  nDofLast = nDofEff;
83  }
84  }
85 
86  if (nDof > 0 && nDof != nDofEff) cout << "WARNING : fixing nDof = " << nDof << ", but computer nDofEff = " << nDofEff << endl;
87  h->Draw();
88  if (fitMax < 0) return h;
89  TF1* fChi2 = fitChi2(*h, "chi2", 0, max, nDof, 0, fitMax);
90  TPaveText* p = new TPaveText(0.6, 0.85, 0.85, 0.65, "NDC");
91  p->AddText(Form("k = %.1f%%", kFactor*100));
92  p->AddText(Form("Fit #chi^{2} = %.1f", fChi2->GetChisquare()));
93  p->AddText(Form("Fit nDOF = %.1f", fChi2->GetParameter(1)));
94  p->Draw();
95  return h;
96 }
97 
98 
99 TF1* DigitMonitor::chi2Func(const char* name, double xMin, double xMax)
100 {
101  TF1* fChi2 = new TF1(name, "[0]*TMath::Power(x/2, [1]/2)/x*TMath::Exp(-x/2)/TMath::Gamma([1]/2)", xMin, xMax);
102  fChi2->SetParameter(0, 100);
103  fChi2->SetParameter(1, 3);
104  fChi2->SetParLimits(1, 1, 7);
105  return fChi2;
106 }
107 
108 
109 TF1* DigitMonitor::fitChi2(TH1D& h, const char* name, double xMin, double xMax,
110  double nDof, double fitMin, double fitMax)
111 {
112  if (fitMin == Definitions::none) fitMin = xMin;
113  if (fitMax == Definitions::none) fitMin = xMax;
114 
115  TF1* fChi2 = chi2Func(name, xMin, xMax);
116  fChi2->SetParameter(0, 4*h.GetMaximum());
117  if (nDof > 0) fChi2->FixParameter(1, nDof);
118 
119  h.Fit(fChi2, "", "", fitMin, fitMax);
120  fChi2->SetRange(xMin, xMax);
121  return fChi2;
122 }
123 
124 
125 TH1D* DigitMonitor::bestChi2Dist(const TString& name, int nBins, double max,
126  double refErrMin, double refErrMax, unsigned int refErrNBins,
127  int lwb, int upb, unsigned int chi2Pars, unsigned int nDof) const
128 {
129  TF1* fChi2 = chi2Func("fChi2", 0, max);
130  if (nDof > 0) fChi2->FixParameter(1, nDof);
131 
132  double delta = (refErrMax - refErrMin)/refErrNBins;
133  TH1D* fitResults = new TH1D("fitResults", "", refErrNBins, refErrMin - delta/2, refErrMax - delta/2);
134 
135  for (unsigned int k = 0; k < refErrNBins; k++) {
136  double refErr = refErrMin + delta*k;
137  TH1D* h = new TH1D(name, "#chi^{2} distribution", nBins, 0, max);
138  UniformShapeErrorGetter kFactorGetter(refErr);
139  CombinedShapeErrorGetter shapeErrorGetter;
140  if (interface().shapeErrorGetter()) shapeErrorGetter.add(*interface().shapeErrorGetter());
141  shapeErrorGetter.add(kFactorGetter);
142  for (unsigned int i = 0; i < nChannels(); i++) {
143  const History* history = cellHistory(i);
144  if (!history) continue;
145  history->setShapeErrorGetter(&shapeErrorGetter);
146  for (unsigned int j = 0; j < history->nData(); j++) {
147  h->Fill(history->chi2(j, lwb, upb, chi2Pars));
148  }
149  }
150  h->Fit(fChi2);
151  fitResults->SetBinContent(k+1, fChi2->GetChisquare());
152  delete h;
153  }
154 
155  TF1* fPol2 = new TF1("fPol2", "[0]*(x - [1])*(x - [1]) + [2]", refErrMin, refErrMax);
156 
157  double center = fitResults->GetBinCenter(fitResults->GetMinimumBin());
158  fPol2->SetParameter(0, (fitResults->GetMaximum() - fitResults->GetMinimum())/center/center);
159  fPol2->SetParameter(1, center);
160  fPol2->SetParameter(2, fitResults->GetMinimum());
161 
162  fitResults->Draw();
163  fitResults->Fit(fPol2);
164 
165  TPaveText* p = new TPaveText(0.6, 0.85, 0.85, 0.65, "NDC");
166  p->AddText(Form("Fit k_{best} = %.1f%%", fPol2->GetParameter(1)*100));
167  p->AddText(Form("Fit #delta k_{best} = %.2f%%", 1/TMath::Sqrt(fPol2->GetParameter(0))*100));
168  p->AddText(Form("Fit #chi^{2}_{min} = %.2f", fPol2->GetParameter(2)));
169  p->Draw();
170 
171  return fitResults;
172 }
173 
174 
175 
176 TH1D* DigitMonitor::gainDist(const TString& name) const
177 {
178  TH1D* h = new TH1D(name, "gain distribution", 3, -0.5, 2.5);
179  h->GetXaxis()->SetBinLabel(1, "Low");
180  h->GetXaxis()->SetBinLabel(2, "Medium");
181  h->GetXaxis()->SetBinLabel(3, "High");
182  for (unsigned int i = 0; i < nChannels(); i++) {
183  const History* history = cellHistory(i);
184  if (!history) continue;
185  for (unsigned int j = 0; j < history->nData(); j++) h->Fill(history->data(j)->gain());
186  }
187  return h;
188 }
189 
190 
191 TH1D* DigitMonitor::layerDist(const TString& name) const
192 {
193  TH1D* h = new TH1D(name, "layer distribution", 4, -0.5, 3.5);
194  h->GetXaxis()->SetBinLabel(1, "PS");
195  h->GetXaxis()->SetBinLabel(2, "Strips");
196  h->GetXaxis()->SetBinLabel(3, "Middle");
197  h->GetXaxis()->SetBinLabel(4, "Back");
198  for (unsigned int i = 0; i < nChannels(); i++) {
199  const History* history = cellHistory(i);
200  if (!history) continue;
201  for (unsigned int j = 0; j < history->nData(); j++) h->Fill(history->cellInfo()->layer());
202  }
203  return h;
204 }
205 
206 
207 TH1D* DigitMonitor::residualDist(unsigned int k, const TString& name, int nBins, double rMin, double rMax, bool norm) const
208 {
209  TH1D* h = new TH1D(name, "Residual distribution", nBins, rMin, rMax);
210  h->GetXaxis()->SetTitle(norm ? "Residual/ADCMax (%)" : "Residual (ADC counts)");
211 
212  for (unsigned int i = 0; i < nChannels(); i++) {
213  const History* history = cellHistory(i);
214  if (!history) continue;
215  for (unsigned int j = 0; j < history->nData(); j++)
216  h->Fill(history->data(j)->delta(k)*(norm ? 100/history->data(j)->adcMax() : 1));
217  }
218  return h;
219 }
220 
221 
222 bool DigitMonitor::residualPlots(CaloId calo, unsigned int layer, CaloGain::CaloGain gain, bool xip, bool ring,
223  double lo, double hi, const TString& fileName)
224 {
225  TH1D* h1[5];
226  TH2D* hEtaPhi[5];
227  TString globalTitle = TString(xip ? "Derivative" : "Offset") + " correction";
228  for (unsigned short k = 0; k < 5; k++) {
229  TString title = globalTitle + " for sampling " + Form("%d", k);
230  h1[k] = new TH1D(Form("h%d", k), title, 100, lo, hi);
231  h1[k]->GetXaxis()->SetTitle("correction");
232  hEtaPhi[k] = Geo::etaPhiHist(calo, layer, Form("hEtaPhi%d", k), title);
233  hEtaPhi[k]->SetMinimum(lo);
234  hEtaPhi[k]->SetMaximum(hi);
235  }
236  FilterParams f;
237  f.addCalo(calo);
238  f.addLayer(layer);
239  TH2D* hAll = new TH2D("hAll", globalTitle + " vs. sampling", 5, 0, 5, 100, lo, hi);
240  hAll->GetXaxis()->SetTitle("sampling");
241  hAll->GetYaxis()->SetTitle("correction");
242  for (unsigned int i = 0; i < nChannels(); i++) {
243  if ((i+1) % 10000 == 0) cout << "Cell # " << i+1 << "/" << nChannels() << endl;
244  const History* history = pass(i, f);
245  if (!history || history->nData() == 0) continue;
246  for (unsigned short k = 0; k < 5; k++) {
247  ShapeErrorType type = (ring ? RingShapeError : CellShapeError);
248  double xi = history->data(0)->xi(k, type, gain, xip);
249  h1[k]->Fill(xi);
250  hEtaPhi[k]->Fill(history->cellInfo()->eta(), history->cellInfo()->phi(), xi);
251  if (Id::matchCalo(calo, HEC) && history->cellInfo()->region() == 1) // For HEC, take care of larger phi size in region 1
252  hEtaPhi[k]->Fill(history->cellInfo()->eta(), history->cellInfo()->phi() - Geo::phiSize(HEC_A, 1, 1), xi);
253  }
254  }
255  for (unsigned int i = 1; i <= 100; i++)
256  for (unsigned short k = 0; k < 5; k++)
257  hAll->SetBinContent(k+1, i, h1[k]->GetBinContent(i));
258  TFile* ff = TFile::Open(fileName + ".root", "RECREATE");
259  for (unsigned short k = 0; k < 5; k++) {
260  h1[k]->Write();
261  hEtaPhi[k]->Write();
262  }
263  hAll->Write();
264  ff->Close();
265  TCanvas* c1 = new TCanvas("c1", Id::str(calo) + " layer " + Form("%d", layer), 1600, 2000);
266  c1->Divide(2,3);
267  for (unsigned short k = 0; k < 5; k++) {
268  c1->cd(k+1);
269  hEtaPhi[k]->Draw("COLZ");
270  }
271  c1->cd(6);
272  hAll->Draw();
273  c1->Print(fileName + ".png");
274  c1->Print(fileName + ".eps");
275  return true;
276 }
277 
278 
279 bool DigitMonitor::residualPlotsRingComp(CaloId calo, unsigned int layer, CaloGain::CaloGain gain, bool xip,
280  double lo, double hi, const TString& fileName)
281 {
282  unsigned int nBins = 100;
283  TH1D* h[5][5];
284  TString globalTitle = TString(xip ? "Derivative" : "Offset") + " correction";
285  for (unsigned short l = 0; l < 5; l++)
286  for (unsigned short k = 0; k < 5; k++)
287  h[l][k] = new TH1D(Form("h%d%d", l, k), "", nBins, l < 2 ? -15 : lo, l < 2 ? 15 : hi);
288  FilterParams f;
289  f.addCalo(calo);
290  f.addLayer(layer);
291  TH2D* hAll[5];
292  hAll[0] = new TH2D("h0All", "#xi_{cell}/#delta#xi vs. sampling", 5, 0, 5, nBins, -15, 15);
293  hAll[1] = new TH2D("hDAll", "(#xi_{cell} - #xi_{ring})/#delta#xi vs. sampling", 5, 0, 5, nBins, -15, 15);
294  hAll[2] = new TH2D("hEAll", globalTitle + " per-cell value vs. sampling", 5, 0, 5, nBins, 0, hi);
295  hAll[3] = new TH2D("hRAll", globalTitle + " per-ring value vs. sampling", 5, 0, 5, nBins, lo, hi);
296  hAll[4] = new TH2D("hCAll", globalTitle + " per-cell error vs. sampling", 5, 0, 5, nBins, lo, hi);
297  TH2D* hCor = new TH2D("hCor", globalTitle + " per-ring vs. per-cell values", nBins, lo, hi, nBins, lo, hi);
298  hCor->GetXaxis()->SetTitle("#xi_{cell}");
299  hCor->GetYaxis()->SetTitle("#xi_{ring}");
300  for (unsigned short l = 0; l < 4; l++) {
301  hAll[l]->GetXaxis()->SetTitle("sampling");
302  hAll[l]->GetYaxis()->SetTitle("significance");
303  }
304  for (unsigned int i = 0; i < nChannels(); i++) {
305  if ((i+1) % 10000 == 0) cout << "Cell # " << i+1 << "/" << nChannels() << endl;
306  const History* history = pass(i, f);
307  if (!history || history->nData() == 0) continue;
308  for (unsigned short k = 0; k < 5; k++) {
309  const ShapeErrorData* cellSED = history->shapeErrorData(gain, CellShapeError);
310  const ShapeErrorData* ringSED = history->shapeErrorData(gain, RingShapeError);
311  if (!cellSED || !ringSED) {
312  if (cellSED) delete cellSED;
313  if (ringSED) delete ringSED;
314  continue;
315  }
316  double cellVal = (xip ? cellSED->xi()(k) : cellSED->xip()(k));
317  double ringVal = (xip ? ringSED->xi()(k) : ringSED->xip()(k));
318  double cellErr = TMath::Sqrt(xip ? cellSED->xiErr()(k,k) : cellSED->xipErr()(k,k));
319  delete cellSED;
320  delete ringSED;
321  h[0][k]->Fill(cellVal/cellErr);
322  h[1][k]->Fill((cellVal - ringVal)/cellErr);
323  h[2][k]->Fill(cellVal);
324  h[3][k]->Fill(ringVal);
325  h[4][k]->Fill(cellErr);
326  hCor->Fill(cellVal, ringVal);
327  }
328  }
329  for (unsigned int i = 1; i <= nBins; i++)
330  for (unsigned short l = 0; l < 5; l++)
331  for (unsigned short k = 0; k < 5; k++)
332  hAll[l]->SetBinContent(k+1, i, h[l][k]->GetBinContent(i));
333  TFile* ff = TFile::Open(fileName + ".root", "RECREATE");
334  for (unsigned short l = 0; l < 5; l++) {
335  for (unsigned short k = 0; k < 5; k++)
336  h[l][k]->Write();
337  hAll[l]->Write();
338  }
339  ff->Close();
340  TCanvas* c1 = new TCanvas("c1", Id::str(calo) + " layer " + Form("%d", layer), 1600, 500);
341  c1->Divide(3,1);
342  c1->cd(1); hCor->Draw();
343  TLatex l; l.SetNDC();
344  TString corrStr = Form("Correlation : %.1f%%", 100*hCor->GetCorrelationFactor());
345  l.DrawLatex(0.2, 0.8, corrStr);
346  c1->cd(2); hAll[0]->Draw("COL");
347  c1->cd(3); hAll[1]->Draw("COL");
348  c1->Print(fileName + ".png");
349  c1->Print(fileName + ".eps");
350  return true;
351 }
352 
353 
354 bool DigitMonitor::residualPlotsGainComp(CaloId calo, unsigned int layer, bool ring, bool xip,
355  double lo, double hi, const TString& fileName)
356 {
357  unsigned int nBins = 100;
358  TH1D* h[5][5];
359  TString globalTitle = TString(xip ? "Derivative" : "Offset") + " correction";
360  for (unsigned short l = 0; l < 5; l++)
361  for (unsigned short k = 0; k < 5; k++)
362  h[l][k] = new TH1D(Form("h%d%d", l, k), "", nBins, l < 2 ? -15 : lo, l < 2 ? 15 : hi);
363  FilterParams f;
364  f.addCalo(calo);
365  f.addLayer(layer);
366  TH2D* hAll[5];
367  hAll[0] = new TH2D("h0All", "#xi_{high}/#delta#xi vs. sampling", 5, 0, 5, nBins, -15, 15);
368  hAll[1] = new TH2D("hDAll", "(#xi_{high} - #xi_{low})/#delta#xi vs. sampling", 5, 0, 5, nBins, -15, 15);
369  hAll[2] = new TH2D("hEAll", globalTitle + " higher-gain value vs. sampling", 5, 0, 5, nBins, 0, hi);
370  hAll[3] = new TH2D("hRAll", globalTitle + " lower-gain value vs. sampling", 5, 0, 5, nBins, lo, hi);
371  hAll[4] = new TH2D("hCAll", globalTitle + " higher-gain error vs. sampling", 5, 0, 5, nBins, lo, hi);
372  TH2D* hCor = new TH2D("hCor", globalTitle + " lower-gain vs. higher-gain values", nBins, lo, hi, nBins, lo, hi);
373  hCor->GetXaxis()->SetTitle("#xi_{high}");
374  hCor->GetYaxis()->SetTitle("#xi_{low}");
375  for (unsigned short l = 0; l < 4; l++) {
376  hAll[l]->GetXaxis()->SetTitle("sampling");
377  hAll[l]->GetYaxis()->SetTitle("significance");
378  }
379 
380  CaloGain::CaloGain hiGain = (Id::matchCalo(calo, HEC) || Id::matchCalo(calo, FCAL) ? CaloGain::LARMEDIUMGAIN : CaloGain::LARHIGHGAIN);
381  CaloGain::CaloGain loGain = (Id::matchCalo(calo, HEC) || Id::matchCalo(calo, FCAL) ? CaloGain::LARLOWGAIN : CaloGain::LARMEDIUMGAIN);
382 
383  for (unsigned int i = 0; i < nChannels(); i++) {
384  if ((i+1) % 10000 == 0) cout << "Cell # " << i+1 << "/" << nChannels() << endl;
385  const History* history = pass(i, f);
386  if (!history || history->nData() == 0) continue;
387  for (unsigned short k = 0; k < 5; k++) {
388  const ShapeErrorData* hgSED = history->shapeErrorData(hiGain, ring ? RingShapeError : CellShapeError);
389  const ShapeErrorData* lgSED = history->shapeErrorData(loGain, ring ? RingShapeError : CellShapeError);
390  if (!hgSED || !lgSED) {
391  if (hgSED) delete hgSED;
392  if (lgSED) delete lgSED;
393  continue;
394  }
395  double hgVal = (xip ? hgSED->xi()(k) : hgSED->xip()(k));
396  double lgVal = (xip ? lgSED->xi()(k) : lgSED->xip()(k));
397  double hgErr = TMath::Sqrt(xip ? hgSED->xiErr()(k,k) : hgSED->xipErr()(k,k));
398  double error = TMath::Sqrt(xip ? hgSED->xiErr()(k,k) + lgSED->xiErr()(k,k) : hgSED->xipErr()(k,k) + lgSED->xipErr()(k,k));
399  delete hgSED;
400  delete lgSED;
401  h[0][k]->Fill(hgVal/hgErr);
402  h[1][k]->Fill((hgVal - lgVal)/error);
403  h[2][k]->Fill(hgVal);
404  h[3][k]->Fill(lgVal);
405  h[4][k]->Fill(hgErr);
406  hCor->Fill(hgVal, lgVal);
407  }
408  }
409  for (unsigned int i = 1; i <= nBins; i++)
410  for (unsigned short l = 0; l < 5; l++)
411  for (unsigned short k = 0; k < 5; k++)
412  hAll[l]->SetBinContent(k+1, i, h[l][k]->GetBinContent(i));
413  TFile* ff = TFile::Open(fileName + ".root", "RECREATE");
414  for (unsigned short l = 0; l < 5; l++) {
415  for (unsigned short k = 0; k < 5; k++)
416  h[l][k]->Write();
417  hAll[l]->Write();
418  }
419  ff->Close();
420  TCanvas* c1 = new TCanvas("c1", Id::str(calo) + " layer " + Form("%d", layer), 1600, 500);
421  c1->Divide(3,1);
422  c1->cd(1); hCor->Draw();
423  TLatex l; l.SetNDC();
424  TString corrStr = Form("Correlation : %.1f%%", 100*hCor->GetCorrelationFactor());
425  l.DrawLatex(0.2, 0.8, corrStr);
426  c1->cd(2); hAll[0]->Draw("COL");
427  c1->cd(3); hAll[1]->Draw("COL");
428  c1->Print(fileName + ".png");
429  c1->Print(fileName + ".eps");
430  return true;
431 }
432 
433 
434 TH1D* DigitMonitor::shapeErrorDist(unsigned int k, const TString& name, int nBins, double rMin, double rMax, double mean) const
435 {
436  TH1D* h = new TH1D(name, "Shape error", nBins, rMin, rMax);
437  h->GetXaxis()->SetTitle("#epsilon^{2}");
438 
439  for (unsigned int i = 0; i < nChannels(); i++) {
440  const History* history = cellHistory(i);
441  if (!history) continue;
442  for (unsigned int j = 0; j < history->nData(); j++) {
443  OFC* ofc = history->ofc(j);
444  double delta = history->data(j)->delta(k) - mean;
445  double adcMax = history->data(j)->adcMax();
446  h->Fill((delta*delta - ofc->Gamma()(k, k))/(adcMax*adcMax));
447  }
448  }
449  cout << sqrt(h->GetMean())*100 << " +/- " << sqrt(h->GetMeanError())*100 << endl;
450  return h;
451 }
452 
453 
454 double DigitMonitor::residualCorr(unsigned int k1, unsigned int k2) const
455 {
456  double sum = 0;
457  int n = 0;
458 
459  for (unsigned int i = 0; i < nChannels(); i++) {
460  const History* history = cellHistory(i);
461  if (!history) continue;
462  for (unsigned int j = 0; j < history->nData(); j++) {
463  double residual1 = history->data(j)->delta(k1);
464  double residual2 = history->data(j)->delta(k2);
465  sum += residual1*residual2;
466  n++;
467  }
468  }
469 
470  return n==0 ? 0 : sum/n;
471 }
472 
473 
474 bool DigitMonitor::residualParams(int lwb, int upb, CovMatrix& k, TVectorD& means) const
475 {
476  CovMatrix sum2N(lwb, upb), sumCN(lwb, upb);
477  TVectorD sum(lwb, upb), sumN(lwb, upb);
478  double sumA = 0, n = 0;
479  for (unsigned int i = 0; i < nChannels(); i++) {
480  const History* history = cellHistory(i);
481  if (!history) continue;
482  for (unsigned int j = 0; j < history->nData(); j++) {
483  if (history->data(j)->adcMax() < 0) continue;
484  OFC* ofc = history->ofc(j);
485  if (!ofc || ofc->hasSameRange(lwb, upb)) { cout << "Invalid index bounds!" << endl; return false; }
486  for (int i1 = lwb; i1 <= upb; i1++) {
487  sum(i1) = sum(i1) + history->data(j)->delta(i1);
488  sumN(i1) = sumN(i1) + history->data(j)->delta(i1)/TMath::Power(history->data(j)->adcMax(), 2);
489  for (int i2 = lwb; i2 <= upb; i2++) {
490  sumCN(i1, i2) = sumCN(i1, i2) + ofc->Gamma()(i1, i2)/TMath::Power(history->data(j)->adcMax(), 2);
491  sum2N(i1, i2) = sum2N(i1, i2) + history->data(j)->delta(i1)*history->data(j)->delta(i2)/TMath::Power(history->data(j)->adcMax(),2);
492  }
493  }
494  sumA += 1/TMath::Power(history->data(j)->adcMax(), 2);
495  n++;
496  }
497  }
498 
499  sumA /= n;
500  k.ResizeTo(lwb, upb);
501  means.ResizeTo(lwb, upb);
502 
503  for (int i1 = lwb; i1 <= upb; i1++) {
504  sum(i1) = sum(i1)/n;
505  sumN(i1) = sumN(i1)/n;
506  means(i1) = sum(i1);
507  }
508 
509  for (int i1 = lwb; i1 <= upb; i1++) {
510  for (int i2 = lwb; i2 <= upb; i2++) {
511  if (n != 0) {
512  sumCN(i1, i2) = sumCN(i1, i2)/n;
513  sum2N(i1, i2) = sum2N(i1, i2)/n;
514  }
515  k(i1, i2) = sum2N(i1, i2) - sumN(i1)*sum(i2) - sum(i1)*sumN(i2)
516  + sum(i1)*sum(i2)*sumA - sumCN(i1,i2);
517  k(i1, i2) = (k(i1, i2) >= 0 ? sqrt(k(i1, i2)) : -sqrt(-k(i1, i2)));
518  }
519  }
520  return true;
521 }
522 
523 
524 int DigitMonitor::combine(SimpleShape*& pshape, SimpleShape*& pref, const TString& selection, bool timeAligned) const
525 {
526  FilterParams f;
527  if (!f.set(selection)) return 0;
528  int n = 0;
529  pshape = pref = nullptr;
530  //don't change existing interface, use these internally
531  std::unique_ptr<SimpleShape> shape;
532  std::unique_ptr<SimpleShape> ref;
533  double maxSum = 0;
534 
535  for (unsigned int i = 0; i < nChannels(); i++) {
536  if ((i+1) % 50000 == 0) cout << "Cell # " << i+1 << "/" << nChannels() << endl;
537  const History* history = pass(i, f);
538  if (!history) continue;
539  for (unsigned int j = 0; j < history->nData(); j++) {
540  if (!f.passEvent(*history->data(j))) continue;
541  maxSum += history->data(j)->maxValue();
542  cout << "Adding pulse (" << n+1 << ") at hash " << i << ", index " << j << ", max = " << maxSum/(n+1) << endl;
543  //cout << i << " " << j << endl;
544  auto thisData = std::make_unique<SimpleShape>(*history->data(j));
545  auto thisRef = std::unique_ptr<SimpleShape>(history->referenceShape(j));
546  if (timeAligned) {
547  if (!SimpleShape::scaleAndShift(thisData, 1, -history->data(j)->ofcTime())) return -1;
548  if (!SimpleShape::scaleAndShift(thisRef, 1, -history->data(j)->ofcTime())) return -1;
549  }
550  if (!SimpleShape::add(shape, *thisData)) return -1;
551  if (!SimpleShape::add(ref, *thisRef)) return -1;
552  n++;
553  }
554  }
555 
556  if (n==0) return -1;
557  if (!SimpleShape::scaleAndShift(shape, 1.0/n)) return -1;
558  if (!SimpleShape::scaleAndShift(ref, 1.0/n)) return -1;
559  //hand these back in existing interface
560  pshape = shape.release();
561  pref = ref.release();
562  return n;
563 }
564 
565 Residuals* DigitMonitor::getResiduals(unsigned int hash, CaloGain::CaloGain gain, double absResTrunc, bool adjust, bool zeroTime) const
566 {
567  const History* history = cellHistory(hash);
568  if (!history) return nullptr;
569  if (adjust) history = history->adjust();
570  Residuals* residuals = history->residuals(gain, absResTrunc, false, zeroTime);
571  if (adjust) delete history;
572  return residuals;
573 }
574 
575 
576 bool DigitMonitor::makeResidualCorrections(const TString& outputFile, short resTrunc, short timeTrunc, double absResTrunc,
577  unsigned int minSize, bool weigh, bool adjust, bool zeroTime) const
578 {
579  TreeShapeErrorGetter* shapeError = new TreeShapeErrorGetter(outputFile, true);
580  std::vector< std::vector<ResidualCalculator> > ringCalcs;
581  std::vector<ResidualCalculator> gainRingCalcs(Geo::nPhiRings());
582  for (unsigned int g = 0; g < CaloGain::LARNGAIN; g++) ringCalcs.push_back(gainRingCalcs);
583 
584  cout << "Processing cells" << endl;
585 
586  for (unsigned int i = 0; i < nChannels(); i++) {
587  if ((i+1) % 10000 == 0) cout << "Processing hash = " << (i+1) << endl;
588  for (unsigned int g = 0; g < CaloGain::LARNGAIN; g++) {
589  Residuals* residuals = getResiduals(i, (CaloGain::CaloGain)g, absResTrunc, adjust, zeroTime);
590  if (residuals && (resTrunc > 0 || timeTrunc > 0)) {
591  Residuals* truncated = residuals->truncate(resTrunc, timeTrunc);
592  delete residuals;
593  residuals = truncated;
594  }
595 
596  if (residuals && residuals->size() < minSize) {
597  delete residuals;
598  residuals = nullptr;
599  }
600 
601  ResidualCalculator* resCalc = nullptr;
602  if (residuals) {
603  resCalc = residuals->calculator(weigh);
604  delete residuals;
605  }
606 
607  if (!resCalc) {
608  shapeError->addCell(ResidualCalculator(), (CaloGain::CaloGain)g);
609  continue;
610  }
611 
612  cout << i << " : Phi ring " << cellHistory(i)->cellInfo()->globalPhiRing() << " : adding " << resCalc->size() << endl;
613  shapeError->addCell(*resCalc, (CaloGain::CaloGain)g);
614  short iring = cellHistory(i)->cellInfo()->globalPhiRing();
615  if (iring >= 0) {
616  ringCalcs[g][iring].append(*resCalc);
617  if (iring == 1142) cout << ringCalcs[g][iring].regresser()->mean(0) << endl;
618  }
619  delete resCalc;
620  }
621  }
622 
623  cout << "Processing phi-symmetric data" << endl;
624 
625  for (unsigned int g = 0; g < CaloGain::LARNGAIN; g++) {
626  for (short i = 0; i < Geo::nPhiRings(); i++) {
627  if (ringCalcs[g][i].size() > 0)
628  cout << "Saving data for gain " << g << ", phi ring " << i << ", n = " << ringCalcs[g][i].size() << ", range = " << ringCalcs[g][i].rangeStr() << endl;
629  shapeError->addRing(ringCalcs[g][i], (CaloGain::CaloGain)g);
630  }
631  }
632 
633  cout << "Done!" << endl;
634 
635  delete shapeError;
636  return true;
637 }
TreeShapeErrorGetter.h
LArSamples::Geo::phiSize
static double phiSize(CaloId calo, short layer, short region, short iPhi=-1)
Definition: Geometry.cxx:344
LArSamples::MonitorBase::dist
TH1D * dist(const DataFuncSet &func, const DataFuncArgs &args, const TString &name, int nBins, double xMin, double xMax, const TString &title="", const TString &xTitle="", const TString &yTitle="", const FilterParams &f=FilterParams()) const
Definition: LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx:69
LArSamples::Geo::etaPhiHist
static TH2D * etaPhiHist(CaloId calo, short layer, const TString &name, const TString &title)
Definition: Geometry.cxx:549
PlotCalibFromCool.norm
norm
Definition: PlotCalibFromCool.py:100
LArSamples::HEC
@ HEC
Definition: CaloId.h:26
WriteBchToCool.adjust
adjust
Definition: WriteBchToCool.py:73
LArSamples::CellInfo::globalPhiRing
short globalPhiRing() const
Definition: CellInfo.cxx:125
mean
void mean(std::vector< double > &bins, std::vector< double > &values, const std::vector< std::string > &files, const std::string &histname, const std::string &tplotname, const std::string &label="")
Definition: dependence.cxx:254
LArSamples::TreeShapeErrorGetter::addCell
int addCell(const ResidualCalculator &calc, CaloGain::CaloGain gain)
Definition: TreeShapeErrorGetter.cxx:131
LArSamples::Id::matchCalo
static bool matchCalo(CaloId id, CaloId idSpec)
Definition: CaloId.cxx:188
LArSamples::FilterParams
Definition: FilterParams.h:50
python.ZdcOnlineRecMonitorConfig.partition
partition
Definition: ZdcOnlineRecMonitorConfig.py:328
ScaledErrorData.h
SimpleShape.h
LArSamples::Data::delta
double delta(short sample) const
Definition: Data.cxx:427
LArSamples::History::cellInfo
const CellInfo * cellInfo() const
Definition: History.h:56
LArSamples::CovMatrix
TMatrixTSym< double > CovMatrix
Definition: LArCalorimeter/LArCafJobs/LArCafJobs/Definitions.h:11
LArSamples::History::referenceShape
SimpleShape * referenceShape(unsigned int k, double adcMax=-1, double time=Definitions::none, bool samplesOnly=false) const
Definition: History.cxx:545
LArSamples::DigitMonitor::residualPlots
bool residualPlots(CaloId calo, unsigned int layer, CaloGain::CaloGain gain=CaloGain::LARHIGHGAIN, bool xip=false, bool ring=false, double lo=-0.05, double hi=0.05, const TString &fileName="residuals")
Definition: DigitMonitor.cxx:222
LArSamples::AbsLArCells::cellHistory
virtual const History * cellHistory(unsigned int i) const
Definition: AbsLArCells.cxx:59
LArSamples::DigitMonitor::bestChi2Dist
TH1D * bestChi2Dist(const TString &name, int nBins, double max, double refErrMin, double refErrMax, unsigned int refErrNBins, int lwb=-1, int upb=-1, unsigned int chi2Pars=DefaultChi2, unsigned int nDof=0) const
Definition: DigitMonitor.cxx:125
LArSamples::Data::adcMax
double adcMax() const
Definition: Data.h:124
max
constexpr double max()
Definition: ap_fixedTest.cxx:33
Residual.h
LArSamples::History
Definition: History.h:35
LArSamples::MinValue
@ MinValue
Definition: LArCalorimeter/LArSamplesMon/LArSamplesMon/MonitorBase.h:27
extractSporadic.c1
c1
Definition: extractSporadic.py:133
LArSamples::DigitMonitor::minValueMap
TH2D * minValueMap(TString name, PartitionId partition) const
Definition: DigitMonitor.cxx:56
LArSamples::MonitorBase::interface
const Interface & interface() const
Definition: LArCalorimeter/LArSamplesMon/LArSamplesMon/MonitorBase.h:74
CaloCondBlobAlgs_fillNoiseFromASCII.gain
gain
Definition: CaloCondBlobAlgs_fillNoiseFromASCII.py:109
UniformShapeErrorGetter.h
Geometry.h
LArSamples::SimpleShape
Definition: SimpleShape.h:25
LArSamples::Id::str
static TString str(CaloId id)
Definition: CaloId.cxx:15
LArSamples::History::residuals
Residuals * residuals(CaloGain::CaloGain gain=CaloGain::LARNGAIN, double absResTrunc=-1, bool correct=true, bool zeroTime=false) const
Definition: History.cxx:613
OFC.h
LArSamples::DigitMonitor::combine
int combine(SimpleShape *&shape, SimpleShape *&ref, const TString &selection="", bool timeAligned=true) const
Definition: DigitMonitor.cxx:524
LArSamples::SimpleShape::scaleAndShift
static bool scaleAndShift(std::unique_ptr< SimpleShape > &s1, double scale, double shift=0)
Definition: SimpleShape.cxx:137
UploadAMITag.l
list l
Definition: UploadAMITag.larcaf.py:157
DigitMonitor.h
LArSamples
Definition: AbsShape.h:24
LArSamples::DigitMonitor::residualParams
bool residualParams(int lwb, int upb, CovMatrix &k, TVectorD &means) const
Definition: DigitMonitor.cxx:474
LArSamples::Data::ofcTime
double ofcTime() const
Definition: Data.h:111
LArSamples::CellInfo::phi
double phi() const
Definition: CellInfo.h:94
LArSamples::HEC_A
@ HEC_A
Definition: CaloId.h:24
LArSamples::ResidualCalculator::size
unsigned int size() const
Definition: LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h:117
LArSamples::ResidualCalculator
Definition: LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h:101
LArSamples::History::ofc
OFC * ofc(unsigned int i, int lwb=-1, int upb=-1, double time=Definitions::none, bool useCorrs=true) const
Definition: History.cxx:207
read_hist_ntuple.h1
h1
Definition: read_hist_ntuple.py:21
CaloGain::LARNGAIN
@ LARNGAIN
Definition: CaloGain.h:19
LArSamples::TreeShapeErrorGetter::addRing
int addRing(const ResidualCalculator &calc, CaloGain::CaloGain gain)
Definition: TreeShapeErrorGetter.cxx:141
LArSamples::DigitMonitor::timeDist
TH1D * timeDist(const TString &name, int nBins, double tMin=-25, double tMax=25) const
Definition: DigitMonitor.cxx:44
LArSamples::DigitMonitor::adcMaxDist
TH1D * adcMaxDist(const TString &name, int nBins, double aMin=0, double aMax=4096) const
Definition: DigitMonitor.cxx:47
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
LArSamples::DigitMonitor::energyDist
TH1D * energyDist(const TString &name, int nBins, double eMax=10000) const
Definition: DigitMonitor.cxx:41
LArSamples::DigitMonitor::residualPlotsRingComp
bool residualPlotsRingComp(CaloId calo, unsigned int layer, CaloGain::CaloGain gain=CaloGain::LARHIGHGAIN, bool xip=false, double lo=-0.05, double hi=0.05, const TString &fileName="residuals")
Definition: DigitMonitor.cxx:279
LArSamples::CombinedShapeErrorGetter::add
void add(const AbsShapeErrorGetter &getter)
Definition: CombinedShapeErrorGetter.h:25
LArSamples::DigitMonitor::getResiduals
Residuals * getResiduals(unsigned int hash, CaloGain::CaloGain gain, double absResTrunc=-1, bool adjust=false, bool zeroTime=false) const
Definition: DigitMonitor.cxx:565
LArSamples::MonitorBase::nChannels
unsigned int nChannels() const
Definition: LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx:50
LArSamples::SimpleShape::add
bool add(unsigned int k, double value, double error)
Definition: SimpleShape.cxx:90
LArSamples::DigitMonitor::noiseDist
TH1D * noiseDist(const TString &name, int nBins, double max) const
Definition: DigitMonitor.cxx:50
LArSamples::ShapeErrorData::xip
const TVectorD & xip() const
Definition: ShapeErrorData.h:41
compareGeometries.outputFile
string outputFile
Definition: compareGeometries.py:25
python.setupRTTAlg.size
int size
Definition: setupRTTAlg.py:39
LArSamples::History::data
const Data * data(unsigned int i) const
Definition: History.cxx:91
LArSamples::DigitMonitor::shapeErrorDist
TH1D * shapeErrorDist(unsigned int k, const TString &name, int nBins, double rMin, double rMax, double mean=0) const
Definition: DigitMonitor.cxx:434
LArSamples::DigitMonitor::makeResidualCorrections
bool makeResidualCorrections(const TString &outputFile, short resTrunc=-1, short timeTrunc=-1, double absResTrunc=-1, unsigned int minSize=0, bool weigh=false, bool adjust=false, bool zeroTime=false) const
Definition: DigitMonitor.cxx:576
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
LArSamples::DigitMonitor::gainDist
TH1D * gainDist(const TString &name) const
Definition: DigitMonitor.cxx:176
LArSamples::MaxValue
@ MaxValue
Definition: LArCalorimeter/LArSamplesMon/LArSamplesMon/MonitorBase.h:27
convertTimingResiduals.sum
sum
Definition: convertTimingResiduals.py:55
lumiFormat.i
int i
Definition: lumiFormat.py:85
D3PDSizeSummary.ff
ff
Definition: D3PDSizeSummary.py:305
h
python.CaloCondTools.g
g
Definition: CaloCondTools.py:15
beamspotman.n
n
Definition: beamspotman.py:727
extractSporadic.h
list h
Definition: extractSporadic.py:96
LArSamples::Definitions::none
static const double none
Definition: LArCalorimeter/LArCafJobs/LArCafJobs/Definitions.h:18
CombinedShapeErrorGetter.h
TRT::Hit::layer
@ layer
Definition: HitInfo.h:79
LArSamples::History::setShapeErrorGetter
void setShapeErrorGetter(const AbsShapeErrorGetter *err) const
Definition: History.h:87
LArSamples::DigitMonitor::residualPlotsGainComp
bool residualPlotsGainComp(CaloId calo, unsigned int layer, bool ring=false, bool xip=false, double lo=-0.05, double hi=0.05, const TString &fileName="residuals")
Definition: DigitMonitor.cxx:354
covarianceTool.title
title
Definition: covarianceTool.py:542
LArSamples::ShapeErrorData::xipErr
const CovMatrix & xipErr() const
Definition: ShapeErrorData.h:43
LArSamples::DigitMonitor::residualDist
TH1D * residualDist(unsigned int k, const TString &name, int nBins, double rMin, double rMax, bool norm=false) const
Definition: DigitMonitor.cxx:207
GetAllXsec.kFactor
kFactor
Definition: GetAllXsec.py:95
LArSamples::ShapeErrorType
ShapeErrorType
Definition: LArCalorimeter/LArCafJobs/LArCafJobs/Definitions.h:22
LArSamples::DigitMonitor::layerDist
TH1D * layerDist(const TString &name) const
Definition: DigitMonitor.cxx:191
hist_file_dump.f
f
Definition: hist_file_dump.py:140
LArSamples::PartitionId
PartitionId
Definition: CaloId.h:29
LArSamples::ShapeErrorData::xiErr
const CovMatrix & xiErr() const
Definition: ShapeErrorData.h:42
LArSamples::History::nData
unsigned int nData() const
Definition: History.h:51
LArSamples::Residuals::calculator
ResidualCalculator * calculator(bool weigh=false) const
Definition: LArCalorimeter/LArSamplesMon/src/Residual.cxx:185
LArSamples::History::chi2
double chi2(int i, int lwb=-1, int upb=-1, int chi2Params=DefaultChi2, ShapeErrorType shapeErrorType=BestShapeError, unsigned int *nDof=0) const
Definition: History.cxx:164
LArSamples::OFC
Definition: OFC.h:27
LArSamples::Residuals
Definition: LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h:72
LArSamples::Data::_ofcTime
double _ofcTime(const DataFuncArgs &) const
Definition: Data.h:192
dumpTgcDigiJitter.nBins
list nBins
Definition: dumpTgcDigiJitter.py:29
LArSamples::AbsLArCells::pass
const History * pass(unsigned int i, const FilterParams &f) const
Definition: AbsLArCells.cxx:99
selection
const std::string selection
Definition: fbtTestBasics.cxx:75
LArSamples::DigitMonitor::maxValueMap
TH2D * maxValueMap(TString name, PartitionId partition) const
Definition: DigitMonitor.cxx:53
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
LArSamples::AbsShape::maxValue
double maxValue(bool withErrors=false) const
Definition: AbsShape.cxx:30
LArSamples::Data::_minValue
double _minValue(const DataFuncArgs &) const
Definition: Data.h:200
LArSamples::Residuals::truncate
Residuals * truncate(double nWidthsRes, double nWidthsTime=-1, unsigned int nMax=0) const
Definition: LArCalorimeter/LArSamplesMon/src/Residual.cxx:143
LArSamples::Data::gain
CaloGain::CaloGain gain() const
Definition: Data.h:85
LArSamples::Data::_maxValue
double _maxValue(const DataFuncArgs &) const
Definition: Data.h:199
LArSamples::Data::_noise
double _noise(const DataFuncArgs &) const
Definition: Data.h:198
LArSamples::DigitMonitor::fitChi2
static TF1 * fitChi2(TH1D &h, const char *name, double xMin, double xMax, double nDof=-1, double fitMin=Definitions::none, double fitMax=Definitions::none)
Definition: DigitMonitor.cxx:109
LArSamples::MonitorBase::partitionMap
TH2D * partitionMap(const DataFuncSet &func, const DataFuncArgs &args, TString name, PartitionId partition, const TString &title="", CombinationType comb=AverageValue, const FilterParams &f=FilterParams()) const
Definition: LArCalorimeter/LArSamplesMon/src/MonitorBase.cxx:113
LArSamples::TreeShapeErrorGetter
Definition: TreeShapeErrorGetter.h:30
CaloGain::LARHIGHGAIN
@ LARHIGHGAIN
Definition: CaloGain.h:18
LArSamples::DigitMonitor::means
TVectorD means(int lwb, int upb)
Definition: DigitMonitor.h:79
LArSamples::CaloId
CaloId
Definition: CaloId.h:21
LArSamples::DigitMonitor::chi2Dist
TH1D * chi2Dist(const TString &name, int nBins, double max, double kFactor=0, double fitMax=-1, int lwb=-1, int upb=-1, unsigned int chi2Pars=DefaultChi2, ShapeErrorType shapeErrorType=BestShapeError, unsigned int nDof=0) const
Definition: DigitMonitor.cxx:60
LArSamples::Data::_energy
double _energy(const DataFuncArgs &) const
Definition: Data.h:191
CaloGain::CaloGain
CaloGain
Definition: CaloGain.h:11
LArSamples::FCAL
@ FCAL
Definition: CaloId.h:26
CaloGain::LARMEDIUMGAIN
@ LARMEDIUMGAIN
Definition: CaloGain.h:18
LArSamples::History::shapeErrorData
const ShapeErrorData * shapeErrorData(CaloGain::CaloGain gain, ShapeErrorType shapeErrorType=BestShapeError, const Residual *res=0) const
Definition: History.cxx:302
CaloCondBlobAlgs_fillNoiseFromASCII.hash
dictionary hash
Definition: CaloCondBlobAlgs_fillNoiseFromASCII.py:108
LArSamples::DigitMonitor::residualCorr
double residualCorr(unsigned int k1, unsigned int k2) const
Definition: DigitMonitor.cxx:454
ref
const boost::regex ref(r_ef)
LArSamples::Residuals::size
unsigned int size() const
Definition: LArCalorimeter/LArSamplesMon/LArSamplesMon/Residual.h:80
LArSamples::History::adjust
History * adjust() const
Definition: History.cxx:259
LArSamples::DataFuncArgs
Definition: Data.h:34
LArSamples::CombinedShapeErrorGetter
Definition: CombinedShapeErrorGetter.h:18
LArSamples::CellInfo::region
short region() const
Definition: CellInfo.h:59
LArSamples::Data::xi
double xi(short sample, ShapeErrorType shapeErrorType=BestShapeError, CaloGain::CaloGain g=CaloGain::UNKNOWNGAIN, bool xip=false) const
Definition: Data.cxx:459
LArSamples::CellShapeError
@ CellShapeError
Definition: LArCalorimeter/LArCafJobs/LArCafJobs/Definitions.h:22
Data.h
LArSamples::DigitMonitor::chi2Func
static TF1 * chi2Func(const char *name, double xMin, double xMax)
Definition: DigitMonitor.cxx:99
LArSamples::CellInfo::eta
double eta() const
Definition: CellInfo.h:93
LArSamples::ShapeErrorData
Definition: ShapeErrorData.h:19
History.h
jobOptions.fileName
fileName
Definition: jobOptions.SuperChic_ALP2.py:39
CaloGain::LARLOWGAIN
@ LARLOWGAIN
Definition: CaloGain.h:18
LArSamples::ShapeErrorData::xi
const TVectorD & xi() const
Definition: ShapeErrorData.h:40
error
Definition: IImpactPoint3dEstimator.h:70
LArSamples::Geo::nPhiRings
static short nPhiRings()
Definition: Geometry.cxx:84
ReadOfcFromCool.ofc
ofc
Definition: ReadOfcFromCool.py:110
Interface.h
LArSamples::Data::_adcMax
double _adcMax(const DataFuncArgs &) const
Definition: Data.h:196
LArSamples::CellInfo::layer
short layer() const
Definition: CellInfo.h:53
LArSamples::RingShapeError
@ RingShapeError
Definition: LArCalorimeter/LArCafJobs/LArCafJobs/Definitions.h:24
fitman.k
k
Definition: fitman.py:528
LArSamples::UniformShapeErrorGetter
Definition: UniformShapeErrorGetter.h:17
Generated on Thu Oct 30 2025 21:08:59 for ATLAS Offline Software by   doxygen 1.8.18