Belle II Software development
OverlapResidualsModule.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/VXDOverlaps/OverlapResidualsModule.h>
10#include <framework/datastore/StoreArray.h>
11#include <framework/gearbox/Unit.h>
12#include <framework/core/Environment.h>
13#include <svd/dataobjects/SVDCluster.h>
14#include <svd/dataobjects/SVDRecoDigit.h>
15#include <svd/dataobjects/SVDShaperDigit.h>
16#include <svd/dataobjects/SVDTrueHit.h>
17#include <vxd/dataobjects/VxdID.h>
18#include <vxd/geometry/GeoCache.h>
19#include <vxd/geometry/SensorInfoBase.h>
20#include <tracking/dataobjects/RecoTrack.h>
21#include <tracking/dataobjects/RecoHitInformation.h>
22#include <genfit/TrackPoint.h>
23#include <Math/Vector3D.h>
24#include <TDirectory.h>
25#include <Math/Boost.h>
26#include <math.h>
27#include <iostream>
28#include <algorithm>
29#include <mdst/dataobjects/Track.h>
30#include <mdst/dataobjects/TrackFitResult.h>
31
32using namespace Belle2;
33
34
35//-----------------------------------------------------------------
36// Register the Module
37//-----------------------------------------------------------------
38REG_MODULE(OverlapResiduals);
39
40//-----------------------------------------------------------------
41// Implementation
42//-----------------------------------------------------------------
43
44OverlapResidualsModule::OverlapResidualsModule() : HistoModule()
45{
46 //Set module properties
47 setDescription("The module studies consecutive hits in overlapping sensors of a same VXD layer, and the differences of their residuals, to monitor the detector alignment.");
48 //Parameter to take only specific RecoTracks as input
49 addParam("recoTracksStoreArrayName", m_recoTracksStoreArrayName, "StoreArray name of the input and output RecoTracks.",
50 m_recoTracksStoreArrayName);
51 //Parameter to produce TTrees storing global information on VXD overlaps
52 addParam("ExpertLevel", m_ExpertLevel,
53 "If set, enables the production of TTrees containing low level information on PXD and SVD overlaps", false);
54}
55
56void OverlapResidualsModule::initialize()
57{
58 StoreArray<RecoTrack> recoTracks(m_recoTracksStoreArrayName);
59 recoTracks.isOptional();
60 //Register histograms (calls back defineHisto)
61 REG_HISTOGRAM
62}
63
64void OverlapResidualsModule::defineHisto()
65{
66 //Create directory to store monitoring histograms
67 TDirectory* oldDir = gDirectory;
68 TDirectory* ResDir = oldDir->mkdir("Monitoring_VXDOverlaps");
69 ResDir->cd();
70 //Define histograms of residuals differences
71 h_U_DeltaRes = new TH1F("h_U_DeltaRes", "Histrogram of residual difference #Delta res_{u} for overlapping hits", 100, -1000, 1000);
72 h_U_DeltaRes->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
73 h_U_DeltaRes->GetYaxis()->SetTitle("counts");
74 h_V_DeltaRes = new TH1F("h_V_DeltaRes", "Histrogram of residual difference #Delta res_{v} for overlapping hits", 100, -1000, 1000);
75 h_V_DeltaRes->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
76 h_V_DeltaRes->GetYaxis()->SetTitle("counts");
77 h_U_DeltaRes_PXD = new TH1F("h_U_DeltaRes_PXD", "Histrogram of residual difference #Delta res_{u} for overlapping PXD hits", 100,
78 -1000, 1000);
79 h_U_DeltaRes_PXD->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
80 h_U_DeltaRes_PXD->GetYaxis()->SetTitle("counts");
81 h_U_DeltaRes_PXD_Lyr1 = new TH1F("h_U_DeltaRes_PXD_Lyr1",
82 "Layer1: histrogram of residual difference #Delta res_{u} for overlapping PXD hits", 100, -1000, 1000);
83 h_U_DeltaRes_PXD_Lyr1->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
84 h_U_DeltaRes_PXD_Lyr1->GetYaxis()->SetTitle("counts");
85 h_U_DeltaRes_PXD_Lyr2 = new TH1F("h_U_DeltaRes_PXD_Lyr2",
86 "Layer 2: hstrogram of residual difference #Delta res_{u} for overlapping PXD hits", 100, -1000, 1000);
87 h_U_DeltaRes_PXD_Lyr2->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
88 h_U_DeltaRes_PXD_Lyr2->GetYaxis()->SetTitle("counts");
89 h_V_DeltaRes_PXD = new TH1F("h_V_DeltaRes_PXD", "Histrogram of residual difference #Delta res_{v} for overlapping PXD hits", 100,
90 -1000, 1000);
91 h_V_DeltaRes_PXD->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
92 h_V_DeltaRes_PXD->GetYaxis()->SetTitle("counts");
93 h_V_DeltaRes_PXD_Lyr1 = new TH1F("h_V_DeltaRes_PXD_Lyr1",
94 "Layer 1: histrogram of residual difference #Delta res_{v} for overlapping PXD hits", 100, -1000, 1000);
95 h_V_DeltaRes_PXD_Lyr1->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
96 h_V_DeltaRes_PXD_Lyr1->GetYaxis()->SetTitle("counts");
97 h_V_DeltaRes_PXD_Lyr2 = new TH1F("h_V_DeltaRes_PXD_Lyr2",
98 "Layer 2: histrogram of residual difference #Delta res_{v} for overlapping PXD hits", 100, -1000, 1000);
99 h_V_DeltaRes_PXD_Lyr2->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
100 h_V_DeltaRes_PXD_Lyr2->GetYaxis()->SetTitle("counts");
101 h_U_DeltaRes_SVD = new TH1F("h_U_DeltaRes_SVD", "Histrogram of residual difference #Delta res_{u} for overlapping SVD hits", 100,
102 -1000, 1000);
103 h_U_DeltaRes_SVD->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
104 h_U_DeltaRes_SVD->GetYaxis()->SetTitle("counts");
105 h_U_DeltaRes_SVD_Lyr3 = new TH1F("h_U_DeltaRes_SVD_Lyr3",
106 "Layer 3: histrogram of residual difference #Delta res_{u} for overlapping SVD hits", 100, -1000, 1000);
107 h_U_DeltaRes_SVD_Lyr3->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
108 h_U_DeltaRes_SVD_Lyr3->GetYaxis()->SetTitle("counts");
109 h_U_DeltaRes_SVD_Lyr4 = new TH1F("h_U_DeltaRes_SVD_Lyr4",
110 "Layer 4: histrogram of residual difference #Delta res_{u} for overlapping SVD hits", 100, -1000, 1000);
111 h_U_DeltaRes_SVD_Lyr4->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
112 h_U_DeltaRes_SVD_Lyr4->GetYaxis()->SetTitle("counts");
113 h_U_DeltaRes_SVD_Lyr5 = new TH1F("h_U_DeltaRes_SVD_Lyr5",
114 "Layer 5: histrogram of residual difference #Delta res_{u} for overlapping SVD hits", 100, -1000, 1000);
115 h_U_DeltaRes_SVD_Lyr5->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
116 h_U_DeltaRes_SVD_Lyr5->GetYaxis()->SetTitle("counts");
117 h_U_DeltaRes_SVD_Lyr6 = new TH1F("h_U_DeltaRes_SVD_Lyr6",
118 "Layer 6: histrogram of residual difference #Delta res_{u} for overlapping SVD hits", 100, -1000, 1000);
119 h_U_DeltaRes_SVD_Lyr6->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
120 h_U_DeltaRes_SVD_Lyr6->GetYaxis()->SetTitle("counts");
121 h_V_DeltaRes_SVD = new TH1F("h_V_DeltaRes_SVD", "Histrogram of residual difference #Delta res_{v} for overlapping SVD hits", 100,
122 -1000, 1000);
123 h_V_DeltaRes_SVD->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
124 h_V_DeltaRes_SVD->GetYaxis()->SetTitle("counts");
125 h_V_DeltaRes_SVD_Lyr3 = new TH1F("h_V_DeltaRes_SVD_Lyr3",
126 "Layer 3: histrogram of residual difference #Delta res_{v} for overlapping SVD hits", 100, -1000, 1000);
127 h_V_DeltaRes_SVD_Lyr3->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
128 h_V_DeltaRes_SVD_Lyr3->GetYaxis()->SetTitle("counts");
129 h_V_DeltaRes_SVD_Lyr4 = new TH1F("h_V_DeltaRes_SVD_Lyr4",
130 "Layer 4: histrogram of residual difference #Delta res_{v} for overlapping SVD hits", 100, -1000, 1000);
131 h_V_DeltaRes_SVD_Lyr4->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
132 h_V_DeltaRes_SVD_Lyr4->GetYaxis()->SetTitle("counts");
133 h_V_DeltaRes_SVD_Lyr5 = new TH1F("h_V_DeltaRes_SVD_Lyr5",
134 "Layer 5: histrogram of residual difference #Delta res_{v} for overlapping SVD hits", 100, -1000, 1000);
135 h_V_DeltaRes_SVD_Lyr5->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
136 h_V_DeltaRes_SVD_Lyr5->GetYaxis()->SetTitle("counts");
137 h_V_DeltaRes_SVD_Lyr6 = new TH1F("h_V_DeltaRes_SVD_Lyr6",
138 "Layer 6: histrogram of residual difference #Delta res_{v} for overlapping SVD hits", 100, -1000, 1000);
139 h_V_DeltaRes_SVD_Lyr6->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
140 h_V_DeltaRes_SVD_Lyr6->GetYaxis()->SetTitle("counts");
141 h_SVDstrips_Mult = new TH1F("h_SVDstrips_Mult", "SVD strips multipicity for SVD clusters in overlapping sensors", 15, 0.5, 15.5);
142 h_SVDstrips_Mult->GetXaxis()->SetTitle("N. of SVD strips contributing to the cluster");
143 h_SVDstrips_Mult->GetYaxis()->SetTitle("counts");
144 //Define 2D histograms: difference of u-residuals vs phi of VXD overlaps for each layer (1 to 6)
145 h_DeltaResUPhi_Lyr1 = new TH2F("h_DeltaResUPhi_Lyr1", "Layer 1: #Delta res_{u} vs #phi", 200, -3.4, 3.4, 100, -200, 200);
146 h_DeltaResUPhi_Lyr1->GetXaxis()->SetTitle("#phi(rad)");
147 h_DeltaResUPhi_Lyr1->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
148 h_DeltaResUPhi_Lyr2 = new TH2F("h_DeltaResUPhi_Lyr2", "Layer 2: #Delta res_{u} vs #phi", 200, -3.4, 3.4, 100, -200, 200);
149 h_DeltaResUPhi_Lyr2->GetXaxis()->SetTitle("#phi(rad)");
150 h_DeltaResUPhi_Lyr2->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
151 h_DeltaResUPhi_Lyr3 = new TH2F("h_DeltaResUPhi_Lyr3", "Layer 3: #Delta res_{u} vs #phi", 200, -3.4, 3.4, 100, -1000, 1000);
152 h_DeltaResUPhi_Lyr3->GetXaxis()->SetTitle("#phi(rad)");
153 h_DeltaResUPhi_Lyr3->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
154 h_DeltaResUPhi_Lyr4 = new TH2F("h_DeltaResUPhi_Lyr4", "Layer 4: #Delta res_{u} vs #phi", 200, -3.4, 3.4, 100, -1000, 1000);
155 h_DeltaResUPhi_Lyr4->GetXaxis()->SetTitle("#phi(rad)");
156 h_DeltaResUPhi_Lyr4->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
157 h_DeltaResUPhi_Lyr5 = new TH2F("h_DeltaResUPhi_Lyr5", "Layer 5: #Delta res_{u} vs #phi", 200, -3.4, 3.4, 100, -1000, 1000);
158 h_DeltaResUPhi_Lyr5->GetXaxis()->SetTitle("#phi(rad)");
159 h_DeltaResUPhi_Lyr5->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
160 h_DeltaResUPhi_Lyr6 = new TH2F("h_DeltaResUPhi_Lyr6", "Layer 6: #Delta res_{u} vs #phi", 200, -3.4, 3.4, 100, -1000, 1000);
161 h_DeltaResUPhi_Lyr6->GetXaxis()->SetTitle("#phi(rad)");
162 h_DeltaResUPhi_Lyr6->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
163 //Define 2D histograms: difference of u-residuals vs z of VXD overlaps for each layer (1 to 6)
164 h_DeltaResUz_Lyr1 = new TH2F("h_DeltaResUz_Lyr1", "Layer 1: #Delta res_{u} vs z", 100, -4, 8, 100, -200, 200);
165 h_DeltaResUz_Lyr1->GetXaxis()->SetTitle("z (cm)");
166 h_DeltaResUz_Lyr1->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
167 h_DeltaResUz_Lyr2 = new TH2F("h_DeltaResUz_Lyr2", "Layer 2: #Delta res_{u} vs z", 100, -10, 15, 100, -200, 200);
168 h_DeltaResUz_Lyr2->GetXaxis()->SetTitle("z (cm)");
169 h_DeltaResUz_Lyr2->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
170 h_DeltaResUz_Lyr3 = new TH2F("h_DeltaResUz_Lyr3", "Layer 3: #Delta res_{u} vs z", 250, -15, 20, 100, -1000, 1000);
171 h_DeltaResUz_Lyr3->GetXaxis()->SetTitle("z (cm)");
172 h_DeltaResUz_Lyr3->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
173 h_DeltaResUz_Lyr4 = new TH2F("h_DeltaResUz_Lyr4", "Layer 4: #Delta res_{u} vs z", 250, -20, 25, 100, -1000, 1000);
174 h_DeltaResUz_Lyr4->GetXaxis()->SetTitle("z (cm)");
175 h_DeltaResUz_Lyr4->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
176 h_DeltaResUz_Lyr5 = new TH2F("h_DeltaResUz_Lyr5", "Layer 5: #Delta res_{u} vs z", 250, -25, 35, 100, -1000, 1000);
177 h_DeltaResUz_Lyr5->GetXaxis()->SetTitle("z (cm)");
178 h_DeltaResUz_Lyr5->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
179 h_DeltaResUz_Lyr6 = new TH2F("h_DeltaResUz_Lyr6", "Layer 6: #Delta res_{u} vs z", 250, -30, 45, 100, -1000, 1000);
180 h_DeltaResUz_Lyr6->GetXaxis()->SetTitle("z (cm)");
181 h_DeltaResUz_Lyr6->GetYaxis()->SetTitle("#Delta res_{u} (#mum)");
182 //Define 2D histograms: difference of u-residuals vs phi of VXD overlaps for each layer (1 to 6)
183 h_DeltaResVPhi_Lyr1 = new TH2F("h_DeltaResVPhi_Lyr1", "Layer 1: #Delta res_{v} vs #phi", 200, -3.4, 3.4, 100, -200, 200);
184 h_DeltaResVPhi_Lyr1->GetXaxis()->SetTitle("#phi(rad)");
185 h_DeltaResVPhi_Lyr1->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
186 h_DeltaResVPhi_Lyr2 = new TH2F("h_DeltaResVPhi_Lyr2", "Layer 2: #Delta res_{v} vs #phi", 200, -3.4, 3.4, 100, -200, 200);
187 h_DeltaResVPhi_Lyr2->GetXaxis()->SetTitle("#phi(rad)");
188 h_DeltaResVPhi_Lyr2->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
189 h_DeltaResVPhi_Lyr3 = new TH2F("h_DeltaResVPhi_Lyr3", "Layer 3: #Delta res_{v} vs #phi", 200, -3.4, 3.4, 100, -1000, 1000);
190 h_DeltaResVPhi_Lyr3->GetXaxis()->SetTitle("#phi(rad)");
191 h_DeltaResVPhi_Lyr3->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
192 h_DeltaResVPhi_Lyr4 = new TH2F("h_DeltaResVPhi_Lyr4", "Layer 4: #Delta res_{v} vs #phi", 200, -3.4, 3.4, 100, -1000, 1000);
193 h_DeltaResVPhi_Lyr4->GetXaxis()->SetTitle("#phi(rad)");
194 h_DeltaResVPhi_Lyr4->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
195 h_DeltaResVPhi_Lyr5 = new TH2F("h_DeltaResVPhi_Lyr5", "Layer 5: #Delta res_{v} vs #phi", 200, -3.4, 3.4, 100, -1000, 1000);
196 h_DeltaResVPhi_Lyr5->GetXaxis()->SetTitle("#phi(rad)");
197 h_DeltaResVPhi_Lyr5->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
198 h_DeltaResVPhi_Lyr6 = new TH2F("h_DeltaResVPhi_Lyr6", "Layer 6: #Delta res_{v} vs #phi", 200, -3.4, 3.4, 100, -1000, 1000);
199 h_DeltaResVPhi_Lyr6->GetXaxis()->SetTitle("#phi(rad)");
200 h_DeltaResVPhi_Lyr6->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
201 //Define 2D histograms: difference of v-residuals vs z of VXD overlaps for each layer (1 to 6)
202 h_DeltaResVz_Lyr1 = new TH2F("h_DeltaResVz_Lyr1", "Layer 1: #Delta res_{v} vs z", 100, -4, 8, 100, -200, 200);
203 h_DeltaResVz_Lyr1->GetXaxis()->SetTitle("z (cm)");
204 h_DeltaResVz_Lyr1->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
205 h_DeltaResVz_Lyr2 = new TH2F("h_DeltaResVz_Lyr2", "Layer 2: #Delta res_{v} vs z", 100, -10, 15, 100, -200, 200);
206 h_DeltaResVz_Lyr2->GetXaxis()->SetTitle("z (cm)");
207 h_DeltaResVz_Lyr2->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
208 h_DeltaResVz_Lyr3 = new TH2F("h_DeltaResVz_Lyr3", "Layer 3: #Delta res_{v} vs z", 250, -15, 20, 100, -1000, 1000);
209 h_DeltaResVz_Lyr3->GetXaxis()->SetTitle("z (cm)");
210 h_DeltaResVz_Lyr3->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
211 h_DeltaResVz_Lyr4 = new TH2F("h_DeltaResVz_Lyr4", "Layer 4: #Delta res_{v} vs z", 250, -20, 25, 100, -1000, 1000);
212 h_DeltaResVz_Lyr4->GetXaxis()->SetTitle("z (cm)");
213 h_DeltaResVz_Lyr4->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
214 h_DeltaResVz_Lyr5 = new TH2F("h_DeltaResVz_Lyr5", "Layer 5: #Delta res_{v} vs z", 250, -25, 35, 100, -1000, 1000);
215 h_DeltaResVz_Lyr5->GetXaxis()->SetTitle("z (cm)");
216 h_DeltaResVz_Lyr5->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
217 h_DeltaResVz_Lyr6 = new TH2F("h_DeltaResVz_Lyr6", "Layer 6: #Delta res_{v} vs z", 250, -30, 45, 100, -1000, 1000);
218 h_DeltaResVz_Lyr6->GetXaxis()->SetTitle("z (cm)");
219 h_DeltaResVz_Lyr6->GetYaxis()->SetTitle("#Delta res_{v} (#mum)");
220 //Restricting to SVD clusters sizes
221 for (int i = 1; i < 5; i++) {
222 //The name is the product of cluster sizes for 2 consecutive hits (maximum size considered is 2)
223 TString h_name_U = "h_U_Cl1Cl2_" + std::to_string(i);
224 TString h_name_V = "h_V_Cl1Cl2_" + std::to_string(i);
225 TString title_U = "#Delta res_{u}: SVDClusterSize_1 x SVDClusterSize_2 = " + std::to_string(i);
226 TString title_V = "#Delta res_{v}: SVDClusterSize_1 x SVDClusterSize_2 = " + std::to_string(i);
227 h_U_Cl1Cl2_DeltaRes[i] = new TH1F(h_name_U, title_U, 100, -1000, 1000);
228 h_U_Cl1Cl2_DeltaRes[i]->GetXaxis()->SetTitle("#Delta res_{u} (#mum)");
229 h_U_Cl1Cl2_DeltaRes[i]->GetYaxis()->SetTitle("counts");
230 h_V_Cl1Cl2_DeltaRes[i] = new TH1F(h_name_V, title_V, 100, -1000, 1000);
231 h_V_Cl1Cl2_DeltaRes[i]->GetXaxis()->SetTitle("#Delta res_{v} (#mum)");
232 h_V_Cl1Cl2_DeltaRes[i]->GetYaxis()->SetTitle("counts");
233 }
234
235 //Special case of ExpertLevel option enabled
236 if (m_ExpertLevel) {
237 //Create directory to store PXD and SVD hitmaps for overlapping hits
238 TDirectory* HMDir = oldDir->mkdir("HitMaps_VXDOverlaps");
239 HMDir->cd();
240 //Define 2D sensor hit-maps for reconstructed hits
241 for (int i = 1; i <= 5; i++) {
242 for (int j = 1; j <= 16; j++) {
243 TString h_name = "h_6" + std::to_string(j) + std::to_string(i);
244 TString title = "Layer:Ladder:Sensor = 6:" + std::to_string(j) + ":" + std::to_string(i);
245 h_Lyr6[j][i] = new TH2F(h_name, title, 100, -2.88, 2.88, 100, -6.14, 6.14);
246 h_Lyr6[j][i]->GetXaxis()->SetTitle("u (cm)");
247 h_Lyr6[j][i]->GetYaxis()->SetTitle("v (cm)");
248 }
249 }
250 for (int i = 1; i <= 4; i++) {
251 for (int j = 1; j <= 12; j++) {
252 TString h_name = "h_5" + std::to_string(j) + std::to_string(i);
253 TString title = "Layer:Ladder:Sensor = 5:" + std::to_string(j) + ":" + std::to_string(i);
254 h_Lyr5[j][i] = new TH2F(h_name, title, 100, -2.88, 2.88, 100, -6.14, 6.14);
255 h_Lyr5[j][i]->GetXaxis()->SetTitle("u (cm)");
256 h_Lyr5[j][i]->GetYaxis()->SetTitle("v (cm)");
257 }
258 }
259 for (int i = 1; i <= 3; i++) {
260 for (int j = 1; j <= 10; j++) {
261 TString h_name = "h_4" + std::to_string(j) + std::to_string(i);
262 TString title = "Layer:Ladder:Sensor = 4:" + std::to_string(j) + ":" + std::to_string(i);
263 h_Lyr4[j][i] = new TH2F(h_name, title, 100, -2.88, 2.88, 100, -6.14, 6.14);
264 h_Lyr4[j][i]->GetXaxis()->SetTitle("u (cm)");
265 h_Lyr4[j][i]->GetYaxis()->SetTitle("v (cm)");
266 }
267 }
268 for (int i = 1; i <= 2; i++) {
269 for (int j = 1; j <= 7; j++) {
270 TString h_name = "h_3" + std::to_string(j) + std::to_string(i);
271 TString title = "Layer:Ladder:Sensor = 3:" + std::to_string(j) + ":" + std::to_string(i);
272 h_Lyr3[j][i] = new TH2F(h_name, title, 100, -1.92, 1.92, 100, -6.14, 6.14);
273 h_Lyr3[j][i]->GetXaxis()->SetTitle("u (cm)");
274 h_Lyr3[j][i]->GetYaxis()->SetTitle("v (cm)");
275 }
276 }
277 for (int i = 1; i <= 2; i++) {
278 for (int j = 1; j <= 12; j++) {
279 TString h_name = "h_2" + std::to_string(j) + std::to_string(i);
280 TString title = "Layer:Ladder:Sensor = 2:" + std::to_string(j) + ":" + std::to_string(i);
281 h_Lyr2[j][i] = new TH2F(h_name, title, 100, -0.625, 0.625, 100, -3.072, 3.072);
282 h_Lyr2[j][i]->GetXaxis()->SetTitle("u (cm)");
283 h_Lyr2[j][i]->GetYaxis()->SetTitle("v (cm)");
284 }
285 }
286 for (int i = 1; i <= 2; i++) {
287 for (int j = 1; j <= 8; j++) {
288 TString h_name = "h_1" + std::to_string(j) + std::to_string(i);
289 TString title = "Layer:Ladder:Sensor = 1:" + std::to_string(j) + ":" + std::to_string(i);
290 h_Lyr1[j][i] = new TH2F(h_name, title, 100, -0.625, 0.625, 100, -2.24, 2.24);
291 h_Lyr1[j][i]->GetXaxis()->SetTitle("u (cm)");
292 h_Lyr1[j][i]->GetYaxis()->SetTitle("v (cm)");
293 }
294 }
295 //Create directory to store PXD and SVD trees
296 TDirectory* TreeDir = oldDir->mkdir("Trees_VXDOverlaps");
297 TreeDir->cd();
298 //Tree for PXD
299 t_PXD = new TTree("t_PXD", "Tree for PXD overlaps");
300 t_PXD->Branch("deltaResU_PXD", &deltaResU_PXD, "deltaResU_PXD/F");
301 t_PXD->Branch("intResU_PXD", &intResU_PXD, "intResU_PXD/F");
302 t_PXD->Branch("intResV_PXD", &intResV_PXD, "intResV_PXD/F");
303 t_PXD->Branch("intU_PXD", &intU_PXD, "intU_PXD/F");
304 t_PXD->Branch("intV_PXD", &intV_PXD, "intV_PXD/F");
305 t_PXD->Branch("intPhi_PXD", &intPhi_PXD, "intPhi_PXD/F");
306 t_PXD->Branch("intZ_PXD", &intZ_PXD, "intZ_PXD/F");
307 t_PXD->Branch("intLayer_PXD", &intLayer_PXD, "intLayer_PXD/i");
308 t_PXD->Branch("intLadder_PXD", &intLadder_PXD, "intLadder_PXD/i");
309 t_PXD->Branch("intSensor_PXD", &intSensor_PXD, "intSensor_PXD/i");
310 t_PXD->Branch("extResU_PXD", &extResU_PXD, "extResU_PXD/F");
311 t_PXD->Branch("extResV_PXD", &extResV_PXD, "extResV_PXD/F");
312 t_PXD->Branch("extU_PXD", &extU_PXD, "extU_PXD/F");
313 t_PXD->Branch("extV_PXD", &extV_PXD, "extV_PXD/F");
314 t_PXD->Branch("extPhi_PXD", &extPhi_PXD, "extPhi_PXD/F");
315 t_PXD->Branch("extZ_PXD", &extZ_PXD, "extZ_PXD/F");
316 t_PXD->Branch("extLayer_PXD", &extLayer_PXD, "extLayer_PXD/i");
317 t_PXD->Branch("extLadder_PXD", &extLadder_PXD, "extLadder_PXD/i");
318 t_PXD->Branch("extSensor_PXD", &extSensor_PXD, "extSensor_PXD/i");
319 //Tree for SVD u overlapping clusters
320 t_SVD_U = new TTree("t_SVD_U", "Tree for SVD u-overlaps");
321 t_SVD_U->Branch("svdDeltaRes", &svdDeltaRes_U, "svdDeltaRes/F");
322 t_SVD_U->Branch("svdTrkPXDHits", &svdTrkPXDHits, "svdTrkPXDHits/i");
323 t_SVD_U->Branch("svdTrkSVDHits", &svdTrkSVDHits, "svdTrkSVDHits/i");
324 t_SVD_U->Branch("svdTrkCDCHits", &svdTrkCDCHits, "svdTrkCDCHits/i");
325 t_SVD_U->Branch("svdTrkd0", &svdTrkd0, "svdTrkd0/F");
326 t_SVD_U->Branch("svdTrkz0", &svdTrkz0, "svdTrkz0/F");
327 t_SVD_U->Branch("svdTrkpT", &svdTrkpT, "svdTrkpT/F");
328 t_SVD_U->Branch("svdTrkpCM", &svdTrkpCM, "svdTrkpCM/F");
329 // Internal ladder variables
330 t_SVD_U->Branch("svdClSNR_int", &svdClSNR_U_int, "svdClSNR_int/F");
331 t_SVD_U->Branch("svdClCharge_int", &svdClCharge_U_int, "svdClCharge_int/F");
332 t_SVD_U->Branch("svdStripCharge_int", &svdStripCharge_U_int);
333 t_SVD_U->Branch("svdStrip6Samples_int", &svdStrip6Samples_U_int);
334 t_SVD_U->Branch("svdClTime_int", &svdClTime_U_int, "svdClTime_int/F");
335 t_SVD_U->Branch("svdStripTime_int", &svdStripTime_U_int);
336 t_SVD_U->Branch("svdStripPosition_int", &svdStripPosition_U_int);
337 t_SVD_U->Branch("svdRes_int", &svdRes_U_int, "svdRes_int/F");
338 t_SVD_U->Branch("svdClIntStrPos_int", &svdClIntStrPos_U_int, "svdClIntStrPos_int/F");
339 t_SVD_U->Branch("svdClPos_int", &svdClPos_U_int, "svdClPos_int/F");
340 t_SVD_U->Branch("svdClPosErr_int", &svdClPosErr_U_int, "svdClPosErr_int/F");
341 t_SVD_U->Branch("svdTruePos_int", &svdTruePos_U_int, "svdTruePos_int/F");
342 t_SVD_U->Branch("svdClPhi_int", &svdClPhi_U_int, "svdClPhi_int/F");
343 t_SVD_U->Branch("svdClZ_int", &svdClZ_U_int, "svdClZ_int/F");
344 t_SVD_U->Branch("svdTrkTraversedLength_int", &svdTrkTraversedLength_U_int, "svdTrkTraversedLength_int/F");
345 t_SVD_U->Branch("svdTrkPos_int", &svdTrkPos_U_int, "svdTrkPos_int/F");
346 t_SVD_U->Branch("svdTrkPosOS_int", &svdTrkPosOS_U_int, "svdTrkPosOS_int/F");
347 t_SVD_U->Branch("svdTrkPosErr_int", &svdTrkPosErr_U_int, "svdTrkPosErr_int/F");
348 t_SVD_U->Branch("svdTrkPosErrOS_int", &svdTrkPosErrOS_U_int, "svdTrkPosErrOS_int/F");
349 t_SVD_U->Branch("svdTrkQoP_int", &svdTrkQoP_U_int, "svdTrkQoP_int/F");
350 t_SVD_U->Branch("svdTrkPrime_int", &svdTrkPrime_U_int, "svdTrkPrime_int/F");
351 t_SVD_U->Branch("svdTrkPrimeOS_int", &svdTrkPrimeOS_U_int, "svdTrkPrimeOS_int/F");
352 t_SVD_U->Branch("svdTrkPosUnbiased_int", &svdTrkPosUnbiased_U_int, "svdTrkPosUnbiased_int/F");
353 t_SVD_U->Branch("svdTrkPosErrUnbiased_int", &svdTrkPosErrUnbiased_U_int, "svdTrkPosErrUnbiased_int/F");
354 t_SVD_U->Branch("svdTrkQoPUnbiased_int", &svdTrkQoPUnbiased_U_int, "svdTrkQoPUnbiased_int/F");
355 t_SVD_U->Branch("svdTrkPrimeUnbiased_int", &svdTrkPrimeUnbiased_U_int, "svdTrkPrimeUnbiased_int/F");
356 t_SVD_U->Branch("svdLayer_int", &svdLayer_U_int, "svdLayer_int/i");
357 t_SVD_U->Branch("svdLadder_int", &svdLadder_U_int, "svdLadder_int/i");
358 t_SVD_U->Branch("svdSensor_int", &svdSensor_U_int, "svdSensor_int/i");
359 t_SVD_U->Branch("svdSize_int", &svdSize_U_int, "svdSize_int/i");
360 // External ladder variables
361 t_SVD_U->Branch("svdClSNR_ext", &svdClSNR_U_ext, "svdClSNR_ext/F");
362 t_SVD_U->Branch("svdClCharge_ext", &svdClCharge_U_ext, "svdClCharge_ext/F");
363 t_SVD_U->Branch("svdStripCharge_ext", &svdStripCharge_U_ext);
364 t_SVD_U->Branch("svdStrip6Samples_ext", &svdStrip6Samples_U_ext);
365 t_SVD_U->Branch("svdClTime_ext", &svdClTime_U_ext, "svdClTime_ext/F");
366 t_SVD_U->Branch("svdStripTime_ext", &svdStripTime_U_ext);
367 t_SVD_U->Branch("svdStripPosition_ext", &svdStripPosition_U_ext);
368 t_SVD_U->Branch("svdRes_ext", &svdRes_U_ext, "svdRes_ext/F");
369 t_SVD_U->Branch("svdClIntStrPos_ext", &svdClIntStrPos_U_ext, "svdClIntStrPos_ext/F");
370 t_SVD_U->Branch("svdClPos_ext", &svdClPos_U_ext, "svdClPos_ext/F");
371 t_SVD_U->Branch("svdClPosErr_ext", &svdClPosErr_U_ext, "svdClPosErr_ext/F");
372 t_SVD_U->Branch("svdTruePos_ext", &svdTruePos_U_ext, "svdTruePos_ext/F");
373 t_SVD_U->Branch("svdClPhi_ext", &svdClPhi_U_ext, "svdClPhi_ext/F");
374 t_SVD_U->Branch("svdClZ_ext", &svdClZ_U_ext, "svdClZ_ext/F");
375 t_SVD_U->Branch("svdTrkTraversedLength_ext", &svdTrkTraversedLength_U_ext, "svdTrkTraversedLength_ext/F");
376 t_SVD_U->Branch("svdTrkPos_ext", &svdTrkPos_U_ext, "svdTrkPos_ext/F");
377 t_SVD_U->Branch("svdTrkPosOS_ext", &svdTrkPosOS_U_ext, "svdTrkPosOS_ext/F");
378 t_SVD_U->Branch("svdTrkPosErr_ext", &svdTrkPosErr_U_ext, "svdTrkPosErr_ext/F");
379 t_SVD_U->Branch("svdTrkPosErrOS_ext", &svdTrkPosErrOS_U_ext, "svdTrkPosErrOS_ext/F");
380 t_SVD_U->Branch("svdTrkQoP_ext", &svdTrkQoP_U_ext, "svdTrkQoP_ext/F");
381 t_SVD_U->Branch("svdTrkPrime_ext", &svdTrkPrime_U_ext, "svdTrkPrime_ext/F");
382 t_SVD_U->Branch("svdTrkPrimeOS_ext", &svdTrkPrimeOS_U_ext, "svdTrkPrimeOS_ext/F");
383 t_SVD_U->Branch("svdTrkPosUnbiased_ext", &svdTrkPosUnbiased_U_ext, "svdTrkPosUnbiased_ext/F");
384 t_SVD_U->Branch("svdTrkPosErrUnbiased_ext", &svdTrkPosErrUnbiased_U_ext, "svdTrkPosErrUnbiased_ext/F");
385 t_SVD_U->Branch("svdTrkQoPUnbiased_ext", &svdTrkQoPUnbiased_U_ext, "svdTrkQoPUnbiased_ext/F");
386 t_SVD_U->Branch("svdTrkPrimeUnbiased_ext", &svdTrkPrimeUnbiased_U_ext, "svdTrkPrimeUnbiased_ext/F");
387 t_SVD_U->Branch("svdLayer_ext", &svdLayer_U_ext, "svdLayer_ext/i");
388 t_SVD_U->Branch("svdLadder_ext", &svdLadder_U_ext, "svdLadder_ext/i");
389 t_SVD_U->Branch("svdSensor_ext", &svdSensor_U_ext, "svdSensor_ext/i");
390 t_SVD_U->Branch("svdSize_ext", &svdSize_U_ext, "svdSize_ext/i");
391 //Tree for SVD v overlapping clusters
392 t_SVD_V = new TTree("t_SVD_V", "Tree for SVD v-overlaps");
393 t_SVD_V->Branch("svdDeltaRes", &svdDeltaRes_V, "svdDeltaRes/F");
394 t_SVD_V->Branch("svdTrkPXDHits", &svdTrkPXDHits, "svdTrkPXDHits/i");
395 t_SVD_V->Branch("svdTrkSVDHits", &svdTrkSVDHits, "svdTrkSVDHits/i");
396 t_SVD_V->Branch("svdTrkCDCHits", &svdTrkCDCHits, "svdTrkCDCHits/i");
397 t_SVD_V->Branch("svdTrkd0", &svdTrkd0, "svdTrkd0/F");
398 t_SVD_V->Branch("svdTrkz0", &svdTrkz0, "svdTrkz0/F");
399 t_SVD_V->Branch("svdTrkpT", &svdTrkpT, "svdTrkpT/F");
400 t_SVD_V->Branch("svdTrkpCM", &svdTrkpCM, "svdTrkpCM/F");
401 // Internal ladder variables
402 t_SVD_V->Branch("svdClSNR_int", &svdClSNR_V_int, "svdClSNR_int/F");
403 t_SVD_V->Branch("svdClCharge_int", &svdClCharge_V_int, "svdClCharge_int/F");
404 t_SVD_V->Branch("svdStripCharge_int", &svdStripCharge_V_int);
405 t_SVD_V->Branch("svdStrip6Samples_int", &svdStrip6Samples_V_int);
406 t_SVD_V->Branch("svdClTime_int", &svdClTime_V_int, "svdClTime_int/F");
407 t_SVD_V->Branch("svdStripTime_int", &svdStripTime_V_int);
408 t_SVD_V->Branch("svdStripPosition_int", &svdStripPosition_V_int);
409 t_SVD_V->Branch("svdRes_int", &svdRes_V_int, "svdRes_int/F");
410 t_SVD_V->Branch("svdClIntStrPos_int", &svdClIntStrPos_V_int, "svdClIntStrPos_int/F");
411 t_SVD_V->Branch("svdClPos_int", &svdClPos_V_int, "svdClPos_int/F");
412 t_SVD_V->Branch("svdClPosErr_int", &svdClPosErr_V_int, "svdClPosErr_int/F");
413 t_SVD_V->Branch("svdTruePos_int", &svdTruePos_V_int, "svdTruePos_int/F");
414 t_SVD_V->Branch("svdClPhi_int", &svdClPhi_V_int, "svdClPhi_int/F");
415 t_SVD_V->Branch("svdClZ_int", &svdClZ_V_int, "svdClZ_int/F");
416 t_SVD_V->Branch("svdTrkTraversedLength_int", &svdTrkTraversedLength_V_int, "svdTrkTraversedLength_int/F");
417 t_SVD_V->Branch("svdTrkPos_int", &svdTrkPos_V_int, "svdTrkPos_int/F");
418 t_SVD_V->Branch("svdTrkPosOS_int", &svdTrkPosOS_V_int, "svdTrkPosOS_int/F");
419 t_SVD_V->Branch("svdTrkPosErr_int", &svdTrkPosErr_V_int, "svdTrkPosErr_int/F");
420 t_SVD_V->Branch("svdTrkPosErrOS_int", &svdTrkPosErrOS_V_int, "svdTrkPosErrOS_int/F");
421 t_SVD_V->Branch("svdTrkQoP_int", &svdTrkQoP_V_int, "svdTrkQoP_int/F");
422 t_SVD_V->Branch("svdTrkPrime_int", &svdTrkPrime_V_int, "svdTrkPrime_int/F");
423 t_SVD_V->Branch("svdTrkPrimeOS_int", &svdTrkPrimeOS_V_int, "svdTrkPrimeOS_int/F");
424 t_SVD_V->Branch("svdTrkPosUnbiased_int", &svdTrkPosUnbiased_V_int, "svdTrkPosUnbiased_int/F");
425 t_SVD_V->Branch("svdTrkPosErrUnbiased_int", &svdTrkPosErrUnbiased_V_int, "svdTrkPosErrUnbiased_int/F");
426 t_SVD_V->Branch("svdTrkQoPUnbiased_int", &svdTrkQoPUnbiased_V_int, "svdTrkQoPUnbiased_int/F");
427 t_SVD_V->Branch("svdTrkPrimeUnbiased_int", &svdTrkPrimeUnbiased_V_int, "svdTrkPrimeUnbiased_int/F");
428 t_SVD_V->Branch("svdLayer_int", &svdLayer_V_int, "svdLayer_int/i");
429 t_SVD_V->Branch("svdLadder_int", &svdLadder_V_int, "svdLadder_int/i");
430 t_SVD_V->Branch("svdSensor_int", &svdSensor_V_int, "svdSensor_int/i");
431 t_SVD_V->Branch("svdSize_int", &svdSize_V_int, "svdSize_int/i");
432 // External ladder variables
433 t_SVD_V->Branch("svdClSNR_ext", &svdClSNR_V_ext, "svdClSNR_ext/F");
434 t_SVD_V->Branch("svdClCharge_ext", &svdClCharge_V_ext, "svdClCharge_ext/F");
435 t_SVD_V->Branch("svdStripCharge_ext", &svdStripCharge_V_ext);
436 t_SVD_V->Branch("svdStrip6Samples_ext", &svdStrip6Samples_V_ext);
437 t_SVD_V->Branch("svdClTime_ext", &svdClTime_V_ext, "svdClTime_ext/F");
438 t_SVD_V->Branch("svdStripTime_ext", &svdStripTime_V_ext);
439 t_SVD_V->Branch("svdStripPosition_ext", &svdStripPosition_V_ext);
440 t_SVD_V->Branch("svdRes_ext", &svdRes_V_ext, "svdRes_ext/F");
441 t_SVD_V->Branch("svdClIntStrPos_ext", &svdClIntStrPos_V_ext, "svdClIntStrPos_ext/F");
442 t_SVD_V->Branch("svdClPos_ext", &svdClPos_V_ext, "svdClPos_ext/F");
443 t_SVD_V->Branch("svdClPosErr_ext", &svdClPosErr_V_ext, "svdClPosErr_ext/F");
444 t_SVD_V->Branch("svdTruePos_ext", &svdTruePos_V_ext, "svdTruePos_ext/F");
445 t_SVD_V->Branch("svdClPhi_ext", &svdClPhi_V_ext, "svdClPhi_ext/F");
446 t_SVD_V->Branch("svdClZ_ext", &svdClZ_V_ext, "svdClZ_ext/F");
447 t_SVD_V->Branch("svdTrkTraversedLength_ext", &svdTrkTraversedLength_V_ext, "svdTrkTraversedLength_ext/F");
448 t_SVD_V->Branch("svdTrkPos_ext", &svdTrkPos_V_ext, "svdTrkPos_ext/F");
449 t_SVD_V->Branch("svdTrkPosOS_ext", &svdTrkPosOS_V_ext, "svdTrkPosOS_ext/F");
450 t_SVD_V->Branch("svdTrkPosErr_ext", &svdTrkPosErr_V_ext, "svdTrkPosErr_ext/F");
451 t_SVD_V->Branch("svdTrkPosErrOS_ext", &svdTrkPosErrOS_V_ext, "svdTrkPosErrOS_ext/F");
452 t_SVD_V->Branch("svdTrkQoP_ext", &svdTrkQoP_V_ext, "svdTrkQoP_ext/F");
453 t_SVD_V->Branch("svdTrkPrime_ext", &svdTrkPrime_V_ext, "svdTrkPrime_ext/F");
454 t_SVD_V->Branch("svdTrkPrimeOS_ext", &svdTrkPrimeOS_V_ext, "svdTrkPrimeOS_ext/F");
455 t_SVD_V->Branch("svdTrkPosUnbiased_ext", &svdTrkPosUnbiased_V_ext, "svdTrkPosUnbiased_ext/F");
456 t_SVD_V->Branch("svdTrkPosErrUnbiased_ext", &svdTrkPosErrUnbiased_V_ext, "svdTrkPosErrUnbiased_ext/F");
457 t_SVD_V->Branch("svdTrkQoPUnbiased_ext", &svdTrkQoPUnbiased_V_ext, "svdTrkQoPUnbiased_ext/F");
458 t_SVD_V->Branch("svdTrkPrimeUnbiased_ext", &svdTrkPrimeUnbiased_V_ext, "svdTrkPrimeUnbiased_ext/F");
459 t_SVD_V->Branch("svdLayer_ext", &svdLayer_V_ext, "svdLayer_ext/i");
460 t_SVD_V->Branch("svdLadder_ext", &svdLadder_V_ext, "svdLadder_ext/i");
461 t_SVD_V->Branch("svdSensor_ext", &svdSensor_V_ext, "svdSensor_ext/i");
462 t_SVD_V->Branch("svdSize_ext", &svdSize_V_ext, "svdSize_ext/i");
463 }
464 //Go back to the default output directory
465 oldDir->cd();
466}
467
468void OverlapResidualsModule::event()
469{
470 bool isMC = Environment::Instance().isMC();
471
472 static VXD::GeoCache& geo = VXD::GeoCache::getInstance();
473 StoreArray<RecoTrack> recoTracks(m_recoTracksStoreArrayName);
474 for (const auto& trk : recoTracks) {
475 if (! trk.wasFitSuccessful()) {
476 continue;
477 }
478 const std::vector<PXDCluster* > pxdClusters = trk.getPXDHitList();
479 const std::vector<SVDCluster* > svdClusters = trk.getSVDHitList();
480 B2DEBUG(29, "FITTED TRACK: NUMBER OF PXD HITS = " << pxdClusters.size() << " NUMBER OF SVD HITS = " << svdClusters.size());
481 //LOOKING FOR 2 CONSECUTIVE PXD HITS IN OVERLAPPING MODULES OF A SAME LAYER
482 for (unsigned int i = 0; i < pxdClusters.size(); i++) {
483 const PXDCluster* pxd_1 = pxdClusters[i];
484 const RecoHitInformation* infoPXD_1 = trk.getRecoHitInformation(pxd_1);
485 if (!infoPXD_1) {
486 continue;
487 }
488 const auto* hitTrackPoint_1 = trk.getCreatedTrackPoint(infoPXD_1);
489 const auto* fittedResult_1 = hitTrackPoint_1->getFitterInfo();
490 if (!fittedResult_1) {
491 continue;
492 }
493 const VxdID pxd_id_1 = pxd_1->getSensorID();
494 const unsigned short pxd_Layer_1 = pxd_id_1.getLayerNumber();
495 const unsigned short pxd_Ladder_1 = pxd_id_1.getLadderNumber();
496 const unsigned short pxd_Sensor_1 = pxd_id_1.getSensorNumber();
497 for (unsigned int l = i + 1; l < pxdClusters.size(); l++) {
498 const PXDCluster* pxd_2 = pxdClusters[l];
499 const RecoHitInformation* infoPXD_2 = trk.getRecoHitInformation(pxd_2);
500 if (!infoPXD_2) {
501 continue;
502 }
503 const auto* hitTrackPoint_2 = trk.getCreatedTrackPoint(infoPXD_2);
504 const auto* fittedResult_2 = hitTrackPoint_2->getFitterInfo();
505 if (!fittedResult_2) {
506 continue;
507 }
508 const VxdID pxd_id_2 = pxd_2->getSensorID();
509 const unsigned short pxd_Layer_2 = pxd_id_2.getLayerNumber();
510 const unsigned short pxd_Ladder_2 = pxd_id_2.getLadderNumber();
511 const unsigned short pxd_Sensor_2 = pxd_id_2.getSensorNumber();
512 if (pxd_Layer_1 == pxd_Layer_2 && ((pxd_Ladder_2 - pxd_Ladder_1) == -1.0 || (pxd_Ladder_2 - pxd_Ladder_1) == 7.0
513 || (pxd_Ladder_2 - pxd_Ladder_1) == 11.0)) {
514 B2DEBUG(29, " ============= 2 hits in a PXD overlap ============= ");
515 const TVectorD resUnBias_PXD_1 = fittedResult_1->getResidual(0, false).getState();
516 const TVectorD resUnBias_PXD_2 = fittedResult_2->getResidual(0, false).getState();
517 const float res_U_1 = resUnBias_PXD_1.GetMatrixArray()[0] * Unit::convertValueToUnit(1.0, "um");
518 const float res_V_1 = resUnBias_PXD_1.GetMatrixArray()[1] * Unit::convertValueToUnit(1.0, "um");
519 const double res_U_2 = resUnBias_PXD_2.GetMatrixArray()[0] * Unit::convertValueToUnit(1.0, "um");
520 const double res_V_2 = resUnBias_PXD_2.GetMatrixArray()[1] * Unit::convertValueToUnit(1.0, "um");
521 const float over_U_PXD = res_U_2 - res_U_1;
522 const float over_V_PXD = res_V_2 - res_V_1;
523 const ROOT::Math::XYZVector pxdLocal_1(pxd_1->getU(), pxd_1->getV(), 0.);
524 const ROOT::Math::XYZVector pxdLocal_2(pxd_2->getU(), pxd_2->getV(), 0.);
525 const VXD::SensorInfoBase& pxdSensor_1 = geo.getSensorInfo(pxd_id_1);
526 const VXD::SensorInfoBase& pxdSensor_2 = geo.getSensorInfo(pxd_id_2);
527 const ROOT::Math::XYZVector& pxdGlobal_1 = pxdSensor_1.pointToGlobal(pxdLocal_1);
528 const ROOT::Math::XYZVector& pxdGlobal_2 = pxdSensor_2.pointToGlobal(pxdLocal_2);
529 double pxdPhi_1 = atan2(pxdGlobal_1.Y(), pxdGlobal_1.X()); // maybe use pxdGlobal_1.Phi() instead
530 double pxdPhi_2 = atan2(pxdGlobal_2.Y(), pxdGlobal_2.X()); // maybe use pxdGlobal_2.Phi() instead
531 double pxdZ_1 = pxdGlobal_1.Z();
532 double pxdZ_2 = pxdGlobal_2.Z();
533 B2DEBUG(29, "PXD: difference of residuals " << over_U_PXD << " " << over_V_PXD);
534 //Fill PXD tree for overlaps if required by the user
535 if (m_ExpertLevel) {
536 deltaResU_PXD = over_U_PXD;
537 intResU_PXD = res_U_1;
538 intResV_PXD = res_V_1;
539 intU_PXD = pxd_1->getU();
540 intV_PXD = pxd_1->getV();
541 intPhi_PXD = pxdPhi_1;
542 intZ_PXD = pxdZ_1;
543 intLayer_PXD = pxd_Layer_1;
544 intLadder_PXD = pxd_Ladder_1;
545 intSensor_PXD = pxd_Sensor_1;
546 extResU_PXD = res_U_2;
547 extResV_PXD = res_V_2;
548 extU_PXD = pxd_2->getU();
549 extV_PXD = pxd_2->getV();
550 extPhi_PXD = pxdPhi_2;
551 extZ_PXD = pxdZ_2;
552 extLayer_PXD = pxd_Layer_2;
553 extLadder_PXD = pxd_Ladder_2;
554 extSensor_PXD = pxd_Sensor_2;
555 t_PXD->Fill();
556 }
557 //Fill histograms of residuals differences with PXD clusters
558 h_U_DeltaRes->Fill(over_U_PXD);
559 h_V_DeltaRes->Fill(over_V_PXD);
560 h_U_DeltaRes_PXD->Fill(over_U_PXD);
561 h_V_DeltaRes_PXD->Fill(over_V_PXD);
562 //Fill sensor hit-maps and 2D histograms with PXD clusters
563 if (pxd_Layer_1 == 1 && pxd_Layer_2 == 1) {
564 if (m_ExpertLevel) {
565 h_Lyr1[pxd_Ladder_1][pxd_Sensor_1]->Fill(pxd_1->getU(), pxd_1->getV());
566 h_Lyr1[pxd_Ladder_2][pxd_Sensor_2]->Fill(pxd_2->getU(), pxd_2->getV());
567 }
568 h_U_DeltaRes_PXD_Lyr1->Fill(over_U_PXD);
569 h_V_DeltaRes_PXD_Lyr1->Fill(over_V_PXD);
570 h_DeltaResUPhi_Lyr1->Fill(pxdPhi_1, over_U_PXD);
571 h_DeltaResUPhi_Lyr1->Fill(pxdPhi_2, over_U_PXD);
572 h_DeltaResVPhi_Lyr1->Fill(pxdPhi_1, over_V_PXD);
573 h_DeltaResVPhi_Lyr1->Fill(pxdPhi_2, over_V_PXD);
574 h_DeltaResVz_Lyr1->Fill(pxdZ_1, over_V_PXD);
575 h_DeltaResVz_Lyr1->Fill(pxdZ_2, over_V_PXD);
576 h_DeltaResUz_Lyr1->Fill(pxdZ_1, over_U_PXD);
577 h_DeltaResUz_Lyr1->Fill(pxdZ_2, over_U_PXD);
578 }
579 if (pxd_Layer_1 == 2 && pxd_Layer_2 == 2) {
580 if (m_ExpertLevel) {
581 h_Lyr2[pxd_Ladder_1][pxd_Sensor_1]->Fill(pxd_1->getU(), pxd_1->getV());
582 h_Lyr2[pxd_Ladder_2][pxd_Sensor_2]->Fill(pxd_2->getU(), pxd_2->getV());
583 }
584 h_U_DeltaRes_PXD_Lyr2->Fill(over_U_PXD);
585 h_V_DeltaRes_PXD_Lyr2->Fill(over_V_PXD);
586 h_DeltaResUPhi_Lyr2->Fill(pxdPhi_1, over_U_PXD);
587 h_DeltaResUPhi_Lyr2->Fill(pxdPhi_2, over_U_PXD);
588 h_DeltaResVPhi_Lyr2->Fill(pxdPhi_1, over_V_PXD);
589 h_DeltaResVPhi_Lyr2->Fill(pxdPhi_2, over_V_PXD);
590 h_DeltaResVz_Lyr2->Fill(pxdZ_1, over_V_PXD);
591 h_DeltaResVz_Lyr2->Fill(pxdZ_2, over_V_PXD);
592 h_DeltaResUz_Lyr2->Fill(pxdZ_1, over_U_PXD);
593 h_DeltaResUz_Lyr2->Fill(pxdZ_2, over_U_PXD);
594 }
595 }
596 }
597
598 }
599
600 //LOOKING FOR 2 CONSECUTIVE SVD HITS IN OVERLAPPING MODULES OF A SAME LAYER
601 RelationVector<Track> theTK = DataStore::getRelationsWithObj<Track>(&trk);
602
603 const TrackFitResult* tfr = theTK[0]->getTrackFitResultWithClosestMass(Const::pion);
604 if (tfr) {
605 svdTrkd0 = tfr->getD0();
606 svdTrkz0 = tfr->getZ0();
607 svdTrkpT = tfr->getMomentum().Rho();
608 ROOT::Math::PxPyPzEVector pStar = tfr->get4Momentum();
609 ROOT::Math::BoostZ boost(3. / 11);
610 pStar = boost(pStar);
611 svdTrkpCM = pStar.P();
612 }
613
614 svdTrkPXDHits = (trk.getPXDHitList()).size();
615 svdTrkSVDHits = (trk.getSVDHitList()).size();
616 svdTrkCDCHits = (trk.getCDCHitList()).size();
617
618 for (unsigned int i = 0; i < svdClusters.size(); i++) {
619 const SVDCluster* svd_1 = svdClusters[i];
620
621 //get true hits, used only if isMC
622 RelationVector<SVDTrueHit> trueHit_1 = DataStore::getRelationsWithObj<SVDTrueHit>(svd_1);
623
624 const RecoHitInformation* infoSVD_1 = trk.getRecoHitInformation(svd_1);
625 if (!infoSVD_1) {
626 continue;
627 }
628 const auto* hitTrackPoint_1 = trk.getCreatedTrackPoint(infoSVD_1);
629 const auto* fittedResult_1 = hitTrackPoint_1->getFitterInfo();
630 if (!fittedResult_1) {
631 continue;
632 }
633 const VxdID svd_id_1 = svd_1->getSensorID();
634 const unsigned short svd_Layer_1 = svd_id_1.getLayerNumber();
635 const unsigned short svd_Ladder_1 = svd_id_1.getLadderNumber();
636 const unsigned short svd_Sensor_1 = svd_id_1.getSensorNumber();
637 for (unsigned int l = i + 1; l < svdClusters.size(); l++) {
638 const SVDCluster* svd_2 = svdClusters[l];
639
640 //get true hits, used only if isMC
641 RelationVector<SVDTrueHit> trueHit_2 = DataStore::getRelationsWithObj<SVDTrueHit>(svd_2);
642
643 const RecoHitInformation* infoSVD_2 = trk.getRecoHitInformation(svd_2);
644 if (!infoSVD_2) {
645 continue;
646 }
647 const auto* hitTrackPoint_2 = trk.getCreatedTrackPoint(infoSVD_2);
648 const auto* fittedResult_2 = hitTrackPoint_2->getFitterInfo();
649 if (!fittedResult_2) {
650 continue;
651 }
652 const VxdID svd_id_2 = svd_2->getSensorID();
653 const unsigned short svd_Layer_2 = svd_id_2.getLayerNumber();
654 const unsigned short svd_Ladder_2 = svd_id_2.getLadderNumber();
655 const unsigned short svd_Sensor_2 = svd_id_2.getSensorNumber();
656 if (svd_Layer_1 == svd_Layer_2 && ((svd_Ladder_2 - svd_Ladder_1) == 1.0 || (svd_Ladder_2 - svd_Ladder_1) == -6.0
657 || (svd_Ladder_2 - svd_Ladder_1) == -9.0 || (svd_Ladder_2 - svd_Ladder_1) == -11.0 || (svd_Ladder_2 - svd_Ladder_1) == -15.0)) {
658 B2DEBUG(29, " ============= 2 hits in a SVD overlap ============= ");
659 const TVectorD resUnBias_SVD_1 = fittedResult_1->getResidual(0, false).getState();
660 const TVectorD resUnBias_SVD_2 = fittedResult_2->getResidual(0, false).getState();
661 genfit::MeasuredStateOnPlane state_unbiased_1 = fittedResult_1->getFittedState(false);
662 genfit::MeasuredStateOnPlane state_unbiased_2 = fittedResult_2->getFittedState(false);
663 const TVectorD& svd_predIntersect_unbiased_1 = state_unbiased_1.getState();
664 const TVectorD& svd_predIntersect_unbiased_2 = state_unbiased_2.getState();
665 const TMatrixDSym& covMatrix_unbiased_1 = state_unbiased_1.getCov();
666 const TMatrixDSym& covMatrix_unbiased_2 = state_unbiased_2.getCov();
667 genfit::MeasuredStateOnPlane state_1 = trk.getMeasuredStateOnPlaneFromRecoHit(infoSVD_1);
668 genfit::MeasuredStateOnPlane state_2 = trk.getMeasuredStateOnPlaneFromRecoHit(infoSVD_2);
669 const TVectorD& svd_predIntersect_1 = state_1.getState();
670 const TVectorD& svd_predIntersect_2 = state_2.getState();
671 const TMatrixDSym& covMatrix_1 = state_1.getCov();
672 const TMatrixDSym& covMatrix_2 = state_2.getCov();
673 //Restricting to consecutive SVD u-clusters
674 if (svd_1->isUCluster() == true && svd_2->isUCluster() == true) {
675 const int strips_1 = svd_1->getSize();
676 h_SVDstrips_Mult->Fill(strips_1);
677 const int strips_2 = svd_2->getSize();
678 h_SVDstrips_Mult->Fill(strips_2);
679 const double res_U_1 = resUnBias_SVD_1.GetMatrixArray()[0] * Unit::convertValueToUnit(1.0, "um");
680 const double res_U_2 = resUnBias_SVD_2.GetMatrixArray()[0] * Unit::convertValueToUnit(1.0, "um");
681 const float over_U_SVD = res_U_2 - res_U_1;
682 const ROOT::Math::XYZVector svdLocal_1(svd_1->getPosition(), svd_predIntersect_1[4], 0.);
683 const ROOT::Math::XYZVector svdLocal_2(svd_2->getPosition(), svd_predIntersect_2[4], 0.);
684 const VXD::SensorInfoBase& svdSensor_1 = geo.getSensorInfo(svd_id_1);
685 const VXD::SensorInfoBase& svdSensor_2 = geo.getSensorInfo(svd_id_2);
686 const ROOT::Math::XYZVector& svdGlobal_1 = svdSensor_1.pointToGlobal(svdLocal_1);
687 const ROOT::Math::XYZVector& svdGlobal_2 = svdSensor_2.pointToGlobal(svdLocal_2);
688 double svdPhi_1 = atan2(svdGlobal_1.Y(), svdGlobal_1.X()); // maybe use svdGlobal_1.Phi() instead
689 double svdPhi_2 = atan2(svdGlobal_2.Y(), svdGlobal_2.X()); // maybe use svdGlobal_2.Phi() instead
690 double svdZ_1 = svdGlobal_1.Z();
691 double svdZ_2 = svdGlobal_2.Z();
692 B2DEBUG(29, "SVD: difference of u-residuals =========> " << over_U_SVD);
693 //Fill SVD tree for u-overlaps if required by the user
694 if (m_ExpertLevel) {
695 svdDeltaRes_U = over_U_SVD;
696 svdRes_U_int = res_U_1;
697 svdClTime_U_int = svd_1->getClsTime();
698 svdClSNR_U_int = svd_1->getSNR();
699 svdClCharge_U_int = svd_1->getCharge();
700 svdClPosErr_U_int = svd_1->getPositionSigma();
701 if (isMC && trueHit_1.size() > 0)
702 svdTruePos_U_int = trueHit_1[0]->getU();
703 else
704 svdTruePos_U_int = -99;
705 svdClPhi_U_int = svdPhi_1;
706 svdClZ_U_int = svdZ_1;
707 svdTrkPos_U_int = svd_predIntersect_1[3];
708 svdTrkPosOS_U_int = svd_predIntersect_1[4];
709 svdTrkPosErr_U_int = sqrt(covMatrix_1[3][3]);
710 svdTrkPosErrOS_U_int = sqrt(covMatrix_1[4][4]);
711 svdTrkQoP_U_int = svd_predIntersect_1[0];
712 svdTrkPrime_U_int = svd_predIntersect_1[1];
713 svdTrkPrimeOS_U_int = svd_predIntersect_1[2];
714 svdTrkTraversedLength_U_int = svdSensor_1.getThickness() * sqrt(1 + svdTrkPrimeOS_U_int * svdTrkPrimeOS_U_int + svdTrkPrime_U_int *
715 svdTrkPrime_U_int);
716 svdTrkPosUnbiased_U_int = svd_predIntersect_unbiased_1[3];
717 svdClPos_U_int = svdRes_U_int / 1e4 + svdTrkPosUnbiased_U_int;
718 svdTrkPosErrUnbiased_U_int = sqrt(covMatrix_unbiased_1[3][3]);
719 svdTrkQoPUnbiased_U_int = svd_predIntersect_unbiased_1[0];
720 svdTrkPrimeUnbiased_U_int = svd_predIntersect_unbiased_1[1];
721 svdLayer_U_int = svd_Layer_1;
722 svdLadder_U_int = svd_Ladder_1;
723 svdSensor_U_int = svd_Sensor_1;
724 svdSize_U_int = strips_1;
725
726 float pitch = 50e-4;
727 float halfLength = 1.92;
728 if (svdLayer_U_int > 3) {
729 pitch = 75e-4;
730 halfLength = 2.88;
731 }
732 svdClIntStrPos_U_int = fmod(svdClPos_U_int + halfLength, pitch) / pitch;
733
734 svdStripCharge_U_int.clear();
735 svdStrip6Samples_U_int.clear();
736 svdStripTime_U_int.clear();
737 svdStripPosition_U_int.clear();
738 //retrieve relations and set strip charges and times
739 RelationVector<SVDRecoDigit> theRecoDigits_1 = DataStore::getRelationsWithObj<SVDRecoDigit>(svd_1);
740 if ((theRecoDigits_1.size() != svdSize_U_int) && (svdSize_U_int != 128)) //virtual cluster
741 B2ERROR(" Inconsistency with cluster size! # recoDigits = " << theRecoDigits_1.size() << " != " << svdSize_U_int <<
742 " cluster size");
743
744 //skip clusters created beacuse of missing APV
745 if (svdSize_U_int < 128)
746 for (unsigned int d = 0; d < svdSize_U_int; d++) {
747 svdStripCharge_U_int.push_back(theRecoDigits_1[d]->getCharge());
748 SVDShaperDigit* ShaperDigit_1 = theRecoDigits_1[d]->getRelated<SVDShaperDigit>();
749 std::array<float, 6> Samples_1 = ShaperDigit_1->getSamples();
750 std::copy(std::begin(Samples_1), std::end(Samples_1), std::back_inserter(svdStrip6Samples_U_int));
751 svdStripTime_U_int.push_back(theRecoDigits_1[d]->getTime());
752 double misalignedStripPos = svdSensor_1.getUCellPosition(theRecoDigits_1[d]->getCellID());
753 //aligned strip pos = misaligned strip - ( misaligned cluster - aligned cluster)
754 svdStripPosition_U_int.push_back(misalignedStripPos - svd_1->getPosition() + svdClPos_U_int);
755 }
756
757 svdRes_U_ext = res_U_2;
758 svdClTime_U_ext = svd_2->getClsTime();
759 svdClSNR_U_ext = svd_2->getSNR();
760 svdClCharge_U_ext = svd_2->getCharge();
761 svdClPosErr_U_ext = svd_2->getPositionSigma();
762 if (isMC && trueHit_2.size() > 0)
763 svdTruePos_U_ext = trueHit_2[0]->getU();
764 else
765 svdTruePos_U_ext = -99;
766 svdClPhi_U_ext = svdPhi_2;
767 svdClZ_U_ext = svdZ_2;
768 svdTrkPos_U_ext = svd_predIntersect_2[3];
769 svdTrkPosOS_U_ext = svd_predIntersect_2[4];
770 svdTrkPosErr_U_ext = sqrt(covMatrix_2[3][3]);
771 svdTrkPosErrOS_U_ext = sqrt(covMatrix_2[4][4]);
772 svdTrkQoP_U_ext = svd_predIntersect_2[0];
773 svdTrkPrime_U_ext = svd_predIntersect_2[1];
774 svdTrkPrimeOS_U_ext = svd_predIntersect_2[2];
775 svdTrkTraversedLength_U_ext = svdSensor_2.getThickness() * sqrt(1 + svdTrkPrimeOS_U_ext * svdTrkPrimeOS_U_ext + svdTrkPrime_U_ext *
776 svdTrkPrime_U_ext);
777 svdTrkPosUnbiased_U_ext = svd_predIntersect_unbiased_2[3];
778 svdClPos_U_ext = svdRes_U_ext / 1e4 + svdTrkPosUnbiased_U_ext;
779 svdTrkPosErrUnbiased_U_ext = sqrt(covMatrix_unbiased_2[3][3]);
780 svdTrkQoPUnbiased_U_ext = svd_predIntersect_unbiased_2[0];
781 svdTrkPrimeUnbiased_U_ext = svd_predIntersect_unbiased_2[1];
782 svdLayer_U_ext = svd_Layer_2;
783 svdLadder_U_ext = svd_Ladder_2;
784 svdSensor_U_ext = svd_Sensor_2;
785 svdSize_U_ext = strips_2;
786
787 svdClIntStrPos_U_ext = fmod(svdClPos_U_ext + halfLength, pitch) / pitch;
788
789 svdStripCharge_U_ext.clear();
790 svdStrip6Samples_U_ext.clear();
791 svdStripTime_U_ext.clear();
792 svdStripPosition_U_ext.clear();
793 //retrieve relations and set strip charges and times
794 RelationVector<SVDRecoDigit> theRecoDigits_2 = DataStore::getRelationsWithObj<SVDRecoDigit>(svd_2);
795 if ((theRecoDigits_2.size() != svdSize_U_ext) && (svdSize_U_ext != 128)) //virtual cluster
796 B2ERROR(" Inconsistency with cluster size! # recoDigits = " << theRecoDigits_2.size() << " != " << svdSize_U_ext <<
797 " cluster size");
798 //skip clusters created beacuse of missing APV
799 if (svdSize_U_ext < 128)
800 for (unsigned int d = 0; d < svdSize_U_ext; d++) {
801 svdStripCharge_U_ext.push_back(theRecoDigits_2[d]->getCharge());
802 SVDShaperDigit* ShaperDigit_2 = theRecoDigits_2[d]->getRelated<SVDShaperDigit>();
803 std::array<float, 6> Samples_2 = ShaperDigit_2->getSamples();
804 std::copy(std::begin(Samples_2), std::end(Samples_2), std::back_inserter(svdStrip6Samples_U_ext));
805 svdStripTime_U_ext.push_back(theRecoDigits_2[d]->getTime());
806 double misalignedStripPos = svdSensor_2.getUCellPosition(theRecoDigits_2[d]->getCellID());
807 //aligned strip pos = misaligned strip - ( misaligned cluster - aligned cluster)
808 svdStripPosition_U_ext.push_back(misalignedStripPos - svd_2->getPosition() + svdClPos_U_ext);
809 }
810 t_SVD_U->Fill();
811 }
812 //Fill histograms of residuals differences with SVD u clusters
813 h_U_DeltaRes->Fill(over_U_SVD);
814 h_U_DeltaRes_SVD->Fill(over_U_SVD);
815 if (svd_1->getSize() < 3 && svd_2->getSize() < 3) { //Consider only clusters with 3 or less strips involved
816 h_U_Cl1Cl2_DeltaRes[svd_1->getSize() * svd_2->getSize()]->Fill(over_U_SVD);
817 }
818 //Fill sensor hit-maps and 2D histograms with SVD u clusters
819 if (svd_Layer_1 == 3 && svd_Layer_2 == 3) {
820 if (m_ExpertLevel) {
821 h_Lyr3[svd_Ladder_1][svd_Sensor_1]->Fill(svd_1->getPosition(), 0.);
822 h_Lyr3[svd_Ladder_2][svd_Sensor_2]->Fill(svd_2->getPosition(), 0.);
823 }
824 h_U_DeltaRes_SVD_Lyr3->Fill(over_U_SVD);
825 h_DeltaResUPhi_Lyr3->Fill(svdPhi_1, over_U_SVD);
826 h_DeltaResUPhi_Lyr3->Fill(svdPhi_2, over_U_SVD);
827 h_DeltaResUz_Lyr3->Fill(svdZ_1, over_U_SVD);
828 h_DeltaResUz_Lyr3->Fill(svdZ_2, over_U_SVD);
829 }
830 if (svd_Layer_1 == 4 && svd_Layer_2 == 4) {
831 if (m_ExpertLevel) {
832 h_Lyr4[svd_Ladder_1][svd_Sensor_1]->Fill(svd_1->getPosition(), 0.);
833 h_Lyr4[svd_Ladder_2][svd_Sensor_2]->Fill(svd_2->getPosition(), 0.);
834 }
835 h_U_DeltaRes_SVD_Lyr4->Fill(over_U_SVD);
836 h_DeltaResUPhi_Lyr4->Fill(svdPhi_1, over_U_SVD);
837 h_DeltaResUPhi_Lyr4->Fill(svdPhi_2, over_U_SVD);
838 h_DeltaResUz_Lyr4->Fill(svdZ_1, over_U_SVD);
839 h_DeltaResUz_Lyr4->Fill(svdZ_2, over_U_SVD);
840 }
841 if (svd_Layer_1 == 5 && svd_Layer_2 == 5) {
842 if (m_ExpertLevel) {
843 h_Lyr5[svd_Ladder_1][svd_Sensor_1]->Fill(svd_1->getPosition(), 0.);
844 h_Lyr5[svd_Ladder_2][svd_Sensor_1]->Fill(svd_1->getPosition(), 0.);
845 }
846 h_U_DeltaRes_SVD_Lyr5->Fill(over_U_SVD);
847 h_DeltaResUPhi_Lyr5->Fill(svdPhi_1, over_U_SVD);
848 h_DeltaResUPhi_Lyr5->Fill(svdPhi_2, over_U_SVD);
849 h_DeltaResUz_Lyr5->Fill(svdZ_1, over_U_SVD);
850 h_DeltaResUz_Lyr5->Fill(svdZ_2, over_U_SVD);
851 }
852 if (svd_Layer_1 == 6 && svd_Layer_2 == 6) {
853 if (m_ExpertLevel) {
854 h_Lyr6[svd_Ladder_1][svd_Sensor_1]->Fill(svd_1->getPosition(), 0.);
855 h_Lyr6[svd_Ladder_2][svd_Sensor_2]->Fill(svd_2->getPosition(), 0.);
856 }
857 h_U_DeltaRes_SVD_Lyr6->Fill(over_U_SVD);
858 h_DeltaResUPhi_Lyr6->Fill(svdPhi_1, over_U_SVD);
859 h_DeltaResUPhi_Lyr6->Fill(svdPhi_2, over_U_SVD);
860 h_DeltaResUz_Lyr6->Fill(svdZ_1, over_U_SVD);
861 h_DeltaResUz_Lyr6->Fill(svdZ_2, over_U_SVD);
862 }
863 }
864 //Restricting to consecutive SVD v-clusters
865 if (svd_1->isUCluster() != true && svd_2->isUCluster() != true) {
866 const int strips_1 = svd_1->getSize();
867 h_SVDstrips_Mult->Fill(strips_1);
868 const int strips_2 = svd_2->getSize();
869 h_SVDstrips_Mult->Fill(strips_2);
870 const double res_V_1 = resUnBias_SVD_1.GetMatrixArray()[0] * Unit::convertValueToUnit(1.0, "um");
871 const double res_V_2 = resUnBias_SVD_2.GetMatrixArray()[0] * Unit::convertValueToUnit(1.0, "um");
872 const float over_V_SVD = res_V_2 - res_V_1;
873 const ROOT::Math::XYZVector svdLocal_1(svd_predIntersect_1[3], svd_1->getPosition(), 0.);
874 const ROOT::Math::XYZVector svdLocal_2(svd_predIntersect_2[3], svd_2->getPosition(), 0.);
875 const VXD::SensorInfoBase& svdSensor_1 = geo.getSensorInfo(svd_id_1);
876 const VXD::SensorInfoBase& svdSensor_2 = geo.getSensorInfo(svd_id_2);
877 const ROOT::Math::XYZVector& svdGlobal_1 = svdSensor_1.pointToGlobal(svdLocal_1);
878 const ROOT::Math::XYZVector& svdGlobal_2 = svdSensor_2.pointToGlobal(svdLocal_2);
879 double svdPhi_1 = atan2(svdGlobal_1.Y(), svdGlobal_1.X()); // maybe use svdGlobal_1.Phi() instead
880 double svdPhi_2 = atan2(svdGlobal_2.Y(), svdGlobal_2.X()); // maybe use svdGlobal_2.Phi() instead
881 double svdZ_1 = svdGlobal_1.Z();
882 double svdZ_2 = svdGlobal_2.Z();
883 B2DEBUG(29, "SVD: difference of v-residuals =========> " << over_V_SVD);
884 //Fill SVD tree for v-overlaps if required by the user
885 if (m_ExpertLevel) {
886 svdDeltaRes_V = over_V_SVD;
887 svdRes_V_int = res_V_1;
888 svdClTime_V_int = svd_1->getClsTime();
889 svdClSNR_V_int = svd_1->getSNR();
890 svdClCharge_V_int = svd_1->getCharge();
891 svdClPosErr_V_int = svd_1->getPositionSigma();
892 if (isMC && trueHit_1.size() > 0)
893 svdTruePos_V_int = trueHit_1[0]->getV();
894 else
895 svdTruePos_V_int = -99;
896 svdClPhi_V_int = svdPhi_1;
897 svdClZ_V_int = svdZ_1;
898 svdTrkPos_V_int = svd_predIntersect_1[4];
899 svdTrkPosOS_V_int = svd_predIntersect_1[3];
900 svdTrkPosErr_V_int = sqrt(covMatrix_1[4][4]);
901 svdTrkPosErrOS_V_int = sqrt(covMatrix_1[3][3]);
902 svdTrkQoP_V_int = svd_predIntersect_1[0];
903 svdTrkPrime_V_int = svd_predIntersect_1[2];
904 svdTrkPrimeOS_V_int = svd_predIntersect_1[1];
905 svdTrkTraversedLength_V_int = svdSensor_1.getThickness() * sqrt(1 + svdTrkPrimeOS_V_int * svdTrkPrimeOS_V_int + svdTrkPrime_V_int *
906 svdTrkPrime_V_int);
907 svdTrkPosUnbiased_V_int = svd_predIntersect_unbiased_1[4];
908 svdClPos_V_int = svdRes_V_int / 1e4 + svdTrkPosUnbiased_V_int;
909 svdTrkPosErrUnbiased_V_int = sqrt(covMatrix_unbiased_1[4][4]);
910 svdTrkQoPUnbiased_V_int = svd_predIntersect_unbiased_1[0];
911 svdTrkPrimeUnbiased_V_int = svd_predIntersect_unbiased_1[2];
912 svdLayer_V_int = svd_Layer_1;
913 svdLadder_V_int = svd_Ladder_1;
914 svdSensor_V_int = svd_Sensor_1;
915 svdSize_V_int = strips_1;
916
917 float pitch = 160e-4;
918 float halfLength = 6.144;
919 if (svdLayer_V_int > 3) {
920 pitch = 240e-4;
921 }
922 svdClIntStrPos_V_int = fmod(svdClPos_V_int + halfLength, pitch) / pitch;
923
924 svdStripCharge_V_int.clear();
925 svdStrip6Samples_V_int.clear();
926 svdStripTime_V_int.clear();
927 svdStripPosition_V_int.clear();
928 //retrieve relations and set strip charges and times
929 RelationVector<SVDRecoDigit> theRecoDigits_1 = DataStore::getRelationsWithObj<SVDRecoDigit>(svd_1);
930 if ((theRecoDigits_1.size() != svdSize_V_int) && (svdSize_V_int != 128)) //virtual cluster
931 B2ERROR(" Inconsistency with cluster size! # recoDigits = " << theRecoDigits_1.size() << " != " << svdSize_V_int <<
932 " cluster size");
933 //skip clusters created beacuse of missing APV
934 if (svdSize_V_int < 128)
935 for (unsigned int d = 0; d < svdSize_V_int; d++) {
936 svdStripCharge_V_int.push_back(theRecoDigits_1[d]->getCharge());
937 SVDShaperDigit* ShaperDigit_1 = theRecoDigits_1[d]->getRelated<SVDShaperDigit>();
938 std::array<float, 6> Samples_1 = ShaperDigit_1->getSamples();
939 std::copy(std::begin(Samples_1), std::end(Samples_1), std::back_inserter(svdStrip6Samples_V_int));
940 svdStripTime_V_int.push_back(theRecoDigits_1[d]->getTime());
941 double misalignedStripPos = svdSensor_1.getVCellPosition(theRecoDigits_1[d]->getCellID());
942 //aligned strip pos = misaligned strip - ( misaligned cluster - aligned cluster)
943 svdStripPosition_V_int.push_back(misalignedStripPos - svd_1->getPosition() + svdClPos_V_int);
944 }
945
946 svdRes_V_ext = res_V_2;
947 svdClTime_V_ext = svd_2->getClsTime();
948 svdClSNR_V_ext = svd_2->getSNR();
949 svdClCharge_V_ext = svd_2->getCharge();
950 svdClPosErr_V_ext = svd_2->getPositionSigma();
951 if (isMC && trueHit_2.size() > 0)
952 svdTruePos_V_ext = trueHit_2[0]->getV();
953 else
954 svdTruePos_V_ext = -99;
955 svdClPhi_V_ext = svdPhi_2;
956 svdClZ_V_ext = svdZ_2;
957 svdTrkPos_V_ext = svd_predIntersect_2[4];
958 svdTrkPosOS_V_ext = svd_predIntersect_2[3];
959 svdTrkPosErr_V_ext = sqrt(covMatrix_2[4][4]);
960 svdTrkPosErrOS_V_ext = sqrt(covMatrix_1[3][3]);
961 svdTrkQoP_V_ext = svd_predIntersect_2[0];
962 svdTrkPrime_V_ext = svd_predIntersect_2[2];
963 svdTrkPrimeOS_V_ext = svd_predIntersect_2[1];
964 svdTrkTraversedLength_V_ext = svdSensor_2.getThickness() * sqrt(1 + svdTrkPrimeOS_V_ext * svdTrkPrimeOS_V_ext + svdTrkPrime_V_ext *
965 svdTrkPrime_V_ext);
966 svdTrkPosUnbiased_V_ext = svd_predIntersect_unbiased_2[4];
967 svdClPos_V_ext = svdRes_V_ext / 1e4 + svdTrkPosUnbiased_V_ext;
968 svdTrkPosErrUnbiased_V_ext = sqrt(covMatrix_unbiased_2[4][4]);
969 svdTrkQoPUnbiased_V_ext = svd_predIntersect_unbiased_2[0];
970 svdTrkPrimeUnbiased_V_ext = svd_predIntersect_unbiased_2[2];
971 svdLayer_V_ext = svd_Layer_2;
972 svdLadder_V_ext = svd_Ladder_2;
973 svdSensor_V_ext = svd_Sensor_2;
974 svdSize_V_ext = strips_2;
975
976 svdClIntStrPos_V_ext = fmod(svdClPos_V_ext + halfLength, pitch) / pitch;
977
978 svdStripCharge_V_ext.clear();
979 svdStrip6Samples_V_ext.clear();
980 svdStripTime_V_ext.clear();
981 svdStripPosition_V_ext.clear();
982 //retrieve relations and set strip charges and times
983 RelationVector<SVDRecoDigit> theRecoDigits_2 = DataStore::getRelationsWithObj<SVDRecoDigit>(svd_2);
984 if ((theRecoDigits_2.size() != svdSize_V_ext) && (svdSize_V_ext != 128)) //virtual cluster
985 B2ERROR(" Inconsistency with cluster size! # recoDigits = " << theRecoDigits_2.size() << " != " << svdSize_V_ext <<
986 " cluster size");
987 //skip clusters created beacuse of missing APV
988 if (svdSize_V_ext < 128)
989 for (unsigned int d = 0; d < svdSize_V_ext; d++) {
990 svdStripCharge_V_ext.push_back(theRecoDigits_2[d]->getCharge());
991 SVDShaperDigit* ShaperDigit_2 = theRecoDigits_2[d]->getRelated<SVDShaperDigit>();
992 std::array<float, 6> Samples_2 = ShaperDigit_2->getSamples();
993 std::copy(std::begin(Samples_2), std::end(Samples_2), std::back_inserter(svdStrip6Samples_V_ext));
994 svdStripTime_V_ext.push_back(theRecoDigits_2[d]->getTime());
995 double misalignedStripPos = svdSensor_2.getVCellPosition(theRecoDigits_2[d]->getCellID());
996 //aligned strip pos = misaligned strip - ( misaligned cluster - aligned cluster)
997 svdStripPosition_V_ext.push_back(misalignedStripPos - svd_2->getPosition() + svdClPos_V_ext);
998 }
999
1000 t_SVD_V->Fill();
1001 }
1002 //Fill histograms of residuals differences with SVD v clusters
1003 h_V_DeltaRes->Fill(over_V_SVD);
1004 h_V_DeltaRes_SVD->Fill(over_V_SVD);
1005 if (svd_1->getSize() < 3 && svd_2->getSize() < 3) { //Consider only clusters with 3 or less strips involved
1006 h_V_Cl1Cl2_DeltaRes[svd_1->getSize() * svd_2->getSize()]->Fill(over_V_SVD);
1007 }
1008 //Fill sensor hit-maps and 2D histograms with SVD v clusters
1009 if (svd_Layer_1 == 3 && svd_Layer_2 == 3) {
1010 if (m_ExpertLevel) {
1011 h_Lyr3[svd_Ladder_1][svd_Sensor_1]->Fill(0., svd_1->getPosition());
1012 h_Lyr3[svd_Ladder_2][svd_Sensor_2]->Fill(0., svd_2->getPosition());
1013 }
1014 h_V_DeltaRes_SVD_Lyr3->Fill(over_V_SVD);
1015 h_DeltaResVPhi_Lyr3->Fill(svdPhi_1, over_V_SVD);
1016 h_DeltaResVPhi_Lyr3->Fill(svdPhi_2, over_V_SVD);
1017 h_DeltaResVz_Lyr3->Fill(svdZ_1, over_V_SVD);
1018 h_DeltaResVz_Lyr3->Fill(svdZ_2, over_V_SVD);
1019 }
1020 if (svd_Layer_1 == 4 && svd_Layer_2 == 4) {
1021 if (m_ExpertLevel) {
1022 h_Lyr4[svd_Ladder_1][svd_Sensor_1]->Fill(0., svd_1->getPosition());
1023 h_Lyr4[svd_Ladder_2][svd_Sensor_2]->Fill(0., svd_2->getPosition());
1024 }
1025 h_V_DeltaRes_SVD_Lyr4->Fill(over_V_SVD);
1026 h_DeltaResVPhi_Lyr4->Fill(svdPhi_1, over_V_SVD);
1027 h_DeltaResVPhi_Lyr4->Fill(svdPhi_2, over_V_SVD);
1028 h_DeltaResVz_Lyr4->Fill(svdZ_1, over_V_SVD);
1029 h_DeltaResVz_Lyr4->Fill(svdZ_2, over_V_SVD);
1030 }
1031 if (svd_Layer_1 == 5 && svd_Layer_2 == 5) {
1032 if (m_ExpertLevel) {
1033 h_Lyr5[svd_Ladder_1][svd_Sensor_1]->Fill(0., svd_1->getPosition());
1034 h_Lyr6[svd_Ladder_2][svd_Sensor_2]->Fill(0., svd_2->getPosition());
1035 }
1036 h_V_DeltaRes_SVD_Lyr5->Fill(over_V_SVD);
1037 h_DeltaResVPhi_Lyr5->Fill(svdPhi_1, over_V_SVD);
1038 h_DeltaResVPhi_Lyr5->Fill(svdPhi_2, over_V_SVD);
1039 h_DeltaResVz_Lyr5->Fill(svdZ_1, over_V_SVD);
1040 h_DeltaResVz_Lyr5->Fill(svdZ_2, over_V_SVD);
1041 }
1042 if (svd_Layer_1 == 6 && svd_Layer_2 == 6) {
1043 if (m_ExpertLevel) {
1044 h_Lyr6[svd_Ladder_1][svd_Sensor_1]->Fill(0., svd_1->getPosition());
1045 h_Lyr6[svd_Ladder_2][svd_Sensor_2]->Fill(0., svd_2->getPosition());
1046 }
1047 h_V_DeltaRes_SVD_Lyr6->Fill(over_V_SVD);
1048 h_DeltaResVPhi_Lyr6->Fill(svdPhi_1, over_V_SVD);
1049 h_DeltaResVPhi_Lyr6->Fill(svdPhi_2, over_V_SVD);
1050 h_DeltaResVz_Lyr6->Fill(svdZ_1, over_V_SVD);
1051 h_DeltaResVz_Lyr6->Fill(svdZ_2, over_V_SVD);
1052 }
1053 }
1054 }
1055 }
1056 }
1057 }
1058}
1059
1060
Belle2::Const::pion
static const ChargedStable pion
charged pion particle
Definition: Const.h:661
Belle2::Environment::isMC
bool isMC() const
Do we have generated, not real data?
Definition: Environment.cc:54
Belle2::Environment::Instance
static Environment & Instance()
Static method to get a reference to the Environment instance.
Definition: Environment.cc:28
Belle2::HistoModule
HistoModule.h is supposed to be used instead of Module.h for the modules with histogram definitions t...
Definition: HistoModule.h:29
Belle2::Module::setDescription
void setDescription(const std::string &description)
Sets the description of the module.
Definition: Module.cc:214
Belle2::OverlapResidualsModule::svdTrkQoP_V_int
float svdTrkQoP_V_int
V internal track q/p.
Definition: OverlapResidualsModule.h:266
Belle2::OverlapResidualsModule::svdTrkPosUnbiased_U_int
float svdTrkPosUnbiased_U_int
U internal unbiased track position.
Definition: OverlapResidualsModule.h:204
Belle2::OverlapResidualsModule::h_DeltaResVPhi_Lyr6
TH2F * h_DeltaResVPhi_Lyr6
2D histogram of DeltaRes_v vs phi of VXD overlaps for layer 6
Definition: OverlapResidualsModule.h:125
Belle2::OverlapResidualsModule::svdTrkQoP_U_int
float svdTrkQoP_U_int
U internal track q/p.
Definition: OverlapResidualsModule.h:201
Belle2::OverlapResidualsModule::h_V_DeltaRes_SVD
TH1F * h_V_DeltaRes_SVD
Histogram of SVD v-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:79
Belle2::OverlapResidualsModule::svdRes_U_int
float svdRes_U_int
U internal residual computed by genfit.
Definition: OverlapResidualsModule.h:185
Belle2::OverlapResidualsModule::svdTrkPosErr_U_int
float svdTrkPosErr_U_int
U internal track position error.
Definition: OverlapResidualsModule.h:199
Belle2::OverlapResidualsModule::svdStrip6Samples_V_ext
std::vector< float > svdStrip6Samples_V_ext
V external 6 samples of the strips of the cluster.
Definition: OverlapResidualsModule.h:308
Belle2::OverlapResidualsModule::svdStripPosition_U_ext
std::vector< float > svdStripPosition_U_ext
U external position of the strips of the cluster.
Definition: OverlapResidualsModule.h:248
Belle2::OverlapResidualsModule::h_DeltaResVz_Lyr1
TH2F * h_DeltaResVz_Lyr1
2D histogram of DeltaRes_v vs z of VXD overlaps for layer 1
Definition: OverlapResidualsModule.h:123
Belle2::OverlapResidualsModule::svdTrkQoPUnbiased_V_ext
float svdTrkQoPUnbiased_V_ext
V external unbiased track q/p.
Definition: OverlapResidualsModule.h:293
Belle2::OverlapResidualsModule::h_DeltaResUz_Lyr3
TH2F * h_DeltaResUz_Lyr3
2D histogram of DeltaRes_u vs z of VXD overlaps for layer 3
Definition: OverlapResidualsModule.h:107
Belle2::OverlapResidualsModule::svdTrkPos_U_ext
float svdTrkPos_U_ext
U external track position.
Definition: OverlapResidualsModule.h:219
Belle2::OverlapResidualsModule::svdTruePos_V_int
float svdTruePos_V_int
V internal true position.
Definition: OverlapResidualsModule.h:258
Belle2::OverlapResidualsModule::svdClPhi_V_ext
float svdClPhi_V_ext
V external cluster global phi.
Definition: OverlapResidualsModule.h:281
Belle2::OverlapResidualsModule::extResV_PXD
float extResV_PXD
tree t_PXD branch extResV_PXD/F
Definition: OverlapResidualsModule.h:169
Belle2::OverlapResidualsModule::svdTrkPosUnbiased_V_ext
float svdTrkPosUnbiased_V_ext
V external unbiased track position.
Definition: OverlapResidualsModule.h:291
Belle2::OverlapResidualsModule::h_U_DeltaRes_SVD_Lyr4
TH1F * h_U_DeltaRes_SVD_Lyr4
Histogram of SVD layer-4 u-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:73
Belle2::OverlapResidualsModule::svdStripCharge_U_ext
std::vector< float > svdStripCharge_U_ext
U external charge of the strips of the cluster.
Definition: OverlapResidualsModule.h:245
Belle2::OverlapResidualsModule::svdLayer_U_ext
unsigned int svdLayer_U_ext
U external layer.
Definition: OverlapResidualsModule.h:237
Belle2::OverlapResidualsModule::svdTrkPrimeUnbiased_U_int
float svdTrkPrimeUnbiased_U_int
U internal unbiased tan of incident angle projected on u,w.
Definition: OverlapResidualsModule.h:207
Belle2::OverlapResidualsModule::extResU_PXD
float extResU_PXD
tree t_PXD branch extResU_PXD/F
Definition: OverlapResidualsModule.h:168
Belle2::OverlapResidualsModule::h_DeltaResUPhi_Lyr2
TH2F * h_DeltaResUPhi_Lyr2
2D histogram of DeltaRes_u vs phi of VXD overlaps for layer 2
Definition: OverlapResidualsModule.h:97
Belle2::OverlapResidualsModule::h_DeltaResUPhi_Lyr5
TH2F * h_DeltaResUPhi_Lyr5
2D histogram of DeltaRes_u vs phi of VXD overlaps for layer 5
Definition: OverlapResidualsModule.h:91
Belle2::OverlapResidualsModule::h_DeltaResUPhi_Lyr3
TH2F * h_DeltaResUPhi_Lyr3
2D histogram of DeltaRes_u vs phi of VXD overlaps for layer 3
Definition: OverlapResidualsModule.h:95
Belle2::OverlapResidualsModule::svdTrkPos_V_int
float svdTrkPos_V_int
V internal track position.
Definition: OverlapResidualsModule.h:262
Belle2::OverlapResidualsModule::h_U_DeltaRes_PXD_Lyr1
TH1F * h_U_DeltaRes_PXD_Lyr1
Histogram of PXD layer-1 u-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:59
Belle2::OverlapResidualsModule::svdTrkTraversedLength_V_int
float svdTrkTraversedLength_V_int
V internal traversed length of the track in the sensor.
Definition: OverlapResidualsModule.h:261
Belle2::OverlapResidualsModule::h_DeltaResUz_Lyr2
TH2F * h_DeltaResUz_Lyr2
2D histogram of DeltaRes_u vs z of VXD overlaps for layer 2
Definition: OverlapResidualsModule.h:109
Belle2::OverlapResidualsModule::svdTrkPosOS_V_ext
float svdTrkPosOS_V_ext
V external track position on the other side.
Definition: OverlapResidualsModule.h:285
Belle2::OverlapResidualsModule::svdClTime_V_ext
float svdClTime_V_ext
V external cluster time.
Definition: OverlapResidualsModule.h:275
Belle2::OverlapResidualsModule::svdRes_V_ext
float svdRes_V_ext
V external residual computed by genfit.
Definition: OverlapResidualsModule.h:276
Belle2::OverlapResidualsModule::svdTrkPosOS_U_ext
float svdTrkPosOS_U_ext
U external track position on the other side.
Definition: OverlapResidualsModule.h:220
Belle2::OverlapResidualsModule::svdTrkPosErrUnbiased_V_int
float svdTrkPosErrUnbiased_V_int
V internal unbiased track position error.
Definition: OverlapResidualsModule.h:270
Belle2::OverlapResidualsModule::svdTrkPrimeOS_U_int
float svdTrkPrimeOS_U_int
U internal tan of incident angle projected on v/u,w (other side)
Definition: OverlapResidualsModule.h:203
Belle2::OverlapResidualsModule::svdStripTime_V_ext
std::vector< float > svdStripTime_V_ext
V external time of the strips of the cluster.
Definition: OverlapResidualsModule.h:309
Belle2::OverlapResidualsModule::h_DeltaResUPhi_Lyr1
TH2F * h_DeltaResUPhi_Lyr1
2D Histogram of DeltaRes_u vs phi of VXD overlaps for layer 1
Definition: OverlapResidualsModule.h:99
Belle2::OverlapResidualsModule::svdLadder_U_ext
unsigned int svdLadder_U_ext
U external ladder.
Definition: OverlapResidualsModule.h:238
Belle2::OverlapResidualsModule::h_DeltaResVz_Lyr4
TH2F * h_DeltaResVz_Lyr4
2D histogram of DeltaRes_v vs z of VXD overlaps for layer 4
Definition: OverlapResidualsModule.h:117
Belle2::OverlapResidualsModule::svdTrkPXDHits
int svdTrkPXDHits
number of PXD hits on the track
Definition: OverlapResidualsModule.h:230
Belle2::OverlapResidualsModule::h_U_DeltaRes_PXD
TH1F * h_U_DeltaRes_PXD
Histogram of PXD u-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:57
Belle2::OverlapResidualsModule::extLayer_PXD
unsigned int extLayer_PXD
tree t_PXD branch extLayer_PXD/i
Definition: OverlapResidualsModule.h:177
Belle2::OverlapResidualsModule::initialize
void initialize() override
Register input and output data.
Definition: OverlapResidualsModule.cc:56
Belle2::OverlapResidualsModule::h_DeltaResVPhi_Lyr5
TH2F * h_DeltaResVPhi_Lyr5
2D histogram of DeltaRes_v vs phi of VXD overlaps for layer 5
Definition: OverlapResidualsModule.h:127
Belle2::OverlapResidualsModule::svdTrkQoP_V_ext
float svdTrkQoP_V_ext
V external track q/p.
Definition: OverlapResidualsModule.h:288
Belle2::OverlapResidualsModule::svdClIntStrPos_V_ext
float svdClIntStrPos_V_ext
V external cluster interstrip position.
Definition: OverlapResidualsModule.h:277
Belle2::OverlapResidualsModule::svdStrip6Samples_U_ext
std::vector< float > svdStrip6Samples_U_ext
U external 6 samples of the strips of the cluster.
Definition: OverlapResidualsModule.h:246
Belle2::OverlapResidualsModule::svdClSNR_V_int
float svdClSNR_V_int
V internal cluster SNR.
Definition: OverlapResidualsModule.h:252
Belle2::OverlapResidualsModule::svdTrkPrimeUnbiased_V_ext
float svdTrkPrimeUnbiased_V_ext
V external unbiased tan of incident angle projected on u,w.
Definition: OverlapResidualsModule.h:294
Belle2::OverlapResidualsModule::svdLadder_V_ext
unsigned int svdLadder_V_ext
V external ladder.
Definition: OverlapResidualsModule.h:300
Belle2::OverlapResidualsModule::h_Lyr1
TH2F * h_Lyr1[9][3]
Sensor hit-map for layer 1 from reconstructed u and v coordinates.
Definition: OverlapResidualsModule.h:153
Belle2::OverlapResidualsModule::svdClSNR_U_ext
float svdClSNR_U_ext
U external cluster SNR.
Definition: OverlapResidualsModule.h:209
Belle2::OverlapResidualsModule::svdLadder_U_int
unsigned int svdLadder_U_int
U internal ladder.
Definition: OverlapResidualsModule.h:234
Belle2::OverlapResidualsModule::svdClSNR_V_ext
float svdClSNR_V_ext
V external cluster SNR.
Definition: OverlapResidualsModule.h:274
Belle2::OverlapResidualsModule::svdLayer_V_ext
unsigned int svdLayer_V_ext
V external layer.
Definition: OverlapResidualsModule.h:299
Belle2::OverlapResidualsModule::svdTrkQoPUnbiased_U_int
float svdTrkQoPUnbiased_U_int
U internal unbiased track q/p.
Definition: OverlapResidualsModule.h:206
Belle2::OverlapResidualsModule::svdTrkPrimeUnbiased_U_ext
float svdTrkPrimeUnbiased_U_ext
U external unbiased tan of incident angle projected on u,w.
Definition: OverlapResidualsModule.h:229
Belle2::OverlapResidualsModule::event
void event() override
Compute the difference of coordinate residuals between two hits in overlapping sensors of a same VXD ...
Definition: OverlapResidualsModule.cc:468
Belle2::OverlapResidualsModule::h_DeltaResVPhi_Lyr3
TH2F * h_DeltaResVPhi_Lyr3
2D histogram of DeltaRes_v vs phi of VXD overlaps for layer 3
Definition: OverlapResidualsModule.h:131
Belle2::OverlapResidualsModule::h_DeltaResVPhi_Lyr4
TH2F * h_DeltaResVPhi_Lyr4
2D histogram of DeltaRes_v vs phi of VXD overlaps for layer 4
Definition: OverlapResidualsModule.h:129
Belle2::OverlapResidualsModule::h_V_DeltaRes_PXD_Lyr1
TH1F * h_V_DeltaRes_PXD_Lyr1
Histogram of PXD layer-1 v-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:65
Belle2::OverlapResidualsModule::extU_PXD
float extU_PXD
tree t_PXD branch extU_PXD/F
Definition: OverlapResidualsModule.h:170
Belle2::OverlapResidualsModule::svdTrkTraversedLength_U_ext
float svdTrkTraversedLength_U_ext
U external traversed length of the track in the sensor.
Definition: OverlapResidualsModule.h:218
Belle2::OverlapResidualsModule::extPhi_PXD
float extPhi_PXD
tree t_PXD branch extPhi/PXD/F
Definition: OverlapResidualsModule.h:172
Belle2::OverlapResidualsModule::svdClZ_U_int
float svdClZ_U_int
U internal cluster global Z.
Definition: OverlapResidualsModule.h:191
Belle2::OverlapResidualsModule::svdTruePos_U_int
float svdTruePos_U_int
U internal true position.
Definition: OverlapResidualsModule.h:189
Belle2::OverlapResidualsModule::h_DeltaResVz_Lyr3
TH2F * h_DeltaResVz_Lyr3
2D histogram of DeltaRes_v vs z of VXD overlaps for layer 3
Definition: OverlapResidualsModule.h:119
Belle2::OverlapResidualsModule::svdRes_U_ext
float svdRes_U_ext
U external residual computed by genfit.
Definition: OverlapResidualsModule.h:211
Belle2::OverlapResidualsModule::extLadder_PXD
unsigned int extLadder_PXD
tree t_PXD branch extLadder_PXD/i
Definition: OverlapResidualsModule.h:178
Belle2::OverlapResidualsModule::svdSensor_V_int
unsigned int svdSensor_V_int
V internal sensor.
Definition: OverlapResidualsModule.h:297
Belle2::OverlapResidualsModule::svdTruePos_V_ext
float svdTruePos_V_ext
V external true position.
Definition: OverlapResidualsModule.h:280
Belle2::OverlapResidualsModule::svdClPos_U_int
float svdClPos_U_int
U internal cluster position.
Definition: OverlapResidualsModule.h:187
Belle2::OverlapResidualsModule::h_U_DeltaRes_SVD_Lyr6
TH1F * h_U_DeltaRes_SVD_Lyr6
Histogram of SVD layer-6 u-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:77
Belle2::OverlapResidualsModule::t_SVD_V
TTree * t_SVD_V
Tree containing global information on SVD v-coordinate overlaps.
Definition: OverlapResidualsModule.h:159
Belle2::OverlapResidualsModule::h_V_DeltaRes_PXD
TH1F * h_V_DeltaRes_PXD
Histogram of PXD v-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:63
Belle2::OverlapResidualsModule::intLayer_PXD
unsigned int intLayer_PXD
tree t_PXD branch intLayer_PXD/i
Definition: OverlapResidualsModule.h:174
Belle2::OverlapResidualsModule::svdClPosErr_U_ext
float svdClPosErr_U_ext
U external cluster position error.
Definition: OverlapResidualsModule.h:214
Belle2::OverlapResidualsModule::t_SVD_U
TTree * t_SVD_U
Tree containing global information on SVD u-coordinate overlaps.
Definition: OverlapResidualsModule.h:157
Belle2::OverlapResidualsModule::svdStripCharge_V_int
std::vector< float > svdStripCharge_V_int
V internal charge of the strips of the cluster.
Definition: OverlapResidualsModule.h:303
Belle2::OverlapResidualsModule::svdClCharge_U_int
float svdClCharge_U_int
U internal cluster charge.
Definition: OverlapResidualsModule.h:182
Belle2::OverlapResidualsModule::svdClCharge_U_ext
float svdClCharge_U_ext
U external cluster charge.
Definition: OverlapResidualsModule.h:208
Belle2::OverlapResidualsModule::svdStrip6Samples_V_int
std::vector< float > svdStrip6Samples_V_int
V internal 6 samples of the strips of the cluster.
Definition: OverlapResidualsModule.h:304
Belle2::OverlapResidualsModule::h_DeltaResUPhi_Lyr6
TH2F * h_DeltaResUPhi_Lyr6
2D histogram of DeltaRes_u vs phi of VXD overlaps for layer 6
Definition: OverlapResidualsModule.h:89
Belle2::OverlapResidualsModule::svdClTime_U_int
float svdClTime_U_int
U internal cluster time.
Definition: OverlapResidualsModule.h:184
Belle2::OverlapResidualsModule::h_DeltaResUz_Lyr5
TH2F * h_DeltaResUz_Lyr5
2D histogram of DeltaRes_u vs z of VXD overlaps for layer 5
Definition: OverlapResidualsModule.h:103
Belle2::OverlapResidualsModule::svdTrkPosUnbiased_V_int
float svdTrkPosUnbiased_V_int
V internal unbiased track position.
Definition: OverlapResidualsModule.h:269
Belle2::OverlapResidualsModule::svdClIntStrPos_U_int
float svdClIntStrPos_U_int
U internal cluster interstrip position.
Definition: OverlapResidualsModule.h:186
Belle2::OverlapResidualsModule::svdTrkQoPUnbiased_V_int
float svdTrkQoPUnbiased_V_int
V internal unbiased track q/p.
Definition: OverlapResidualsModule.h:271
Belle2::OverlapResidualsModule::svdClTime_U_ext
float svdClTime_U_ext
U external cluster time.
Definition: OverlapResidualsModule.h:210
Belle2::OverlapResidualsModule::svdStripCharge_U_int
std::vector< float > svdStripCharge_U_int
U internal charge of the strips of the cluster.
Definition: OverlapResidualsModule.h:241
Belle2::OverlapResidualsModule::svdLadder_V_int
unsigned int svdLadder_V_int
V internal ladder.
Definition: OverlapResidualsModule.h:296
Belle2::OverlapResidualsModule::svdSize_V_int
unsigned int svdSize_V_int
V internal size.
Definition: OverlapResidualsModule.h:298
Belle2::OverlapResidualsModule::h_DeltaResVz_Lyr5
TH2F * h_DeltaResVz_Lyr5
2D histogram of DeltaRes_v vs z of VXD overlaps for layer 5
Definition: OverlapResidualsModule.h:115
Belle2::OverlapResidualsModule::svdLayer_V_int
unsigned int svdLayer_V_int
V internal layer.
Definition: OverlapResidualsModule.h:295
Belle2::OverlapResidualsModule::svdTrkPosErr_V_ext
float svdTrkPosErr_V_ext
V external track position error.
Definition: OverlapResidualsModule.h:286
Belle2::OverlapResidualsModule::svdTrkPrime_U_ext
float svdTrkPrime_U_ext
U external tan of incident angle projected on u,w.
Definition: OverlapResidualsModule.h:224
Belle2::OverlapResidualsModule::svdDeltaRes_V
float svdDeltaRes_V
V difference between external and internal residual.
Definition: OverlapResidualsModule.h:250
Belle2::OverlapResidualsModule::h_U_Cl1Cl2_DeltaRes
TH1F * h_U_Cl1Cl2_DeltaRes[5]
Histogram of SVD differences of u-coordinate residuals grouped by clusters sizes.
Definition: OverlapResidualsModule.h:139
Belle2::OverlapResidualsModule::svdTrkPosErrUnbiased_U_ext
float svdTrkPosErrUnbiased_U_ext
U external unbiased track position error.
Definition: OverlapResidualsModule.h:227
Belle2::OverlapResidualsModule::h_DeltaResVz_Lyr6
TH2F * h_DeltaResVz_Lyr6
2D histogram of DeltaRes_v vs z of VXD overlaps for layer 6
Definition: OverlapResidualsModule.h:113
Belle2::OverlapResidualsModule::extV_PXD
float extV_PXD
tree t_PXD branch extV_PXD/F
Definition: OverlapResidualsModule.h:171
Belle2::OverlapResidualsModule::svdTrkPosOS_U_int
float svdTrkPosOS_U_int
U internal track position on the other side.
Definition: OverlapResidualsModule.h:198
Belle2::OverlapResidualsModule::svdTrkQoPUnbiased_U_ext
float svdTrkQoPUnbiased_U_ext
U external unbiased track q/p.
Definition: OverlapResidualsModule.h:228
Belle2::OverlapResidualsModule::h_Lyr4
TH2F * h_Lyr4[11][4]
Sensor hit-map for layer 4 from reconstructed u and v coordinates.
Definition: OverlapResidualsModule.h:147
Belle2::OverlapResidualsModule::h_V_DeltaRes_SVD_Lyr6
TH1F * h_V_DeltaRes_SVD_Lyr6
Histogram of SVD layer-6 v-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:87
Belle2::OverlapResidualsModule::svdSensor_V_ext
unsigned int svdSensor_V_ext
V external sensor.
Definition: OverlapResidualsModule.h:301
Belle2::OverlapResidualsModule::svdTrkCDCHits
int svdTrkCDCHits
number of PXD hits on the track
Definition: OverlapResidualsModule.h:232
Belle2::OverlapResidualsModule::svdClPosErr_V_int
float svdClPosErr_V_int
V internal cluster position error.
Definition: OverlapResidualsModule.h:257
Belle2::OverlapResidualsModule::h_V_DeltaRes_SVD_Lyr3
TH1F * h_V_DeltaRes_SVD_Lyr3
Histogram of SVD layer-3 v-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:81
Belle2::OverlapResidualsModule::h_V_Cl1Cl2_DeltaRes
TH1F * h_V_Cl1Cl2_DeltaRes[5]
Histogram of SVD differences of v-coordinate residuals grouped by clusters sizes.
Definition: OverlapResidualsModule.h:141
Belle2::OverlapResidualsModule::svdTrkPosErr_U_ext
float svdTrkPosErr_U_ext
U external track position error.
Definition: OverlapResidualsModule.h:221
Belle2::OverlapResidualsModule::svdTrkPrime_U_int
float svdTrkPrime_U_int
U internal tan of incident angle projected on u,w.
Definition: OverlapResidualsModule.h:202
Belle2::OverlapResidualsModule::svdRes_V_int
float svdRes_V_int
V internal residual computed by genfit.
Definition: OverlapResidualsModule.h:254
Belle2::OverlapResidualsModule::svdTrkPrimeOS_V_ext
float svdTrkPrimeOS_V_ext
V external tan of incident angle projected on v/u,w (other side)
Definition: OverlapResidualsModule.h:290
Belle2::OverlapResidualsModule::h_V_DeltaRes_PXD_Lyr2
TH1F * h_V_DeltaRes_PXD_Lyr2
Histogram of PXD layer-2 v-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:67
Belle2::OverlapResidualsModule::h_U_DeltaRes_SVD_Lyr5
TH1F * h_U_DeltaRes_SVD_Lyr5
Histogram of SVD layer-5 u-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:75
Belle2::OverlapResidualsModule::svdTrkPosErr_V_int
float svdTrkPosErr_V_int
V internal track position error.
Definition: OverlapResidualsModule.h:264
Belle2::OverlapResidualsModule::h_V_DeltaRes_SVD_Lyr5
TH1F * h_V_DeltaRes_SVD_Lyr5
Histogram of SVD layer-5 v-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:85
Belle2::OverlapResidualsModule::svdClPosErr_V_ext
float svdClPosErr_V_ext
V external cluster position error.
Definition: OverlapResidualsModule.h:279
Belle2::OverlapResidualsModule::deltaResU_PXD
float deltaResU_PXD
tree t_PXD branch deltaResU/F
Definition: OverlapResidualsModule.h:161
Belle2::OverlapResidualsModule::h_Lyr2
TH2F * h_Lyr2[13][3]
Sensor hit-map for layer 2 from reconstructed u and v coordinates.
Definition: OverlapResidualsModule.h:151
Belle2::OverlapResidualsModule::svdDeltaRes_U
float svdDeltaRes_U
U difference between external and internal residual.
Definition: OverlapResidualsModule.h:181
Belle2::OverlapResidualsModule::intU_PXD
float intU_PXD
tree t_PXD branch intU_PXD/F
Definition: OverlapResidualsModule.h:164
Belle2::OverlapResidualsModule::extSensor_PXD
unsigned int extSensor_PXD
tree t_PXD branch extSensor_PXD/i
Definition: OverlapResidualsModule.h:179
Belle2::OverlapResidualsModule::svdStripPosition_V_ext
std::vector< float > svdStripPosition_V_ext
V external position of the strips of the cluster.
Definition: OverlapResidualsModule.h:310
Belle2::OverlapResidualsModule::svdClIntStrPos_V_int
float svdClIntStrPos_V_int
V internal cluster interstrip position.
Definition: OverlapResidualsModule.h:255
Belle2::OverlapResidualsModule::svdTrkPosErrUnbiased_U_int
float svdTrkPosErrUnbiased_U_int
U internal unbiased track position error.
Definition: OverlapResidualsModule.h:205
Belle2::OverlapResidualsModule::h_V_DeltaRes_SVD_Lyr4
TH1F * h_V_DeltaRes_SVD_Lyr4
Histogram of SVD layer-4 v-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:83
Belle2::OverlapResidualsModule::svdTrkPosErrUnbiased_V_ext
float svdTrkPosErrUnbiased_V_ext
V external unbiased track position error.
Definition: OverlapResidualsModule.h:292
Belle2::OverlapResidualsModule::svdStripCharge_V_ext
std::vector< float > svdStripCharge_V_ext
V external charge of the strips of the cluster.
Definition: OverlapResidualsModule.h:307
Belle2::OverlapResidualsModule::svdStripTime_U_ext
std::vector< float > svdStripTime_U_ext
U external time of the strips of the cluster.
Definition: OverlapResidualsModule.h:247
Belle2::OverlapResidualsModule::h_DeltaResUPhi_Lyr4
TH2F * h_DeltaResUPhi_Lyr4
2D histogram of DeltaRes_u vs phi of VXD overlaps for layer 4
Definition: OverlapResidualsModule.h:93
Belle2::OverlapResidualsModule::intZ_PXD
float intZ_PXD
tree t_PXD branch intZ_PXD/F
Definition: OverlapResidualsModule.h:167
Belle2::OverlapResidualsModule::svdClPos_V_ext
float svdClPos_V_ext
V external cluster position.
Definition: OverlapResidualsModule.h:278
Belle2::OverlapResidualsModule::svdTrkPosErrOS_V_int
float svdTrkPosErrOS_V_int
V internal track position error on the other side.
Definition: OverlapResidualsModule.h:265
Belle2::OverlapResidualsModule::svdStripPosition_U_int
std::vector< float > svdStripPosition_U_int
U internal position of the strips of the cluster.
Definition: OverlapResidualsModule.h:244
Belle2::OverlapResidualsModule::m_ExpertLevel
bool m_ExpertLevel
Expert level switch.
Definition: OverlapResidualsModule.h:45
Belle2::OverlapResidualsModule::svdSensor_U_int
unsigned int svdSensor_U_int
U internal sensor.
Definition: OverlapResidualsModule.h:235
Belle2::OverlapResidualsModule::svdStrip6Samples_U_int
std::vector< float > svdStrip6Samples_U_int
U internal 6 samples of the strips of the cluster.
Definition: OverlapResidualsModule.h:242
Belle2::OverlapResidualsModule::svdClPosErr_U_int
float svdClPosErr_U_int
U internal cluster position error.
Definition: OverlapResidualsModule.h:188
Belle2::OverlapResidualsModule::svdClIntStrPos_U_ext
float svdClIntStrPos_U_ext
U external cluster interstrip position.
Definition: OverlapResidualsModule.h:212
Belle2::OverlapResidualsModule::svdClTime_V_int
float svdClTime_V_int
V internal cluster time.
Definition: OverlapResidualsModule.h:253
Belle2::OverlapResidualsModule::h_DeltaResVPhi_Lyr2
TH2F * h_DeltaResVPhi_Lyr2
2D histogram of DeltaRes_v vs phi of VXD overlaps for layer 2
Definition: OverlapResidualsModule.h:133
Belle2::OverlapResidualsModule::h_U_DeltaRes
TH1F * h_U_DeltaRes
Histogram of VXD ( PXD + SVD ) differences of u-coordinate residuals.
Definition: OverlapResidualsModule.h:53
Belle2::OverlapResidualsModule::svdTrkPrime_V_ext
float svdTrkPrime_V_ext
V external tan of incident angle projected on u,w.
Definition: OverlapResidualsModule.h:289
Belle2::OverlapResidualsModule::svdTrkTraversedLength_U_int
float svdTrkTraversedLength_U_int
U internal traversed length of the track in the sensor.
Definition: OverlapResidualsModule.h:196
Belle2::OverlapResidualsModule::h_DeltaResUz_Lyr4
TH2F * h_DeltaResUz_Lyr4
2D histogram of DeltaRes_u vs z of VXD overlaps for layer 4
Definition: OverlapResidualsModule.h:105
Belle2::OverlapResidualsModule::h_Lyr5
TH2F * h_Lyr5[13][5]
Sensor hit-map for layer 5 from reconstructed u and v coordinates.
Definition: OverlapResidualsModule.h:145
Belle2::OverlapResidualsModule::h_SVDstrips_Mult
TH1F * h_SVDstrips_Mult
Histogram of SVD strips multiplicity.
Definition: OverlapResidualsModule.h:137
Belle2::OverlapResidualsModule::h_DeltaResUz_Lyr6
TH2F * h_DeltaResUz_Lyr6
2D histogram of DeltaRes_u vs z of VXD overlaps for layer 6
Definition: OverlapResidualsModule.h:101
Belle2::OverlapResidualsModule::extZ_PXD
float extZ_PXD
tree t_PXD branch extZ_PXD/F
Definition: OverlapResidualsModule.h:173
Belle2::OverlapResidualsModule::h_Lyr3
TH2F * h_Lyr3[8][3]
Sensor hit-map for layer 3 from reconstructed u and v coordinates.
Definition: OverlapResidualsModule.h:149
Belle2::OverlapResidualsModule::svdClZ_V_int
float svdClZ_V_int
V internal cluster global Z.
Definition: OverlapResidualsModule.h:260
Belle2::OverlapResidualsModule::svdTrkPrime_V_int
float svdTrkPrime_V_int
V internal tan of incident angle projected on u,w.
Definition: OverlapResidualsModule.h:267
Belle2::OverlapResidualsModule::svdTrkTraversedLength_V_ext
float svdTrkTraversedLength_V_ext
V external traversed length of the track in the sensor.
Definition: OverlapResidualsModule.h:283
Belle2::OverlapResidualsModule::svdClCharge_V_ext
float svdClCharge_V_ext
V external cluster charge.
Definition: OverlapResidualsModule.h:273
Belle2::OverlapResidualsModule::svdSize_U_int
unsigned int svdSize_U_int
U internal size.
Definition: OverlapResidualsModule.h:236
Belle2::OverlapResidualsModule::h_DeltaResVPhi_Lyr1
TH2F * h_DeltaResVPhi_Lyr1
2D histogram of DeltaRes_v vs phi of VXD overlaps for layer 1
Definition: OverlapResidualsModule.h:135
Belle2::OverlapResidualsModule::svdStripTime_V_int
std::vector< float > svdStripTime_V_int
V internal time of the strips of the cluster.
Definition: OverlapResidualsModule.h:305
Belle2::OverlapResidualsModule::svdClPhi_U_int
float svdClPhi_U_int
U internal cluster global phi.
Definition: OverlapResidualsModule.h:190
Belle2::OverlapResidualsModule::h_DeltaResUz_Lyr1
TH2F * h_DeltaResUz_Lyr1
2D histogram of DeltaRes_u vs z of VXD overlaps for layer 1
Definition: OverlapResidualsModule.h:111
Belle2::OverlapResidualsModule::svdClSNR_U_int
float svdClSNR_U_int
U internal cluster SNR.
Definition: OverlapResidualsModule.h:183
Belle2::OverlapResidualsModule::svdLayer_U_int
unsigned int svdLayer_U_int
U internal layer.
Definition: OverlapResidualsModule.h:233
Belle2::OverlapResidualsModule::svdClCharge_V_int
float svdClCharge_V_int
V internal cluster charge.
Definition: OverlapResidualsModule.h:251
Belle2::OverlapResidualsModule::svdSensor_U_ext
unsigned int svdSensor_U_ext
U external sensor.
Definition: OverlapResidualsModule.h:239
Belle2::OverlapResidualsModule::intV_PXD
float intV_PXD
tree t_PXD branch intV_PXD/F
Definition: OverlapResidualsModule.h:165
Belle2::OverlapResidualsModule::h_Lyr6
TH2F * h_Lyr6[17][6]
Sensor hit-map for layer 6 from reconstructed u and v coordinates.
Definition: OverlapResidualsModule.h:143
Belle2::OverlapResidualsModule::svdStripPosition_V_int
std::vector< float > svdStripPosition_V_int
V internal position of the strips of the cluster.
Definition: OverlapResidualsModule.h:306
Belle2::OverlapResidualsModule::svdTrkPrimeOS_U_ext
float svdTrkPrimeOS_U_ext
U external tan of incident angle projected on v/u,w (other side)
Definition: OverlapResidualsModule.h:225
Belle2::OverlapResidualsModule::intSensor_PXD
unsigned int intSensor_PXD
tree t_PXD branch intSensor_PXD/i
Definition: OverlapResidualsModule.h:176
Belle2::OverlapResidualsModule::svdStripTime_U_int
std::vector< float > svdStripTime_U_int
U internal time of the strips of the cluster.
Definition: OverlapResidualsModule.h:243
Belle2::OverlapResidualsModule::svdTrkPrimeUnbiased_V_int
float svdTrkPrimeUnbiased_V_int
V internal unbiased tan of incident angle projected on u,w.
Definition: OverlapResidualsModule.h:272
Belle2::OverlapResidualsModule::h_DeltaResVz_Lyr2
TH2F * h_DeltaResVz_Lyr2
2D histogram of DeltaRes_v vs z of VXD overlaps for layer 2
Definition: OverlapResidualsModule.h:121
Belle2::OverlapResidualsModule::intLadder_PXD
unsigned int intLadder_PXD
tree t_PXD branch intLadder_PXD/i
Definition: OverlapResidualsModule.h:175
Belle2::OverlapResidualsModule::svdClPos_V_int
float svdClPos_V_int
V internal cluster position.
Definition: OverlapResidualsModule.h:256
Belle2::OverlapResidualsModule::m_recoTracksStoreArrayName
std::string m_recoTracksStoreArrayName
StoreArray name of the input and output RecoTracks.
Definition: OverlapResidualsModule.h:47
Belle2::OverlapResidualsModule::svdTrkPos_V_ext
float svdTrkPos_V_ext
V external track position.
Definition: OverlapResidualsModule.h:284
Belle2::OverlapResidualsModule::svdClZ_U_ext
float svdClZ_U_ext
U external cluster global Z.
Definition: OverlapResidualsModule.h:217
Belle2::OverlapResidualsModule::svdTrkSVDHits
int svdTrkSVDHits
number of PXD hits on the track
Definition: OverlapResidualsModule.h:231
Belle2::OverlapResidualsModule::intPhi_PXD
float intPhi_PXD
tree t_PXD branch intPhi_PXD/F
Definition: OverlapResidualsModule.h:166
Belle2::OverlapResidualsModule::svdTrkPrimeOS_V_int
float svdTrkPrimeOS_V_int
V internal tan of incident angle projected on v/u,w (other side)
Definition: OverlapResidualsModule.h:268
Belle2::OverlapResidualsModule::svdTrkPosErrOS_U_ext
float svdTrkPosErrOS_U_ext
U external track position error on the other side.
Definition: OverlapResidualsModule.h:222
Belle2::OverlapResidualsModule::t_PXD
TTree * t_PXD
Tree containing global information on PXD overlaps.
Definition: OverlapResidualsModule.h:155
Belle2::OverlapResidualsModule::svdTrkPosErrOS_U_int
float svdTrkPosErrOS_U_int
U internal track position error on the other side.
Definition: OverlapResidualsModule.h:200
Belle2::OverlapResidualsModule::svdTrkQoP_U_ext
float svdTrkQoP_U_ext
U external track q/p.
Definition: OverlapResidualsModule.h:223
Belle2::OverlapResidualsModule::svdClZ_V_ext
float svdClZ_V_ext
V external cluster global Z.
Definition: OverlapResidualsModule.h:282
Belle2::OverlapResidualsModule::svdClPhi_V_int
float svdClPhi_V_int
V internal cluster global phi.
Definition: OverlapResidualsModule.h:259
Belle2::OverlapResidualsModule::h_V_DeltaRes
TH1F * h_V_DeltaRes
Histogram of VXD ( PXD + SVD ) differences of v-coordinate residuals.
Definition: OverlapResidualsModule.h:55
Belle2::OverlapResidualsModule::svdTrkPosOS_V_int
float svdTrkPosOS_V_int
V internal track position on the other side.
Definition: OverlapResidualsModule.h:263
Belle2::OverlapResidualsModule::svdSize_V_ext
unsigned int svdSize_V_ext
V external size.
Definition: OverlapResidualsModule.h:302
Belle2::OverlapResidualsModule::intResV_PXD
float intResV_PXD
tree t_PXD branch intResV/PXD/F
Definition: OverlapResidualsModule.h:163
Belle2::OverlapResidualsModule::h_U_DeltaRes_SVD_Lyr3
TH1F * h_U_DeltaRes_SVD_Lyr3
Histogram of SVD layer-3 u-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:71
Belle2::OverlapResidualsModule::h_U_DeltaRes_PXD_Lyr2
TH1F * h_U_DeltaRes_PXD_Lyr2
Histogram of PXD layer-2 u-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:61
Belle2::OverlapResidualsModule::svdClPhi_U_ext
float svdClPhi_U_ext
U external cluster global phi.
Definition: OverlapResidualsModule.h:216
Belle2::OverlapResidualsModule::h_U_DeltaRes_SVD
TH1F * h_U_DeltaRes_SVD
Histogram of SVD u-coordinate differences of residuals.
Definition: OverlapResidualsModule.h:69
Belle2::OverlapResidualsModule::svdClPos_U_ext
float svdClPos_U_ext
U external cluster position.
Definition: OverlapResidualsModule.h:213
Belle2::OverlapResidualsModule::svdTruePos_U_ext
float svdTruePos_U_ext
U external true position.
Definition: OverlapResidualsModule.h:215
Belle2::OverlapResidualsModule::svdTrkPosErrOS_V_ext
float svdTrkPosErrOS_V_ext
V external track position error on the other side.
Definition: OverlapResidualsModule.h:287
Belle2::OverlapResidualsModule::defineHisto
void defineHisto() override
Histogram definitions such as TH1(), TH2(), TNtuple(), TTree()....
Definition: OverlapResidualsModule.cc:64
Belle2::OverlapResidualsModule::svdTrkPosUnbiased_U_ext
float svdTrkPosUnbiased_U_ext
U external unbiased track position.
Definition: OverlapResidualsModule.h:226
Belle2::OverlapResidualsModule::svdTrkPos_U_int
float svdTrkPos_U_int
U internal track position.
Definition: OverlapResidualsModule.h:197
Belle2::OverlapResidualsModule::svdSize_U_ext
unsigned int svdSize_U_ext
U external size.
Definition: OverlapResidualsModule.h:240
Belle2::OverlapResidualsModule::intResU_PXD
float intResU_PXD
tree t_PXD branch intResU_PXD/F
Definition: OverlapResidualsModule.h:162
Belle2::PXDCluster
The PXD Cluster class This class stores all information about reconstructed PXD clusters The position...
Definition: PXDCluster.h:30
Belle2::PXDCluster::getV
float getV() const
Get v coordinate of hit position.
Definition: PXDCluster.h:136
Belle2::PXDCluster::getSensorID
VxdID getSensorID() const
Get the sensor ID.
Definition: PXDCluster.h:126
Belle2::PXDCluster::getU
float getU() const
Get u coordinate of hit position.
Definition: PXDCluster.h:131
Belle2::RecoHitInformation
This class stores additional information to every CDC/SVD/PXD hit stored in a RecoTrack.
Definition: RecoHitInformation.h:38
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::SVDCluster
The SVD Cluster class This class stores all information about reconstructed SVD clusters.
Definition: SVDCluster.h:29
Belle2::SVDCluster::getClsTime
float getClsTime() const
Get average of waveform maximum times of cluster strip signals.
Definition: SVDCluster.h:134
Belle2::SVDCluster::getSNR
float getSNR() const
Get cluster SNR.
Definition: SVDCluster.h:159
Belle2::SVDCluster::getSize
unsigned short getSize() const
Get cluster size.
Definition: SVDCluster.h:154
Belle2::SVDCluster::getCharge
float getCharge() const
Get collected charge.
Definition: SVDCluster.h:144
Belle2::SVDCluster::getSensorID
VxdID getSensorID() const
Get the sensor ID.
Definition: SVDCluster.h:102
Belle2::SVDCluster::isUCluster
bool isUCluster() const
Get the direction of strips.
Definition: SVDCluster.h:110
Belle2::SVDCluster::getPosition
float getPosition(double v=0) const
Get the coordinate of reconstructed hit.
Definition: SVDCluster.h:117
Belle2::SVDCluster::getPositionSigma
float getPositionSigma() const
Get the error of the reconstructed hit coordinate.
Definition: SVDCluster.h:129
Belle2::SVDShaperDigit
The SVD ShaperDigit class.
Definition: SVDShaperDigit.h:32
Belle2::SVDShaperDigit::getSamples
APVFloatSamples getSamples() const
Get array of samples.
Definition: SVDShaperDigit.h:119
Belle2::StoreAccessorBase::isOptional
bool isOptional(const std::string &name="")
Tell the DataStore about an optional input.
Definition: StoreAccessorBase.h:93
Belle2::StoreArray
Accessor to arrays stored in the data store.
Definition: StoreArray.h:113
Belle2::TrackFitResult
Values of the result of a track fit with a given particle hypothesis.
Definition: TrackFitResult.h:49
Belle2::TrackFitResult::get4Momentum
ROOT::Math::PxPyPzEVector get4Momentum() const
Getter for the 4Momentum at the closest approach of the track in the r/phi projection.
Definition: TrackFitResult.h:125
Belle2::TrackFitResult::getD0
double getD0() const
Getter for d0.
Definition: TrackFitResult.h:179
Belle2::TrackFitResult::getZ0
double getZ0() const
Getter for z0.
Definition: TrackFitResult.h:201
Belle2::TrackFitResult::getMomentum
ROOT::Math::XYZVector getMomentum() const
Getter for vector of momentum at closest approach of track in r/phi projection.
Definition: TrackFitResult.h:116
Belle2::VXD::GeoCache
Class to faciliate easy access to sensor information of the VXD like coordinate transformations or pi...
Definition: GeoCache.h:39
Belle2::VXD::GeoCache::getSensorInfo
const SensorInfoBase & getSensorInfo(Belle2::VxdID id) const
Return a referecne to the SensorInfo of a given SensorID.
Definition: GeoCache.cc:67
Belle2::VXD::GeoCache::getInstance
static GeoCache & getInstance()
Return a reference to the singleton instance.
Definition: GeoCache.cc:214
Belle2::VXD::SensorInfoBase
Base class to provide Sensor Information for PXD and SVD.
Definition: SensorInfoBase.h:29
Belle2::VXD::SensorInfoBase::getVCellPosition
double getVCellPosition(int vID) const
Return the position of a specific strip/pixel in v direction.
Definition: SensorInfoBase.h:180
Belle2::VXD::SensorInfoBase::getUCellPosition
double getUCellPosition(int uID, int vID=-1) const
Return the position of a specific strip/pixel in u direction.
Definition: SensorInfoBase.h:168
Belle2::VXD::SensorInfoBase::pointToGlobal
ROOT::Math::XYZVector pointToGlobal(const ROOT::Math::XYZVector &local, bool reco=false) const
Convert a point from local to global coordinates.
Definition: SensorInfoBase.h:369
Belle2::VXD::SensorInfoBase::getThickness
double getThickness() const
Return the thickness of the sensor.
Definition: SensorInfoBase.h:110
Belle2::VxdID
Class to uniquely identify a any structure of the PXD and SVD.
Definition: VxdID.h:33
Belle2::VxdID::getSensorNumber
baseType getSensorNumber() const
Get the sensor id.
Definition: VxdID.h:100
Belle2::VxdID::getLadderNumber
baseType getLadderNumber() const
Get the ladder id.
Definition: VxdID.h:98
Belle2::VxdID::getLayerNumber
baseType getLayerNumber() const
Get the layer id.
Definition: VxdID.h:96
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::Unit::convertValueToUnit
static double convertValueToUnit(double value, const std::string &unitString)
Converts a floating point value from the standard framework unit to the given unit.
Definition: UnitConst.cc:139
REG_MODULE
#define REG_MODULE(moduleName)
Register the given module (without 'Module' suffix) with the framework.
Definition: Module.h:650
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