Skip to content

Commit 7bba362

Browse files
giovannimalfattoregiovannimalfattore
authored
[PWGLF] LightNucleiTask - Add calibration friendly mode and new TPC cuts (#1898)
* [PWGLF] LightNucleiTask - Fix TPC cuts * [PWGLF] LightNucleiTask - Add calibration friendly mode and new TPC cuts
1 parent fc0ad64 commit 7bba362

File tree

1 file changed

+36
-33
lines changed

1 file changed

+36
-33
lines changed

‎PWGLF/Tasks/LFNucleiBATask.cxx‎

Lines changed: 36 additions & 33 deletions
Original file line numberDiff line numberDiff line change
@@ -58,7 +58,8 @@ struct LFNucleiBATask {
5858
ConfigurableAxis binsdEdx{"binsdEdx", {1000, 0.f, 1000.f}, ""};
5959
Configurable<bool> enableEvTimeSplitting{"enableEvTimeSplitting", true, "Flag to enable histograms splitting depending on the Event Time used"};
6060
Configurable<bool> enableDebug{"enableDebug", false, "Flag to enable histograms for debug"};
61-
Configurable<bool> enableStrongCut{"enableStrongCut", false, "Flag to change | NSigma TPC(nucl)| < nSigmaTPC --> NOT | NSigma TPC(p)| > 3"};
61+
Configurable<bool> enableStrongCut{"enableStrongCut", false, "Flag to change | NSigma TPC(nucl)| < nSigmaTPC --> NOT | NSigma TPC(p)| > nStrongCut"};
62+
Configurable<float> nsigmaTPCStrongCut{"nsigmaTPCStrongCut", 3.f, "Value of the Nsigma TPC (Strong) proton cut, if enabled"};
6263
Configurable<int> useHasTRDConfig{"useHasTRDConfig", 0, "No selections on TRD (0); With TRD (1); Without TRD (2)"};
6364

6465
static constexpr int PDGPion = 211;
@@ -311,7 +312,8 @@ struct LFNucleiBATask {
311312
}
312313

313314
// Bethe-Bloch TPC distribution and Beta vs pT TOF distribution
314-
histos.add<TH2>("tracks/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {8000, 0.0, 1E3}});
315+
histos.add<TH2>("tracks/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {dedxAxis}});
316+
315317
histos.add<TH2>("tracks/h2TOFbetaVsP", "TOF #beta vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TOF #beta", HistType::kTH2F, {{1000, -5.f, 5.f}, {1200, 0.0, 1.2}});
316318
histos.add<TH2>("tracks/h2TOFbetaVsP_debug", "TOF #beta vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TOF #beta", HistType::kTH2F, {{1000, -5.f, 5.f}, {1200, 0.0, 1.2}});
317319
histos.add<TH2>("tracks/h2withTPCProtonCutTOFbetaVsP", "TOF #beta (with N_{#sigma(TPC(p))} cut) vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TOF #beta (with N_{#sigma(TPC(p))} cut)", HistType::kTH2F, {{1000, -5.f, 5.f}, {1200, 0.0, 1.2}});
@@ -394,10 +396,10 @@ struct LFNucleiBATask {
394396
// TOF EvTime Splitting plots
395397
if (enableEvTimeSplitting) {
396398
// Bethe-Bloch TPC distribution - TOF EvTime Splitted
397-
histos.add<TH2>("tracks/evtime/fill/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {8000, 0.0, 1E3}});
398-
histos.add<TH2>("tracks/evtime/tof/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {8000, 0.0, 1E3}});
399-
histos.add<TH2>("tracks/evtime/ft0/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {8000, 0.0, 1E3}});
400-
histos.add<TH2>("tracks/evtime/ft0tof/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {8000, 0.0, 1E3}});
399+
histos.add<TH2>("tracks/evtime/fill/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {dedxAxis}});
400+
histos.add<TH2>("tracks/evtime/tof/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {dedxAxis}});
401+
histos.add<TH2>("tracks/evtime/ft0/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {dedxAxis}});
402+
histos.add<TH2>("tracks/evtime/ft0tof/h2TPCsignVsTPCmomentum", "TPC <-dE/dX> vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TPC <-dE/dx> (a.u.)", HistType::kTH2F, {{1000, -5.f, 5.f}, {dedxAxis}});
401403

402404
// Beta vs pT TOF distribution - TOF EvTime Splitted
403405
histos.add<TH2>("tracks/evtime/fill/h2TOFbetaVsP", "TOF #beta vs #it{p}/Z; Signed #it{p} (GeV/#it{c}); TOF #beta", HistType::kTH2F, {{1000, -5.f, 5.f}, {1200, 0.0, 1.2}});
@@ -695,7 +697,7 @@ struct LFNucleiBATask {
695697
histos.fill(HIST("tracks/hDCAxyVsPt"), track.pt(), track.dcaXY());
696698
histos.fill(HIST("tracks/hDCAzVsPt"), track.pt(), track.dcaZ());
697699

698-
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCcut) {
700+
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCPr) {
699701
if (track.sign() > 0) {
700702
histos.fill(HIST("tracks/proton/dca/hDCAxyVsPtProton"), track.pt(), track.dcaXY());
701703
histos.fill(HIST("tracks/proton/dca/hDCAzVsPtProton"), track.pt(), track.dcaZ());
@@ -704,7 +706,7 @@ struct LFNucleiBATask {
704706
histos.fill(HIST("tracks/proton/dca/hDCAzVsPtantiProton"), track.pt(), track.dcaZ());
705707
}
706708
}
707-
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCcut) {
709+
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCDe) {
708710
if (track.sign() > 0) {
709711
histos.fill(HIST("tracks/deuteron/dca/hDCAxyVsPtDeuteron"), track.pt(), track.dcaXY());
710712
histos.fill(HIST("tracks/deuteron/dca/hDCAzVsPtDeuteron"), track.pt(), track.dcaZ());
@@ -713,7 +715,7 @@ struct LFNucleiBATask {
713715
histos.fill(HIST("tracks/deuteron/dca/hDCAzVsPtantiDeuteron"), track.pt(), track.dcaZ());
714716
}
715717
}
716-
if (std::abs(track.tpcNSigmaTr()) < nsigmaTPCcut) {
718+
if (std::abs(track.tpcNSigmaTr()) < nsigmaTPCTr) {
717719
if (track.sign() > 0) {
718720
histos.fill(HIST("tracks/triton/dca/hDCAxyVsPtTriton"), track.pt(), track.dcaXY());
719721
histos.fill(HIST("tracks/triton/dca/hDCAzVsPtTriton"), track.pt(), track.dcaZ());
@@ -722,7 +724,7 @@ struct LFNucleiBATask {
722724
histos.fill(HIST("tracks/triton/dca/hDCAzVsPtantiTriton"), track.pt(), track.dcaZ());
723725
}
724726
}
725-
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCcut) {
727+
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCHe) {
726728
if (track.sign() > 0) {
727729
histos.fill(HIST("tracks/helium/dca/hDCAxyVsPtHelium"), track.pt(), track.dcaXY());
728730
histos.fill(HIST("tracks/helium/dca/hDCAzVsPtHelium"), track.pt(), track.dcaZ());
@@ -731,7 +733,7 @@ struct LFNucleiBATask {
731733
histos.fill(HIST("tracks/helium/dca/hDCAzVsPtantiHelium"), track.pt(), track.dcaZ());
732734
}
733735
}
734-
if (std::abs(track.tpcNSigmaAl()) < nsigmaTPCcut) {
736+
if (std::abs(track.tpcNSigmaAl()) < nsigmaTPCAl) {
735737
if (track.sign() > 0) {
736738
histos.fill(HIST("tracks/alpha/dca/hDCAxyVsPtAlpha"), track.pt(), track.dcaXY());
737739
histos.fill(HIST("tracks/alpha/dca/hDCAzVsPtAlpha"), track.pt(), track.dcaZ());
@@ -1133,36 +1135,36 @@ struct LFNucleiBATask {
11331135

11341136
// PID
11351137

1136-
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCcut && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Proton))) < yCut) {
1138+
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCPr && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Proton))) < yCut) {
11371139
if (track.sign() > 0) {
11381140
histos.fill(HIST("tracks/proton/h1ProtonSpectra"), track.pt());
11391141
} else {
11401142
histos.fill(HIST("tracks/proton/h1antiProtonSpectra"), track.pt());
11411143
}
11421144
}
11431145

1144-
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaDe()) < nsigmaTPCcut)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= 3))) && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Deuteron))) < yCut) {
1146+
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaDe()) < nsigmaTPCDe)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= nsigmaTPCStrongCut))) && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Deuteron))) < yCut) {
11451147
if (track.sign() > 0) {
11461148
histos.fill(HIST("tracks/deuteron/h1DeuteronSpectra"), track.pt());
11471149
} else {
11481150
histos.fill(HIST("tracks/deuteron/h1antiDeuteronSpectra"), track.pt());
11491151
}
11501152
}
1151-
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaTr()) < nsigmaTPCcut)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= 3))) && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Triton))) < yCut) {
1153+
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaTr()) < nsigmaTPCTr)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= nsigmaTPCStrongCut))) && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Triton))) < yCut) {
11521154
if (track.sign() > 0) {
11531155
histos.fill(HIST("tracks/triton/h1TritonSpectra"), track.pt());
11541156
} else {
11551157
histos.fill(HIST("tracks/triton/h1antiTritonSpectra"), track.pt());
11561158
}
11571159
}
1158-
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaHe()) < nsigmaTPCcut)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= 3))) && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3))) < yCut) {
1160+
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaHe()) < nsigmaTPCHe)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= nsigmaTPCStrongCut))) && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3))) < yCut) {
11591161
if (track.sign() > 0) {
11601162
histos.fill(HIST("tracks/helium/h1HeliumSpectra"), track.pt());
11611163
} else {
11621164
histos.fill(HIST("tracks/helium/h1antiHeliumSpectra"), track.pt());
11631165
}
11641166
}
1165-
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaAl()) < nsigmaTPCcut)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= 3))) && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Alpha))) < yCut) {
1167+
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaAl()) < nsigmaTPCAl)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= nsigmaTPCStrongCut))) && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Alpha))) < yCut) {
11661168
if (track.sign() > 0) {
11671169
histos.fill(HIST("tracks/alpha/h1AlphaSpectra"), track.pt());
11681170
} else {
@@ -1190,21 +1192,21 @@ struct LFNucleiBATask {
11901192
if (track.tpcNSigmaPr() + 3 > ProtonNsigmaTPCcutInTOFbeta) {
11911193
histos.fill(HIST("tracks/h2withTPCProtonCutTOFbetaVsP"), track.p() / (1.f * track.sign()), track.beta());
11921194
}
1193-
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCcut && track.sign() > 0)
1195+
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCPr && track.sign() > 0)
11941196
histos.fill(HIST("tracks/proton/h2ProtonTOFbetaVsP"), track.p(), track.beta());
1195-
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCcut && track.sign() < 0)
1197+
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCPr && track.sign() < 0)
11961198
histos.fill(HIST("tracks/proton/h2antiProtonTOFbetaVsP"), track.p(), track.beta());
1197-
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCcut && track.sign() > 0)
1199+
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCDe && track.sign() > 0)
11981200
histos.fill(HIST("tracks/deuteron/h2DeuteronTOFbetaVsP"), track.p(), track.beta());
1199-
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCcut && track.sign() < 0)
1201+
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCDe && track.sign() < 0)
12001202
histos.fill(HIST("tracks/deuteron/h2antiDeuteronTOFbetaVsP"), track.p(), track.beta());
1201-
if (std::abs(track.tpcNSigmaTr()) < nsigmaTPCcut && track.sign() > 0)
1203+
if (std::abs(track.tpcNSigmaTr()) < nsigmaTPCTr && track.sign() > 0)
12021204
histos.fill(HIST("tracks/triton/h2TritonTOFbetaVsP"), track.p(), track.beta());
1203-
if (std::abs(track.tpcNSigmaTr()) < nsigmaTPCcut && track.sign() < 0)
1205+
if (std::abs(track.tpcNSigmaTr()) < nsigmaTPCTr && track.sign() < 0)
12041206
histos.fill(HIST("tracks/triton/h2antiTritonTOFbetaVsP"), track.p(), track.beta());
1205-
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCcut && track.sign() > 0)
1207+
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCHe && track.sign() > 0)
12061208
histos.fill(HIST("tracks/helium/h2HeliumTOFbetaVsP"), track.p(), track.beta());
1207-
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCcut && track.sign() < 0)
1209+
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCHe && track.sign() < 0)
12081210
histos.fill(HIST("tracks/helium/h2antiHeliumTOFbetaVsP"), track.p(), track.beta());
12091211

12101212
if (track.sign() > 0) {
@@ -1252,7 +1254,7 @@ struct LFNucleiBATask {
12521254
}
12531255
}
12541256

1255-
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCcut && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Proton))) < yCut) {
1257+
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCPr && TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Proton))) < yCut) {
12561258
if (track.sign() > 0) {
12571259
histos.fill(HIST("tracks/proton/h2TOFmass2ProtonVsPt"), massTOF * massTOF - fMassProton * fMassProton, track.pt());
12581260
histos.fill(HIST("tracks/proton/h2ProtonTOFExpSignalDiffVsPtCut"), track.pt(), track.tofExpSignalDiffPr());
@@ -1294,7 +1296,7 @@ struct LFNucleiBATask {
12941296
}
12951297
}
12961298

1297-
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaDe()) < nsigmaTPCcut)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= 3))) && (TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Deuteron))) < yCut)) {
1299+
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaDe()) < nsigmaTPCDe)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= nsigmaTPCStrongCut))) && (TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Deuteron))) < yCut)) {
12981300
if (track.sign() > 0) {
12991301
histos.fill(HIST("tracks/deuteron/h2TOFmass2DeuteronVsPt"), massTOF * massTOF - fMassDeuteron * fMassDeuteron, track.pt());
13001302
histos.fill(HIST("tracks/deuteron/h2DeuteronTOFExpSignalDiffVsPtCut"), track.pt(), track.tofExpSignalDiffDe());
@@ -1336,7 +1338,8 @@ struct LFNucleiBATask {
13361338
}
13371339
}
13381340
}
1339-
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaHe()) < nsigmaTPCcut)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= 3))) && (TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3))) < yCut)) {
1341+
1342+
if ((((!enableStrongCut) && (std::abs(track.tpcNSigmaHe()) < nsigmaTPCHe)) || ((enableStrongCut) && (std::abs(track.tpcNSigmaPr()) >= nsigmaTPCStrongCut))) && (TMath::Abs(track.rapidity(o2::track::PID::getMass2Z(o2::track::PID::Helium3))) < yCut)) {
13401343
if (track.sign() > 0) {
13411344
histos.fill(HIST("tracks/helium/h2TOFmass2HeliumVsPt"), massTOF * massTOF - fMassHelium * fMassHelium, track.pt());
13421345
histos.fill(HIST("tracks/helium/h2HeliumTOFExpSignalDiffVsPtCut"), track.pt(), track.tofExpSignalDiffHe());
@@ -1370,7 +1373,7 @@ struct LFNucleiBATask {
13701373
histos.fill(HIST("tracks/proton/h1ProtonSpectraTrue"), track.pt());
13711374
histos.fill(HIST("tracks/proton/dca/hDCAxyVsPtProtonTrue"), track.pt(), track.dcaXY());
13721375
histos.fill(HIST("tracks/proton/dca/hDCAzVsPtProtonTrue"), track.pt(), track.dcaZ());
1373-
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCcut) {
1376+
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCPr) {
13741377
histos.fill(HIST("tracks/proton/h1ProtonSpectraTrueWPID"), track.pt());
13751378
}
13761379
if (isPhysPrim) {
@@ -1393,7 +1396,7 @@ struct LFNucleiBATask {
13931396
histos.fill(HIST("tracks/proton/h1antiProtonSpectraTrue"), track.pt());
13941397
histos.fill(HIST("tracks/proton/dca/hDCAxyVsPtantiProtonTrue"), track.pt(), track.dcaXY());
13951398
histos.fill(HIST("tracks/proton/dca/hDCAzVsPtantiProtonTrue"), track.pt(), track.dcaZ());
1396-
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCcut) {
1399+
if (std::abs(track.tpcNSigmaPr()) < nsigmaTPCPr) {
13971400
histos.fill(HIST("tracks/proton/h1antiProtonSpectraTrueWPID"), track.pt());
13981401
}
13991402
if (isPhysPrim) {
@@ -1417,7 +1420,7 @@ struct LFNucleiBATask {
14171420
histos.fill(HIST("tracks/deuteron/h1DeuteronSpectraTrue"), track.pt());
14181421
histos.fill(HIST("tracks/deuteron/dca/hDCAxyVsPtDeuteronTrue"), track.pt(), track.dcaXY());
14191422
histos.fill(HIST("tracks/deuteron/dca/hDCAzVsPtDeuteronTrue"), track.pt(), track.dcaZ());
1420-
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCcut) {
1423+
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCDe) {
14211424
histos.fill(HIST("tracks/proton/h1DeuteronSpectraTrueWPID"), track.pt());
14221425
}
14231426
if (isPhysPrim) {
@@ -1439,7 +1442,7 @@ struct LFNucleiBATask {
14391442
histos.fill(HIST("tracks/deuteron/h1antiDeuteronSpectraTrue"), track.pt());
14401443
histos.fill(HIST("tracks/deuteron/dca/hDCAxyVsPtantiDeuteronTrue"), track.pt(), track.dcaXY());
14411444
histos.fill(HIST("tracks/deuteron/dca/hDCAzVsPtantiDeuteronTrue"), track.pt(), track.dcaZ());
1442-
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCcut) {
1445+
if (std::abs(track.tpcNSigmaDe()) < nsigmaTPCDe) {
14431446
histos.fill(HIST("tracks/proton/h1antiDeuteronSpectraTrueWPID"), track.pt());
14441447
}
14451448
if (isPhysPrim) {
@@ -1501,7 +1504,7 @@ struct LFNucleiBATask {
15011504
histos.fill(HIST("tracks/helium/h1HeliumSpectraTrue"), track.pt());
15021505
histos.fill(HIST("tracks/helium/dca/hDCAxyVsPtHeliumTrue"), track.pt(), track.dcaXY());
15031506
histos.fill(HIST("tracks/helium/dca/hDCAzVsPtHeliumTrue"), track.pt(), track.dcaZ());
1504-
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCcut) {
1507+
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCHe) {
15051508
histos.fill(HIST("tracks/proton/h1HeliumSpectraTrueWPID"), track.pt());
15061509
}
15071510
if (isPhysPrim) {
@@ -1524,7 +1527,7 @@ struct LFNucleiBATask {
15241527
histos.fill(HIST("tracks/helium/h1antiHeliumSpectraTrue"), track.pt());
15251528
histos.fill(HIST("tracks/helium/dca/hDCAxyVsPtantiHeliumTrue"), track.pt(), track.dcaXY());
15261529
histos.fill(HIST("tracks/helium/dca/hDCAzVsPtantiHeliumTrue"), track.pt(), track.dcaZ());
1527-
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCcut) {
1530+
if (std::abs(track.tpcNSigmaHe()) < nsigmaTPCHe) {
15281531
histos.fill(HIST("tracks/proton/h1antiHeliumSpectraTrueWPID"), track.pt());
15291532
}
15301533
if (isPhysPrim) {

0 commit comments

Comments
 (0)