Loading [MathJax]/extensions/tex2jax.js
ATLAS Offline Software
•All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Properties Friends Macros Modules Pages
NewVrtSecInclusiveTool.cxx
Go to the documentation of this file.
1 /*
2  Copyright (C) 2002-2025 CERN for the benefit of the ATLAS collaboration
3 */
7 
8 
9 // Header include
10 #include "NewVrtSecInclusiveTool/NewVrtSecInclusiveTool.h"
11 #include "VxSecVertex/VxSecVertexInfo.h"
12 #include "PathResolver/PathResolver.h"
13 #include "CxxUtils/checker_macros.h"
14 
15 #include "GaudiKernel/ITHistSvc.h"
16 #include "GaudiKernel/ConcurrencyFlags.h"
17 #include "TH1.h"
18 #include "TH2.h"
19 #include "TTree.h"
20 #include "TMath.h"
21 #include "TFile.h"
22 #include "MVAUtils/BDT.h"
23 #include "xAODEgamma/ElectronxAODHelpers.h"
24 
25 
26 namespace Rec {
27 
28 //
29 //Constructor--------------------------------------------------------------
30 NewVrtSecInclusiveTool::NewVrtSecInclusiveTool(const std::string& type,
31  const std::string& name,
32  const IInterface* parent):
33  AthAlgTool(type,name,parent),
34  m_SV2T_BDT(nullptr),
35  m_instanceName(name),
36  m_is_selected("is_selected"),
37  m_is_svtrk_final("is_svtrk_final"),
38  m_pt_wrtSV("pt_wrtSV"),
39  m_eta_wrtSV("eta_wrtSV"),
40  m_phi_wrtSV("phi_wrtSV"),
41  m_d0_wrtSV("d0_wrtSV"),
42  m_z0_wrtSV("z0_wrtSV"),
43  m_errP_wrtSV("errP_wrtSV"),
44  m_errd0_wrtSV("errd0_wrtSV"),
45  m_errz0_wrtSV("errz0_wrtSV"),
46  m_chi2_toSV("chi2_toSV")
47  {
48 //
49 // Declare additional interface
50 //
51  declareInterface< IVrtInclusive >(this);
52 //
53  m_massPi = Trk::ParticleMasses::mass[Trk::pion];
54  m_massP = Trk::ParticleMasses::mass[Trk::proton];
55  m_massE = Trk::ParticleMasses::mass[Trk::electron];
56  m_massK0 = Trk::ParticleMasses::mass[Trk::k0];
57  m_massLam = 1115.683 ;
58 
59  }
60 
61 //Destructor---------------------------------------------------------------
62  NewVrtSecInclusiveTool::~NewVrtSecInclusiveTool(){
63  ATH_MSG_DEBUG("NewVrtSecInclusiveTool destructor called");
64  }
65 
66 //Initialize---------------------------------------------------------------
67  StatusCode NewVrtSecInclusiveTool::initialize(){
68  ATH_MSG_DEBUG( "Initialising NewVrtSecInclusiveTool" );
69  ATH_CHECK( m_extrapolator.retrieve() );
70  ATH_CHECK( m_trackToVertexTool.retrieve() );
71  ATH_CHECK(m_beamSpotKey.initialize());
72  ATH_CHECK( m_fitSvc.retrieve() );
73  ATH_MSG_DEBUG("NewVrtSecInclusiveTool TrkVKalVrtFitter found");
74 
75  m_is_selected = SG::AuxElement::Decorator<char>("is_selected"+m_augString);
76  m_is_svtrk_final = SG::AuxElement::Decorator<char>("is_svtrk_final"+m_augString);
77  m_pt_wrtSV = SG::AuxElement::Decorator<float>("pt_wrtSV"+m_augString);
78  m_eta_wrtSV = SG::AuxElement::Decorator<float>("eta_wrtSV"+m_augString);
79  m_phi_wrtSV = SG::AuxElement::Decorator<float>("phi_wrtSV"+m_augString);
80  m_d0_wrtSV = SG::AuxElement::Decorator<float>("d0_wrtSV"+m_augString);
81  m_z0_wrtSV = SG::AuxElement::Decorator<float>("z0_wrtSV"+m_augString);
82  m_errP_wrtSV = SG::AuxElement::Decorator<float>("errP_wrtSV"+m_augString);
83  m_errd0_wrtSV = SG::AuxElement::Decorator<float>("errd0_wrtSV"+m_augString);
84  m_errz0_wrtSV = SG::AuxElement::Decorator<float>("errz0_wrtSV"+m_augString);
85  m_chi2_toSV = SG::AuxElement::Decorator<float>("chi2_toSV"+m_augString);
86 
87 //------------------------------------------
88 //
89  if(m_fillHist){
90  if (Gaudi::Concurrency::ConcurrencyFlags::numThreads() > 1) {
91  ATH_MSG_FATAL("Filling histograms not supported in MT jobs.");
92  return StatusCode::FAILURE;
93  }
94 
95  SmartIF<ITHistSvc> hist_root{service("THistSvc")};
96  if( !hist_root ) ATH_MSG_DEBUG("Could not find THistSvc service");
97  else ATH_MSG_DEBUG("NewVrtSecInclusiveTool Histograms found");
98  std::string histDir;
99  histDir="/NVSI/stat/MultiSVrt"+m_instanceName+"/";
100 
101  m_h = std::make_unique<Hists>();
102  ATH_CHECK( m_h->book (*hist_root, histDir) );
103 
104  m_w_1 = 1.;
105  }
106 
107 //--------------------------------------------------------
108  //std::string fileName="NewVrtSecInclusiveTool/Fake2TrVertexReject.MVA.v01.root"; ///For local calibration file
109  //std::string rootFilePath = PathResolver::find_file(fileName, "DATAPATH"); ///
110  std::string rootFilePath = PathResolver::find_calib_file("NewVrtSecInclusiveTool/"+m_calibFileName);
111  TFile* rootFile = TFile::Open(rootFilePath.c_str(), "READ");
112  if (!rootFile) {
113  ATH_MSG_FATAL("Could not retrieve root file: " << m_calibFileName);
114  return StatusCode::FAILURE;
115  }
116  TTree * training = (TTree*)rootFile->Get("BDT");
117  m_SV2T_BDT = std::make_unique<MVAUtils::BDT>(training);
118 //--------------------------------------------------------
119  return StatusCode::SUCCESS;
120 
121  }
122 
123 
124  StatusCode NewVrtSecInclusiveTool::Hists::book (ITHistSvc& histSvc,
125  const std::string& histDir)
126  {
127  m_hb_massPiPi = new TH1D("massPiPi"," mass PiPi",200,0., 4000.);
128  m_hb_massPiPi1 = new TH1D("massPiPi1"," mass PiPi",200,0., 4000.);
129  m_hb_massPPi = new TH1D("massPPi"," massPPi", 100,1000., 1250.);
130  m_hb_massEE = new TH1D("massEE"," massEE", 100,0., 200.);
131  m_hb_nvrt2 = new TH1D("nvrt2"," vertices2", 50,0., 50.);
132  m_hb_ratio = new TH1D("ratio"," ratio", 51,0., 1.02);
133  m_hb_totmass = new TH1D("totmass"," totmass", 250,0., 10000.);
134  m_hb_impact = new TH1D("impact", " impact", 100,0., 20.);
135  m_hb_impactR = new TH1D("impactR"," impactR", 400,-30., 70.);
136  m_hb_impactZ = new TH1D("impactZ"," impactZ", 100,-30., 70.);
137  m_hb_impactRZ = new TH2D("impactRZ"," impactRZ", 40,-10., 10., 60, -30.,30. );
138  m_hb_trkD0 = new TH1D("trkD0"," d0 of tracks", 100, 0., 10.);
139  m_hb_trkZ = new TH1D("trkZ"," Z of tracks", 120,-60., 60.);
140  m_hb_r2d = new TH1D("r2interact","Interaction radius 2tr selected", 150,0., 150.);
141  m_hb_ntrksel = new TH1F("NTrkSel","Number of selected tracks", 200,0., 200.);
142  m_hb_ntrkInput = new TH1F("NTrkInput","Number of provided tracks", 200,0., 1000.);
143  m_hb_trkSelect = new TH1F("TrkSelect","Track selection efficiency", 15,0., 15.);
144  m_hb_signif3D = new TH1D("signif3D"," Signif3D for initial 2tr vertex", 140,-20., 50.);
145  m_hb_sig3DTot = new TH1D("sig3dcommon"," Signif3D for common vertex", 140,-20., 50.);
146  m_hb_sig3D1tr = new TH1D("sig3D1tr","Signif3D for 1tr vertices", 140,-20., 50.);
147  m_hb_sig3D2tr = new TH1D("sig3D2tr","Signif3D for 2tr single vertex", 140,-20., 50.);
148  m_hb_sig3DNtr = new TH1D("sig3DNtr","Signif3D for many-tr single vertex", 140,-20., 50.);
149  m_hb_goodvrtN = new TH1F("goodvrtN","Number of good vertices", 20,0., 20.);
150  m_hb_goodvrt1N = new TH1F("goodvrt1N","Number of good 1-track vertices", 20,0., 20.);
151  m_hb_distVV = new TH1D("distvv","Vertex-Vertex dist", 100,0., 20.);
152  m_hb_diffPS = new TH1D("diffPS","Primary-Secondary assoc", 200,-20., 20.);
153  m_hb_rawVrtN = new TH1F("rawVrtN","Number of raw vertices multivertex case", 20, 0., 20.);
154  m_hb_cosSVMom = new TH1F("cosSVMom","SV-PV vs SV momentum ", 100, 0., 1.);
155  m_hb_etaSV = new TH1F("etaSV"," Eta of SV-PV ", 100, -5., 5.);
156  m_hb_fakeSVBDT = new TH1F("fakeSVBDT"," BDT for fake SV rejection", 100, -1., 1.);
157 
158  ATH_CHECK( histSvc.regHist(histDir+"massPiPi", m_hb_massPiPi) );
159  ATH_CHECK( histSvc.regHist(histDir+"massPiPi1", m_hb_massPiPi1) );
160  ATH_CHECK( histSvc.regHist(histDir+"massPPi", m_hb_massPPi) );
161  ATH_CHECK( histSvc.regHist(histDir+"massEE", m_hb_massEE ) );
162  ATH_CHECK( histSvc.regHist(histDir+"nvrt2", m_hb_nvrt2) );
163  ATH_CHECK( histSvc.regHist(histDir+"ratio", m_hb_ratio) );
164  ATH_CHECK( histSvc.regHist(histDir+"totmass", m_hb_totmass) );
165  ATH_CHECK( histSvc.regHist(histDir+"impact", m_hb_impact) );
166  ATH_CHECK( histSvc.regHist(histDir+"impactR", m_hb_impactR) );
167  ATH_CHECK( histSvc.regHist(histDir+"impactZ", m_hb_impactZ) );
168  ATH_CHECK( histSvc.regHist(histDir+"impactRZ", m_hb_impactRZ) );
169  ATH_CHECK( histSvc.regHist(histDir+"trkD0", m_hb_trkD0) );
170  ATH_CHECK( histSvc.regHist(histDir+"trkZ", m_hb_trkZ) );
171  ATH_CHECK( histSvc.regHist(histDir+"r2interact",m_hb_r2d) );
172  ATH_CHECK( histSvc.regHist(histDir+"NTrkSel", m_hb_ntrksel) );
173  ATH_CHECK( histSvc.regHist(histDir+"NTrkInput", m_hb_ntrkInput) );
174  ATH_CHECK( histSvc.regHist(histDir+"TrkSelect", m_hb_trkSelect) );
175  ATH_CHECK( histSvc.regHist(histDir+"signif3D", m_hb_signif3D) );
176  ATH_CHECK( histSvc.regHist(histDir+"sig3dcommon", m_hb_sig3DTot) );
177  ATH_CHECK( histSvc.regHist(histDir+"sig3D1tr", m_hb_sig3D1tr) );
178  ATH_CHECK( histSvc.regHist(histDir+"sig3D2tr", m_hb_sig3D2tr) );
179  ATH_CHECK( histSvc.regHist(histDir+"sig3DNtr", m_hb_sig3DNtr) );
180  ATH_CHECK( histSvc.regHist(histDir+"goodvrtN", m_hb_goodvrtN) );
181  ATH_CHECK( histSvc.regHist(histDir+"goodvrt1N", m_hb_goodvrt1N) );
182  ATH_CHECK( histSvc.regHist(histDir+"distVV", m_hb_distVV) );
183  ATH_CHECK( histSvc.regHist(histDir+"diffPS", m_hb_diffPS) );
184  ATH_CHECK( histSvc.regHist(histDir+"rawVrtN", m_hb_rawVrtN) );
185  ATH_CHECK( histSvc.regHist(histDir+"cosSVMom", m_hb_cosSVMom) );
186  ATH_CHECK( histSvc.regHist(histDir+"etaSV", m_hb_etaSV) );
187  ATH_CHECK( histSvc.regHist(histDir+"fakeSVBDT", m_hb_fakeSVBDT) );
188 
189  m_tuple = new TTree("Vertices","Vertices");
190  ATH_CHECK( histSvc.regTree(histDir, m_tuple) );
191  m_curTup=new DevTuple();
192  m_tuple->Branch("ntrk", &m_curTup->nTrk, "ntrk/I");
193  m_tuple->Branch("pttrk", &m_curTup->pttrk, "pttrk[ntrk]/F");
194  m_tuple->Branch("d0trk", &m_curTup->d0trk, "d0trk[ntrk]/F");
195  m_tuple->Branch("Sig3D", &m_curTup->Sig3D, "Sig3D[ntrk]/F");
196  m_tuple->Branch("dRdZrat", &m_curTup->dRdZrat, "dRdZrat[ntrk]/F");
197  m_tuple->Branch("idHF", &m_curTup->idHF, "idHF[ntrk]/I");
198  m_tuple->Branch("trkTRT", &m_curTup->trkTRT, "trkTRT[ntrk]/I");
199  m_tuple->Branch("etatrk", &m_curTup->etatrk, "etatrk[ntrk]/F");
200 
201  m_tuple->Branch("n2Vrt", &m_curTup->n2Vrt, "n2Vrt/I");
202  m_tuple->Branch("VrtTrkHF", &m_curTup->VrtTrkHF, "VrtTrkHF[n2Vrt]/I");
203  m_tuple->Branch("VrtTrkI", &m_curTup->VrtTrkI, "VrtTrkI[n2Vrt]/I");
204  m_tuple->Branch("VrtCh", &m_curTup->VrtCh, "VrtCh[n2Vrt]/I");
205  m_tuple->Branch("VrtDist2D", &m_curTup->VrtDist2D, "VrtDist2D[n2Vrt]/F");
206  m_tuple->Branch("VrtSig3D", &m_curTup->VrtSig3D, "VrtSig3D[n2Vrt]/F");
207  m_tuple->Branch("VrtSig2D", &m_curTup->VrtSig2D, "VrtSig2D[n2Vrt]/F");
208  m_tuple->Branch("VrtM", &m_curTup->VrtM, "VrtM[n2Vrt]/F");
209  m_tuple->Branch("VrtZ", &m_curTup->VrtZ, "VrtZ[n2Vrt]/F");
210  m_tuple->Branch("VrtPt", &m_curTup->VrtPt, "VrtPt[n2Vrt]/F");
211  m_tuple->Branch("VrtEta", &m_curTup->VrtEta, "VrtEta[n2Vrt]/F");
212  m_tuple->Branch("VrtIBL", &m_curTup->VrtIBL, "VrtIBL[n2Vrt]/I");
213  m_tuple->Branch("VrtBL", &m_curTup->VrtBL, "VrtBL[n2Vrt]/I");
214  m_tuple->Branch("VrtCosSPM", &m_curTup->VrtCosSPM, "VrtCosSPM[n2Vrt]/F");
215  m_tuple->Branch("VMinPtT", &m_curTup->VMinPtT, "VMinPtT[n2Vrt]/F");
216  m_tuple->Branch("VMinS3DT", &m_curTup->VMinS3DT, "VMinS3DT[n2Vrt]/F");
217  m_tuple->Branch("VMaxS3DT", &m_curTup->VMaxS3DT, "VMaxS3DT[n2Vrt]/F");
218  m_tuple->Branch("VrtProb", &m_curTup->VrtProb, "VrtProb[n2Vrt]/F");
219  m_tuple->Branch("VrtHR1", &m_curTup->VrtHR1, "VrtHR1[n2Vrt]/F");
220  m_tuple->Branch("VrtHR2", &m_curTup->VrtHR2, "VrtHR2[n2Vrt]/F");
221  m_tuple->Branch("VrtBDT", &m_curTup->VrtBDT, "VrtBDT[n2Vrt]/F");
222  m_tuple->Branch("VrtDZ", &m_curTup->VrtDZ, "VrtDZ[n2Vrt]/F");
223  m_tuple->Branch("VrtDisk", &m_curTup->VrtDisk, "VrtDisk[n2Vrt]/I");
224  m_tuple->Branch("VSigMat", &m_curTup->VSigMat, "VSigMat[n2Vrt]/F");
225 
226  m_tuple->Branch("nNVrt", &m_curTup->nNVrt, "nNVrt/I");
227  m_tuple->Branch("NVrtTrk", &m_curTup->NVrtTrk, "NVrtTrk[nNVrt]/I");
228  m_tuple->Branch("NVrtTrkHF", &m_curTup->NVrtTrkHF, "NVrtTrkHF[nNVrt]/I");
229  m_tuple->Branch("NVrtTrkI", &m_curTup->NVrtTrkI, "NVrtTrkI[nNVrt]/I");
230  m_tuple->Branch("NVrtCh", &m_curTup->NVrtCh, "NVrtCh[nNVrt]/I");
231  m_tuple->Branch("NVrtDist2D", &m_curTup->NVrtDist2D, "NVrtDist2D[nNVrt]/F");
232  m_tuple->Branch("NVrtSig3D", &m_curTup->NVrtSig3D, "NVrtSig3D[nNVrt]/F");
233  m_tuple->Branch("NVrtSig2D", &m_curTup->NVrtSig2D, "NVrtSig2D[nNVrt]/F");
234  m_tuple->Branch("NVrtM", &m_curTup->NVrtM, "NVrtM[nNVrt]/F");
235  m_tuple->Branch("NVrtPt", &m_curTup->NVrtPt, "NVrtPt[nNVrt]/F");
236  m_tuple->Branch("NVrtEta", &m_curTup->NVrtEta, "NVrtEta[nNVrt]/F");
237  m_tuple->Branch("NVrtIBL", &m_curTup->NVrtIBL, "NVrtIBL[nNVrt]/I");
238  m_tuple->Branch("NVrtBL", &m_curTup->NVrtBL, "NVrtBL[nNVrt]/I");
239  m_tuple->Branch("NVrtCosSPM", &m_curTup->NVrtCosSPM, "NVrtCosSPM[nNVrt]/F");
240  m_tuple->Branch("NVMinPtT", &m_curTup->NVMinPtT, "NVMinPtT[nNVrt]/F");
241  m_tuple->Branch("NVMinS3DT", &m_curTup->NVMinS3DT, "NVMinS3DT[nNVrt]/F");
242  m_tuple->Branch("NVrtProb", &m_curTup->NVrtProb, "NVrtProb[nNVrt]/F");
243  m_tuple->Branch("NVrtBDT", &m_curTup->NVrtBDT, "NVrtBDT[nNVrt]/F");
244  m_tuple->Branch("NVrtHR1", &m_curTup->NVrtHR1, "NVrtHR1[nNVrt]/F");
245  m_tuple->Branch("NVrtHR2", &m_curTup->NVrtHR2, "NVrtHR2[nNVrt]/F");
246 
247  return StatusCode::SUCCESS;
248  }
249 
250 
251  StatusCode NewVrtSecInclusiveTool::finalize()
252  {
253  ATH_MSG_DEBUG("NewVrtSecInclusiveTool finalize()");
254  return StatusCode::SUCCESS;
255  }
256 
257 
258 
259  std::unique_ptr<Trk::VxSecVertexInfo> NewVrtSecInclusiveTool::findAllVertices (
260  const std::vector<const xAOD::TrackParticle*> & inpTrk,
261  const xAOD::Vertex & primVrt ) const
262  {
263  std::vector<xAOD::Vertex*> listVrtSec(0);
264 
265  if(m_fillHist) {
266  Hists& h = getHists();
267  h.m_curTup->nTrk=0;
268  h.m_curTup->n2Vrt=0;
269  h.m_curTup->nNVrt=0;
270  };
271 
272  workVectorArrxAOD * tmpVectxAOD=new workVectorArrxAOD();
273  tmpVectxAOD->inpTrk.resize(inpTrk.size());
274  std::copy(inpTrk.begin(),inpTrk.end(), tmpVectxAOD->inpTrk.begin());
275  SG::ReadCondHandle<InDet::BeamSpotData> beamSpotHandle { m_beamSpotKey };
276  tmpVectxAOD->beamX=beamSpotHandle->beamPos().x();
277  tmpVectxAOD->beamY=beamSpotHandle->beamPos().y();
278  tmpVectxAOD->beamZ=beamSpotHandle->beamPos().z();
279  tmpVectxAOD->tanBeamTiltX=tan(beamSpotHandle->beamTilt(0));
280  tmpVectxAOD->tanBeamTiltY=tan(beamSpotHandle->beamTilt(1));
281 
282  compatibilityGraph_t compatibilityGraph;
283  listVrtSec = getVrtSecMulti(tmpVectxAOD,primVrt,compatibilityGraph);
284 
285  for (const auto trk : tmpVectxAOD->listSelTracks) {
286  // Mark the track as selected
287  m_is_selected(*trk) = true;
288  if (trk->trackFitter() == xAOD::GaussianSumFilter) m_is_selected(*xAOD::EgammaHelpers::getOriginalTrackParticleFromGSF(trk)) = true;
289  }
290  delete tmpVectxAOD;
291 
292  for (const auto &vrt : listVrtSec) {
293  for (const auto &trk : vrt->trackParticleLinks()) {
294 
295  // Mark the track as a final track
296  m_is_svtrk_final(**trk) = true;
297 
298  // Get the perigee of the track at the vertex
299  ATH_MSG_VERBOSE(" > " << __FUNCTION__ << ": > Track index " << (*trk)->index() << ": Get the perigee of the track at the vertex." );
300  auto sv_perigee = m_trackToVertexTool->perigeeAtVertex(Gaudi::Hive::currentContext(), **trk, vrt->position() );
301  if( !sv_perigee ) {
302  ATH_MSG_WARNING(" > " << __FUNCTION__ << ": > Track index " << (*trk)->index() << ": Failed in obtaining the SV perigee!" );
303  }
304 
305  float qOverP_wrtSV = sv_perigee ? sv_perigee->parameters() [Trk::qOverP] : -FLT_MAX;
306  float theta_wrtSV = sv_perigee ? sv_perigee->parameters() [Trk::theta] : -FLT_MAX;
307  float p_wrtSV = sv_perigee ? 1.0 / std::abs( qOverP_wrtSV ) : -FLT_MAX;
308  float pt_wrtSV = sv_perigee ? p_wrtSV * sin( theta_wrtSV ) : -FLT_MAX;
309  float eta_wrtSV = sv_perigee ? -log( tan( theta_wrtSV/2. ) ) : -FLT_MAX;
310  float phi_wrtSV = sv_perigee ? sv_perigee->parameters() [Trk::phi] : -FLT_MAX;
311  float d0_wrtSV = sv_perigee ? sv_perigee->parameters() [Trk::d0] : -FLT_MAX;
312  float z0_wrtSV = sv_perigee ? sv_perigee->parameters() [Trk::z0] : -FLT_MAX;
313  float errd0_wrtSV = sv_perigee ? std::sqrt((*sv_perigee->covariance())( Trk::d0, Trk::d0 )) : -FLT_MAX;
314  float errz0_wrtSV = sv_perigee ? std::sqrt((*sv_perigee->covariance())( Trk::z0, Trk::z0 )) : -FLT_MAX;
315  float errP_wrtSV = sv_perigee ? std::sqrt((*sv_perigee->covariance())( Trk::qOverP, Trk::qOverP )) : -FLT_MAX;
316 
317  m_pt_wrtSV(**trk) = pt_wrtSV;
318  m_eta_wrtSV(**trk) = eta_wrtSV;
319  m_phi_wrtSV(**trk) = phi_wrtSV;
320  m_d0_wrtSV(**trk) = d0_wrtSV;
321  m_z0_wrtSV(**trk) = z0_wrtSV;
322  m_errP_wrtSV(**trk) = errP_wrtSV;
323  m_errd0_wrtSV(**trk) = errd0_wrtSV;
324  m_errz0_wrtSV(**trk) = errz0_wrtSV;
325  }
326  }
327 
328  std::unique_ptr<Trk::VxSecVertexInfo> res = std::make_unique<Trk::VxSecVertexInfo>(Trk::VxSecVertexInfo(listVrtSec));
329 
330  if(m_fillHist){
331  Hists& h = getHists();
332  h.m_tuple->Fill();
333  };
334 
335  lockDecorations (inpTrk);
336 
337  return res;
338  }
339 
340 
341  void NewVrtSecInclusiveTool::lockDecorations (const std::vector<const xAOD::TrackParticle*> & inpTrk) const
342  {
343  // We may have track from several containers. Use this to keep track
344  // of which ones we've processed. We only expect a few distinct ones,
345  // so just use a vector.
346  std::vector<const SG::AuxVectorData*> containers;
347  containers.reserve (16);
348 
349  for (const xAOD::TrackParticle* t : inpTrk) {
350  const SG::AuxVectorData* c = t->container();
351  if (c && std::find (containers.begin(), containers.end(), c) == containers.end())
352  {
353  containers.push_back (c);
354  // Ok because we just made these decorations.
355  SG::AuxVectorData* c_nc ATLAS_THREAD_SAFE = const_cast<SG::AuxVectorData*> (c);
356  c_nc->lockDecoration (m_is_selected.auxid());
357  c_nc->lockDecoration (m_is_svtrk_final.auxid());
358  c_nc->lockDecoration (m_pt_wrtSV.auxid());
359  c_nc->lockDecoration (m_eta_wrtSV.auxid());
360  c_nc->lockDecoration (m_phi_wrtSV.auxid());
361  c_nc->lockDecoration (m_d0_wrtSV.auxid());
362  c_nc->lockDecoration (m_z0_wrtSV.auxid());
363  c_nc->lockDecoration (m_errP_wrtSV.auxid());
364  c_nc->lockDecoration (m_errd0_wrtSV.auxid());
365  c_nc->lockDecoration (m_errz0_wrtSV.auxid());
366  c_nc->lockDecoration (m_chi2_toSV.auxid());
367  }
368  }
369  }
370 
371 
372  NewVrtSecInclusiveTool::Hists&
373  NewVrtSecInclusiveTool::getHists() const
374  {
375  // We earlier checked that no more than one thread is being used.
376  Hists* h ATLAS_THREAD_SAFE = m_h.get();
377  return *h;
378  }
379 
380 
381 } // end Rec namespace
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtBDT
float NVrtBDT[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:285
Rec::NewVrtSecInclusiveTool::~NewVrtSecInclusiveTool
virtual ~NewVrtSecInclusiveTool()
Definition: NewVrtSecInclusiveTool.cxx:62
PathResolver::find_calib_file
static std::string find_calib_file(const std::string &logical_file_name)
Definition: PathResolver.cxx:384
Rec::NewVrtSecInclusiveTool::Hists::book
StatusCode book(ITHistSvc &histSvc, const std::string &histDir)
Definition: NewVrtSecInclusiveTool.cxx:124
BDT.h
Rec::NewVrtSecInclusiveTool::DevTuple::pttrk
float pttrk[maxNTrk]
Definition: NewVrtSecInclusiveTool.h:235
Rec::NewVrtSecInclusiveTool::DevTuple::NVMinS3DT
float NVMinS3DT[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:279
Rec::workVectorArrxAOD::beamY
double beamY
Definition: GNNVertexFitterTool.h:43
Trk::proton
@ proton
Definition: ParticleHypothesis.h:31
ATH_MSG_FATAL
#define ATH_MSG_FATAL(x)
Definition: AthMsgStreamMacros.h:34
Trk::k0
@ k0
Definition: ParticleHypothesis.h:35
Rec::NewVrtSecInclusiveTool::Hists::m_hb_trkZ
TH1D * m_hb_trkZ
Definition: NewVrtSecInclusiveTool.h:120
Rec::NewVrtSecInclusiveTool::DevTuple::VSigMat
float VSigMat[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:265
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtHR2
float NVrtHR2[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:287
Rec::NewVrtSecInclusiveTool::DevTuple::VrtBDT
float VrtBDT[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:256
Rec::NewVrtSecInclusiveTool::Hists
Definition: NewVrtSecInclusiveTool.h:105
SG::ReadCondHandle
Definition: ReadCondHandle.h:44
Rec::NewVrtSecInclusiveTool::getVrtSecMulti
std::vector< xAOD::Vertex * > getVrtSecMulti(workVectorArrxAOD *inpParticlesxAOD, const xAOD::Vertex &primVrt, compatibilityGraph_t &compatibilityGraph) const
Definition: VrtSecMulti.cxx:29
find
std::string find(const std::string &s)
return a remapped string
Definition: hcg.cxx:135
Rec::NewVrtSecInclusiveTool::Hists::m_hb_trkD0
TH1D * m_hb_trkD0
Definition: NewVrtSecInclusiveTool.h:119
Rec::NewVrtSecInclusiveTool::Hists::m_hb_ntrksel
TH1F * m_hb_ntrksel
Definition: NewVrtSecInclusiveTool.h:121
Rec::NewVrtSecInclusiveTool::m_eta_wrtSV
SG::AuxElement::Decorator< float > m_eta_wrtSV
Definition: NewVrtSecInclusiveTool.h:211
Rec::NewVrtSecInclusiveTool::m_h
std::unique_ptr< Hists > m_h
Definition: NewVrtSecInclusiveTool.h:141
Rec::NewVrtSecInclusiveTool::m_instanceName
std::string m_instanceName
Definition: NewVrtSecInclusiveTool.h:206
Rec::NewVrtSecInclusiveTool::NewVrtSecInclusiveTool
NewVrtSecInclusiveTool(const std::string &type, const std::string &name, const IInterface *parent)
Definition: NewVrtSecInclusiveTool.cxx:30
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtSig2D
float NVrtSig2D[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:283
Rec::NewVrtSecInclusiveTool::Hists::m_hb_diffPS
TH1D * m_hb_diffPS
Definition: NewVrtSecInclusiveTool.h:132
Rec::workVectorArrxAOD::listSelTracks
std::vector< const xAOD::TrackParticle * > listSelTracks
Definition: GNNVertexFitterTool.h:41
Rec::NewVrtSecInclusiveTool::Hists::m_hb_signif3D
TH1D * m_hb_signif3D
Definition: NewVrtSecInclusiveTool.h:126
Rec::NewVrtSecInclusiveTool::DevTuple::n2Vrt
int n2Vrt
Definition: NewVrtSecInclusiveTool.h:242
Rec::NewVrtSecInclusiveTool::DevTuple::VrtZ
float VrtZ[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:253
Rec::NewVrtSecInclusiveTool::Hists::m_hb_sig3D2tr
TH1D * m_hb_sig3D2tr
Definition: NewVrtSecInclusiveTool.h:134
Rec::NewVrtSecInclusiveTool::DevTuple::NVMinPtT
float NVMinPtT[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:278
Rec::NewVrtSecInclusiveTool::DevTuple::VMinPtT
float VMinPtT[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:262
ElectronxAODHelpers.h
Trk::z0
@ z0
Definition: ParamDefs.h:64
Rec::NewVrtSecInclusiveTool::DevTuple::VMinS3DT
float VMinS3DT[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:263
Rec::workVectorArrxAOD::tanBeamTiltX
double tanBeamTiltX
Definition: GNNVertexFitterTool.h:45
Rec::NewVrtSecInclusiveTool::DevTuple::VrtEta
float VrtEta[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:255
SG::AuxVectorData::lockDecoration
void lockDecoration(SG::auxid_t auxid)
Explicitly lock a decoration.
Definition: AuxVectorData.cxx:687
Rec::NewVrtSecInclusiveTool::DevTuple::VrtProb
float VrtProb[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:257
read_hist_ntuple.t
t
Definition: read_hist_ntuple.py:5
ATH_MSG_VERBOSE
#define ATH_MSG_VERBOSE(x)
Definition: AthMsgStreamMacros.h:28
Rec::NewVrtSecInclusiveTool::m_trackToVertexTool
ToolHandle< Reco::ITrackToVertex > m_trackToVertexTool
Definition: NewVrtSecInclusiveTool.h:197
Rec::NewVrtSecInclusiveTool::m_is_svtrk_final
SG::AuxElement::Decorator< char > m_is_svtrk_final
Definition: NewVrtSecInclusiveTool.h:209
Rec::NewVrtSecInclusiveTool::DevTuple
Definition: NewVrtSecInclusiveTool.h:231
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtPt
float NVrtPt[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:275
Rec::NewVrtSecInclusiveTool::Hists::m_hb_trkSelect
TH1F * m_hb_trkSelect
Definition: NewVrtSecInclusiveTool.h:123
Rec::NewVrtSecInclusiveTool::DevTuple::idHF
int idHF[maxNTrk]
Definition: NewVrtSecInclusiveTool.h:240
python.CaloAddPedShiftConfig.type
type
Definition: CaloAddPedShiftConfig.py:42
Rec::NewVrtSecInclusiveTool::Hists::m_hb_etaSV
TH1F * m_hb_etaSV
Definition: NewVrtSecInclusiveTool.h:138
Rec::NewVrtSecInclusiveTool::m_beamSpotKey
SG::ReadCondHandleKey< InDet::BeamSpotData > m_beamSpotKey
Definition: NewVrtSecInclusiveTool.h:193
Rec::NewVrtSecInclusiveTool::DevTuple::nTrk
int nTrk
Definition: NewVrtSecInclusiveTool.h:234
Rec::NewVrtSecInclusiveTool::Hists::m_hb_fakeSVBDT
TH1F * m_hb_fakeSVBDT
Definition: NewVrtSecInclusiveTool.h:139
Rec::NewVrtSecInclusiveTool::finalize
StatusCode finalize()
Definition: NewVrtSecInclusiveTool.cxx:251
Rec::NewVrtSecInclusiveTool::m_fitSvc
ToolHandle< Trk::TrkVKalVrtFitter > m_fitSvc
Definition: NewVrtSecInclusiveTool.h:196
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtDist2D
float NVrtDist2D[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:281
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtIBL
int NVrtIBL[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:272
Rec::NewVrtSecInclusiveTool::Hists::m_hb_impactRZ
TH2D * m_hb_impactRZ
Definition: NewVrtSecInclusiveTool.h:118
Rec::NewVrtSecInclusiveTool::Hists::m_hb_goodvrtN
TH1F * m_hb_goodvrtN
Definition: NewVrtSecInclusiveTool.h:129
Rec::NewVrtSecInclusiveTool::DevTuple::VrtM
float VrtM[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:252
Rec::NewVrtSecInclusiveTool::Hists::m_hb_totmass
TH1D * m_hb_totmass
Definition: NewVrtSecInclusiveTool.h:115
Rec::NewVrtSecInclusiveTool::DevTuple::trkTRT
int trkTRT[maxNTrk]
Definition: NewVrtSecInclusiveTool.h:241
SG::Decorator
Helper class to provide type-safe access to aux data.
Definition: Decorator.h:59
Rec::NewVrtSecInclusiveTool::Hists::m_hb_sig3DTot
TH1D * m_hb_sig3DTot
Definition: NewVrtSecInclusiveTool.h:128
Rec::NewVrtSecInclusiveTool::m_z0_wrtSV
SG::AuxElement::Decorator< float > m_z0_wrtSV
Definition: NewVrtSecInclusiveTool.h:214
IsoCloseByCorrectionTest.containers
containers
Associate the close-by pflow objects and the calorimeter clusters.
Definition: IsoCloseByCorrectionTest.py:82
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtTrkHF
int NVrtTrkHF[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:269
Rec::NewVrtSecInclusiveTool::m_calibFileName
Gaudi::Property< std::string > m_calibFileName
Definition: NewVrtSecInclusiveTool.h:189
Rec::NewVrtSecInclusiveTool::lockDecorations
void lockDecorations(const std::vector< const xAOD::TrackParticle * > &inpTrk) const
Definition: NewVrtSecInclusiveTool.cxx:341
Rec::NewVrtSecInclusiveTool::Hists::m_hb_ntrkInput
TH1F * m_hb_ntrkInput
Definition: NewVrtSecInclusiveTool.h:122
Rec
Name: MuonSpContainer.h Package : offline/Reconstruction/MuonIdentification/muonEvent.
Definition: FakeTrackBuilder.h:10
Rec::NewVrtSecInclusiveTool::m_d0_wrtSV
SG::AuxElement::Decorator< float > m_d0_wrtSV
Definition: NewVrtSecInclusiveTool.h:213
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtM
float NVrtM[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:274
Trk::theta
@ theta
Definition: ParamDefs.h:66
EL::StatusCode
::StatusCode StatusCode
StatusCode definition for legacy code.
Definition: PhysicsAnalysis/D3PDTools/EventLoop/EventLoop/StatusCode.h:22
Trk::electron
@ electron
Definition: ParticleHypothesis.h:27
extractSporadic.h
list h
Definition: extractSporadic.py:97
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition: AthMsgStreamMacros.h:29
Rec::NewVrtSecInclusiveTool::DevTuple::VrtDist2D
float VrtDist2D[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:249
Rec::NewVrtSecInclusiveTool::m_phi_wrtSV
SG::AuxElement::Decorator< float > m_phi_wrtSV
Definition: NewVrtSecInclusiveTool.h:212
xAOD::GaussianSumFilter
@ GaussianSumFilter
Tracks from Gaussian Sum Filter.
Definition: TrackingPrimitives.h:50
Rec::NewVrtSecInclusiveTool::m_massPi
double m_massPi
Definition: NewVrtSecInclusiveTool.h:201
Trk::pion
@ pion
Definition: ParticleHypothesis.h:29
Rec::NewVrtSecInclusiveTool::m_massP
double m_massP
Definition: NewVrtSecInclusiveTool.h:202
res
std::pair< std::vector< unsigned int >, bool > res
Definition: JetGroupProductTest.cxx:11
Rec::NewVrtSecInclusiveTool::m_fillHist
Gaudi::Property< bool > m_fillHist
Definition: NewVrtSecInclusiveTool.h:173
Rec::NewVrtSecInclusiveTool::DevTuple::VrtSig3D
float VrtSig3D[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:250
test_pyathena.parent
parent
Definition: test_pyathena.py:15
MuonSegmentReaderConfig.histSvc
histSvc
Definition: MuonSegmentReaderConfig.py:96
VxSecVertexInfo.h
Rec::workVectorArrxAOD::tanBeamTiltY
double tanBeamTiltY
Definition: GNNVertexFitterTool.h:46
ATH_CHECK
#define ATH_CHECK
Definition: AthCheckMacros.h:40
makeComparison.rootFile
rootFile
Definition: makeComparison.py:27
Rec::NewVrtSecInclusiveTool::DevTuple::d0trk
float d0trk[maxNTrk]
Definition: NewVrtSecInclusiveTool.h:236
drawFromPickle.tan
tan
Definition: drawFromPickle.py:36
Rec::NewVrtSecInclusiveTool::DevTuple::VrtTrkHF
int VrtTrkHF[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:243
Rec::NewVrtSecInclusiveTool::m_augString
Gaudi::Property< std::string > m_augString
Definition: NewVrtSecInclusiveTool.h:199
Rec::NewVrtSecInclusiveTool::Hists::m_tuple
TTree * m_tuple
Definition: NewVrtSecInclusiveTool.h:107
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtBL
int NVrtBL[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:273
Rec::NewVrtSecInclusiveTool::Hists::m_hb_sig3DNtr
TH1D * m_hb_sig3DNtr
Definition: NewVrtSecInclusiveTool.h:135
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtCh
int NVrtCh[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:271
Rec::NewVrtSecInclusiveTool::m_massK0
double m_massK0
Definition: NewVrtSecInclusiveTool.h:204
Rec::NewVrtSecInclusiveTool::DevTuple::VrtSig2D
float VrtSig2D[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:251
Rec::NewVrtSecInclusiveTool::Hists::m_hb_impactZ
TH1D * m_hb_impactZ
Definition: NewVrtSecInclusiveTool.h:124
Trk::ParticleMasses::mass
constexpr double mass[PARTICLEHYPOTHESES]
the array of masses
Definition: ParticleHypothesis.h:53
Rec::NewVrtSecInclusiveTool::m_errP_wrtSV
SG::AuxElement::Decorator< float > m_errP_wrtSV
Definition: NewVrtSecInclusiveTool.h:215
PathResolver.h
Rec::workVectorArrxAOD
Definition: GNNVertexFitterTool.h:40
Rec::NewVrtSecInclusiveTool::DevTuple::VMaxS3DT
float VMaxS3DT[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:264
name
std::string name
Definition: Control/AthContainers/Root/debug.cxx:228
Trk::d0
@ d0
Definition: ParamDefs.h:63
Rec::NewVrtSecInclusiveTool::DevTuple::VrtTrkI
int VrtTrkI[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:244
Rec::NewVrtSecInclusiveTool::findAllVertices
std::unique_ptr< Trk::VxSecVertexInfo > findAllVertices(const std::vector< const xAOD::TrackParticle * > &inputTracks, const xAOD::Vertex &primaryVertex) const final
Definition: NewVrtSecInclusiveTool.cxx:259
Rec::NewVrtSecInclusiveTool::Hists::m_hb_impactR
TH1D * m_hb_impactR
Definition: NewVrtSecInclusiveTool.h:117
Rec::NewVrtSecInclusiveTool::m_is_selected
SG::AuxElement::Decorator< char > m_is_selected
Definition: NewVrtSecInclusiveTool.h:208
SG::CondHandleKey::initialize
StatusCode initialize(bool used=true)
Trk::VxSecVertexInfo
Definition: VxSecVertexInfo.h:63
Rec::workVectorArrxAOD::inpTrk
std::vector< const xAOD::TrackParticle * > inpTrk
Definition: NewVrtSecInclusiveTool.h:70
Rec::NewVrtSecInclusiveTool::DevTuple::VrtCh
int VrtCh[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:245
xAOD::EgammaHelpers::getOriginalTrackParticleFromGSF
const xAOD::TrackParticle * getOriginalTrackParticleFromGSF(const xAOD::TrackParticle *trkPar)
Helper function for getting the "Original" Track Particle (i.e before GSF) via the GSF Track Particle...
Definition: ElectronxAODHelpers.cxx:22
Rec::NewVrtSecInclusiveTool::m_chi2_toSV
SG::AuxElement::Decorator< float > m_chi2_toSV
Definition: NewVrtSecInclusiveTool.h:218
Rec::NewVrtSecInclusiveTool::Hists::m_hb_massPPi
TH1D * m_hb_massPPi
Definition: NewVrtSecInclusiveTool.h:111
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtTrkI
int NVrtTrkI[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:270
Rec::NewVrtSecInclusiveTool::Hists::m_hb_ratio
TH1D * m_hb_ratio
Definition: NewVrtSecInclusiveTool.h:114
Rec::NewVrtSecInclusiveTool::DevTuple::dRdZrat
float dRdZrat[maxNTrk]
Definition: NewVrtSecInclusiveTool.h:239
Rec::NewVrtSecInclusiveTool::Hists::m_hb_massPiPi
TH1D * m_hb_massPiPi
Definition: NewVrtSecInclusiveTool.h:109
Rec::NewVrtSecInclusiveTool::Hists::m_hb_massEE
TH1D * m_hb_massEE
Definition: NewVrtSecInclusiveTool.h:112
NewVrtSecInclusiveTool.h
Rec::NewVrtSecInclusiveTool::Hists::m_hb_distVV
TH1D * m_hb_distVV
Definition: NewVrtSecInclusiveTool.h:131
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtHR1
float NVrtHR1[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:286
Rec::NewVrtSecInclusiveTool::DevTuple::VrtPt
float VrtPt[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:254
Rec::NewVrtSecInclusiveTool::m_errd0_wrtSV
SG::AuxElement::Decorator< float > m_errd0_wrtSV
Definition: NewVrtSecInclusiveTool.h:216
Rec::NewVrtSecInclusiveTool::DevTuple::etatrk
float etatrk[maxNTrk]
Definition: NewVrtSecInclusiveTool.h:237
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtEta
float NVrtEta[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:276
Rec::NewVrtSecInclusiveTool::DevTuple::nNVrt
int nNVrt
Definition: NewVrtSecInclusiveTool.h:267
Rec::NewVrtSecInclusiveTool::compatibilityGraph_t
boost::adjacency_list< boost::listS, boost::vecS, boost::undirectedS > compatibilityGraph_t
Definition: NewVrtSecInclusiveTool.h:308
Rec::NewVrtSecInclusiveTool::m_w_1
double m_w_1
Definition: NewVrtSecInclusiveTool.h:103
Rec::NewVrtSecInclusiveTool::initialize
StatusCode initialize()
Definition: NewVrtSecInclusiveTool.cxx:67
Rec::NewVrtSecInclusiveTool::m_extrapolator
ToolHandle< Trk::IExtrapolator > m_extrapolator
Definition: NewVrtSecInclusiveTool.h:195
Rec::NewVrtSecInclusiveTool::Hists::m_curTup
DevTuple * m_curTup
Definition: NewVrtSecInclusiveTool.h:108
Rec::NewVrtSecInclusiveTool::Hists::m_hb_sig3D1tr
TH1D * m_hb_sig3D1tr
Definition: NewVrtSecInclusiveTool.h:133
Rec::NewVrtSecInclusiveTool::DevTuple::VrtHR2
float VrtHR2[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:259
Rec::NewVrtSecInclusiveTool::Hists::m_hb_rawVrtN
TH1F * m_hb_rawVrtN
Definition: NewVrtSecInclusiveTool.h:136
Rec::NewVrtSecInclusiveTool::DevTuple::VrtCosSPM
float VrtCosSPM[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:261
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtCosSPM
float NVrtCosSPM[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:277
h
xAOD::Vertex_v1
Class describing a Vertex.
Definition: Vertex_v1.h:42
python.BackTrackingConfig.numThreads
int numThreads
Definition: BackTrackingConfig.py:61
ATH_MSG_WARNING
#define ATH_MSG_WARNING(x)
Definition: AthMsgStreamMacros.h:32
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtTrk
int NVrtTrk[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:268
Rec::NewVrtSecInclusiveTool::getHists
Hists & getHists() const
Definition: NewVrtSecInclusiveTool.cxx:373
Trk::qOverP
@ qOverP
perigee
Definition: ParamDefs.h:67
python.CaloCondTools.log
log
Definition: CaloCondTools.py:20
Rec::NewVrtSecInclusiveTool::m_pt_wrtSV
SG::AuxElement::Decorator< float > m_pt_wrtSV
Definition: NewVrtSecInclusiveTool.h:210
Rec::NewVrtSecInclusiveTool::DevTuple::VrtIBL
int VrtIBL[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:246
Rec::NewVrtSecInclusiveTool::DevTuple::Sig3D
float Sig3D[maxNTrk]
Definition: NewVrtSecInclusiveTool.h:238
Rec::NewVrtSecInclusiveTool::DevTuple::VrtBL
int VrtBL[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:247
Rec::NewVrtSecInclusiveTool::Hists::m_hb_massPiPi1
TH1D * m_hb_massPiPi1
Definition: NewVrtSecInclusiveTool.h:110
Rec::NewVrtSecInclusiveTool::Hists::m_hb_cosSVMom
TH1F * m_hb_cosSVMom
Definition: NewVrtSecInclusiveTool.h:137
Trk::phi
@ phi
Definition: ParamDefs.h:75
Rec::NewVrtSecInclusiveTool::m_errz0_wrtSV
SG::AuxElement::Decorator< float > m_errz0_wrtSV
Definition: NewVrtSecInclusiveTool.h:217
python.TrigEgammaMonitorHelper.TH1F
def TH1F(name, title, nxbins, bins_par2, bins_par3=None, path='', **kwargs)
Definition: TrigEgammaMonitorHelper.py:24
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtProb
float NVrtProb[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:284
Rec::NewVrtSecInclusiveTool::DevTuple::NVrtSig3D
float NVrtSig3D[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:282
SG::AuxVectorData
Manage lookup of vectors of auxiliary data.
Definition: AuxVectorData.h:168
calibdata.copy
bool copy
Definition: calibdata.py:27
ATLAS_THREAD_SAFE
#define ATLAS_THREAD_SAFE
Definition: checker_macros.h:211
Rec::NewVrtSecInclusiveTool::Hists::m_hb_nvrt2
TH1D * m_hb_nvrt2
Definition: NewVrtSecInclusiveTool.h:113
xAOD::TrackParticle_v1
Class describing a TrackParticle.
Definition: TrackParticle_v1.h:43
Rec::NewVrtSecInclusiveTool::DevTuple::VrtDZ
float VrtDZ[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:260
Rec::workVectorArrxAOD::beamX
double beamX
Definition: GNNVertexFitterTool.h:42
drawFromPickle.sin
sin
Definition: drawFromPickle.py:36
AthAlgTool
Definition: AthAlgTool.h:26
Rec::NewVrtSecInclusiveTool::m_SV2T_BDT
std::unique_ptr< MVAUtils::BDT > m_SV2T_BDT
Definition: NewVrtSecInclusiveTool.h:191
Rec::NewVrtSecInclusiveTool::DevTuple::VrtDisk
int VrtDisk[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:248
Rec::NewVrtSecInclusiveTool::Hists::m_hb_impact
TH1D * m_hb_impact
Definition: NewVrtSecInclusiveTool.h:116
checker_macros.h
Define macros for attributes used to control the static checker.
Rec::NewVrtSecInclusiveTool::DevTuple::VrtHR1
float VrtHR1[maxNVrt]
Definition: NewVrtSecInclusiveTool.h:258
Rec::NewVrtSecInclusiveTool::Hists::m_hb_goodvrt1N
TH1F * m_hb_goodvrt1N
Definition: NewVrtSecInclusiveTool.h:130
python.compressB64.c
def c
Definition: compressB64.py:93
Rec::NewVrtSecInclusiveTool::m_massE
double m_massE
Definition: NewVrtSecInclusiveTool.h:203
Rec::workVectorArrxAOD::beamZ
double beamZ
Definition: GNNVertexFitterTool.h:44
Rec::NewVrtSecInclusiveTool::m_massLam
double m_massLam
Definition: NewVrtSecInclusiveTool.h:205
Rec::NewVrtSecInclusiveTool::Hists::m_hb_r2d
TH1D * m_hb_r2d
Definition: NewVrtSecInclusiveTool.h:125
Generated on Fri Mar 14 2025 21:16:07 for ATLAS Offline Software by   doxygen 1.8.18