ATLAS Offline Software
Loading...
Searching...
No Matches
SecVertexTruthMatchAlg.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#include "TrackingAnalysisAlgorithms/SecVertexTruthMatchAlg.h"
6#include "InDetSecVtxTruthMatchTool/InDetSecVtxTruthMatchTool.h"
7#include "xAODTruth/TruthParticle.h"
8
9#include "TMath.h"
10#include "TH1.h"
11#include "TEfficiency.h"
12
13const float GeV = 1000.;
14
15namespace CP {
16
17 SecVertexTruthMatchAlg::SecVertexTruthMatchAlg( const std::string& name, ISvcLocator* svcLoc )
18 : EL::AnaAlgorithm( name, svcLoc ) {}
19
20 StatusCode SecVertexTruthMatchAlg::initialize() {
21
22 // Initializing Keys
23 ATH_CHECK(m_secVtxContainerKey.initialize());
24 ATH_CHECK(m_truthVtxContainerKey.initialize());
25 ATH_CHECK(m_trackParticleContainerKey.initialize());
26
27 // Retrieving the tool
28 ATH_CHECK(m_matchTool.retrieve());
29
30 if(m_writeHistograms) {
31 std::vector<std::string> recoTypes{"All", "Matched", "Merged", "Fake", "Split", "Other"};
32 std::vector<std::string> truthTypes{"Inclusive", "Reconstructable", "Accepted", "Seeded", "Reconstructed", "ReconstructedSplit"};
33 // define SM origin categories if SM origin tracking is enabled
34 std::vector<std::string> smOriginTypes;
35 if(m_doSMOrigin) {
36 smOriginTypes = {"FakeOrigin", "Pileup", "KshortDecay", "StrangeMesonDecay", "LambdaDecay",
37 "StrangeBaryonDecay", "TauDecay", "GammaConversion", "OtherDecay",
38 "HadronicInteraction", "OtherSecondary", "BHadronDecay", "DHadronDecay",
39 "Fragmentation", "OtherOrigin", "Signal"};
40 }
41
42 //determine histogram ranges depending if we are in normal or MuSA mode
43 //total bin counts stay the same for simplicity -- MuSA has less precision in general
44 float maxX = m_doMuSA ? 8000 : 500;
45 float maxY = m_doMuSA ? 10000 : 500;
46 float maxZ = m_doMuSA ? 10000 : 1500;
47 float maxLxy = m_doMuSA ? 8000 : 500;
48 float maxR = m_doMuSA ? 8000 : 600;
49 float mind0 = m_doMuSA ? 2000 : 100;
50 float maxd0 = m_doMuSA ? 2000 : 100;
51 float maxTrackd0 = m_doMuSA ? 3000 : 300;
52 float maxTrackz0 = m_doMuSA ? 5000 : 500;
53 float maxErrd0 = m_doMuSA ? 300 : 30;
54 float maxErrz0 = m_doMuSA ? 500 : 50;
55
56 float maxResR = m_doMuSA ? 2000 : 20;
57 float maxResZ = m_doMuSA ? 2000 : 20;
58
59 ANA_CHECK (book(TH1F("RecoVertex/matchType", "Vertex Match Type", 65, -0.5, 64.5)));
60 if (m_doSMOrigin) {
61 ANA_CHECK (book(TH1F("RecoVertex/smOriginType", "Vertex SM Origin Type", 65537, -0.5, 65536.5)));
62 }
63
64
65 for(const auto& recoType : recoTypes) {
66 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_x").c_str(), "Reco vertex x [mm]", 1000, -maxX, maxX)));
67 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_y").c_str(), "Reco vertex y [mm]", 1000, -maxY, maxY)));
68 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_z").c_str(), "Reco vertex z [mm]", 1000, -maxZ, maxZ)));
69 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Lxy").c_str(), "Reco vertex L_{xy} [mm]", 500, 0, maxLxy)));
70 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_pT").c_str(), "Reco vertex p_{T} [GeV]", 100, 0, 100)));
71 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_eta").c_str(), "Reco vertex #eta", 100, -5, 5)));
72 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_phi").c_str(), "Reco vertex #phi", 100, -TMath::Pi(), TMath::Pi())));
73 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_mass").c_str(), "Reco vertex mass [GeV]", 500, 0, 100)));
74 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_mu").c_str(), "Reco vertex Red. Mass [GeV]", 500, 0, 100)));
75 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_chi2").c_str(), "Reco vertex recoChi2", 100, 0, 10)));
76 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_dir").c_str(), "Reco vertex recoDirection", 100, -1, 1)));
77 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_charge").c_str(), "Reco vertex recoCharge", 20, -10, 10)));
78 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_H").c_str(), "Reco vertex H [GeV]", 100, 0, 100)));
79 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_HT").c_str(), "Reco vertex Mass [GeV]", 100, 0, 100)));
80 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_minOpAng").c_str(), "Reco vertex minOpAng", 100, -1, 1)));
81 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_maxOpAng").c_str(), "Reco vertex maxOpAng", 100, -1, 1)));
82 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_maxdR").c_str(), "Reco vertex maxDR", 100, 0, 10)));
83 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_mind0").c_str(), "Reco vertex min d0 [mm]", 100, 0, mind0)));
84 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_maxd0").c_str(), "Reco vertex max d0 [mm]", 100, 0, maxd0)));
85 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_ntrk").c_str(), "Reco vertex n tracks", 30, 0, 30)));
86
87 // tracks
88 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_qOverP").c_str(), "Reco track qOverP ", 100, 0, .01)));
89 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_theta").c_str(), "Reco track theta ", 64, 0, 3.2)));
90 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_E").c_str(), "Reco track E ", 100, 0, 100)));
91 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_M").c_str(), "Reco track M ", 100, 0, 10)));
92 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_Pt").c_str(), "Reco track Pt ", 100, 0, 100)));
93 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_Px").c_str(), "Reco track Px ", 100, 0, 100)));
94 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_Py").c_str(), "Reco track Py ", 100, 0, 100)));
95 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_Pz").c_str(), "Reco track Pz ", 100, 0, 100)));
96 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_Eta").c_str(), "Reco track Eta ", 100, -5, 5)));
97 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_Phi").c_str(), "Reco track Phi ", 63, -3.2, 3.2)));
98 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_D0").c_str(), "Reco track D0 ", 300, -maxTrackd0, maxTrackd0)));
99 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_Z0").c_str(), "Reco track Z0 ", 500, -maxTrackz0, maxTrackz0)));
100 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_errD0").c_str(), "Reco track errD0 ", 300, 0, maxErrd0)));
101 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_errZ0").c_str(), "Reco track errZ0 ", 500, 0, maxErrz0)));
102 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_Chi2").c_str(), "Reco track Chi2 ", 100, 0, 10)));
103 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_nDoF").c_str(), "Reco track nDoF ", 100, 0, 100)));
104 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_Trk_charge").c_str(), "Reco track charge ", 3, -1.5, 1.5)));
105
106 // truth matching -- don't book for non-matched vertices
107 if ( recoType != "All" and recoType != "Fake" ) {
108 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_positionRes_R").c_str(), "Position resolution for vertices matched to truth decays", 400, -maxResR, maxResR)));
109 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_positionRes_Z").c_str(), "Position resolution for vertices matched to truth decays", 400, -maxResZ, maxResZ)));
110 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_matchScore_weight").c_str(), "Vertex Match Score (weight)", 101, 0, 1.01)));
111 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_matchScore_pt").c_str(), "Vertex Match Score (pT)", 101, 0, 1.01)));
112 ANA_CHECK (book(TH1F(("RecoVertex/" + recoType + "_matchedTruthID").c_str(), "Vertex Truth Match ID", 100, 0, 100)));
113 }
114 }
115
116 // do reco vertices by SM origin if enabled -- NOTE these types are not exclusive (d decays will also be b decays in a cascade etc)
117 if (m_doSMOrigin) {
118 for(const auto& smOriginType : smOriginTypes) {
119 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_x").c_str(), "Reco vertex x [mm]", 1000, -maxX, maxX)));
120 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_y").c_str(), "Reco vertex y [mm]", 1000, -maxY, maxY)));
121 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_z").c_str(), "Reco vertex z [mm]", 1000, -maxZ, maxZ)));
122 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Lxy").c_str(), "Reco vertex L_{xy} [mm]", 500, 0, maxLxy)));
123 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_pT").c_str(), "Reco vertex p_{T} [GeV]", 100, 0, 100)));
124 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_eta").c_str(), "Reco vertex #eta", 100, -5, 5)));
125 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_phi").c_str(), "Reco vertex #phi", 100, -TMath::Pi(), TMath::Pi())));
126 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_mass").c_str(), "Reco vertex mass [GeV]", 500, 0, 100)));
127 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_mu").c_str(), "Reco vertex Red. Mass [GeV]", 500, 0, 100)));
128 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_chi2").c_str(), "Reco vertex recoChi2", 100, 0, 10)));
129 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_dir").c_str(), "Reco vertex recoDirection", 100, -1, 1)));
130 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_charge").c_str(), "Reco vertex recoCharge", 20, -10, 10)));
131 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_H").c_str(), "Reco vertex H [GeV]", 100, 0, 100)));
132 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_HT").c_str(), "Reco vertex Mass [GeV]", 100, 0, 100)));
133 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_minOpAng").c_str(), "Reco vertex minOpAng", 100, -1, 1)));
134 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_maxOpAng").c_str(), "Reco vertex maxOpAng", 100, -1, 1)));
135 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_maxdR").c_str(), "Reco vertex maxDR", 100, 0, 10)));
136 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_mind0").c_str(), "Reco vertex min d0 [mm]", 100, 0, mind0)));
137 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_maxd0").c_str(), "Reco vertex max d0 [mm]", 100, 0, maxd0)));
138 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_ntrk").c_str(), "Reco vertex n tracks", 30, 0, 30)));
139
140 // tracks
141 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_qOverP").c_str(), "Reco track qOverP ", 100, 0, .01)));
142 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_theta").c_str(), "Reco track theta ", 64, 0, 3.2)));
143 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_E").c_str(), "Reco track E ", 100, 0, 100)));
144 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_M").c_str(), "Reco track M ", 100, 0, 10)));
145 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_Pt").c_str(), "Reco track Pt ", 100, 0, 100)));
146 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_Px").c_str(), "Reco track Px ", 100, 0, 100)));
147 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_Py").c_str(), "Reco track Py ", 100, 0, 100)));
148 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_Pz").c_str(), "Reco track Pz ", 100, 0, 100)));
149 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_Eta").c_str(), "Reco track Eta ", 100, -5, 5)));
150 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_Phi").c_str(), "Reco track Phi ", 63, -3.2, 3.2)));
151 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_D0").c_str(), "Reco track D0 ", 300, -maxTrackd0, maxTrackd0)));
152 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_Z0").c_str(), "Reco track Z0 ", 500, -maxTrackz0, maxTrackz0)));
153 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_errD0").c_str(), "Reco track errD0 ", 300, 0, maxErrd0)));
154 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_errZ0").c_str(), "Reco track errZ0 ", 500, 0, maxErrz0)));
155 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_Chi2").c_str(), "Reco track Chi2 ", 100, 0, 10)));
156 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_nDoF").c_str(), "Reco track nDoF ", 100, 0, 100)));
157 ANA_CHECK (book(TH1F(("RecoVertex/" + smOriginType + "_Trk_charge").c_str(), "Reco track charge ", 3, -1.5, 1.5)));
158
159
160 }
161 }
162
163
164 for(const auto& truthType : truthTypes) {
165 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_x").c_str(), "Truth vertex x [mm]", 1000, -maxX, maxX)));
166 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_y").c_str(), "Truth vertex y [mm]", 500, -maxY, maxY)));
167 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_z").c_str(), "Truth vertex z [mm]", 500, -maxZ, maxZ)));
168 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_R").c_str(), "Truth vertex r [mm]", 6000, 0, maxR)));
169 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_distFromPV").c_str(), "Truth vertex distFromPV [mm]", 600, 0, maxR)));
170 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_Eta").c_str(), "Truth vertex Eta", 100, -5, 5)));
171 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_Phi").c_str(), "Truth vertex Phi", 64, -3.2, 3.2)));
172 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_Ntrk_out").c_str(), "Truth vertex n outgoing tracks", 100, 0, 100)));
173 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_Parent_E").c_str(), "Reco track E", 100, 0, 100)));
174 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_Parent_M").c_str(), "Reco track M", 500, 0, 500)));
175 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_Parent_Pt").c_str(), "Reco track Pt", 100, 0, 100)));
176 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_Parent_Eta").c_str(), "Reco track Eta", 100, -5, 5)));
177 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_Parent_Phi").c_str(), "Reco track Phi", 63, -3.2, 3.2)));
178 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_Parent_charge").c_str(), "Reco track charge", 3, -1, 1)));
179 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_ParentProdX").c_str(), "truthParentProd vertex x [mm]", 500, -500, 500)));
180 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_ParentProdY").c_str(), "truthParentProd vertex y [mm]", 500, -500, 500)));
181 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_ParentProdZ").c_str(), "truthParentProd vertex z [mm]", 500, -500, 500)));
182 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_ParentProdR").c_str(), "truthParentProd vertex r [mm]", 6000, 0, 600)));
183 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_ParentProddistFromPV").c_str(), "truthParentProd vertex distFromPV [mm]", 500, 0, 500)));
184 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_ParentProdEta").c_str(), "truthParentProd vertex Eta", 100, -5, 5)));
185 ANA_CHECK (book(TH1F(("TruthVertex/" + truthType + "_ParentProdPhi").c_str(), "truthParentProd vertex Phi", 64, -3.2, 3.2)));
186 }
187 // now add the efficiencies
188 // Define two different bin arrays - one for standard mode, one for MuSA
189 Double_t standard_bins[] = {0.0, 0.1, 0.2, 0.3, 0.4, 0.5, 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 15, 20, 25, 30, 35, 40, 50, 60, 70, 80, 90, 100, 125, 150, 200, 300, 500};
190 Double_t muSA_bins[] = {0.0, 1, 5, 10, 20, 50, 100, 200, 300, 500, 750, 1000, 1500, 2000, 2500, 3000, 4000, 5000, 6000, 7000, 8000};
191
192 // Determine which bins to use based on whether we're in MuSA mode
193 Double_t* bins = m_doMuSA ? muSA_bins : standard_bins;
194 size_t nbins = m_doMuSA ? sizeof(muSA_bins)/sizeof(muSA_bins[0])-1 : sizeof(standard_bins)/sizeof(standard_bins[0])-1;
195
196 ANA_CHECK (book(TEfficiency("Acceptance", "Acceptance", nbins, bins)));
197 ANA_CHECK (book(TEfficiency("eff_seed", "Seed efficiency", nbins, bins)));
198 ANA_CHECK (book(TEfficiency("eff_core", "Core efficiency", nbins, bins)));
199 ANA_CHECK (book(TEfficiency("eff_total", "Total efficiency", nbins, bins)));
200
201 }
202
203
204 return StatusCode::SUCCESS;
205 }
206
207 StatusCode SecVertexTruthMatchAlg::execute() {
208
209 //Retrieve the vertices:
210 SG::ReadHandle<xAOD::VertexContainer> recoVertexContainer(m_secVtxContainerKey);
211 SG::ReadHandle<xAOD::TruthVertexContainer> truthVertexContainer(m_truthVtxContainerKey);
212 SG::ReadHandle<xAOD::TrackParticleContainer> trackParticleContainer(m_trackParticleContainerKey);
213
214 std::vector<const xAOD::Vertex*> recoVerticesToMatch;
215 std::vector<const xAOD::TruthVertex*> truthVerticesToMatch;
216
217 for(const auto recoVertex : *recoVertexContainer) {
218 xAOD::VxType::VertexType vtxType = static_cast<xAOD::VxType::VertexType>( recoVertex->vertexType() );
219
220 if(vtxType != xAOD::VxType::SecVtx ){
221 ATH_MSG_DEBUG("Vertex not labeled as secondary");
222 continue;
223 }
224 recoVerticesToMatch.push_back(recoVertex);
225 }
226
227 for(const auto truthVertex : *truthVertexContainer) {
228 if(truthVertex->nIncomingParticles() != 1) {
229 continue;
230 }
231 const xAOD::TruthParticle* truthPart = truthVertex->incomingParticle(0);
232 if(not truthPart) {
233 continue;
234 }
235 if(std::find(m_targetPDGIDs.begin(), m_targetPDGIDs.end(), std::abs(truthPart->pdgId())) == m_targetPDGIDs.end()) {
236 continue;
237 }
238 if(truthVertex->nOutgoingParticles() < 2) {
239 continue;
240 }
241 truthVerticesToMatch.push_back(truthVertex);
242 }
243
244 //pass to the tool for decoration:
245 ATH_CHECK( m_matchTool->matchVertices( recoVerticesToMatch, truthVerticesToMatch, trackParticleContainer.cptr() ) );
246
247 if(m_writeHistograms) {
248 static const xAOD::Vertex::Decorator<int> matchTypeDecor("vertexMatchType");
249 for(const auto& secVtx : recoVerticesToMatch) {
250 int matchTypeBitset = matchTypeDecor(*secVtx);
251 hist("RecoVertex/matchType")->Fill(matchTypeBitset);
252
253 if(InDetSecVtxTruthMatchUtils::isMatched(matchTypeBitset)) {
254 fillRecoHistograms(secVtx, "Matched");
255 }
256 if(InDetSecVtxTruthMatchUtils::isMerged(matchTypeBitset)) {
257 fillRecoHistograms(secVtx, "Merged");
258 }
259 if(InDetSecVtxTruthMatchUtils::isFake(matchTypeBitset)) {
260 fillRecoHistograms(secVtx, "Fake");
261 }
262 if(InDetSecVtxTruthMatchUtils::isSplit(matchTypeBitset)) {
263 fillRecoHistograms(secVtx, "Split");
264 }
265 if(InDetSecVtxTruthMatchUtils::isOther(matchTypeBitset)) {
266 fillRecoHistograms(secVtx, "Other");
267 }
268 fillRecoHistograms(secVtx, "All");
269 }
270
271 static const xAOD::Vertex::Decorator<int> truthTypeDecor("truthVertexMatchType");
272 for(const auto& truthVtx : truthVerticesToMatch) {
273 int truthTypeBitset = truthTypeDecor(*truthVtx);
274 if(InDetSecVtxTruthMatchUtils::isReconstructable(truthTypeBitset)) {
275 fillTruthHistograms(truthVtx, "Reconstructable");
276
277 // fill efficiencies
278 efficiency("Acceptance")->Fill(InDetSecVtxTruthMatchUtils::isAccepted(truthTypeBitset), truthVtx->perp());
279 efficiency("eff_total")->Fill(InDetSecVtxTruthMatchUtils::isReconstructed(truthTypeBitset), truthVtx->perp());
280 }
281 if(InDetSecVtxTruthMatchUtils::isAccepted(truthTypeBitset)) {
282 fillTruthHistograms(truthVtx, "Accepted");
283 efficiency("eff_seed")->Fill(InDetSecVtxTruthMatchUtils::isSeeded(truthTypeBitset), truthVtx->perp());
284 }
285 if(InDetSecVtxTruthMatchUtils::isSeeded(truthTypeBitset)) {
286 fillTruthHistograms(truthVtx, "Seeded");
287 efficiency("eff_core")->Fill(InDetSecVtxTruthMatchUtils::isReconstructed(truthTypeBitset), truthVtx->perp());
288 }
289 if(InDetSecVtxTruthMatchUtils::isReconstructed(truthTypeBitset)) {
290 fillTruthHistograms(truthVtx, "Reconstructed");
291 }
292 if(InDetSecVtxTruthMatchUtils::isReconstructedSplit(truthTypeBitset)) {
293 fillTruthHistograms(truthVtx, "ReconstructedSplit");
294 }
295 fillTruthHistograms(truthVtx, "Inclusive");
296
297 }
298 if (m_doSMOrigin) {
299 static const xAOD::Vertex::Decorator<int> originTypeDecor("vertexMatchOriginType");
300
301 static const std::map<InDetSecVtxTruthMatchUtils::VertexMatchOriginType, std::string> originTypeMap = {
302 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::FakeOrigin, "FakeOrigin"},
303 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::Pileup, "Pileup"},
304 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::KshortDecay, "KshortDecay"},
305 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::StrangeMesonDecay, "StrangeMesonDecay"},
306 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::LambdaDecay, "LambdaDecay"},
307 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::StrangeBaryonDecay, "StrangeBaryonDecay"},
308 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::TauDecay, "TauDecay"},
309 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::GammaConversion, "GammaConversion"},
310 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::OtherDecay, "OtherDecay"},
311 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::HadronicInteraction, "HadronicInteraction"},
312 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::OtherSecondary, "OtherSecondary"},
313 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::BHadronDecay, "BHadronDecay"},
314 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::DHadronDecay, "DHadronDecay"},
315 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::Fragmentation, "Fragmentation"},
316 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::OtherOrigin, "OtherOrigin"},
317 {InDetSecVtxTruthMatchUtils::VertexMatchOriginType::Signal, "Signal"},
318 };
319
320 for(const auto& secVtx : recoVerticesToMatch) {
321 int smOriginTypeBitset = originTypeDecor(*secVtx);
322 hist("RecoVertex/smOriginType")->Fill(smOriginTypeBitset);
323
324 for(const auto& entry : originTypeMap) {
325 if(InDetSecVtxTruthMatchUtils::isOriginType(smOriginTypeBitset, entry.first)) {
326 fillOriginHistograms(secVtx, entry.second);
327 }
328 }
329 }
330 }
331
332 }
333
334 return StatusCode::SUCCESS;
335
336 }
337 void SecVertexTruthMatchAlg::fillRecoHistograms(const xAOD::Vertex* secVtx, const std::string& matchType) {
338
339 // set of accessors for tracks and weights
340 xAOD::Vertex::ConstAccessor<xAOD::Vertex::TrackParticleLinks_t> trkAcc("trackParticleLinks");
341 xAOD::Vertex::ConstAccessor<std::vector<float> > weightAcc("trackWeights");
342
343 // set of decorators for truth matching info
344 const xAOD::Vertex::Decorator<std::vector<InDetSecVtxTruthMatchUtils::VertexTruthMatchInfo> > matchInfoDecor("truthVertexMatchingInfos");
345
346 TVector3 reco_pos(secVtx->x(), secVtx->y(), secVtx->z());
347 float Lxy = reco_pos.Perp();
348
349 size_t ntracks;
350 const xAOD::Vertex::TrackParticleLinks_t & trkParts = trkAcc( *secVtx );
351 ntracks = trkParts.size();
352
353 TLorentzVector sumP4(0,0,0,0);
354 double H = 0.0;
355 double HT = 0.0;
356 int charge = 0;
357 double minOpAng = -1.0* 1.e10;
358 double maxOpAng = 1.0* 1.e10;
359 double minD0 = 1.0* 1.e10;
360 double maxD0 = 0.0;
361 double maxDR = 0.0;
362
363 xAOD::TrackParticle::ConstAccessor< std::vector< float > > accCovMatrixDiag( "definingParametersCovMatrixDiag" );
364
365 ATH_MSG_DEBUG("Loop over tracks");
366 for(size_t t = 0; t < ntracks; t++){
367 if(!trkParts[t].isValid()){
368 ATH_MSG_DEBUG("Track " << t << " is bad!");
369 continue;
370 }
371 const xAOD::TrackParticle & trk = **trkParts[t];
372
373 double trk_d0 = std::abs(trk.definingParameters()[0]);
374 double trk_z0 = std::abs(trk.definingParameters()[1]);
375
376 if(trk_d0 < minD0){ minD0 = trk_d0; }
377 if(trk_d0 > maxD0){ maxD0 = trk_d0; }
378
379 TLorentzVector vv;
380 // TODO: use values computed w.r.t SV
381 vv.SetPtEtaPhiM(trk.pt(),trk.eta(), trk.phi0(), trk.m());
382 sumP4 += vv;
383 H += vv.Vect().Mag();
384 HT += vv.Pt();
385
386 TLorentzVector v_minus_iv(0,0,0,0);
387 for(size_t j = 0; j < ntracks; j++){
388 if (j == t){ continue; }
389 if(!trkParts[j].isValid()){
390 ATH_MSG_DEBUG("Track " << j << " is bad!");
391 continue;
392 }
393
394 const xAOD::TrackParticle & trk_2 = **trkParts[j];
395
396 TLorentzVector tmp;
397 // TODO: use values computed w.r.t. SV
398 tmp.SetPtEtaPhiM(trk_2.pt(),trk_2.eta(), trk_2.phi0(), trk_2.m());
399 v_minus_iv += tmp;
400
401 if( j > t ) {
402 double tm = vv * tmp / ( vv.Mag() * tmp.Mag() );
403 if( minOpAng < tm ) minOpAng = tm;
404 if( maxOpAng > tm ) maxOpAng = tm;
405 }
406 }
407 double DR = vv.DeltaR(v_minus_iv);
408 if( DR > maxDR ){ maxDR = DR;}
409
410 charge += trk.charge();
411
412 xAOD::TrackParticle::ConstAccessor<float> Trk_Chi2("chiSquared");
413 xAOD::TrackParticle::ConstAccessor<float> Trk_nDoF("numberDoF");
414
415 if ( Trk_Chi2.isAvailable(trk) && Trk_Chi2(trk) && Trk_nDoF.isAvailable(trk) && Trk_nDoF(trk) ) {
416 hist("RecoVertex/" + matchType + "_Trk_Chi2")->Fill(Trk_Chi2(trk) / Trk_nDoF(trk));
417 hist("RecoVertex/" + matchType + "_Trk_nDoF")->Fill(Trk_nDoF(trk));
418 }
419 hist("RecoVertex/" + matchType + "_Trk_D0")->Fill(trk_d0);
420 hist("RecoVertex/" + matchType + "_Trk_Z0")->Fill(trk_z0);
421 hist("RecoVertex/" + matchType + "_Trk_theta")->Fill(trk.definingParameters()[3]);
422 hist("RecoVertex/" + matchType + "_Trk_qOverP")->Fill(trk.definingParameters()[4]);
423 hist("RecoVertex/" + matchType + "_Trk_Eta")->Fill(trk.eta());
424 hist("RecoVertex/" + matchType + "_Trk_Phi")->Fill(trk.phi0());
425 hist("RecoVertex/" + matchType + "_Trk_E")->Fill(trk.e() / GeV);
426 hist("RecoVertex/" + matchType + "_Trk_M")->Fill(trk.m() / GeV);
427 hist("RecoVertex/" + matchType + "_Trk_Pt")->Fill(trk.pt() / GeV);
428 hist("RecoVertex/" + matchType + "_Trk_Px")->Fill(trk.p4().Px() / GeV);
429 hist("RecoVertex/" + matchType + "_Trk_Py")->Fill(trk.p4().Py() / GeV);
430 hist("RecoVertex/" + matchType + "_Trk_Pz")->Fill(trk.p4().Pz() / GeV);
431 hist("RecoVertex/" + matchType + "_Trk_charge")->Fill(trk.charge());
432 hist("RecoVertex/" + matchType + "_Trk_errD0")->Fill(trk.definingParametersCovMatrix()(0,0));
433 hist("RecoVertex/" + matchType + "_Trk_errZ0")->Fill(trk.definingParametersCovMatrix()(1,1));
434
435 } // end loop over tracks
436
437 const double dir = sumP4.Vect().Dot( reco_pos ) / sumP4.Vect().Mag() / reco_pos.Mag();
438
439 xAOD::Vertex::ConstAccessor<float> Chi2("chiSquared");
440 xAOD::Vertex::ConstAccessor<float> nDoF("numberDoF");
441
442 hist("RecoVertex/" + matchType + "_x")->Fill(secVtx->x());
443 hist("RecoVertex/" + matchType + "_y")->Fill(secVtx->y());
444 hist("RecoVertex/" + matchType + "_z")->Fill(secVtx->z());
445 hist("RecoVertex/" + matchType + "_Lxy")->Fill(Lxy);
446 hist("RecoVertex/" + matchType + "_ntrk")->Fill(ntracks);
447 hist("RecoVertex/" + matchType + "_pT")->Fill(sumP4.Pt() / GeV);
448 hist("RecoVertex/" + matchType + "_eta")->Fill(sumP4.Eta());
449 hist("RecoVertex/" + matchType + "_phi")->Fill(sumP4.Phi());
450 hist("RecoVertex/" + matchType + "_mass")->Fill(sumP4.M() / GeV);
451 hist("RecoVertex/" + matchType + "_mu")->Fill(sumP4.M()/maxDR / GeV);
452 hist("RecoVertex/" + matchType + "_chi2")->Fill(Chi2(*secVtx)/nDoF(*secVtx));
453 hist("RecoVertex/" + matchType + "_dir")->Fill(dir);
454 hist("RecoVertex/" + matchType + "_charge")->Fill(charge);
455 hist("RecoVertex/" + matchType + "_H")->Fill(H / GeV);
456 hist("RecoVertex/" + matchType + "_HT")->Fill(HT / GeV);
457 hist("RecoVertex/" + matchType + "_minOpAng")->Fill(minOpAng);
458 hist("RecoVertex/" + matchType + "_maxOpAng")->Fill(maxOpAng);
459 hist("RecoVertex/" + matchType + "_mind0")->Fill(minD0);
460 hist("RecoVertex/" + matchType + "_maxd0")->Fill(maxD0);
461 hist("RecoVertex/" + matchType + "_maxdR")->Fill(maxDR);
462
463 std::vector<InDetSecVtxTruthMatchUtils::VertexTruthMatchInfo> truthmatchinfo;
464 truthmatchinfo = matchInfoDecor(*secVtx);
465
466 // This includes all matched vertices, including splits
467 if (matchType != "All" and matchType != "Fake") {
468 if(not truthmatchinfo.empty()){
469 float matchScore_weight = std::get<1>(truthmatchinfo.at(0));
470 float matchScore_pt = std::get<2>(truthmatchinfo.at(0));
471
472 ATH_MSG_DEBUG("Match Score and probability: " << matchScore_weight << " " << matchScore_pt/0.01);
473
474 const ElementLink<xAOD::TruthVertexContainer>& truthVertexLink = std::get<0>(truthmatchinfo.at(0));
475 const xAOD::TruthVertex& truthVtx = **truthVertexLink ;
476
477 hist("RecoVertex/" + matchType + "_positionRes_R")->Fill(Lxy - truthVtx.perp());
478 hist("RecoVertex/" + matchType + "_positionRes_Z")->Fill(secVtx->z() - truthVtx.z());
479 hist("RecoVertex/" + matchType + "_matchScore_weight")->Fill(matchScore_weight);
480 hist("RecoVertex/" + matchType + "_matchScore_pt")->Fill(matchScore_pt);
481 }
482 }
483 }
484
485 void SecVertexTruthMatchAlg::fillTruthHistograms(const xAOD::TruthVertex* truthVtx, const std::string& truthType) {
486
487 hist("TruthVertex/" + truthType + "_x")->Fill(truthVtx->x());
488 hist("TruthVertex/" + truthType + "_y")->Fill(truthVtx->y());
489 hist("TruthVertex/" + truthType + "_z")->Fill(truthVtx->z());
490 hist("TruthVertex/" + truthType + "_R")->Fill(truthVtx->perp());
491 hist("TruthVertex/" + truthType + "_Eta")->Fill(truthVtx->eta());
492 hist("TruthVertex/" + truthType + "_Phi")->Fill(truthVtx->phi());
493 hist("TruthVertex/" + truthType + "_Ntrk_out")->Fill(truthVtx->nOutgoingParticles());
494
495 ATH_MSG_DEBUG("Plotting truth parent");
496 const xAOD::TruthParticle& truthPart = *truthVtx->incomingParticle(0);
497
498 hist("TruthVertex/" + truthType + "_Parent_E")->Fill(truthPart.e() / GeV);
499 hist("TruthVertex/" + truthType + "_Parent_M")->Fill(truthPart.m() / GeV);
500 hist("TruthVertex/" + truthType + "_Parent_Pt")->Fill(truthPart.pt() / GeV);
501 hist("TruthVertex/" + truthType + "_Parent_Phi")->Fill(truthPart.phi());
502 hist("TruthVertex/" + truthType + "_Parent_Eta")->Fill(truthPart.eta());
503 hist("TruthVertex/" + truthType + "_Parent_charge")->Fill(truthPart.charge());
504
505 ATH_MSG_DEBUG("Plotting truth prod vtx");
506 if(truthPart.hasProdVtx()){
507 const xAOD::TruthVertex & vertex = *truthPart.prodVtx();
508
509 hist("TruthVertex/" + truthType + "_ParentProdX")->Fill(vertex.x());
510 hist("TruthVertex/" + truthType + "_ParentProdY")->Fill(vertex.y());
511 hist("TruthVertex/" + truthType + "_ParentProdZ")->Fill(vertex.z());
512 hist("TruthVertex/" + truthType + "_ParentProdR")->Fill(vertex.perp());
513 hist("TruthVertex/" + truthType + "_ParentProdEta")->Fill(vertex.eta());
514 hist("TruthVertex/" + truthType + "_ParentProdPhi")->Fill(vertex.phi());
515 }
516 // m_truthInclusive_r->Fill(truthVtx.perp());
517
518 // if(matchTypeDecor(truthVtx) >= RECONSTRUCTABLE){
519 // m_truthReconstructable_r->Fill(truthVtx.perp());
520 // }
521 // if(matchTypeDecor(truthVtx) >= ACCEPTED){
522 // m_truthAccepted_r->Fill(truthVtx.perp());
523 // }
524 // if(matchTypeDecor(truthVtx) >= SEEDED){
525 // m_truthSeeded_r->Fill(truthVtx.perp());
526 // }
527 // if(matchTypeDecor(truthVtx) >= RECONSTRUCTED){
528 // m_truthReconstructed_r->Fill(truthVtx.perp());
529 // }
530 // if(matchTypeDecor(truthVtx) >= RECONSTRUCTEDSPLIT){
531 // m_truthSplit_r->Fill(truthVtx.perp());
532 // }
533 //
534 }
535
536 void SecVertexTruthMatchAlg::fillOriginHistograms(const xAOD::Vertex* secVtx, const std::string& originType) {
537
538 // set of accessors for tracks and weights
539 xAOD::Vertex::ConstAccessor<xAOD::Vertex::TrackParticleLinks_t> trkAcc("trackParticleLinks");
540 xAOD::Vertex::ConstAccessor<std::vector<float> > weightAcc("trackWeights");
541
542 TVector3 reco_pos(secVtx->x(), secVtx->y(), secVtx->z());
543 float Lxy = reco_pos.Perp();
544
545 size_t ntracks;
546 const xAOD::Vertex::TrackParticleLinks_t & trkParts = trkAcc( *secVtx );
547 ntracks = trkParts.size();
548
549 TLorentzVector sumP4(0,0,0,0);
550 double H = 0.0;
551 double HT = 0.0;
552 int charge = 0;
553 double minOpAng = -1.0* 1.e10;
554 double maxOpAng = 1.0* 1.e10;
555 double minD0 = 1.0* 1.e10;
556 double maxD0 = 0.0;
557 double maxDR = 0.0;
558
559 xAOD::TrackParticle::ConstAccessor< std::vector< float > > accCovMatrixDiag( "definingParametersCovMatrixDiag" );
560
561 // Loop over tracks to calculate derived quantities
562 for(size_t t = 0; t < ntracks; t++){
563 if(!trkParts[t].isValid()){
564 continue;
565 }
566 const xAOD::TrackParticle & trk = **trkParts[t];
567
568 double trk_d0 = std::abs(trk.definingParameters()[0]);
569 double trk_z0 = std::abs(trk.definingParameters()[1]);
570
571 if(trk_d0 < minD0){ minD0 = trk_d0; }
572 if(trk_d0 > maxD0){ maxD0 = trk_d0; }
573
574 TLorentzVector vv;
575 vv.SetPtEtaPhiM(trk.pt(), trk.eta(), trk.phi0(), trk.m());
576 sumP4 += vv;
577 H += vv.Vect().Mag();
578 HT += vv.Pt();
579
580 TLorentzVector v_minus_iv(0,0,0,0);
581 for(size_t j = 0; j < ntracks; j++){
582 if (j == t){ continue; }
583 if(!trkParts[j].isValid()){
584 continue;
585 }
586
587 const xAOD::TrackParticle & trk_2 = **trkParts[j];
588
589 TLorentzVector tmp;
590 tmp.SetPtEtaPhiM(trk_2.pt(), trk_2.eta(), trk_2.phi0(), trk_2.m());
591 v_minus_iv += tmp;
592
593 if( j > t ) {
594 double tm = vv * tmp / ( vv.Mag() * tmp.Mag() );
595 if( minOpAng < tm ) minOpAng = tm;
596 if( maxOpAng > tm ) maxOpAng = tm;
597 }
598 }
599 double DR = vv.DeltaR(v_minus_iv);
600 if( DR > maxDR ){ maxDR = DR;}
601
602 charge += trk.charge();
603
604 // Fill track-level histograms
605 xAOD::TrackParticle::ConstAccessor<float> Trk_Chi2("chiSquared");
606 xAOD::TrackParticle::ConstAccessor<float> Trk_nDoF("numberDoF");
607
608 if ( Trk_Chi2.isAvailable(trk) && Trk_Chi2(trk) && Trk_nDoF.isAvailable(trk) && Trk_nDoF(trk) ) {
609 hist("RecoVertex/" + originType + "_Trk_Chi2")->Fill(Trk_Chi2(trk) / Trk_nDoF(trk));
610 hist("RecoVertex/" + originType + "_Trk_nDoF")->Fill(Trk_nDoF(trk));
611 }
612 hist("RecoVertex/" + originType + "_Trk_D0")->Fill(trk_d0);
613 hist("RecoVertex/" + originType + "_Trk_Z0")->Fill(trk_z0);
614 hist("RecoVertex/" + originType + "_Trk_theta")->Fill(trk.definingParameters()[3]);
615 hist("RecoVertex/" + originType + "_Trk_qOverP")->Fill(trk.definingParameters()[4]);
616 hist("RecoVertex/" + originType + "_Trk_Eta")->Fill(trk.eta());
617 hist("RecoVertex/" + originType + "_Trk_Phi")->Fill(trk.phi0());
618 hist("RecoVertex/" + originType + "_Trk_E")->Fill(trk.e() / GeV);
619 hist("RecoVertex/" + originType + "_Trk_M")->Fill(trk.m() / GeV);
620 hist("RecoVertex/" + originType + "_Trk_Pt")->Fill(trk.pt() / GeV);
621 hist("RecoVertex/" + originType + "_Trk_Px")->Fill(trk.p4().Px() / GeV);
622 hist("RecoVertex/" + originType + "_Trk_Py")->Fill(trk.p4().Py() / GeV);
623 hist("RecoVertex/" + originType + "_Trk_Pz")->Fill(trk.p4().Pz() / GeV);
624 hist("RecoVertex/" + originType + "_Trk_charge")->Fill(trk.charge());
625 hist("RecoVertex/" + originType + "_Trk_errD0")->Fill(trk.definingParametersCovMatrix()(0,0));
626 hist("RecoVertex/" + originType + "_Trk_errZ0")->Fill(trk.definingParametersCovMatrix()(1,1));
627 }
628
629 const double dir = sumP4.Vect().Dot( reco_pos ) / sumP4.Vect().Mag() / reco_pos.Mag();
630
631 xAOD::Vertex::ConstAccessor<float> Chi2("chiSquared");
632 xAOD::Vertex::ConstAccessor<float> nDoF("numberDoF");
633
634 // Fill vertex-level histograms
635 hist("RecoVertex/" + originType + "_x")->Fill(secVtx->x());
636 hist("RecoVertex/" + originType + "_y")->Fill(secVtx->y());
637 hist("RecoVertex/" + originType + "_z")->Fill(secVtx->z());
638 hist("RecoVertex/" + originType + "_Lxy")->Fill(Lxy);
639 hist("RecoVertex/" + originType + "_ntrk")->Fill(ntracks);
640 hist("RecoVertex/" + originType + "_pT")->Fill(sumP4.Pt() / GeV);
641 hist("RecoVertex/" + originType + "_eta")->Fill(sumP4.Eta());
642 hist("RecoVertex/" + originType + "_phi")->Fill(sumP4.Phi());
643 hist("RecoVertex/" + originType + "_mass")->Fill(sumP4.M() / GeV);
644 hist("RecoVertex/" + originType + "_mu")->Fill(sumP4.M()/maxDR / GeV);
645 hist("RecoVertex/" + originType + "_chi2")->Fill(Chi2(*secVtx)/nDoF(*secVtx));
646 hist("RecoVertex/" + originType + "_dir")->Fill(dir);
647 hist("RecoVertex/" + originType + "_charge")->Fill(charge);
648 hist("RecoVertex/" + originType + "_H")->Fill(H / GeV);
649 hist("RecoVertex/" + originType + "_HT")->Fill(HT / GeV);
650 hist("RecoVertex/" + originType + "_minOpAng")->Fill(minOpAng);
651 hist("RecoVertex/" + originType + "_maxOpAng")->Fill(maxOpAng);
652 hist("RecoVertex/" + originType + "_mind0")->Fill(minD0);
653 hist("RecoVertex/" + originType + "_maxd0")->Fill(maxD0);
654 hist("RecoVertex/" + originType + "_maxdR")->Fill(maxDR);
655 }
656
657} // namespace CP
ATH_CHECK
#define ATH_CHECK
Evaluate an expression and check for errors.
Definition AthCheckMacros.h:40
ATH_MSG_DEBUG
#define ATH_MSG_DEBUG(x)
Definition AthMsgStreamMacros.h:29
isValid
bool isValid(const T &p)
Av: we implement here an ATLAS-sepcific convention: all particles which are 99xxxxx are fine.
Definition AtlasPID.h:878
ANA_CHECK
#define ANA_CHECK(EXP)
check whether the given expression was successful
Definition Control/AthToolSupport/AsgMessaging/AsgMessaging/MessageCheck.h:324
H
#define H(x, y, z)
Definition MD5.cxx:114
bins
static const std::vector< std::string > bins
Definition MuonTriggerSFFilesTest.cxx:22
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
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
AthHistogramming::efficiency
TEfficiency * efficiency(const std::string &effName, const std::string &tDir="", const std::string &stream="")
Simplify the retrieval of registered TEfficiency.
Definition AthHistogramming.cxx:250
CP::SecVertexTruthMatchAlg::fillOriginHistograms
void fillOriginHistograms(const xAOD::Vertex *secVtx, const std::string &originType)
Definition SecVertexTruthMatchAlg.cxx:536
CP::SecVertexTruthMatchAlg::m_doMuSA
Gaudi::Property< bool > m_doMuSA
Definition SecVertexTruthMatchAlg.h:55
CP::SecVertexTruthMatchAlg::m_writeHistograms
Gaudi::Property< bool > m_writeHistograms
Definition SecVertexTruthMatchAlg.h:52
CP::SecVertexTruthMatchAlg::m_targetPDGIDs
Gaudi::Property< std::vector< int > > m_targetPDGIDs
Definition SecVertexTruthMatchAlg.h:50
CP::SecVertexTruthMatchAlg::m_doSMOrigin
Gaudi::Property< bool > m_doSMOrigin
Definition SecVertexTruthMatchAlg.h:56
CP::SecVertexTruthMatchAlg::initialize
virtual StatusCode initialize() override
Definition SecVertexTruthMatchAlg.cxx:20
CP::SecVertexTruthMatchAlg::fillRecoHistograms
void fillRecoHistograms(const xAOD::Vertex *secVtx, const std::string &matchType)
Definition SecVertexTruthMatchAlg.cxx:337
CP::SecVertexTruthMatchAlg::execute
virtual StatusCode execute() override
Definition SecVertexTruthMatchAlg.cxx:207
CP::SecVertexTruthMatchAlg::m_trackParticleContainerKey
SG::ReadHandleKey< xAOD::TrackParticleContainer > m_trackParticleContainerKey
Definition SecVertexTruthMatchAlg.h:47
CP::SecVertexTruthMatchAlg::m_secVtxContainerKey
SG::ReadHandleKey< xAOD::VertexContainer > m_secVtxContainerKey
Definition SecVertexTruthMatchAlg.h:41
CP::SecVertexTruthMatchAlg::m_matchTool
ToolHandle< IInDetSecVtxTruthMatchTool > m_matchTool
Definition SecVertexTruthMatchAlg.h:54
CP::SecVertexTruthMatchAlg::SecVertexTruthMatchAlg
SecVertexTruthMatchAlg(const std::string &name, ISvcLocator *svcLoc)
Regular Algorithm constructor.
Definition SecVertexTruthMatchAlg.cxx:17
CP::SecVertexTruthMatchAlg::fillTruthHistograms
void fillTruthHistograms(const xAOD::TruthVertex *truthVtx, const std::string &truthType)
Definition SecVertexTruthMatchAlg.cxx:485
CP::SecVertexTruthMatchAlg::m_truthVtxContainerKey
SG::ReadHandleKey< xAOD::TruthVertexContainer > m_truthVtxContainerKey
Definition SecVertexTruthMatchAlg.h:44
EL::AnaAlgorithm::AnaAlgorithm
AnaAlgorithm(const std::string &name, ISvcLocator *pSvcLocator)
constructor with parameters
Definition AnaAlgorithm.cxx:40
SG::ConstAccessor::isAvailable
bool isAvailable(const ELT &e) const
Test to see if this variable exists in the store.
SG::ConstAuxElement::Decorator
SG::Decorator< T, ALLOC > Decorator
class to provide type-safe access to aux data.
Definition AuxElement.h:135
SG::ConstAuxElement::ConstAccessor
SG::ConstAccessor< T, ALLOC > ConstAccessor
Helper class to provide type-safe access to aux data.
Definition AuxElement.h:131
SG::ReadHandle::cptr
const_pointer_type cptr()
Dereference the pointer.
xAOD::TrackParticle_v1::definingParameters
DefiningParameters_t definingParameters() const
Returns a SVector of the Perigee track parameters.
Definition TrackParticle_v1.cxx:163
xAOD::TrackParticle_v1::m
virtual double m() const override final
The invariant mass of the particle..
Definition TrackParticle_v1.cxx:85
xAOD::TrackParticle_v1::p4
virtual FourMom_t p4() const override final
The full 4-momentum of the particle.
Definition TrackParticle_v1.cxx:122
xAOD::TrackParticle_v1::definingParametersCovMatrix
const ParametersCovMatrix_t definingParametersCovMatrix() const
Returns the 5x5 symmetric matrix containing the defining parameters covariance matrix.
Definition TrackParticle_v1.cxx:238
xAOD::TrackParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition TrackParticle_v1.cxx:75
xAOD::TrackParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition TrackParticle_v1.cxx:79
xAOD::TrackParticle_v1::charge
float charge() const
Returns the charge.
Definition TrackParticle_v1.cxx:143
xAOD::TrackParticle_v1::phi0
float phi0() const
Returns the parameter, which has range to .
Definition TrackParticle_v1.cxx:150
xAOD::TrackParticle_v1::e
virtual double e() const override final
The total energy of the particle.
Definition TrackParticle_v1.cxx:111
xAOD::TruthParticle_v1::m
virtual double m() const override final
The mass of the particle.
xAOD::TruthParticle_v1::pdgId
int pdgId() const
PDG ID code.
xAOD::TruthParticle_v1::hasProdVtx
bool hasProdVtx() const
Check for a production vertex on this particle.
Definition TruthParticle_v1.cxx:69
xAOD::TruthParticle_v1::e
virtual double e() const override final
The total energy of the particle.
xAOD::TruthParticle_v1::pt
virtual double pt() const override final
The transverse momentum ( ) of the particle.
Definition TruthParticle_v1.cxx:161
xAOD::TruthParticle_v1::prodVtx
const TruthVertex_v1 * prodVtx() const
The production vertex of this particle.
Definition TruthParticle_v1.cxx:75
xAOD::TruthParticle_v1::charge
double charge() const
Physical charge.
xAOD::TruthParticle_v1::eta
virtual double eta() const override final
The pseudorapidity ( ) of the particle.
Definition TruthParticle_v1.cxx:169
xAOD::TruthParticle_v1::phi
virtual double phi() const override final
The azimuthal angle ( ) of the particle.
Definition TruthParticle_v1.cxx:176
xAOD::TruthVertex_v1::z
float z() const
Vertex longitudinal distance along the beam line form the origin.
xAOD::TruthVertex_v1::eta
float eta() const
Vertex pseudorapidity.
Definition TruthVertex_v1.cxx:171
xAOD::TruthVertex_v1::y
float y() const
Vertex y displacement.
xAOD::TruthVertex_v1::phi
float phi() const
Vertex azimuthal angle.
Definition TruthVertex_v1.cxx:177
xAOD::TruthVertex_v1::incomingParticle
const TruthParticle_v1 * incomingParticle(size_t index) const
Get one of the incoming particles.
Definition TruthVertex_v1.cxx:70
xAOD::TruthVertex_v1::perp
float perp() const
Vertex transverse distance from the beam line.
Definition TruthVertex_v1.cxx:164
xAOD::TruthVertex_v1::nOutgoingParticles
size_t nOutgoingParticles() const
Get the number of outgoing particles.
xAOD::TruthVertex_v1::x
float x() const
Vertex x displacement.
xAOD::Vertex_v1::z
float z() const
Returns the z position.
xAOD::Vertex_v1::TrackParticleLinks_t
std::vector< ElementLink< xAOD::TrackParticleContainer > > TrackParticleLinks_t
Type for the associated track particles.
Definition Vertex_v1.h:128
xAOD::Vertex_v1::y
float y() const
Returns the y position.
xAOD::Vertex_v1::x
float x() const
Returns the x position.
CP
Select isolated Photons, Electrons and Muons.
Definition Control/xAODRootAccess/xAODRootAccess/TEvent.h:27
EL
This module defines the arguments passed from the BATCH driver to the BATCH worker.
Definition AsgComponentFactories.h:16
InDetSecVtxTruthMatchUtils::isOriginType
bool isOriginType(int matchInfo, VertexMatchOriginType type)
Definition InDetSecVtxTruthMatchTool.h:117
xAOD::VxType::VertexType
VertexType
Vertex types.
Definition TrackingPrimitives.h:570
xAOD::VxType::SecVtx
@ SecVtx
Secondary vertex.
Definition TrackingPrimitives.h:573
xAOD::TruthVertex
TruthVertex_v1 TruthVertex
Typedef to implementation.
Definition TruthVertex.h:15
xAOD::TrackParticle
TrackParticle_v1 TrackParticle
Reference the current persistent version:
Definition Event/xAOD/xAODTracking/xAODTracking/TrackParticle.h:13
xAOD::Vertex
Vertex_v1 Vertex
Define the latest version of the vertex class.
Definition Event/xAOD/xAODTracking/xAODTracking/Vertex.h:16
xAOD::TruthParticle
TruthParticle_v1 TruthParticle
Typedef to implementation.
Definition Event/xAOD/xAODTruth/xAODTruth/TruthParticle.h:15
xAOD::accCovMatrixDiag
static const SG::AuxElement::Accessor< std::vector< float > > accCovMatrixDiag("definingParametersCovMatrixDiag")
Generated on for ATLAS Offline Software by doxygen 1.14.0