Skip to content

Commit ed6e83e

Browse files
mgignacmgignac
authored
Merge pull request #219 from JeffersonLab/mgignac-trksmearing-2021
Track Momentum Smearing: 2021
2 parents b4a5238 + 7b748d5 commit ed6e83e

17 files changed

Lines changed: 1331 additions & 85 deletions

File tree

‎CMakeLists.txt‎

Lines changed: 1 addition & 1 deletion
Original file line numberDiff line numberDiff line change
@@ -65,7 +65,7 @@ include(MacroModule)
6565
# import macro for declaring external dependencies
6666
include(MacroExtDeps)
6767

68-
set(MODULES event utils analysis processing processors)
68+
set(MODULES event analysis utils processing processors)
6969

7070
# build each module in the list
7171
foreach(module ${MODULES})

‎analysis/data/smearing/trackSmearing_2021.json‎

Lines changed: 13 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,13 @@
1+
{
2+
"pSmearing": {
3+
"inclusive": 0.16219805756419403,
4+
"top": 0.2008113904741088,
5+
"bot": 0.14066369819431446
6+
},
7+
"z0Smearing": {
8+
"inclusive": 0.0489615408651479,
9+
"top": 0.049285536755409885,
10+
"bot": 0.04589554481673371
11+
},
12+
"relSmearing": false
13+
}

‎analysis/plotconfigs/fee_smearing/feeSmearing_2021.json‎

Lines changed: 304 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1,304 @@
1+
{
2+
"n_tracks_h" : {
3+
"bins" : 10,
4+
"minX" : 0,
5+
"maxX" : 10,
6+
"xtitle" : "N_{tracks}",
7+
"ytitle" : "Events"
8+
},
9+
"nHits_2d_h" : {
10+
"bins" : 15,
11+
"minX" : 0,
12+
"maxX" : 15,
13+
"xtitle" : "N_{2Dhits}",
14+
"ytitle" : "Tracks"
15+
},
16+
"nShared_h" : {
17+
"bins" : 8,
18+
"minX" : -0.5,
19+
"maxX" : 7.5,
20+
"xtitle" : "N_{shared}^{hits}",
21+
"ytitle" : "Tracks"
22+
},
23+
"d0_h" : {
24+
"bins" : 200,
25+
"minX" : -10,
26+
"maxX" : 10,
27+
"xtitle" : "d_{0} [mm]",
28+
"ytitle" : "Tracks"
29+
},
30+
"Phi_h" : {
31+
"bins" : 100,
32+
"minX" : -0.5,
33+
"maxX" : 0.5,
34+
"xtitle" : "phi_{0}",
35+
"ytitle" : "Tracks"
36+
},
37+
"Omega_h" : {
38+
"bins" : 100,
39+
"minX" : -0.002,
40+
"maxX" : 0.002,
41+
"xtitle" : "#omega",
42+
"ytitle" : "Tracks"
43+
},
44+
"TanLambda_h" : {
45+
"bins" : 200,
46+
"minX" : -0.3,
47+
"maxX" : 0.3,
48+
"xtitle" : "tan(#lambda)",
49+
"ytitle" : "Tracks"
50+
},
51+
"Z0_h" : {
52+
"bins" : 100,
53+
"minX" : -0.5,
54+
"maxX" : 0.5,
55+
"xtitle" : "z_0 [mm]",
56+
"ytitle" : "Tracks"
57+
},
58+
"top_track_z0_h" : {
59+
"bins" : 100,
60+
"minX" : -0.5,
61+
"maxX" : 0.5,
62+
"xtitle" : "z_0 [mm]",
63+
"ytitle" : "Tracks"
64+
},
65+
"bot_track_z0_h" : {
66+
"bins" : 100,
67+
"minX" : -0.5,
68+
"maxX" : 0.5,
69+
"xtitle" : "z_0 [mm]",
70+
"ytitle" : "Tracks"
71+
},
72+
"track_xpos_h" : {
73+
"bins" : 200,
74+
"minX" : -3.0,
75+
"maxX" : 3.0,
76+
"xtitle" : "Ref X [mm]",
77+
"ytitle" : "Tracks"
78+
},
79+
"track_ypos_h" : {
80+
"bins" : 200,
81+
"minX" : -0.5,
82+
"maxX" : 0.5,
83+
"xtitle" : "Ref Y [mm]",
84+
"ytitle" : "Tracks"
85+
},
86+
"time_h" : {
87+
"bins" : 200,
88+
"minX" : -40,
89+
"maxX" : 40,
90+
"xtitle" : "track time [ns]",
91+
"ytitle" : "Tracks"
92+
},
93+
"chi2_h" : {
94+
"bins" : 200,
95+
"minX" : 0,
96+
"maxX" : 100,
97+
"xtitle" : "track #chi^{2}",
98+
"ytitle" : "Tracks"
99+
},
100+
"p_h" : {
101+
"bins" : 200,
102+
"minX" : 0,
103+
"maxX" : 6,
104+
"xtitle" : "p_{e^{-}} [GeV]",
105+
"ytitle" : "Tracks"
106+
},
107+
"p_top_h" : {
108+
"bins" : 200,
109+
"minX" : 0,
110+
"maxX" : 6,
111+
"xtitle" : "p_{e^{-}} [GeV]",
112+
"ytitle" : "Tracks (top)"
113+
},
114+
"p_bot_h" : {
115+
"bins" : 200,
116+
"minX" : 0,
117+
"maxX" : 6,
118+
"xtitle" : "p_{e^{-}} [GeV]",
119+
"ytitle" : "Tracks (bottom)"
120+
},
121+
"p_vs_nHits_hh" : {
122+
"binsX" : 7,
123+
"minX" : 8,
124+
"maxX" : 15,
125+
"binsY" : 200,
126+
"minY" : 0,
127+
"maxY" : 6,
128+
"xtitle" : "Nhits",
129+
"ytitle" : "p_{e^{-}} [GeV]"
130+
},
131+
"p_vs_nHits_top_hh" : {
132+
"binsX" : 7,
133+
"minX" : 8,
134+
"maxX" : 15,
135+
"binsY" : 200,
136+
"minY" : 0,
137+
"maxY" : 6,
138+
"xtitle" : "Nhits",
139+
"ytitle" : "top p_{e^{-}} [GeV]"
140+
},
141+
"p_vs_nHits_bot_hh" : {
142+
"binsX" : 7,
143+
"minX" : 8,
144+
"maxX" : 15,
145+
"binsY" : 200,
146+
"minY" : 0,
147+
"maxY" : 6,
148+
"xtitle" : "Nhits",
149+
"ytitle" : "bot p_{e^{-}} [GeV]"
150+
},
151+
"p_vs_TanLambda_hh" : {
152+
"binsX" : 100,
153+
"minX" : -0.1,
154+
"maxX" : 0.1,
155+
"binsY" : 200,
156+
"minY" : 0,
157+
"maxY" : 6,
158+
"xtitle" : "tan(#lambda)",
159+
"ytitle" : "p_{e^{-}} [GeV]"
160+
},
161+
"p_vs_Phi_hh" : {
162+
"binsX" : 100,
163+
"minX" : -0.2,
164+
"maxX" : 0.2,
165+
"binsY" : 200,
166+
"minY" : 0,
167+
"maxY" : 6,
168+
"xtitle" : "#phi_{0}",
169+
"ytitle" : "p_{e^{-}} [GeV]"
170+
},
171+
"p_vs_TanLambda_Phi_hhh" : {
172+
"binsX" : 100,
173+
"minX" : -0.2,
174+
"maxX" : 0.2,
175+
"binsY" : 100,
176+
"minY" : -0.1,
177+
"maxY" : 0.1,
178+
"binsZ" : 200,
179+
"minZ" : 0,
180+
"maxZ" : 6,
181+
"xtitle" : "#phi_{0}",
182+
"ytitle" : "tan(#lambda)",
183+
"ztitle" : "p_{e^{-}} [GeV]"
184+
},
185+
"p_vs_TanLambda_nHits_hhh" : {
186+
"binsX" : 100,
187+
"minX" : -0.2,
188+
"maxX" : 0.2,
189+
"binsY" : 7,
190+
"minY" : 8,
191+
"maxY" : 15,
192+
"binsZ" : 200,
193+
"minZ" : 0,
194+
"maxZ" : 6,
195+
"xtitle" : "#phi_{0}",
196+
"ytitle" : "tan(#lambda)",
197+
"ztitle" : "p_{e^{-}} [GeV]"
198+
},
199+
"chi2ndf_h" : {
200+
"bins" : 200,
201+
"minX" : 0,
202+
"maxX" : 30,
203+
"xtitle" : "track #chi^{2} / ndf",
204+
"ytitle" : "Tracks"
205+
},
206+
"d0_vs_p_hh" : {
207+
"binsX" : 100,
208+
"minX" : 0,
209+
"maxX" : 6,
210+
"binsY" : 200,
211+
"minY" : -10,
212+
"maxY" : 10,
213+
"xtitle" : "trk p [GeV]",
214+
"ytitle" : "trk d0 [mm]"
215+
},
216+
"d0_vs_phi0_hh" : {
217+
"binsX" : 100,
218+
"minX" : -0.4,
219+
"maxX" : 0.4,
220+
"binsY" : 200,
221+
"minY" : -10,
222+
"maxY" : 10,
223+
"xtitle" : "trk #phi_{0}",
224+
"ytitle" : "trk d0 [mm]"
225+
},
226+
227+
"d0_vs_tanlambda_hh" : {
228+
"binsX" : 200,
229+
"minX" : -0.1,
230+
"maxX" : 0.1,
231+
"binsY" : 200,
232+
"minY" : -10,
233+
"maxY" : 10,
234+
"xtitle" : "trk(#lamdba)",
235+
"ytitle" : "trk d0 [mm]"
236+
},
237+
238+
"z0_vs_p_hh" : {
239+
"binsX" : 100,
240+
"minX" : 0,
241+
"maxX" : 6,
242+
"binsY" : 200,
243+
"minY" : -5,
244+
"maxY" : 5,
245+
"xtitle" : "trk p [GeV]",
246+
"ytitle" : "trk z0 [mm]"
247+
},
248+
249+
"z0_vs_phi0_hh" : {
250+
"binsX" : 100,
251+
"minX" : -0.4,
252+
"maxX" : 0.4,
253+
"binsY" : 200,
254+
"minY" : -5,
255+
"maxY" : 5,
256+
"xtitle" : "trk #phi_{0}",
257+
"ytitle" : "trk z0 [mm]"
258+
},
259+
260+
"z0_vs_tanlambda_hh" : {
261+
"binsX" : 200,
262+
"minX" : -0.1,
263+
"maxX" : 0.1,
264+
"binsY" : 200,
265+
"minY" : -5,
266+
"maxY" : 5,
267+
"xtitle" : "trk(#lamdba)",
268+
"ytitle" : "trk z0 [mm]"
269+
},
270+
"tanlambda_vs_phi0_hh" : {
271+
"binsX" : 100,
272+
"minX" : -0.4,
273+
"maxX" : 0.4,
274+
"binsY" : 100,
275+
"minY" : -0.2,
276+
"maxY" : 0.2,
277+
"xtitle" : "trk #phi_{0}",
278+
"ytitle" : "trk tan(#lambda)"
279+
},
280+
"xpos_at_ecal_h" : {
281+
"bins" : 400,
282+
"minX" : -100,
283+
"maxX" : 100,
284+
"xtitle" : "track x pos [mm]",
285+
"ytitle" : "Tracks"
286+
},
287+
"ypos_at_ecal_h" : {
288+
"bins" : 400,
289+
"minX" : -100,
290+
"maxX" : 100,
291+
"xtitle" : "track y pos [mm]",
292+
"ytitle" : "Tracks"
293+
},
294+
"xypos_at_ecal_hh" : {
295+
"binsX" : 100,
296+
"minX" : -200,
297+
"maxX" : 200,
298+
"binsY" : 100,
299+
"minY" : -200,
300+
"maxY" : 200,
301+
"xtitle" : "trk ecal x [mm]",
302+
"ytitle" : "trk ecal y [mm]"
303+
}
304+
}

‎analysis/selections/trackHit/trackHitAna.json‎

Lines changed: 2 additions & 2 deletions
Original file line numberDiff line numberDiff line change
@@ -10,9 +10,9 @@
1010
"info" : "#chi^{2}/ndf <= 10.0"
1111
},
1212
"p_gt" : {
13-
"cut" : 0.4,
13+
"cut" : 2.0,
1414
"id" : 2,
15-
"info" : "p > 0.4"
15+
"info" : "p > 2.0"
1616
},
1717
"p_lt" : {
1818
"cut" : 6.0,

‎analysis/src/TrackHistos.cxx‎

Lines changed: 6 additions & 3 deletions
Original file line numberDiff line numberDiff line change
@@ -256,11 +256,14 @@ void TrackHistos::Fill1DTrack(Track* track, float weight, const std::string& trk
256256
Fill1DHisto(trkname+"zpos_at_ecal_h",track->getPositionAtEcal().at(2) ,weight);
257257

258258
//Top vs Bot
259-
if(track->getTanLambda() > 0.0)
259+
if(track->getTanLambda() > 0.0){
260260
Fill1DHisto(trkname+"top_track_z0_h", track->getZ0(), weight);
261-
else
261+
Fill1DHisto(trkname+"p_top_h" ,track->getP() , weight);
262+
}
263+
else{
262264
Fill1DHisto(trkname+"bot_track_z0_h", track->getZ0(), weight);
263-
265+
Fill1DHisto(trkname+"p_bot_h" ,track->getP() , weight);
266+
}
264267
//Track param errors
265268
Fill1DHisto(trkname+"d0_err_h" ,track->getD0Err() ,weight);
266269
Fill1DHisto(trkname+"Phi_err_h" ,track->getPhiErr() ,weight);

‎event/include/Track.h‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -90,6 +90,7 @@ class Track : public TObject {
9090
double getOmega () const {return omega_;}
9191
double getTanLambda() const {return tan_lambda_;}
9292
double getZ0 () const {return z0_;}
93+
void setZ0 (const double z0) {z0_ = z0;}
9394

9495
/** Set the covariance matrix **/
9596
void setCov(const std::vector<float>& cov) {cov_ = cov;}

‎plotUtils/.gitignore‎

Lines changed: 1 addition & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -0,0 +1 @@
1+
*root

0 commit comments

Comments
 (0)