Belle II Software  release-08-01-10
EffPlotsModule.cc
1 /**************************************************************************
2  * basf2 (Belle II Analysis Software Framework) *
3  * Author: The Belle II Collaboration *
4  * *
5  * See git log for contributors and copyright holders. *
6  * This file is licensed under LGPL-3.0, see LICENSE.md. *
7  **************************************************************************/
8 
9 #include <tracking/modules/trackingPerformanceEvaluation/EffPlotsModule.h>
10 
11 #include <framework/datastore/RelationVector.h>
12 
13 #include <framework/geometry/BFieldManager.h>
14 
15 #include <framework/gearbox/Const.h>
16 
17 #include <tracking/dataobjects/MCParticleInfo.h>
18 #include <tracking/dataobjects/RecoTrack.h>
19 #include <tracking/dataobjects/V0ValidationVertex.h>
20 
21 #include <root/TObject.h>
22 
23 #include <vector>
24 
25 using namespace Belle2;
26 
27 //-----------------------------------------------------------------
28 // Register the Module
29 //-----------------------------------------------------------------
30 
31 REG_MODULE(EffPlots);
32 
33 EffPlotsModule::EffPlotsModule() :
34  Module()
35 {
36 
37  setDescription("This module makes some plots related to V0 and saves them into a root file. For the efficiency plots: _noGeoAcc -> normalized to MCParticles; _withGeoAcc -> normalized to RecoTracks.");
38 
39  addParam("outputFileName", m_rootFileName, "Name of output root file.",
40  std::string("EffPlots_output.root"));
41  addParam("V0sName", m_V0sName, "Name of V0 collection.", std::string("V0ValidationVertexs"));
42  addParam("MCParticlesName", m_MCParticlesName, "Name of MC Particle collection.", std::string(""));
43  addParam("TrackColName", m_TrackColName,
44  "Belle2::Track collection name (input). Note that the V0s use "
45  "pointers indices into these arrays, so all hell may break loose "
46  "if you change this.", std::string(""));
47  addParam("RecoTracksName", m_RecoTracksName, "Name of RecoTrack collection.", std::string("RecoTracks"));
48  addParam("MCRecoTracksName", m_MCRecoTracksName, "Name of MCRecoTrack collection.", std::string("MCRecoTracks"));
49 
50  addParam("V0sType", m_V0sType, "Type of V0 to perform plots. Default is Lambda0, alternatively Ks", std::string("Lambda0"));
51  addParam("AllHistograms", m_allHistograms, "Create output for all histograms, not only efficiencies.", bool(false));
52  addParam("GeometricalAccettance", m_geometricalAccettance, "Create output for geometrical accettance.", bool(false));
53 }
54 
55 EffPlotsModule::~EffPlotsModule()
56 {
57 
58 }
59 
60 void EffPlotsModule::initialize()
61 {
62  m_MCParticles.isRequired(m_MCParticlesName);
63  StoreArray<V0ValidationVertex> v0ValidationVertices;
64  v0ValidationVertices.isRequired(m_V0sName);
65 
66  StoreArray<TrackFitResult> trackFitResults(m_TFRColName);
67  trackFitResults.isRequired();
68 
69  StoreArray<Track> tracks(m_TrackColName);
70  tracks.isRequired();
71 
72  StoreArray<RecoTrack> recoTracks;
73  recoTracks.isRequired(m_MCRecoTracksName);
74 
75  //create list of histograms to be saved in the rootfile
76  m_histoList = new TList;
77  m_histoList_MCParticles = new TList;
78  m_histoList_RecoTracks = new TList;
79  m_histoList_Tracks = new TList;
80  m_histoList_Efficiencies = new TList;
81  m_histoList_GA = new TList;
82  m_histoList_check = new TList;
83 
84  //set the ROOT File
85  m_rootFilePtr = new TFile(m_rootFileName.c_str(), "RECREATE");
86 
87  Double_t bins_pt_new[25 + 1] = {0., 0.05, 0.1, 0.15, 0.2, 0.25, 0.3, 0.4, 0.5, 0.6, 0.7, 0.8, 0.9, 1, 1.2, 1.4, 1.6, 1.8, 2, 2.2, 2.4, 2.6, 2.8, 3., 3.2, 3.4};
88 
89  Double_t bins_theta[10 + 1];
90  Double_t width_theta = TMath::Pi() / 10;
91  for (unsigned int bin = 0; bin < 10 + 1; bin++)
92  bins_theta[bin] = bin * width_theta;
93 
94  Double_t bins_phi[14 + 1];
95  Double_t width_phi = 2 * TMath::Pi() / 14;
96  for (unsigned int bin = 0; bin < 14 + 1; bin++)
97  bins_phi[bin] = - TMath::Pi() + bin * width_phi;
98 
99  Double_t bins_costheta[20 + 1];
100  Double_t width_cosTheta = 2. / 20;
101  for (unsigned int bin1 = 0; bin1 < 20 + 1; bin1++)
102  bins_costheta[bin1] = - 1 + bin1 * width_cosTheta;
103 
104  //create histograms
105 
106  //------------------------------------------------------------------//
107  // MC PARTICLES //
108  //------------------------------------------------------------------//
109 
110  // MC dau0
111  m_h1_MC_dau0_d0 = createHistogram1D("h1MCdau0D0", "d0 dau_{0}", 100, -10, 10, "d0_{dau_{0}}", m_histoList_MCParticles);
112  m_h1_MC_dau0_z0 = createHistogram1D("h1MCdau0Z0", "z0 dau_{0}", 100, -10, 10, "z0_{dau_{0}}", m_histoList_MCParticles);
113  m_h1_MC_dau0_RMother = createHistogram1D("h1MCdau0RMother", "dau_{0}, R mother", 200, 0, 20, "R mother", m_histoList_MCParticles);
114  m_h3_MC_dau0 = createHistogram3D("h3MCdau0", "entry per MC dau_{0}",
115  25, bins_pt_new, "p_{t} (GeV/c)",
116  10, bins_theta, "#theta",
117  14, bins_phi, "#phi", m_histoList_MCParticles);
118  m_h1_MC_dau0_pt = createHistogram1D("h1MCdau0Pt", "dau_{0}, p_{T}", 25, bins_pt_new, "p_{T} (GeV/c)", m_histoList_MCParticles);
119  m_h1_MC_dau0_pz = createHistogram1D("h1MCdau0Pz", "dau_{0}, p_{z}", 25, bins_pt_new, "p_{z} (GeV/c)", m_histoList_MCParticles);
120  m_h1_MC_dau0_p = createHistogram1D("h1MCdau0P", "dau_{0}, p", 25, bins_pt_new, "p (GeV/c)", m_histoList_MCParticles);
121  m_h1_MC_dau0_theta = createHistogram1D("h1MCdau0Theta", "dau_{0}, #theta", 10, bins_theta, "#theta",
122  m_histoList_MCParticles);
123  m_h1_MC_dau0_costheta = createHistogram1D("h1MCdau0CosTheta", "dau_{0}, cos#theta", 20, bins_costheta, "cos#theta",
124  m_histoList_MCParticles);
125  m_h1_MC_dau0_Mother_cosAngle = createHistogram1D("h1MCdau0MothercosAngle", "cos#theta_{mother,dau_{0}}", 20,
126  bins_costheta, "cos#theta", m_histoList_MCParticles);
127  m_h1_MC_dau0_phi = createHistogram1D("h1MCdau0Phi", "dau_{0}, #phi", 14, bins_phi, "#phi", m_histoList_MCParticles);
128  m_h1_MC_dau0_phi_BW = createHistogram1D("h1MCdau0PhiBW", "dau_{0}, #phi", 14, bins_phi, "#phi", m_histoList_MCParticles);
129  m_h1_MC_dau0_phi_barrel = createHistogram1D("h1MCdau0Phibarrel", "dau_{0}, #phi", 14, bins_phi, "#phi",
130  m_histoList_MCParticles);
131  m_h1_MC_dau0_phi_FW = createHistogram1D("h1MCdau0PhiFW", "dau_{0}, #phi", 14, bins_phi, "#phi", m_histoList_MCParticles);
132 
133  m_h1_MC_dau0_phiMother_total = createHistogram1D("h1MCdau0phiMothertotal", "dau_{0}, #phi_{mother}", 14, bins_phi,
134  "#phi_{mother}", m_histoList_MCParticles);
135  m_h1_MC_dau0_phiMother_BW = createHistogram1D("h1MCdau0phiMotherBW", "dau_{0}, #phi_{mother}, BW", 14, bins_phi,
136  "#phi_{mother} BW", m_histoList_MCParticles);
137  m_h1_MC_dau0_phiMother_barrel = createHistogram1D("h1MCdau0phiMotherbarrel", "dau_{0}, #phi_{mother}, barrel", 14, bins_phi,
138  "#phi_{mother} barrel", m_histoList_MCParticles);
139  m_h1_MC_dau0_phiMother_FW = createHistogram1D("h1MCdau0phiMotherFW", "dau_{0}, #phi_{mother}, FW", 14, bins_phi,
140  "#phi_{mother} FW", m_histoList_MCParticles);
141 
142  m_h1_MC_dau0_thetaMother = createHistogram1D("h1MCdau0ThetaMother", "dau_{0}, #theta_{mother}", 10, bins_theta,
143  "#theta_{mother}", m_histoList_MCParticles);
144  m_h1_MC_dau0_ptMother = createHistogram1D("h1MCdau0PtMother", "dau_{0}, p_{T,mother}", 25, bins_pt_new, "p_{T,mother} (GeV/c)",
145  m_histoList_MCParticles);
146 
147  m_h2_MC_dau0_2D = createHistogram2D("h2MCdau0", "entry per MC dau_{0}",
148  10, bins_theta, "#theta",
149  25, bins_pt_new, "p_{t} (GeV/c)", m_histoList_MCParticles);
150  m_h2_MC_dau0_2D_BP = createHistogram2D("h2MCdau0BP", "entry per MC dau_{0}, beam pipe",
151  10, bins_theta, "#theta",
152  25, bins_pt_new, "p_{t} (GeV/c)", m_histoList_MCParticles);
153 
154  m_h2_MC_dau0_2DMother = createHistogram2D("h2MCdau0Mother", "entry mother per MC dau_{0}",
155  10, bins_theta, "#theta_{mother}",
156  25, bins_pt_new, "p_{T,mother} (GeV/c)", m_histoList_MCParticles);
157 
158  m_h2_MC_dau0_pVScostheta = createHistogram2D("h2MCdau0pVScostheta", "p_{CM} VS cos(#theta)_{CM}, dau_{0}",
159  20, -1., 1.,
160  "cos(#theta)_{CM}",
161  50, 0., 5.,
162  "p_{CM} [GeV]", m_histoList_MCParticles);
163 
164  m_h1_MC_dau0_PDG = createHistogram1D("h1MCdau0PDG", "PDG code, dau_{0}", 4600, -2300, 2300, "PDG", m_histoList_check);
165 
166  //MC dau1
167  m_h1_MC_dau1_d0 = createHistogram1D("h1MCdau1D0", "d0 dau_{1}", 100, -10, 10, "d0_{dau_{1}}", m_histoList_MCParticles);
168  m_h1_MC_dau1_z0 = createHistogram1D("h1MCdau1Z0", "z0 dau_{1}", 100, -10, 10, "z0_{dau_{1}}", m_histoList_MCParticles);
169  m_h1_MC_dau1_RMother = createHistogram1D("h1MCdau1RMother", "dau_{1}, R mother", 200, 0, 20, "R mother", m_histoList_MCParticles);
170  m_h3_MC_dau1 = createHistogram3D("h3MCdau1", "entry per MC dau_{1}",
171  25, bins_pt_new, "p_{t} (GeV/c)",
172  10, bins_theta, "#theta",
173  14, bins_phi, "#phi", m_histoList_MCParticles);
174  m_h1_MC_dau1_pt = createHistogram1D("h1MCdau1Pt", "dau_{1}, p_{T}", 25, bins_pt_new, "p_{T} (GeV/c)", m_histoList_MCParticles);
175  m_h1_MC_dau1_pz = createHistogram1D("h1MCdau1Pz", "dau_{1}, p_{z}", 25, bins_pt_new, "p_{z} (GeV/c)", m_histoList_MCParticles);
176  m_h1_MC_dau1_p = createHistogram1D("h1MCdau1P", "dau_{1}, ", 25, bins_pt_new, "p (GeV/c)", m_histoList_MCParticles);
177  m_h1_MC_dau1_theta = createHistogram1D("h1MCdau1Theta", "dau_{1}, #theta", 10, bins_theta, "#theta",
178  m_histoList_MCParticles);
179  m_h1_MC_dau1_costheta = createHistogram1D("h1MCdau1CosTheta", "dau_{1}, cos#theta", 20, bins_costheta, "cos#theta",
180  m_histoList_MCParticles);
181  m_h1_MC_dau1_phi = createHistogram1D("h1MCdau1Phi", "dau_{1}, #phi", 14, bins_phi, "#phi", m_histoList_MCParticles);
182  m_h1_MC_dau1_phi_BW = createHistogram1D("h1MCdau1PhiBW", "dau_{1}, #phi", 14, bins_phi, "#phi", m_histoList_MCParticles);
183  m_h1_MC_dau1_phi_barrel = createHistogram1D("h1MCdau1Phibarrel", "dau_{1}, #phi", 14, bins_phi, "#phi",
184  m_histoList_MCParticles);
185  m_h1_MC_dau1_phi_FW = createHistogram1D("h1MCdau1PhiFW", "dau_{1}, #phi", 14, bins_phi, "#phi", m_histoList_MCParticles);
186 
187  m_h1_MC_dau1_Mother_cosAngle = createHistogram1D("h1MCdau1MothercosAngle", "cos#theta_{mother,p}", 20, bins_costheta,
188  "cos#theta", m_histoList_MCParticles);
189 
190  m_h1_MC_dau1_phiMother_total = createHistogram1D("h1MCdau1phiMothertotal", "dau_{1}, #phi_{mother}", 14, bins_phi,
191  "#phi_{mother}", m_histoList_MCParticles);
192  m_h1_MC_dau1_phiMother_BW = createHistogram1D("h1MCdau1phiMotherBW", "dau_{1}, #phi_{mother}, BW", 14, bins_phi,
193  "#phi_{mother} BW", m_histoList_MCParticles);
194  m_h1_MC_dau1_phiMother_barrel = createHistogram1D("h1MCdau1phiMotherbarrel", "dau_{1}, #phi_{mother}, barrel", 14, bins_phi,
195  "#phi_{mother} barrel", m_histoList_MCParticles);
196  m_h1_MC_dau1_phiMother_FW = createHistogram1D("h1MCdau1phiMotherFW", "dau_{1}, #phi_{mother}, FW", 14, bins_phi,
197  "#phi_{mother} FW", m_histoList_MCParticles);
198 
199  m_h1_MC_dau1_thetaMother = createHistogram1D("h1MCdau1ThetaMother", "dau_{1}, #theta_{mother}", 10, bins_theta,
200  "#theta_{mother}", m_histoList_MCParticles);
201  m_h1_MC_dau1_ptMother = createHistogram1D("h1MCdau1PtMother", "dau_{1}, p_{T,mother}", 25, bins_pt_new, "p_{T,mother} (GeV/c)",
202  m_histoList_MCParticles);
203 
204  m_h2_MC_dau1_2D = createHistogram2D("h2MCdau1", "entry per MC dau_{1}",
205  10, bins_theta, "#theta",
206  25, bins_pt_new, "p_{t} (GeV/c)", m_histoList_MCParticles);
207 
208  m_h2_MC_dau1_2D_BP = createHistogram2D("h2MCdau1BP", "entry per MC dau_{1}, beam pipe",
209  10, bins_theta, "#theta",
210  25, bins_pt_new, "p_{t} (GeV/c)", m_histoList_MCParticles);
211 
212  m_h2_MC_dau1_2DMother = createHistogram2D("h2MCdau1Mother", "entry mother per MC dau_{1}",
213  10, bins_theta, "#theta_{mother}",
214  25, bins_pt_new, "p_{T,mother} (GeV/c)", m_histoList_MCParticles);
215 
216  m_h2_MC_dau1_pVScostheta = createHistogram2D("h2MCdau1pVScostheta", "p_{CM} VS cos(#theta)_{CM}, dau_{1}",
217  20, -1., 1.,
218  "cos(#theta)_{CM}",
219  50, 0., 5.,
220  "p_{CM} [GeV]", m_histoList_MCParticles);
221 
222  m_h1_MC_dau1_PDG = createHistogram1D("h1MCdau1PDG", "PDG code, dau_{1}", 4600, -2300, 2300, "PDG", m_histoList_check);
223 
224  //MC mother
225  m_h1_MC_Mother_RMother = createHistogram1D("h1MCMotherRMother", "mother, R mother", 200, 0, 20, "R mother",
226  m_histoList_MCParticles);
227  m_h3_MC_Mother = createHistogram3D("h3MCMother", "entry per MCmother",
228  25, bins_pt_new, "p_{t} (GeV/c)",
229  10, bins_theta, "#theta",
230  14, bins_phi, "#phi", m_histoList_MCParticles);
231  m_h1_MC_Mother_pt = createHistogram1D("h1MCMotherPt", "mother, p_{T}", 25, bins_pt_new, "p_{T} (GeV/c)",
232  m_histoList_MCParticles);
233  m_h1_MC_Mother_pz = createHistogram1D("h1MCMotherPz", "mother, p_{z}", 25, bins_pt_new, "p_{z} (GeV/c)",
234  m_histoList_MCParticles);
235  m_h1_MC_Mother_p = createHistogram1D("h1MCMotherP", "mother, p", 25, bins_pt_new, "p (GeV/c)", m_histoList_MCParticles);
236  m_h1_MC_Mother_theta = createHistogram1D("h1MCMotherTheta", "mother, #theta", 10, bins_theta, "#theta",
237  m_histoList_MCParticles);
238  m_h1_MC_Mother_costheta = createHistogram1D("h1MCMotherCosTheta", "mother, cos#theta", 20, bins_costheta, "cos#theta",
239  m_histoList_MCParticles);
240  m_h1_MC_Mother_phi = createHistogram1D("h1MCMotherPhi", "mother, #phi", 14, bins_phi, "#phi", m_histoList_MCParticles);
241  m_h1_MC_Mother_phi_BW = createHistogram1D("h1MCMotherPhiBW", "mother, #phi", 14, bins_phi, "#phi", m_histoList_MCParticles);
242  m_h1_MC_Mother_phi_barrel = createHistogram1D("h1MCMotherPhibarrel", "mother, #phi", 14, bins_phi, "#phi",
243  m_histoList_MCParticles);
244  m_h1_MC_Mother_phi_FW = createHistogram1D("h1MCMotherPhiFW", "mother, #phi", 14, bins_phi, "#phi", m_histoList_MCParticles);
245 
246  m_h2_MC_Mother_2D = createHistogram2D("h2MCMother", "entry per MCmother",
247  10, bins_theta, "#theta",
248  25, bins_pt_new, "p_{t} (GeV/c)", m_histoList_MCParticles);
249 
250  m_h2_MC_Mother_2D_BP = createHistogram2D("h2MCMotherBP", "entry per MCmother, beam pipe",
251  10, bins_theta, "#theta",
252  25, bins_pt_new, "p_{t} (GeV/c)", m_histoList_MCParticles);
253 
254  m_h2_MC_Mother_pVScostheta = createHistogram2D("h2MCMotherpVScostheta", "p_{CM} VS cos(#theta)_{CM}, mother",
255  20, -1., 1.,
256  "cos(#theta)_{CM}",
257  50, 0., 5.,
258  "p_{CM} [GeV]", m_histoList_MCParticles);
259 
260  m_h1_MC_Mother_PDG = createHistogram1D("h1MCMotherPDG", "PDG code, mother", 6400, -3200, 3200, "PDG", m_histoList_check);
261 
262 
263  //------------------------------------------------------------------//
264  // TRACKS //
265  //------------------------------------------------------------------//
266 
267  //track dau0
268  m_h1_track_dau0_d0 = (TH1F*)duplicateHistogram("h1trackdau0D0", "d0 dau_{0}", m_h1_MC_dau0_d0, m_histoList_Tracks);
269  m_h1_track_dau0_z0 = (TH1F*)duplicateHistogram("h1trackdau0Z0", "z0 dau_{0}", m_h1_MC_dau0_z0, m_histoList_Tracks);
270  m_h1_track_dau0_RMother = (TH1F*)duplicateHistogram("h1trackdau0RMother", "dau_{0}, R mother", m_h1_MC_dau0_RMother,
271  m_histoList_Tracks);
272  m_h3_track_dau0 = (TH3F*)duplicateHistogram("h3trackdau0", "entry per track dau_{0}", m_h3_MC_dau0, m_histoList_Tracks);
273  m_h1_track_dau0_pt = (TH1F*)duplicateHistogram("h1trackdau0Pt", "p_{T} dau_{0}", m_h1_MC_dau0_pt, m_histoList_Tracks);
274  m_h1_track_dau0_pz = (TH1F*)duplicateHistogram("h1trackdau0Pz", "p_{z} dau_{0}", m_h1_MC_dau0_pz, m_histoList_Tracks);
275  m_h1_track_dau0_p = (TH1F*)duplicateHistogram("h1trackdau0P", "p dau_{0}", m_h1_MC_dau0_p, m_histoList_Tracks);
276  m_h1_track_dau0_theta = (TH1F*)duplicateHistogram("h1trackdau0Theta", "#theta dau_{0}", m_h1_MC_dau0_theta, m_histoList_Tracks);
277  m_h1_track_dau0_costheta = (TH1F*)duplicateHistogram("h1trackdau0CosTheta", "cos#theta dau_{0}", m_h1_MC_dau0_costheta,
278  m_histoList_Tracks);
279  m_h1_track_dau0_Mother_cosAngle = (TH1F*)duplicateHistogram("h1trackdau0MothercosAngle", "#alpha_{mother,dau_{0}}",
280  m_h1_MC_dau0_Mother_cosAngle, m_histoList_Tracks);
281  m_h1_track_dau0_phi = (TH1F*)duplicateHistogram("h1trackdau0Phi", "#phi dau_{0}", m_h1_MC_dau0_phi, m_histoList_Tracks);
282  m_h1_track_dau0_phi_BW = (TH1F*)duplicateHistogram("h1trackdau0PhiBW", "#phi dau_{0}", m_h1_MC_dau0_phi_BW, m_histoList_Tracks);
283  m_h1_track_dau0_phi_barrel = (TH1F*)duplicateHistogram("h1trackdau0Phibarrel", "#phi dau_{0}", m_h1_MC_dau0_phi_barrel,
284  m_histoList_Tracks);
285  m_h1_track_dau0_phi_FW = (TH1F*)duplicateHistogram("h1trackdau0PhiFW", "#phi dau_{0}", m_h1_MC_dau0_phi_FW, m_histoList_Tracks);
286 
287  m_h1_track_dau0_phiMother_total = (TH1F*)duplicateHistogram("h1trackdau0PhiMothertotal", "dau_{0}, #phi_{mother}",
288  m_h1_MC_dau0_phiMother_total, m_histoList_Tracks);
289  m_h1_track_dau0_phiMother_BW = (TH1F*)duplicateHistogram("h1trackdau0PhiMotherBW", "dau_{0}, #phi_{mother}, BW",
290  m_h1_MC_dau0_phiMother_BW, m_histoList_Tracks);
291  m_h1_track_dau0_phiMother_barrel = (TH1F*)duplicateHistogram("h1trackdau0PhiMotherbarrel", "dau_{0}, #phi_{mother}, barrel",
292  m_h1_MC_dau0_phiMother_barrel, m_histoList_Tracks);
293  m_h1_track_dau0_phiMother_FW = (TH1F*)duplicateHistogram("h1trackdau0PhiMotherFW", "dau_{0}, #phi_{mother}, FW",
294  m_h1_MC_dau0_phiMother_FW, m_histoList_Tracks);
295 
296  m_h1_track_dau0_thetaMother = (TH1F*)duplicateHistogram("h1trackdau0ThetaMother", "#theta_{mother} dau_{0}",
297  m_h1_MC_dau0_thetaMother, m_histoList_Tracks);
298  m_h1_track_dau0_ptMother = (TH1F*)duplicateHistogram("h1trackdau0PtMother", "p_{T,mother} dau_{0}", m_h1_MC_dau0_ptMother,
299  m_histoList_Tracks);
300 
301  m_h2_track_dau0_2D = (TH2F*)duplicateHistogram("h2trackdau02D", "p_{T} VS #theta, dau_{0}", m_h2_MC_dau0_2D, m_histoList_Tracks);
302  m_h2_track_dau0_2D_BP = (TH2F*)duplicateHistogram("h2trackdau02dBP", "p_{T} VS #theta, dau_{0} BP", m_h2_MC_dau0_2D_BP,
303  m_histoList_Tracks);
304  m_h2_track_dau0_2DMother = (TH2F*)duplicateHistogram("h2trackdau02DMother", "p_{T,mother} VS #theta_{mother}, dau_{0}",
305  m_h2_MC_dau0_2DMother, m_histoList_Tracks);
306 
307  m_h2_track_dau0_pVScostheta = (TH2F*)duplicateHistogram("h2trackdau0pVScostheta", "p VS cos(#theta), dau_{0}",
308  m_h2_MC_dau0_pVScostheta, m_histoList_Tracks);
309 
310 
311  //track dau1
312  m_h1_track_dau1_d0 = (TH1F*)duplicateHistogram("h1trackdau1D0", "d0 dau_{1}", m_h1_MC_dau1_d0, m_histoList_Tracks);
313  m_h1_track_dau1_z0 = (TH1F*)duplicateHistogram("h1trackdau1Z0", "z0 dau_{1}", m_h1_MC_dau1_z0, m_histoList_Tracks);
314  m_h1_track_dau1_RMother = (TH1F*)duplicateHistogram("h1trackdau1RMother", "p, R mother", m_h1_MC_dau1_RMother, m_histoList_Tracks);
315  m_h3_track_dau1 = (TH3F*)duplicateHistogram("h3trackdau1", "entry per track dau_{1}", m_h3_MC_dau1, m_histoList_Tracks);
316  m_h1_track_dau1_pt = (TH1F*)duplicateHistogram("h1trackdau1Pt", "p_{T} dau_{1}", m_h1_MC_dau1_pt, m_histoList_Tracks);
317  m_h1_track_dau1_pz = (TH1F*)duplicateHistogram("h1trackdau1Pz", "p_{z} dau_{1}", m_h1_MC_dau1_pz, m_histoList_Tracks);
318  m_h1_track_dau1_p = (TH1F*)duplicateHistogram("h1trackdau1P", "p dau_{1}", m_h1_MC_dau1_p, m_histoList_Tracks);
319  m_h1_track_dau1_theta = (TH1F*)duplicateHistogram("h1trackdau1Theta", "#theta dau_{1}", m_h1_MC_dau1_theta, m_histoList_Tracks);
320  m_h1_track_dau1_costheta = (TH1F*)duplicateHistogram("h1trackdau1CosTheta", "cos#theta dau_{1}", m_h1_MC_dau1_costheta,
321  m_histoList_Tracks);
322  m_h1_track_dau1_Mother_cosAngle = (TH1F*)duplicateHistogram("h1trackdau1MothercosAngle", "#alpha_{mother,p}",
323  m_h1_MC_dau1_Mother_cosAngle, m_histoList_Tracks);
324  m_h1_track_dau1_phi = (TH1F*)duplicateHistogram("h1trackdau1Phi", "#phi dau_{1}", m_h1_MC_dau1_phi, m_histoList_Tracks);
325  m_h1_track_dau1_phi_BW = (TH1F*)duplicateHistogram("h1trackdau1PhiBW", "#phi dau_{1}", m_h1_MC_dau1_phi_BW, m_histoList_Tracks);
326  m_h1_track_dau1_phi_barrel = (TH1F*)duplicateHistogram("h1trackdau1Phibarrel", "#phi dau_{1}", m_h1_MC_dau1_phi_barrel,
327  m_histoList_Tracks);
328  m_h1_track_dau1_phi_FW = (TH1F*)duplicateHistogram("h1trackdau1PhiFW", "#phi dau_{1}", m_h1_MC_dau1_phi_FW, m_histoList_Tracks);
329 
330  m_h1_track_dau1_phiMother_total = (TH1F*)duplicateHistogram("h1trackdau1PhiMothertotal", "dau_{1}, #phi_{mother}",
331  m_h1_MC_dau1_phiMother_total, m_histoList_Tracks);
332  m_h1_track_dau1_phiMother_BW = (TH1F*)duplicateHistogram("h1trackdau1PhiMotherBW", "dau_{1}, #phi_{mother}, BW",
333  m_h1_MC_dau1_phiMother_BW, m_histoList_Tracks);
334  m_h1_track_dau1_phiMother_barrel = (TH1F*)duplicateHistogram("h1trackdau1PhiMotherbarrel", "dau_{1}, #phi_{mother}, barrel",
335  m_h1_MC_dau1_phiMother_barrel, m_histoList_Tracks);
336  m_h1_track_dau1_phiMother_FW = (TH1F*)duplicateHistogram("h1trackdau1PhiMotherFW", "dau_{1}, #phi_{mother}, FW",
337  m_h1_MC_dau1_phiMother_FW, m_histoList_Tracks);
338 
339  m_h1_track_dau1_thetaMother = (TH1F*)duplicateHistogram("h1trackdau1ThetaMother", "#theta_{mother} dau_{1}",
340  m_h1_MC_dau1_thetaMother, m_histoList_Tracks);
341  m_h1_track_dau1_ptMother = (TH1F*)duplicateHistogram("h1trackdau1PtMother", "p_{T,mother} dau_{1}", m_h1_MC_dau1_ptMother,
342  m_histoList_Tracks);
343 
344  m_h2_track_dau1_2D = (TH2F*)duplicateHistogram("h2trackdau12D", "p_{T} VS #theta, dau_{1}", m_h2_MC_dau1_2D, m_histoList_Tracks);
345  m_h2_track_dau1_2D_BP = (TH2F*)duplicateHistogram("h2trackdau12dBP", "p_{T} VS #theta, p BDAU_{1}", m_h2_MC_dau1_2D_BP,
346  m_histoList_Tracks);
347  m_h2_track_dau1_2DMother = (TH2F*)duplicateHistogram("h2trackdau12DMother", "p_{T,mother} VS #theta_{mother}, dau_{1}",
348  m_h2_MC_dau1_2DMother, m_histoList_Tracks);
349 
350  m_h2_track_dau1_pVScostheta = (TH2F*)duplicateHistogram("h2trackdau1pVScostheta", "p VS cos(#theta), dau_{1}",
351  m_h2_MC_dau1_pVScostheta, m_histoList_Tracks);
352 
353  //V0
354  m_h1_V0_RMother = (TH1F*)duplicateHistogram("h1V0RMother", "mother, R mother", m_h1_MC_Mother_RMother, m_histoList_Tracks);
355  m_h3_V0 = (TH3F*)duplicateHistogram("h3V0", "entry per V0", m_h3_MC_Mother, m_histoList_Tracks);
356  m_h1_V0_pt = (TH1F*)duplicateHistogram("h1V0Pt", "p_{T} mother", m_h1_MC_Mother_pt, m_histoList_Tracks);
357  m_h1_V0_pz = (TH1F*)duplicateHistogram("h1V0Pz", "p_{z} mother", m_h1_MC_Mother_pz, m_histoList_Tracks);
358  m_h1_V0_p = (TH1F*)duplicateHistogram("h1V0P", "p mother", m_h1_MC_Mother_p, m_histoList_Tracks);
359  m_h1_V0_theta = (TH1F*)duplicateHistogram("h1V0Theta", "#theta mother", m_h1_MC_Mother_theta, m_histoList_Tracks);
360  m_h1_V0_costheta = (TH1F*)duplicateHistogram("h1V0CosTheta", "cos#theta mother", m_h1_MC_Mother_costheta, m_histoList_Tracks);
361  m_h1_V0_phi = (TH1F*)duplicateHistogram("h1V0Phi", "#phi mother", m_h1_MC_Mother_phi, m_histoList_Tracks);
362  m_h1_V0_phi_BW = (TH1F*)duplicateHistogram("h1V0PhiBW", "#phi mother", m_h1_MC_Mother_phi_BW, m_histoList_Tracks);
363  m_h1_V0_phi_barrel = (TH1F*)duplicateHistogram("h1V0Phibarrel", "#phi mother", m_h1_MC_Mother_phi_barrel, m_histoList_Tracks);
364  m_h1_V0_phi_FW = (TH1F*)duplicateHistogram("h1V0PhiFW", "#phi mother", m_h1_MC_Mother_phi_FW, m_histoList_Tracks);
365 
366  m_h2_V0_Mother_2D = (TH2F*)duplicateHistogram("h2V0Mother2D", "p_{T} VS #theta, mother", m_h2_MC_Mother_2D, m_histoList_Tracks);
367  m_h2_V0_Mother_2D_BP = (TH2F*)duplicateHistogram("h2V0Mother2dBP", "p_{T} VS #theta, mother BP", m_h2_MC_Mother_2D_BP,
368  m_histoList_Tracks);
369 
370  m_h2_V0_Mother_pVScostheta = (TH2F*)duplicateHistogram("h2V0MotherpVScostheta", "p VS cos(#theta), mother",
371  m_h2_MC_Mother_pVScostheta, m_histoList_Tracks);
372 
373 
374  //------------------------------------------------------------------//
375  // RECO TRACKS //
376  //------------------------------------------------------------------//
377 
378  //RecoTrack dau0
379  m_h1_RecoTrack_dau0_d0 = (TH1F*)duplicateHistogram("h1RecoTrackdau0D0", "d0 dau_{0}", m_h1_MC_dau0_d0, m_histoList_RecoTracks);
380  m_h1_RecoTrack_dau0_z0 = (TH1F*)duplicateHistogram("h1RecoTrackdau0Z0", "z0 dau_{0}", m_h1_MC_dau0_z0, m_histoList_RecoTracks);
381  m_h1_RecoTrack_dau0_RMother = (TH1F*)duplicateHistogram("h1RecoTrackdau0RMother", "dau_{0}, R mother", m_h1_MC_dau0_RMother,
382  m_histoList_RecoTracks);
383  m_h3_RecoTrack_dau0 = (TH3F*)duplicateHistogram("h3RecoTrackdau0", "entry per RecoTrack dau_{0}", m_h3_MC_dau0,
384  m_histoList_RecoTracks);
385  m_h1_RecoTrack_dau0_pt = (TH1F*)duplicateHistogram("h1RecoTrackdau0Pt", "p_{T} dau_{0}", m_h1_MC_dau0_pt, m_histoList_RecoTracks);
386  m_h1_RecoTrack_dau0_pz = (TH1F*)duplicateHistogram("h1RecoTrackdau0Pz", "p_{z} dau_{0}", m_h1_MC_dau0_pz, m_histoList_RecoTracks);
387  m_h1_RecoTrack_dau0_p = (TH1F*)duplicateHistogram("h1RecoTrackdau0P", "p dau_{0}", m_h1_MC_dau0_p, m_histoList_RecoTracks);
388  m_h1_RecoTrack_dau0_theta = (TH1F*)duplicateHistogram("h1RecoTrackdau0Theta", "#theta dau_{0}", m_h1_MC_dau0_theta,
389  m_histoList_RecoTracks);
390  m_h1_RecoTrack_dau0_costheta = (TH1F*)duplicateHistogram("h1RecoTrackdau0CosTheta", "cos#theta dau_{0}", m_h1_MC_dau0_costheta,
391  m_histoList_RecoTracks);
392  m_h1_RecoTrack_dau0_Mother_cosAngle = (TH1F*)duplicateHistogram("h1RecoTrackdau0MothercosAngle", "#alpha_{mother,dau_{0}}",
393  m_h1_MC_dau0_Mother_cosAngle, m_histoList_RecoTracks);
394  m_h1_RecoTrack_dau0_phi = (TH1F*)duplicateHistogram("h1RecoTrackdau0Phi", "#phi dau_{0}", m_h1_MC_dau0_phi, m_histoList_RecoTracks);
395  m_h1_RecoTrack_dau0_phi_BW = (TH1F*)duplicateHistogram("h1RecoTrackdau0PhiBW", "#phi dau_{0}", m_h1_MC_dau0_phi_BW,
396  m_histoList_RecoTracks);
397  m_h1_RecoTrack_dau0_phi_barrel = (TH1F*)duplicateHistogram("h1RecoTrackdau0Phibarrel", "#phi dau_{0}", m_h1_MC_dau0_phi_barrel,
398  m_histoList_RecoTracks);
399  m_h1_RecoTrack_dau0_phi_FW = (TH1F*)duplicateHistogram("h1RecoTrackdau0PhiFW", "#phi dau_{0}", m_h1_MC_dau0_phi_FW,
400  m_histoList_RecoTracks);
401 
402  m_h1_RecoTrack_dau0_phiMother_total = (TH1F*)duplicateHistogram("h1RecoTrackdau0PhiMothertotal", "dau_{0}, #phi_{mother}",
403  m_h1_MC_Mother_phi, m_histoList_RecoTracks);
404  m_h1_RecoTrack_dau0_phiMother_BW = (TH1F*)duplicateHistogram("h1RecoTrackdau0PhiMotherBW", "dau_{0}, #phi_{mother}, BW",
405  m_h1_MC_dau0_phiMother_BW, m_histoList_RecoTracks);
406  m_h1_RecoTrack_dau0_phiMother_barrel = (TH1F*)duplicateHistogram("h1RecoTrackdau0PhiMotherbarrel", "dau_{0}, #phi_{mother}, barrel",
407  m_h1_MC_dau0_phiMother_barrel, m_histoList_RecoTracks);
408  m_h1_RecoTrack_dau0_phiMother_FW = (TH1F*)duplicateHistogram("h1RecoTrackdau0PhiMotherFW", "dau_{0}, #phi_{mother}, FW",
409  m_h1_MC_dau0_phiMother_FW, m_histoList_RecoTracks);
410 
411  m_h1_RecoTrack_dau0_thetaMother = (TH1F*)duplicateHistogram("h1RecoTrackdau0ThetaMother", "#theta_{mother} dau_{0}",
412  m_h1_MC_dau0_theta, m_histoList_RecoTracks);
413  m_h1_RecoTrack_dau0_ptMother = (TH1F*)duplicateHistogram("h1RecoTrackdau0PtMother", "p_{T,mother} dau_{0}", m_h1_MC_dau0_pt,
414  m_histoList_RecoTracks);
415 
416  m_h2_RecoTrack_dau0_2D = (TH2F*)duplicateHistogram("h2RecoTrackdau02D", "p_{T} VS #theta, dau_{0}", m_h2_MC_dau0_2D,
417  m_histoList_RecoTracks);
418  m_h2_RecoTrack_dau0_2D_BP = (TH2F*)duplicateHistogram("h2RecoTrackdau02dBP", "p_{T} VS #theta, dau_{0} BP", m_h2_MC_dau0_2D_BP,
419  m_histoList_RecoTracks);
420  m_h2_RecoTrack_dau0_2DMother = (TH2F*)duplicateHistogram("h2RecoTrackdau02DMother", "p_{T,mother} VS #theta_{mother}, dau_{0}",
421  m_h2_MC_dau0_2DMother, m_histoList_RecoTracks);
422 
423  m_h2_RecoTrack_dau0_pVScostheta = (TH2F*)duplicateHistogram("h2RecoTrackdau0pVscostheta", "p VS cos(#theta), dau_{0}",
424  m_h2_MC_dau0_pVScostheta, m_histoList_RecoTracks);
425 
426  //RecoTrack dau1
427  m_h1_RecoTrack_dau1_d0 = (TH1F*)duplicateHistogram("h1RecoTrackdau1D0", "d0 dau_{1}", m_h1_MC_dau1_d0, m_histoList_RecoTracks);
428  m_h1_RecoTrack_dau1_z0 = (TH1F*)duplicateHistogram("h1RecoTrackdau1Z0", "z0 dau_{1}", m_h1_MC_dau1_z0, m_histoList_RecoTracks);
429  m_h1_RecoTrack_dau1_RMother = (TH1F*)duplicateHistogram("h1RecoTrackdau1RMother", "dau_{1}, R mother", m_h1_MC_dau1_RMother,
430  m_histoList_RecoTracks);
431  m_h3_RecoTrack_dau1 = (TH3F*)duplicateHistogram("h3RecoTrackdau1", "entry per RecoTrack dau_{1}", m_h3_MC_dau1,
432  m_histoList_RecoTracks);
433  m_h1_RecoTrack_dau1_pt = (TH1F*)duplicateHistogram("h1RecoTrackdau1Pt", "p_{T} dau_{1}", m_h1_MC_dau1_pt, m_histoList_RecoTracks);
434  m_h1_RecoTrack_dau1_pz = (TH1F*)duplicateHistogram("h1RecoTrackdau1Pz", "p_{z} dau_{1}", m_h1_MC_dau1_pz, m_histoList_RecoTracks);
435  m_h1_RecoTrack_dau1_p = (TH1F*)duplicateHistogram("h1RecoTrackdau1P", "p dau_{1}", m_h1_MC_dau1_p, m_histoList_RecoTracks);
436  m_h1_RecoTrack_dau1_theta = (TH1F*)duplicateHistogram("h1RecoTrackdau1Theta", "#theta dau_{1}", m_h1_MC_dau1_theta,
437  m_histoList_RecoTracks);
438  m_h1_RecoTrack_dau1_costheta = (TH1F*)duplicateHistogram("h1RecoTrackdau1CosTheta", "cos#theta dau_{1}", m_h1_MC_dau1_costheta,
439  m_histoList_RecoTracks);
440  m_h1_RecoTrack_dau1_Mother_cosAngle = (TH1F*)duplicateHistogram("h1RecoTrackdau1MothercosAngle", "#alpha_{mother,p}",
441  m_h1_MC_dau1_Mother_cosAngle, m_histoList_RecoTracks);
442  m_h1_RecoTrack_dau1_phi = (TH1F*)duplicateHistogram("h1RecoTrackdau1Phi", "#phi dau_{1}", m_h1_MC_dau1_phi, m_histoList_RecoTracks);
443  m_h1_RecoTrack_dau1_phi_BW = (TH1F*)duplicateHistogram("h1RecoTrackdau1PhiBW", "#phi dau_{1}", m_h1_MC_dau1_phi_BW,
444  m_histoList_RecoTracks);
445  m_h1_RecoTrack_dau1_phi_barrel = (TH1F*)duplicateHistogram("h1RecoTrackdau1Phibarrel", "#phi dau_{1}", m_h1_MC_dau1_phi_barrel,
446  m_histoList_RecoTracks);
447  m_h1_RecoTrack_dau1_phi_FW = (TH1F*)duplicateHistogram("h1RecoTrackdau1PhiFW", "#phi dau_{1}", m_h1_MC_dau1_phi_FW,
448  m_histoList_RecoTracks);
449 
450  m_h1_RecoTrack_dau1_phiMother_total = (TH1F*)duplicateHistogram("h1RecoTrackdau1PhiMothertotal", "dau_{1}, #phi_{mother}",
451  m_h1_MC_dau1_phiMother_total, m_histoList_RecoTracks);
452  m_h1_RecoTrack_dau1_phiMother_BW = (TH1F*)duplicateHistogram("h1RecoTrackdau1PhiMotherBW", "dau_{1}, #phi_{mother}, BW",
453  m_h1_MC_dau1_phiMother_BW, m_histoList_RecoTracks);
454  m_h1_RecoTrack_dau1_phiMother_barrel = (TH1F*)duplicateHistogram("h1RecoTrackdau1PhiMotherbarrel", "dau_{1}, #phi_{mother}, barrel",
455  m_h1_MC_dau1_phiMother_barrel, m_histoList_RecoTracks);
456  m_h1_RecoTrack_dau1_phiMother_FW = (TH1F*)duplicateHistogram("h1RecoTrackdau1PhiMotherFW", "dau_{1}, #phi_{mother}, FW",
457  m_h1_MC_dau1_phiMother_FW, m_histoList_RecoTracks);
458 
459  m_h1_RecoTrack_dau1_thetaMother = (TH1F*)duplicateHistogram("h1RecoTrackdau1ThetaMother", "#theta_{mother} dau_{1}",
460  m_h1_MC_dau1_theta, m_histoList_RecoTracks);
461  m_h1_RecoTrack_dau1_ptMother = (TH1F*)duplicateHistogram("h1RecoTrackdau1PtMother", "p_{T,mother} dau_{1}", m_h1_MC_dau1_pt,
462  m_histoList_RecoTracks);
463 
464  m_h2_RecoTrack_dau1_2D = (TH2F*)duplicateHistogram("h2RecoTrackdau12D", "p_{T} VS #theta, dau_{1}", m_h2_MC_dau1_2D,
465  m_histoList_RecoTracks);
466  m_h2_RecoTrack_dau1_2D_BP = (TH2F*)duplicateHistogram("h2RecoTrackdau12dBP", "p_{T} VS #theta, p BP", m_h2_MC_dau1_2D_BP,
467  m_histoList_RecoTracks);
468  m_h2_RecoTrack_dau1_2DMother = (TH2F*)duplicateHistogram("h2RecoTrackdau12DMother", "p_{T,mother} VS #theta_{mother}, dau_{1}",
469  m_h2_MC_dau1_2DMother, m_histoList_RecoTracks);
470 
471  m_h2_RecoTrack_dau1_pVScostheta = (TH2F*)duplicateHistogram("h2RecoTrackdau1pVscostheta", "p VS cos(#theta), dau_{1}",
472  m_h2_MC_dau1_pVScostheta, m_histoList_RecoTracks);
473 
474  //RecoTrack Mother
475  m_h1_RecoTrack_Mother_RMother = (TH1F*)duplicateHistogram("h1RecoTrack_MotherRMother", "mother, R mother", m_h1_MC_Mother_RMother,
476  m_histoList_RecoTracks);
477  m_h3_RecoTrack_Mother = (TH3F*)duplicateHistogram("h3RecoTrack_Mother", "entry per RecoTrack_Mother", m_h3_MC_Mother,
478  m_histoList_RecoTracks);
479  m_h1_RecoTrack_Mother_pt = (TH1F*)duplicateHistogram("h1RecoTrack_MotherPt", "p_{T} mother", m_h1_MC_dau0_pt,
480  m_histoList_RecoTracks);
481  m_h1_RecoTrack_Mother_pz = (TH1F*)duplicateHistogram("h1RecoTrack_MotherPz", "p_{z} mother", m_h1_MC_dau0_pz,
482  m_histoList_RecoTracks);
483  m_h1_RecoTrack_Mother_p = (TH1F*)duplicateHistogram("h1RecoTrack_MotherP", "p mother", m_h1_MC_dau0_p, m_histoList_RecoTracks);
484  m_h1_RecoTrack_Mother_theta = (TH1F*)duplicateHistogram("h1RecoTrack_MotherTheta", "#theta mother", m_h1_MC_dau0_theta,
485  m_histoList_RecoTracks);
486  m_h1_RecoTrack_Mother_costheta = (TH1F*)duplicateHistogram("h1RecoTrack_MotherCosTheta", "cos#theta mother", m_h1_MC_dau0_costheta,
487  m_histoList_RecoTracks);
488  m_h1_RecoTrack_Mother_phi = (TH1F*)duplicateHistogram("h1RecoTrack_MotherPhi", "#phi mother", m_h1_MC_dau0_phi,
489  m_histoList_RecoTracks);
490  m_h1_RecoTrack_Mother_phi_BW = (TH1F*)duplicateHistogram("h1RecoTrack_MotherPhiBW", "#phi mother", m_h1_MC_dau0_phi_BW,
491  m_histoList_RecoTracks);
492  m_h1_RecoTrack_Mother_phi_barrel = (TH1F*)duplicateHistogram("h1RecoTrack_MotherPhibarrel", "#phi mother", m_h1_MC_dau0_phi_barrel,
493  m_histoList_RecoTracks);
494  m_h1_RecoTrack_Mother_phi_FW = (TH1F*)duplicateHistogram("h1RecoTrack_MotherPhiFW", "#phi mother", m_h1_MC_dau0_phi_FW,
495  m_histoList_RecoTracks);
496 
497  m_h2_RecoTrack_Mother_2D = (TH2F*)duplicateHistogram("h2RecoTrack_Mother2D", "p_{T} VS #theta, mother", m_h2_MC_Mother_2D,
498  m_histoList_RecoTracks);
499  m_h2_RecoTrack_Mother_2D_BP = (TH2F*)duplicateHistogram("h2RecoTrack_Mother2dBP", "p_{T} VS #theta, mother BP",
500  m_h2_MC_Mother_2D_BP, m_histoList_RecoTracks);
501 
502  m_h2_RecoTrack_Mother_pVScostheta = (TH2F*)duplicateHistogram("h2RecoTrackMotherpVscostheta", "p VS cos(#theta), mother",
503  m_h2_MC_Mother_pVScostheta, m_histoList_RecoTracks);
504 
505 }
506 
507 void EffPlotsModule::beginRun()
508 {
509 
510 }
511 
512 void EffPlotsModule::event()
513 {
514 
515  ROOT::Math::XYZVector magField = BFieldManager::getField(0, 0, 0) / Unit::T;
516 
517  B2DEBUG(29, "+++++ 1. loop on MCParticles");
518  for (const MCParticle& mcParticle : m_MCParticles) {
519 
520  //------------------------------------------------------------------//
521  // MC PARTICLES //
522  //------------------------------------------------------------------//
523 
524  if (m_V0sType == "Lambda0") {
525  if (! isLambda0(mcParticle))
526  continue;
527 
528  std::vector< MCParticle* > MCPart_dau = mcParticle.getDaughters();
529 
530  if (abs(MCPart_dau[0]->getPDG()) == Const::pion.getPDGCode() && abs(MCPart_dau[1]->getPDG()) == Const::proton.getPDGCode()) {
531  m_MCDaughter0 = MCPart_dau[0];
532  m_MCDaughter1 = MCPart_dau[1];
533  } else if (abs(MCPart_dau[0]->getPDG()) == Const::proton.getPDGCode() && abs(MCPart_dau[1]->getPDG()) == Const::pion.getPDGCode()) {
534  m_MCDaughter0 = MCPart_dau[1];
535  m_MCDaughter1 = MCPart_dau[0];
536  } else B2INFO("Lambda daughters != pi & p");
537  }
538 
539  else if (m_V0sType == "Ks") {
540  if (! isK_Short(mcParticle))
541  continue;
542 
543  std::vector< MCParticle* > MCPart_dau = mcParticle.getDaughters();
544 
545  if (MCPart_dau[0]->getPDG() == Const::pion.getPDGCode() && MCPart_dau[1]->getPDG() == -Const::pion.getPDGCode()) {
546  m_MCDaughter0 = MCPart_dau[0];
547  m_MCDaughter1 = MCPart_dau[1];
548  } else if (MCPart_dau[0]->getPDG() == -Const::pion.getPDGCode() && MCPart_dau[1]->getPDG() == Const::pion.getPDGCode()) {
549  m_MCDaughter0 = MCPart_dau[1];
550  m_MCDaughter1 = MCPart_dau[0];
551  } else B2INFO("Ks daughters != pi+ & pi-");
552  }
553 
554  MCParticleInfo mcParticleInfo(mcParticle, magField);
555  MCParticleInfo mcParticleInfo_dau0(*m_MCDaughter0, magField);
556  MCParticleInfo mcParticleInfo_dau1(*m_MCDaughter1, magField);
557 
558  const ROOT::Math::XYZVector& MC_vtx = mcParticle.getDecayVertex();
559 
560  float MC_transDist = sqrt(MC_vtx.X() * MC_vtx.X() + MC_vtx.Y() * MC_vtx.Y());
561  float MC_pt = mcParticle.getMomentum().Rho();
562  float MC_p = mcParticle.getMomentum().R();
563  float MC_phi = mcParticle.getMomentum().Phi();
564  float MC_theta = mcParticle.getMomentum().Theta();
565  float MC_costheta = cos(mcParticle.getMomentum().Theta());
566 
567  m_h1_MC_dau0_d0->Fill(mcParticleInfo_dau0.getD0());
568  m_h1_MC_dau0_z0->Fill(mcParticleInfo_dau0.getZ0());
569  m_h1_MC_dau0_RMother->Fill(MC_transDist);
570  m_h3_MC_dau0->Fill(m_MCDaughter0->getMomentum().Rho(), m_MCDaughter0->getMomentum().Theta(), m_MCDaughter0->getMomentum().Phi());
571  m_h1_MC_dau0_pt->Fill(m_MCDaughter0->getMomentum().Rho());
572  m_h1_MC_dau0_pz->Fill(m_MCDaughter0->getMomentum().z());
573  m_h1_MC_dau0_p->Fill(m_MCDaughter0->getMomentum().R());
574  m_h1_MC_dau0_phi->Fill(m_MCDaughter0->getMomentum().Phi());
575  m_h1_MC_dau0_theta->Fill(m_MCDaughter0->getMomentum().Theta());
576  m_h1_MC_dau0_costheta->Fill(cos(m_MCDaughter0->getMomentum().Theta()));
577  m_h1_MC_dau0_Mother_cosAngle->Fill(mcParticle.getMomentum().Dot(m_MCDaughter0->getMomentum()) / mcParticle.getMomentum().R() /
578  m_MCDaughter0->getMomentum().R());
579 
580  m_h1_MC_dau0_thetaMother->Fill(MC_theta);
581  m_h1_MC_dau0_ptMother->Fill(MC_pt);
582 
583  m_h1_MC_dau0_phiMother_total->Fill(MC_phi);
584  m_h2_MC_dau0_2D->Fill(m_MCDaughter0->getMomentum().Theta(), m_MCDaughter0->getMomentum().Rho());
585  m_h2_MC_dau0_2DMother->Fill(MC_theta, MC_pt);
586  m_h2_MC_dau0_pVScostheta->Fill(cos(m_MCDaughter0->getMomentum().Theta()), m_MCDaughter0->getMomentum().R());
587  m_h1_MC_dau0_PDG->Fill(m_MCDaughter0->getPDG());
588 
589  m_h1_MC_dau1_d0->Fill(mcParticleInfo_dau1.getD0());
590  m_h1_MC_dau1_z0->Fill(mcParticleInfo_dau1.getZ0());
591  m_h1_MC_dau1_RMother->Fill(MC_transDist);
592  m_h3_MC_dau1->Fill(m_MCDaughter1->getMomentum().Rho(), m_MCDaughter1->getMomentum().Theta(), m_MCDaughter1->getMomentum().Phi());
593  m_h1_MC_dau1_pt->Fill(m_MCDaughter1->getMomentum().Rho());
594  m_h1_MC_dau1_pz->Fill(m_MCDaughter1->getMomentum().z());
595  m_h1_MC_dau1_p->Fill(m_MCDaughter1->getMomentum().R());
596  m_h1_MC_dau1_phi->Fill(m_MCDaughter1->getMomentum().Phi());
597  m_h1_MC_dau1_theta->Fill(m_MCDaughter1->getMomentum().Theta());
598  m_h1_MC_dau1_costheta->Fill(cos(m_MCDaughter1->getMomentum().Theta()));
599  m_h1_MC_dau1_Mother_cosAngle->Fill(mcParticle.getMomentum().Dot(m_MCDaughter1->getMomentum()) / mcParticle.getMomentum().R() /
600  m_MCDaughter1->getMomentum().R());
601 
602  m_h1_MC_dau1_thetaMother->Fill(MC_theta);
603  m_h1_MC_dau1_ptMother->Fill(MC_pt);
604 
605  m_h1_MC_dau1_phiMother_total->Fill(MC_phi);
606  m_h2_MC_dau1_2D->Fill(m_MCDaughter1->getMomentum().Theta(), m_MCDaughter1->getMomentum().Rho());
607  m_h2_MC_dau1_2DMother->Fill(MC_theta, MC_pt);
608  m_h2_MC_dau1_pVScostheta->Fill(cos(m_MCDaughter1->getMomentum().Theta()), m_MCDaughter1->getMomentum().R());
609  m_h1_MC_dau1_PDG->Fill(m_MCDaughter1->getPDG());
610 
611  m_h1_MC_Mother_RMother->Fill(MC_transDist);
612  m_h3_MC_Mother->Fill(MC_pt, MC_theta, MC_phi);
613  m_h1_MC_Mother_pt->Fill(MC_pt);
614  m_h1_MC_Mother_pz->Fill(mcParticle.getMomentum().z());
615  m_h1_MC_Mother_p->Fill(mcParticle.getMomentum().R());
616  m_h1_MC_Mother_phi->Fill(MC_phi);
617  m_h1_MC_Mother_theta->Fill(MC_theta);
618  m_h1_MC_Mother_costheta->Fill(cos(mcParticle.getMomentum().Theta()));
619  m_h2_MC_Mother_2D->Fill(MC_theta, MC_pt);
620  m_h2_MC_Mother_pVScostheta->Fill(MC_costheta, MC_p);
621  m_h1_MC_Mother_PDG->Fill(mcParticle.getPDG());
622 
623  //beam pipe
624  if (MC_transDist < 1.) {
625  m_h2_MC_dau0_2D_BP->Fill(m_MCDaughter0->getMomentum().Theta(), m_MCDaughter0->getMomentum().Rho());
626  m_h2_MC_dau1_2D_BP->Fill(m_MCDaughter1->getMomentum().Theta(), m_MCDaughter1->getMomentum().Rho());
627  m_h2_MC_Mother_2D_BP->Fill(MC_theta, MC_pt);
628  }
629 
630 
631  if (MC_theta > (120 * TMath::Pi() / 180.)) { //BW
632  m_h1_MC_dau0_phiMother_BW->Fill(MC_phi);
633  m_h1_MC_dau1_phiMother_BW->Fill(MC_phi);
634  m_h1_MC_dau0_phi_BW->Fill(m_MCDaughter0->getMomentum().Phi());
635  m_h1_MC_dau1_phi_BW->Fill(m_MCDaughter1->getMomentum().Phi());
636  m_h1_MC_Mother_phi_BW->Fill(MC_phi);
637  } else if (MC_theta < (30. * TMath::Pi() / 180.)) { //FW, theta < 30)
638  m_h1_MC_dau0_phiMother_FW->Fill(MC_phi);
639  m_h1_MC_dau1_phiMother_FW->Fill(MC_phi);
640  m_h1_MC_dau0_phi_FW->Fill(m_MCDaughter0->getMomentum().Phi());
641  m_h1_MC_dau1_phi_FW->Fill(m_MCDaughter1->getMomentum().Phi());
642  m_h1_MC_Mother_phi_FW->Fill(MC_phi);
643  } else { //barrel
644  m_h1_MC_dau0_phiMother_barrel->Fill(MC_phi);
645  m_h1_MC_dau1_phiMother_barrel->Fill(MC_phi);
646  m_h1_MC_dau0_phi_barrel->Fill(m_MCDaughter0->getMomentum().Phi());
647  m_h1_MC_dau1_phi_barrel->Fill(m_MCDaughter1->getMomentum().Phi());
648  m_h1_MC_Mother_phi_barrel->Fill(MC_phi);
649  }
650 
651  //------------------------------------------------------------------//
652  // MC RECO TRACKS //
653  //------------------------------------------------------------------//
654 
655  RelationVector<RecoTrack> MCRecoTracks_MCdau0 =
656  DataStore::getRelationsWithObj<RecoTrack>(m_MCDaughter0, m_MCRecoTracksName); //oppure &m_MCDaughter0;
657 
658  RelationVector<RecoTrack> MCRecoTracks_MCdau1 =
659  DataStore::getRelationsWithObj<RecoTrack>(m_MCDaughter1, m_MCRecoTracksName);
660 
661  if (MCRecoTracks_MCdau0.size() > 0) {
662  m_h1_RecoTrack_dau0_d0->Fill(mcParticleInfo_dau0.getD0());
663  m_h1_RecoTrack_dau0_z0->Fill(mcParticleInfo_dau0.getZ0());
664  m_h1_RecoTrack_dau0_RMother->Fill(MC_transDist);
665  m_h3_RecoTrack_dau0->Fill(m_MCDaughter0->getMomentum().Rho(), m_MCDaughter0->getMomentum().Theta(),
666  m_MCDaughter0->getMomentum().Phi());
667  m_h1_RecoTrack_dau0_pt->Fill(m_MCDaughter0->getMomentum().Rho());
668  m_h1_RecoTrack_dau0_pz->Fill(m_MCDaughter0->getMomentum().z());
669  m_h1_RecoTrack_dau0_p->Fill(m_MCDaughter0->getMomentum().R());
670  m_h1_RecoTrack_dau0_phi->Fill(m_MCDaughter0->getMomentum().Phi());
671  m_h1_RecoTrack_dau0_theta->Fill(m_MCDaughter0->getMomentum().Theta());
672  m_h1_RecoTrack_dau0_costheta->Fill(cos(m_MCDaughter0->getMomentum().Theta()));
673  m_h1_RecoTrack_dau0_Mother_cosAngle->Fill(mcParticle.getMomentum().Dot(m_MCDaughter0->getMomentum()) /
674  mcParticle.getMomentum().R() /
675  m_MCDaughter0->getMomentum().R());
676 
677  m_h1_RecoTrack_dau0_thetaMother->Fill(MC_theta);
678  m_h1_RecoTrack_dau0_ptMother->Fill(MC_pt);
679 
680  m_h1_RecoTrack_dau0_phiMother_total->Fill(MC_phi);
681  m_h2_RecoTrack_dau0_2D->Fill(m_MCDaughter0->getMomentum().Theta(), m_MCDaughter0->getMomentum().Rho());
682  m_h2_RecoTrack_dau0_2DMother->Fill(MC_theta, MC_pt);
683  m_h2_RecoTrack_dau0_pVScostheta->Fill(cos(m_MCDaughter0->getMomentum().Theta()), m_MCDaughter0->getMomentum().R());
684 
685  if (MC_transDist < 1.) {
686  m_h2_RecoTrack_dau0_2D_BP->Fill(m_MCDaughter0->getMomentum().Theta(), m_MCDaughter0->getMomentum().Rho());
687  }
688 
689  if (m_MCDaughter0->getMomentum().Theta() > (120 * TMath::Pi() / 180.)) m_h1_RecoTrack_dau0_phi_BW->Fill(
690  m_MCDaughter0->getMomentum().Phi());
691  else if (m_MCDaughter0->getMomentum().Theta() < (30. * TMath::Pi() / 180.)) m_h1_RecoTrack_dau0_phi_FW->Fill(
692  m_MCDaughter0->getMomentum().Phi());
693  else m_h1_RecoTrack_dau0_phi_barrel->Fill(m_MCDaughter0->getMomentum().Phi());
694 
695  if (MC_theta > (120 * TMath::Pi() / 180.)) m_h1_RecoTrack_dau0_phiMother_BW->Fill(MC_phi);
696  else if (MC_theta < (30. * TMath::Pi() / 180.)) m_h1_RecoTrack_dau0_phiMother_FW->Fill(MC_phi);
697  else m_h1_RecoTrack_dau0_phiMother_barrel->Fill(MC_phi);
698 
699  //V0 candidates
700  if (MCRecoTracks_MCdau1.size() > 0) {
701  m_h1_RecoTrack_Mother_RMother->Fill(MC_transDist);
702  m_h3_RecoTrack_Mother->Fill(MC_pt, MC_theta, MC_phi);
703  m_h1_RecoTrack_Mother_pt->Fill(MC_pt);
704  m_h1_RecoTrack_Mother_pz->Fill(mcParticle.getMomentum().z());
705  m_h1_RecoTrack_Mother_p->Fill(mcParticle.getMomentum().R());
706  m_h1_RecoTrack_Mother_phi->Fill(MC_phi);
707  m_h1_RecoTrack_Mother_theta->Fill(MC_theta);
708  m_h1_RecoTrack_Mother_costheta->Fill(cos(mcParticle.getMomentum().Theta()));
709 
710  m_h1_RecoTrack_Mother_pt->Fill(MC_pt);
711  m_h2_RecoTrack_Mother_2D->Fill(MC_theta, MC_pt);
712 
713  m_h2_RecoTrack_Mother_pVScostheta->Fill(MC_costheta, MC_p);
714 
715  if (MC_transDist < 1.) m_h2_RecoTrack_Mother_2D_BP->Fill(MC_theta, MC_pt);
716 
717  if (MC_theta > (120 * TMath::Pi() / 180.)) //BW
718  m_h1_RecoTrack_Mother_phi_BW->Fill(MC_phi);
719  else if (MC_theta < (30. * TMath::Pi() / 180.))//FW
720  m_h1_RecoTrack_Mother_phi_FW->Fill(MC_phi);
721  else
722  m_h1_RecoTrack_Mother_phi_barrel->Fill(MC_phi);
723  }
724  }
725 
726  if (MCRecoTracks_MCdau1.size() > 0) {
727  m_h1_RecoTrack_dau1_d0->Fill(mcParticleInfo_dau1.getD0());
728  m_h1_RecoTrack_dau1_z0->Fill(mcParticleInfo_dau1.getZ0());
729  m_h1_RecoTrack_dau1_RMother->Fill(MC_transDist);
730  m_h3_RecoTrack_dau1->Fill(m_MCDaughter1->getMomentum().Rho(), m_MCDaughter1->getMomentum().Theta(),
731  m_MCDaughter1->getMomentum().Phi());
732  m_h1_RecoTrack_dau1_pt->Fill(m_MCDaughter1->getMomentum().Rho());
733  m_h1_RecoTrack_dau1_pz->Fill(m_MCDaughter1->getMomentum().z());
734  m_h1_RecoTrack_dau1_p->Fill(m_MCDaughter1->getMomentum().R());
735  m_h1_RecoTrack_dau1_phi->Fill(m_MCDaughter1->getMomentum().Phi());
736  m_h1_RecoTrack_dau1_theta->Fill(m_MCDaughter1->getMomentum().Theta());
737  m_h1_RecoTrack_dau1_costheta->Fill(cos(m_MCDaughter1->getMomentum().Theta()));
738  m_h1_RecoTrack_dau1_Mother_cosAngle->Fill(mcParticle.getMomentum().Dot(m_MCDaughter1->getMomentum()) /
739  mcParticle.getMomentum().R() /
740  m_MCDaughter1->getMomentum().R());
741 
742  m_h1_RecoTrack_dau1_thetaMother->Fill(MC_theta);
743  m_h1_RecoTrack_dau1_ptMother->Fill(MC_pt);
744 
745  m_h1_RecoTrack_dau1_phiMother_total->Fill(MC_phi);
746  m_h2_RecoTrack_dau1_2D->Fill(m_MCDaughter1->getMomentum().Theta(), m_MCDaughter1->getMomentum().Rho());
747  m_h2_RecoTrack_dau1_2DMother->Fill(MC_theta, MC_pt);
748  m_h2_RecoTrack_dau1_pVScostheta->Fill(cos(m_MCDaughter1->getMomentum().Theta()), m_MCDaughter1->getMomentum().R());
749 
750  if (MC_transDist < 1.) {
751  m_h2_RecoTrack_dau1_2D_BP->Fill(m_MCDaughter1->getMomentum().Theta(), m_MCDaughter1->getMomentum().Rho());
752 
753  if (m_MCDaughter1->getMomentum().Theta() > (120 * TMath::Pi() / 180.)) m_h1_RecoTrack_dau1_phi_BW->Fill(
754  m_MCDaughter1->getMomentum().Phi());
755  else if (m_MCDaughter1->getMomentum().Theta() < (30. * TMath::Pi() / 180.)) m_h1_RecoTrack_dau1_phi_FW->Fill(
756  m_MCDaughter1->getMomentum().Phi());
757  else m_h1_RecoTrack_dau1_phi_barrel->Fill(m_MCDaughter1->getMomentum().Phi());
758 
759  if (MC_theta > (120 * TMath::Pi() / 180.)) m_h1_RecoTrack_dau1_phiMother_BW->Fill(MC_phi);
760  else if (MC_theta < (30. * TMath::Pi() / 180.)) m_h1_RecoTrack_dau1_phiMother_FW->Fill(MC_phi);
761  else m_h1_RecoTrack_dau1_phiMother_barrel->Fill(MC_phi);
762  }
763  }
764 
765  //------------------------------------------------------------------//
766  // TRACKS //
767  //------------------------------------------------------------------//
768 
769  Track* Track_dau0ToMCParticle = m_MCDaughter0->getRelated<Track>();
770  Track* Track_dau1ToMCParticle = m_MCDaughter1->getRelated<Track>();
771 
772  if (Track_dau0ToMCParticle) {
773  m_h1_track_dau0_d0->Fill(mcParticleInfo_dau0.getD0());
774  m_h1_track_dau0_z0->Fill(mcParticleInfo_dau0.getZ0());
775  m_h1_track_dau0_RMother->Fill(MC_transDist);
776  m_h3_track_dau0->Fill(m_MCDaughter0->getMomentum().Rho(), m_MCDaughter0->getMomentum().Theta(), m_MCDaughter0->getMomentum().Phi());
777  m_h1_track_dau0_pt->Fill(m_MCDaughter0->getMomentum().Rho());
778  m_h1_track_dau0_pz->Fill(m_MCDaughter0->getMomentum().z());
779  m_h1_track_dau0_p->Fill(m_MCDaughter0->getMomentum().R());
780  m_h1_track_dau0_phi->Fill(m_MCDaughter0->getMomentum().Phi());
781  m_h1_track_dau0_theta->Fill(m_MCDaughter0->getMomentum().Theta());
782  m_h1_track_dau0_costheta->Fill(cos(m_MCDaughter0->getMomentum().Theta()));
783  m_h1_track_dau0_Mother_cosAngle->Fill(mcParticle.getMomentum().Dot(m_MCDaughter0->getMomentum()) / mcParticle.getMomentum().R() /
784  m_MCDaughter0->getMomentum().R());
785 
786  m_h1_track_dau0_thetaMother->Fill(MC_theta);
787  m_h1_track_dau0_ptMother->Fill(MC_pt);
788 
789  m_h1_track_dau0_phiMother_total->Fill(MC_phi);
790 
791  m_h2_track_dau0_2D->Fill(m_MCDaughter0->getMomentum().Theta(), m_MCDaughter0->getMomentum().Rho());
792  m_h2_track_dau0_2DMother->Fill(MC_theta, MC_pt);
793  m_h2_track_dau0_pVScostheta->Fill(cos(m_MCDaughter0->getMomentum().Theta()), m_MCDaughter0->getMomentum().R());
794 
795  if (MC_transDist < 1.)
796  m_h2_track_dau0_2D_BP->Fill(m_MCDaughter0->getMomentum().Theta(), m_MCDaughter0->getMomentum().Rho());
797 
798  if (MC_theta > (120 * TMath::Pi() / 180.)) { //BW
799  m_h1_track_dau0_phiMother_BW->Fill(MC_phi);
800  m_h1_track_dau0_phi_BW->Fill(m_MCDaughter0->getMomentum().Phi());
801  } else if (MC_theta < (30 * TMath::Pi() / 180.)) { //FW
802  m_h1_track_dau0_phiMother_FW->Fill(MC_phi);
803  m_h1_track_dau0_phi_FW->Fill(m_MCDaughter0->getMomentum().Phi());
804  } else { //barrel
805  m_h1_track_dau0_phiMother_barrel->Fill(MC_phi);
806  m_h1_track_dau0_phi_barrel->Fill(m_MCDaughter0->getMomentum().Phi());
807  }
808  }
809 
810  if (Track_dau1ToMCParticle) {
811  m_h1_track_dau1_d0->Fill(mcParticleInfo_dau1.getD0());
812  m_h1_track_dau1_z0->Fill(mcParticleInfo_dau1.getZ0());
813  m_h1_track_dau1_RMother->Fill(MC_transDist);
814  m_h3_track_dau1->Fill(m_MCDaughter1->getMomentum().Rho(), m_MCDaughter1->getMomentum().Theta(), m_MCDaughter1->getMomentum().Phi());
815  m_h1_track_dau1_pt->Fill(m_MCDaughter1->getMomentum().Rho());
816  m_h1_track_dau1_pz->Fill(m_MCDaughter1->getMomentum().z());
817  m_h1_track_dau1_p->Fill(m_MCDaughter1->getMomentum().R());
818  m_h1_track_dau1_phi->Fill(m_MCDaughter1->getMomentum().Phi());
819  m_h1_track_dau1_theta->Fill(m_MCDaughter1->getMomentum().Theta());
820  m_h1_track_dau1_costheta->Fill(cos(m_MCDaughter1->getMomentum().Theta()));
821  m_h1_track_dau1_Mother_cosAngle->Fill(mcParticle.getMomentum().Dot(m_MCDaughter1->getMomentum()) / mcParticle.getMomentum().R() /
822  m_MCDaughter1->getMomentum().R());
823 
824  m_h1_track_dau1_thetaMother->Fill(MC_theta);
825  m_h1_track_dau1_ptMother->Fill(MC_pt);
826 
827  m_h1_track_dau1_phiMother_total->Fill(MC_phi);
828 
829  m_h2_track_dau1_2D->Fill(m_MCDaughter1->getMomentum().Theta(), m_MCDaughter1->getMomentum().Rho());
830  m_h2_track_dau1_2DMother->Fill(MC_theta, MC_pt);
831  m_h2_track_dau1_pVScostheta->Fill(cos(m_MCDaughter1->getMomentum().Theta()), m_MCDaughter1->getMomentum().R());
832 
833  if (MC_transDist < 1.)
834  m_h2_track_dau1_2D_BP->Fill(m_MCDaughter1->getMomentum().Theta(), m_MCDaughter1->getMomentum().Rho());
835 
836  if (MC_theta > (120 * TMath::Pi() / 180.)) { //BW
837  m_h1_track_dau1_phiMother_BW->Fill(MC_phi);
838  m_h1_track_dau1_phi_BW->Fill(m_MCDaughter1->getMomentum().Phi());
839  } else if (MC_theta < (30 * TMath::Pi() / 180.)) { //FW
840  m_h1_track_dau1_phiMother_FW->Fill(MC_phi);
841  m_h1_track_dau1_phi_FW->Fill(m_MCDaughter1->getMomentum().Phi());
842  } else { //barrel
843  m_h1_track_dau1_phiMother_barrel->Fill(MC_phi);
844  m_h1_track_dau1_phi_barrel->Fill(m_MCDaughter1->getMomentum().Phi());
845  }
846  }
847 
848  int nMatchedDau = nMatchedDaughters(mcParticle);
849 
850  //V0: proceed only in case the MCParticle daughters have one associated reconstructed track:
851  if (nMatchedDau != 2)
852  continue;
853 
854  int pdgCode = mcParticle.getPDG();
855  B2DEBUG(29, "MCParticle has PDG code " << pdgCode);
856 
857  RelationVector<V0ValidationVertex> V0s_toMCParticle =
858  DataStore::getRelationsWithObj<V0ValidationVertex>(&mcParticle, m_V0sName);
859 
860  if (V0s_toMCParticle.size() > 0) {
861 
862  m_h1_V0_RMother->Fill(MC_transDist);
863  m_h3_V0->Fill(mcParticleInfo.getPt(), mcParticleInfo.getPtheta(), mcParticleInfo.getPphi());
864  m_h1_V0_pt->Fill(mcParticleInfo.getPt());
865  m_h1_V0_pz->Fill(mcParticleInfo.getPz());
866  m_h1_V0_p->Fill(mcParticleInfo.getP());
867  m_h1_V0_phi->Fill(mcParticleInfo.getPphi());
868  m_h1_V0_theta->Fill(mcParticleInfo.getPtheta());
869  m_h1_V0_costheta->Fill(cos(mcParticle.getMomentum().Theta()));
870  m_h2_V0_Mother_2D->Fill(MC_theta, MC_p);
871  m_h2_V0_Mother_pVScostheta->Fill(MC_costheta, MC_p);
872 
873  if (MC_transDist < 1.) m_h2_V0_Mother_2D_BP->Fill(MC_theta, MC_pt);
874 
875  if (MC_theta > (120 * TMath::Pi() / 180.)) //BW
876  m_h1_V0_phi_BW->Fill(mcParticleInfo.getPphi());
877 
878  else if (MC_theta < (30 * TMath::Pi() / 180.))//FW
879  m_h1_V0_phi_FW->Fill(mcParticleInfo.getPphi());
880 
881  else //barrel
882  m_h1_V0_phi_barrel->Fill(mcParticleInfo.getPphi());
883  }
884  }
885 }
886 
887 void EffPlotsModule::endRun()
888 {
889  double track_dau0 = m_h1_track_dau0_RMother->GetEntries();
890  double RecoTrack_dau0 = m_h1_RecoTrack_dau0_RMother->GetEntries();
891  double MC_dau0 = m_h1_MC_dau0_RMother->GetEntries();
892  double eff_dau0_noGA = track_dau0 / MC_dau0;
893  double effErr_dau0_noGA = sqrt(eff_dau0_noGA * (1 - eff_dau0_noGA)) / sqrt(MC_dau0);
894  double eff_dau0_withGA = track_dau0 / RecoTrack_dau0;
895  double effErr_dau0_withGA = sqrt(eff_dau0_withGA * (1 - eff_dau0_withGA)) / sqrt(RecoTrack_dau0);
896 
897  double track_dau1 = m_h1_track_dau1_RMother->GetEntries();
898  double RecoTrack_dau1 = m_h1_RecoTrack_dau1_RMother->GetEntries();
899  double MC_dau1 = m_h1_MC_dau1_RMother->GetEntries();
900  double eff_dau1_noGA = track_dau1 / MC_dau1;
901  double effErr_dau1_noGA = sqrt(eff_dau1_noGA * (1 - eff_dau1_noGA)) / sqrt(MC_dau1);
902  double eff_dau1_withGA = track_dau1 / RecoTrack_dau1;
903  double effErr_dau1_withGA = sqrt(eff_dau1_withGA * (1 - eff_dau1_withGA)) / sqrt(RecoTrack_dau1);
904 
905  double track_Mother = m_h1_V0_RMother->GetEntries();
906  double RecoTrack_Mother = m_h1_RecoTrack_Mother_RMother->GetEntries();
907  double MC_Mother = m_h1_MC_Mother_RMother->GetEntries();
908  double eff_Mother_noGA = track_Mother / MC_Mother;
909  double effErr_Mother_noGA = sqrt(eff_Mother_noGA * (1 - eff_Mother_noGA)) / sqrt(MC_Mother);
910  double eff_Mother_withGA = track_Mother / RecoTrack_Mother;
911  double effErr_Mother_withGA = sqrt(eff_Mother_withGA * (1 - eff_Mother_withGA)) / sqrt(RecoTrack_Mother);
912 
913 
914 
915 
916  B2INFO("");
917  B2INFO("~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~");
918  B2INFO("~ V0 Finding Performance Evaluation ~ SHORT SUMMARY ~");
919  B2INFO("");
920  B2INFO(" + overall, normalized to MC particles (_noGeoAcc):");
921  B2INFO("");
922  B2INFO(" efficiency dau0 = (" << eff_dau0_noGA * 100 << " +/- " << effErr_dau0_noGA * 100 << ")% ");
923  B2INFO(" efficiency dau1 = (" << eff_dau1_noGA * 100 << " +/- " << effErr_dau1_noGA * 100 << ")% ");
924  B2INFO(" efficiency Mother = (" << eff_Mother_noGA * 100 << " +/- " << effErr_Mother_noGA * 100 << ")% ");
925  B2INFO("");
926  B2INFO(" + overall, normalized to RecoTracks(_withGeoAcc):");
927  B2INFO("");
928  B2INFO(" efficiency dau0 = (" << eff_dau0_withGA * 100 << " +/- " << effErr_dau0_withGA * 100 << ")% ");
929  B2INFO(" efficiency dau1 = (" << eff_dau1_withGA * 100 << " +/- " << effErr_dau1_withGA * 100 << ")% ");
930  B2INFO(" efficiency Mother = (" << eff_Mother_withGA * 100 << " +/- " << effErr_Mother_withGA * 100 << ")% ");
931  B2INFO("~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~");
932 
933 }
934 
935 void EffPlotsModule::terminate()
936 {
937  //------------------------------------------------------------------//
938  // EFFICIENCIES //
939  //------------------------------------------------------------------//
940 
941  //dau0
942  TH1F* h_eff_dau0_d0 = effPlot1D(m_h1_MC_dau0_d0, m_h1_RecoTrack_dau0_d0, m_h1_track_dau0_d0, "h_eff_dau0_d0",
943  "efficiency VS d0, dau_{0}", m_histoList_Efficiencies);
944  if (h_eff_dau0_d0->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_d0");
945 
946  TH1F* h_eff_dau0_z0 = effPlot1D(m_h1_MC_dau0_z0, m_h1_RecoTrack_dau0_z0, m_h1_track_dau0_z0, "h_eff_dau0_z0",
947  "efficiency VS z0, dau_{0}", m_histoList_Efficiencies);
948  if (h_eff_dau0_z0->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_z0");
949 
950  TH1F* h_eff_dau0_RMother = effPlot1D(m_h1_MC_dau0_RMother, m_h1_RecoTrack_dau0_RMother, m_h1_track_dau0_RMother,
951  "h_eff_dau0_RMother", "efficiency VS R_{mother}, dau_{0}", m_histoList_Efficiencies);
952  if (h_eff_dau0_RMother->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_RMother");
953 
954  TH1F* h_eff_dau0_thetaMother = effPlot1D(m_h1_MC_dau0_thetaMother, m_h1_RecoTrack_dau0_thetaMother, m_h1_track_dau0_thetaMother,
955  "h_eff_dau0_thetaMother", "efficiency VS #theta_{mother}, dau_{0}", m_histoList_Efficiencies);
956  if (h_eff_dau0_thetaMother->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_thetaMother");
957 
958  TH1F* h_eff_dau0_phiMother_total = effPlot1D(m_h1_MC_dau0_phiMother_total, m_h1_RecoTrack_dau0_phiMother_total,
959  m_h1_track_dau0_phiMother_total, "h_eff_dau0_phiMother_total", "efficiency VS #phi_{mother}, dau_{0}", m_histoList_Efficiencies);
960  if (h_eff_dau0_phiMother_total->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_phiMother_total");
961 
962  TH1F* h_eff_dau0_phiMother_BW = effPlot1D(m_h1_MC_dau0_phiMother_BW, m_h1_RecoTrack_dau0_phiMother_BW, m_h1_track_dau0_phiMother_BW,
963  "h_eff_dau0_phiMother_BW", "efficiency VS #phi_{mother}, dau_{0} BW", m_histoList_Efficiencies);
964  if (h_eff_dau0_phiMother_BW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_phiMother_BW");
965 
966  TH1F* h_eff_dau0_phiMother_barrel = effPlot1D(m_h1_MC_dau0_phiMother_barrel, m_h1_RecoTrack_dau0_phiMother_barrel,
967  m_h1_track_dau0_phiMother_barrel, "h_eff_dau0_phiMother_barrel", "efficiency VS #phi_{mother}, dau_{0} barrel",
968  m_histoList_Efficiencies);
969  if (h_eff_dau0_phiMother_barrel->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_phiMother_barrel");
970 
971  TH1F* h_eff_dau0_phiMother_FW = effPlot1D(m_h1_MC_dau0_phiMother_FW, m_h1_RecoTrack_dau0_phiMother_FW, m_h1_track_dau0_phiMother_FW,
972  "h_eff_dau0_phiMother_FW", "efficiency VS #phi_{mother}, dau_{0} FW", m_histoList_Efficiencies);
973  if (h_eff_dau0_phiMother_FW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_phiMother_FW");
974 
975  TH1F* h_eff_dau0_phi_total = effPlot1D(m_h1_MC_dau0_phi, m_h1_RecoTrack_dau0_phi, m_h1_track_dau0_phi, "h_eff_dau0_phi_total",
976  "efficiency VS #phi_{dau_{0}}, dau_{0}", m_histoList_Efficiencies);
977  if (h_eff_dau0_phi_total->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_phi_total");
978 
979  TH1F* h_eff_dau0_phi_BW = effPlot1D(m_h1_MC_dau0_phi_BW, m_h1_RecoTrack_dau0_phi_BW, m_h1_track_dau0_phi_BW, "h_eff_dau0_phi_BW",
980  "efficiency VS #phi_{dau_{0}}, dau_{0} BW", m_histoList_Efficiencies);
981  if (h_eff_dau0_phi_BW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_phi_BW");
982 
983  TH1F* h_eff_dau0_phi_barrel = effPlot1D(m_h1_MC_dau0_phi_barrel, m_h1_RecoTrack_dau0_phi_barrel, m_h1_track_dau0_phi_barrel,
984  "h_eff_dau0_phi_barrel", "efficiency VS #phi_{dau_{0}}, dau_{0} barrel", m_histoList_Efficiencies);
985  if (h_eff_dau0_phi_barrel->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_phi_barrel");
986 
987  TH1F* h_eff_dau0_phi_FW = effPlot1D(m_h1_MC_dau0_phi_FW, m_h1_RecoTrack_dau0_phi_FW, m_h1_track_dau0_phi_FW, "h_eff_dau0_phi_FW",
988  "efficiency VS #phi_{dau_{0}}, dau_{0} FW", m_histoList_Efficiencies);
989  if (h_eff_dau0_phi_FW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau0_phi_FW");
990 
991  TH2F* h2_effMap_dau0_pVScostheta = effPlot2D(m_h2_MC_dau0_pVScostheta, m_h2_RecoTrack_dau0_pVScostheta, m_h2_track_dau0_pVScostheta,
992  "h2_effMap_dau0_pVScostheta", "efficiency map, p_{CM} VS cos(#theta)_{CM} dau_{0}", m_histoList_Efficiencies);
993  if (h2_effMap_dau0_pVScostheta->GetEntries() == 0) B2WARNING("Empty histogram h2_effMap_dau0_pVScostheta");
994 
995  //dau1
996  TH1F* h_eff_dau1_d0 = effPlot1D(m_h1_MC_dau1_d0, m_h1_RecoTrack_dau1_d0, m_h1_track_dau1_d0, "h_eff_dau1_d0", "efficiency VS d0, p",
997  m_histoList_Efficiencies);
998  if (h_eff_dau1_d0->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_d0");
999 
1000  TH1F* h_eff_dau1_z0 = effPlot1D(m_h1_MC_dau1_z0, m_h1_RecoTrack_dau1_z0, m_h1_track_dau1_z0, "h_eff_dau1_z0", "efficiency VS z0, p",
1001  m_histoList_Efficiencies);
1002  if (h_eff_dau1_z0->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_z0");
1003 
1004  TH1F* h_eff_dau1_RMother = effPlot1D(m_h1_MC_dau1_RMother, m_h1_RecoTrack_dau1_RMother, m_h1_track_dau1_RMother,
1005  "h_eff_dau1_RMother", "efficiency VS R_{mother}, p", m_histoList_Efficiencies);
1006  if (h_eff_dau1_RMother->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_RMother");
1007 
1008  TH1F* h_eff_dau1_thetaMother = effPlot1D(m_h1_MC_dau1_thetaMother, m_h1_RecoTrack_dau1_thetaMother, m_h1_track_dau1_thetaMother,
1009  "h_eff_dau1_thetaMother", "efficiency VS #theta_{mother}, p", m_histoList_Efficiencies);
1010  if (h_eff_dau1_thetaMother->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_thetaMother");
1011 
1012  TH1F* h_eff_dau1_phiMother_total = effPlot1D(m_h1_MC_dau1_phiMother_total, m_h1_RecoTrack_dau1_phiMother_total,
1013  m_h1_track_dau1_phiMother_total, "h_eff_dau1_phiMother_total", "efficiency VS #phi_{mother}, p", m_histoList_Efficiencies);
1014  if (h_eff_dau1_phiMother_total->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_phiMother_total");
1015 
1016  TH1F* h_eff_dau1_phiMother_BW = effPlot1D(m_h1_MC_dau1_phiMother_BW, m_h1_RecoTrack_dau1_phiMother_BW, m_h1_track_dau1_phiMother_BW,
1017  "h_eff_dau1_phiMother_BW", "efficiency VS #phi_{mother}, p BW", m_histoList_Efficiencies);
1018  if (h_eff_dau1_phiMother_BW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_phiMother_BW");
1019 
1020  TH1F* h_eff_dau1_phiMother_barrel = effPlot1D(m_h1_MC_dau1_phiMother_barrel, m_h1_RecoTrack_dau1_phiMother_barrel,
1021  m_h1_track_dau1_phiMother_barrel, "h_eff_dau1_phiMother_barrel", "efficiency VS #phi_{mother}, p barrel", m_histoList_Efficiencies);
1022  if (h_eff_dau1_phiMother_barrel->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_phiMother_barrel");
1023 
1024  TH1F* h_eff_dau1_phiMother_FW = effPlot1D(m_h1_MC_dau1_phiMother_FW, m_h1_RecoTrack_dau1_phiMother_FW, m_h1_track_dau1_phiMother_FW,
1025  "h_eff_dau1_phiMother_FW", "efficiency VS #phi_{mother}, p FW", m_histoList_Efficiencies);
1026  if (h_eff_dau1_phiMother_FW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_phiMother_FW");
1027 
1028  TH1F* h_eff_dau1_phi_total = effPlot1D(m_h1_MC_dau1_phi, m_h1_RecoTrack_dau1_phi, m_h1_track_dau1_phi, "h_eff_dau1_phi_total",
1029  "efficiency VS #phi_{p}, p", m_histoList_Efficiencies);
1030  if (h_eff_dau1_phi_total->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_phi_total");
1031 
1032  TH1F* h_eff_dau1_phi_BW = effPlot1D(m_h1_MC_dau1_phi_BW, m_h1_RecoTrack_dau1_phi_BW, m_h1_track_dau1_phi_BW, "h_eff_dau1_phi_BW",
1033  "efficiency VS #phi_{p}, p BW", m_histoList_Efficiencies);
1034  if (h_eff_dau1_phi_BW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_phi_BW");
1035 
1036  TH1F* h_eff_dau1_phi_barrel = effPlot1D(m_h1_MC_dau1_phi_barrel, m_h1_RecoTrack_dau1_phi_barrel, m_h1_track_dau1_phi_barrel,
1037  "h_eff_dau1_phi_barrel", "efficiency VS #phi_{p}, p barrel", m_histoList_Efficiencies);
1038  if (h_eff_dau1_phi_barrel->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_phi_barrel");
1039 
1040  TH1F* h_eff_dau1_phi_FW = effPlot1D(m_h1_MC_dau1_phi_FW, m_h1_RecoTrack_dau1_phi_FW, m_h1_track_dau1_phi_FW, "h_eff_dau1_phi_FW",
1041  "efficiency VS #phi_{p}, p FW", m_histoList_Efficiencies);
1042  if (h_eff_dau1_phi_FW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_dau1_phi_FW");
1043 
1044  TH2F* h2_effMap_dau1_pVScostheta = effPlot2D(m_h2_MC_dau1_pVScostheta, m_h2_RecoTrack_dau1_pVScostheta, m_h2_track_dau1_pVScostheta,
1045  "h2_effMap_dau1_pVScostheta", "efficiency map, p_{CM} VS cos(#theta)_{CM} dau_{1}", m_histoList_Efficiencies);
1046  if (h2_effMap_dau1_pVScostheta->GetEntries() == 0) B2WARNING("Empty histogram h2_effMap_dau1_pVScostheta");
1047 
1048  //mother
1049  TH1F* h_eff_Mother_RMother = effPlot1D(m_h1_MC_Mother_RMother, m_h1_RecoTrack_Mother_RMother, m_h1_V0_RMother,
1050  "h_eff_Mother_RMother", "efficiency VS R_{mother}, mother", m_histoList_Efficiencies);
1051  if (h_eff_Mother_RMother->GetEntries() == 0) B2WARNING("Empty histogram h_eff_Mother_RMother");
1052 
1053  TH1F* h_eff_Mother_theta = effPlot1D(m_h1_MC_Mother_theta, m_h1_RecoTrack_Mother_theta, m_h1_V0_theta, "h_eff_Mother_theta",
1054  "efficiency VS #theta_{mother}, mother", m_histoList_Efficiencies);
1055  if (h_eff_Mother_theta->GetEntries() == 0) B2WARNING("Empty histogram h_eff_Mother_theta");
1056 
1057  TH1F* h_eff_Mother_phi_total = effPlot1D(m_h1_MC_Mother_phi, m_h1_RecoTrack_Mother_phi, m_h1_V0_phi, "h_eff_Mother_phi_total",
1058  "efficiency VS #phi_{mother}, mother", m_histoList_Efficiencies);
1059  if (h_eff_Mother_phi_total->GetEntries() == 0) B2WARNING("Empty histogram h_eff_Mother_phi_total");
1060 
1061  TH1F* h_eff_Mother_phiMother_BW = effPlot1D(m_h1_MC_Mother_phi_BW, m_h1_RecoTrack_Mother_phi_BW, m_h1_V0_phi_BW,
1062  "h_eff_Mother_phi_BW", "efficiency VS #phi_{mother}, mother BW", m_histoList_Efficiencies);
1063  if (h_eff_Mother_phiMother_BW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_Mother_phiMother_BW");
1064 
1065  TH1F* h_eff_Mother_phiMother_barrel = effPlot1D(m_h1_MC_Mother_phi_barrel, m_h1_RecoTrack_Mother_phi_barrel, m_h1_V0_phi_barrel,
1066  "h_eff_Mother_phi_barrel", "efficiency VS #phi_{mother}, mother barrel", m_histoList_Efficiencies);
1067  if (h_eff_Mother_phiMother_barrel->GetEntries() == 0) B2WARNING("Empty histogram h_eff_Mother_phiMother_barrel");
1068 
1069  TH1F* h_eff_Mother_phiMother_FW = effPlot1D(m_h1_MC_Mother_phi_FW, m_h1_RecoTrack_Mother_phi_FW, m_h1_V0_phi_FW,
1070  "h_eff_Mother_phi_FW", "efficiency VS #phi_{mother}, mother FW", m_histoList_Efficiencies);
1071  if (h_eff_Mother_phiMother_FW->GetEntries() == 0) B2WARNING("Empty histogram h_eff_Mother_phiMother_FW");
1072 
1073  TH2F* h2_effMap_Mother_pVScostheta = effPlot2D(m_h2_MC_Mother_pVScostheta, m_h2_RecoTrack_Mother_pVScostheta,
1074  m_h2_V0_Mother_pVScostheta, "h2_effMap_Mother_pVScostheta", "efficiency map, p_{CM} VS cos(#theta)_{CM} mother",
1075  m_histoList_Efficiencies);
1076  if (h2_effMap_Mother_pVScostheta->GetEntries() == 0) B2WARNING("Empty histogram h2_effMap_Mother_pVScostheta");
1077 
1078  //eff_V0Finder [ = eff_mother / (eff_dau0 * eff_dau1)]
1079  TH1F* h_eff_V0Finder_RMother_withGeoAcc = V0FinderEff(h_eff_dau0_RMother, h_eff_dau1_RMother, h_eff_Mother_RMother,
1080  "h_eff_V0Finder_RMother_withGeoAcc", "efficiency VS R_{mother}, V0Finder", m_histoList_Efficiencies);
1081  if (h_eff_V0Finder_RMother_withGeoAcc->GetEntries() == 0) B2WARNING("Empty histogram h_eff_V0Finder_RMother_withGeoAcc");
1082 
1083  TH1F* h_eff_V0Finder_theta_withGeoAcc = V0FinderEff(h_eff_dau0_thetaMother, h_eff_dau1_thetaMother, h_eff_Mother_theta,
1084  "h_eff_V0Finder_theta_withGeoAcc", "efficiency VS #theta_{mother}, V0Finder", m_histoList_Efficiencies);
1085  if (h_eff_V0Finder_theta_withGeoAcc->GetEntries() == 0) B2WARNING("Empty histogram h_eff_V0Finder_theta_withGeoAcc");
1086 
1087  TH1F* h_eff_V0Finder_phi_total_withGeoAcc = V0FinderEff(h_eff_dau0_phiMother_total, h_eff_dau1_phiMother_total,
1088  h_eff_Mother_phi_total, "h_eff_V0Finder_phi_total_withGeoAcc", "efficiency VS #phi_{mother}, V0Finder", m_histoList_Efficiencies);
1089  if (h_eff_V0Finder_phi_total_withGeoAcc->GetEntries() == 0) B2WARNING("Empty histogram h_eff_V0Finder_phi_total_withGeoAcc");
1090 
1091  TH1F* h_eff_V0Finder_phiMother_BW_withGeoAcc = V0FinderEff(h_eff_dau0_phiMother_BW, h_eff_dau1_phiMother_BW,
1092  h_eff_Mother_phiMother_BW, "h_eff_V0Finder_phi_BW_withGeoAcc", "efficiency VS #phi_{mother}, V0Finder BW",
1093  m_histoList_Efficiencies);
1094  if (h_eff_V0Finder_phiMother_BW_withGeoAcc->GetEntries() == 0) B2WARNING("Empty histogram h_eff_V0Finder_phiMother_BW_withGeoAcc");
1095 
1096  TH1F* h_eff_V0Finder_phiMother_barrel_withGeoAcc = V0FinderEff(h_eff_dau0_phiMother_barrel, h_eff_dau1_phiMother_barrel,
1097  h_eff_Mother_phiMother_barrel, "h_eff_V0Finder_phi_barrel_withGeoAcc", "efficiency VS #phi_{mother}, V0Finder barrel",
1098  m_histoList_Efficiencies);
1099  if (h_eff_V0Finder_phiMother_barrel_withGeoAcc->GetEntries() == 0)
1100  B2WARNING("Empty histogram h_eff_V0Finder_phiMother_barrel_withGeoAcc");
1101 
1102  TH1F* h_eff_V0Finder_phiMother_FW_withGeoAcc = V0FinderEff(h_eff_dau0_phiMother_FW, h_eff_dau1_phiMother_FW,
1103  h_eff_Mother_phiMother_FW, "h_eff_V0Finder_phi_FW_withGeoAcc", "efficiency VS #phi_{mother}, V0Finder FW",
1104  m_histoList_Efficiencies);
1105  if (h_eff_V0Finder_phiMother_FW_withGeoAcc->GetEntries() == 0) B2WARNING("Empty histogram h_eff_V0Finder_phiMother_FW_withGeoAcc");
1106 
1107 
1108  //-------------------------------------------------------------------------//
1109  // GEOMETRICAL ACCETTANCE //
1110  //-------------------------------------------------------------------------//
1111 
1112  if (m_geometricalAccettance) {
1113  //dau0
1114  TH1F* h_dau0_geoAcc_theta = geoAcc1D(m_h1_MC_dau0_theta, m_h1_RecoTrack_dau0_theta, "h_dau0_geoAcc_theta",
1115  "geometrical acceptance VS #theta_{dau_{0}}, dau_{0}", m_histoList_GA);
1116  if (h_dau0_geoAcc_theta->GetEntries() == 0) B2WARNING("Empty histogram h_dau0_geoAcc_theta");
1117 
1118  TH1F* h_dau0_geoAcc_phi = geoAcc1D(m_h1_MC_dau0_phi, m_h1_RecoTrack_dau0_phi, "h_dau0_geoAcc_phi",
1119  "geometrical acceptance VS #phi_{dau_{0}}, dau_{0}", m_histoList_GA);
1120  if (h_dau0_geoAcc_phi->GetEntries() == 0) B2WARNING("Empty histogram h_dau0_geoAcc_phi");
1121 
1122  TH1F* h_dau0_geoAcc_pt = geoAcc1D(m_h1_MC_dau0_pt, m_h1_RecoTrack_dau0_pt, "h_dau0_geoAcc_pt",
1123  "geometrical acceptance VS p_{T,dau_{0}}, dau_{0}", m_histoList_GA);
1124  if (h_dau0_geoAcc_pt->GetEntries() == 0) B2WARNING("Empty histogram h_dau0_geoAcc_pt");
1125 
1126  TH1F* h_dau0_geoAcc_thetaMother = geoAcc1D(m_h1_MC_dau0_thetaMother, m_h1_RecoTrack_dau0_thetaMother, "h_dau0_geoAcc_thetaMother",
1127  "geometrical acceptance VS #theta_{mother}, dau_{0}", m_histoList_GA);
1128  if (h_dau0_geoAcc_thetaMother->GetEntries() == 0) B2WARNING("Empty histogram h_dau0_geoAcc_thetaMother");
1129 
1130  TH1F* h_dau0_geoAcc_phiMother = geoAcc1D(m_h1_MC_dau0_phiMother_total, m_h1_RecoTrack_dau0_phiMother_total,
1131  "h_dau0_geoAcc_phiMother", "geometrical acceptance VS #phi_{mother}, dau_{0}", m_histoList_GA);
1132  if (h_dau0_geoAcc_phiMother->GetEntries() == 0) B2WARNING("Empty histogram h_dau0_geoAcc_phiMother");
1133 
1134  TH1F* h_dau0_geoAcc_ptMother = geoAcc1D(m_h1_MC_dau0_ptMother, m_h1_RecoTrack_dau0_ptMother, "h_dau0_geoAcc_ptMother",
1135  "geometrical acceptance VS #p_{T,mother}, dau_{0}", m_histoList_GA);
1136  if (h_dau0_geoAcc_ptMother->GetEntries() == 0) B2WARNING("Empty histogram h_dau0_geoAcc_ptMother");
1137 
1138  //dau1
1139  TH1F* h_dau1_geoAcc_theta = geoAcc1D(m_h1_MC_dau1_theta, m_h1_RecoTrack_dau1_theta, "h_dau1_geoAcc_theta",
1140  "geometrical acceptance VS #theta_{p}, p", m_histoList_GA);
1141  if (h_dau1_geoAcc_theta->GetEntries() == 0) B2WARNING("Empty histogram h_dau1_geoAcc_theta");
1142 
1143  TH1F* h_dau1_geoAcc_phi = geoAcc1D(m_h1_MC_dau1_phi, m_h1_RecoTrack_dau1_phi, "h_dau1_geoAcc_phi",
1144  "geometrical acceptance VS #phi_{p}, p", m_histoList_GA);
1145  if (h_dau1_geoAcc_phi->GetEntries() == 0) B2WARNING("Empty histogram h_dau1_geoAcc_phi");
1146 
1147  TH1F* h_dau1_geoAcc_pt = geoAcc1D(m_h1_MC_dau1_pt, m_h1_RecoTrack_dau1_pt, "h_dau1_geoAcc_pt",
1148  "geometrical acceptance VS p_{T,p}, p", m_histoList_GA);
1149  if (h_dau1_geoAcc_pt->GetEntries() == 0) B2WARNING("Empty histogram h_dau1_geoAcc_pt");
1150 
1151  TH1F* h_dau1_geoAcc_thetaMother = geoAcc1D(m_h1_MC_dau1_thetaMother, m_h1_RecoTrack_dau1_thetaMother, "h_dau1_geoAcc_thetaMother",
1152  "geometrical acceptance VS #theta_{mother}, p", m_histoList_GA);
1153  if (h_dau1_geoAcc_thetaMother->GetEntries() == 0) B2WARNING("Empty histogram h_dau1_geoAcc_thetaMother");
1154 
1155  TH1F* h_dau1_geoAcc_phiMother = geoAcc1D(m_h1_MC_dau1_phiMother_total, m_h1_RecoTrack_dau1_phiMother_total,
1156  "h_dau1_geoAcc_phiMother", "geometrical acceptance VS #phi_{mother}, p", m_histoList_GA);
1157  if (h_dau1_geoAcc_phiMother->GetEntries() == 0) B2WARNING("Empty histogram h_dau1_geoAcc_phiMother");
1158 
1159  TH1F* h_dau1_geoAcc_ptMother = geoAcc1D(m_h1_MC_dau1_ptMother, m_h1_RecoTrack_dau1_ptMother, "h_dau1_geoAcc_ptMother",
1160  "geometrical acceptance VS #p_{T,mother}, p", m_histoList_GA);
1161  if (h_dau1_geoAcc_ptMother->GetEntries() == 0) B2WARNING("Empty histogram h_dau1_geoAcc_ptMother");
1162 
1163  //mother
1164  TH1F* h_Mother_geoAcc_theta = geoAcc1D(m_h1_MC_Mother_theta, m_h1_RecoTrack_Mother_theta, "h_Mother_geoAcc_theta",
1165  "geometrical acceptance VS #theta_{mother}, mother", m_histoList_GA);
1166  if (h_Mother_geoAcc_theta->GetEntries() == 0) B2WARNING("Empty histogram h_Mother_geoAcc_theta");
1167 
1168  TH1F* h_Mother_geoAcc_phi = geoAcc1D(m_h1_MC_Mother_phi, m_h1_RecoTrack_Mother_phi, "h_Mother_geoAcc_phi",
1169  "geometrical acceptance VS #phi_{mother}, mother", m_histoList_GA);
1170  if (h_Mother_geoAcc_phi->GetEntries() == 0) B2WARNING("Empty histogram h_Mother_geoAcc_phi");
1171 
1172  TH1F* h_Mother_geoAcc_pt = geoAcc1D(m_h1_MC_Mother_pt, m_h1_RecoTrack_Mother_pt, "h_Mother_geoAcc_pt",
1173  "geometrical acceptance VS p_{T,mother}, mother", m_histoList_GA);
1174  if (h_Mother_geoAcc_pt->GetEntries() == 0) B2WARNING("Empty histogram h_Mother_geoAcc_pt");
1175 
1176  //2D
1177  TH2F* h2_dau0_geoAcc2D = geoAcc2D(m_h2_MC_dau0_2D, m_h2_RecoTrack_dau0_2D, "h2_dau0_geoAcc2D",
1178  "geometrical acceptance dau_{0}, p_{T,dau_{0}} VS #theta_{dau_{0}}", m_histoList_GA);
1179  if (h2_dau0_geoAcc2D->GetEntries() == 0) B2WARNING("Empty histogram h2_dau0_geoAcc2D");
1180 
1181  TH2F* h2_dau1_geoAcc2D = geoAcc2D(m_h2_MC_dau1_2D, m_h2_RecoTrack_dau1_2D, "h2_dau1_geoAcc2D",
1182  "geometrical acceptance p, p_{T,p} VS #theta_{p}", m_histoList_GA);
1183  if (h2_dau1_geoAcc2D->GetEntries() == 0) B2WARNING("Empty histogram h2_dau1_geoAcc2D");
1184 
1185  TH2F* h2_Mother_geoAcc2D = geoAcc2D(m_h2_MC_Mother_2D, m_h2_RecoTrack_Mother_2D, "h2_Mother_geoAcc2D",
1186  "geometrical acceptance mother, p_{T,mother} VS #theta_{mother}", m_histoList_GA);
1187  if (h2_Mother_geoAcc2D->GetEntries() == 0) B2WARNING("Empty histogram h2_Mother_geoAcc2D");
1188 
1189  TH2F* h2_dau0Mother_geoAcc2D = geoAcc2D(m_h2_MC_dau0_2DMother, m_h2_RecoTrack_dau0_2DMother, "h2_dau0Mother_geoAcc2D",
1190  "geometrical acceptance dau_{0}, p_{T,mother} VS #theta_{mother}", m_histoList_GA);
1191  if (h2_dau0Mother_geoAcc2D->GetEntries() == 0) B2WARNING("Empty histogram h2_dau0Mother_geoAcc2D");
1192 
1193  TH2F* h2_dau1Mother_geoAcc2D = geoAcc2D(m_h2_MC_dau1_2DMother, m_h2_RecoTrack_dau1_2DMother, "h2_dau1Mother_geoAcc2D",
1194  "geometrical acceptance p, p_{T,mother} VS #theta_{mother}", m_histoList_GA);
1195  if (h2_dau1Mother_geoAcc2D->GetEntries() == 0) B2WARNING("Empty histogram h2_dau1Mother_geoAcc2D");
1196 
1197  //particles coming from inside the Beam Pipe
1198  TH2F* h2_dau0_geoAcc2D_BP = geoAcc2D(m_h2_MC_dau0_2D_BP, m_h2_RecoTrack_dau0_2D_BP, "h2_dau0_geoAcc2D_BP",
1199  "geometrical acceptance dau_{0}, p_{T,dau_{0}} VS #theta_{dau_{0}}, BP", m_histoList_GA);
1200  if (h2_dau0_geoAcc2D_BP->GetEntries() == 0) B2WARNING("Empty histogram h2_dau0_geoAcc2D_BP");
1201 
1202  TH2F* h2_dau1_geoAcc2D_BP = geoAcc2D(m_h2_MC_dau1_2D_BP, m_h2_RecoTrack_dau1_2D_BP, "h2_dau1_geoAcc2D_BP",
1203  "geometrical acceptance p, p_{T,p} VS #theta_{p}, BP", m_histoList_GA);
1204  if (h2_dau1_geoAcc2D_BP->GetEntries() == 0) B2WARNING("Empty histogram h2_dau1_geoAcc2D_BP");
1205 
1206  TH2F* h2_Mother_geoAcc2D_BP = geoAcc2D(m_h2_MC_Mother_2D_BP, m_h2_RecoTrack_Mother_2D_BP, "h2_Mother_geoAcc2D_BP",
1207  "geometrical acceptance mother, p_{T,mother} VS #theta_{mother}, BP", m_histoList_GA);
1208  if (h2_Mother_geoAcc2D_BP->GetEntries() == 0) B2WARNING("Empty histogram h2_Mother_geoAcc2D_BP");
1209  }
1210 
1211  //---------------------------------------------------------------------------------
1212 
1213  //write histograms on the output root file
1214  if (m_rootFilePtr != nullptr) {
1215  m_rootFilePtr->cd();
1216 
1217  TIter nextHMC(m_histoList_MCParticles);
1218  TIter nextHTC(m_histoList_RecoTracks);
1219  TIter nextHT(m_histoList_Tracks);
1220  TIter nextHEff(m_histoList_Efficiencies);
1221  TIter nextHGA(m_histoList_GA);
1222  TIter nextHCheck(m_histoList_check);
1223 
1224  TObject* obj;
1225 
1226  if (m_allHistograms) {
1227  TDirectory* MCdir = m_rootFilePtr->mkdir("MCParticles");
1228  MCdir->cd();
1229  while ((obj = nextHMC()))
1230  obj->Write("", BIT(2), 0);
1231 
1232  TDirectory* TCdir = m_rootFilePtr->mkdir("RecoTracks");
1233  TCdir->cd();
1234  while ((obj = nextHTC()))
1235  obj->Write("", BIT(2), 0);
1236 
1237  TDirectory* Trdir = m_rootFilePtr->mkdir("Tracks");
1238  Trdir->cd();
1239  while ((obj = nextHT()))
1240  obj->Write("", BIT(2), 0);
1241 
1242  TDirectory* check = m_rootFilePtr->mkdir("check");
1243  check->cd();
1244  while ((obj = nextHCheck()))
1245  obj->Write("", BIT(2), 0);
1246  }
1247 
1248  TDirectory* Effdir = m_rootFilePtr->mkdir("Efficiencies");
1249  Effdir->cd();
1250  while ((obj = nextHEff()))
1251  obj->Write("", BIT(2), 0);
1252 
1253  if (m_geometricalAccettance) {
1254  TDirectory* geoAcc = m_rootFilePtr->mkdir("geoAccettance");
1255  geoAcc->cd();
1256  while ((obj = nextHGA()))
1257  obj->Write("", BIT(2), 0);
1258  }
1259 
1260 
1261  m_rootFilePtr->Close();
1262  }
1263 }
1264 
1265 bool EffPlotsModule::isK_Short(const MCParticle& the_mcParticle)
1266 {
1267 
1268  bool isK_S0 = false;
1269  if (abs(the_mcParticle.getPDG()) == Const::Kshort.getPDGCode())
1270  isK_S0 = true;
1271 
1272  bool twoChargedProngs = false;
1273 
1274  if (the_mcParticle.getDaughters().size() == 2 &&
1275  ((the_mcParticle.getDaughters()[0]->getPDG() == Const::pion.getPDGCode()
1276  && the_mcParticle.getDaughters()[1]->getPDG() == -Const::pion.getPDGCode()) ||
1277  (the_mcParticle.getDaughters()[0]->getPDG() == -Const::pion.getPDGCode()
1278  && the_mcParticle.getDaughters()[1]->getPDG() == Const::pion.getPDGCode())))
1279  twoChargedProngs = true;
1280 
1281  return (isK_S0 && twoChargedProngs);
1282 
1283 }
1284 
1285 bool EffPlotsModule::isLambda0(const MCParticle& the_mcParticle)
1286 {
1287 
1288  bool isLambda = false;
1289  if (abs(the_mcParticle.getPDG()) == Const::Lambda.getPDGCode())
1290  isLambda = true;
1291 
1292  bool twoChargedProngs = false;
1293 
1294  if (the_mcParticle.getDaughters().size() == 2 &&
1295  ((the_mcParticle.getDaughters()[0]->getPDG() == Const::pion.getPDGCode()
1296  && the_mcParticle.getDaughters()[1]->getPDG() == -Const::proton.getPDGCode()) ||
1297  (the_mcParticle.getDaughters()[0]->getPDG() == -Const::pion.getPDGCode()
1298  && the_mcParticle.getDaughters()[1]->getPDG() == Const::proton.getPDGCode()) ||
1299  (the_mcParticle.getDaughters()[0]->getPDG() == Const::proton.getPDGCode()
1300  && the_mcParticle.getDaughters()[1]->getPDG() == -Const::pion.getPDGCode()) ||
1301  (the_mcParticle.getDaughters()[0]->getPDG() == -Const::proton.getPDGCode()
1302  && the_mcParticle.getDaughters()[1]->getPDG() == Const::pion.getPDGCode())))
1303  twoChargedProngs = true;
1304 
1305  return (isLambda && twoChargedProngs);
1306 
1307 }
1308 
1309 int EffPlotsModule::nMatchedDaughters(const MCParticle& the_mcParticle)
1310 {
1311 
1312  int nMatchedDau = 0;
1313 
1314  std::vector< MCParticle* > MCPart_dau = the_mcParticle.getDaughters();
1315 
1316  bool first = false;
1317  bool second = false;
1318 
1319  RelationVector<Track> Tracks_fromMCParticle_0 = DataStore::getRelationsWithObj<Track>(MCPart_dau[0]);
1320  if (Tracks_fromMCParticle_0.size() > 0)
1321  first = true;
1322 
1323  RelationVector<Track> Tracks_fromMCParticle_1 = DataStore::getRelationsWithObj<Track>(MCPart_dau[1]);
1324  if (Tracks_fromMCParticle_1.size() > 0)
1325  second = true;
1326 
1327 
1328  if (first)
1329  nMatchedDau++;
1330 
1331  if (second)
1332  nMatchedDau++;
1333 
1334 
1335  return nMatchedDau;
1336 
1337 }
Belle2::Const::ParticleType::getPDGCode
int getPDGCode() const
PDG code.
Definition: Const.h:464
Belle2::Const::Lambda
static const ParticleType Lambda
Lambda particle.
Definition: Const.h:670
Belle2::Const::pion
static const ChargedStable pion
charged pion particle
Definition: Const.h:652
Belle2::Const::proton
static const ChargedStable proton
proton particle
Definition: Const.h:654
Belle2::Const::Kshort
static const ParticleType Kshort
K^0_S particle.
Definition: Const.h:668
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_phiMother_FW
TH1F * m_h1_RecoTrack_dau0_phiMother_FW
histogram of RecoTrack daughter 0's mother's phi (forward region)
Definition: EffPlotsModule.h:264
Belle2::EffPlotsModule::m_h1_MC_dau1_Mother_cosAngle
TH1F * m_h1_MC_dau1_Mother_cosAngle
histogram of MCParticle daughter 1's and mother's cos(opening-angle)
Definition: EffPlotsModule.h:151
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_costheta
TH1F * m_h1_RecoTrack_dau1_costheta
histogram of RecoTrack daughter 1's cos(theta)
Definition: EffPlotsModule.h:284
Belle2::EffPlotsModule::m_h1_V0_pz
TH1F * m_h1_V0_pz
histogram of V0 mother's pz
Definition: EffPlotsModule.h:234
Belle2::EffPlotsModule::m_geometricalAccettance
bool m_geometricalAccettance
true to create output for the geometrical acceptance
Definition: EffPlotsModule.h:99
Belle2::EffPlotsModule::m_h1_V0_phi
TH1F * m_h1_V0_phi
histogram of V0 mother's phi
Definition: EffPlotsModule.h:236
Belle2::EffPlotsModule::m_h3_track_dau1
TH3F * m_h3_track_dau1
histogram of Track daughter 1's pt vs theta vs phi
Definition: EffPlotsModule.h:209
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_pz
TH1F * m_h1_RecoTrack_dau0_pz
histogram of RecoTrack daughter 0's pz
Definition: EffPlotsModule.h:252
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_phi_BW
TH1F * m_h1_RecoTrack_dau0_phi_BW
histogram of RecoTrack daughter 0's phi (backward region)
Definition: EffPlotsModule.h:255
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_theta
TH1F * m_h1_RecoTrack_Mother_theta
histogram of RecoTrack mother's theta
Definition: EffPlotsModule.h:306
Belle2::EffPlotsModule::m_h1_track_dau0_phi_barrel
TH1F * m_h1_track_dau0_phi_barrel
histogram of Track daughter 0's phi (barrel region)
Definition: EffPlotsModule.h:190
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_p
TH1F * m_h1_RecoTrack_dau1_p
histogram of RecoTrack daughter 1's p
Definition: EffPlotsModule.h:278
Belle2::EffPlotsModule::m_TFRColName
std::string m_TFRColName
name of the TFR dataobjects collection
Definition: EffPlotsModule.h:92
Belle2::EffPlotsModule::m_h1_track_dau1_phi_barrel
TH1F * m_h1_track_dau1_phi_barrel
histogram of Track daughter 1's phi (barrel region)
Definition: EffPlotsModule.h:215
Belle2::EffPlotsModule::m_h2_RecoTrack_Mother_pVScostheta
TH2F * m_h2_RecoTrack_Mother_pVScostheta
histogram of RecoTrack mother's p vs cos(theta)
Definition: EffPlotsModule.h:310
Belle2::EffPlotsModule::m_h1_track_dau0_ptMother
TH1F * m_h1_track_dau0_ptMother
histogram of Track daughter 0's mother's pt
Definition: EffPlotsModule.h:200
Belle2::EffPlotsModule::isK_Short
bool isK_Short(const MCParticle &the_mcParticle)
determine if the MCParticle is a K-short
Definition: EffPlotsModule.cc:1265
Belle2::EffPlotsModule::m_h1_MC_dau0_PDG
TH1F * m_h1_MC_dau0_PDG
histogram of MCParticle daughter 0's PDG code
Definition: EffPlotsModule.h:136
Belle2::EffPlotsModule::m_h1_track_dau1_z0
TH1F * m_h1_track_dau1_z0
histogram of Track daughter 1's z0
Definition: EffPlotsModule.h:207
Belle2::EffPlotsModule::m_h1_track_dau1_pt
TH1F * m_h1_track_dau1_pt
histogram of Track daughter 1's pt
Definition: EffPlotsModule.h:210
Belle2::EffPlotsModule::m_h1_track_dau1_ptMother
TH1F * m_h1_track_dau1_ptMother
histogram of Track daughter 1's mother's pt
Definition: EffPlotsModule.h:225
Belle2::EffPlotsModule::m_h1_MC_dau1_PDG
TH1F * m_h1_MC_dau1_PDG
histogram of MCParticle daughter 1's PDG code
Definition: EffPlotsModule.h:162
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_phiMother_BW
TH1F * m_h1_RecoTrack_dau0_phiMother_BW
histogram of RecoTrack daughter 0's mother's phi (backward region)
Definition: EffPlotsModule.h:262
Belle2::EffPlotsModule::m_h1_V0_pt
TH1F * m_h1_V0_pt
histogram of V0 mother's pt
Definition: EffPlotsModule.h:233
Belle2::EffPlotsModule::m_h1_track_dau0_RMother
TH1F * m_h1_track_dau0_RMother
histogram of Track daughter 0's RMother
Definition: EffPlotsModule.h:183
Belle2::EffPlotsModule::m_h2_V0_Mother_2D
TH2F * m_h2_V0_Mother_2D
histogram of V0 mother's pt vs theta
Definition: EffPlotsModule.h:242
Belle2::EffPlotsModule::m_h1_track_dau1_pz
TH1F * m_h1_track_dau1_pz
histogram of Track daughter 1's pz
Definition: EffPlotsModule.h:211
Belle2::EffPlotsModule::m_h1_track_dau0_phiMother_barrel
TH1F * m_h1_track_dau0_phiMother_barrel
histogram of Track daughter 0's mother's phi (barrel region)
Definition: EffPlotsModule.h:197
Belle2::EffPlotsModule::m_h1_track_dau0_costheta
TH1F * m_h1_track_dau0_costheta
histogram of Track daughter 0's cos(theta)
Definition: EffPlotsModule.h:193
Belle2::EffPlotsModule::m_h2_track_dau1_2DMother
TH2F * m_h2_track_dau1_2DMother
histogram of Track daughter 1's mother's pt vs theta
Definition: EffPlotsModule.h:228
Belle2::EffPlotsModule::m_h1_MC_dau0_z0
TH1F * m_h1_MC_dau0_z0
histogram of MCParticle daughter 0's z0
Definition: EffPlotsModule.h:113
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_z0
TH1F * m_h1_RecoTrack_dau1_z0
histogram of RecoTrack daughter 1's z0
Definition: EffPlotsModule.h:273
Belle2::EffPlotsModule::m_h1_MC_dau1_ptMother
TH1F * m_h1_MC_dau1_ptMother
histogram of MCParticle daughter 1's mother's pt
Definition: EffPlotsModule.h:157
Belle2::EffPlotsModule::initialize
void initialize() override
Initializer.
Definition: EffPlotsModule.cc:60
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_phiMother_total
TH1F * m_h1_RecoTrack_dau1_phiMother_total
histogram of RecoTrack daughter 1's mother's phi
Definition: EffPlotsModule.h:286
Belle2::EffPlotsModule::m_h2_MC_dau1_2D_BP
TH2F * m_h2_MC_dau1_2D_BP
histogram of MCParticle daughter 1's pt vs theta (beam pipe)
Definition: EffPlotsModule.h:159
Belle2::EffPlotsModule::m_h2_RecoTrack_Mother_2D_BP
TH2F * m_h2_RecoTrack_Mother_2D_BP
histogram of RecoTrack mother's pt vs theta (beam pipe)
Definition: EffPlotsModule.h:309
Belle2::EffPlotsModule::m_h2_MC_dau1_pVScostheta
TH2F * m_h2_MC_dau1_pVScostheta
histogram of MCParticle daughter 1's p vs cos(theta)
Definition: EffPlotsModule.h:161
Belle2::EffPlotsModule::m_h1_V0_phi_barrel
TH1F * m_h1_V0_phi_barrel
histogram of V0 mother's phi (barrel region)
Definition: EffPlotsModule.h:238
Belle2::EffPlotsModule::m_h3_RecoTrack_Mother
TH3F * m_h3_RecoTrack_Mother
histogram of RecoTrack mother's pt vs theta vs phi
Definition: EffPlotsModule.h:298
Belle2::EffPlotsModule::m_h1_track_dau0_pz
TH1F * m_h1_track_dau0_pz
histogram of Track daughter 0's pz
Definition: EffPlotsModule.h:186
Belle2::EffPlotsModule::m_h1_MC_dau1_phi_BW
TH1F * m_h1_MC_dau1_phi_BW
histogram of MCParticle daughter 1's phi (backward region)
Definition: EffPlotsModule.h:146
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_Mother_cosAngle
TH1F * m_h1_RecoTrack_dau1_Mother_cosAngle
histogram of RecoTrack daughter 1's and mother's cos(opening-angle)
Definition: EffPlotsModule.h:285
Belle2::EffPlotsModule::m_h2_track_dau0_2D
TH2F * m_h2_track_dau0_2D
histogram of Track daughter 0's pt vs theta
Definition: EffPlotsModule.h:201
Belle2::EffPlotsModule::m_h2_MC_Mother_pVScostheta
TH2F * m_h2_MC_Mother_pVScostheta
histogram of MCParticle mother's p vs cos(theta)
Definition: EffPlotsModule.h:177
Belle2::EffPlotsModule::m_h1_V0_phi_FW
TH1F * m_h1_V0_phi_FW
histogram of V0 mother's phi (forward region)
Definition: EffPlotsModule.h:239
Belle2::EffPlotsModule::m_h1_track_dau0_phiMother_total
TH1F * m_h1_track_dau0_phiMother_total
histogram of Track daughter 0's mother's phi
Definition: EffPlotsModule.h:195
Belle2::EffPlotsModule::m_h1_track_dau1_phi
TH1F * m_h1_track_dau1_phi
histogram of Track daughter 1's phi
Definition: EffPlotsModule.h:213
Belle2::EffPlotsModule::m_h1_V0_RMother
TH1F * m_h1_V0_RMother
histogram of V0 mother's RMother
Definition: EffPlotsModule.h:231
Belle2::EffPlotsModule::m_h1_track_dau0_theta
TH1F * m_h1_track_dau0_theta
histogram of Track daughter 0's theta
Definition: EffPlotsModule.h:192
Belle2::EffPlotsModule::m_h1_MC_dau0_thetaMother
TH1F * m_h1_MC_dau0_thetaMother
histogram of MCParticle daughter 0's mother's theta
Definition: EffPlotsModule.h:130
Belle2::EffPlotsModule::m_h1_MC_Mother_pt
TH1F * m_h1_MC_Mother_pt
histogram of MCParticle mother's pt
Definition: EffPlotsModule.h:166
Belle2::EffPlotsModule::event
void event() override
This method is called for each event.
Definition: EffPlotsModule.cc:512
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_phi
TH1F * m_h1_RecoTrack_dau1_phi
histogram of RecoTrack daughter 1's phi
Definition: EffPlotsModule.h:279
Belle2::EffPlotsModule::m_h2_V0_Mother_2D_BP
TH2F * m_h2_V0_Mother_2D_BP
histogram of V0 mother's pt vs theta (beam pipe)
Definition: EffPlotsModule.h:243
Belle2::EffPlotsModule::m_h3_RecoTrack_dau1
TH3F * m_h3_RecoTrack_dau1
histogram of RecoTrack daughter 1's pt vs theta vs phi
Definition: EffPlotsModule.h:275
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_phi_FW
TH1F * m_h1_RecoTrack_dau1_phi_FW
histogram of RecoTrack daughter 1's phi (forward region)
Definition: EffPlotsModule.h:282
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_d0
TH1F * m_h1_RecoTrack_dau0_d0
list of histograms filled per RecoTracks found in the event
Definition: EffPlotsModule.h:247
Belle2::EffPlotsModule::m_h1_MC_Mother_phi
TH1F * m_h1_MC_Mother_phi
histogram of MCParticle mother's phi
Definition: EffPlotsModule.h:169
Belle2::EffPlotsModule::m_h1_V0_costheta
TH1F * m_h1_V0_costheta
histogram of V0 mother's cos(theta)
Definition: EffPlotsModule.h:241
Belle2::EffPlotsModule::m_h1_MC_dau1_RMother
TH1F * m_h1_MC_dau1_RMother
histogram of MCParticle daughter 1's RMother
Definition: EffPlotsModule.h:140
Belle2::EffPlotsModule::m_h1_track_dau1_Mother_cosAngle
TH1F * m_h1_track_dau1_Mother_cosAngle
histogram of Track daughter 1's and mother's cos(opening-angle)
Definition: EffPlotsModule.h:219
Belle2::EffPlotsModule::m_h1_MC_Mother_phi_FW
TH1F * m_h1_MC_Mother_phi_FW
histogram of MCParticle mother's phi (forward region)
Definition: EffPlotsModule.h:172
Belle2::EffPlotsModule::m_h2_RecoTrack_dau1_pVScostheta
TH2F * m_h2_RecoTrack_dau1_pVScostheta
histogram of RecoTrack daughter 1's p vs cos(theta)
Definition: EffPlotsModule.h:295
Belle2::EffPlotsModule::m_h1_MC_dau0_phi
TH1F * m_h1_MC_dau0_phi
histogram of MCParticle daughter 0's phi
Definition: EffPlotsModule.h:119
Belle2::EffPlotsModule::m_allHistograms
bool m_allHistograms
true to create all histograms
Definition: EffPlotsModule.h:98
Belle2::EffPlotsModule::endRun
void endRun() override
This method is called if the current run ends.
Definition: EffPlotsModule.cc:887
Belle2::EffPlotsModule::m_h2_track_dau0_pVScostheta
TH2F * m_h2_track_dau0_pVScostheta
histogram of Track daughter 0's p vs cos(theta)
Definition: EffPlotsModule.h:204
Belle2::EffPlotsModule::m_h2_RecoTrack_dau0_2D_BP
TH2F * m_h2_RecoTrack_dau0_2D_BP
histogram of RecoTrack daughter 0's pt vs theta (beam pipe)
Definition: EffPlotsModule.h:268
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_phi_BW
TH1F * m_h1_RecoTrack_dau1_phi_BW
histogram of RecoTrack daughter 1's phi (backward region)
Definition: EffPlotsModule.h:280
Belle2::EffPlotsModule::m_h1_MC_dau1_thetaMother
TH1F * m_h1_MC_dau1_thetaMother
histogram of MCParticle daughter 1's mother's theta
Definition: EffPlotsModule.h:156
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_phiMother_total
TH1F * m_h1_RecoTrack_dau0_phiMother_total
histogram of RecoTrack daughter 0's mother's phi
Definition: EffPlotsModule.h:261
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_phi_barrel
TH1F * m_h1_RecoTrack_Mother_phi_barrel
histogram of RecoTrack mother's phi (barrel region)
Definition: EffPlotsModule.h:304
Belle2::EffPlotsModule::m_h1_MC_dau0_pt
TH1F * m_h1_MC_dau0_pt
histogram of MCParticle daughter 0's pt
Definition: EffPlotsModule.h:116
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_ptMother
TH1F * m_h1_RecoTrack_dau0_ptMother
histogram of RecoTrack daughter 0's mother's pt
Definition: EffPlotsModule.h:266
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_phi_FW
TH1F * m_h1_RecoTrack_dau0_phi_FW
histogram of RecoTrack daughter 0's phi (forward region)
Definition: EffPlotsModule.h:257
Belle2::EffPlotsModule::terminate
void terminate() override
This method is called at the end of the event processing.
Definition: EffPlotsModule.cc:935
Belle2::EffPlotsModule::m_h1_MC_dau0_costheta
TH1F * m_h1_MC_dau0_costheta
histogram of MCParticle daughter 0's cos(theta)
Definition: EffPlotsModule.h:124
Belle2::EffPlotsModule::m_histoList_check
TList * m_histoList_check
list of histograms of MCParticle mother and daughters
Definition: EffPlotsModule.h:109
Belle2::EffPlotsModule::m_h1_MC_dau1_theta
TH1F * m_h1_MC_dau1_theta
histogram of MCParticle daughter 1's theta
Definition: EffPlotsModule.h:149
Belle2::EffPlotsModule::m_h1_track_dau1_p
TH1F * m_h1_track_dau1_p
histogram of Track daughter 1's p
Definition: EffPlotsModule.h:212
Belle2::EffPlotsModule::m_h1_MC_Mother_PDG
TH1F * m_h1_MC_Mother_PDG
histogram of MCParticle mother's PDG code
Definition: EffPlotsModule.h:178
Belle2::EffPlotsModule::m_h1_MC_dau0_Mother_cosAngle
TH1F * m_h1_MC_dau0_Mother_cosAngle
histogram of MCParticle daughter 0's and mother's cos(opening-angle)
Definition: EffPlotsModule.h:125
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_theta
TH1F * m_h1_RecoTrack_dau0_theta
histogram of RecoTrack daughter 0's theta
Definition: EffPlotsModule.h:258
Belle2::EffPlotsModule::m_TrackColName
std::string m_TrackColName
name of the Tracks dataobjects collection
Definition: EffPlotsModule.h:93
Belle2::EffPlotsModule::m_h3_MC_dau0
TH3F * m_h3_MC_dau0
histogram of MCParticle daughter 0's pt vs theta vs phi
Definition: EffPlotsModule.h:115
Belle2::EffPlotsModule::m_h1_track_dau0_pt
TH1F * m_h1_track_dau0_pt
histogram of Track daughter 0's pt
Definition: EffPlotsModule.h:185
Belle2::EffPlotsModule::m_h2_MC_Mother_2D_BP
TH2F * m_h2_MC_Mother_2D_BP
histogram of MCParticle mother's pt vs theta (beam pipe)
Definition: EffPlotsModule.h:176
Belle2::EffPlotsModule::m_h1_MC_dau1_p
TH1F * m_h1_MC_dau1_p
histogram of MCParticle daughter 1's p
Definition: EffPlotsModule.h:144
Belle2::EffPlotsModule::m_h1_track_dau0_Mother_cosAngle
TH1F * m_h1_track_dau0_Mother_cosAngle
histogram of Track daughter 0's and mother's cos(opening-angle)
Definition: EffPlotsModule.h:194
Belle2::EffPlotsModule::m_h1_track_dau0_z0
TH1F * m_h1_track_dau0_z0
histogram of Track daughter 0's z0
Definition: EffPlotsModule.h:182
Belle2::EffPlotsModule::m_h1_MC_dau0_RMother
TH1F * m_h1_MC_dau0_RMother
histogram of MCParticle daughter 0's RMother
Definition: EffPlotsModule.h:114
Belle2::EffPlotsModule::m_h1_V0_phi_BW
TH1F * m_h1_V0_phi_BW
histogram of V0 mother's phi (backward region)
Definition: EffPlotsModule.h:237
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_pz
TH1F * m_h1_RecoTrack_dau1_pz
histogram of RecoTrack daughter 1's pz
Definition: EffPlotsModule.h:277
Belle2::EffPlotsModule::m_histoList_Efficiencies
TList * m_histoList_Efficiencies
list of histograms of efficiencies
Definition: EffPlotsModule.h:107
Belle2::EffPlotsModule::m_h3_V0
TH3F * m_h3_V0
histogram of V0 mother's pt vs theta vs phi
Definition: EffPlotsModule.h:232
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_ptMother
TH1F * m_h1_RecoTrack_dau1_ptMother
histogram of RecoTrack daughter 1's mother's pt
Definition: EffPlotsModule.h:291
Belle2::EffPlotsModule::m_h1_MC_dau1_pt
TH1F * m_h1_MC_dau1_pt
histogram of MCParticle daughter 1's pt
Definition: EffPlotsModule.h:142
Belle2::EffPlotsModule::m_h1_track_dau0_p
TH1F * m_h1_track_dau0_p
histogram of Track daughter 0's p
Definition: EffPlotsModule.h:187
Belle2::EffPlotsModule::m_h2_RecoTrack_Mother_2D
TH2F * m_h2_RecoTrack_Mother_2D
histogram of RecoTrack mother's pt vs theta
Definition: EffPlotsModule.h:308
Belle2::EffPlotsModule::nMatchedDaughters
int nMatchedDaughters(const MCParticle &the_mcParticle)
get the number of matched daughters of the MCParticle
Definition: EffPlotsModule.cc:1309
Belle2::EffPlotsModule::m_V0sType
std::string m_V0sType
type (as a string) of the selected V0
Definition: EffPlotsModule.h:94
Belle2::EffPlotsModule::m_h1_MC_Mother_p
TH1F * m_h1_MC_Mother_p
histogram of MCParticle mother's p
Definition: EffPlotsModule.h:168
Belle2::EffPlotsModule::m_h2_MC_dau1_2D
TH2F * m_h2_MC_dau1_2D
histogram of MCParticle daughter 1's pt vs theta
Definition: EffPlotsModule.h:158
Belle2::EffPlotsModule::m_h3_MC_dau1
TH3F * m_h3_MC_dau1
histogram of MCParticle daughter 1's pt vs theta vs phi
Definition: EffPlotsModule.h:141
Belle2::EffPlotsModule::m_h2_track_dau1_pVScostheta
TH2F * m_h2_track_dau1_pVScostheta
histogram of Track daughter 1's p vs cos(theta)
Definition: EffPlotsModule.h:229
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_pz
TH1F * m_h1_RecoTrack_Mother_pz
histogram of RecoTrack mother's pz
Definition: EffPlotsModule.h:300
Belle2::EffPlotsModule::m_h1_MC_Mother_phi_BW
TH1F * m_h1_MC_Mother_phi_BW
histogram of MCParticle mother's phi (backward region)
Definition: EffPlotsModule.h:170
Belle2::EffPlotsModule::m_h1_V0_p
TH1F * m_h1_V0_p
histogram of V0 mother's p
Definition: EffPlotsModule.h:235
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_thetaMother
TH1F * m_h1_RecoTrack_dau0_thetaMother
histogram of RecoTrack daughter 0's mother's theta
Definition: EffPlotsModule.h:265
Belle2::EffPlotsModule::m_V0sName
std::string m_V0sName
name of the V0s dataobjects collection
Definition: EffPlotsModule.h:89
Belle2::EffPlotsModule::m_histoList_GA
TList * m_histoList_GA
list of histograms of geometric acceptance
Definition: EffPlotsModule.h:108
Belle2::EffPlotsModule::m_h1_MC_dau0_phi_BW
TH1F * m_h1_MC_dau0_phi_BW
histogram of MCParticle daughter 0's phi (backward region)
Definition: EffPlotsModule.h:120
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_RMother
TH1F * m_h1_RecoTrack_dau0_RMother
histogram of RecoTrack daughter 0's RMother
Definition: EffPlotsModule.h:249
Belle2::EffPlotsModule::m_h2_track_dau1_2D
TH2F * m_h2_track_dau1_2D
histogram of Track daughter 1's pt vs theta
Definition: EffPlotsModule.h:226
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_phi_FW
TH1F * m_h1_RecoTrack_Mother_phi_FW
histogram of RecoTrack mother's phi (forward region)
Definition: EffPlotsModule.h:305
Belle2::EffPlotsModule::beginRun
void beginRun() override
Called when entering a new run.
Definition: EffPlotsModule.cc:507
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_costheta
TH1F * m_h1_RecoTrack_Mother_costheta
histogram of RecoTrack mother's cos(theta)
Definition: EffPlotsModule.h:307
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_z0
TH1F * m_h1_RecoTrack_dau0_z0
histogram of RecoTrack daughter 0's z0
Definition: EffPlotsModule.h:248
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_RMother
TH1F * m_h1_RecoTrack_Mother_RMother
histogram of RecoTrack mother's RMother
Definition: EffPlotsModule.h:297
Belle2::EffPlotsModule::m_h2_track_dau0_2D_BP
TH2F * m_h2_track_dau0_2D_BP
histogram of Track daughter 0's pt vs theta (beam pipe)
Definition: EffPlotsModule.h:202
Belle2::EffPlotsModule::m_h1_MC_dau0_theta
TH1F * m_h1_MC_dau0_theta
histogram of MCParticle daughter 0's theta
Definition: EffPlotsModule.h:123
Belle2::EffPlotsModule::m_h1_MC_dau1_phiMother_FW
TH1F * m_h1_MC_dau1_phiMother_FW
histogram of MCParticle daughter 1's mother's phi (forward region)
Definition: EffPlotsModule.h:155
Belle2::EffPlotsModule::m_h2_RecoTrack_dau0_2DMother
TH2F * m_h2_RecoTrack_dau0_2DMother
histogram of RecoTrack daughter 0's mother's pt vs theta
Definition: EffPlotsModule.h:269
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_phi_barrel
TH1F * m_h1_RecoTrack_dau1_phi_barrel
histogram of RecoTrack daughter 1's phi (barrel region)
Definition: EffPlotsModule.h:281
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_theta
TH1F * m_h1_RecoTrack_dau1_theta
histogram of RecoTrack daughter 1's theta
Definition: EffPlotsModule.h:283
Belle2::EffPlotsModule::m_h2_RecoTrack_dau1_2D
TH2F * m_h2_RecoTrack_dau1_2D
histogram of RecoTrack daughter 1's pt vs theta
Definition: EffPlotsModule.h:292
Belle2::EffPlotsModule::m_h3_RecoTrack_dau0
TH3F * m_h3_RecoTrack_dau0
histogram of RecoTrack daughter 0's pt vs theta vs phi
Definition: EffPlotsModule.h:250
Belle2::EffPlotsModule::m_h2_RecoTrack_dau0_2D
TH2F * m_h2_RecoTrack_dau0_2D
histogram of RecoTrack daughter 0's pt vs theta
Definition: EffPlotsModule.h:267
Belle2::EffPlotsModule::m_MCDaughter0
MCParticle * m_MCDaughter0
daughter 0 of a decayed MCParticle
Definition: EffPlotsModule.h:101
Belle2::EffPlotsModule::m_h1_MC_dau0_phiMother_barrel
TH1F * m_h1_MC_dau0_phiMother_barrel
histogram of MCParticle daughter 0's mother's phi (barrel region)
Definition: EffPlotsModule.h:128
Belle2::EffPlotsModule::m_RecoTracksName
std::string m_RecoTracksName
name of the RecoTracks dataobjects collection
Definition: EffPlotsModule.h:90
Belle2::EffPlotsModule::m_h1_MC_dau0_ptMother
TH1F * m_h1_MC_dau0_ptMother
histogram of MCParticle daughter 0's mother's pt
Definition: EffPlotsModule.h:131
Belle2::EffPlotsModule::m_h1_track_dau0_d0
TH1F * m_h1_track_dau0_d0
list of histograms filled per Tracks/V0 found in the event
Definition: EffPlotsModule.h:181
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_phiMother_barrel
TH1F * m_h1_RecoTrack_dau0_phiMother_barrel
histogram of RecoTrack daughter 0's mother's phi (barrel region)
Definition: EffPlotsModule.h:263
Belle2::EffPlotsModule::m_h1_track_dau1_costheta
TH1F * m_h1_track_dau1_costheta
histogram of Track daughter 1's cos(theta)
Definition: EffPlotsModule.h:218
Belle2::EffPlotsModule::m_h2_RecoTrack_dau1_2D_BP
TH2F * m_h2_RecoTrack_dau1_2D_BP
histogram of RecoTrack daughter 1's pt vs theta (beam pipe)
Definition: EffPlotsModule.h:293
Belle2::EffPlotsModule::m_h1_MC_dau1_phiMother_BW
TH1F * m_h1_MC_dau1_phiMother_BW
histogram of MCParticle daughter 1's mother's phi (backward region)
Definition: EffPlotsModule.h:153
Belle2::EffPlotsModule::m_histoList_MCParticles
TList * m_histoList_MCParticles
list of histograms for MCParticles
Definition: EffPlotsModule.h:104
Belle2::EffPlotsModule::m_h1_MC_dau1_d0
TH1F * m_h1_MC_dau1_d0
histogram of MCParticle daughter 1's d0
Definition: EffPlotsModule.h:138
Belle2::EffPlotsModule::m_h1_MC_dau1_phi_FW
TH1F * m_h1_MC_dau1_phi_FW
histogram of MCParticle daughter 1's phi (forward region)
Definition: EffPlotsModule.h:148
Belle2::EffPlotsModule::m_h1_MC_dau1_phiMother_total
TH1F * m_h1_MC_dau1_phiMother_total
histogram of MCParticle daughter 1's mother's phi
Definition: EffPlotsModule.h:152
Belle2::EffPlotsModule::m_h2_MC_dau0_2DMother
TH2F * m_h2_MC_dau0_2DMother
histogram of MCParticle daughter 0's mother's pt vs theta
Definition: EffPlotsModule.h:134
Belle2::EffPlotsModule::m_h1_MC_dau0_phi_FW
TH1F * m_h1_MC_dau0_phi_FW
histogram of MCParticle daughter 0's phi (forward region)
Definition: EffPlotsModule.h:122
Belle2::EffPlotsModule::EffPlotsModule
EffPlotsModule()
Constructor.
Definition: EffPlotsModule.cc:33
Belle2::EffPlotsModule::m_MCParticlesName
std::string m_MCParticlesName
user-defined parameters
Definition: EffPlotsModule.h:88
Belle2::EffPlotsModule::m_h1_track_dau0_phi_BW
TH1F * m_h1_track_dau0_phi_BW
histogram of Track daughter 0's phi (backward region)
Definition: EffPlotsModule.h:189
Belle2::EffPlotsModule::m_h1_MC_dau0_phiMother_BW
TH1F * m_h1_MC_dau0_phiMother_BW
histogram of MCParticle daughter 0's mother's phi (backward region)
Definition: EffPlotsModule.h:127
Belle2::EffPlotsModule::m_h2_track_dau0_2DMother
TH2F * m_h2_track_dau0_2DMother
histogram of Track daughter 0's mother's pt vs theta
Definition: EffPlotsModule.h:203
Belle2::EffPlotsModule::m_h1_MC_dau0_pz
TH1F * m_h1_MC_dau0_pz
histogram of MCParticle daughter 0's pz
Definition: EffPlotsModule.h:117
Belle2::EffPlotsModule::m_h1_track_dau0_phiMother_FW
TH1F * m_h1_track_dau0_phiMother_FW
histogram of Track daughter 0's mother's phi (forward region)
Definition: EffPlotsModule.h:198
Belle2::EffPlotsModule::m_h1_MC_Mother_phi_barrel
TH1F * m_h1_MC_Mother_phi_barrel
histogram of MCParticle mother's phi (barrel region)
Definition: EffPlotsModule.h:171
Belle2::EffPlotsModule::isLambda0
bool isLambda0(const MCParticle &the_mcParticle)
determine if the MCParticle is a Lambda0
Definition: EffPlotsModule.cc:1285
Belle2::EffPlotsModule::m_h1_MC_dau1_costheta
TH1F * m_h1_MC_dau1_costheta
histogram of MCParticle daughter 1's cos(theta)
Definition: EffPlotsModule.h:150
Belle2::EffPlotsModule::m_h2_RecoTrack_dau1_2DMother
TH2F * m_h2_RecoTrack_dau1_2DMother
histogram of RecoTrack daughter 1's mother's pt vs theta
Definition: EffPlotsModule.h:294
Belle2::EffPlotsModule::m_h1_track_dau1_theta
TH1F * m_h1_track_dau1_theta
histogram of Track daughter 1's theta
Definition: EffPlotsModule.h:217
Belle2::EffPlotsModule::m_h3_MC_Mother
TH3F * m_h3_MC_Mother
histogram of MCParticle mother's pt vs theta vs phi
Definition: EffPlotsModule.h:165
Belle2::EffPlotsModule::m_h1_MC_dau1_phi_barrel
TH1F * m_h1_MC_dau1_phi_barrel
histogram of MCParticle daughter 1's phi (barrel region)
Definition: EffPlotsModule.h:147
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_pt
TH1F * m_h1_RecoTrack_Mother_pt
histogram of RecoTrack mother's pt
Definition: EffPlotsModule.h:299
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_phiMother_FW
TH1F * m_h1_RecoTrack_dau1_phiMother_FW
histogram of RecoTrack daughter 1's mother's phi (forward region)
Definition: EffPlotsModule.h:289
Belle2::EffPlotsModule::m_h3_track_dau0
TH3F * m_h3_track_dau0
histogram of Track daughter 0's pt vs theta vs phi
Definition: EffPlotsModule.h:184
Belle2::EffPlotsModule::m_h1_track_dau1_phi_FW
TH1F * m_h1_track_dau1_phi_FW
histogram of Track daughter 1's phi (forward region)
Definition: EffPlotsModule.h:216
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_phiMother_BW
TH1F * m_h1_RecoTrack_dau1_phiMother_BW
histogram of RecoTrack daughter 1's mother's phi (backward region)
Definition: EffPlotsModule.h:287
Belle2::EffPlotsModule::m_h2_MC_dau0_pVScostheta
TH2F * m_h2_MC_dau0_pVScostheta
histogram of MCParticle daughter 0's p vs cos(theta)
Definition: EffPlotsModule.h:135
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_phi
TH1F * m_h1_RecoTrack_dau0_phi
histogram of RecoTrack daughter 0's phi
Definition: EffPlotsModule.h:254
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_costheta
TH1F * m_h1_RecoTrack_dau0_costheta
histogram of RecoTrack daughter 0's cos(theta)
Definition: EffPlotsModule.h:259
Belle2::EffPlotsModule::m_h2_MC_dau0_2D
TH2F * m_h2_MC_dau0_2D
histogram of MCParticle daughter 0's pt vs theta
Definition: EffPlotsModule.h:132
Belle2::EffPlotsModule::m_h1_MC_dau0_p
TH1F * m_h1_MC_dau0_p
histogram of MCParticle daughter 0's p
Definition: EffPlotsModule.h:118
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_d0
TH1F * m_h1_RecoTrack_dau1_d0
histogram of RecoTrack daughter 1's d0
Definition: EffPlotsModule.h:272
Belle2::EffPlotsModule::m_h2_MC_dau1_2DMother
TH2F * m_h2_MC_dau1_2DMother
histogram of MCParticle daughter 1's mother's pt vs theta
Definition: EffPlotsModule.h:160
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_thetaMother
TH1F * m_h1_RecoTrack_dau1_thetaMother
histogram of RecoTrack daughter 1's mother's theta
Definition: EffPlotsModule.h:290
Belle2::EffPlotsModule::m_h1_MC_dau0_phiMother_FW
TH1F * m_h1_MC_dau0_phiMother_FW
histogram of MCParticle daughter 0's mother's phi (forward region)
Definition: EffPlotsModule.h:129
Belle2::EffPlotsModule::m_h2_V0_Mother_pVScostheta
TH2F * m_h2_V0_Mother_pVScostheta
histogram of V0 mother's p vs cos(theta)
Definition: EffPlotsModule.h:244
Belle2::EffPlotsModule::m_h1_MC_Mother_RMother
TH1F * m_h1_MC_Mother_RMother
histogram of MCParticle mother's RMother
Definition: EffPlotsModule.h:164
Belle2::EffPlotsModule::m_h1_track_dau1_phiMother_FW
TH1F * m_h1_track_dau1_phiMother_FW
histogram of Track daughter 1's mother's phi (forward region)
Definition: EffPlotsModule.h:223
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_phi_BW
TH1F * m_h1_RecoTrack_Mother_phi_BW
histogram of RecoTrack mother's phi (backward region)
Definition: EffPlotsModule.h:303
Belle2::EffPlotsModule::m_h1_MC_Mother_pz
TH1F * m_h1_MC_Mother_pz
histogram of MCParticle mother's pz
Definition: EffPlotsModule.h:167
Belle2::EffPlotsModule::m_h1_track_dau0_thetaMother
TH1F * m_h1_track_dau0_thetaMother
histogram of Track daughter 0's mother's theta
Definition: EffPlotsModule.h:199
Belle2::EffPlotsModule::m_h1_MC_dau1_phiMother_barrel
TH1F * m_h1_MC_dau1_phiMother_barrel
histogram of MCParticle daughter 1's mother's phi (barrel region)
Definition: EffPlotsModule.h:154
Belle2::EffPlotsModule::m_h1_MC_dau1_z0
TH1F * m_h1_MC_dau1_z0
histogram of MCParticle daughter 1's z0
Definition: EffPlotsModule.h:139
Belle2::EffPlotsModule::m_h1_V0_theta
TH1F * m_h1_V0_theta
histogram of V0 mother's theta
Definition: EffPlotsModule.h:240
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_RMother
TH1F * m_h1_RecoTrack_dau1_RMother
histogram of RecoTrack daughter 1's RMother
Definition: EffPlotsModule.h:274
Belle2::EffPlotsModule::m_h1_track_dau1_RMother
TH1F * m_h1_track_dau1_RMother
histogram of Track daughter 1's RMother
Definition: EffPlotsModule.h:208
Belle2::EffPlotsModule::m_histoList_RecoTracks
TList * m_histoList_RecoTracks
list of histograms for RecoTracks
Definition: EffPlotsModule.h:105
Belle2::EffPlotsModule::m_h1_MC_dau0_phiMother_total
TH1F * m_h1_MC_dau0_phiMother_total
histogram of MCParticle daughter 0's mother's phi
Definition: EffPlotsModule.h:126
Belle2::EffPlotsModule::m_h1_MC_Mother_costheta
TH1F * m_h1_MC_Mother_costheta
histogram of MCParticle mother's cos(theta)
Definition: EffPlotsModule.h:174
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_phi
TH1F * m_h1_RecoTrack_Mother_phi
histogram of RecoTrack mother's phi
Definition: EffPlotsModule.h:302
Belle2::EffPlotsModule::m_h1_track_dau0_phi
TH1F * m_h1_track_dau0_phi
histogram of Track daughter 0's phi
Definition: EffPlotsModule.h:188
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_pt
TH1F * m_h1_RecoTrack_dau1_pt
histogram of RecoTrack daughter 1's pt
Definition: EffPlotsModule.h:276
Belle2::EffPlotsModule::m_h1_track_dau1_phiMother_total
TH1F * m_h1_track_dau1_phiMother_total
histogram of Track daughter 1's mother's phi
Definition: EffPlotsModule.h:220
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_p
TH1F * m_h1_RecoTrack_dau0_p
histogram of RecoTrack daughter 0's p
Definition: EffPlotsModule.h:253
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_Mother_cosAngle
TH1F * m_h1_RecoTrack_dau0_Mother_cosAngle
histogram of RecoTrack daughter 0's and mother's cos(opening-angle)
Definition: EffPlotsModule.h:260
Belle2::EffPlotsModule::m_h1_track_dau0_phiMother_BW
TH1F * m_h1_track_dau0_phiMother_BW
histogram of Track daughter 0's mother's phi (backward region)
Definition: EffPlotsModule.h:196
Belle2::EffPlotsModule::m_h1_MC_dau0_d0
TH1F * m_h1_MC_dau0_d0
list of histograms filled per MCParticle found in the event
Definition: EffPlotsModule.h:112
Belle2::EffPlotsModule::m_MCRecoTracksName
std::string m_MCRecoTracksName
name of the MCRecoTracks dataobjects collection
Definition: EffPlotsModule.h:91
Belle2::EffPlotsModule::m_MCParticles
StoreArray< MCParticle > m_MCParticles
MCParticle StoreArray.
Definition: EffPlotsModule.h:96
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_pt
TH1F * m_h1_RecoTrack_dau0_pt
histogram of RecoTrack daughter 0's pt
Definition: EffPlotsModule.h:251
Belle2::EffPlotsModule::~EffPlotsModule
~EffPlotsModule()
Destructor.
Definition: EffPlotsModule.cc:55
Belle2::EffPlotsModule::m_h2_MC_dau0_2D_BP
TH2F * m_h2_MC_dau0_2D_BP
histogram of MCParticle daughter 0's pt vs theta (beam pipe)
Definition: EffPlotsModule.h:133
Belle2::EffPlotsModule::m_h1_track_dau1_phiMother_barrel
TH1F * m_h1_track_dau1_phiMother_barrel
histogram of Track daughter 1's mother's phi (barrel region)
Definition: EffPlotsModule.h:222
Belle2::EffPlotsModule::m_h1_RecoTrack_dau0_phi_barrel
TH1F * m_h1_RecoTrack_dau0_phi_barrel
histogram of RecoTrack daughter 0's phi (barrel region)
Definition: EffPlotsModule.h:256
Belle2::EffPlotsModule::m_h1_MC_Mother_theta
TH1F * m_h1_MC_Mother_theta
histogram of MCParticle mother's theta
Definition: EffPlotsModule.h:173
Belle2::EffPlotsModule::m_h1_track_dau1_phiMother_BW
TH1F * m_h1_track_dau1_phiMother_BW
histogram of Track daughter 1's mother's phi (backward region)
Definition: EffPlotsModule.h:221
Belle2::EffPlotsModule::m_h2_RecoTrack_dau0_pVScostheta
TH2F * m_h2_RecoTrack_dau0_pVScostheta
histogram of RecoTrack daughter 0's p vs cos(theta)
Definition: EffPlotsModule.h:270
Belle2::EffPlotsModule::m_MCDaughter1
MCParticle * m_MCDaughter1
daughter 1 of a decayed MCParticle
Definition: EffPlotsModule.h:102
Belle2::EffPlotsModule::m_histoList_Tracks
TList * m_histoList_Tracks
list of histograms for Tracks
Definition: EffPlotsModule.h:106
Belle2::EffPlotsModule::m_h1_RecoTrack_dau1_phiMother_barrel
TH1F * m_h1_RecoTrack_dau1_phiMother_barrel
histogram of RecoTrack daughter 1's mother's phi (barrel region)
Definition: EffPlotsModule.h:288
Belle2::EffPlotsModule::m_h2_track_dau1_2D_BP
TH2F * m_h2_track_dau1_2D_BP
histogram of Track daughter 1's pt vs theta (beam pipe)
Definition: EffPlotsModule.h:227
Belle2::EffPlotsModule::m_h1_MC_dau1_phi
TH1F * m_h1_MC_dau1_phi
histogram of MCParticle daughter 1's phi
Definition: EffPlotsModule.h:145
Belle2::EffPlotsModule::m_h2_MC_Mother_2D
TH2F * m_h2_MC_Mother_2D
histogram of MCParticle mother's pt vs theta
Definition: EffPlotsModule.h:175
Belle2::EffPlotsModule::m_h1_track_dau1_d0
TH1F * m_h1_track_dau1_d0
histogram of Track daughter 1's d0
Definition: EffPlotsModule.h:206
Belle2::EffPlotsModule::m_h1_track_dau0_phi_FW
TH1F * m_h1_track_dau0_phi_FW
histogram of Track daughter 0's phi (forward region)
Definition: EffPlotsModule.h:191
Belle2::EffPlotsModule::m_h1_MC_dau0_phi_barrel
TH1F * m_h1_MC_dau0_phi_barrel
histogram of MCParticle daughter 0's phi (barrel region)
Definition: EffPlotsModule.h:121
Belle2::EffPlotsModule::m_h1_MC_dau1_pz
TH1F * m_h1_MC_dau1_pz
histogram of MCParticle daughter 1's pz
Definition: EffPlotsModule.h:143
Belle2::EffPlotsModule::m_h1_track_dau1_phi_BW
TH1F * m_h1_track_dau1_phi_BW
histogram of Track daughter 1's phi (backward region)
Definition: EffPlotsModule.h:214
Belle2::EffPlotsModule::m_h1_RecoTrack_Mother_p
TH1F * m_h1_RecoTrack_Mother_p
histogram of RecoTrack mother's p
Definition: EffPlotsModule.h:301
Belle2::EffPlotsModule::m_h1_track_dau1_thetaMother
TH1F * m_h1_track_dau1_thetaMother
histogram of Track daughter 1's mother's theta
Definition: EffPlotsModule.h:224
Belle2::MCParticleInfo
This struct is used by the TrackingPerformanceEvaluation Module to save information of reconstructed ...
Definition: MCParticleInfo.h:22
Belle2::MCParticleInfo::getPt
double getPt()
Getter for transverse momentum.
Definition: MCParticleInfo.h:36
Belle2::MCParticleInfo::getPtheta
double getPtheta()
Getter for theta of momentum vector.
Definition: MCParticleInfo.h:50
Belle2::MCParticleInfo::getZ0
double getZ0()
Getter for Z0.
Definition: MCParticleInfo.cc:73
Belle2::MCParticleInfo::getPz
double getPz()
Getter for z component of momentum.
Definition: MCParticleInfo.h:34
Belle2::MCParticleInfo::getPphi
double getPphi()
Getter for phi of momentum vector.
Definition: MCParticleInfo.h:52
Belle2::MCParticleInfo::getD0
double getD0()
Getter for D0.
Definition: MCParticleInfo.cc:24
Belle2::MCParticleInfo::getP
double getP()
Getter for magnitut of momentum.
Definition: MCParticleInfo.h:38
Belle2::MCParticle
A Class to store the Monte Carlo particle information.
Definition: MCParticle.h:32
Belle2::MCParticle::getDaughters
std::vector< Belle2::MCParticle * > getDaughters() const
Get vector of all daughter particles, empty vector if none.
Definition: MCParticle.cc:52
Belle2::MCParticle::getPDG
int getPDG() const
Return PDG code of particle.
Definition: MCParticle.h:112
Belle2::MCParticle::getMomentum
ROOT::Math::XYZVector getMomentum() const
Return momentum.
Definition: MCParticle.h:198
Belle2::Module
Base class for Modules.
Definition: Module.h:72
Belle2::Module::setDescription
void setDescription(const std::string &description)
Sets the description of the module.
Definition: Module.cc:214
Belle2::PerformanceEvaluationBaseClass::createHistogram1D
TH1F * createHistogram1D(const char *name, const char *title, Int_t nbins, Double_t min, Double_t max, const char *xtitle, TList *histoList=nullptr)
Create a 1D histogram and add it to the TList of 1D-histograms.
Definition: PerformanceEvaluationBaseClass.cc:26
Belle2::PerformanceEvaluationBaseClass::effPlot2D
TH2F * effPlot2D(TH2F *h2_den, TH2F *h2_num, const char *name, const char *title, bool geo_accettance, TList *histoList=nullptr)
Create a 2D efficiency histogram and add it to the TList of 2D-histograms.
Definition: PerformanceEvaluationBaseClass.cc:440
Belle2::PerformanceEvaluationBaseClass::effPlot1D
TH1F * effPlot1D(TH1F *h1_den, TH1F *h1_num, const char *name, const char *title, bool geo_accettance, TList *histoList=nullptr)
Create a 1D efficiency histogram and add it to the TList of 1D-histograms.
Definition: PerformanceEvaluationBaseClass.cc:332
Belle2::PerformanceEvaluationBaseClass::duplicateHistogram
TH1 * duplicateHistogram(const char *newname, const char *newtitle, TH1 *h, TList *histoList=nullptr)
Make a copy of a 1D histogram and add it to the TList of 1D-histograms.
Definition: PerformanceEvaluationBaseClass.cc:137
Belle2::PerformanceEvaluationBaseClass::V0FinderEff
TH1F * V0FinderEff(TH1F *h1_dau0, TH1F *h1_dau1, TH1F *h1_Mother, const char *name, const char *title, TList *histoList=nullptr)
Create a 1D efficiency histogram for V0 finding and add it to the TList of 1D-histograms.
Definition: PerformanceEvaluationBaseClass.cc:644
Belle2::PerformanceEvaluationBaseClass::m_histoList
TList * m_histoList
List of performance-evaluation histograms.
Definition: PerformanceEvaluationBaseClass.h:38
Belle2::PerformanceEvaluationBaseClass::geoAcc1D
TH1F * geoAcc1D(TH1F *h1_den, TH1F *h1_num, const char *name, const char *title, TList *histoList=nullptr)
Create a 1D efficiency histogram for geometric acceptance and add it to the TList of 1D-histograms.
Definition: PerformanceEvaluationBaseClass.cc:580
Belle2::PerformanceEvaluationBaseClass::createHistogram3D
TH3F * createHistogram3D(const char *name, const char *title, Int_t nbinsX, Double_t minX, Double_t maxX, const char *titleX, Int_t nbinsY, Double_t minY, Double_t maxY, const char *titleY, Int_t nbinsZ, Double_t minZ, Double_t maxZ, const char *titleZ, TList *histoList=nullptr)
Create a 3D histogram and add it to the TList of 3D-histograms.
Definition: PerformanceEvaluationBaseClass.cc:93
Belle2::PerformanceEvaluationBaseClass::createHistogram2D
TH2F * createHistogram2D(const char *name, const char *title, Int_t nbinsX, Double_t minX, Double_t maxX, const char *titleX, Int_t nbinsY, Double_t minY, Double_t maxY, const char *titleY, TList *histoList=nullptr)
Create a 2D histogram and add it to the TList of 2D-histograms.
Definition: PerformanceEvaluationBaseClass.cc:56
Belle2::PerformanceEvaluationBaseClass::geoAcc2D
TH2F * geoAcc2D(TH2F *h2_den, TH2F *h2_num, const char *name, const char *title, TList *histoList=nullptr)
Create a 2D efficiency histogram for geometric acceptance and add it to the TList of 2D-histograms.
Definition: PerformanceEvaluationBaseClass.cc:605
Belle2::PerformanceEvaluationBaseClass::m_rootFilePtr
TFile * m_rootFilePtr
pointer at root file used for storing histograms
Definition: PerformanceEvaluationBaseClass.h:140
Belle2::RelationVector
Class for type safe access to objects that are referred to in relations.
Definition: RelationVector.h:67
Belle2::RelationVector::size
size_t size() const
Get number of relations.
Definition: RelationVector.h:88
Belle2::RelationsInterface::getRelated
T * getRelated(const std::string &name="", const std::string &namedRelation="") const
Get the object to or from which this object has a relation.
Definition: RelationsObject.h:278
Belle2::StoreAccessorBase::isRequired
bool isRequired(const std::string &name="")
Ensure this array/object has been registered previously.
Definition: StoreAccessorBase.h:78
Belle2::StoreArray
Accessor to arrays stored in the data store.
Definition: StoreArray.h:113
Belle2::Track
Class that bundles various TrackFitResults.
Definition: Track.h:25
Belle2::Unit::T
static const double T
[tesla]
Definition: Unit.h:120
Belle2::REG_MODULE
REG_MODULE(arichBtest)
Register the Module.
Belle2::Module::addParam
void addParam(const std::string &name, T &paramVariable, const std::string &description, const T &defaultValue)
Adds a new parameter to the module.
Definition: Module.h:560
Belle2::BFieldManager::getField
static void getField(const double *pos, double *field)
return the magnetic field at a given position.
Definition: BFieldManager.h:91
Belle2::sqrt
double sqrt(double a)
sqrt for double
Definition: beamHelpers.h:28
Belle2
Abstract base class for different kinds of events.
Definition: MillepedeAlgorithm.h:17