Skip to content

Commit 99cb57c

Browse files
abylinkinabylinkin
authored
[PWGUD] Added deuterons and TOF into personal PID task (#8644)
1 parent 733c47d commit 99cb57c

File tree

1 file changed

+138
-0
lines changed

1 file changed

+138
-0
lines changed

‎PWGUD/Tasks/sgPIDAnalyzer.cxx‎

Lines changed: 138 additions & 0 deletions
Original file line numberDiff line numberDiff line change
@@ -73,28 +73,86 @@ struct sgPIDAnalyzer {
7373
histos.add("TPC/pTPC_Pi_Ka", "Positive TPC Pi vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
7474
histos.add("TPC/pTPC_Pi_Pr", "Positive TPC Pi vs Pr", {HistType::kTH2F, {ptBins, nSigmaBins}});
7575
histos.add("TPC/pTPC_Pi_El", "Positive TPC Pi vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
76+
histos.add("TPC/pTPC_Pi_De", "Positive TPC Pi vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
7677
histos.add("TPC/pTPC_Ka_Pi", "Positive TPC Ka vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
7778
histos.add("TPC/pTPC_Ka_Pr", "Positive TPC Ka vs Pr", {HistType::kTH2F, {ptBins, nSigmaBins}});
7879
histos.add("TPC/pTPC_Ka_El", "Positive TPC Ka vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
80+
histos.add("TPC/pTPC_Ka_De", "Positive TPC Ka vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
7981
histos.add("TPC/pTPC_Pr_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
8082
histos.add("TPC/pTPC_Pr_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
8183
histos.add("TPC/pTPC_Pr_El", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
84+
histos.add("TPC/pTPC_Pr_De", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
8285
histos.add("TPC/pTPC_El_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
8386
histos.add("TPC/pTPC_El_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
8487
histos.add("TPC/pTPC_El_Pr", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
88+
histos.add("TPC/pTPC_El_De", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
89+
histos.add("TPC/pTPC_De_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
90+
histos.add("TPC/pTPC_De_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
91+
histos.add("TPC/pTPC_De_Pr", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
92+
histos.add("TPC/pTPC_De_El", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
8593

8694
histos.add("TPC/nTPC_Pi_Ka", "Positive TPC Pi vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
8795
histos.add("TPC/nTPC_Pi_Pr", "Positive TPC Pi vs Pr", {HistType::kTH2F, {ptBins, nSigmaBins}});
8896
histos.add("TPC/nTPC_Pi_El", "Positive TPC Pi vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
97+
histos.add("TPC/nTPC_Pi_De", "Positive TPC Pi vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
8998
histos.add("TPC/nTPC_Ka_Pi", "Positive TPC Ka vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
9099
histos.add("TPC/nTPC_Ka_Pr", "Positive TPC Ka vs Pr", {HistType::kTH2F, {ptBins, nSigmaBins}});
91100
histos.add("TPC/nTPC_Ka_El", "Positive TPC Ka vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
101+
histos.add("TPC/nTPC_Ka_De", "Positive TPC Ka vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
92102
histos.add("TPC/nTPC_Pr_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
93103
histos.add("TPC/nTPC_Pr_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
94104
histos.add("TPC/nTPC_Pr_El", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
105+
histos.add("TPC/nTPC_Pr_De", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
95106
histos.add("TPC/nTPC_El_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
96107
histos.add("TPC/nTPC_El_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
97108
histos.add("TPC/nTPC_El_Pr", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
109+
histos.add("TPC/nTPC_El_De", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
110+
histos.add("TPC/nTPC_De_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
111+
histos.add("TPC/nTPC_De_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
112+
histos.add("TPC/nTPC_De_Pr", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
113+
histos.add("TPC/nTPC_De_El", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
114+
115+
histos.add("TPC/pTOF_Pi_Ka", "Positive TPC Pi vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
116+
histos.add("TPC/pTOF_Pi_Pr", "Positive TPC Pi vs Pr", {HistType::kTH2F, {ptBins, nSigmaBins}});
117+
histos.add("TPC/pTOF_Pi_El", "Positive TPC Pi vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
118+
histos.add("TPC/pTOF_Pi_De", "Positive TPC Pi vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
119+
histos.add("TPC/pTOF_Ka_Pi", "Positive TPC Ka vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
120+
histos.add("TPC/pTOF_Ka_Pr", "Positive TPC Ka vs Pr", {HistType::kTH2F, {ptBins, nSigmaBins}});
121+
histos.add("TPC/pTOF_Ka_El", "Positive TPC Ka vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
122+
histos.add("TPC/pTOF_Ka_De", "Positive TPC Ka vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
123+
histos.add("TPC/pTOF_Pr_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
124+
histos.add("TPC/pTOF_Pr_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
125+
histos.add("TPC/pTOF_Pr_El", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
126+
histos.add("TPC/pTOF_Pr_De", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
127+
histos.add("TPC/pTOF_El_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
128+
histos.add("TPC/pTOF_El_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
129+
histos.add("TPC/pTOF_El_Pr", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
130+
histos.add("TPC/pTOF_El_De", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
131+
histos.add("TPC/pTOF_De_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
132+
histos.add("TPC/pTOF_De_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
133+
histos.add("TPC/pTOF_De_Pr", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
134+
histos.add("TPC/pTOF_De_El", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
135+
136+
histos.add("TPC/nTOF_Pi_Ka", "Positive TPC Pi vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
137+
histos.add("TPC/nTOF_Pi_Pr", "Positive TPC Pi vs Pr", {HistType::kTH2F, {ptBins, nSigmaBins}});
138+
histos.add("TPC/nTOF_Pi_El", "Positive TPC Pi vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
139+
histos.add("TPC/nTOF_Pi_De", "Positive TPC Pi vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
140+
histos.add("TPC/nTOF_Ka_Pi", "Positive TPC Ka vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
141+
histos.add("TPC/nTOF_Ka_Pr", "Positive TPC Ka vs Pr", {HistType::kTH2F, {ptBins, nSigmaBins}});
142+
histos.add("TPC/nTOF_Ka_El", "Positive TPC Ka vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
143+
histos.add("TPC/nTOF_Ka_De", "Positive TPC Ka vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
144+
histos.add("TPC/nTOF_Pr_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
145+
histos.add("TPC/nTOF_Pr_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
146+
histos.add("TPC/nTOF_Pr_El", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
147+
histos.add("TPC/nTOF_Pr_De", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
148+
histos.add("TPC/nTOF_El_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
149+
histos.add("TPC/nTOF_El_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
150+
histos.add("TPC/nTOF_El_Pr", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
151+
histos.add("TPC/nTOF_El_De", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
152+
histos.add("TPC/nTOF_De_Pi", "Positive TPC Pr vs Pi", {HistType::kTH2F, {ptBins, nSigmaBins}});
153+
histos.add("TPC/nTOF_De_Ka", "Positive TPC Pr vs Ka", {HistType::kTH2F, {ptBins, nSigmaBins}});
154+
histos.add("TPC/nTOF_De_Pr", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
155+
histos.add("TPC/nTOF_De_El", "Positive TPC Pr vs El", {HistType::kTH2F, {ptBins, nSigmaBins}});
98156

99157
histos.add("TOF/pPi", "Positive TPC Pi vs TOF Pi vs pt", {HistType::kTH3F, {ptBins, ntofBins, ntofBins}});
100158
histos.add("TOF/nPi", "Negative TPC Pi vs TOF Pi vs pt", {HistType::kTH3F, {ptBins, ntofBins, ntofBins}});
@@ -139,21 +197,31 @@ struct sgPIDAnalyzer {
139197
histos.fill(HIST("TPC/pTPC_Ka_Pi"), track.pt(), track.tpcka());
140198
histos.fill(HIST("TPC/pTPC_Pr_Pi"), track.pt(), track.tpcpr());
141199
histos.fill(HIST("TPC/pTPC_El_Pi"), track.pt(), track.tpcel());
200+
histos.fill(HIST("TPC/pTPC_De_Pi"), track.pt(), track.tpcde());
142201
}
143202
if (std::abs(track.tpcka()) < 1) {
144203
histos.fill(HIST("TPC/pTPC_Pi_Ka"), track.pt(), track.tpcpi());
145204
histos.fill(HIST("TPC/pTPC_Pr_Ka"), track.pt(), track.tpcpr());
146205
histos.fill(HIST("TPC/pTPC_El_Ka"), track.pt(), track.tpcel());
206+
histos.fill(HIST("TPC/pTPC_De_Ka"), track.pt(), track.tpcde());
147207
}
148208
if (std::abs(track.tpcpr()) < 1) {
149209
histos.fill(HIST("TPC/pTPC_Pi_Pr"), track.pt(), track.tpcpi());
150210
histos.fill(HIST("TPC/pTPC_Ka_Pr"), track.pt(), track.tpcka());
151211
histos.fill(HIST("TPC/pTPC_El_Pr"), track.pt(), track.tpcel());
212+
histos.fill(HIST("TPC/pTPC_De_Pr"), track.pt(), track.tpcde());
152213
}
153214
if (std::abs(track.tpcel()) < 1) {
154215
histos.fill(HIST("TPC/pTPC_Pi_El"), track.pt(), track.tpcpi());
155216
histos.fill(HIST("TPC/pTPC_Ka_El"), track.pt(), track.tpcka());
156217
histos.fill(HIST("TPC/pTPC_Pr_El"), track.pt(), track.tpcpr());
218+
histos.fill(HIST("TPC/pTPC_De_El"), track.pt(), track.tpcde());
219+
}
220+
if (std::abs(track.tpcde()) < 1) {
221+
histos.fill(HIST("TPC/pTPC_Pi_De"), track.pt(), track.tpcpi());
222+
histos.fill(HIST("TPC/pTPC_Ka_De"), track.pt(), track.tpcka());
223+
histos.fill(HIST("TPC/pTPC_Pr_De"), track.pt(), track.tpcpr());
224+
histos.fill(HIST("TPC/pTPC_El_De"), track.pt(), track.tpcel());
157225
}
158226
} else {
159227
histos.fill(HIST("TPC/nTPC_Pi"), track.pt(), track.tpcpi());
@@ -169,21 +237,31 @@ struct sgPIDAnalyzer {
169237
histos.fill(HIST("TPC/nTPC_Ka_Pi"), track.pt(), track.tpcka());
170238
histos.fill(HIST("TPC/nTPC_Pr_Pi"), track.pt(), track.tpcpr());
171239
histos.fill(HIST("TPC/nTPC_El_Pi"), track.pt(), track.tpcel());
240+
histos.fill(HIST("TPC/nTPC_De_Pi"), track.pt(), track.tpcde());
172241
}
173242
if (std::abs(track.tpcka()) < 1) {
174243
histos.fill(HIST("TPC/nTPC_Pi_Ka"), track.pt(), track.tpcpi());
175244
histos.fill(HIST("TPC/nTPC_Pr_Ka"), track.pt(), track.tpcpr());
176245
histos.fill(HIST("TPC/nTPC_El_Ka"), track.pt(), track.tpcel());
246+
histos.fill(HIST("TPC/nTPC_De_Ka"), track.pt(), track.tpcde());
177247
}
178248
if (std::abs(track.tpcpr()) < 1) {
179249
histos.fill(HIST("TPC/nTPC_Pi_Pr"), track.pt(), track.tpcpi());
180250
histos.fill(HIST("TPC/nTPC_Ka_Pr"), track.pt(), track.tpcka());
181251
histos.fill(HIST("TPC/nTPC_El_Pr"), track.pt(), track.tpcel());
252+
histos.fill(HIST("TPC/nTPC_De_Pr"), track.pt(), track.tpcde());
182253
}
183254
if (std::abs(track.tpcel()) < 1) {
184255
histos.fill(HIST("TPC/nTPC_Pi_El"), track.pt(), track.tpcpi());
185256
histos.fill(HIST("TPC/nTPC_Ka_El"), track.pt(), track.tpcka());
186257
histos.fill(HIST("TPC/nTPC_Pr_El"), track.pt(), track.tpcpr());
258+
histos.fill(HIST("TPC/nTPC_De_El"), track.pt(), track.tpcde());
259+
}
260+
if (std::abs(track.tpcde()) < 1) {
261+
histos.fill(HIST("TPC/nTPC_Pi_De"), track.pt(), track.tpcpi());
262+
histos.fill(HIST("TPC/nTPC_Ka_De"), track.pt(), track.tpcka());
263+
histos.fill(HIST("TPC/nTPC_Pr_De"), track.pt(), track.tpcpr());
264+
histos.fill(HIST("TPC/nTPC_El_De"), track.pt(), track.tpcel());
187265
}
188266
}
189267
} else {
@@ -197,6 +275,36 @@ struct sgPIDAnalyzer {
197275
histos.fill(HIST("TOF/pHe"), track.pt(), track.tpcpr(), track.tofhe());
198276
histos.fill(HIST("TOF/pAl"), track.pt(), track.tpcel(), track.tofal());
199277
histos.fill(HIST("TOF/pMu"), track.pt(), track.tpcel(), track.tofmu());
278+
if (std::abs(track.tofpi()) < 1) {
279+
histos.fill(HIST("TOF/pTOF_Ka_Pi"), track.pt(), track.tofka());
280+
histos.fill(HIST("TOF/pTOF_Pr_Pi"), track.pt(), track.tofpr());
281+
histos.fill(HIST("TOF/pTOF_El_Pi"), track.pt(), track.tofel());
282+
histos.fill(HIST("TPF/pTOF_De_Pi"), track.pt(), track.tofde());
283+
}
284+
if (std::abs(track.tofka()) < 1) {
285+
histos.fill(HIST("TOF/pTOF_Pi_Ka"), track.pt(), track.tofpi());
286+
histos.fill(HIST("TOF/pTOF_Pr_Ka"), track.pt(), track.tofpr());
287+
histos.fill(HIST("TOF/pTOF_El_Ka"), track.pt(), track.tofel());
288+
histos.fill(HIST("TOF/pTOF_De_Ka"), track.pt(), track.tofde());
289+
}
290+
if (std::abs(track.tofpr()) < 1) {
291+
histos.fill(HIST("TOF/pTOF_Pi_Pr"), track.pt(), track.tofpi());
292+
histos.fill(HIST("TOF/pTOF_Ka_Pr"), track.pt(), track.tofka());
293+
histos.fill(HIST("TOF/pTOF_El_Pr"), track.pt(), track.tofel());
294+
histos.fill(HIST("TOF/pTOF_De_Pr"), track.pt(), track.tofde());
295+
}
296+
if (std::abs(track.tofel()) < 1) {
297+
histos.fill(HIST("TOF/pTOF_Pi_El"), track.pt(), track.tofpi());
298+
histos.fill(HIST("TOF/pTOF_Ka_El"), track.pt(), track.tofka());
299+
histos.fill(HIST("TOF/pTOF_Pr_El"), track.pt(), track.tofpr());
300+
histos.fill(HIST("TOF/pTOF_De_El"), track.pt(), track.tofde());
301+
}
302+
if (std::abs(track.tofde()) < 1) {
303+
histos.fill(HIST("TOF/pTOF_Pi_De"), track.pt(), track.tofpi());
304+
histos.fill(HIST("TOF/pTOF_Ka_De"), track.pt(), track.tofka());
305+
histos.fill(HIST("TOF/pTOF_Pr_De"), track.pt(), track.tofpr());
306+
histos.fill(HIST("TOF/pTOF_El_De"), track.pt(), track.tofel());
307+
}
200308
} else {
201309
histos.fill(HIST("TOF/nPi"), track.pt(), track.tpcpi(), track.tofpi());
202310
histos.fill(HIST("TOF/nKa"), track.pt(), track.tpcka(), track.tofka());
@@ -207,6 +315,36 @@ struct sgPIDAnalyzer {
207315
histos.fill(HIST("TOF/nHe"), track.pt(), track.tpcpr(), track.tofhe());
208316
histos.fill(HIST("TOF/nAl"), track.pt(), track.tpcel(), track.tofal());
209317
histos.fill(HIST("TOF/nMu"), track.pt(), track.tpcel(), track.tofmu());
318+
if (std::abs(track.tofpi()) < 1) {
319+
histos.fill(HIST("TOF/nTOF_Ka_Pi"), track.pt(), track.tofka());
320+
histos.fill(HIST("TOF/nTOF_Pr_Pi"), track.pt(), track.tofpr());
321+
histos.fill(HIST("TOF/nTOF_El_Pi"), track.pt(), track.tofel());
322+
histos.fill(HIST("TPF/nTOF_De_Pi"), track.pt(), track.tofde());
323+
}
324+
if (std::abs(track.tofka()) < 1) {
325+
histos.fill(HIST("TOF/nTOF_Pi_Ka"), track.pt(), track.tofpi());
326+
histos.fill(HIST("TOF/nTOF_Pr_Ka"), track.pt(), track.tofpr());
327+
histos.fill(HIST("TOF/nTOF_El_Ka"), track.pt(), track.tofel());
328+
histos.fill(HIST("TOF/nTOF_De_Ka"), track.pt(), track.tofde());
329+
}
330+
if (std::abs(track.tofpr()) < 1) {
331+
histos.fill(HIST("TOF/nTOF_Pi_Pr"), track.pt(), track.tofpi());
332+
histos.fill(HIST("TOF/nTOF_Ka_Pr"), track.pt(), track.tofka());
333+
histos.fill(HIST("TOF/nTOF_El_Pr"), track.pt(), track.tofel());
334+
histos.fill(HIST("TOF/nTOF_De_Pr"), track.pt(), track.tofde());
335+
}
336+
if (std::abs(track.tofel()) < 1) {
337+
histos.fill(HIST("TOF/nTOF_Pi_El"), track.pt(), track.tofpi());
338+
histos.fill(HIST("TOF/nTOF_Ka_El"), track.pt(), track.tofka());
339+
histos.fill(HIST("TOF/nTOF_Pr_El"), track.pt(), track.tofpr());
340+
histos.fill(HIST("TOF/nTOF_De_El"), track.pt(), track.tofde());
341+
}
342+
if (std::abs(track.tofde()) < 1) {
343+
histos.fill(HIST("TOF/nTOF_Pi_De"), track.pt(), track.tofpi());
344+
histos.fill(HIST("TOF/nTOF_Ka_De"), track.pt(), track.tofka());
345+
histos.fill(HIST("TOF/nTOF_Pr_De"), track.pt(), track.tofpr());
346+
histos.fill(HIST("TOF/nTOF_El_De"), track.pt(), track.tofel());
347+
}
210348
}
211349
}
212350
}

0 commit comments

Comments
 (0)