diff --git a/PWGLF/Tasks/Nuspex/AntiNucleiTask.cxx b/PWGLF/Tasks/Nuspex/AntiNucleiTask.cxx index 3cf9fcb8758..96a16bd79e9 100644 --- a/PWGLF/Tasks/Nuspex/AntiNucleiTask.cxx +++ b/PWGLF/Tasks/Nuspex/AntiNucleiTask.cxx @@ -69,7 +69,7 @@ struct AntiNucleiTask { ConfigurableAxis phiAxis{"phiAxis", {70, 0.f, 7.f}, "#phi"}; ConfigurableAxis zVtxAxis{"zVtxAxis", {100, -20.f, 20.f}, "Primary Vertex z (cm)"}; ConfigurableAxis nSigmaAxis{"nSigmaAxis", {50, -5.f, 5.f}, "N_{#sigma}"}; - ConfigurableAxis ptAxis{"ptAxis", {100, -5.0f, 5.0f}, "p_{T} (GeV/c)"}; + ConfigurableAxis ptAxis{"ptAxis", {200, -10.0f, 10.0f}, "p_{T} (GeV/c)"}; ConfigurableAxis centAxis{"centAxis", {100, 0, 100.0f}, "Centrality"}; ConfigurableAxis momAxis{"momAxis", {5.e2, 0.f, 5.f}, "momentum axis binning"}; ConfigurableAxis tpcAxis{"tpcAxis", {4.e2, 0.f, 4.e3f}, "tpc signal axis binning"}; @@ -140,14 +140,14 @@ struct AntiNucleiTask { double expBethe{tpc::BetheBlochAleph(static_cast(track.tpcInnerParam()), cfgBetheBlochParams->get("p0"), cfgBetheBlochParams->get("p1"), cfgBetheBlochParams->get("p2"), cfgBetheBlochParams->get("p3"), cfgBetheBlochParams->get("p4"))}; double expSigma{expBethe * cfgBetheBlochParams->get("resolution")}; - float tpcNSigmaDeuteron = static_cast((track.tpcSignal() - expBethe) / expSigma); + float tpcNSigma = static_cast((track.tpcSignal() - expBethe) / expSigma); - float pt = track.sign() > 0 ? track.pt() : -track.pt(); + float pt = track.sign() > 0 ? 2 * track.pt() : -2 * track.pt(); // Filling histograms with track data before applying any cuts. histos.fill(HIST("RawEta"), track.eta()); histos.fill(HIST("RawPhi"), track.phi()); histos.fill(HIST("RawPt"), pt); - histos.fill(HIST("RawtpcNSigma"), collision.centFT0C(), pt, tpcNSigmaDeuteron); + histos.fill(HIST("RawtpcNSigma"), collision.centFT0C(), pt, tpcNSigma); histos.fill(HIST("RawtofNSigma"), collision.centFT0C(), pt, track.tofNSigmaDe()); // If the track is good, fill the "after cuts" histograms. @@ -155,10 +155,10 @@ struct AntiNucleiTask { histos.fill(HIST("Eta"), track.eta()); histos.fill(HIST("Phi"), track.phi()); histos.fill(HIST("Pt"), pt); - histos.fill(HIST("tpcNSigma"), collision.centFT0C(), pt, tpcNSigmaDeuteron); + histos.fill(HIST("tpcNSigma"), collision.centFT0C(), pt, tpcNSigma); histos.fill(HIST("TpcSignal"), track.tpcInnerParam(), track.tpcSignal()); - if (std::abs(tpcNSigmaDeuteron) < tpcNSigmaCut) { + if (std::abs(tpcNSigma) < tpcNSigmaCut) { histos.fill(HIST("tofNSigma"), collision.centFT0C(), pt, track.tofNSigmaDe()); } }