Skip to content

Commit dfdceaa

Browse files
alcalivafgrosa
alcaliva
authored and
fgrosa
committed
added histos for ka,pr,de (#5263)
1 parent 9be44df commit dfdceaa

1 file changed

Lines changed: 218 additions & 9 deletions

File tree

‎PWGLF/Tasks/Nuspex/AntimatterAbsorptionHMPID.cxx‎

Lines changed: 218 additions & 9 deletions
Original file line numberDiff line numberDiff line change
@@ -91,6 +91,72 @@ struct AntimatterAbsorptionHMPID {
9191
registryDA.add("Pi_Neg_ClsSize_8cm", "Pi_Neg_ClsSize_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
9292
registryDA.add("Pi_Neg_ClsSize_4cm", "Pi_Neg_ClsSize_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
9393
registryDA.add("Pi_Neg_momentum", "Pi_Neg_momentum", HistType::kTH2F, {{100, 0.0, 3.0, "#it{p}_{vtx} (GeV/#it{c})"}, {100, 0.0, 3.0, "#it{p}_{mhpid} (GeV/#it{c})"}});
94+
95+
// Kaon Pos
96+
registryDA.add("incomingKa_Pos_8cm", "incomingKa_Pos_8cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
97+
registryDA.add("incomingKa_Pos_4cm", "incomingKa_Pos_4cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
98+
registryDA.add("survivingKa_Pos_8cm", "survivingKa_Pos_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
99+
registryDA.add("survivingKa_Pos_4cm", "survivingKa_Pos_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
100+
registryDA.add("Ka_Pos_Q_8cm", "Ka_Pos_Q_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
101+
registryDA.add("Ka_Pos_Q_4cm", "Ka_Pos_Q_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
102+
registryDA.add("Ka_Pos_ClsSize_8cm", "Ka_Pos_ClsSize_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
103+
registryDA.add("Ka_Pos_ClsSize_4cm", "Ka_Pos_ClsSize_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
104+
registryDA.add("Ka_Pos_momentum", "Ka_Pos_momentum", HistType::kTH2F, {{100, 0.0, 3.0, "#it{p}_{vtx} (GeV/#it{c})"}, {100, 0.0, 3.0, "#it{p}_{mhpid} (GeV/#it{c})"}});
105+
106+
// Kaon Neg
107+
registryDA.add("incomingKa_Neg_8cm", "incomingKa_Neg_8cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
108+
registryDA.add("incomingKa_Neg_4cm", "incomingKa_Neg_4cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
109+
registryDA.add("survivingKa_Neg_8cm", "survivingKa_Neg_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
110+
registryDA.add("survivingKa_Neg_4cm", "survivingKa_Neg_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
111+
registryDA.add("Ka_Neg_Q_8cm", "Ka_Neg_Q_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
112+
registryDA.add("Ka_Neg_Q_4cm", "Ka_Neg_Q_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
113+
registryDA.add("Ka_Neg_ClsSize_8cm", "Ka_Neg_ClsSize_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
114+
registryDA.add("Ka_Neg_ClsSize_4cm", "Ka_Neg_ClsSize_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
115+
registryDA.add("Ka_Neg_momentum", "Ka_Neg_momentum", HistType::kTH2F, {{100, 0.0, 3.0, "#it{p}_{vtx} (GeV/#it{c})"}, {100, 0.0, 3.0, "#it{p}_{mhpid} (GeV/#it{c})"}});
116+
117+
// Proton Pos
118+
registryDA.add("incomingPr_Pos_8cm", "incomingPr_Pos_8cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
119+
registryDA.add("incomingPr_Pos_4cm", "incomingPr_Pos_4cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
120+
registryDA.add("survivingPr_Pos_8cm", "survivingPr_Pos_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
121+
registryDA.add("survivingPr_Pos_4cm", "survivingPr_Pos_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
122+
registryDA.add("Pr_Pos_Q_8cm", "Pr_Pos_Q_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
123+
registryDA.add("Pr_Pos_Q_4cm", "Pr_Pos_Q_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
124+
registryDA.add("Pr_Pos_ClsSize_8cm", "Pr_Pos_ClsSize_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
125+
registryDA.add("Pr_Pos_ClsSize_4cm", "Pr_Pos_ClsSize_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
126+
registryDA.add("Pr_Pos_momentum", "Pr_Pos_momentum", HistType::kTH2F, {{100, 0.0, 3.0, "#it{p}_{vtx} (GeV/#it{c})"}, {100, 0.0, 3.0, "#it{p}_{mhpid} (GeV/#it{c})"}});
127+
128+
// Proton Neg
129+
registryDA.add("incomingPr_Neg_8cm", "incomingPr_Neg_8cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
130+
registryDA.add("incomingPr_Neg_4cm", "incomingPr_Neg_4cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
131+
registryDA.add("survivingPr_Neg_8cm", "survivingPr_Neg_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
132+
registryDA.add("survivingPr_Neg_4cm", "survivingPr_Neg_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
133+
registryDA.add("Pr_Neg_Q_8cm", "Pr_Neg_Q_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
134+
registryDA.add("Pr_Neg_Q_4cm", "Pr_Neg_Q_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
135+
registryDA.add("Pr_Neg_ClsSize_8cm", "Pr_Neg_ClsSize_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
136+
registryDA.add("Pr_Neg_ClsSize_4cm", "Pr_Neg_ClsSize_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
137+
registryDA.add("Pr_Neg_momentum", "Pr_Neg_momentum", HistType::kTH2F, {{100, 0.0, 3.0, "#it{p}_{vtx} (GeV/#it{c})"}, {100, 0.0, 3.0, "#it{p}_{mhpid} (GeV/#it{c})"}});
138+
139+
// Deuteron Pos
140+
registryDA.add("incomingDe_Pos_8cm", "incomingDe_Pos_8cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
141+
registryDA.add("incomingDe_Pos_4cm", "incomingDe_Pos_4cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
142+
registryDA.add("survivingDe_Pos_8cm", "survivingDe_Pos_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
143+
registryDA.add("survivingDe_Pos_4cm", "survivingDe_Pos_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
144+
registryDA.add("De_Pos_Q_8cm", "De_Pos_Q_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
145+
registryDA.add("De_Pos_Q_4cm", "De_Pos_Q_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
146+
registryDA.add("De_Pos_ClsSize_8cm", "De_Pos_ClsSize_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
147+
registryDA.add("De_Pos_ClsSize_4cm", "De_Pos_ClsSize_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
148+
registryDA.add("De_Pos_momentum", "De_Pos_momentum", HistType::kTH2F, {{100, 0.0, 3.0, "#it{p}_{vtx} (GeV/#it{c})"}, {100, 0.0, 3.0, "#it{p}_{mhpid} (GeV/#it{c})"}});
149+
150+
// Deuteron Neg
151+
registryDA.add("incomingDe_Neg_8cm", "incomingDe_Neg_8cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
152+
registryDA.add("incomingDe_Neg_4cm", "incomingDe_Neg_4cm", HistType::kTH1F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}});
153+
registryDA.add("survivingDe_Neg_8cm", "survivingDe_Neg_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
154+
registryDA.add("survivingDe_Neg_4cm", "survivingDe_Neg_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {300, 0.0, 30.0, "#Delta R (cm)"}});
155+
registryDA.add("De_Neg_Q_8cm", "De_Neg_Q_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
156+
registryDA.add("De_Neg_Q_4cm", "De_Neg_Q_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Q (ADC)"}});
157+
registryDA.add("De_Neg_ClsSize_8cm", "De_Neg_ClsSize_8cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
158+
registryDA.add("De_Neg_ClsSize_4cm", "De_Neg_ClsSize_4cm", HistType::kTH2F, {{290, 0.1, 3.0, "#it{p} (GeV/#it{c})"}, {200, 0.0, 2000.0, "Cls size"}});
159+
registryDA.add("De_Neg_momentum", "De_Neg_momentum", HistType::kTH2F, {{100, 0.0, 3.0, "#it{p}_{vtx} (GeV/#it{c})"}, {100, 0.0, 3.0, "#it{p}_{mhpid} (GeV/#it{c})"}});
94160
}
95161

96162
// Single-Track Selection
@@ -103,8 +169,6 @@ struct AntimatterAbsorptionHMPID {
103169
return false;
104170
if (!track.hasTOF())
105171
return false;
106-
// if (!track.has_hmpid())
107-
// return false;
108172
if (!track.passedITSRefit())
109173
return false;
110174
if (!track.passedTPCRefit())
@@ -152,6 +216,54 @@ struct AntimatterAbsorptionHMPID {
152216
return true;
153217
}
154218

219+
// Particle Identification (Kaons)
220+
template <typename kaonCandidate>
221+
bool passedKaonSelection(const kaonCandidate& track)
222+
{
223+
if (track.tpcNSigmaKa() < nsigmaTPCMin)
224+
return false;
225+
if (track.tpcNSigmaKa() > nsigmaTPCMax)
226+
return false;
227+
if (track.tofNSigmaKa() < nsigmaTOFMin)
228+
return false;
229+
if (track.tofNSigmaKa() > nsigmaTOFMax)
230+
return false;
231+
232+
return true;
233+
}
234+
235+
// Particle Identification (Protons)
236+
template <typename protonCandidate>
237+
bool passedProtonSelection(const protonCandidate& track)
238+
{
239+
if (track.tpcNSigmaPr() < nsigmaTPCMin)
240+
return false;
241+
if (track.tpcNSigmaPr() > nsigmaTPCMax)
242+
return false;
243+
if (track.tofNSigmaPr() < nsigmaTOFMin)
244+
return false;
245+
if (track.tofNSigmaPr() > nsigmaTOFMax)
246+
return false;
247+
248+
return true;
249+
}
250+
251+
// Particle Identification (Deuterons)
252+
template <typename deuteronCandidate>
253+
bool passedDeuteronSelection(const deuteronCandidate& track)
254+
{
255+
if (track.tpcNSigmaDe() < nsigmaTPCMin)
256+
return false;
257+
if (track.tpcNSigmaDe() > nsigmaTPCMax)
258+
return false;
259+
if (track.tofNSigmaDe() < nsigmaTOFMin)
260+
return false;
261+
if (track.tofNSigmaDe() > nsigmaTOFMax)
262+
return false;
263+
264+
return true;
265+
}
266+
155267
// Info for TPC PID
156268
using PidInfoTPC = soa::Join<aod::pidTPCLfFullPi, aod::pidTPCLfFullKa,
157269
aod::pidTPCLfFullPr, aod::pidTPCLfFullDe,
@@ -198,11 +310,6 @@ struct AntimatterAbsorptionHMPID {
198310
registryDA.fill(HIST("hmpidXYneg"), hmpid.hmpidXMip(), hmpid.hmpidYMip());
199311
}
200312

201-
// Particle Identification
202-
bool passedPionSel = false;
203-
if (passedPionSelection(track))
204-
passedPionSel = true;
205-
206313
// Absorber
207314
bool hmpidAbs8cm = true;
208315
bool hmpidAbs4cm = true;
@@ -213,7 +320,7 @@ struct AntimatterAbsorptionHMPID {
213320
float dr = sqrt(dx * dx + dy * dy);
214321

215322
// Fill Histograms for Positive Pions
216-
if (passedPionSel && track.sign() > 0) {
323+
if (passedPionSelection(track) && track.sign() > 0) {
217324

218325
if (hmpidAbs8cm) {
219326
registryDA.fill(HIST("incomingPi_Pos_8cm"), hmpid.hmpidMom());
@@ -230,7 +337,7 @@ struct AntimatterAbsorptionHMPID {
230337
}
231338

232339
// Fill Histograms for Negative Pions
233-
if (passedPionSel && track.sign() < 0) {
340+
if (passedPionSelection(track) && track.sign() < 0) {
234341

235342
if (hmpidAbs8cm) {
236343
registryDA.fill(HIST("incomingPi_Neg_8cm"), hmpid.hmpidMom());
@@ -245,6 +352,108 @@ struct AntimatterAbsorptionHMPID {
245352
registryDA.fill(HIST("Pi_Neg_ClsSize_4cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
246353
}
247354
}
355+
356+
// Fill Histograms for Positive Kaons
357+
if (passedKaonSelection(track) && track.sign() > 0) {
358+
359+
if (hmpidAbs8cm) {
360+
registryDA.fill(HIST("incomingKa_Pos_8cm"), hmpid.hmpidMom());
361+
registryDA.fill(HIST("survivingKa_Pos_8cm"), hmpid.hmpidMom(), dr);
362+
registryDA.fill(HIST("Ka_Pos_Q_8cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
363+
registryDA.fill(HIST("Ka_Pos_ClsSize_8cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
364+
}
365+
if (hmpidAbs4cm) {
366+
registryDA.fill(HIST("incomingKa_Pos_4cm"), hmpid.hmpidMom());
367+
registryDA.fill(HIST("survivingKa_Pos_4cm"), hmpid.hmpidMom(), dr);
368+
registryDA.fill(HIST("Ka_Pos_Q_4cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
369+
registryDA.fill(HIST("Ka_Pos_ClsSize_4cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
370+
}
371+
}
372+
373+
// Fill Histograms for Negative Kaons
374+
if (passedKaonSelection(track) && track.sign() < 0) {
375+
376+
if (hmpidAbs8cm) {
377+
registryDA.fill(HIST("incomingKa_Neg_8cm"), hmpid.hmpidMom());
378+
registryDA.fill(HIST("survivingKa_Neg_8cm"), hmpid.hmpidMom(), dr);
379+
registryDA.fill(HIST("Ka_Neg_Q_8cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
380+
registryDA.fill(HIST("Ka_Neg_ClsSize_8cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
381+
}
382+
if (hmpidAbs4cm) {
383+
registryDA.fill(HIST("incomingKa_Neg_4cm"), hmpid.hmpidMom());
384+
registryDA.fill(HIST("survivingKa_Neg_4cm"), hmpid.hmpidMom(), dr);
385+
registryDA.fill(HIST("Ka_Neg_Q_4cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
386+
registryDA.fill(HIST("Ka_Neg_ClsSize_4cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
387+
}
388+
}
389+
390+
// Fill Histograms for Positive Protons
391+
if (passedProtonSelection(track) && track.sign() > 0) {
392+
393+
if (hmpidAbs8cm) {
394+
registryDA.fill(HIST("incomingPr_Pos_8cm"), hmpid.hmpidMom());
395+
registryDA.fill(HIST("survivingPr_Pos_8cm"), hmpid.hmpidMom(), dr);
396+
registryDA.fill(HIST("Pr_Pos_Q_8cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
397+
registryDA.fill(HIST("Pr_Pos_ClsSize_8cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
398+
}
399+
if (hmpidAbs4cm) {
400+
registryDA.fill(HIST("incomingPr_Pos_4cm"), hmpid.hmpidMom());
401+
registryDA.fill(HIST("survivingPr_Pos_4cm"), hmpid.hmpidMom(), dr);
402+
registryDA.fill(HIST("Pr_Pos_Q_4cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
403+
registryDA.fill(HIST("Pr_Pos_ClsSize_4cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
404+
}
405+
}
406+
407+
// Fill Histograms for Negative Protons
408+
if (passedProtonSelection(track) && track.sign() < 0) {
409+
410+
if (hmpidAbs8cm) {
411+
registryDA.fill(HIST("incomingPr_Neg_8cm"), hmpid.hmpidMom());
412+
registryDA.fill(HIST("survivingPr_Neg_8cm"), hmpid.hmpidMom(), dr);
413+
registryDA.fill(HIST("Pr_Neg_Q_8cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
414+
registryDA.fill(HIST("Pr_Neg_ClsSize_8cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
415+
}
416+
if (hmpidAbs4cm) {
417+
registryDA.fill(HIST("incomingPr_Neg_4cm"), hmpid.hmpidMom());
418+
registryDA.fill(HIST("survivingPr_Neg_4cm"), hmpid.hmpidMom(), dr);
419+
registryDA.fill(HIST("Pr_Neg_Q_4cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
420+
registryDA.fill(HIST("Pr_Neg_ClsSize_4cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
421+
}
422+
}
423+
424+
// Fill Histograms for Positive Deuterons
425+
if (passedDeuteronSelection(track) && track.sign() > 0) {
426+
427+
if (hmpidAbs8cm) {
428+
registryDA.fill(HIST("incomingDe_Pos_8cm"), hmpid.hmpidMom());
429+
registryDA.fill(HIST("survivingDe_Pos_8cm"), hmpid.hmpidMom(), dr);
430+
registryDA.fill(HIST("De_Pos_Q_8cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
431+
registryDA.fill(HIST("De_Pos_ClsSize_8cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
432+
}
433+
if (hmpidAbs4cm) {
434+
registryDA.fill(HIST("incomingDe_Pos_4cm"), hmpid.hmpidMom());
435+
registryDA.fill(HIST("survivingDe_Pos_4cm"), hmpid.hmpidMom(), dr);
436+
registryDA.fill(HIST("De_Pos_Q_4cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
437+
registryDA.fill(HIST("De_Pos_ClsSize_4cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
438+
}
439+
}
440+
441+
// Fill Histograms for Negative Deuterons
442+
if (passedDeuteronSelection(track) && track.sign() < 0) {
443+
444+
if (hmpidAbs8cm) {
445+
registryDA.fill(HIST("incomingDe_Neg_8cm"), hmpid.hmpidMom());
446+
registryDA.fill(HIST("survivingDe_Neg_8cm"), hmpid.hmpidMom(), dr);
447+
registryDA.fill(HIST("De_Neg_Q_8cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
448+
registryDA.fill(HIST("De_Neg_ClsSize_8cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
449+
}
450+
if (hmpidAbs4cm) {
451+
registryDA.fill(HIST("incomingDe_Neg_4cm"), hmpid.hmpidMom());
452+
registryDA.fill(HIST("survivingDe_Neg_4cm"), hmpid.hmpidMom(), dr);
453+
registryDA.fill(HIST("De_Neg_Q_4cm"), hmpid.hmpidMom(), hmpid.hmpidQMip());
454+
registryDA.fill(HIST("De_Neg_ClsSize_4cm"), hmpid.hmpidMom(), hmpid.hmpidClusSize());
455+
}
456+
}
248457
}
249458
}
250459
PROCESS_SWITCH(AntimatterAbsorptionHMPID, processData, "process data", true);

0 commit comments

Comments
 (0)