diff --git a/PWGLF/Utils/nucleiUtils.h b/PWGLF/Utils/nucleiUtils.h index 10720befcb4..f75dab12a74 100644 --- a/PWGLF/Utils/nucleiUtils.h +++ b/PWGLF/Utils/nucleiUtils.h @@ -12,9 +12,6 @@ #ifndef PWGLF_UTILS_NUCLEIUTILS_H_ #define PWGLF_UTILS_NUCLEIUTILS_H_ -#include "Framework/HistogramRegistry.h" -#include "Framework/HistogramSpec.h" - #include "Common/DataModel/Centrality.h" #include "Common/DataModel/EventSelection.h" #include "Common/DataModel/PIDResponse.h" @@ -22,10 +19,12 @@ #include "Common/TableProducer/PID/pidTOFBase.h" #include "DataFormatsTPC/BetheBlochAleph.h" +#include "Framework/HistogramRegistry.h" +#include "Framework/HistogramSpec.h" +#include #include #include -#include using namespace o2; using namespace o2::framework; @@ -435,7 +434,8 @@ class PidManager template float getBetaTOF(const Ttrack& track) { - if (!track.hasTOF()) return -999.f; + if (!track.hasTOF()) + return -999.f; float beta = o2::pid::tof::Beta::GetBeta(track); return std::min(1.f - 1.e-6f, std::max(1.e-4f, beta)); /// sometimes beta > 1 or < 0, to be checked } @@ -443,7 +443,8 @@ class PidManager template float getMassTOF(const Ttrack& track) { - if (!track.hasTOF()) return -999.f; + if (!track.hasTOF()) + return -999.f; const float charge{1.f + static_cast(mSpecies == Species::kHe || mSpecies == Species::kAl)}; const float beta = getBetaTOF(track); return track.tpcInnerParam() * charge * std::sqrt(1.f / (beta * beta) - 1.f); @@ -452,15 +453,22 @@ class PidManager template float getNSigmaTOF(const Ttrack& track) { - if (!track.hasTOF()) return -999.f; + if (!track.hasTOF()) + return -999.f; switch (mSpecies) { - case Species::kPr: return track.tofNSigmaPr(); - case Species::kDe: return track.tofNSigmaDe(); - case Species::kTr: return track.tofNSigmaTr(); - case Species::kHe: return track.tofNSigmaHe(); - case Species::kAl: return track.tofNSigmaAl(); - default: return -999.f; + case Species::kPr: + return track.tofNSigmaPr(); + case Species::kDe: + return track.tofNSigmaDe(); + case Species::kTr: + return track.tofNSigmaTr(); + case Species::kHe: + return track.tofNSigmaHe(); + case Species::kAl: + return track.tofNSigmaAl(); + default: + return -999.f; } } @@ -497,7 +505,8 @@ class PidManager // TPC float getExpectedTPCsignal(const float p) { - if (!mUseTpcCentralCalibration) return -999.f; + if (!mUseTpcCentralCalibration) + return -999.f; float pScaled = p * mMomScaling[0] + mMomScaling[1]; float betaGamma = pScaled / masses[mSpecies]; @@ -526,12 +535,18 @@ class PidManager float getNSigmaTPCcentral(const Ttrack& track) { switch (mSpecies) { - case Species::kPr: return track.tpcNSigmaPr(); - case Species::kDe: return track.tpcNSigmaDe(); - case Species::kTr: return track.tpcNSigmaTr(); - case Species::kHe: return track.tpcNSigmaHe(); - case Species::kAl: return track.tpcNSigmaAl(); - default: return -999.f; + case Species::kPr: + return track.tpcNSigmaPr(); + case Species::kDe: + return track.tpcNSigmaDe(); + case Species::kTr: + return track.tpcNSigmaTr(); + case Species::kHe: + return track.tpcNSigmaHe(); + case Species::kAl: + return track.tpcNSigmaAl(); + default: + return -999.f; } }