ATLAS Offline Software
EFTrackingSmearingAlg.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
3 */
4 
5 
6 #include "EFTrackingSmearingAlg.h"
7 #include "xAODTracking/VertexContainer.h"
8 #include "AthContainers/ConstAccessor.h"
9 
10 #include "TH1.h"
11 
12 #ifndef XAOD_ANALYSIS
13 #include "EventPrimitives/EventPrimitivesHelpers.h"
14 #include "TrkEventPrimitives/ParamDefs.h"
15 #else
16 
17 namespace Amg
18 {
19  // expression template to evaluate the required size of the compressed matrix at compile time
20  template<int N>
21  struct CalculateCompressedSize {
22  static const int value = N + CalculateCompressedSize<N - 1>::value;
23  };
24 
25  template<>
26  struct CalculateCompressedSize<1> {
27  static const int value = 1;
28  };
29 
30  template<int N>
31  inline void compress(const AmgSymMatrix(N)& covMatrix, std::vector<float>& vec ) {
32  vec.reserve(CalculateCompressedSize<N>::value);
33  for (unsigned int i = 0; i < N ; ++i)
34  for (unsigned int j = 0; j <= i; ++j)
35  vec.push_back(covMatrix(i,j));
36  }
37  inline void compress(const MatrixX& covMatrix, std::vector<float>& vec) {
38  int rows = covMatrix.rows();
39  for (int i = 0; i < rows; ++i) {
40  for (int j = 0; j <= i; ++j) {
41  vec.push_back(covMatrix(i, j));
42  }
43  }
44  }
45 }
46 
47 namespace Trk
48 {
49  enum ParamDefs {
50  // Enums for LocalParameters - LocalPosition/
51  loc1=0, loc2=1,
52  locX=0, locY=1 ,
53  locRPhi=0, locPhiR=0, locZ=1,
54  locR=0, locPhi=1,
55  iPhi=0, iEta=1,
56  distPhi=0, distEta=1,
57  driftRadius=0,
58  // Enums for const Amg::Vector3D & GlobalMomentum /
59  x=0, y=1, z=2,
60  px=0, py=1, pz=2,
61  // Enums for Perigee //
62  d0=0, z0=1, phi0=2, theta=3, qOverP=4,
63  /* Enums for TrackState on Surfaces
64  The first two enums in the TrackParameters refer to the local Frame, i.e.
65  - LocalCartesian for AtanArbitraryPlane
66  - LocalCylindrical for AtaCylinder (includes line)
67  - LocalPolar for AtaDisc
68  The other three enums are standard \f$(\phi, \eta, \frac{q}{p_{T}})\f$
69  */
70  phi=2,
72  u=0, v=1,
73 
75  trkMass=5
76 
77  };
78 
79 }
80 #endif
81 
82 
83 // from Tracking/TrkEventCnv/TrkEventTPCnv/TrkEventTPCnv/helpers/EigenHelpers.h
84 namespace EigenHelpers
85 {
86  template <class VECTOR, class COVARIANCE>
87  inline static void eigenMatrixToVector(VECTOR& vec, COVARIANCE& cov, const char* ) {
88  Amg::compress(cov, vec);
89  }
90 }
91 
92 
93 
94 EFTrackingSmearingAlg::EFTrackingSmearingAlg( const std::string& name, ISvcLocator* pSvcLocator )
95 : AthHistogramAlgorithm( name, pSvcLocator ){}
96 
97 
98 EFTrackingSmearingAlg::~EFTrackingSmearingAlg() {}
99 
100 
101 StatusCode EFTrackingSmearingAlg::initialize() {
102  ATH_MSG_INFO ("Initializing " << name() << "...");
103 
104  ATH_CHECK( m_inputTrackParticleKey.initialize() );
105  ATH_CHECK( m_outputTrackParticleKey.initialize() );
106 
107  ATH_CHECK( m_inputTruthParticleKey.initialize() );
108  ATH_CHECK( m_outputTruthParticleKey.initialize() );
109 
110  // Decoration keys
111  m_d0DecoratorKey = m_outputTruthParticleKey.key() + ".d0";
112  m_z0DecoratorKey = m_outputTruthParticleKey.key() + ".z0";
113  m_ptDecoratorKey = m_outputTruthParticleKey.key() + ".pt";
114 
115  ATH_CHECK(m_d0DecoratorKey.initialize());
116  ATH_CHECK(m_z0DecoratorKey.initialize());
117  ATH_CHECK(m_ptDecoratorKey.initialize());
118 
119  ATH_MSG_INFO("########## EFTrackingSmearingAlg Configurations are ########## ");
120  ATH_MSG_INFO("------- InputTrackParticleKey: "<<m_inputTrackParticleKey.key());
121  ATH_MSG_INFO("------- OutputTrackParticleKey: "<<m_outputTrackParticleKey.key());
122  ATH_MSG_INFO("------- InputTruthParticleKey: " <<m_inputTruthParticleKey.key());
123  ATH_MSG_INFO("------- OutputTruthParticleKey: "<<m_outputTruthParticleKey.key());
124  ATH_MSG_INFO("------- inputTracksPtCut [GeV]: "<<m_inputTracksPtCut);
125  ATH_MSG_INFO("------- outputTracksPtCut [GeV]: "<<m_outputTracksPtCut);
126  ATH_MSG_INFO("------- SmearingSigma: "<<m_SigmaScaleFactor);
127  ATH_MSG_INFO("------- trackEfficiency: "<<m_smearedTrackEfficiency);
128  ATH_MSG_INFO("------- parameterizedEfficiency: "<<m_parameterizedTrackEfficiency);
129  ATH_MSG_INFO("------- parameterizedEfficiency LRT: "<<m_parameterizedTrackEfficiency_LRT);
130  ATH_MSG_INFO("------- parameterizedEfficiency LRT high d0 cut: "<<m_smearedTrackEfficiency_d0high_LRT );
131  ATH_MSG_INFO("------- parameterizedEfficiency LRT low d0 cut: "<<m_smearedTrackEfficiency_d0low_LRT );
132  ATH_MSG_INFO("------- UseResolutionPtCutOff: "<<m_UseResolutionPtCutOff);
133  ATH_MSG_INFO("------- SetResolutionPtCutOff: "<<m_SetResolutionPtCutOff);
134  ATH_MSG_INFO("------- EnableMonitoring:" <<m_enableMonitoring);
135  ATH_MSG_INFO("------- SmearTruthParticle:"<< m_smearTruthParticle);
136 
137  ATH_MSG_INFO("------- IncludeDuplicatesAndFakes: "<<m_EnableFakes);
138  ATH_MSG_INFO("------- RandomSeed: " <<m_RandomSeed);
139  ATH_MSG_INFO("------- UseCoinToss: " <<m_UseCoinToss);
140  ATH_MSG_INFO("------- FakeKillerEnable: "<<m_FakeKillerEnable);
141  ATH_MSG_INFO("------- IncludeFakesInResolutionCalculation: "<<m_IncludeFakesInResolutionCalculation);
142  ATH_MSG_INFO("########## EFTrackingSmearingAlg Configurations: That's it. ########## ");
143 
144  std::string smearerName;
145  if (m_smearTruthParticle){
146  ATH_MSG_INFO("Will output new truth container with name "<<m_outputTruthParticleKey.key());
147  smearerName = m_outputTruthParticleKey.key()+"_smearer";
148  }
149  else {
150  ATH_MSG_INFO("Will output new track container with name "<<m_outputTrackParticleKey.key());
151  smearerName = m_outputTrackParticleKey.key()+"_smearer";
152  }
153 
154  // configure the Smearer
155  m_mySmearer = std::make_unique<FakeTrackSmearer>(smearerName, m_RandomSeed, msgLvl (MSG::DEBUG));
156  m_mySmearer->SetInputTracksPtCut(m_inputTracksPtCut);
157  m_mySmearer->SetOutputTracksPtCut(m_outputTracksPtCut);
158  m_mySmearer->SetTrackingEfficiency(m_smearedTrackEfficiency);
159  m_mySmearer->SetParameterizedEfficiency(m_parameterizedTrackEfficiency);
160  m_mySmearer->SetParameterizedEfficiency_LRT(m_parameterizedTrackEfficiency_LRT);
161  m_mySmearer->SetParameterizedEfficiency_highd0_LRT(m_smearedTrackEfficiency_d0high_LRT);
162  m_mySmearer->SetParameterizedEfficiency_lowd0_LRT(m_smearedTrackEfficiency_d0low_LRT);
163 
164  m_mySmearer->SetSigmaScaleFactor(m_SigmaScaleFactor.value());
165  m_mySmearer->UseResolutionPtCutOff(m_UseResolutionPtCutOff.value());
166  m_mySmearer->SetResolutionPtCutOff(m_SetResolutionPtCutOff.value());
167 
168  m_mySmearer->EnableFakes(m_EnableFakes.value());
169  m_mySmearer->UseCoinToss(m_UseCoinToss.value());
170  m_mySmearer->FakeKillerEnable(m_FakeKillerEnable.value());
171  m_mySmearer->IncludeFakesInResolutionCalculation(m_IncludeFakesInResolutionCalculation.value());
172 
173  if (m_enableMonitoring) {
174  // store the smearing functions
175  TF1 *d0res_eta = m_mySmearer->d0res_eta;
176  TF1 *z0res_eta = m_mySmearer->z0res_eta;
177  TF1 *curvres_eta = m_mySmearer->curvres_eta;
178  TF1 *d0res_pt = m_mySmearer->d0res_pt;
179  TF1 *z0res_pt = m_mySmearer->z0res_pt;
180  TF1 *curvres_pt = m_mySmearer->curvres_pt;
181  TF1 *effLRT_d0 = m_mySmearer->effLRT_d0;
182 
183  CHECK(book(new TH1F("d0res_function_vs_eta","#eta of track (p_{T}=10GeV);#eta",100, 0.0,4.0)));
184  CHECK(book(new TH1F("z0res_function_vs_eta","#eta of track (p_{T}=10GeV);#eta",100, 0.0,4.0)));
185  CHECK(book(new TH1F("curvres_function_vs_eta","#eta of track (p_{T}=10GeV);#eta",100, 0.0,4.0)));
186  CHECK(book(new TH1F("d0res_function_vs_pt","p_{T} of track (#eta=1);p_{T} [GeV]",100, 1.0,200.0)));
187  CHECK(book(new TH1F("z0res_function_vs_pt","p_{T} of track (#eta=1);p_{T} [GeV]",100, 1.0,200.0)));
188  CHECK(book(new TH1F("curvres_function_vs_pt","p_{T} of track (#eta=1);p_{T} [GeV]",100, 1.0,200.0)));
189  CHECK(book(new TH1F("effLRT_function_vs_d0","d_{0} of track;d_{0} [mm]",100, 0.001,600.0)));
190  hist("d0res_function_vs_eta")->Add(d0res_eta);
191  hist("z0res_function_vs_eta")->Add(z0res_eta);
192  hist("curvres_function_vs_eta")->Add(curvres_eta);
193  hist("d0res_function_vs_pt")->Add(d0res_pt);
194  hist("z0res_function_vs_pt")->Add(z0res_pt);
195  hist("curvres_function_vs_pt")->Add(curvres_pt);
196  hist("effLRT_function_vs_d0")->Add(effLRT_d0);
197  // book historgams
198  CHECK(book_histograms());
199  }
200  return StatusCode::SUCCESS;
201 }
202 
203 
204 StatusCode EFTrackingSmearingAlg::smearTruthParticles(const EventContext& ctx) {
205 
206  SG::ReadHandle<xAOD::TruthParticleContainer> inputTruth_handle( m_inputTruthParticleKey, ctx );
207  const xAOD::TruthParticleContainer* inputTruth = inputTruth_handle.cptr();
208  if (not inputTruth) {
209  ATH_MSG_FATAL("Unable to retrieve input truth particle");
210  return StatusCode::FAILURE;
211  }
212  SG::WriteHandle<xAOD::TruthParticleContainer> outputTruth_handle( m_outputTruthParticleKey, ctx );
213  ATH_CHECK( outputTruth_handle.record( std::make_unique<xAOD::TruthParticleContainer>(), std::make_unique<xAOD::TruthParticleAuxContainer>() ) );
214  auto outputTruth = outputTruth_handle.ptr();
215 
216  // create decorators
217  SG::WriteDecorHandle<xAOD::TruthParticleContainer, float > d0Decorator(m_d0DecoratorKey, ctx);
218  SG::WriteDecorHandle<xAOD::TruthParticleContainer, float > z0Decorator(m_z0DecoratorKey, ctx);
219  SG::WriteDecorHandle<xAOD::TruthParticleContainer, float > ptDecorator(m_ptDecoratorKey, ctx);
220 
221 
222  // clear the smearear
223  m_mySmearer->Clear();
224 
225  static const SG::ConstAccessor<float> ptAcc("pt");
226  static const SG::ConstAccessor<float> d0Acc("d0");
227  static const SG::ConstAccessor<float> z0Acc("z0");
228  static const SG::ConstAccessor<float> thetaAcc("theta");
229 
230  int n_input_tracks=0;
231  int n_output_tracks=0;
232  int n_output_broad_tracks=0;
233  int n_output_narrow_tracks=0;
234  ATH_MSG_DEBUG ("Found "<<inputTruth->size()<< " input truth particles");
235  for ( const auto* part : *inputTruth )
236  {
237  double pt = part->pt();// MeV
238  float theta = thetaAcc(*part);
239  float z0 = z0Acc(*part);
240  float d0 = d0Acc(*part);
241  float eta = part->eta();
242  float phi = part->phi();
243  if (part->isNeutral()) continue;
244  if (pt <=0.) continue;
245 
246  ATH_MSG_DEBUG ("===> New Truth: "
247  <<" curv=" << 1./part->pt()
248  <<" phi=" << part->phi()
249  <<" eta=" << part->eta()
250  <<" d0=" << d0Acc(*part)
251  <<" z0=" << z0Acc(*part)
252  <<" pT=" << part->pt()
253  <<" PDGID=" << part->pdgId()
254  <<" status=" << part->status()
255  );
256  if (part->parent()) ATH_MSG_DEBUG (" parent status=" << part->parent()->pdgId());
257 
258  if (std::abs(pt)/1000. > m_inputTracksPtCut) //GeV cut
259  {
260  n_input_tracks++;
261  if (m_enableMonitoring) {
262  hist("track_input_eta")->Fill(eta);
263  hist("track_input_theta")->Fill(theta);
264  hist("track_input_pt" )->Fill(pt/1000.);
265  hist("track_input_phi")->Fill(phi);
266  hist("track_input_z0" )->Fill(z0);
267  hist("track_input_d0" )->Fill(d0);
268  }
269  double qoverPt = part->charge()*1000./pt; //this must be in GeV
270  m_mySmearer->AddTrack(d0,z0,qoverPt,eta,phi); // smearing here
271  n_output_tracks += m_mySmearer->GetNTracks();
272 
273  ATH_MSG_DEBUG ("Looping on output tracks #"<< m_mySmearer->GetNTracks());
274  for (const auto& otrack : m_mySmearer->Tracks)
275  {
276  xAOD::TruthParticle * newtrk = new xAOD::TruthParticle(*part);
277  outputTruth->push_back(newtrk);
278  *newtrk = *part;
279  auto newpt = part->pt();
280  if (m_SigmaScaleFactor !=0) {
281  // set the decorators
282  d0Decorator(*newtrk) = otrack.d0();
283  z0Decorator(*newtrk) = otrack.z0();
284  ptDecorator(*newtrk) = otrack.pt()*1000.; //MeV
285  //TrackParticle has already ::pt(), so the smeared value is in the decorator
286  // and can be accessed by
287  newpt = ptAcc(*newtrk);
288  }
289  if (newpt==0.) continue;
290  ATH_MSG_DEBUG ("Smeared Truth: "
291  <<" curv=" << 1./newpt
292  <<" phi=" << newtrk->phi()
293  <<" eta=" << newtrk->eta()
294  <<" d0=" << d0Acc(*newtrk)
295  <<" z0=" << z0Acc(*newtrk)
296  <<" pT=" << newpt
297  <<" PDGID=" << newtrk->pdgId()
298  <<" status=" << newtrk->status()
299  );
300  if (newtrk->parent()) ATH_MSG_DEBUG (" parent status=" << newtrk->parent()->pdgId());
301 
302  if (m_enableMonitoring) {
303  hist("track_output_eta")->Fill(otrack.eta());
304  hist("track_output_theta")->Fill(otrack.theta());
305  hist("track_output_pt" )->Fill(part->charge()*otrack.pt() );
306  hist("track_output_phi")->Fill(otrack.phi());
307  hist("track_output_z0" )->Fill(otrack.z0() );
308  hist("track_output_d0" )->Fill(otrack.d0() );
309 
310  hist("track_outputcoll_eta")->Fill(newtrk->eta());
311  hist("track_outputcoll_theta")->Fill(thetaAcc(*newtrk));
312  hist("track_outputcoll_pt" )->Fill(part->charge()* newpt/1000.);
313  hist("track_outputcoll_phi")->Fill(newtrk->phi());
314  hist("track_outputcoll_z0" )->Fill(z0Acc(*newtrk));
315  hist("track_outputcoll_d0" )->Fill(d0Acc(*newtrk));
316 
317  hist("track_delta_eta")->Fill(newtrk->eta() - part->eta());
318  hist("track_delta_pt") ->Fill((newpt - part->pt())/1000.);
319  hist("track_delta_crv")->Fill(newtrk->charge()*1000./newpt - ((part->charge()*1000./part->pt())));
320  hist("track_delta_phi")->Fill(newtrk->phi() - part->phi());
321  hist("track_delta_z0" )->Fill(z0Acc(*newtrk) - z0Acc(*part));
322  hist("track_delta_d0" )->Fill(d0Acc(*newtrk) - d0Acc(*part));
323  }
324  } // end of loop
325  }
326  m_mySmearer->Clear(); // clear the smearer after each input track
327  }
328 
329  ATH_MSG_DEBUG ("End of loop track #"<<n_input_tracks<<" ---> "<<" "<< n_output_tracks
330  <<" "<<n_output_narrow_tracks<<" "<<n_output_broad_tracks);
331  if (m_enableMonitoring) {
332  hist("n_input_tracks")->Fill(n_input_tracks);
333  hist("n_output_tracks")->Fill(n_output_tracks);
334  hist("n_output_narrow_tracks")->Fill(n_output_narrow_tracks);
335  hist("n_output_broad_tracks")->Fill(n_output_broad_tracks);
336  }
337  return StatusCode::SUCCESS;
338 }
339 
340 
341 
342 StatusCode EFTrackingSmearingAlg::execute() {
343  ATH_MSG_DEBUG ("Executing " << name() << "...");
344 
345  auto ctx = getContext() ;
346  if (m_smearTruthParticle)
347  return smearTruthParticles(ctx);
348 
349  SG::ReadHandle<xAOD::TrackParticleContainer> inputTracks_handle( m_inputTrackParticleKey, ctx );
350  const xAOD::TrackParticleContainer* inputTracks = inputTracks_handle.cptr();
351  if (not inputTracks) {
352  ATH_MSG_FATAL("Unable to retrieve input ID tacks");
353  return StatusCode::FAILURE;
354  }
355 
356  SG::WriteHandle<xAOD::TrackParticleContainer> outputTracks_handle( m_outputTrackParticleKey, ctx );
357  ATH_CHECK( outputTracks_handle.record( std::make_unique<xAOD::TrackParticleContainer>(), std::make_unique<xAOD::TrackParticleAuxContainer>() ) );
358  auto outputTracks = outputTracks_handle.ptr();
359 
360  // clear the smearear
361  m_mySmearer->Clear();
362 
363  //int trackno=0;
364  int n_input_tracks=0;
365  int n_output_tracks=0;
366  int n_output_broad_tracks=0;
367  int n_output_narrow_tracks=0;
368  ATH_MSG_DEBUG ("Found "<<inputTracks->size()<< " input tracks");
369  for ( const auto* trk : *inputTracks )
370  {
371  // get Cov matrix of input track
372  xAOD::ParametersCovMatrix_t trkcov = trk->definingParametersCovMatrix();
373  auto trkcovvec = trk->definingParametersCovMatrixVec();
374  double theta=trk->theta();
375  double pt = std::sin(theta)/std::abs(trk->qOverP());// MeV
376 
377  ATH_MSG_DEBUG ("===> New Track: "
378  <<" curv=" << 1./trk->pt()
379  <<" phi=" << trk->phi0()
380  <<" eta=" << trk->eta()
381  <<" d0=" << trk->d0()
382  <<" z0=" << trk->z0()
383  <<" pT=" << trk->pt()
384  <<" cov_d0=" << trkcov(Trk::d0,Trk::d0)
385  <<" cov_z0=" << trkcov(Trk::z0,Trk::z0)
386  <<" sigma_d0=" << std::sqrt(std::abs(trkcov(Trk::d0,Trk::d0)))
387  <<" sigma_z0=" << std::sqrt(std::abs(trkcov(Trk::z0,Trk::z0))) );
388 
389 
390 
391  if (std::abs(pt) > m_inputTracksPtCut)
392  {
393  n_input_tracks++;
394  if (m_enableMonitoring) {
395  hist("track_input_eta")->Fill(trk->eta());
396  hist("track_input_theta")->Fill(trk->theta());
397  hist("track_input_pt" )->Fill(pt/1000.);
398  hist("track_input_phi")->Fill(trk->phi0());
399  hist("track_input_z0" )->Fill(trk->z0());
400  hist("track_input_d0" )->Fill(trk->d0());
401 
402  hist("track_input_sigma_theta") ->Fill(std::sqrt(trkcov(Trk::theta,Trk::theta)));
403  hist("track_input_sigma_qOverP")->Fill(std::sqrt(trkcov(Trk::qOverP,Trk::qOverP)));
404  hist("track_input_sigma_phi") ->Fill(std::sqrt(trkcov(Trk::phi0,Trk::phi0)));
405  hist("track_input_sigma_z0") ->Fill(std::sqrt(trkcov(Trk::z0,Trk::z0)));
406  hist("track_input_sigma_d0") ->Fill(std::sqrt(trkcov(Trk::d0,Trk::d0)));
407  }
408 
409  // get Cov matrix of input track
410  auto trkcovvec = trk->definingParametersCovMatrixVec();
411  double qoverPt = trk->charge()*1000./pt; //this must be in GeV
412  m_mySmearer->AddTrack(trk->d0(),trk->z0(),qoverPt,trk->eta(),trk->phi0());
413  n_output_tracks += m_mySmearer->GetNTracks();
414 
415  ATH_MSG_DEBUG ("Looping on output tracks #"<< m_mySmearer->GetNTracks());
416  for (const auto& otrack : m_mySmearer->Tracks)
417  {
418  xAOD::TrackParticle * newtrk = new xAOD::TrackParticle(*trk);
419  outputTracks->push_back(newtrk);
420  *newtrk = *trk;
421 
422  double sintheta=std::sin(otrack.theta());
423  trkcov = trk->definingParametersCovMatrix();
424  auto newtrkcov = trkcov;
425  if (m_SigmaScaleFactor !=0) {
426  // modify the cov matrix
427  ATH_MSG_DEBUG ("Setting parameters in covariance vector");
428  for (int ii=0;ii<5;ii++) for (int jj=0;jj<5;jj++) {
429  // correct the CM to be consistent with the smearing ofthe parameters
430  newtrkcov(ii,jj)=trkcov(ii,jj) * (std::pow(m_SigmaScaleFactor,2) + 1);
431  }
432  trkcovvec.clear();
433  EigenHelpers::eigenMatrixToVector(trkcovvec,newtrkcov,"");
434  newtrk->setDefiningParametersCovMatrixVec(trkcovvec);
435  ATH_MSG_DEBUG ("Setting parameters covariance");
436  newtrk->setDefiningParametersCovMatrix(newtrkcov);
437 
438  /* old method stays for reference, since it's not yet validated
439  for (int ii=0;ii<5;ii++) for (int jj=0;jj<5;jj++) if (ii!=jj) trkcov(ii,jj)=0.0;
440  trkcov(Trk::d0 ,Trk::d0 ) = otrack.sigma_d0() * otrack.sigma_d0() ;
441  trkcov(Trk::z0 ,Trk::z0 ) = otrack.sigma_z0() * otrack.sigma_z0() ;
442  trkcov(Trk::theta ,Trk::theta ) = otrack.sigma_theta()* otrack.sigma_theta();
443  trkcov(Trk::phi0 ,Trk::phi0 ) = otrack.sigma_phi() * otrack.sigma_phi() ;
444  trkcov(Trk::qOverP,Trk::qOverP) = otrack.sigma_curv()/1000. * otrack.sigma_curv()/1000. *sintheta*sintheta;
445 
446  trkcovvec.clear();
447  EigenHelpers::eigenMatrixToVector(trkcovvec,trkcov,"");
448  newtrk->setDefiningParametersCovMatrixVec(trkcovvec);
449  ATH_MSG_DEBUG ("Setting parameters covariance");
450  newtrk->setDefiningParametersCovMatrix(trkcov);
451  */
452 
453  //(float d0, float z0, float phi0, float theta, float qOverP)
454  newtrk->setDefiningParameters(
455  otrack.d0(),
456  otrack.z0(),
457  otrack.phi(),
458  otrack.theta(),
459  (otrack.curv()*sintheta/1000.) //oneOverp in MeV
460  );
461  }
462  xAOD::ParametersCovMatrix_t trkcov_out = newtrk->definingParametersCovMatrix();
463 
464 
465  ATH_MSG_DEBUG ("Smeared Track: "
466  <<" curv=" << 1./newtrk->pt()
467  <<" phi=" << newtrk->phi()
468  <<" eta=" << newtrk->eta()
469  <<" d0=" << newtrk->d0()
470  <<" z0=" << newtrk->z0()
471  <<" pT=" << newtrk->pt()
472  <<" cov_d0=" << trkcov_out(Trk::d0,Trk::d0)
473  <<" cov_z0=" << trkcov_out(Trk::z0,Trk::z0)
474  <<" sigma_d0=" << std::sqrt(std::abs(trkcov_out(Trk::d0,Trk::d0)))
475  <<" sigma_z0=" << std::sqrt(std::abs(trkcov_out(Trk::z0,Trk::z0))) );
476 
477 
478  if (m_enableMonitoring) {
479  hist("track_output_eta")->Fill(otrack.eta());
480  hist("track_output_theta")->Fill(otrack.theta());
481  hist("track_output_pt" )->Fill(trk->charge()*otrack.pt() );
482  hist("track_output_phi")->Fill(otrack.phi());
483  hist("track_output_z0" )->Fill(otrack.z0() );
484  hist("track_output_d0" )->Fill(otrack.d0() );
485 
486  hist("track_outputcoll_eta")->Fill(newtrk->eta());
487  hist("track_outputcoll_theta")->Fill(newtrk->theta());
488  hist("track_outputcoll_pt" )->Fill(trk->charge()*std::sin(newtrk->theta())/(1000.0*newtrk->qOverP()));
489  hist("track_outputcoll_phi")->Fill(newtrk->phi0());
490  hist("track_outputcoll_z0" )->Fill(newtrk->z0());
491  hist("track_outputcoll_d0" )->Fill(newtrk->d0());
492 
493  hist("track_delta_eta")->Fill(newtrk->eta() - trk->eta());
494  hist("track_delta_pt")->Fill((newtrk->pt() - trk->pt())/1000.);
495  hist("track_delta_crv")->Fill((1000.0*newtrk->qOverP()/std::sin(theta))-((1000.0*trk->qOverP())/std::sin(theta)));
496  hist("track_delta_phi")->Fill(newtrk->phi() - trk->phi0());
497  hist("track_delta_z0" )->Fill(newtrk->z0() - trk->z0());
498  hist("track_delta_d0" )->Fill(newtrk->d0() - trk->d0());
499 
500  auto trkcov_original = trk->definingParametersCovMatrix();
501  hist("track_delta_sigma_theta") ->Fill(std::sqrt(trkcov_out(Trk::theta,Trk::theta)) - std::sqrt(trkcov_original(Trk::theta,Trk::theta)));
502  hist("track_delta_sigma_qOverP")->Fill(std::sqrt(trkcov_out(Trk::qOverP,Trk::qOverP)) - std::sqrt(trkcov_original(Trk::qOverP,Trk::qOverP)));
503  hist("track_delta_sigma_phi") ->Fill(std::sqrt(trkcov_out(Trk::phi0,Trk::phi0)) - std::sqrt(trkcov_original(Trk::phi0,Trk::phi0)));
504  hist("track_delta_sigma_z0") ->Fill(std::sqrt(trkcov_out(Trk::z0,Trk::z0)) - std::sqrt(trkcov_original(Trk::z0,Trk::z0)));
505  hist("track_delta_sigma_d0") ->Fill(std::sqrt(trkcov_out(Trk::d0,Trk::d0)) - std::sqrt(trkcov_original(Trk::d0,Trk::d0)));
506 
507 
508  hist("track_outputcoll_sigma_theta") ->Fill(std::sqrt(trkcov(Trk::theta,Trk::theta)));
509  hist("track_outputcoll_sigma_qOverP")->Fill(std::sqrt(trkcov(Trk::qOverP,Trk::qOverP)));
510  hist("track_outputcoll_sigma_phi") ->Fill(std::sqrt(trkcov(Trk::phi0,Trk::phi0)));
511  hist("track_outputcoll_sigma_z0") ->Fill(std::sqrt(trkcov(Trk::z0,Trk::z0)));
512  hist("track_outputcoll_sigma_d0") ->Fill(std::sqrt(trkcov(Trk::d0,Trk::d0)));
513 
514  }
515  // do we need this hack?
516  // hack!!!
517  // Chi2: 1 -> "Broad Fake"
518  // 0 -> "core tracks"
519  //
520  // NDF: index of parent track in input track collection
521  //newtrk->setFitQuality(otrack.FakeFlags(),trackno);
522 
523 
524  }
525 
526  }
527  m_mySmearer->Clear(); // clear teh smearer after each input track
528  //trackno++;
529  }
530 
531  ATH_MSG_DEBUG ("End of loop track #"<<n_input_tracks<<" ---> "<<" "<< n_output_tracks
532  <<" "<<n_output_narrow_tracks<<" "<<n_output_broad_tracks);
533  if (m_enableMonitoring) {
534  hist("n_input_tracks")->Fill(n_input_tracks);
535  hist("n_output_tracks")->Fill(n_output_tracks);
536  hist("n_output_narrow_tracks")->Fill(n_output_narrow_tracks);
537  hist("n_output_broad_tracks")->Fill(n_output_broad_tracks);
538  }
539  return StatusCode::SUCCESS;
540 }
541 
542 
543 
544 
545 StatusCode EFTrackingSmearingAlg::book_histograms()
546 {
547  double b_eta[3]={50,-5,5};
548  CHECK(book(new TH1F("track_input_eta","#eta of input tracks",b_eta[0],b_eta[1],b_eta[2])));
549  CHECK(book(new TH1F("track_input_theta","#Theta of input tracks",50,0.0,4.0)));
550  CHECK(book(new TH1F("track_input_pt" ,"p_T of input tracks",50.0,-10.0,10.0)));
551  CHECK(book(new TH1F("track_input_phi","#phi of input tracks",50,-6.28,6.28)));
552  CHECK(book(new TH1F("track_input_z0" ,"z_0 of input tracks",100,-50,50)));
553  CHECK(book(new TH1F("track_input_d0" ,"d_0 of input tracks",50,-5,5)));
554 
555  CHECK(book(new TH1F("track_input_sigma_theta","#sigma_{#Theta} of input tracks" ,50,0.0,0.001)));
556  CHECK(book(new TH1F("track_input_sigma_qOverP" ,"#sigma_{q/P} of input tracks" ,50.0,0.0,0.0001)));
557  CHECK(book(new TH1F("track_input_sigma_phi" ,"#sigma_{#phi} of input tracks" ,50,0.0,0.002)));
558  CHECK(book(new TH1F("track_input_sigma_z0" ,"#sigma_{z_0} of input tracks" ,100,0.0,5)));
559  CHECK(book(new TH1F("track_input_sigma_d0" ,"#sigma_{d_0} of input tracks" ,50,0.0,5)));
560 
561  int maxtracks=500;
562  CHECK(book(new TH1F("n_input_tracks" ,"Number of input tracks" ,maxtracks,0,maxtracks)));
563  CHECK(book(new TH1F("n_output_tracks" ,"Number of output tracks" ,maxtracks,0,maxtracks)));
564  CHECK(book(new TH1F("n_output_narrow_tracks","Number of output tracks (narrow)" ,maxtracks,0,maxtracks)));
565  CHECK(book(new TH1F("n_output_broad_tracks" ,"Number of output tracks (broad)" ,maxtracks,0,maxtracks)));
566 
567  CHECK(book(new TH1F("track_output_eta","#eta of output tracks",50,-5,5)));
568  CHECK(book(new TH1F("track_output_theta","#Theta of output tracks",50,0.0,4.0)));
569  CHECK(book(new TH1F("track_output_pt" ,"p_T of output tracks [GeV]",50,-10,10)));
570  CHECK(book(new TH1F("track_output_phi","#phi of output tracks",50,-6.28,6.28)));
571  CHECK(book(new TH1F("track_output_z0" ,"z_0 of output tracks",100,-50,50)));
572  CHECK(book(new TH1F("track_output_d0" ,"d_0 of output tracks",50,-5,5)));
573 
574  CHECK(book(new TH1F("track_outputcoll_eta","#eta of output tracks collection",50,-5,5)));
575  CHECK(book(new TH1F("track_outputcoll_theta","#Theta of output tracks collection",50,0.0,4.0)));
576  CHECK(book(new TH1F("track_outputcoll_pt" ,"p_T of output tracks collection [GeV]",50,-10,10)));
577  CHECK(book(new TH1F("track_outputcoll_phi","#phi of output tracks collection",50,-6.28,6.28)));
578  CHECK(book(new TH1F("track_outputcoll_z0" ,"z_0 of output tracks collection",100,-50,50)));
579  CHECK(book(new TH1F("track_outputcoll_d0" ,"d_0 of output tracks collection",50,-5,5)));
580 
581  CHECK(book(new TH1F("track_outputcoll_sigma_theta","#sigma_{#Theta} of output tracks collection" ,50,0.0,0.001)));
582  CHECK(book(new TH1F("track_outputcoll_sigma_qOverP" ,"#sigma_{q/P} of output tracks collection",50.0,0.0,0.0001)));
583  CHECK(book(new TH1F("track_outputcoll_sigma_phi" ,"#sigma_{#phi} of output tracks collection" ,50,0.0,0.002)));
584  CHECK(book(new TH1F("track_outputcoll_sigma_z0" ,"#sigma_{z_0} of output tracks collection" ,100,0.0,5)));
585  CHECK(book(new TH1F("track_outputcoll_sigma_d0" ,"#sigma_{d_0} of output tracks collection" ,50,0.0,5)));
586 
587  CHECK(book(new TH1F("track_delta_sigma_theta", "#sigma_{#Theta} of output tracks collection",100,-0.001,0.001)));
588  CHECK(book(new TH1F("track_delta_sigma_qOverP", "#sigma_{q/P} of output tracks collection",100,-0.0001,0.0001)));
589  CHECK(book(new TH1F("track_delta_sigma_phi", "#sigma_{#phi} of output tracks collection" ,100,-0.002,0.002)));
590  CHECK(book(new TH1F("track_delta_sigma_z0", "#sigma_{z_0} of output tracks collection" ,100,-1.,1.)));
591  CHECK(book(new TH1F("track_delta_sigma_d0", "#sigma_{d_0} of output tracks collection" ,100,-0.5,0.5)));
592 
593  CHECK(book(new TH1F("track_delta_eta","tracks #Delta #eta",50,-1,1)));
594  CHECK(book(new TH1F("track_delta_pt","tracks #Delta #pt [GeV]",50,-2.,2.)));
595  CHECK(book(new TH1F("track_delta_crv" ,"tracks #Delta crv",50,-1,1)));
596  CHECK(book(new TH1F("track_delta_phi","tracks #Delta #phi",50,-0.5,0.5)));
597  CHECK(book(new TH1F("track_delta_z0" ,"tracks #Delta z_0 ",100,-10,10)));
598  CHECK(book(new TH1F("track_delta_d0" ,"tracks #Delta d_0 ",50,-2,2)));
599  return StatusCode::SUCCESS;
600 }
601 
602 
603 
LArG4FSStartPointFilter.part
part
Definition: LArG4FSStartPointFilter.py:21
EFTrackingSmearingAlg::smearTruthParticles
StatusCode smearTruthParticles(const EventContext &ctx)
Definition: EFTrackingSmearingAlg.cxx:204
xAOD::TrackParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition: TrackParticle_v1.cxx:73
xAOD::TruthParticle_v1::parent
const TruthParticle_v1 * parent(size_t i=0) const
Retrieve the i-th mother (TruthParticle) of this TruthParticle.
Definition: TruthParticle_v1.cxx:131
Trk::y
@ y
Definition: ParamDefs.h:56
Trk::py
@ py
Definition: ParamDefs.h:60
EFTrackingSmearingAlg::m_smearedTrackEfficiency_d0low_LRT
DoubleProperty m_smearedTrackEfficiency_d0low_LRT
Definition: EFTrackingSmearingAlg.h:60
EFTrackingSmearingAlg::m_outputTracksPtCut
DoubleProperty m_outputTracksPtCut
Definition: EFTrackingSmearingAlg.h:52
EFTrackingSmearingAlg::m_EnableFakes
BooleanProperty m_EnableFakes
Definition: EFTrackingSmearingAlg.h:63
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
Amg::MatrixX
Eigen::Matrix< double, Eigen::Dynamic, Eigen::Dynamic > MatrixX
Dynamic Matrix - dynamic allocation.
Definition: EventPrimitives.h:27
EFTrackingSmearingAlg::m_z0DecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_z0DecoratorKey
Definition: EFTrackingSmearingAlg.h:43
Amg::compress
void compress(const AmgSymMatrix(N) &covMatrix, std::vector< float > &vec)
Definition: EventPrimitivesHelpers.h:56
Trk::z
@ z
global position (cartesian)
Definition: ParamDefs.h:57
phi
Scalar phi() const
phi method
Definition: AmgMatrixBasePlugin.h:67
ATH_MSG_INFO
#define ATH_MSG_INFO(x)
Definition: AthMsgStreamMacros.h:31
Trk::locX
@ locX
Definition: ParamDefs.h:37
Trk::locY
@ locY
local cartesian
Definition: ParamDefs.h:38
SG::ReadHandle::cptr
const_pointer_type cptr()
Dereference the pointer.
AthHistogramming::book
StatusCode book(const TH1 &hist, const std::string &tDir="", const std::string &stream="")
Simplify the booking and registering (into THistSvc) of histograms.
Definition: AthHistogramming.h:303
EFTrackingSmearingAlg::book_histograms
StatusCode book_histograms()
Definition: EFTrackingSmearingAlg.cxx:545
EFTrackingSmearingAlg::m_IncludeFakesInResolutionCalculation
BooleanProperty m_IncludeFakesInResolutionCalculation
Definition: EFTrackingSmearingAlg.h:65
eta
Scalar eta() const
pseudorapidity method
Definition: AmgMatrixBasePlugin.h:83
SG::ReadHandle
Definition: StoreGate/StoreGate/ReadHandle.h:70
EFTrackingSmearingAlg::m_smearTruthParticle
BooleanProperty m_smearTruthParticle
Definition: EFTrackingSmearingAlg.h:56
xAOD::TrackParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition: TrackParticle_v1.cxx:77
EFTrackingSmearingAlg::m_ptDecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_ptDecoratorKey
Definition: EFTrackingSmearingAlg.h:44
Trk::ParamDefs
ParamDefs
Definition: ParamDefs.h:32
EventPrimitivesHelpers.h
theta
Scalar theta() const
theta method
Definition: AmgMatrixBasePlugin.h:75
Trk::locRPhi
@ locRPhi
Definition: ParamDefs.h:40
Amg::CalculateCompressedSize
constexpr int CalculateCompressedSize(int n)
Definition: EventPrimitivesHelpers.h:51
plotBeamSpotVxVal.cov
cov
Definition: plotBeamSpotVxVal.py:201
test_pyathena.pt
pt
Definition: test_pyathena.py:11
xAOD::TrackParticle_v1::z0
float z0() const
Returns the parameter.
xAOD::TrackParticle_v1::setDefiningParameters
void setDefiningParameters(float d0, float z0, float phi0, float theta, float qOverP)
Set the defining parameters.
Definition: TrackParticle_v1.cxx:177
Trk::z0
@ z0
Definition: ParamDefs.h:64
AthCommonMsg< Algorithm >::msgLvl
bool msgLvl(const MSG::Level lvl) const
Definition: AthCommonMsg.h:30
Trk::loc2
@ loc2
generic first and second local coordinate
Definition: ParamDefs.h:35
athena.value
value
Definition: athena.py:124
JetTiledMap::N
@ N
Definition: TiledEtaPhiMap.h:44
vec
std::vector< size_t > vec
Definition: CombinationsGeneratorTest.cxx:12
SG::ConstAccessor< float >
EFTrackingSmearingAlg::m_mySmearer
std::unique_ptr< FakeTrackSmearer > m_mySmearer
Definition: EFTrackingSmearingAlg.h:72
Trk::locR
@ locR
Definition: ParamDefs.h:44
EFTrackingSmearingAlg::m_outputTruthParticleKey
SG::WriteHandleKey< xAOD::TruthParticleContainer > m_outputTruthParticleKey
Definition: EFTrackingSmearingAlg.h:38
SG::VarHandleKey::key
const std::string & key() const
Return the StoreGate ID for the referenced object.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:141
Trk::pz
@ pz
global momentum (cartesian)
Definition: ParamDefs.h:61
ParamDefs.h
Trk::u
@ u
Enums for curvilinear frames.
Definition: ParamDefs.h:77
AmgSymMatrix
#define AmgSymMatrix(dim)
Definition: EventPrimitives.h:50
xAOD::TrackParticle_v1::d0
float d0() const
Returns the parameter.
EFTrackingSmearingAlg::m_d0DecoratorKey
SG::WriteDecorHandleKey< xAOD::TruthParticleContainer > m_d0DecoratorKey
Definition: EFTrackingSmearingAlg.h:42
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition: Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
Trk::iEta
@ iEta
(old readout) will be skipped
Definition: ParamDefs.h:48
EFTrackingSmearingAlg::m_UseResolutionPtCutOff
BooleanProperty m_UseResolutionPtCutOff
Definition: EFTrackingSmearingAlg.h:49
EFTrackingSmearingAlg::m_inputTrackParticleKey
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_inputTrackParticleKey
Definition: EFTrackingSmearingAlg.h:29
Trk::locZ
@ locZ
local cylindrical
Definition: ParamDefs.h:42
Trk::distEta
@ distEta
readout for silicon
Definition: ParamDefs.h:51
python.TrigInDetConfig.inputTracks
inputTracks
Definition: TrigInDetConfig.py:183
lumiFormat.i
int i
Definition: lumiFormat.py:85
EFTrackingSmearingAlg::m_inputTruthParticleKey
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_inputTruthParticleKey
Definition: EFTrackingSmearingAlg.h:35
Trk::theta
@ theta
Definition: ParamDefs.h:66
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
SG::WriteHandle::ptr
pointer_type ptr()
Dereference the pointer.
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
xAOD::covMatrix
covMatrix
Definition: TrackMeasurement_v1.cxx:19
xAOD::TruthParticle_v1
Class describing a truth particle in the MC record.
Definition: TruthParticle_v1.h:37
Trk::driftRadius
@ driftRadius
trt, straws
Definition: ParamDefs.h:53
SG::WriteDecorHandle
Handle class for adding a decoration to an object.
Definition: StoreGate/StoreGate/WriteDecorHandle.h:100
xAOD::TruthParticle
TruthParticle_v1 TruthParticle
Typedef to implementation.
Definition: Event/xAOD/xAODTruth/xAODTruth/TruthParticle.h:15
TRT::Track::d0
@ d0
Definition: InnerDetector/InDetCalibEvent/TRT_CalibData/TRT_CalibData/TrackInfo.h:62
Trk::px
@ px
Definition: ParamDefs.h:59
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
CHECK
#define CHECK(...)
Evaluate an expression and check for errors.
Definition: Control/AthenaKernel/AthenaKernel/errorcheck.h:422
EFTrackingSmearingAlg::m_outputTrackParticleKey
SG::WriteHandleKey< xAOD::TrackParticleContainer > m_outputTrackParticleKey
Definition: EFTrackingSmearingAlg.h:32
beamspotnt.rows
list rows
Definition: bin/beamspotnt.py:1112
SG::VarHandleKey::initialize
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
Definition: AthToolSupport/AsgDataHandles/Root/VarHandleKey.cxx:103
TRT::Track::z0
@ z0
Definition: InnerDetector/InDetCalibEvent/TRT_CalibData/TRT_CalibData/TrackInfo.h:63
AthHistogramAlgorithm
Definition: AthHistogramAlgorithm.h:32
DataVector
Derived DataVector<T>.
Definition: DataVector.h:794
Trk::distPhi
@ distPhi
Definition: ParamDefs.h:50
EigenHelpers::eigenMatrixToVector
static void eigenMatrixToVector(VECTOR &vec, COVARIANCE &cov, const char *)
Helper fn to get raw data (vec<float>) from the covariance.
Definition: EigenHelpers.h:30
EigenHelpers
Definition: EFTrackingSmearingAlg.cxx:85
xAOD::TrackParticle_v1::phi0
float phi0() const
Returns the parameter, which has range to .
Definition: TrackParticle_v1.cxx:158
Trk
Ensure that the ATLAS eigen extensions are properly loaded.
Definition: FakeTrackBuilder.h:9
xAOD::TrackParticle_v1::qOverP
float qOverP() const
Returns the parameter.
xAOD::TrackParticle_v1::definingParametersCovMatrix
const ParametersCovMatrix_t definingParametersCovMatrix() const
Returns the 5x5 symmetric matrix containing the defining parameters covariance matrix.
Definition: TrackParticle_v1.cxx:246
Amg
Definition of ATLAS Math & Geometry primitives (Amg)
Definition: AmgStringHelpers.h:19
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
Trk::d0
@ d0
Definition: ParamDefs.h:63
EFTrackingSmearingAlg::m_parameterizedTrackEfficiency
BooleanProperty m_parameterizedTrackEfficiency
Definition: EFTrackingSmearingAlg.h:54
Trk::trkMass
@ trkMass
Extended perigee: mass.
Definition: ParamDefs.h:81
Trk::locPhi
@ locPhi
local polar
Definition: ParamDefs.h:45
xAOD::TruthParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition: TruthParticle_v1.cxx:174
xAOD::TrackParticle_v1::setDefiningParametersCovMatrix
void setDefiningParametersCovMatrix(const ParametersCovMatrix_t &cov)
Set the defining parameters covariance matrix.
Definition: TrackParticle_v1.cxx:213
EFTrackingSmearingAlg::m_inputTracksPtCut
DoubleProperty m_inputTracksPtCut
Definition: EFTrackingSmearingAlg.h:51
EFTrackingSmearingAlg::initialize
virtual StatusCode initialize() override
Definition: EFTrackingSmearingAlg.cxx:101
Trk::iPhi
@ iPhi
Definition: ParamDefs.h:47
EFTrackingSmearingAlg::m_smearedTrackEfficiency
DoubleProperty m_smearedTrackEfficiency
Definition: EFTrackingSmearingAlg.h:53
EFTrackingSmearingAlg::m_RandomSeed
LongLongProperty m_RandomSeed
Definition: EFTrackingSmearingAlg.h:70
xAOD::TrackParticle_v1::setDefiningParametersCovMatrixVec
void setDefiningParametersCovMatrixVec(const std::vector< float > &cov)
Definition: TrackParticle_v1.cxx:460
SG::WriteHandle
Definition: StoreGate/StoreGate/WriteHandle.h:76
EFTrackingSmearingAlg::m_UseCoinToss
BooleanProperty m_UseCoinToss
Definition: EFTrackingSmearingAlg.h:67
xAOD::TruthParticle_v1::status
int status() const
Status code.
VertexContainer.h
SG::WriteHandle::record
StatusCode record(std::unique_ptr< T > data)
Record a const object to the store.
VECTOR
TVectorD VECTOR
Definition: sincos_poly.cxx:20
Trk::locPhiR
@ locPhiR
Definition: ParamDefs.h:41
xAOD::TruthParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle.
Definition: TruthParticle_v1.cxx:181
EFTrackingSmearingAlg::~EFTrackingSmearingAlg
virtual ~EFTrackingSmearingAlg()
Definition: EFTrackingSmearingAlg.cxx:98
EFTrackingSmearingAlg::execute
virtual StatusCode execute() override
Definition: EFTrackingSmearingAlg.cxx:342
DEBUG
#define DEBUG
Definition: page_access.h:11
Trk::qOverP
@ qOverP
perigee
Definition: ParamDefs.h:67
AthHistogramming::hist
TH1 * hist(const std::string &histName, const std::string &tDir="", const std::string &stream="")
Simplify the retrieval of registered histograms of any type.
Definition: AthHistogramming.cxx:198
Trk::phi
@ phi
Definition: ParamDefs.h:75
python.TrigEgammaMonitorHelper.TH1F
def TH1F(name, title, nxbins, bins_par2, bins_par3=None, path='', **kwargs)
Definition: TrigEgammaMonitorHelper.py:24
EFTrackingSmearingAlg::m_smearedTrackEfficiency_d0high_LRT
DoubleProperty m_smearedTrackEfficiency_d0high_LRT
Definition: EFTrackingSmearingAlg.h:59
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
ConstAccessor.h
Helper class to provide constant type-safe access to aux data.
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
Trk::x
@ x
Definition: ParamDefs.h:55
EFTrackingSmearingAlg.h
EFTrackingSmearingAlg::m_enableMonitoring
BooleanProperty m_enableMonitoring
Definition: EFTrackingSmearingAlg.h:55
EFTrackingSmearingAlg::m_SigmaScaleFactor
DoubleProperty m_SigmaScaleFactor
Definition: EFTrackingSmearingAlg.h:48
xAOD::TruthParticle_v1::pdgId
int pdgId() const
PDG ID code.
Trk::loc1
@ loc1
Definition: ParamDefs.h:34
EFTrackingSmearingAlg::EFTrackingSmearingAlg
EFTrackingSmearingAlg(const std::string &name, ISvcLocator *pSvcLocator)
Definition: EFTrackingSmearingAlg.cxx:94
EFTrackingSmearingAlg::m_parameterizedTrackEfficiency_LRT
BooleanProperty m_parameterizedTrackEfficiency_LRT
Definition: EFTrackingSmearingAlg.h:58
pow
constexpr int pow(int base, int exp) noexcept
Definition: ap_fixedTest.cxx:15
xAOD::TrackParticle_v1::theta
float theta() const
Returns the parameter, which has range 0 to .
EFTrackingSmearingAlg::m_FakeKillerEnable
BooleanProperty m_FakeKillerEnable
Definition: EFTrackingSmearingAlg.h:64
EFTrackingSmearingAlg::m_SetResolutionPtCutOff
DoubleProperty m_SetResolutionPtCutOff
Definition: EFTrackingSmearingAlg.h:50
Trk::phi0
@ phi0
Definition: ParamDefs.h:65
xAOD::TruthParticle_v1::charge
double charge() const
Physical charge.
Trk::v
@ v
Definition: ParamDefs.h:78
xAOD::TrackParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle (has range to .)
Generated on Mon Dec 23 2024 21:09:46 for ATLAS Offline Software by   doxygen 1.8.18