ATLAS Offline Software
Functions | Variables
CalibrateIBL.cxx File Reference
#include "CxxUtils/checker_macros.h"
#include "TFile.h"
#include "TDirectory.h"
#include "TDirectoryFile.h"
#include "TKey.h"
#include "TH2D.h"
#include "TH2F.h"
#include "TProfile.h"
#include "TGraph.h"
#include "TGraphErrors.h"
#include "TMultiGraph.h"
#include "TMath.h"
#include "TF1.h"
#include "TStyle.h"
#include "TVectorT.h"
#include "THStack.h"
#include "TRandom3.h"
#include "ChargeCalibration/common/PixelMapping.h"
#include "PathResolver/PathResolver.h"
#include <cmath>
#include <fstream>
#include <iostream>
#include <iterator>
#include <string>
#include <sstream>
#include <utility>
#include <map>
#include <array>

Go to the source code of this file.

Functions

double funcDisp (double *x, double *par)
 
double ChrgfromToT (float tot, double par[10])
 
double ToTfromChrg (float chrg, double par[5])
 
double funcRation5 (double *x, double *par)
 
double funcDispConvoluted (double *x, double *par)
 
int iblCalib (const std::string &InDir, const std::string &THRscan, const std::string &scan, const std::string &scanLowCharge)
 
void printError ()
 
int main (int argc, char *argv[])
 

Variables

 ATLAS_NO_CHECK_FILE_THREAD_SAFETY
 

Function Documentation

◆ ChrgfromToT()

double ChrgfromToT ( float  tot,
double  par[10] 
)

Definition at line 69 of file CalibrateIBL.cxx.

70 {
71  int lr = (tot <= 6. ? 0 : 5);
72  double chrg = 1400.;
73 
74  double num = par[0 + lr] + tot * (par[1 + lr] + tot * (par[2 + lr] + tot * par[3 + lr]));
75  double denom = 1. + tot * par[4 + lr];
76  if (denom != 0.)
77  chrg = num / denom;
78 
79  return chrg;
80 }

◆ funcDisp()

double funcDisp ( double *  x,
double *  par 
)

Definition at line 62 of file CalibrateIBL.cxx.

63 {
64  double ret = 9.9e10;
65  ret = par[0] * 0.1 * (1. - exp(-(x[0] - 1400.) * (x[0] - 1400.) / (1000. * par[1])) + par[2] * 0.0001 * sqrt(x[0]));
66  return ret;
67 }

◆ funcDispConvoluted()

double funcDispConvoluted ( double *  x,
double *  par 
)

Definition at line 102 of file CalibrateIBL.cxx.

103 {
104  double ret = 9.9e09;
105  ret = -par[0] * (x[0] - par[1]) + par[2];
106  return ret;
107 }

◆ funcRation5()

double funcRation5 ( double *  x,
double *  par 
)

Definition at line 92 of file CalibrateIBL.cxx.

93 {
94  double ret = 9.9e10;
95  double num = par[0] + x[0] * (par[1] + x[0] * (par[2] + x[0] * par[3]));
96  double denom = 1. + x[0] * par[4];
97  if (denom != 0.0)
98  ret = num / denom;
99  return ret;
100 }

◆ iblCalib()

int iblCalib ( const std::string &  InDir,
const std::string &  THRscan,
const std::string &  scan,
const std::string &  scanLowCharge 
) 
    TF1 f1Disp( "f1Disp", funcDisp_old, chrgAbaciIBL[1], chrgAbaciIBL[nchargeIBL-1], 2 );
    grTotSig->Fit( &f1Disp,"MRQ" );
    Double_t parP0 = 1.*f1Disp.GetParameter(0);
    Double_t parP1 = 1.*f1Disp.GetParameter(1);
    Double_t parP2 = 0. ;

    TF1 f1Disp( "f1Disp", funcDisp, chrgAbaciIBL[1], chrgAbaciIBL[nchargeIBL-1], 3 );
    f1Disp.SetParLimits( 0, 3, 7 ) ;
    f1Disp.SetParLimits( 1, 300., 1000. ) ;
    f1Disp.SetParLimits( 2, -5, 5 ) ;

    grTotSig->Fit( &f1Disp,"MRQ" );
    Double_t parP0 = 1.*f1Disp.GetParameter(0);
    Double_t parP1 = 1.*f1Disp.GetParameter(1);
    Double_t parP2 = 1.*f1Disp.GetParameter(2);

paraToTRes[sfe][0] = parP0; paraToTRes[sfe][1] = parP1; paraToTRes[sfe][2] = parP2;

Definition at line 109 of file CalibrateIBL.cxx.

110 {
111  // creating the object for the pixel mapping
112  PixelMapping pixmap(PathResolver::find_file("PixelCalibAlgs/mapping.csv", "DATAPATH"));
113 
114  constexpr bool run3 = true; // new series of injection charges (22) and new format
115  const float badRDfracCut = 0.2; // Cut on the fraction of bad read out per frontend
116  float HDCshift = 2;
117 
118  // Shall we try this later ?
119  constexpr int finerToT = 1; // smaller bin size for ToT if doing fit
120  constexpr bool moreFE = false;
121  constexpr int npsFEs = moreFE ? 8 : 2;
122 
123  gStyle->SetOptFit(1100);
124 
125 #if defined(DEMOXCHECK)
126 
127  std::cout << " DEMOXCHECK will run ... " << std::endl;
128 
129  // For debugging a small area of pixels will print their readout and injection
130  int XcheckCharge[2] = {13, 14}, XcheckToT = 10;
131  int XcheckPhi[2] = {57, 63}, XcheckEta[2] = {23, 29};
132  TString XcheckModule = "LI_S11_A_M3_A6";
133 
134  // fill some histograms for monitoring
135  int ToTfill = 3 * finerToT, CHRGfill = 3;
136  const int HDC = 2; // specifying output
137 
138 #endif
139 
140  Double_t tmpRCf_H[4] = {0.868295, 0.67523, 0.732279, 0.678458}; // initial values from history calibration
141  Double_t tmpRCf_T[4] = {1.02844, 1.10788, 1.07204, 1.11015}; // initial values from history calibration
142  Double_t reverseCF_H[4], reverseCF_T[4];
143  for (int t = 0; t < 4; t++)
144  {
145  reverseCF_H[t] = tmpRCf_H[t];
146  reverseCF_T[t] = tmpRCf_T[t];
147  }
148 
149  Double_t RDentries[33] = {0, 10494.8, 14605.6, 62952.6, 8663.57, 17320.2, 20417, 49508.9,
150  21646.8, 31503.9, 14955.6, 16290.7, 11617.5, 16461.2, 11980.5, 17210.5,
151  12773.7, 18574, 13930.9, 19777, 14629.4, 20619.5, 15191.6, 20883,
152  15411.1, 19198, 12097.4, 12109.7, 8881.96, 17235.1, 20802.3, 11265.8,
153  0.};
154 
155  std::vector<int> reversedModules = {182, 211, 213, 229, 223, 231, 246, 280, 295, 350, 357, 360, 371, 403, 416, 423};
156 
157  // outputs : something@somewhere as one like ...
158  TString Outdir = "./";
159  string spec = "ToTbin" + to_string(finerToT) + "_FrtEnd" + to_string(npsFEs) + "_";
160 
161 #if defined(DEMOXCHECK)
162  spec += "DEMOXCHECK_";
163 #endif
164 
165  TString StrFileName = spec + scan;
166  TString rootFileName = Outdir + "/TotChargeCalib_" + StrFileName + ".root";
167  TString logFileName = Outdir + "/ChargeCalib_" + StrFileName + ".log";
168  TString dbFileName = Outdir + "/ChargeCalib_" + StrFileName + ".TXT";
169  if (!run3)
170  {
171  rootFileName = Outdir + "/TotChargeCalib_run2_" + StrFileName + ".root";
172  logFileName = Outdir + "/ChargeCalib_run2_" + StrFileName + ".log";
173  dbFileName = Outdir + "/ChargeCalib_run2_" + StrFileName + ".TXT";
174  }
175 
176  // fetch pictures from directory : ToT/ROD_L1_Sxx/LI_Sy_AC_Mz_ACk
177  ofstream logout(logFileName);
178  ofstream txtDB(dbFileName);
179 
180  TString inThrFile = "";
181  TString inTotFile = "";
182  TString rodPath = "";
183  TString inTotFileAux = "";
184  TString rodPathAux = "";
185 
186  std::cout << " Running IBL calibration analysis ... " << endl;
187 
188  inThrFile = InDir + THRscan + ".root";
189  inTotFile = InDir + scan + ".root";
190  rodPath = inTotFile + ":/" + scan.substr(scan.find("_")+1);
191  inTotFileAux = InDir + scanLowCharge + ".root";
192  rodPathAux = inTotFileAux + ":/" + scanLowCharge.substr(scanLowCharge.find("_")+1);
193 
194  // historical record for run2 fomat
195  if (!run3)
196  {
197  inThrFile = "calib2018/IBL/THR_SCAN_S000083686.root";
198  inTotFile = "calib2018/IBL/TOT_SCAN_S000083690.root";
199  rodPath = inTotFile + ":/S000083690";
200  }
201 
202  // selecting Q threshold
203  int nrow = 336, ncol = 160; // y-axis, x-axis respectively
204 
205  const Int_t numChrgs = 22;
206  Double_t IBLchrgs[numChrgs] = {1400., 1500., 1750., 2000., 2500., 3000., 3500., 4000., 5000., 6000., 8000., 10000.,
207  12000., 14000., 16000., 18000., 20000., 22000., 24000., 26000., 28000., 30000.};
208 
209  constexpr Int_t nchargeIBL = (run3 ? 22 : 19); //run3 is constexpr
210 
211  int skip = numChrgs - nchargeIBL;
212  Double_t chrgAbaciIBL[nchargeIBL], chargeErrArrIBL[nchargeIBL];
213  for (int c = 0; c < nchargeIBL; c++)
214  {
215  chrgAbaciIBL[c] = IBLchrgs[c + skip];
216  chargeErrArrIBL[c] = 0.;
217  }
218 
219  Double_t chrgsbins[nchargeIBL + 1];
220  for (int c = 1; c < nchargeIBL; c++)
221  chrgsbins[c] = 0.5 * (chrgAbaciIBL[c - 1] + chrgAbaciIBL[c]);
222 
223  chrgsbins[0] = chrgAbaciIBL[0] - 0.5 * (chrgAbaciIBL[1] - chrgAbaciIBL[0]);
224  chrgsbins[nchargeIBL] = chrgAbaciIBL[nchargeIBL - 1] + 0.5 * (chrgAbaciIBL[nchargeIBL - 1] - chrgAbaciIBL[nchargeIBL - 2]);
225 
226  // please note extra +1 for ending bins.
227  constexpr Int_t nToTibl = 16 * finerToT + 1;
228 
229  Double_t totAbaci[nToTibl], totbins[nToTibl + 1];
230 
231  totbins[0] = -0.25;
232  totbins[1] = 0.5;
233  totAbaci[0] = 0.0;
234  totAbaci[1] = 1.;
235  for (int t = 2; t < nToTibl; t++)
236  {
237  totAbaci[t] = totAbaci[1] + (t - 1) / (1. * finerToT);
238  totbins[t] = 0.5 * (totAbaci[t] + totAbaci[t - 1]);
239  }
240  totbins[nToTibl] = totAbaci[nToTibl - 1] + 0.5 * (totbins[nToTibl - 1] - totbins[nToTibl - 2]);
241 
242  Double_t totErrArr[nToTibl];
243  for (int t = 0; t < nToTibl; t++)
244  totErrArr[t] = 0.;
245 
246 #if defined(DEMOXCHECK)
247  for (int t = 0; t < nToTibl; t++)
248  logout << " totAbaci : " << t << " " << totbins[t] << " < " << totAbaci[t] << " > " << totbins[t + 1] << endl;
249 
250  if (finerToT == 1)
251  {
252  for (int t = 1; t < nToTibl; t++)
253  {
254  RDentries[t] = RDentries[2 * t - 1] + RDentries[2 * t];
255  logout << " finer RDentries : " << t << " " << 2 * t - 1 << " " << 2 * t << " " << RDentries[t] << endl;
256  }
257  }
258 #endif
259 
260  map<string, map<string, map<string, float>>> pcdMap;
261 
262  TFile roFile(rootFileName, "RECREATE");
263  TDirectory *roThrDir = roFile.mkdir("Threshold");
264  TDirectory *roTotDir = roFile.mkdir("ToT");
265 
266  // ****************************** Threshold IBL ******************************
267 
268  Double_t THR_avg[2][npsFEs], ThrSig_avg[2][npsFEs];
269  array<std::unique_ptr<TH2D> , npsFEs> h2_Thr{};
270  array<std::unique_ptr<TH2D> , npsFEs> h2_ThrSig{};
271 
272  std::multimap<float, TString, std::greater<float>> badThr_Order;
273 
274  for (int sfe = 0; sfe < npsFEs; sfe++)
275  {
276  stringstream ss;
277  ss << sfe;
278 
279  string knowModule = "LI_S06_C_M1_C1";
280  string idx = "I" + ss.str();
281 
282  pcdMap[knowModule][idx]["ThrNorm"] = -42.;
283  pcdMap[knowModule][idx]["ThrRmsNorm"] = -42.;
284  pcdMap[knowModule][idx]["ThrSigNorm"] = -42.;
285  pcdMap[knowModule][idx]["ThrLong"] = -42.;
286  pcdMap[knowModule][idx]["ThrRmsLong"] = -42.;
287  pcdMap[knowModule][idx]["ThrSigLong"] = -42.;
288 
289  for (int nl = 0; nl < 2; nl++)
290  THR_avg[nl][sfe] = ThrSig_avg[nl][sfe] = -99.;
291 
292  idx = "Threshold" + ss.str();
293  h2_Thr[sfe] = std::make_unique<TH2D>(idx.c_str(), " ", 2, 0, 2, 200, 0, 5000);
294  idx = "ThresholdSig" + ss.str();
295  h2_ThrSig[sfe] = std::make_unique<TH2D>(idx.c_str(), " ", 2, 0, 2, 200, 0, 500);
296  }
297 
298  if (inThrFile.Length() > 0)
299  {
300  std::cout << endl
301  << "INFO =>> [IBL] threshold scan analysis..." << endl;
302 
303  TFile riThrFile(inThrFile, "READ");
304  TString chi2HistName = "SCURVE_CHI2";
305  TString thrHistName = "SCURVE_MEAN";
306  TString sigHistName = "SCURVE_SIGMA";
307 
308  std::unique_ptr<TH1F> h1dChi2 = std::make_unique<TH1F>("h1dChi2", "", 200, 0, 1);
309  std::unique_ptr<TH1F> h1dThr = std::make_unique<TH1F>("h1dThr" , "", 200, 0, 5000);
310  std::unique_ptr<TH1F> h1dSig = std::make_unique<TH1F>("h1dSig" , "", 200, 0, 500);
311 
312  TDirectoryFile *scanDir = (TDirectoryFile *)((TKey *)riThrFile.GetListOfKeys()->First())->ReadObj();
313  TList *rodKeyList = (TList *)scanDir->GetListOfKeys();
314  TIter rodItr(rodKeyList);
315  TKey *rodKey;
316 
317  while ((rodKey = (TKey *)rodItr()))
318  {
319  TString rodName(rodKey->GetName());
320  TDirectoryFile *rodDir = (TDirectoryFile *)rodKey->ReadObj();
321  TList *modKeyList = (TList *)rodDir->GetListOfKeys();
322  TIter modItr(modKeyList);
323  TKey *modKey;
324 
325  TDirectory *dirRod = roThrDir->mkdir(rodName);
326 
327  while ((modKey = (TKey *)modItr()))
328  {
329  TString modName(modKey->GetName());
330  string modStr(modKey->GetName());
331  TDirectory *dirMod = dirRod->mkdir(modName);
332 
333  TString chi2HistDirPath = modName + "/" + chi2HistName + "/A0/B0";
334 
335  TDirectory *chi2HistDir = (TDirectory *)rodDir->Get(chi2HistDirPath);
336  if (chi2HistDir == NULL)
337  {
338  cout << " Warning : NULL dir " << endl;
339  continue;
340  }
341 
342  std::unique_ptr<TH2D> h2dChi2 (static_cast<TH2D*>(static_cast<TKey*>(chi2HistDir->GetListOfKeys()->First())->ReadObj()));
343  if (h2dChi2 == NULL)
344  {
345  cout << " Warning : NULL dir " << endl;
346  continue;
347  }
348  TString thrHistDirPath = modName + "/" + thrHistName + "/A0/B0";
349  TDirectoryFile *thrHistDir = (TDirectoryFile *)rodDir->Get(thrHistDirPath);
350  std::unique_ptr<TH2D> h2dThr(static_cast<TH2D*>(static_cast<TKey*>(thrHistDir->GetListOfKeys()->First())->ReadObj()));
351  TString sigHistDirPath = modName + "/" + sigHistName + "/A0/B0";
352 
353  TDirectoryFile *sigHistDir = (TDirectoryFile *)rodDir->Get(sigHistDirPath);
354  std::unique_ptr<TH2D> h2dSig (static_cast<TH2D*>(static_cast<TKey*>(sigHistDir->GetListOfKeys()->First())->ReadObj()));
355 
356  array<std::unique_ptr<TH1F>, npsFEs> h1_ThrNorm{};
357  array<std::unique_ptr<TH1F>, npsFEs> h1_ThrSigNorm{};
358  array<std::unique_ptr<TH1F>, npsFEs> h1_ThrLong{};
359  array<std::unique_ptr<TH1F>, npsFEs> h1_ThrSigLong{};
360  array<float, npsFEs> IlledThr{};
361 
362  for (int sfe = 0; sfe < npsFEs; sfe++)
363  {
364  stringstream ss;
365  ss << sfe;
366 
367  string hname = modStr + "ThrNorm" + ss.str();
368  h1_ThrNorm[sfe] = std::make_unique<TH1F>(hname.c_str(), "", 200, 0, 5000);
369  hname = modStr + "ThrSigNorm" + ss.str();
370  h1_ThrSigNorm[sfe] = std::make_unique<TH1F>(hname.c_str(), "", 200, 0, 500);
371  hname = modStr + "ThrLong" + ss.str();
372  h1_ThrLong[sfe] = std::make_unique<TH1F>(hname.c_str(), "", 200, 0, 5000);
373  hname = modStr + "ThrSigLong" + ss.str();
374  h1_ThrSigLong[sfe] = std::make_unique<TH1F>(hname.c_str(), "", 200, 0, 500);
375  }
376 
377  for (int col = 1; col <= ncol; col++)
378  {
379  for (int row = 1; row <= nrow; row++)
380  {
381  float chi2 = h2dChi2->GetBinContent(col, row);
382  float thr = h2dThr->GetBinContent(col, row);
383  float sig = h2dSig->GetBinContent(col, row);
384 
385  h1dChi2->Fill(chi2);
386  h1dThr->Fill(thr);
387  h1dSig->Fill(sig);
388  // Shall we measure/count the goodness of THR scan here ?
389  bool filled = true;
390  if (thr == 0 || thr > 10000 || sig == 0 || sig > 1000 || chi2 > 0.5 || chi2 <= 0)
391  {
392  filled = false;
393  }
394  int circ = -1;
395  if (filled)
396  {
397  if (col == 1 || col == ncol / 2 || col == ncol / 2 + 1 || col == ncol)
398  {
399  if (moreFE)
400  {
401  circ = (int)(row / 84.);
402  if (col > 80)
403  {
404  circ = 7 - circ;
405  }
406  }
407  else
408  {
409  if (col <= ncol / 2)
410  circ = 0;
411  else
412  circ = 1;
413  }
414 
415  h1_ThrLong[circ]->Fill(thr);
416  h1_ThrSigLong[circ]->Fill(sig);
417  }
418  else
419  {
420  if (moreFE)
421  {
422  circ = (int)(row / 84.);
423  if (col > 80)
424  {
425  circ = 7 - circ;
426  }
427  }
428  else
429  {
430  if (col <= ncol / 2)
431  circ = 0;
432  else
433  circ = 1;
434  }
435 
436  h1_ThrNorm[circ]->Fill(thr);
437  h1_ThrSigNorm[circ]->Fill(sig);
438  }
439  }
440  else
441  {
442  // Note: what is the intention? circ is -1 here
443  continue;
444  }
445  }
446  }
447  h2dSig.reset();
448  h2dThr.reset();
449  h2dChi2.reset();
450 
451  for (int sfe = 0; sfe < npsFEs; sfe++)
452  {
453  stringstream ss;
454  ss << sfe;
455  string feName = "I" + ss.str();
456 
457  float THRnorm = h1_ThrNorm[sfe]->GetMean();
458  bool valid = THRnorm > 100.;
459  pcdMap[modStr][feName]["ThrNorm"] = THRnorm;
460  pcdMap[modStr][feName]["ThrRmsNorm"] = h1_ThrNorm[sfe]->GetRMS();
461  float THRnormSig = h1_ThrSigNorm[sfe]->GetMean();
462  pcdMap[modStr][feName]["ThrSigNorm"] = THRnormSig;
463  float THRlong = h1_ThrLong[sfe]->GetMean();
464  pcdMap[modStr][feName]["ThrLong"] = THRlong;
465  pcdMap[modStr][feName]["ThrRmsLong"] = h1_ThrLong[sfe]->GetRMS();
466  float THRlongSig = h1_ThrSigLong[sfe]->GetMean();
467  pcdMap[modStr][feName]["ThrSigLong"] = THRlongSig;
468 
469  float blank = 100. * IlledThr[sfe] * npsFEs / (1. * ncol * nrow);
470 
471  if (blank > 0.01)
472  badThr_Order.insert(std::pair<float, TString>(blank, (TString)(modName + " : " + blank)));
473 
474  if (valid)
475  {
476  h2_Thr[sfe]->Fill(0.5, THRnorm, 1);
477  h2_Thr[sfe]->Fill(1.5, THRlong, 1);
478  h2_ThrSig[sfe]->Fill(0.5, THRnormSig, 1);
479  h2_ThrSig[sfe]->Fill(1.5, THRlongSig, 1);
480  }
481  dirMod->WriteTObject(h1_ThrNorm[sfe].get());
482  dirMod->WriteTObject(h1_ThrLong[sfe].get());
483 
484  h1_ThrNorm[sfe].reset();
485  h1_ThrSigNorm[sfe].reset();
486  h1_ThrLong[sfe].reset();
487  h1_ThrSigLong[sfe].reset();
488  }
489  }
490  }
491  roThrDir->WriteTObject(h1dChi2.get());
492  h1dChi2.reset();
493  roThrDir->WriteTObject(h1dThr.get());
494  h1dThr.reset();
495  roThrDir->WriteTObject(h1dSig.get());
496  h1dSig.reset();
497  }
498 
499  for (int sfe = 0; sfe < npsFEs; sfe++)
500  {
501  stringstream ss;
502  ss << sfe;
503  string hName = "thrNorm" + ss.str();
504  THR_avg[0][sfe] = h2_Thr[sfe]->ProjectionY(hName.c_str(), 1, 1)->GetMean();
505  hName = "thrLong" + ss.str();
506  THR_avg[1][sfe] = h2_Thr[sfe]->ProjectionY(hName.c_str(), 2, 2)->GetMean();
507  hName = "thrSigNorm" + ss.str();
508  ThrSig_avg[0][sfe] = h2_ThrSig[sfe]->ProjectionY(hName.c_str(), 1, 1)->GetMean();
509  hName = "thrSigLong" + ss.str();
510  ThrSig_avg[1][sfe] = h2_ThrSig[sfe]->ProjectionY(hName.c_str(), 2, 2)->GetMean();
511  h2_Thr[sfe].reset();
512  h2_ThrSig[sfe].reset();
513  }
514  // ****************************** End - Threshold IBL ******************************
515 
517  /************************* IBL ToT ***************************/
519 
520  if (inTotFile.Length() == 0 || inTotFileAux.Length() == 0)
521  {
522  logout << " Missing ToT file from calib scan as input " << endl;
523  return 0;
524  }
525 
526  std::cout << endl
527  << "INFO =>> [IBL] tot calib analysis..." << endl;
528 
529  TFile riTotFile(inTotFile, "READ");
530  TFile riTotFileAux(inTotFileAux, "READ");
531 
532  TString totHistName = "TOT_MEAN";
533  TString totSigHistName = "TOT_SIGMA";
534 
535  TDirectoryFile *scanDir = (TDirectoryFile *)((TKey *)riTotFile.GetListOfKeys()->First())->ReadObj();
536  TList *rodKeyList = (TList *)scanDir->GetListOfKeys();
537  TIter rodItr(rodKeyList);
538  TKey *rodKey;
539 
540  std::map<float, std::pair<vector<TString>, vector<Double_t>>> ModuDataToPrint;
541 
542 #if defined(DEMOXCHECK)
543  vector<TH1F *> h1_ChrgEntry(nToTibl, nullptr);
544 
545  vector<TH1F *> h1d_totSprdAll(nToTibl-1, nullptr);
546 
547  for (int t = 0; t < nToTibl; t++)
548  {
549  stringstream tt;
550  tt << t;
551 
552  string hname = "ChrgEntries_ToT_" + tt.str();
553  h1_ChrgEntry[t] = new TH1F(hname.c_str(), "RD entries ", 100, 0., 100000);
554 
555  hname = "totSprdsFll_ToT_" + tt.str();
556  if (t < nToTibl - 1)
557  h1d_totSprdAll[t] = new TH1F(hname.c_str(), "ToT spread ", 36, 0.2, 1.1);
558  }
559 #endif
560 
561  constexpr Int_t totFE = 14 * 16 * npsFEs; // 16 modules on each of 14 ReadOutDisk
562  using Row = Double_t[nchargeIBL];
563  using Row2 = Double_t[nToTibl];
564  //avoid stack overflow, create on heap
565  std::unique_ptr<Row[]> TotArray{ new Row[totFE]{} };
566  std::unique_ptr<Row[]> TotErrArray{ new Row[totFE]{} };
567  std::unique_ptr<Row[]> TotSigArray { new Row[totFE]{} };
568  std::unique_ptr<Row[]> TotSigErrArray{ new Row[totFE]{} };
569  std::unique_ptr<Row2[]> ChrgArray{ new Row2[totFE]{} };
570  std::unique_ptr<Row[]> ChrgErrArray{ new Row[totFE]{} };
571 
572  gRandom = new TRandom3(2203);
573  using Column = Double_t[totFE];
574  std::unique_ptr<Column[]> occuPhiEta{new Column[nchargeIBL]{}};
575 
576  Int_t cntRod = 0;
577  std::map<float, TString> devChrg_Order;
578  std::map<float, TString> devToT_Order;
579 
580  std::multimap<float, TString, std::greater<float>> badModules_Order;
581  std::multimap<float, TString, std::greater<float>> badModules_Order_detailed;
582 
583  while ((rodKey = (TKey *)rodItr()))
584  {
585  TString rodName(rodKey->GetName());
586  string rodStr(rodKey->GetName());
587  TDirectoryFile *rodDir = (TDirectoryFile *)rodKey->ReadObj();
588 
589  TString path = rodDir->GetPath();
590  TString pathAux = path.ReplaceAll(rodPath, rodPathAux);
591 
592  TDirectoryFile *rodDirAux = (TDirectoryFile *)riTotFileAux.Get(pathAux);
593  if (rodDirAux == NULL)
594  {
595  std::cout<< " Fail to get the rodPath in Aux: "<< pathAux << endl;
596  logout << " Fail to get the rodPath in Aux " << endl;
597  continue;
598  }
599 
600  TDirectory *dirRod = roTotDir->mkdir(rodName);
601 
602  TList *modKeyList = (TList *)rodDir->GetListOfKeys();
603  TIter modItr(modKeyList);
604 
605  TKey *modKey;
606 
607  int cntMod = 0;
608  string modNames[16 * npsFEs];
609 
610  float occuChrgs[nToTibl][16 * npsFEs];
611 
612 #if defined(DEMOXCHECK)
613  TString feName_maxDevChrg = "", feName_maxDevToT = "";
614  float maxDevChrg = -9., maxDevToT = -9., avgDevChrg = 0., avgDevToT = 0.;
615 
616  vector<TH1F *> h1d_totSprd(nToTibl - 1, nullptr);
617 
618  for (int t = 0; t < nToTibl - 1; t++)
619  {
620  stringstream tt;
621  tt << t + 1;
622 
623  string prfmodname = "ToT_Sprd_" + tt.str();
624  h1d_totSprd[t] = new TH1F(prfmodname.c_str(), "ToT spread ", 28, 0.3, 1.0);
625  prfmodname = rodStr + "ToT spread @ ToT_" + tt.str();
626  h1d_totSprd[t]->SetTitle(prfmodname.c_str());
627  }
628 #endif
629 
630  while ((modKey = (TKey *)modItr()))
631  {
632  TString modName(modKey->GetName());
633  string modStr(modKey->GetName());
634 
635  TDirectory *dirMod = dirRod->mkdir(modName);
636  bool ibl3Dfe0 = false, ibl3Dfe1 = false;
637 
638  int hashID = -1, hashIDL = -1, hashIDR = -1;
639  int bec = -1;
640  int layer = -1;
641  int phi_module = -1;
642  int eta_module = -1;
643  int eta_moduleL = -1, eta_moduleR = -1;
644  pixmap.mapping(modStr, &hashID, &bec, &layer, &phi_module, &eta_module);
645  if (hashID == -1)
646  {
647  pixmap.mapping(modStr + "_1", &hashID, &bec, &layer, &phi_module, &eta_module);
648  if (hashID != -1)
649  {
650  ibl3Dfe0 = true;
651  hashIDL = hashID;
652  eta_moduleL = eta_module;
653  }
654 
655  pixmap.mapping(modStr + "_2", &hashID, &bec, &layer, &phi_module, &eta_module);
656  if (hashID != -1)
657  {
658  ibl3Dfe1 = true;
659  hashIDR = hashID;
660  eta_moduleR = eta_module;
661  }
662  }
663  if (hashID == -1){
664  std::cout<< "negative hash ID\n";
665  continue;
666  }
667 
668 
669  if (pcdMap.find(modStr) == pcdMap.end()){
670  std::cout<< "MOD not found "<< modStr << "\n";
671  continue;
672  }
673 
674 
675  // a stack of Charge-ToT data for ToT-charge converting
676  array<std::unique_ptr<TH2D>, npsFEs> h2d_XchrgYtot{};
677  array<std::unique_ptr<TH2D>, npsFEs> h2d_XchrgYToTSig{};
678 
679  for (int sfe = 0; sfe < npsFEs; sfe++)
680  {
681  stringstream ss;
682  ss << sfe;
683 
684  modNames[cntMod * npsFEs + sfe] = modStr + "-" + ss.str();
685 
686  string prfmodname = modStr + "ToTvsChrg_FE" + ss.str();
687  h2d_XchrgYtot[sfe] = std::make_unique<TH2D>(prfmodname.c_str(), "ToT vs Charge", nchargeIBL, chrgsbins, nToTibl, totbins);
688 
689  prfmodname = modStr + "ToT_Sig_vsChrg_FE" + ss.str();
690  h2d_XchrgYToTSig[sfe] = std::make_unique<TH2D>(prfmodname.c_str(), "ToT Sig vs Charge", nchargeIBL, chrgsbins, 100, 0., 1.);
691  }
692 
693  // Fill ToT-Chrg for pixels per module
694  for (int c = 0; c < nchargeIBL; c++)
695  {
696  TString totHistDirPath = modName + "/" + totHistName + "/A0/B0/C";
697  totHistDirPath += std::to_string(c);
698 
699  std::unique_ptr<TH2F> h2dTot;
700  TDirectoryFile* totHistDir(static_cast<TDirectoryFile*> (rodDir->Get(totHistDirPath)));
701  if(!totHistDir){
702  std::cout<<" Missing totHistDir in : " << totHistDirPath << endl;
703  abort();
704  }
705  else {
706  h2dTot.reset(static_cast<TH2F*> ((static_cast<TKey*>(totHistDir->GetListOfKeys()->First()))->ReadObj()));
707  if (h2dTot) h2dTot->SetDirectory(0);
708  }
709  if (not h2dTot){
710  std::cout<<" Unrecoverable error in : " <<__LINE__ <<" of CalibrateIBL.cxx\n";
711  abort();
712  }
713 
714  unique_ptr<TH2F> h2dTotAux;
715  TDirectoryFile* totHistDirAux(static_cast<TDirectoryFile*>(rodDirAux->Get(totHistDirPath)));
716  if (!totHistDirAux)
717  {
718  std::cout<<" Missing totHistDir in : " << totHistDirPath << endl;
719  logout << " Missing totHistDir in : " << totHistDirPath << endl;
720  abort();
721  }
722  else {
723  h2dTotAux.reset(static_cast<TH2F*> ((static_cast<TKey*>(totHistDirAux->GetListOfKeys()->First())->ReadObj())));
724 
725  if(h2dTotAux) h2dTotAux->SetDirectory(0);
726  }
727  if (not h2dTotAux){
728  std::cout<<" Unrecoverable error in : " <<__LINE__ <<" of CalibrateIBL.cxx\n";
729  abort();
730  }
731  TString totSigHistDirPath = modName + "/" + totSigHistName + "/A0/B0/C";
732  totSigHistDirPath += std::to_string(c);
733  TDirectoryFile* totSigHistDir (static_cast<TDirectoryFile*>(rodDir->Get(totSigHistDirPath)));
734  unique_ptr<TH2F> h2dTotSig (static_cast<TH2F*> ((static_cast<TKey*>(totSigHistDir->GetListOfKeys()->First()))->ReadObj()));
735  h2dTotSig->SetDirectory(0);
736 
737  // Now one fill the calib-data from all pixles as basics.
738  for (int ieta = 0; ieta < ncol; ieta++)
739  {
740  for (int iphi = 0; iphi < nrow; iphi++)
741  {
742  int circ = -1;
743  if (moreFE)
744  {
745  // virtual FE with 80*84 pixels, so that sFE0, sFE1, sFE7, sFE6 will make up the physical 1'st FE
746  // while sFE2, sFE3, sFE5, sFE4 make up the physical 2'nd.
747  circ = (int)(iphi / 84.); // sFE0, sFE1, sFE2, sFE3
748  if (ieta > 80)
749  {
750  circ = 7 - circ;
751  } // sFE7, sFE6, sFE5, sFE4
752  }
753  else
754  {
755  if (ieta < ncol / 2)
756  circ = 0;
757  else
758  circ = 1;
759  }
760 
761  float tot = h2dTot->GetBinContent(ieta + 1, iphi + 1) + HDCshift;
762  float totAux = h2dTotAux->GetBinContent(ieta + 1, iphi + 1);
763 
764  // use ToT from higher HDC as default
765  float filltot = tot;
766 
767  // use ToT from lower HDC for SMALL charges, dirty in studying
768  // May 4, 2022, please keep below strategy, it give reasonable merging till now
769 
770  if (tot <= 5 || c <= 5)
771  {
772  if ((c == 0 && tot - totAux > 0.8) || (c >= 1 && c <= 3 && tot - totAux > 0.5) || (c > 4 && tot - totAux > 0.3))
773  filltot = totAux;
774  else
775  filltot = 0.5 * (totAux + tot);
776  if (c == 0 && totAux > 2.)
777  filltot = 1.;
778  }
779 
780  if (c > 19)
781  filltot = 0.5 * (totAux + tot + 2.);
782 
783  if (filltot > 0.){
784  h2d_XchrgYtot[circ]->Fill(chrgAbaciIBL[c], filltot, 1.);
785  }
786 
787 
788  float err = h2dTotSig->GetBinContent(ieta + 1, iphi + 1);
789  h2d_XchrgYToTSig[circ]->Fill(chrgAbaciIBL[c], err, 1.);
790  if (err == 0.)
791  err = 0.05;
792 
793 #if defined(DEMOXCHECK)
794  if (modName == XcheckModule && (iphi >= XcheckPhi[0] && iphi <= XcheckPhi[1]) && (ieta >= XcheckEta[0] && ieta <= XcheckEta[1]) && (c == XcheckCharge[0] || c == XcheckCharge[1] || (int)(tot + 0.5) == XcheckToT || (int)(totAux + 0.5) == XcheckToT))
795  logout << " tot = " << tot << ", totAux = " << totAux << " @injection " << chrgAbaciIBL[c] << " " << err << endl;
796 #endif
797  }
798  } // all pixels done over a module
799 
800  h2dTotSig.reset();
801  h2dTotAux.reset();
802  h2dTot.reset();
803 
804  } // end looping over charges
805  float modHash = -1;
806  if (!(ibl3Dfe0 || ibl3Dfe1))
807  {
808  modHash = hashID;
809  }
810 
811  array<TDirectory *, npsFEs> dirFE;
812 
813  Int_t idxMod = cntRod * 16 * npsFEs + cntMod * npsFEs;
814 
815  for (int sfe = 0; sfe < npsFEs; sfe++)
816  {
817  stringstream ss;
818  ss << sfe;
819 
820  if (ibl3Dfe0 && sfe < npsFEs / 2)
821  {
822  modHash = hashIDL;
823  }
824  else if (ibl3Dfe1 && sfe >= npsFEs / 2)
825  {
826  modHash = hashIDR;
827  }
828  else
829  modHash += (sfe < npsFEs / 2 ? 0 : 1) * 0.8;
830 
831  string prfmodname = "DirFE_" + ss.str();
832  dirFE[sfe] = dirMod->mkdir(prfmodname.c_str());
833 
834  // standard deviation will be adopted for ERROR
835  prfmodname = modStr + "_profile_Tot_FE" + ss.str();
836  std::unique_ptr<TProfile> prfl_TotsFE(h2d_XchrgYtot[sfe]->ProfileX(prfmodname.c_str(), 1, -1, "s"));
837  prfl_TotsFE->SetTitle(prfmodname.c_str());
838  prfmodname = modStr + "_profile_Chrg_FE" + ss.str();
839  // shanly : seeming free from ROOT-7770 issue. 2021 Sep 17
840  std::unique_ptr<TProfile> prfl_ChrgsFE(h2d_XchrgYtot[sfe]->ProfileY(prfmodname.c_str(), 1, -1, "s"));
841  prfl_ChrgsFE->SetTitle(prfmodname.c_str());
842 
843  Double_t TotArr[nchargeIBL], TotErrArr[nchargeIBL];
844  Double_t ChrgArr[nToTibl], ChrgErrArr[nToTibl];
845 
846  Int_t idxFE = idxMod + sfe;
847  // ToT distribution along nchargeIBL bins charges
848  int dumbFE = 0;
849  for (int c = 0; c < nchargeIBL; c++)
850  {
851  TotArr[c] = TotArray[idxFE][c] = prfl_TotsFE->GetBinContent(c + 1);
852  TotErrArr[c] = TotErrArray[idxFE][c] = prfl_TotsFE->GetBinError(c + 1);
853 
854  occuPhiEta[c][idxFE] = prfl_TotsFE->GetBinEntries(c + 1) * npsFEs / (1. * ncol * nrow);
855 
856  if ((1. - occuPhiEta[c][idxFE]) > badRDfracCut)
857  dumbFE++;
858 
859  stringstream sc;
860  sc << c;
861  prfmodname = modStr + "_ToTsig_FE" + ss.str() + "Chrg" + sc.str();
862 
863  std::unique_ptr<TH1D> h_ToTsig(h2d_XchrgYToTSig[sfe]->ProjectionY(prfmodname.c_str(), c + 1, c + 1));
864  TotSigArray[idxFE][c] = TMath::Sqrt(h_ToTsig->GetMean() * h_ToTsig->GetMean() + h_ToTsig->GetRMS() * h_ToTsig->GetRMS());
865  TotSigErrArray[idxFE][c] = TMath::Sqrt(h_ToTsig->GetMeanError() * h_ToTsig->GetMeanError() + h_ToTsig->GetRMSError() * h_ToTsig->GetRMSError());
866 
867  h_ToTsig.reset();
868  }
869  prfl_TotsFE.reset();
870 
871  if (dumbFE > 5)
872  {
873  for (int c = 0; c < nchargeIBL; c++)
874  {
875  TotArr[c] = TotArray[idxFE][c] = 0.5 * (TotArray[idxFE - 1][c] + TotArray[idxFE - 2][c]);
876  TotErrArr[c] = TotErrArray[idxFE][c] = 0.5 * (TotErrArray[idxFE - 1][c] + TotErrArray[idxFE - 2][c]);
877 
878  TotSigArray[idxFE][c] = 0.5 * (TotSigArray[idxFE - 1][c] + TotSigArray[idxFE - 2][c]);
879  TotSigErrArray[idxFE][c] = 0.5 * (TotSigErrArray[idxFE - 1][c] + TotSigErrArray[idxFE - 2][c]);
880 
881  }
882  }
883  // Charges distribution along 16 bins ToT
884  Double_t errToT_overChrg[nToTibl];
885  for (int t = 0; t < nToTibl; t++)
886  {
887  ChrgArr[t] = ChrgArray[idxFE][t] = prfl_ChrgsFE->GetBinContent(t + 1);
888  ChrgErrArr[t] = ChrgErrArray[idxFE][t] = prfl_ChrgsFE->GetBinError(t + 1);
889 
890  occuChrgs[t][cntMod * npsFEs + sfe] = prfl_ChrgsFE->GetBinEntries(t + 1) / RDentries[t];
891 
892  stringstream st;
893  st << t;
894  prfmodname = modStr + "_Chrg_FE" + ss.str() + "ToT" + st.str();
895  std::unique_ptr<TH1D> h_chrg(h2d_XchrgYtot[sfe]->ProjectionX(prfmodname.c_str(), t + 1, t + 1));
896  st.str("");
897  st << totAbaci[t];
898  prfmodname = modStr + " Charge_FE" + ss.str() + " @ ToT : " + st.str();
899  h_chrg->SetTitle(prfmodname.c_str());
900 
901  Double_t scl = h_chrg->Integral();
902  errToT_overChrg[t] = 0.;
903  if (scl > 0.)
904  for (int b = 0; b < nchargeIBL; b++)
905  {
906  errToT_overChrg[t] += TotErrArr[b] * (h_chrg->GetBinContent(b + 1)) / scl;
907  }
908  if (errToT_overChrg[t] == 0.)
909  errToT_overChrg[t] = 1.1;
910 
911  h_chrg.reset();
912  }
913  prfl_ChrgsFE.reset();
914 #if defined(DEMOXCHECK)
915  for (int b = 1; b < nToTibl; b++)
916  {
917  double err = errToT_overChrg[b];
918  h1d_totSprd[b - 1]->Fill(err);
919  h1d_totSprdAll[b - 1]->Fill(err);
920  }
921 #endif
922 
923  // try a correction upon dumb FE
924  if (dumbFE > 5)
925  {
926  for (int t = 0; t < nToTibl; t++)
927  {
928  ChrgArr[t] = ChrgArray[idxFE][t] = 0.5 * (ChrgArray[idxFE - 1][t] + ChrgArray[idxFE - 2][t]);
929  ChrgErrArr[t] = ChrgErrArray[idxFE][t] = 0.5 * (ChrgErrArray[idxFE - 1][t] + ChrgErrArray[idxFE - 2][t]);
930  }
931  }
932  bool xoticMod = std::find(reversedModules.begin(), reversedModules.end(), int(modHash)) != reversedModules.end();
933 
934  // when charges' order are distorted along ToTs
935  bool reverseH = false, reverseT = false;
936  for (int t = 4; t >= 1; t--)
937  if (ChrgArr[t] > ChrgArr[t + 1])
938  {
939  reverseH = true;
940  break;
941  }
942  for (int t = 4; t >= 1; t--)
943  {
944  int tl = nToTibl - t;
945  if (ChrgArr[tl] < ChrgArr[tl - 1])
946  {
947  reverseT = true;
948  break;
949  }
950  }
951 
952  if (reverseH)
953  {
954 #if defined(DEMOXCHECK)
955  logout << " Head reversed : " << int(modHash) << " " << modName << " " << sfe << " " << ChrgArr[1] << " "
956  << ChrgArr[2] << " " << ChrgArr[3] << " " << ChrgArr[4] << " " << ChrgArr[5] << std::endl;
957 #endif
958  // try a correction with some averaged factors
959  for (int t = 4; t >= 1; t--)
960  if (ChrgArr[t] > ChrgArr[t + 1] && (occuChrgs[t][cntMod * npsFEs + sfe] < 0.5 || xoticMod))
961  ChrgArr[t] = ChrgArray[idxFE][t] = reverseCF_H[t - 1] * ChrgArr[t + 1];
962  }
963  else // prepare the average factors along ToTs
964  {
965  for (int t = 4; t >= 1; t--)
966  if (ChrgArr[t + 1] != 0.)
967  reverseCF_H[t - 1] = (reverseCF_H[t - 1] == tmpRCf_H[t - 1] ? ChrgArr[t] / ChrgArr[t + 1] : 0.5 * (reverseCF_H[t - 1] + ChrgArr[t] / ChrgArr[t + 1]));
968  }
969 
970  if (reverseT)
971  {
972 #if defined(DEMOXCHECK)
973  logout << " Tail reversed : " << int(modHash) << " " << modName << " " << sfe << " " << ChrgArr[nToTibl - 5] << " "
974  << ChrgArr[nToTibl - 4] << " " << ChrgArr[nToTibl - 3] << " " << ChrgArr[nToTibl - 2] << " " << ChrgArr[nToTibl - 1] << std::endl;
975 #endif
976  // try a correction with some averaged factors
977  for (int t = 4; t >= 1; t--)
978  {
979  int tl = nToTibl - t;
980  if (ChrgArr[tl] < ChrgArr[tl - 1] && (occuChrgs[tl][cntMod * npsFEs + sfe] < 0.5 || xoticMod))
981  ChrgArr[tl] = ChrgArray[idxFE][tl] = reverseCF_T[t - 1] * ChrgArr[tl - 1];
982  }
983  }
984  else // prepare the average factors along ToTs
985  {
986  for (int t = 4; t >= 1; t--)
987  {
988  int tl = nToTibl - t;
989  if (ChrgArr[tl - 1] != 0.)
990  reverseCF_T[t - 1] = (reverseCF_T[t - 1] == tmpRCf_T[t - 1] ? ChrgArr[tl] / ChrgArr[tl - 1] : 0.5 * (reverseCF_T[t - 1] + ChrgArr[tl] / ChrgArr[tl - 1]));
991  }
992  }
993 
994  string gername = modStr + "_grToTsig_FE_" + ss.str();
995  std::unique_ptr<TGraphErrors> grTotSig = std::make_unique<TGraphErrors>(nchargeIBL, chrgAbaciIBL, TotSigArray[sfe], chargeErrArrIBL, TotSigErrArray[sfe]);
996  grTotSig->SetTitle(gername.c_str());
997  grTotSig->SetName(gername.c_str());
998 
999  std::unique_ptr<TGraph> grTotSprd = std::make_unique<TGraph>(nToTibl, totAbaci, errToT_overChrg);
1000  gername = "ToTsprd";
1001  grTotSprd->SetName(gername.c_str());
1002  gername = modStr + " ToT spread over FrontEnd @ FE " + ss.str();
1003  grTotSprd->SetTitle(gername.c_str());
1004  grTotSprd->SetLineColor(4);
1005  grTotSprd->SetMarkerStyle(4);
1006  grTotSprd->SetMarkerSize(1.);
1007 
1008  TF1 ToTres("ToTres", funcDispConvoluted, 0.8, 16.5, 3);
1009  ToTres.SetParameter(1,13.);
1010  ToTres.SetParameter(2,0.65);
1011  ToTres.SetParLimits(1, 10., 16.);
1012  ToTres.SetParLimits(2, 0.3, 1.0);
1013 
1014  grTotSprd->Fit(&ToTres, "MRQ");
1015  Double_t parP0 = ToTres.GetParameter(0);
1016  Double_t parP1 = ToTres.GetParameter(1);
1017 
1018  dirFE[sfe]->WriteTObject(grTotSprd.get());
1019  grTotSprd.reset();
1020  // Lets' keep the bad 2-para fit untill the dataBase will be updated someday
1043  gername = modStr + "_grToT_FE_" + ss.str();
1044  std::unique_ptr<TGraphErrors> grTot = std::make_unique<TGraphErrors>(nchargeIBL, chrgAbaciIBL, TotArr, chargeErrArrIBL, TotErrArr);
1045  grTot->SetTitle(gername.c_str());
1046  grTot->SetName(gername.c_str());
1047  TF1 *f1ToTfromCharge = new TF1("ToTfromCharge", funcRation5, chrgAbaciIBL[0] - 100., chrgAbaciIBL[nchargeIBL - 1] + 300., 5);
1048 
1049  grTot->Fit(f1ToTfromCharge, "MRQ");
1050 
1051  std::unique_ptr<TGraphErrors> grChrg = std::make_unique<TGraphErrors>(nToTibl, totAbaci, ChrgArr, totErrArr, ChrgErrArr);
1052 
1053  gername = modStr + "_grChrg_FE_" + ss.str();
1054  grChrg->SetName(gername.c_str());
1055  grChrg->SetTitle(gername.c_str());
1056  grChrg->SetLineColor(3);
1057  grChrg->SetLineWidth(3);
1058 
1059  TF1 *f1ChargefromToTLeft = new TF1("ChargefromToTL", funcRation5, 1., 7.5, 5);
1060  f1ChargefromToTLeft->SetLineColor(2);
1061  grChrg->Fit(f1ChargefromToTLeft, "MRQ");
1062  if (f1ChargefromToTLeft->GetChisquare() / 4. > 2.){
1063  logout << "bad fit Left... " << std::endl;
1064  }
1065 
1066  TF1 *f1ChargefromToTRight = new TF1("ChargefromToTR", funcRation5, 5.5, 16., 5);
1067  f1ChargefromToTRight->SetLineColor(4);
1068  grChrg->Fit(f1ChargefromToTRight, "MRQ+");
1069 
1070  dirFE[sfe]->WriteTObject(grTot.get());
1071  dirFE[sfe]->WriteTObject(grTotSig.get());
1072  dirFE[sfe]->WriteTObject(h2d_XchrgYtot[sfe].get());
1073 
1074  h2d_XchrgYtot[sfe].reset();
1075  h2d_XchrgYToTSig[sfe].reset();
1076  grTot.reset();
1077  grTotSig.reset();
1078  grChrg.reset();
1079 
1080  // prepare a map sor sorting
1081  string Idx = "I" + ss.str();
1082  vector<TString> modName2prt;
1083  modName2prt.push_back(modName);
1084  modName2prt.push_back((TString)(Idx));
1085  // unfortunately some tuning is necessary to satisfy dataBase reauirements
1086  // 5 for module ID, 4 for threshold, nToTibl - 1 for charges at ToT, 3 for ResToT
1087  std::vector<Double_t> prtAux;
1088  prtAux.reserve(5 + 4 + nToTibl - 1 + 3); // skip ToT = 0
1089  // @ 0 ... 4
1090  prtAux.push_back(modHash);
1091  prtAux.push_back(bec);
1092  prtAux.push_back(layer);
1093  prtAux.push_back(phi_module);
1094  prtAux.push_back(eta_module);
1095 
1096  modName2prt[0] = (modStr + "_" + ss.str()).c_str();
1097  if (ibl3Dfe0 && sfe < npsFEs / 2)
1098  {
1099  if (!moreFE)
1100  modName2prt[0] = modName + "_0";
1101  prtAux[0] = hashIDL;
1102  prtAux[4] = eta_moduleL;
1103  }
1104  if (ibl3Dfe1 && sfe >= npsFEs / 2)
1105  {
1106  if (!moreFE)
1107  modName2prt[0] = modName + "_1";
1108  prtAux[0] = hashIDR;
1109  prtAux[4] = eta_moduleR;
1110  }
1111 
1112  // @ 5 ... 8
1113  float ThrNorm = pcdMap[modStr][Idx]["ThrNorm"];
1114  if (ThrNorm == 0. || ThrNorm == -42.)
1115  {
1116  // retrieve something from average when it's missing in calibration
1117  prtAux.push_back(THR_avg[0][sfe]);
1118  prtAux.push_back(ThrSig_avg[0][sfe]);
1119  prtAux.push_back(THR_avg[1][sfe]);
1120  prtAux.push_back(ThrSig_avg[1][sfe]);
1121  }
1122  else
1123  {
1124  prtAux.push_back(pcdMap[modStr][Idx]["ThrNorm"]);
1125  prtAux.push_back(pcdMap[modStr][Idx]["ThrSigNorm"]);
1126  prtAux.push_back(pcdMap[modStr][Idx]["ThrLong"]);
1127  prtAux.push_back(pcdMap[modStr][Idx]["ThrSigLong"]);
1128  }
1129  // since this version, the goodness will be skipped
1130  // @ 10 ... nToTibl + 8, with ToT = 0 skipped
1131  for (int t = 1; t < nToTibl; t++)
1132  prtAux.push_back(ChrgArr[t]);
1133 
1134  // @ nToTibl + 9, nToTibl + 10
1135  prtAux.push_back(parP0);
1136  prtAux.push_back(parP1);
1137 
1138  std::pair<vector<TString>, vector<Double_t>> payloadDB = std::pair<vector<TString>, vector<Double_t>>(modName2prt, prtAux);
1139  ModuDataToPrint.insert(std::pair<float, std::pair<vector<TString>, vector<Double_t>>>(modHash, payloadDB));
1140  } // end 1'st loop over FrontEnds
1141 
1142  cntMod++;
1143  } // end of loop over Mods
1144 
1145 #if defined(DEMOXCHECK)
1146 
1147  TDirectory *dirToTSprd = dirRod->mkdir("ToT_Spreads");
1148  for (int b = 0; b < nToTibl - 1; b++)
1149  {
1150  Int_t upper = h1d_totSprd[b]->FindLastBinAbove(1, 1);
1151  h1d_totSprd[b]->SetAxisRange(h1d_totSprd[b]->GetBinLowEdge(
1152  h1d_totSprd[b]->FindFirstBinAbove(1, 1)),
1153  h1d_totSprd[b]->GetBinLowEdge(upper) + h1d_totSprd[b]->GetBinWidth(upper), "X");
1154 
1155  dirToTSprd->WriteTObject(h1d_totSprd[b]);
1156  delete h1d_totSprd[b];
1157  }
1158  h1d_totSprd.clear();
1159 
1160  TString hn = "";
1161 
1162  TH2F *h2_badChrgs = new TH2F(rodName + "ChargeOccupancy", "Occupancy ", nchargeIBL, 0, nchargeIBL, 16 * npsFEs, 0, 16 * npsFEs);
1163  hn = "Occupancy along charges @ " + rodName;
1164  h2_badChrgs->SetTitle(hn);
1165 
1166  TH2F *h2_badToT = new TH2F(rodName + "ToTOccupancy", "Occupancy ", nToTibl - 1, 0.5, nToTibl - 0.5, 16 * npsFEs, 0, 16 * npsFEs);
1167  hn = "Occupancy along ToT @ " + rodName;
1168  h2_badToT->SetTitle(hn);
1169 
1170  bool fillChrg = false, fillToT = false;
1171 
1172  for (int sfe = 0; sfe < 16 * npsFEs; sfe++)
1173  {
1174  float PhiEta = 0.;
1175  for (int c = 0; c < nchargeIBL; c++)
1176  {
1177  float occu = occuPhiEta[c][cntRod * 16 * npsFEs + sfe];
1178  float nt = 10. - floor(10. * occu + 0.5);
1179  PhiEta += 0.1 * nt;
1180  TString modchrg = (TString)(modNames[sfe] + "_Chrg_" + c);
1181  if (nt > 0.01)
1182  badModules_Order_detailed.insert(std::pair<float, TString>(1. - occu, modchrg));
1183  if (nt > 5.)
1184  continue;
1185  if (nt < 0.1)
1186  nt = 0.01;
1187 
1188  h2_badChrgs->SetBinContent(c + 1, sfe + 1, nt);
1189  fillChrg = true;
1190  }
1191  PhiEta /= nchargeIBL;
1192  if (PhiEta > 0.01)
1193  badModules_Order.insert(std::pair<float, TString>(PhiEta,
1194  (TString)(modNames[sfe])));
1195 
1196  for (int t = 1; t < nToTibl; t++)
1197  {
1198  float nt = floor(10. * occuChrgs[t][sfe] + 0.5);
1199  if (nt > 5.)
1200  continue;
1201  h2_badToT->SetBinContent(t + 1, sfe + 1, nt);
1202  fillToT = true;
1203  }
1204 
1205  TString shortName = modNames[sfe].substr(9, 13).c_str();
1206  h2_badChrgs->GetYaxis()->SetBinLabel(sfe + 1, shortName);
1207  h2_badToT->GetYaxis()->SetBinLabel(sfe + 1, shortName);
1208  h2_badChrgs->SetStats(kFALSE);
1209  h2_badToT->SetStats(kFALSE);
1210  h2_badChrgs->SetTickLength(0.01, "Y");
1211  h2_badToT->SetTickLength(0.01, "Y");
1212  }
1213 
1214  h2_badChrgs->SetOption("TEXT");
1215  h2_badToT->SetOption("TEXT");
1216 
1217  if (fillChrg)
1218  {
1219  dirRod->WriteTObject(h2_badChrgs, "", "Overwrite");
1220  }
1221 
1222  if (fillToT)
1223  {
1224  dirRod->WriteTObject(h2_badToT, "", "WriteDelete");
1225  }
1226 
1227  delete h2_badChrgs;
1228  delete h2_badToT;
1229 
1230 #endif
1231 
1232  cntRod++;
1233  } // end of loop over Rods
1234 
1235 #if defined(DEMOXCHECK)
1236 
1237  for (int t = 0; t < 4; t++)
1238  {
1239  if (t == 0)
1240  logout << " correction factors for reversed charges : ";
1241  logout << reverseCF_H[t] << ", ";
1242  if (t == 3)
1243  logout << endl;
1244  }
1245  for (int t = 0; t < 4; t++)
1246  {
1247  if (t == 0)
1248  logout << " correction factors for reversed charges : ";
1249  logout << reverseCF_T[t] << ", ";
1250  if (t == 3)
1251  logout << endl;
1252  }
1253 
1254  for (int t = 0; t < nToTibl; t++)
1255  {
1256  if (t == 0)
1257  logout << "RD Entries per ToT bin ";
1258  logout << h1_ChrgEntry[t]->GetMean() << ", ";
1259  if (t == nToTibl - 1)
1260  logout << std::endl;
1261  roTotDir->WriteTObject(h1_ChrgEntry[t]);
1262  delete h1_ChrgEntry[t];
1263  }
1264  h1_ChrgEntry.clear();
1265 
1266  for (std::map<float, TString>::const_iterator itr = devChrg_Order.begin();
1267  itr != devChrg_Order.end(); ++itr)
1268  logout << " Charge dev order : " << itr->second << " : " << itr->first << endl;
1269 
1270  for (std::map<float, TString>::const_iterator itr = devToT_Order.begin();
1271  itr != devToT_Order.end(); ++itr)
1272  logout << " ToT dev order : " << itr->second << " : " << itr->first << endl;
1273 
1274  logout << "ToTRes (spread) & its RMS @ ToT from 1 to 16 : " << endl;
1275  for (int t = 0; t < nToTibl - 1; t++)
1276  {
1277  if (t == 0)
1278  logout << "[ ";
1279  logout << "[ " << h1d_totSprdAll[t]->GetMean() << " , " << h1d_totSprdAll[t]->GetStdDev() << " ] ";
1280  if (t < nToTibl - 2)
1281  logout << " , " << endl;
1282  if (t == nToTibl - 2)
1283  logout << " ] " << endl;
1284  delete h1d_totSprdAll[t];
1285  }
1286  h1d_totSprdAll.clear();
1287 #endif
1288 
1289  logout << " goto print " << std::endl;
1290 
1291  // Comment out for information only
1292  // txtDB << "####################################################################################################" << endl;
1293  // txtDB << "=========== physics meaning of each column ====================" << endl;
1294  // txtDB << "hashId " << "Threshold " << "ThresholdRMS " << "ThresholdLong " << "ThresholdRMSLong " << "Chrg_ToT1 "
1295  // << "Chrg_ToT2 " << "Chrg_ToT3 " << "Chrg_ToT4 " << "Chrg_ToT5 " << "Chrg_ToT6 " << "Chrg_ToT7 " << "Chrg_ToT8 " << "Chrg_ToT9 "
1296  // << "Chrg_ToT10 " << "Chrg_ToT11 " << "Chrg_ToT12 " << "Chrg_ToT13 " << "Chrg_ToT14 " << "Chrg_ToT15 " << "Chrg_ToT16 "
1297  // << "ToT_sprd " << " ToT_sprd_var" << endl;
1298 
1299  // txtDB << "######################## TXT for IBL calibration ____Start : _> ###########################" << endl;
1300 
1301  // now print the dataBase payloads in a TXT format
1302  for (std::map<float, std::pair<vector<TString>, vector<Double_t>>>::const_iterator itr = ModuDataToPrint.begin();
1303  itr != ModuDataToPrint.end(); ++itr)
1304  {
1305  int hash = floor(itr->first);
1306  std::pair<vector<TString>, vector<Double_t>> payload = itr->second;
1307  vector<Double_t> fe = payload.second;
1308 
1309  txtDB << hash << " ";
1310 
1311  for (unsigned int t = 5; t < fe.size(); t++)
1312  {
1313  if (t < fe.size() - 2)
1314  {
1315  txtDB << (int)(fe[t]) << " ";
1316 #if defined(DEMOXCHECK)
1317  if (fe[t] > fe[t + 1] && t >= 9 && t <= 23)
1318  logout << " Reversed charge ! " << t - 5 << " " << fe[t] << " " << fe[t + 1] << std::endl;
1319 #endif
1320  }
1321  else
1322  {
1323  float p01 = abs(fe[t]);
1324  if (t < fe.size() - 1)
1325  txtDB << 0.001 * (int)(p01 * 1000.) << " ";
1326  else
1327  txtDB << 0.000000001 * (int)(p01 * 1000000000.) << " ";
1328  }
1329  }
1330  txtDB << endl;
1331  }
1332  // txtDB << "####################################################################################################" << endl;
1333  // txtDB << "######################## TXT for IBL calibration <___ : End _____ ###########################" << endl;
1334  // txtDB << "####################################################################################################" << endl;
1335 
1336  std::cout << " Please find the file : " << dbFileName << " for dataBase payload " << endl;
1337 
1338  // statistics for BAD frontends, only for ELOG reports
1339  logout << " modules lacking in RD during Threshold scan : " << std::endl;
1340  for (std::multimap<float, TString, std::greater<float>>::const_iterator itr = badThr_Order.begin(); itr != badThr_Order.end(); ++itr)
1341  logout << " " << itr->second << " " << (itr->first) * 100. << "%" << endl;
1342 
1343  logout << " modules lacking in RD during ToT scan : " << std::endl;
1344  for (std::multimap<float, TString, std::greater<float>>::const_iterator itr = badModules_Order.begin(); itr != badModules_Order.end(); ++itr)
1345  logout << " " << itr->second << " " << (itr->first) * 100. << "%" << endl;
1346 
1347  logout << " modules lacking in RD at certain charges during ToT scan : " << std::endl;
1348  for (std::multimap<float, TString, std::greater<float>>::const_iterator itr = badModules_Order_detailed.begin();
1349  itr != badModules_Order_detailed.end(); ++itr)
1350  logout << " " << itr->second << " " << (itr->first) * 100. << "%" << endl;
1351 
1352  roFile.Close();
1353  txtDB.close();
1354  logout.close();
1355 
1356  return 0;
1357 }

◆ main()

int main ( int  argc,
char *  argv[] 
)

Definition at line 1365 of file CalibrateIBL.cxx.

1365  {
1366 
1367  bool saveInfo = false;
1368  std::string THR = "THR";
1369  std::string TOT_HISDIS = "TOT_HISDIS";
1370  std::string TOT_LOWQ = "TOT_LOWQ";
1371  std::string dpath = "directory_path";
1372 
1373  for(int i=1; i<argc; i++){
1374  std::string aux(argv[i]);
1375  if (THR.compare(aux.substr(0,aux.find("="))) == 0) THR = aux.substr(aux.find("=")+1);
1376  else if(TOT_HISDIS.compare(aux.substr(0,aux.find("="))) == 0) TOT_HISDIS = aux.substr(aux.find("=")+1);
1377  else if(TOT_LOWQ.compare(aux.substr(0,aux.find("="))) == 0) TOT_LOWQ = aux.substr(aux.find("=")+1);
1378  else if(dpath.compare(aux.substr(0,aux.find("="))) == 0) dpath = aux.substr(aux.find("=")+1);
1379  else if(aux.compare("saveInfo") == 0) saveInfo = true;
1380  else{
1381  printError();
1382  return 1;
1383  }
1384  }
1385 
1386  printf("%-14s = %s\n","Directory path",dpath.c_str());
1387  printf("%-14s = %s.root\n","THR",THR.c_str());
1388  printf("%-14s = %s.root\n","TOT_HISDIS",TOT_HISDIS.c_str());
1389  printf("%-14s = %s.root\n","TOT_LOWQ",TOT_LOWQ.c_str());
1390  printf("%-14s = %s\n\n\n" ,"Save root file",saveInfo ? "True" : "False" );
1391 
1392  bool correctArgc = (THR.compare("THR") == 0) or (TOT_HISDIS.compare("TOT_HISDIS") == 0) or (TOT_LOWQ.compare("TOT_LOWQ") == 0) or (dpath.compare("directory_path") == 0);
1393 
1394  if(correctArgc){
1395  printf("Cannot continue, one arguments is incorrect or not filled correctly...\n");
1396  printf("Helper below:\n**********************\n\n");
1397  printError();
1398  return 1;
1399  }
1400 
1401  iblCalib(dpath, THR, TOT_HISDIS, TOT_LOWQ);
1402  return 0;
1403 
1404 }

◆ printError()

void printError ( )

Definition at line 1359 of file CalibrateIBL.cxx.

1359  {
1360  printf("ERROR - Argument not expected or wrongly set:\n\n");
1361  printf("Valid format is: ./IBLCalibration THR=SCAN_Sxxxxxxxxx TOT_HISDIS=SCAN_Sxxxxxxxxx TOT_LOWQ=SCAN_Sxxxxxxxxx directory_path=path/to/file/\n");
1362  printf("\n\t i.e: ./IBLCalibration THR=SCAN_S000087719 TOT_HISDIS=SCAN_S000087717 TOT_LOWQ=SCAN_S000087710 directory_path=/eos/user/x/xxxx/\n");
1363 }

◆ ToTfromChrg()

double ToTfromChrg ( float  chrg,
double  par[5] 
)

Definition at line 82 of file CalibrateIBL.cxx.

83 {
84  double tot = 0.;
85  double num = par[0] + chrg * (par[1] + chrg * (par[2] + chrg * par[3]));
86  double denom = 1. + chrg * par[4];
87  if (denom != 0.0)
88  tot = num / denom;
89  return tot;
90 }

Variable Documentation

◆ ATLAS_NO_CHECK_FILE_THREAD_SAFETY

ATLAS_NO_CHECK_FILE_THREAD_SAFETY

Definition at line 24 of file CalibrateIBL.cxx.

printError
void printError()
Definition: CalibrateIBL.cxx:1359
athena.path
path
python interpreter configuration --------------------------------------—
Definition: athena.py:128
Merge_LAr_blobs.filled
list filled
Definition: Merge_LAr_blobs.py:43
PowhegControl_ttHplus_NLO.ss
ss
Definition: PowhegControl_ttHplus_NLO.py:83
find
std::string find(const std::string &s)
return a remapped string
Definition: hcg.cxx:135
ParticleGun_SamplingFraction.bec
int bec
Definition: ParticleGun_SamplingFraction.py:89
dqt_zlumi_pandas.hname
string hname
Definition: dqt_zlumi_pandas.py:280
CaloCondBlobAlgs_fillNoiseFromASCII.spec
spec
Definition: CaloCondBlobAlgs_fillNoiseFromASCII.py:46
SH::scanDir
void scanDir(SampleHandler &sh, DiskList &list, const std::string &pattern, const std::string &samplePattern, const std::string &samplePostfix)
effects: scan the given directory and add all subdirectories as samples that contain root files.
Definition: ToolsDiscovery.cxx:45
GetMean
TGraphErrors * GetMean(TH2F *histo)
Definition: TRTCalib_makeplots.cxx:3913
keylayer_zslicemap.row
row
Definition: keylayer_zslicemap.py:155
upper
int upper(int c)
Definition: LArBadChannelParser.cxx:49
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
drawFromPickle.exp
exp
Definition: drawFromPickle.py:36
python.TrigEgammaMonitorHelper.TH2F
def TH2F(name, title, nxbins, bins_par2, bins_par3, bins_par4, bins_par5=None, bins_par6=None, path='', **kwargs)
Definition: TrigEgammaMonitorHelper.py:45
x
#define x
AthenaPoolTestRead.sc
sc
Definition: AthenaPoolTestRead.py:27
calibdata.valid
list valid
Definition: calibdata.py:44
IBLCalibrationConfig.thr
thr
Definition: IBLCalibrationConfig.py:39
dqt_zlumi_pandas.err
err
Definition: dqt_zlumi_pandas.py:183
lumiFormat.i
int i
Definition: lumiFormat.py:85
LArCellNtuple.argv
argv
Definition: LArCellNtuple.py:152
python.getProblemFolderFromLogs.st
st
Definition: getProblemFolderFromLogs.py:68
TRT::Hit::layer
@ layer
Definition: HitInfo.h:79
chi2
double chi2(TH1 *h0, TH1 *h1)
Definition: comparitor.cxx:525
python.BuildSignatureFlags.sig
sig
Definition: BuildSignatureFlags.py:237
DQHistogramMergeRegExp.argc
argc
Definition: DQHistogramMergeRegExp.py:19
dqt_zlumi_pandas.run3
string run3
Definition: dqt_zlumi_pandas.py:41
compute_lumi.denom
denom
Definition: compute_lumi.py:76
pix::PixelMapping
Definition: PixelMapping.h:18
trigbs_pickEvents.num
num
Definition: trigbs_pickEvents.py:76
PixelCalibrationConfig.tot
tot
Definition: PixelCalibrationConfig.py:28
createCoolChannelIdFile.par
par
Definition: createCoolChannelIdFile.py:28
ActsTrk::to_string
std::string to_string(const DetectorType &type)
Definition: GeometryDefs.h:34
plotBeamSpotMon.b
b
Definition: plotBeamSpotMon.py:76
funcRation5
double funcRation5(double *x, double *par)
Definition: CalibrateIBL.cxx:92
PixelModuleFeMask_create_db.payload
string payload
Definition: PixelModuleFeMask_create_db.py:69
python.CaloAddPedShiftConfig.int
int
Definition: CaloAddPedShiftConfig.py:45
scan
void scan(TDirectory *td=0, int depth=0)
Definition: listroot.cxx:440
CaloCondBlobAlgs_fillNoiseFromASCII.hash
dictionary hash
Definition: CaloCondBlobAlgs_fillNoiseFromASCII.py:108
PathResolver::find_file
static std::string find_file(const std::string &logical_file_name, const std::string &search_path)
Definition: PathResolver.cxx:221
BchCleanup.modName
modName
Definition: BchCleanup.py:201
get
T * get(TKey *tobj)
get a TObject* from a TKey* (why can't a TObject be a TKey?)
Definition: hcg.cxx:127
LArNewCalib_DelayDump_OFC_Cali.idx
idx
Definition: LArNewCalib_DelayDump_OFC_Cali.py:69
DiTauMassTools::MMCFitMethod::shortName
const std::string shortName[MAX]
Definition: PhysicsAnalysis/TauID/DiTauMassTools/DiTauMassTools/HelperFunctions.h:48
dqt_zlumi_alleff_HIST.tl
tl
Definition: dqt_zlumi_alleff_HIST.py:73
python.TrigEgammaMonitorHelper.TH1F
def TH1F(name, title, nxbins, bins_par2, bins_par3=None, path='', **kwargs)
Definition: TrigEgammaMonitorHelper.py:24
beamspotnt.nt
def nt
Definition: bin/beamspotnt.py:1062
skip
bool skip
Definition: TrigGlobEffCorrValidation.cxx:190
TileDCSDataPlotter.tt
tt
Definition: TileDCSDataPlotter.py:876
python.compressB64.c
def c
Definition: compressB64.py:93
funcDispConvoluted
double funcDispConvoluted(double *x, double *par)
Definition: CalibrateIBL.cxx:102
iblCalib
int iblCalib(const std::string &InDir, const std::string &THRscan, const std::string &scan, const std::string &scanLowCharge)
Definition: CalibrateIBL.cxx:109
Generated on Sun Oct 26 2025 21:23:45 for ATLAS Offline Software by   doxygen 1.8.18