ATLAS Offline Software
TRTCalibrator.cxx
Go to the documentation of this file.
1 
2 /*
3  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
4 */
5 
6 /********************************************************************
7 
8  NAME: TRTCalibrator.cxx
9 PACKAGE: TRTCalibTools
10 
11 AUTHORS: Johan Lundquist
12 CREATED: 27-03-2009
13 
14 PURPOSE: Tool for calibrating the TRT
15 
16 ********************************************************************/
17 
18 #include "TRTCalibrator.h"
19 #include "Calibrator.h"
20 #include "TRT_CalibData/TrackInfo.h"
21 #include "TrkTrack/Track.h"
22 #include "TrkParameters/TrackParameters.h"
23 #include "TrkEventPrimitives/LocalParameters.h"
24 #include "EventPrimitives/EventPrimitives.h"
25 
26 #include "TrkRIO_OnTrack/RIO_OnTrack.h"
27 #include "InDetRIO_OnTrack/TRT_DriftCircleOnTrack.h"
28 #include "InDetPrepRawData/TRT_DriftCircle.h"
29 #include "TrkPrepRawData/PrepRawData.h"
30 #include "AtlasDetDescr/AtlasDetectorID.h"
31 #include "InDetIdentifier/TRT_ID.h"
32 #include "TRT_ConditionsData/RtRelation.h"
33 #include "TRT_ConditionsData/BasicRtRelation.h"
34 #include "TRT_ConditionsData/DinesRtRelation.h"
35 #include "TRT_ConditionsData/BinnedRtRelation.h"
36 #include "CommissionEvent/ComTime.h"
37 #include "TRT_ReadoutGeometry/TRT_DetectorManager.h"
38 
39 #include <TNtuple.h>
40 #include <TFile.h>
41 #include <fstream>
42 
43 
44 TRTCalibrator::TRTCalibrator(const std::string& type, const std::string& name, const IInterface* parent) :
45  AthAlgTool(type, name, parent),
46  m_DetID(nullptr),
47  m_TRTID(nullptr),
48  m_trtmanager(nullptr)
49 {
50  declareInterface<ITRTCalibrator>(this);
51 }
52 
53 
54 StatusCode TRTCalibrator::initialize()
55 {
56  ATH_MSG_INFO( "initialize()" );
57 
58  m_trtmanager=nullptr;
59 
60  if ((detStore()->retrieve(m_DetID,"AtlasID")).isFailure()) {
61  ATH_MSG_FATAL( "Problem retrieving ATLASDetectorID helper" );
62  return StatusCode::FAILURE;
63  }
64  if ((detStore()->retrieve(m_trtmanager,m_ntrtmanager)).isFailure()) {
65  ATH_MSG_FATAL ( "Could not get TRT_DetectorManager" );
66  return StatusCode::FAILURE;
67  }
68  if ((detStore()->retrieve(m_TRTID)).isFailure()) {
69  ATH_MSG_FATAL( "Problem retrieving TRTID helper" );
70  return StatusCode::FAILURE;
71  }
72  if(m_trtcaldbTool.retrieve().isFailure()) {
73  ATH_MSG_FATAL( "Could not get TRTCalDbTool !" );
74  return StatusCode::FAILURE;
75  }
76  if(StatusCode::SUCCESS!=m_neighbourSvc.retrieve() ) {
77  ATH_MSG_FATAL("Could not get TRTStrawNeighbourTool !");
78  return StatusCode::FAILURE;
79  }
80 
81 
82  // The tool to get the argon status:
83 
84  if (m_TRTStrawSummaryTool.retrieve().isFailure() ) {
85  ATH_MSG_ERROR ("Failed to retrieve StrawStatus Summary " << m_TRTStrawSummaryTool);
86  return StatusCode::FAILURE;
87  } else {
88  ATH_MSG_INFO( "Retrieved tool " << m_TRTStrawSummaryTool );
89  }
90 
91 
92 
93  /*
94  Pre-define standard calibration configuration
95 
96  The string-vectors "CalibrateRt" ("CalibrateT0") contains the names of the levels for
97  which R-t (T0) calibration should be performed
98 
99  The string-vector "PrintLog" contains the names of the levels for which a status line
100  should be printed to the output log
101 
102  The string-vectors "PrintT0Out" and "PrintRtOut" contains the names of the levels for
103  text files containing the calibration constants should be generated
104 
105  This standard configuration divides the Barrel in 9 parts:
106 
107  "all": a single calibration for all hits in the whole TRT
108  "-1" & "1": calibration of the barrel C- & A-side
109  "X_Y": calibration of layer Y in barrel X-side
110  */
111 
112  //the whole TRT
113  m_config["all"].selections["Detector"].insert(CALIB_NONE); // don't go deeper than TRT level
114  m_config["all"].selections["Layer"].insert(CALIB_NONE);
115  m_config["all"].selections["Module"].insert(CALIB_NONE);
116  m_config["all"].selections["Board"].insert(CALIB_NONE);
117  m_config["all"].selections["Chip"].insert(CALIB_NONE);
118  m_config["all"].selections["Straw"].insert(CALIB_NONE);
119  m_config["all"].CalibrateRt.emplace_back("TRT") ; //only calibrate R-t and T0 on TRT level
120  m_config["all"].CalibrateT0.emplace_back("TRT") ;
121  m_config["all"].PrintLog.emplace_back("TRT") ;
122  m_config["all"].PrintT0Out.emplace_back("TRT") ;
123  m_config["all"].PrintRtOut.emplace_back("TRT") ;
124 
125  //common for all barrel/ec
126  m_config["subdet"].CalibrateRt.emplace_back("TRT") ; //calibrate R-t and T0 on TRT and barrel/ec level
127  m_config["subdet"].CalibrateRt.emplace_back("Detector") ;
128  m_config["subdet"].CalibrateT0=m_config["subdet"].CalibrateRt;
129  m_config["subdet"].NoHistograms.emplace_back("TRT") ; //don't write root histogram on TRT level
130  m_config["subdet"].PrintLog=m_config["subdet"].CalibrateRt;
131  m_config["subdet"].PrintT0Out.emplace_back("Detector") ;
132  m_config["subdet"].PrintRtOut.emplace_back("Detector") ;
133 
134  //individual barrel/ec
135  m_config["1"]=m_config["subdet"];
136  m_config["1"].selections["Detector"].insert(1);
137  m_config["1"].selections["Layer"].insert(CALIB_NONE); // don't go deeper than barrel/ec level
138  m_config["1"].selections["Module"].insert(CALIB_NONE);
139  m_config["1"].selections["Board"].insert(CALIB_NONE);
140  m_config["1"].selections["Chip"].insert(CALIB_NONE);
141  m_config["1"].selections["Straw"].insert(CALIB_NONE);
142  m_config["-1"]=m_config["subdet"];
143  m_config["-1"].selections["Detector"].insert(-1);
144  m_config["-1"].selections["Layer"].insert(CALIB_NONE);
145  m_config["-1"].selections["Module"].insert(CALIB_NONE);
146  m_config["-1"].selections["Board"].insert(CALIB_NONE);
147  m_config["-1"].selections["Chip"].insert(CALIB_NONE);
148  m_config["-1"].selections["Straw"].insert(CALIB_NONE);
149 
150  //common for all layers
151  m_config["layer"].CalibrateRt.emplace_back("TRT") ; //calibrate R-t and T0 on all levels (except Straw level)
152  m_config["layer"].CalibrateRt.emplace_back("Detector") ;
153  m_config["layer"].CalibrateRt.emplace_back("Layer") ;
154  m_config["layer"].CalibrateRt.emplace_back("Module") ;
155  m_config["layer"].CalibrateRt.emplace_back("Board") ;
156  m_config["layer"].CalibrateRt.emplace_back("Chip") ;
157  m_config["layer"].CalibrateT0=m_config["layer"].CalibrateRt;
158  m_config["layer"].NoHistograms.emplace_back("TRT") ; //don't write root histogram on TRT, barrel/ec, straw level
159  m_config["layer"].NoHistograms.emplace_back("Detector") ;
160  m_config["layer"].NoHistograms.emplace_back("Straw") ;
161  m_config["layer"].PrintLog=m_config["layer"].CalibrateRt;
162  m_config["layer"].PrintT0Out.emplace_back("Layer") ;
163  m_config["layer"].PrintT0Out.emplace_back("Module") ;
164  m_config["layer"].PrintT0Out.emplace_back("Board") ;
165  m_config["layer"].PrintT0Out.emplace_back("Chip") ;
166  m_config["layer"].PrintT0Out.emplace_back("Straw") ; //print txt output for each straw
167  m_config["layer"].PrintRtOut.emplace_back("Layer") ;
168  m_config["layer"].PrintRtOut.emplace_back("Module") ;
169  m_config["layer"].PrintRtOut.emplace_back("Board") ;
170  m_config["layer"].PrintRtOut.emplace_back("Chip") ;
171  m_config["layer"].PrintRtOut.emplace_back("Straw") ; //print txt output for each straw
172 
173  //individual layers
174  m_config["1_0"]=m_config["layer"];
175  m_config["1_0"].selections["Detector"].insert(1);
176  m_config["1_0"].selections["Layer"].insert(0);
177  m_config["1_1"]=m_config["layer"];
178  m_config["1_1"].selections["Detector"].insert(1);
179  m_config["1_1"].selections["Layer"].insert(1);
180  m_config["1_2"]=m_config["layer"];
181  m_config["1_2"].selections["Detector"].insert(1);
182  m_config["1_2"].selections["Layer"].insert(2);
183  m_config["-1_0"]=m_config["layer"];
184  m_config["-1_0"].selections["Detector"].insert(-1);
185  m_config["-1_0"].selections["Layer"].insert(0);
186  m_config["-1_1"]=m_config["layer"];
187  m_config["-1_1"].selections["Detector"].insert(-1);
188  m_config["-1_1"].selections["Layer"].insert(1);
189  m_config["-1_2"]=m_config["layer"];
190  m_config["-1_2"].selections["Detector"].insert(-1);
191  m_config["-1_2"].selections["Layer"].insert(2);
192 
193  return StatusCode::SUCCESS;
194 }
195 
196 
197 StatusCode TRTCalibrator::finalize()
198 {
199  return StatusCode::SUCCESS;
200 }
201 
202 
203 std::string TRTCalibrator::SubLev(std::string pat, int lev){
204 
205  std::string sublev;
206 
207  for (int i=0;i<lev;i++){
208  pat.replace(pat.find('_'),1,"");
209  sublev = pat.substr(0,pat.find('_'));
210  pat.replace(0,pat.find('_'),"");
211  }
212  return sublev;
213 }
214 
215 void TRTCalibrator::DumpStrawData(int isid){
216 
217  std::map<std::string,int> strawmap;
218  Identifier ident=(Identifier)(isid);
219 
220 
221 
222  std::string strawkey=std::string(Form("%i %i %i %i %i",
223  (int)m_TRTID->barrel_ec(ident),
224  (int)m_TRTID->layer_or_wheel(ident) ,
225  (int)m_TRTID->phi_module(ident) ,
226  (int)m_TRTID->straw_layer(ident) ,
227  (int)m_TRTID->straw(ident)));
228  strawmap[strawkey]=isid;
229 }
230 
231 
232 bool TRTCalibrator::IsSubLev(const std::string& key, int lev, const std::string& sublev){
233 
234  std::string sl=sublev;
235 
236  if (sl.compare(SubLev(key,lev))==0) return true;
237 
238  while(sl.compare(sl.substr(sl.find(',')+1))!=0){
239  if (SubLev(key,lev).compare(sl.substr(0,sl.find(',')))==0) return true;
240  sl=sl.substr(sl.find(',')+1);
241  if (SubLev(key,lev).compare(sl.substr(0,sl.find(',')))==0) return true;
242  }
243  return false;
244 
245 }
246 
247 int TRTCalibrator::GetSubLevels(const std::string& key, int lev, std::set<int>* levels){
248 
249  std::string sl=SubLev(key,lev);
250 
251  if (sl.compare("-")==0) {
252  levels->insert(-4);
253  return 999;
254  }
255  if (sl.compare("*")==0) {
256  levels->insert(-3);
257  return 999;
258  }
259 
260  if(sl.find(',')!=std::string::npos){
261  int value;
262  while(sl.find(',')!=std::string::npos){
263  sscanf(sl.substr(0,sl.find(',')).data(),"%i",&value);
264  levels->insert(value);
265  sl=sl.substr(sl.find(',')+1);
266  }
267  sscanf(sl.data(),"%i",&value);
268  levels->insert(value);
269  return levels->size();
270  }
271 
272  if(sl.find('t')!=std::string::npos){
273  int min,max;
274  sscanf(sl.substr(0,sl.find('t')).data(),"%i",&min);
275  ATH_MSG_INFO( "min=" << min );
276  sl=sl.substr(sl.find('t')+1);
277  sscanf(sl.data(),"%i",&max);
278  ATH_MSG_INFO("max=" << max);
279  for (int imod=min; imod<=max; imod++){
280  levels->insert(imod);
281  }
282  return levels->size();
283  }
284 
285  int value;
286  sscanf(sl.data(),"%i",&value);
287  levels->insert(value);
288  return levels->size();
289 
290 }
291 
292 void TRTCalibrator::MakeBDKeys(int bec, int low, int mod, int brd, int chp, int strawid){
293  //make map keys for the board-chip view
294 
295  char ckey[100];
296  sprintf(ckey,"_%i_%i_%i_%i_%i_%i",bec,low,mod,brd,chp,strawid); m_Skey=std::string(ckey);
297  sprintf(ckey,"_%i_%i_%i_%i_%i",bec,low,mod,brd,chp); m_Ckey=std::string(ckey);
298  sprintf(ckey,"_%i_%i_%i_%i",bec,low,mod,brd); m_Bkey=std::string(ckey);
299  sprintf(ckey,"_%i_%i_%i",bec,low,mod); m_Mkey=std::string(ckey);
300  sprintf(ckey,"_%i_%i",bec,low); m_Lkey=std::string(ckey);
301  sprintf(ckey,"_%i",bec); m_Dkey=std::string(ckey);
302  sprintf(ckey,"_all"); m_Tkey=std::string(ckey);
303 
304  sprintf(ckey,"_%i_%i",abs(bec),low); m_Lkey_acc=std::string(ckey);
305  sprintf(ckey,"_%i",abs(bec)); m_Dkey_acc=std::string(ckey);
306 
307 }
308 
309 bool TRTCalibrator::calibrate ATLAS_NOT_THREAD_SAFE () {
310 
311  float sid;
312  Identifier ident;
313  const InDetDD::TRT_BaseElement* strawelement;
314  const InDetDD::TRT_BarrelElement* barrelelement;
315  const InDetDD::TRT_EndcapElement* endcapelement;
316  int nTRThist=0, ndethist=0, nlayhist=0, nmodhist=0, nbrdhist=0, nchphist=0, nstwhist=0;
317  int nTRThistAr=0, ndethistAr=0, nlayhistAr=0, nmodhistAr=0, nbrdhistAr=0, nchphistAr=0, nstwhistAr=0;
318  databundle hitdata;
319 
320  std::string infile=m_hittuple;
321  std::string outfile="calibout";
322 
323  std::string detname="Detector";
324  std::string layname="Layer";
325  if (!m_SplitBarrel) {
326  detname="WholeBarrel";
327  layname="WholeBarrelLayer";
328  }
329 
330  //create Calibrator objects for each sub-level
331  Calibrator TRT(0,"TRT",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
332  Calibrator Detector(1,detname.data(),m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
333  Calibrator Layer(2,layname.data(),m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
334  Calibrator Module(3,"Module",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
335  Calibrator Board(4,"Board",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
336  Calibrator Chip(5,"Chip",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
337  Calibrator Straw(6,"Straw",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
338 
339 
340  // make a map of calibrators (so that they can be looped over)
341  std::map<std::string,Calibrator*> calibrators;
342  calibrators["TRT"]=&TRT;
343  calibrators["Detector"]=&Detector;
344  calibrators["Layer"]=&Layer;
345  calibrators["Module"]=&Module;
346  calibrators["Board"]=&Board;
347  calibrators["Chip"]=&Chip;
348  calibrators["Straw"]=&Straw;
349 
350  //use default configurations unless "user" option is used
351  if(m_calsub!="user"){
352  m_doRt=m_config[m_calsub].CalibrateRt;
353  m_doT0=m_config[m_calsub].CalibrateT0;
354  m_doRes=m_config[m_calsub].FitResidual;
355  m_beQuiet=m_config[m_calsub].NoHistograms;
356  m_doOutPrint=m_config[m_calsub].PrintT0Out;
357  m_doRtPrint=m_config[m_calsub].PrintRtOut;
358  m_doLogPrint=m_config[m_calsub].PrintLog;
359  }
360 
361  //configure the calibrators
362  for (const std::string& s : m_doRt) calibrators[s]->dort=true;
363  for (const std::string& s : m_doT0) calibrators[s]->dot0=true;
364  for (const std::string& s : m_doRes) calibrators[s]->dores=true;
365  for (const std::string& s : m_beQuiet) calibrators[s]->bequiet=true;
366  for (const std::string& s : m_doOutPrint) calibrators[s]->printt0=true;
367  for (const std::string& s : m_doRtPrint) calibrators[s]->printrt=true;
368  for (const std::string& s : m_doLogPrint) calibrators[s]->printlog=true;
369  for (const std::string& s : m_useBoardRef) calibrators[s]->usebref=true;
370  for (std::pair<const std::string, Calibrator*>& p : calibrators) {
371  Calibrator* calib = p.second;
372  calib->usep0=m_useP0;
373  calib->floatp3=m_floatP3;
374  calib->useshortstraws=m_DoShortStrawCorrection;
375 
376  if (m_doRt.size()==0) calib->nort=true;
377  if (m_doT0.size()==0) calib->not0=true;
378 
379  std::set<int> selset;
380  GetSubLevels(m_selstring,calib->level+1,&selset); //get the set of selections for level "level+1" from the m_selstring ...
381  for (std::set<int>::iterator imod=selset.begin(); imod!=selset.end(); ++imod){
382  calib->selection = selset;//... and configure the calibrator with them
383  }
384  }
385 
386  // FOR ARGON:
387  Calibrator TRT_Ar (0,"TRT_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
388  Calibrator Detector_Ar(1,"Detector_Ar",m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
389  Calibrator Layer_Ar (2,"Layer_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
390  Calibrator Module_Ar (3,"Module_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
391  Calibrator Board_Ar (4,"Board_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
392  Calibrator Chip_Ar (5,"Chip_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
393  Calibrator Straw_Ar (6,"Straw_Ar" ,m_mint0,m_minrt,m_rtrel,m_rtbinning,m_t0offset);
394 
395  // make a map of calibrators (so that they can be looped over)
396  std::map<std::string,Calibrator*> calibratorsAr;
397  calibratorsAr["TRT"] =&TRT_Ar;
398  calibratorsAr["Detector"] =&Detector_Ar;
399  calibratorsAr["Layer"] =&Layer_Ar;
400  calibratorsAr["Module"] =&Module_Ar;
401  calibratorsAr["Board"] =&Board_Ar;
402  calibratorsAr["Chip"] =&Chip_Ar;
403  calibratorsAr["Straw"] =&Straw_Ar;
404 
405  //configure the calibrators
406  for (const std::string& s : m_doRt) calibratorsAr[s]->dort=true;
407  for (const std::string& s : m_doT0) calibratorsAr[s]->dot0=true;
408  for (const std::string& s : m_doRes) calibratorsAr[s]->dores=true;
409  for (const std::string& s : m_beQuiet) calibratorsAr[s]->bequiet=true;
410  for (const std::string& s : m_doOutPrint) calibratorsAr[s]->printt0=true;
411  for (const std::string& s : m_doRtPrint) calibratorsAr[s]->printrt=true;
412  for (const std::string& s : m_doLogPrint) calibratorsAr[s]->printlog=true;
413  for (const std::string& s : m_useBoardRef) calibratorsAr[s]->usebref=true;
414  for (std::pair<const std::string, Calibrator*>& p : calibratorsAr) {
415  Calibrator* calib = p.second;
416  calib->usep0=m_useP0;
417  calib->floatp3=m_floatP3;
418  calib->useshortstraws=m_DoShortStrawCorrection;
419 
420  if (m_doRt.size()==0) calib->nort=true;
421  if (m_doT0.size()==0) calib->not0=true;
422 
423  std::set<int> selset;
424  GetSubLevels(m_selstring,calib->level+1,&selset); //get the set of selections for level "level+1" from the m_selstring ...
425  for (std::set<int>::iterator imod=selset.begin(); imod!=selset.end(); ++imod){
426  calib->selection = selset;//... and configure the calibrator with them
427  }
428  }
429 
430 
431 
432  m_options = '_' + TRT.GetOptString() + '_' + Detector.GetOptString() + '_' + Layer.GetOptString() + '_' + Module.GetOptString() + '_' + Board.GetOptString() + '_' + Chip.GetOptString() + '_' + Straw.GetOptString();
433  if(m_calsub!="user") m_selstring = '_' + TRT.GetSelString() + '_' + Detector.GetSelString() + '_' + Layer.GetSelString() + '_' + Module.GetSelString() + '_' + Board.GetSelString() + '_' + Chip.GetSelString() + '_' + Straw.GetSelString();
434 
435  //print some info
436  ATH_MSG_INFO(" ");
437  ATH_MSG_INFO( "INPUT FILE : " << infile );
438  ATH_MSG_INFO( "OUTPUT FILE : " << outfile );
439  ATH_MSG_INFO( "SELECTION STRING : " << m_selstring );
440  ATH_MSG_INFO( "OPTION STRING : " << m_options );
441  ATH_MSG_INFO( "RT RELATION : " << m_rtrel );
442  ATH_MSG_INFO( "MIN STATISTICS : RT=" << m_minrt << ", T0=" << m_mint0 );
443  ATH_MSG_INFO( "T0 OFFSET : " << m_t0offset );
444  ATH_MSG_INFO(" ");
445  ATH_MSG_INFO( TRT.PrintInfo() );
446  ATH_MSG_INFO( Detector.PrintInfo() );
447  ATH_MSG_INFO( Layer.PrintInfo() );
448  ATH_MSG_INFO( Module.PrintInfo() );
449  ATH_MSG_INFO( Board.PrintInfo() );
450  ATH_MSG_INFO( Chip.PrintInfo() );
451  ATH_MSG_INFO( Straw.PrintInfo() );
452  ATH_MSG_INFO(" ");
453 
454  if (m_DoArXenonSep){
455  //print some info
456  ATH_MSG_INFO(" ");
457  ATH_MSG_INFO( "Setup for the ARGON Straws: " );
458  ATH_MSG_INFO( "INPUT FILE : " << infile );
459  ATH_MSG_INFO( "OUTPUT FILE : " << outfile );
460  ATH_MSG_INFO( "SELECTION STRING : " << m_selstring );
461  ATH_MSG_INFO( "OPTION STRING : " << m_options );
462  ATH_MSG_INFO( "RT RELATION : " << m_rtrel );
463  ATH_MSG_INFO( "MIN STATISTICS : RT=" << m_minrt << ", T0=" << m_mint0 );
464  ATH_MSG_INFO( "T0 OFFSET : " << m_t0offset );
465  ATH_MSG_INFO( " ");
466  ATH_MSG_INFO( TRT_Ar.PrintInfo() );
467  ATH_MSG_INFO( Detector_Ar.PrintInfo() );
468  ATH_MSG_INFO( Layer_Ar.PrintInfo() );
469  ATH_MSG_INFO( Module_Ar.PrintInfo() );
470  ATH_MSG_INFO( Board_Ar.PrintInfo() );
471  ATH_MSG_INFO( Chip_Ar.PrintInfo() );
472  ATH_MSG_INFO( Straw_Ar.PrintInfo() );
473  ATH_MSG_INFO(" ");
474  }
475  //read the chip reference t0 values from finedelays.txt
476  int rbrd, rchp, rdet, dum;
477  float rt0;
478  std::map<std::string,float> reft0map;
479  std::ifstream t0ref("finedelays.txt",std::ios::in);
480  if(t0ref.is_open()) ATH_MSG_INFO(" Opened finedelays.txt ");
481  for(int iref=0;iref<208;iref++){
482  t0ref >> rbrd >> rchp >> rdet >> dum >> rt0;
483  reft0map[std::string(Form("_%i_%i_%i",rdet,rbrd,rchp))]=rt0;
484  }
485 
486  t0ref.close();
487 
488 
489  //set type of rt-relation
490  bool isdines = (m_rtrel.value()).find("dines")!=std::string::npos;
491  bool isbinned = (m_rtrel.value()).find("binned")!=std::string::npos;
492  int rtint;
493  if (isdines) rtint=2;
494  else if (isbinned) rtint=1;
495  else rtint=0;
496  ATH_MSG_INFO(" Rt relation is: " << rtint);
497 
498  int npop,isid;
499  std::map<std::string,epdata> ephasemap;
500 
501  // open the output histogram file
502  std::unique_ptr<TFile> histfile(TFile::Open("calibout.root","RECREATE"));
503 
504  if(!histfile || histfile->IsZombie()) {
505  ATH_MSG_ERROR ("Failed to open calibout.root ");
506  return false;
507  } else {
508  ATH_MSG_INFO ("Opened calibout.root ");
509  }
510 
511  std::ifstream myFile (infile.data(), std::ios::in | std::ios::binary);
512  ATH_MSG_INFO( " Opened " << infile << " as binary histogram file " );
513  int ihist=0;
514  int ihistAr=0;
515 
516  while(true){
517 
518  //read a binary histogram
519  //dangerous cast reproduces pre-existing c style cast
520  myFile.read (reinterpret_cast<char*>(&npop),sizeof(int)); //number of populated bins
521  if (myFile.eof()) break;
522  int* chist=new int[2*npop+2]; //the histogram
523  //dangerous cast reproduces pre-existing c style cast
524  if (npop>0) myFile.read (reinterpret_cast<char*>(chist+2), sizeof(int)*2*npop);
525  //dangerous cast reproduces pre-existing c style cast
526  myFile.read (reinterpret_cast<char*>(&isid),sizeof(int)); //the straw id
527  sid = (float)isid;
528  if(sid<0) continue;
529 
530  chist[0]=npop;
531  chist[1]=isid;
532  //get the straw address (barrel/ec, layer, module, ...) based on the straw identifier
533  ident=(Identifier)(isid);
534  int chip=0,board=-1;
535  m_neighbourSvc->getChip(ident,chip);
536  if(abs(m_TRTID->barrel_ec(ident))<2){
537  board=m_neighbourSvc->chipToBoardBarrel(chip,m_TRTID->layer_or_wheel(ident));
538  }
539  else if (chip<12) {
540  board=0;
541  }
542  else {
543  chip=chip-20;
544  board=1;
545  }
546  if (m_SplitBarrel) hitdata.det=(int)m_TRTID->barrel_ec(ident);
547  else hitdata.det=abs((int)m_TRTID->barrel_ec(ident));
548  if(hitdata.det!=1 && hitdata.det!=-1 && hitdata.det!=2 && hitdata.det!=-2) {
549  ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " detector decoded to be : " << hitdata.det);
550  continue;
551  }
552  hitdata.lay=(int)m_TRTID->layer_or_wheel(ident);
553  if( (hitdata.det==1 || hitdata.det==-1) && (hitdata.lay<0 || hitdata.lay>2) ) {
554  ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " barrel layer decoded to be : " << hitdata.lay);
555  continue;
556  } else if( (hitdata.det==2 || hitdata.det==-2) && (hitdata.lay<0 || hitdata.lay>13) ) {
557  ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " endcap layer decoded to be : " << hitdata.lay);
558  continue;
559  }
560  hitdata.mod=(int)m_TRTID->phi_module(ident);
561  if( hitdata.mod < 0 || hitdata.mod > 31 ) {
562  ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " phi-module decoded to be : " << hitdata.lay);
563  continue;
564  }
565 
566  hitdata.brd=board;
567  hitdata.chp=chip;
568  hitdata.stl=(int)m_TRTID->straw_layer(ident);
569  if( hitdata.stl < 0 || hitdata.stl > m_TRTID->straw_layer_max(ident) ) {
570  ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " straw-layer decoded to be : " << hitdata.stl);
571  continue;
572  }
573 
574  hitdata.stw=(int)m_TRTID->straw(ident);
575  if( hitdata.stw < 0 || hitdata.stw > m_TRTID->straw_max(ident) ) {
576  ATH_MSG_WARNING( " Invalid Identifier read : " << ident << " straw decoded to be : " << hitdata.stw);
577  continue;
578  }
579 
580  hitdata.sid=isid;
581 
582  //get the old rt parameters based on straw identifier
583  const float* pcal ;
584  std::vector<float> rvalues ;
585  const float defaultpcal[] = {0,0,0,0} ;
586  if (isdines){
587  const TRTCond::DinesRtRelation* rtr = dynamic_cast<const TRTCond::DinesRtRelation*>(m_trtcaldbTool->getRtRelation(ident)) ;
588  pcal = rtr ? rtr->cal() : defaultpcal ;
589  }
590 
591  else {
592  const TRTCond::BasicRtRelation* rtr = dynamic_cast<const TRTCond::BasicRtRelation*>(m_trtcaldbTool->getRtRelation(ident)) ;
593  pcal = rtr ? rtr->cal() : defaultpcal ;
594  }
595  hitdata.rtpar[0]=pcal[0];
596  hitdata.rtpar[1]=pcal[1];
597  hitdata.rtpar[2]=pcal[2];
598  hitdata.rtpar[3]=pcal[3];
599 
600 
601  //build map keys
602  MakeBDKeys(hitdata.det, hitdata.lay, hitdata.mod, hitdata.brd, hitdata.chp, hitdata.sid);
603 
604  //make the level hierachy dictionaries
605  m_trt.t[m_Tkey].d[m_Dkey].l[m_Lkey].m[m_Mkey].b[m_Bkey].c[m_Ckey].s[m_Skey].z=0;
606  m_trt_acc.t[m_Tkey].d[m_Dkey_acc].l[m_Lkey_acc].m[m_Mkey].b[m_Bkey].c[m_Ckey].s[m_Skey].z=0;
607 
608  //populate the hit data structure to be added
609 
610 
611  hitdata.ievt=0;
612  hitdata.tres=0;
613  hitdata.weight=1.0;
614  hitdata.res=0;
615  hitdata.t=0;
616  hitdata.r=0;
617  hitdata.t0=m_trtcaldbTool->getT0(ident);
618 
619  strawelement = m_trtmanager->getElement(ident);
620  if(hitdata.det==1 || hitdata.det==-1) {
621  barrelelement=static_cast<const InDetDD::TRT_BarrelElement*>(strawelement);
622  hitdata.x=(barrelelement->center(ident)).x();
623  hitdata.y=(barrelelement->center(ident)).y();
624  hitdata.z=(barrelelement->center(ident)).z();
625  }else{
626  endcapelement=static_cast<const InDetDD::TRT_EndcapElement*>(strawelement);
627  hitdata.x=(endcapelement->center(ident)).x();
628  hitdata.y=(endcapelement->center(ident)).y();
629  hitdata.z=(endcapelement->center(ident)).z();
630  }
631 
632  //in the short straws corrections, autodetect if it was applied on the previous step
633  if ( m_DoShortStrawCorrection && hitdata.lay==0 && hitdata.stl<9){
634  //If correction was done in ctes in db (readed), undo the correction: This fixes problem on t0 averaging....
635  double t0test1=m_trtcaldbTool->getT0((Identifier)301998432);
636  double t0test2=m_trtcaldbTool->getT0((Identifier)302001504);
637  if (t0test1 != t0test2) hitdata.t0+=0.75; //short straw compensation
638  }
639 
640  //the reference t0 from finedelays
641  hitdata.rt0=reft0map[std::string(Form("_%i_%i_%i",hitdata.det,hitdata.brd,hitdata.chp))];
642 
643  // Prepare for Xe-Ar mixed conditions:
644  int isArgonStraw = 0;
645  if (m_TRTStrawSummaryTool->getStatusHT(ident, Gaudi::Hive::currentContext()) != TRTCond::StrawStatus::Good) {
646  isArgonStraw = 1;
647  }
648 
649 
650  //add histogram to the Calibrators (A and C side separated)
651  if(!m_DoArXenonSep){
652 
653  nTRThist += TRT.AddHit(m_Tkey,hitdata,chist,true);
654  if (Detector.CheckSelection(hitdata.det)) { //only add the histogram if it is in the selection
655  if (m_SplitBarrel) ndethist += Detector.AddHit(m_Dkey,hitdata,chist,true);
656  else ndethist += Detector.AddHit(m_Dkey_acc,hitdata,chist,true);
657  if (Layer.CheckSelection(hitdata.lay)) {
658  if (m_SplitBarrel) nlayhist += Layer.AddHit(m_Lkey,hitdata,chist,true);
659  else { nlayhist += Layer.AddHit(m_Lkey_acc,hitdata,chist,true); continue;}
660  if (Module.CheckSelection(hitdata.mod)) {
661  nmodhist += Module.AddHit(m_Mkey,hitdata,chist,true);
662  if (Board.CheckSelection(hitdata.brd)) {
663  nbrdhist += Board.AddHit(m_Bkey,hitdata,chist,true);
664  if (Chip.CheckSelection(hitdata.chp)) {
665  nchphist += Chip.AddHit(m_Ckey,hitdata,chist,true);
666  if (Straw.CheckSelection(hitdata.stw))
667  nstwhist += Straw.AddHit(m_Skey,hitdata,chist,true);
668  }
669  }
670  }
671  }
672  } // Here it closes
673  ihist++;
674  } else
675  if(isArgonStraw==0){ // Separate Ar and Xe in endcaps. Here Xe
676  nTRThist += TRT.AddHit(m_Tkey,hitdata,chist,true);
677  if (Detector.CheckSelection(hitdata.det)) { //only add the histogram if it is in the selection
678  if (m_SplitBarrel) ndethist += Detector.AddHit(m_Dkey,hitdata,chist,true);
679  else ndethist += Detector.AddHit(m_Dkey_acc,hitdata,chist,true);
680  if (Layer.CheckSelection(hitdata.lay)) {
681  if (m_SplitBarrel) nlayhist += Layer.AddHit(m_Lkey,hitdata,chist,true);
682  else { nlayhist += Layer.AddHit(m_Lkey_acc,hitdata,chist,true); continue;}
683  if (Module.CheckSelection(hitdata.mod)) {
684  nmodhist += Module.AddHit(m_Mkey,hitdata,chist,true);
685  if (Board.CheckSelection(hitdata.brd)) {
686  nbrdhist += Board.AddHit(m_Bkey,hitdata,chist,true);
687  if (Chip.CheckSelection(hitdata.chp)) {
688  nchphist += Chip.AddHit(m_Ckey,hitdata,chist,true);
689  if (Straw.CheckSelection(hitdata.stw)) {
690  nstwhist += Straw.AddHit(m_Skey,hitdata,chist,true);
691  }
692  }
693  }
694  }
695  }
696  } // Here it closes
697  ihist++;
698  } else { // ARGON HITS
699  nTRThistAr += TRT_Ar.AddHit(m_Tkey,hitdata,chist,true);
700  if (Detector_Ar.CheckSelection(hitdata.det)) { //only add the histogram if it is in the selection
701  if (m_SplitBarrel) ndethistAr += Detector_Ar.AddHit(m_Dkey,hitdata,chist,true);
702  else ndethistAr += Detector_Ar.AddHit(m_Dkey_acc,hitdata,chist,true);
703  if (Layer_Ar.CheckSelection(hitdata.lay)) {
704  if (m_SplitBarrel) nlayhistAr += Layer_Ar.AddHit(m_Lkey,hitdata,chist,true);
705  else { nlayhistAr += Layer_Ar.AddHit(m_Lkey_acc,hitdata,chist,true); continue;}
706  if (Module_Ar.CheckSelection(hitdata.mod)) {
707  nmodhistAr += Module_Ar.AddHit(m_Mkey,hitdata,chist,true);
708  if (Board_Ar.CheckSelection(hitdata.brd)) {
709  nbrdhistAr += Board_Ar.AddHit(m_Bkey,hitdata,chist,true);
710  if (Chip_Ar.CheckSelection(hitdata.chp)) {
711  nchphistAr += Chip_Ar.AddHit(m_Ckey,hitdata,chist,true);
712  if (Straw_Ar.CheckSelection(hitdata.stw)) {
713  nstwhistAr += Straw_Ar.AddHit(m_Skey,hitdata,chist,true);
714  }
715  }
716  }
717  }
718  }
719  } // Here it closes
720  ihistAr++;
721  } // END ARGON HITS
722 
723 
724  if ((ihist%10000==9999) | (ihist==m_nevents-1)){
725  ATH_MSG_INFO( Form("%7i HISTOGRAMS READ, UNITS ADDED: %i %i %2i %3i %3i %4i %6i",ihist+1,nTRThist, ndethist, nlayhist, nmodhist, nbrdhist, nchphist, nstwhist) );
726  if(m_DoArXenonSep) ATH_MSG_INFO( Form("%7i Ar HISTOGRAMS READ, UNITS ADDED: %i %i %2i %3i %3i %4i %6i",ihistAr+1,nTRThistAr, ndethistAr, nlayhistAr, nmodhistAr, nbrdhistAr, nchphistAr, nstwhistAr) );
727  }
728 
729  ihist++;
730  delete [] chist;
731  } // Finish Straw Loop
732 
733  ATH_MSG_INFO( Form("%7i HISTOGRAMS READ, UNITS ADDED: %i %i %2i %3i %3i %4i %6i",ihist+1,nTRThist, ndethist, nlayhist, nmodhist, nbrdhist, nchphist, nstwhist) );
734  if(m_DoArXenonSep) ATH_MSG_INFO( Form("%7i Ar HISTOGRAMS READ, UNITS ADDED: %i %i %2i %3i %3i %4i %6i",ihistAr+1,nTRThistAr, ndethistAr, nlayhistAr, nmodhistAr, nbrdhistAr, nchphistAr, nstwhistAr) );
735 
736  TRT.UpdateOldConstants();
737  Detector.UpdateOldConstants();
738  Layer.UpdateOldConstants();
739  Module.UpdateOldConstants();
740  Board.UpdateOldConstants();
741  Chip.UpdateOldConstants();
742 
743  if (ihistAr>0 ){
744  TRT_Ar.UpdateOldConstants();
745  Detector_Ar.UpdateOldConstants();
746  Layer_Ar.UpdateOldConstants();
747  Module_Ar.UpdateOldConstants();
748  Board_Ar.UpdateOldConstants();
749  Chip_Ar.UpdateOldConstants();
750  }
751 
752  myFile.close();
753 
754 
755  ATH_MSG_INFO(" ");
756  ATH_MSG_INFO( TRT.PrintStat() );
757  ATH_MSG_INFO( Detector.PrintStat() );
758  ATH_MSG_INFO( Layer.PrintStat() );
759  ATH_MSG_INFO( Module.PrintStat() );
760  ATH_MSG_INFO( Board.PrintStat() );
761  ATH_MSG_INFO( Chip.PrintStat() );
762  ATH_MSG_INFO( Straw.PrintStat() );
763  ATH_MSG_INFO(" ");
764 
765 
766  if(m_DoArXenonSep){
767  ATH_MSG_INFO(" ");
768  ATH_MSG_INFO( TRT_Ar.PrintStat() );
769  ATH_MSG_INFO( Detector_Ar.PrintStat() );
770  ATH_MSG_INFO( Layer_Ar.PrintStat() );
771  ATH_MSG_INFO( Module_Ar.PrintStat() );
772  ATH_MSG_INFO( Board_Ar.PrintStat() );
773  ATH_MSG_INFO( Chip_Ar.PrintStat() );
774  ATH_MSG_INFO( Straw_Ar.PrintStat() );
775  ATH_MSG_INFO(" ");
776  }
777 
778 
779 
780  //caldata startdata(true,tbins,rbins);
781  caldata startdata(true,55,100);
782  //caldata startdata(true,64,64);
783  startdata.t0=5.0;
784  std::map<std::string,TDirectory*> dirmap;
785  TDirectory* trtdir=gDirectory;
786  TDirectory* detdir=gDirectory;
787  TDirectory* laydir=gDirectory;
788  TDirectory* moddir=gDirectory;
789  TDirectory* brddir=gDirectory;
790  TDirectory* chpdir=gDirectory;
791 
792  std::ofstream rtcalfile(Form("%s_rt.txt",outfile.data()),std::ios::out);
793  std::ofstream binrtcalfile(Form("%s_binrt.txt",outfile.data()),std::ios::out);
794  std::ofstream t0calfile(Form("%s_t0.txt",outfile.data()),std::ios::out);
795 
796 
797  //if (m_SplitBarrel){
798 
799  for (std::pair<const std::string, BDdetector>& pt : m_trt.t) {
800 
801  if (TRT.Skip()) break;
802  if (TRT.HasKey(pt.first)) {
803  trtdir = TRT.Calibrate(histfile.get(),pt.first,SubLev(m_options,1),&startdata);
804  if (TRT.printt0) t0calfile << Form("-3 -1 -1 -1 -1 : %e %e",TRT.data[pt.first].t0,TRT.data[pt.first].t0err) << std::endl;
805  if (TRT.printrt) rtcalfile << Form("-3 -1 -1 -1 -1 : %i %e %e %e %e",rtint,TRT.data[pt.first].rtpar[0],TRT.data[pt.first].rtpar[1],TRT.data[pt.first].rtpar[2],TRT.data[pt.first].rtpar[3]) << std::endl;
806  }
807  for (std::pair<const std::string, BDlayer>& pd : pt.second.d) {
808 
809  if(Detector.Skip()) break;
810  if(Detector.HasKey(pd.first)){
811  detdir = Detector.Calibrate(trtdir,pd.first,SubLev(m_options,2),&TRT.data[pt.first]);
812  }
813  for (std::pair<const std::string, BDmodule>& pl : pd.second.l) {
814 
815  if(Layer.Skip()) break;
816  if(Layer.HasKey(pl.first)){
817  laydir = Layer.Calibrate(detdir,pl.first,SubLev(m_options,3),&Detector.data[pd.first]);
818  if (Layer.printt0) t0calfile << Form("%i %i -1 -1 -1 : %e %e",Layer.data[pl.first].det,Layer.data[pl.first].lay,Layer.data[pl.first].t0,Layer.data[pl.first].t0err) << std::endl;
819  if (Layer.printrt) rtcalfile << Form("%i %i -1 -1 -1 : %i %e %e %e %e",Layer.data[pl.first].det,Layer.data[pl.first].lay,rtint,Layer.data[pl.first].rtpar[0],Layer.data[pl.first].rtpar[1],Layer.data[pl.first].rtpar[2],Layer.data[pl.first].rtpar[3]) << std::endl;
820  if (!m_SplitBarrel) {
821  if (Layer.printt0) t0calfile << Form("%i %i -1 -1 -1 : %e %e",-Layer.data[pl.first].det,Layer.data[pl.first].lay,Layer.data[pl.first].t0,Layer.data[pl.first].t0err) << std::endl;
822  if (Layer.printrt) rtcalfile << Form("%i %i -1 -1 -1 : %i %e %e %e %e",-Layer.data[pl.first].det,Layer.data[pl.first].lay,rtint,Layer.data[pl.first].rtpar[0],Layer.data[pl.first].rtpar[1],Layer.data[pl.first].rtpar[2],Layer.data[pl.first].rtpar[3]) << std::endl;
823  }
824  }
825  for (std::pair<const std::string, BDboard>& pm : pl.second.m) {
826 
827  if(Module.Skip()) break;
828  if(Module.HasKey(pm.first)){
829  moddir = Module.Calibrate(laydir,pm.first,SubLev(m_options,4),&Layer.data[pl.first]);
830  if (Module.printt0) t0calfile << Form("%i %i %i -1 -1 : %e %e",Module.data[pm.first].det,Module.data[pm.first].lay,Module.data[pm.first].mod,Module.data[pm.first].t0,Module.data[pm.first].t0err) << std::endl;
831  if (Module.printrt) rtcalfile << Form("%i %i %i -1 -1 : %i %e %e %e %e",Module.data[pm.first].det,Module.data[pm.first].lay,Module.data[pm.first].mod,rtint,Module.data[pm.first].rtpar[0],Module.data[pm.first].rtpar[1],Module.data[pm.first].rtpar[2],Module.data[pm.first].rtpar[3]) << std::endl;
832  }
833  for (std::pair<const std::string, BDchip>& pb : pm.second.b) {
834 
835  if(Board.Skip()) break;
836  if(Board.HasKey(pb.first)){
837  brddir = Board.Calibrate(moddir,pb.first,SubLev(m_options,5),&Module.data[pm.first]);
838  }
839  for (std::pair<const std::string, BDstraw>& pc : pb.second.c) {
840 
841  if(Chip.Skip()) break;
842  if(Chip.HasKey(pc.first)){
843  chpdir = Chip.Calibrate(brddir,pc.first,SubLev(m_options,6),&Board.data[pb.first]);
844  }
845  for (std::pair<const std::string, BDzero>& ps : pc.second.s) {
846 
847  if(Straw.Skip()) break;
848  if(Straw.HasKey(ps.first)){
849  Straw.Calibrate(chpdir,ps.first,SubLev(m_options,7),&Chip.data[pc.first]);
850  if (Straw.printt0) t0calfile << Form("%i %i %i %i %i : %e %e",Straw.data[ps.first].det,Straw.data[ps.first].lay,Straw.data[ps.first].mod,Straw.data[ps.first].stl,Straw.data[ps.first].stw,Straw.data[ps.first].t0,Straw.data[ps.first].t0err) << std::endl;
851  if (Straw.printrt) rtcalfile << Form("%i %i %i %i %i : %i %e %e %e %e",Straw.data[ps.first].det,Straw.data[ps.first].lay,Straw.data[ps.first].mod,Straw.data[ps.first].stl,Straw.data[ps.first].stw,rtint,Straw.data[ps.first].rtpar[0],Straw.data[ps.first].rtpar[1],Straw.data[ps.first].rtpar[2],Straw.data[ps.first].rtpar[3]) << std::endl;
852 
853 
854  }
855  }
856  }
857  }
858  }
859  }
860  }
861  }
862 
863  std::map<std::string,TDirectory*> dirArmap;
864  TDirectory* trtdirAr=gDirectory;
865  TDirectory* detdirAr=gDirectory;
866  TDirectory* laydirAr=gDirectory;
867  TDirectory* moddirAr=gDirectory;
868  TDirectory* brddirAr=gDirectory;
869  TDirectory* chpdirAr=gDirectory;
870 
871  if (m_DoArXenonSep){
872 
873  for (std::pair<const std::string, BDdetector>& pt : m_trt.t) {
874 
875  if (TRT_Ar.Skip()) break;
876  if (TRT_Ar.HasKey(pt.first)) {
877  trtdirAr = TRT_Ar.Calibrate(histfile.get(),pt.first,SubLev(m_options,1),&startdata);
878  }
879  for (std::pair<const std::string, BDlayer>& pd : pt.second.d) {
880 
881  if(Detector_Ar.Skip()) break;
882  if(Detector_Ar.HasKey(pd.first)){
883  detdirAr = Detector_Ar.Calibrate(trtdirAr,pd.first,SubLev(m_options,2),&TRT_Ar.data[pt.first]);
884  }
885  for (std::pair<const std::string, BDmodule>& pl : pd.second.l) {
886 
887  if(Layer_Ar.Skip()) break;
888  if(Layer_Ar.HasKey(pl.first)){
889  laydirAr = Layer_Ar.Calibrate(detdirAr,pl.first,SubLev(m_options,3),&Detector_Ar.data[pd.first]);
890 
891  }
892  for (std::pair<const std::string, BDboard>& pm : pl.second.m) {
893 
894  if(Module_Ar.Skip()) break;
895  if(Module_Ar.HasKey(pm.first)){
896  moddirAr = Module_Ar.Calibrate(laydirAr,pm.first,SubLev(m_options,4),&Layer_Ar.data[pl.first]);
897  if (Module_Ar.printt0) t0calfile << Form("%i %i %i -1 -1 : %e %e",Module_Ar.data[pm.first].det,Module_Ar.data[pm.first].lay,Module_Ar.data[pm.first].mod,Module_Ar.data[pm.first].t0,Module_Ar.data[pm.first].t0err) << std::endl;
898  if (Layer_Ar.printrt) rtcalfile << Form("%i %i %i -1 -1 : %i %e %e %e %e",Module_Ar.data[pm.first].det,Module_Ar.data[pm.first].lay,Module_Ar.data[pm.first].mod,rtint,Module_Ar.data[pm.first].rtpar[0],Module_Ar.data[pm.first].rtpar[1],Module_Ar.data[pm.first].rtpar[2],Module_Ar.data[pm.first].rtpar[3]) << std::endl;
899  }
900  for (std::pair<const std::string, BDchip>& pb : pm.second.b) {
901  if(Board_Ar.Skip()) break;
902  if(Board_Ar.HasKey(pb.first)){
903  brddirAr = Board_Ar.Calibrate(moddirAr,pb.first,SubLev(m_options,5),&Module_Ar.data[pm.first]);
904  }
905  for (std::pair<const std::string, BDstraw>& pc : pb.second.c) {
906 
907  if(Chip_Ar.Skip()) break;
908  if(Chip_Ar.HasKey(pc.first)){
909  chpdirAr = Chip_Ar.Calibrate(brddirAr,pc.first,SubLev(m_options,6),&Board_Ar.data[pb.first]);
910  }
911  for (std::pair<const std::string, BDzero>& ps : pc.second.s) {
912 
913  if(Straw_Ar.Skip()) break;
914  if(Straw_Ar.HasKey(ps.first)){
915  Straw_Ar.Calibrate(chpdirAr,ps.first,SubLev(m_options,7),&Chip_Ar.data[pc.first]);
916  if (Straw_Ar.printt0) t0calfile << Form("%i %i %i %i %i : %e %e",Straw_Ar.data[ps.first].det,Straw_Ar.data[ps.first].lay,Straw_Ar.data[ps.first].mod,Straw_Ar.data[ps.first].stl,Straw_Ar.data[ps.first].stw,Straw_Ar.data[ps.first].t0,Straw_Ar.data[ps.first].t0err) << std::endl;
917  if (Straw_Ar.printrt) rtcalfile << Form("%i %i %i %i %i : %i %e %e %e %e",Straw_Ar.data[ps.first].det,Straw_Ar.data[ps.first].lay,Straw_Ar.data[ps.first].mod,Straw_Ar.data[ps.first].stl,Straw_Ar.data[ps.first].stw,rtint,Straw_Ar.data[ps.first].rtpar[0],Straw_Ar.data[ps.first].rtpar[1],Straw_Ar.data[ps.first].rtpar[2],Straw_Ar.data[ps.first].rtpar[3]) << std::endl;
918  }
919  }
920  }
921  }
922  }
923  }
924  }
925  }
926 
927 }// Close for ARGON
928 
929 
930 
931 
932  //make ntuples for each calibrator
933  if (!TRT.bequiet) TRT.WriteStat(histfile.get());
934  if (!Detector.bequiet) Detector.WriteStat(histfile.get());
935  if (!Layer.bequiet) Layer.WriteStat(histfile.get());
936  if (!Module.bequiet) Module.WriteStat(histfile.get());
937  if (!Board.bequiet) Board.WriteStat(histfile.get());
938  Chip.WriteStat(histfile.get());
939  Straw.WriteStat(histfile.get());
940 
941  TRT.DumpConstants();
942  Detector.DumpConstants();
943  Layer.DumpConstants();
944  Module.DumpConstants();
945  Board.DumpConstants();
946  Chip.DumpConstants();
947  Straw.DumpConstants();
948 
949  rtcalfile.close();
950  binrtcalfile.close();
951  t0calfile.close();
952 
953  //AR
954  //make ntuples for each calibrator
955  if(m_DoArXenonSep){
956  if (!TRT_Ar.bequiet) TRT_Ar.WriteStat(histfile.get());
957  if (!Detector_Ar.bequiet) Detector_Ar.WriteStat(histfile.get());
958  if (!Layer_Ar.bequiet) Layer_Ar.WriteStat(histfile.get());
959  if (!Module_Ar.bequiet) Module_Ar.WriteStat(histfile.get());
960  if (!Board_Ar.bequiet) Board_Ar.WriteStat(histfile.get());
961  Chip_Ar.WriteStat(histfile.get());
962  Straw_Ar.WriteStat(histfile.get());
963 
964  TRT_Ar.DumpConstants();
965  Detector_Ar.DumpConstants();
966  Layer_Ar.DumpConstants();
967  Module_Ar.DumpConstants();
968  Board_Ar.DumpConstants();
969  Chip_Ar.DumpConstants();
970  Straw_Ar.DumpConstants();
971  }
972 
973  histfile->ls();
974  ATH_MSG_INFO( "writing out calibout.root");
975  histfile->Write();
976 
977  return true;
978 
979 }
980 
981 bool TRTCalibrator::fill(const Trk::Track* aTrack, TRT::TrackInfo* output) {
982  if (aTrack==nullptr) {};
983  if (output==nullptr) {};
984  return true;
985 }
986 
987 
python.PyKernel.retrieve
def retrieve(aClass, aKey=None)
Definition: PyKernel.py:110
xAOD::iterator
JetConstituentVector::iterator iterator
Definition: JetConstituentVector.cxx:68
TRTCond::DinesRtRelation
Definition: DinesRtRelation.h:26
Calibrator::printt0
bool printt0
if true a t0 entry in the calibration output file is prined for each sub-module in this sub-level
Definition: Calibrator.h:349
TRTCalibrator.h
databundle::det
int det
detector (barrel -1 or 1, or end-cap -2 or 2)
Definition: Calibrator.h:79
TRT
@ TRT
Definition: RegSelEnums.h:26
Amg::compare
std::pair< int, int > compare(const AmgSymMatrix(N) &m1, const AmgSymMatrix(N) &m2, double precision=1e-9, bool relative=false)
compare two matrices, returns the indices of the first element that fails the condition,...
Definition: EventPrimitivesHelpers.h:109
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
DinesRtRelation.h
InDetDD::TRT_BarrelElement
Definition: TRT_BarrelElement.h:43
TrackParameters.h
run.infile
string infile
Definition: run.py:13
ComTime.h
keylayer_zslicemap.pb
pb
Definition: keylayer_zslicemap.py:188
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
find
std::string find(const std::string &s)
return a remapped string
Definition: hcg.cxx:135
TRT_DetectorManager.h
Trk::Track
The ATLAS Track class.
Definition: Tracking/TrkEvent/TrkTrack/TrkTrack/Track.h:73
ParticleGun_SamplingFraction.bec
int bec
Definition: ParticleGun_SamplingFraction.py:89
TRTCalibrator::finalize
virtual StatusCode finalize() override
Definition: TRTCalibrator.cxx:197
max
constexpr double max()
Definition: ap_fixedTest.cxx:33
RtRelation.h
abstract base class for rt-relations
Calibrator::UpdateOldConstants
int UpdateOldConstants()
...
Definition: Calibrator.cxx:529
Calibrator
A structure to contain data associated with the calibration of a certain sub-module.
Definition: Calibrator.h:178
min
constexpr double min()
Definition: ap_fixedTest.cxx:26
Calibrator::HasKey
bool HasKey(const std::string &) const
...
Definition: Calibrator.cxx:452
TRTCalibrator::m_TRTID
const TRT_ID * m_TRTID
Definition: TRTCalibrator.h:155
TRTCalibrator::m_Dkey
std::string m_Dkey
Definition: TRTCalibrator.h:162
TRTCalibrator::m_Dkey_acc
std::string m_Dkey_acc
Definition: TRTCalibrator.h:162
TRT_ID.h
This is an Identifier helper class for the TRT subdetector. This class is a factory for creating comp...
test_pyathena.pt
pt
Definition: test_pyathena.py:11
python.AthDsoLogger.out
out
Definition: AthDsoLogger.py:70
caldata
A structure to contain data associated with the calibration of a certain sub-module.
Definition: Calibrator.h:104
databundle::lay
int lay
layer
Definition: Calibrator.h:80
InDetDD::TRT_EndcapElement
Definition: TRT_EndcapElement.h:43
Calibrator::CheckSelection
bool CheckSelection(int)
Checks if a given sub-module is selected for calibration.
Definition: Calibrator.cxx:457
Calibrator::printrt
bool printrt
if true an rt entry in the calibration output file is prined for each sub-module in this sub-level
Definition: Calibrator.h:350
athena.value
value
Definition: athena.py:124
makeDTCalibBlob_pickPhase.pd
pd
Definition: makeDTCalibBlob_pickPhase.py:342
TRTCalibrator::m_DetID
const AtlasDetectorID * m_DetID
Definition: TRTCalibrator.h:154
BasicRtRelation.h
databundle::brd
int brd
board
Definition: Calibrator.h:82
TRTCalibrator::m_config
std::map< std::string, conf > m_config
Definition: TRTCalibrator.h:187
x
#define x
TRT::TrackInfo
Definition: InnerDetector/InDetCalibEvent/TRT_CalibData/TRT_CalibData/TrackInfo.h:83
PrepRawData.h
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
python.RingerConstants.Layer
Layer
Definition: RingerConstants.py:42
AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore
const ServiceHandle< StoreGateSvc > & detStore() const
The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
Definition: AthCommonDataStore.h:95
Track.h
Calibrator::Skip
bool Skip()
...
Definition: Calibrator.cxx:462
databundle::stw
int stw
straw number (within the strawlayer)
Definition: Calibrator.h:85
Calibrator::PrintInfo
std::string PrintInfo()
Definition: Calibrator.cxx:467
databundle::mod
int mod
phi module
Definition: Calibrator.h:81
AtlasDetectorID.h
This class provides an interface to generate or decode an identifier for the upper levels of the dete...
TRTCond::DinesRtRelation::cal
const float * cal() const
return to calibration constants
Definition: DinesRtRelation.h:76
TRTCalibrator::m_ntrtmanager
Gaudi::Property< std::string > m_ntrtmanager
Definition: TRTCalibrator.h:143
TRT_ID::straw
int straw(const Identifier &id) const
Definition: TRT_ID.h:896
python.utils.AtlRunQueryDQUtils.p
p
Definition: AtlRunQueryDQUtils.py:209
TrackInfo.h
ATH_MSG_ERROR
#define ATH_MSG_ERROR(x)
Definition: AthMsgStreamMacros.h:33
maskDeadModules.mod
mod
Definition: maskDeadModules.py:36
TRT
Definition: HitInfo.h:33
ATLAS_NOT_THREAD_SAFE
bool TRTCalibrator::calibrate ATLAS_NOT_THREAD_SAFE()
Install fatal handler with default options.
Definition: TRTCalibrator.cxx:309
Calibrator::GetOptString
std::string GetOptString() const
Creates a string summarizing what is being done at this sub-level.
Definition: Calibrator.cxx:485
ParseInputs.gDirectory
gDirectory
Definition: Final2012/ParseInputs.py:133
lumiFormat.i
int i
Definition: lumiFormat.py:85
z
#define z
python.SystemOfUnits.ps
float ps
Definition: SystemOfUnits.py:150
TRTCond::BasicRtRelation
Definition: BasicRtRelation.h:26
databundle::x
float x
straw x position
Definition: Calibrator.h:96
caldata::t0
float t0
the new t0
Definition: Calibrator.h:126
TRTCalibrator::MakeBDKeys
void MakeBDKeys(int, int, int, int, int, int)
makes the key strings for a sub-module
Definition: TRTCalibrator.cxx:292
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
BinnedRtRelation.h
databundle::res
float res
space residual
Definition: Calibrator.h:90
Calibrator::DumpConstants
void DumpConstants()
Definition: Calibrator.cxx:1225
TRTCalibrator::m_TRTStrawSummaryTool
ToolHandle< ITRT_StrawStatusSummaryTool > m_TRTStrawSummaryTool
Definition: TRTCalibrator.h:116
TRTCalibrator::m_Skey
std::string m_Skey
Definition: TRTCalibrator.h:162
test_pyathena.parent
parent
Definition: test_pyathena.py:15
jet::CompCategory::Detector
@ Detector
Definition: UncertaintyEnum.h:19
databundle::ievt
int ievt
event number
Definition: Calibrator.h:87
TRTCalibrator::SubLev
std::string SubLev(std::string, int)
return sub-level i from a selection string with the form _X_X_X_X_X_X_X
Definition: TRTCalibrator.cxx:203
TRT_DriftCircleOnTrack.h
TRT_ID::barrel_ec
int barrel_ec(const Identifier &id) const
Values of different levels (failure returns 0)
Definition: TRT_ID.h:860
TRT_ID::straw_layer
int straw_layer(const Identifier &id) const
Definition: TRT_ID.h:887
TRT_ID::layer_or_wheel
int layer_or_wheel(const Identifier &id) const
Definition: TRT_ID.h:878
TRT_DriftCircle.h
databundle::weight
float weight
histogram weight (currently set to 1)
Definition: Calibrator.h:89
TRTCond::StrawStatus::Good
@ Good
Definition: StrawStatus.h:18
merge.output
output
Definition: merge.py:16
dso-stats.pat
pat
Definition: dso-stats.py:37
Calibrator::WriteStat
void WriteStat(TDirectory *)
Creates an ntuple with entries containing data associated with the sub-modules in a sub level.
Definition: Calibrator.cxx:1180
EventPrimitives.h
xAOD::Straw
@ Straw
Definition: TrackingPrimitives.h:558
databundle::rtpar
std::array< float, 4 > rtpar
rt-parameters used in the reconstruction
Definition: Calibrator.h:95
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:240
RIO_OnTrack.h
TRTCalibrator::CALIB_NONE
@ CALIB_NONE
Definition: TRTCalibrator.h:189
LArConditions2Ntuple.ckey
string ckey
Definition: LArConditions2Ntuple.py:175
Calibrator::PrintStat
std::string PrintStat()
Prints some sub-level statistics.
Definition: Calibrator.cxx:478
databundle
A structure to contain hit data.
Definition: Calibrator.h:77
TRT::Hit::ident
@ ident
Definition: HitInfo.h:77
TRT_ID::phi_module
int phi_module(const Identifier &id) const
Definition: TRT_ID.h:869
LocalParameters.h
TRTCalibrator::m_trtmanager
const InDetDD::TRT_DetectorManager * m_trtmanager
Definition: TRTCalibrator.h:156
databundle::chp
int chp
chip
Definition: Calibrator.h:83
PlotSFuncertainty.calib
calib
Definition: PlotSFuncertainty.py:110
TRTCalibrator::m_neighbourSvc
ServiceHandle< ITRT_StrawNeighbourSvc > m_neighbourSvc
Definition: TRTCalibrator.h:117
databundle::z
float z
straw z position
Definition: Calibrator.h:98
Calibrator::AddHit
int AddHit(const std::string &, const databundle &, int *, bool)
Adds hits to a sub-module either in the form of a single straw hit or a histogram containing all the ...
Definition: Calibrator.cxx:962
Calibrator.h
Calibrator::GetSelString
std::string GetSelString()
Creates a string summarizing which modules are being calibrated at this sub-level.
Definition: Calibrator.cxx:506
databundle::rt0
float rt0
reference t0 (offset from board mean)
Definition: Calibrator.h:94
y
#define y
TRTCalibrator::m_Mkey
std::string m_Mkey
Definition: TRTCalibrator.h:162
python.CaloAddPedShiftConfig.int
int
Definition: CaloAddPedShiftConfig.py:45
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
TRTCalibrator::m_Bkey
std::string m_Bkey
Definition: TRTCalibrator.h:162
TRTCalibrator::m_Tkey
std::string m_Tkey
Definition: TRTCalibrator.h:162
TRTCalibrator::m_Ckey
std::string m_Ckey
Definition: TRTCalibrator.h:162
InDetDD::TRT_BaseElement::center
virtual const Amg::Vector3D & center() const override final
Element Surface: center of a straw layer.
databundle::stl
int stl
straw-layer
Definition: Calibrator.h:84
databundle::r
float r
drift radius from r(t) relation
Definition: Calibrator.h:92
python.SystemOfUnits.s
float s
Definition: SystemOfUnits.py:147
databundle::sid
int sid
straw ID
Definition: Calibrator.h:86
TRTCalibrator::initialize
virtual StatusCode initialize() override
Pre-define standard calibration configuration.
Definition: TRTCalibrator.cxx:54
TRTCalibrator::DumpStrawData
void DumpStrawData(int)
returns true if the ...
Definition: TRTCalibrator.cxx:215
TRTCalibrator::TRTCalibrator
TRTCalibrator(const std::string &type, const std::string &name, const IInterface *parent)
Constructor.
Definition: TRTCalibrator.cxx:44
TRTCond::BasicRtRelation::cal
const float * cal() const
return to calibration constants
Definition: BasicRtRelation.h:76
TRTCalibrator::fill
virtual bool fill(const Trk::Track *aTrack, TRT::TrackInfo *output) override
Definition: TRTCalibrator.cxx:981
TRTCalibrator::GetSubLevels
int GetSubLevels(const std::string &, int, std::set< int > *)
fills a set of integers with the indexes af the sub-modules that are conatined in a selection
Definition: TRTCalibrator.cxx:247
TRTCalibrator::m_trtcaldbTool
ToolHandle< ITRT_CalDbTool > m_trtcaldbTool
Definition: TRTCalibrator.h:115
AthAlgTool
Definition: AthAlgTool.h:26
Calibrator::bequiet
bool bequiet
if true no histograms are written to the root file
Definition: Calibrator.h:347
databundle::tres
float tres
time residual
Definition: Calibrator.h:88
TRTCalibrator::m_Lkey
std::string m_Lkey
Definition: TRTCalibrator.h:162
databundle::t0
float t0
t0 value used in reconstruction
Definition: Calibrator.h:93
TRTCalibrator::IsSubLev
bool IsSubLev(const std::string &, int, const std::string &)
returns true if the ...
Definition: TRTCalibrator.cxx:232
PrepareReferenceFile.outfile
outfile
Definition: PrepareReferenceFile.py:42
module_driven_slicing.histfile
histfile
Definition: module_driven_slicing.py:570
databundle::y
float y
straw y position
Definition: Calibrator.h:97
python.SystemOfUnits.pc
float pc
Definition: SystemOfUnits.py:114
TRTCalibrator::m_Lkey_acc
std::string m_Lkey_acc
Definition: TRTCalibrator.h:162
Calibrator::data
std::map< std::string, caldata > data
A map between the sub-module identifier string and the calibration data structure (caldata)
Definition: Calibrator.h:339
extractSporadic.myFile
myFile
Definition: extractSporadic.py:86
InDetDD::TRT_BaseElement
Definition: TRT_BaseElement.h:52
databundle::t
float t
raw time
Definition: Calibrator.h:91
python.LArMinBiasAlgConfig.float
float
Definition: LArMinBiasAlgConfig.py:65
mapkey::key
key
Definition: TElectronEfficiencyCorrectionTool.cxx:37
Identifier
Definition: IdentifierFieldParser.cxx:14
Generated on Mon Sep 1 2025 21:23:08 for ATLAS Offline Software by   doxygen 1.8.18