Skip to content

Commit 7367b98

Browse files
jtanwar2212jtanwar2212
authored
Update nucleitpcpbpb.cxx
1 parent 44786a3 commit 7367b98

File tree

1 file changed

+11
-11
lines changed

1 file changed

+11
-11
lines changed

‎PWGLF/Tasks/Nuspex/nucleitpcpbpb.cxx‎

Lines changed: 11 additions & 11 deletions
Original file line numberDiff line numberDiff line change
@@ -67,7 +67,7 @@ constexpr double betheBlochDefault[nParticles][nBetheParams]{
6767
{-126.557359, -0.858569, 1.111643, 1.210323, 2.656374, 0.09}}; // alpha
6868
const int nTrkSettings = 19;
6969
static const std::vector<std::string> trackPIDsettingsNames{"useBBparams", "minITSnCls", "minITSnClscos", "minTPCnCls", "maxTPCchi2", "minTPCchi2", "maxITSchi2", "minRigidity", "maxRigidity", "maxTPCnSigma", "TOFrequiredabove", "minTOFmass", "maxTOFmass", "maxDcaXY", "maxDcaZ", "minITSclsSize", "maxITSclsSize", "minTPCnClsCrossedRows", "minReqClusterITSib"};
70-
constexpr double trackPIDsettings[nParticles][nTrkSettings]{
70+
constexpr double kTrackPIDSettings[nParticles][nTrkSettings]{
7171
{0, 0, 4, 60, 4.0, 0.5, 100, 0.15, 1.2, 2.5, -1, 0, 100, 2., 2., 0., 1000, 70, 1},
7272
{1, 0, 4, 70, 4.0, 0.5, 100, 0.20, 4.0, 3.0, -1, 0, 100, 2., 2., 0., 1000, 70, 1},
7373
{1, 0, 4, 70, 4.0, 0.5, 100, 0.50, 5.0, 3.0, -1, 0, 100, 2., 2., 0., 1000, 70, 1},
@@ -137,7 +137,7 @@ struct NucleitpcPbPb {
137137
Configurable<bool> cfgRequireMCposZ{"cfgRequireMCposZ", true, "Require beta plot"};
138138

139139
Configurable<LabeledArray<double>> cfgBetheBlochParams{"cfgBetheBlochParams", {betheBlochDefault[0], nParticles, nBetheParams, particleNames, betheBlochParNames}, "TPC Bethe-Bloch parameterisation for light nuclei"};
140-
Configurable<LabeledArray<double>> cfgTrackPIDsettings{"cfgTrackPIDsettings", {trackPIDsettings[0], nParticles, nTrkSettings, particleNames, trackPIDsettingsNames}, "track selection and PID criteria"};
140+
Configurable<LabeledArray<double>> cfgTrackPIDsettings{"cfgTrackPIDsettings", {kTrackPIDSettings[0], nParticles, nTrkSettings, particleNames, trackPIDsettingsNames}, "track selection and PID criteria"};
141141
Configurable<bool> cfgFillDeDxWithCut{"cfgFillDeDxWithCut", true, "Fill with cut beth bloch"};
142142
Configurable<bool> cfgFillnsigma{"cfgFillnsigma", false, "Fill n-sigma histograms"};
143143
Configurable<bool> cfgFillmass{"cfgFillmass", false, "Fill mass histograms"};
@@ -177,8 +177,8 @@ struct NucleitpcPbPb {
177177
int mRunNumber, occupancy;
178178
float dBz, momn;
179179
TRandom3 rand;
180-
float He3 = 4;
181-
float He4 = 5;
180+
float he3 = 4;
181+
float he4 = 5;
182182
//----------------------------------------------------------------------------------------------------------------
183183
void init(InitContext const&)
184184
{
@@ -295,7 +295,7 @@ struct NucleitpcPbPb {
295295
float ptMomn;
296296
setTrackParCov(track, mTrackParCov);
297297
mTrackParCov.setPID(track.pidForTracking());
298-
ptMomn = (i == He3 || i == He4) ? 2 * mTrackParCov.getPt() : mTrackParCov.getPt();
298+
ptMomn = (i == he3 || i == he4) ? 2 * mTrackParCov.getPt() : mTrackParCov.getPt();
299299
bool insideDCAxy = (std::abs(track.dcaXY()) <= (cfgTrackPIDsettings->get(i, "maxDcaXY") * (0.0105f + 0.0350f / std::pow(ptMomn, 1.1f)))); // o2-linter: disable=magic-number (To be checked)
300300
if ((!(insideDCAxy) || std::abs(track.dcaZ()) > DCAzSigma(ptMomn, cfgTrackPIDsettings->get(i, "maxDcaZ"))) && cfgDCAwithptRequire)
301301
continue;
@@ -319,7 +319,7 @@ struct NucleitpcPbPb {
319319
if (!track.hasTOF() && cfgTrackPIDsettings->get(i, "TOFrequiredabove") < 1)
320320
continue;
321321
float beta{o2::pid::tof::Beta::GetBeta(track)};
322-
float charge{1.f + static_cast<float>(i == He3 || i == He4)};
322+
float charge{1.f + static_cast<float>(i == he3 || i == he4)};
323323
float tofMasses = getRigidity(track) * charge * std::sqrt(1.f / (beta * beta) - 1.f);
324324
if ((getRigidity(track) > cfgTrackPIDsettings->get(i, "TOFrequiredabove") && (tofMasses < cfgTrackPIDsettings->get(i, "minTOFmass") || tofMasses > cfgTrackPIDsettings->get(i, "maxTOFmass"))) && cfgmassRequire)
325325
continue;
@@ -489,7 +489,7 @@ struct NucleitpcPbPb {
489489
float ptMomn;
490490
setTrackParCov(track, mTrackParCov);
491491
mTrackParCov.setPID(track.pidForTracking());
492-
ptMomn = (i == He3 || i == He4) ? 2 * mTrackParCov.getPt() : mTrackParCov.getPt();
492+
ptMomn = (i == he3 || i == he4) ? 2 * mTrackParCov.getPt() : mTrackParCov.getPt();
493493
bool insideDCAxy = (std::abs(track.dcaXY()) <= (cfgTrackPIDsettings->get(i, "maxDcaXY") * (0.0105f + 0.0350f / std::pow(ptMomn, 1.1f)))); // o2-linter: disable=magic-number (To be checked)
494494
if ((!(insideDCAxy) || std::abs(track.dcaZ()) > DCAzSigma(ptMomn, cfgTrackPIDsettings->get(i, "maxDcaZ"))) && cfgDCAwithptRequire)
495495
continue;
@@ -512,7 +512,7 @@ struct NucleitpcPbPb {
512512
if (!track.hasTOF() && cfgTrackPIDsettings->get(i, "TOFrequiredabove") < 1)
513513
continue;
514514
float beta{o2::pid::tof::Beta::GetBeta(track)};
515-
float charge{1.f + static_cast<float>(i == He3 || i == He4)};
515+
float charge{1.f + static_cast<float>(i == he3 || i == he4)};
516516
float tofMasses = getRigidity(track) * charge * std::sqrt(1.f / (beta * beta) - 1.f);
517517
if ((getRigidity(track) > cfgTrackPIDsettings->get(i, "TOFrequiredabove") && (tofMasses < cfgTrackPIDsettings->get(i, "minTOFmass") || tofMasses > cfgTrackPIDsettings->get(i, "maxTOFmass"))) && cfgmassRequire)
518518
continue;
@@ -639,7 +639,7 @@ struct NucleitpcPbPb {
639639
float ptMomn;
640640
setTrackParCov(track, mTrackParCov);
641641
mTrackParCov.setPID(track.pidForTracking());
642-
ptMomn = (i == He3 || i == He4) ? 2 * mTrackParCov.getPt() : mTrackParCov.getPt();
642+
ptMomn = (i == he3 || i == he4) ? 2 * mTrackParCov.getPt() : mTrackParCov.getPt();
643643
if (track.sign() > 0) {
644644
hNsigmaPt[2 * species]->Fill(ptMomn, tpcNsigma);
645645
}
@@ -659,7 +659,7 @@ struct NucleitpcPbPb {
659659
float beta{o2::pid::tof::Beta::GetBeta(track)};
660660
if (beta <= 0.f || beta >= 1.f)
661661
return;
662-
float charge = (species == He3 || species == He4) ? 2.f : 1.f;
662+
float charge = (species == he3 || species == he4) ? 2.f : 1.f;
663663
float p = getRigidity(track); // assuming this is the momentum from inner TPC
664664
float massTOF = p * charge * std::sqrt(1.f / (beta * beta) - 1.f);
665665
// get PDG mass
@@ -668,7 +668,7 @@ struct NucleitpcPbPb {
668668
float ptMomn;
669669
setTrackParCov(track, mTrackParCov);
670670
mTrackParCov.setPID(track.pidForTracking());
671-
ptMomn = (species == He3 || species == He4) ? 2 * mTrackParCov.getPt() : mTrackParCov.getPt();
671+
ptMomn = (species == he3 || species == he4) ? 2 * mTrackParCov.getPt() : mTrackParCov.getPt();
672672
if (track.sign() > 0) {
673673
hmass[2 * species]->Fill(ptMomn, massDiff * massDiff);
674674
} else if (track.sign() < 0) {

0 commit comments

Comments
 (0)