Skip to content

Commit 788bb0e

Browse files
author
Arkaprabha Saha
committed
updated the pT calculation for helium3 and fixed some variable names
1 parent 383d819 commit 788bb0e

File tree

1 file changed

+5
-5
lines changed

1 file changed

+5
-5
lines changed

PWGLF/Tasks/Nuspex/AntiNucleiTask.cxx

Lines changed: 5 additions & 5 deletions
Original file line numberDiff line numberDiff line change
@@ -140,25 +140,25 @@ struct AntiNucleiTask {
140140

141141
double expBethe{tpc::BetheBlochAleph(static_cast<double>(track.tpcInnerParam()), cfgBetheBlochParams->get("p0"), cfgBetheBlochParams->get("p1"), cfgBetheBlochParams->get("p2"), cfgBetheBlochParams->get("p3"), cfgBetheBlochParams->get("p4"))};
142142
double expSigma{expBethe * cfgBetheBlochParams->get("resolution")};
143-
float tpcNSigmaDeuteron = static_cast<float>((track.tpcSignal() - expBethe) / expSigma);
143+
float tpcNSigma = static_cast<float>((track.tpcSignal() - expBethe) / expSigma);
144144

145-
float pt = track.sign() > 0 ? track.pt() : -track.pt();
145+
float pt = track.sign() > 0 ? 2 * track.pt() : -2 * track.pt();
146146
// Filling histograms with track data before applying any cuts.
147147
histos.fill(HIST("RawEta"), track.eta());
148148
histos.fill(HIST("RawPhi"), track.phi());
149149
histos.fill(HIST("RawPt"), pt);
150-
histos.fill(HIST("RawtpcNSigma"), collision.centFT0C(), pt, tpcNSigmaDeuteron);
150+
histos.fill(HIST("RawtpcNSigma"), collision.centFT0C(), pt, tpcNSigma);
151151
histos.fill(HIST("RawtofNSigma"), collision.centFT0C(), pt, track.tofNSigmaDe());
152152

153153
// If the track is good, fill the "after cuts" histograms.
154154
if (isGoodTrack(track)) {
155155
histos.fill(HIST("Eta"), track.eta());
156156
histos.fill(HIST("Phi"), track.phi());
157157
histos.fill(HIST("Pt"), pt);
158-
histos.fill(HIST("tpcNSigma"), collision.centFT0C(), pt, tpcNSigmaDeuteron);
158+
histos.fill(HIST("tpcNSigma"), collision.centFT0C(), pt, tpcNSigma);
159159
histos.fill(HIST("TpcSignal"), track.tpcInnerParam(), track.tpcSignal());
160160

161-
if (std::abs(tpcNSigmaDeuteron) < tpcNSigmaCut) {
161+
if (std::abs(tpcNSigma) < tpcNSigmaCut) {
162162
histos.fill(HIST("tofNSigma"), collision.centFT0C(), pt, track.tofNSigmaDe());
163163
}
164164
}

0 commit comments

Comments
 (0)