Skip to content

Commit 1c1c116

Browse files
giovannimalfattoregiovannimalfattore
committed
NucleiTask - Add mothers - Small fix
1 parent bc3cd60 commit 1c1c116

File tree

1 file changed

+28
-27
lines changed

1 file changed

+28
-27
lines changed

‎PWGLF/Tasks/Nuspex/LFNucleiBATask.cxx‎

Lines changed: 28 additions & 27 deletions
Original file line numberDiff line numberDiff line change
@@ -217,7 +217,8 @@ struct LFNucleiBATask {
217217
// ITS to TPC - Fake hit loop
218218
static constexpr int kFakeLoop = 10; // Fixed O2Linter error
219219
// TPC low/high momentum range
220-
static constexpr float cfgTpcClasses[] = {0.5f, 0.1f};
220+
static constexpr float kCfgTpcClasses[] = {0.5f, 0.1f};
221+
static constexpr float kCfgKaonCut = 5.f;
221222

222223
// Weak-decay flag
223224
static constexpr int kProcessWeakDecay = 4;
@@ -974,20 +975,20 @@ struct LFNucleiBATask {
974975
histos.add<TH1>("tracks/deuteron/dca/before/hNumMothers", "N mothers per particle; N mothers;counts", HistType::kTH1I, {{7, 1.0, 8.0}});
975976
histos.add<TH3>("tracks/deuteron/dca/before/hMomTrueMaterial", "MC mothers;mother index;mother type; mother #it{p}_{T}", HistType::kTH3F, {{kMaxNumMom + 2, -0.5, static_cast<double>(kMaxNumMom) + 1.5}, {kNumMotherlist + 2, -1.5, static_cast<double>(kNumMotherlist) + 0.5}, {250, 0.0, 10.0}});
976977

977-
std::shared_ptr<TH3> hTemp_De = histos.get<TH3>(HIST("tracks/deuteron/dca/before/hMomTrueMaterial"));
978-
TH3* hPDG_De = hTemp_De.get();
978+
std::shared_ptr<TH3> hTempDe = histos.get<TH3>(HIST("tracks/deuteron/dca/before/hMomTrueMaterial"));
979+
TH3* hPDG_De = hTempDe.get();
979980

980-
TAxis* axPDG_De = hPDG_De->GetXaxis();
981+
TAxis* axPdgDe = hPDG_De->GetXaxis();
981982
for (int i = 0; i <= kMaxNumMom; i++) {
982-
axPDG_De->SetBinLabel(i + 1, Form("%d", i));
983+
axPdgDe->SetBinLabel(i + 1, Form("%d", i));
983984
}
984-
axPDG_De->SetBinLabel(kMaxNumMom + 2, ">=5");
985+
axPdgDe->SetBinLabel(kMaxNumMom + 2, ">=5");
985986

986-
TAxis* ayPDG_De = hPDG_De->GetYaxis();
987-
ayPDG_De->SetBinLabel(1, "undef.");
988-
ayPDG_De->SetBinLabel(2, "other");
987+
TAxis* ayPdgDe = hPDG_De->GetYaxis();
988+
ayPdgDe->SetBinLabel(1, "undef.");
989+
ayPdgDe->SetBinLabel(2, "other");
989990
for (int i = 0; i < kNumMotherlist; i++) {
990-
ayPDG_De->SetBinLabel(i + 3, kMotherNames[i]);
991+
ayPdgDe->SetBinLabel(i + 3, kMotherNames[i]);
991992
}
992993

993994
histos.add<TH2>("tracks/deuteron/dca/before/hDCAxyVsPtDeuteronTrueTransport", "DCAxy vs Pt (d); #it{p}_{T} (GeV/#it{c}); DCAxy (cm)", HistType::kTH2F, {{ptAxis}, {dcaxyAxis}});
@@ -1191,20 +1192,20 @@ struct LFNucleiBATask {
11911192
histos.add<TH3>("tracks/helium/dca/before/hMomTrueMaterial", "MC mothers;mother index;mother type; mother #it{p}_{T}", HistType::kTH3F, {{kMaxNumMom + 2, -0.5, static_cast<double>(kMaxNumMom) + 1.5}, {kNumMotherlist + 2, -1.5, static_cast<double>(kNumMotherlist) + 0.5}, {250, 0.0, 10.0}});
11921193

11931194
// Fix for getting TH3 pointer
1194-
std::shared_ptr<TH3> hTemp_He = histos.get<TH3>(HIST("tracks/helium/dca/before/hMomTrueMaterial"));
1195-
TH3* hPDG_He = hTemp_He.get();
1195+
std::shared_ptr<TH3> hTempHe = histos.get<TH3>(HIST("tracks/helium/dca/before/hMomTrueMaterial"));
1196+
TH3* hPDG_He = hTempHe.get();
11961197

1197-
TAxis* axPDG_He = hPDG_He->GetXaxis();
1198+
TAxis* axPdgHe = hPDG_He->GetXaxis();
11981199
for (int i = 0; i <= kMaxNumMom; i++) {
1199-
axPDG_He->SetBinLabel(i + 1, Form("%d", i));
1200+
axPdgHe->SetBinLabel(i + 1, Form("%d", i));
12001201
}
1201-
axPDG_He->SetBinLabel(kMaxNumMom + 2, ">=5");
1202+
axPdgHe->SetBinLabel(kMaxNumMom + 2, ">=5");
12021203

1203-
TAxis* ayPDG_He = hPDG_He->GetYaxis();
1204-
ayPDG_He->SetBinLabel(1, "undef.");
1205-
ayPDG_He->SetBinLabel(2, "other");
1204+
TAxis* ayPdgHe = hPDG_He->GetYaxis();
1205+
ayPdgHe->SetBinLabel(1, "undef.");
1206+
ayPdgHe->SetBinLabel(2, "other");
12061207
for (int i = 0; i < kNumMotherlist; i++) {
1207-
ayPDG_He->SetBinLabel(i + 3, kMotherNames[i]);
1208+
ayPdgHe->SetBinLabel(i + 3, kMotherNames[i]);
12081209
}
12091210

12101211
histos.add<TH2>("tracks/helium/dca/before/hDCAxyVsPtHeliumTrueTransport", "DCAxy vs Pt (He); #it{p}_{T} (GeV/#it{c}); DCAxy (cm)", HistType::kTH2F, {{ptZHeAxis}, {dcaxyAxis}});
@@ -3976,26 +3977,26 @@ struct LFNucleiBATask {
39763977
debugHistos.fill(HIST("debug/qa/h2TPCncrVsPtPos"), track.tpcInnerParam(), track.tpcNClsCrossedRows());
39773978
debugHistos.fill(HIST("debug/qa/h2TPCncrVsTPCsignalPos"), track.tpcSignal(), track.tpcNClsCrossedRows());
39783979

3979-
if (track.tpcInnerParam() < cfgTpcClasses[0]) {
3980+
if (track.tpcInnerParam() < kCfgTpcClasses[0]) {
39803981
debugHistos.fill(HIST("debug/qa/h1TPCncrLowPPos"), track.tpcNClsCrossedRows());
39813982
}
3982-
if ((track.tpcInnerParam() >= cfgTpcClasses[0]) && (track.tpcInnerParam() < cfgTpcClasses[1])) {
3983+
if ((track.tpcInnerParam() >= kCfgTpcClasses[0]) && (track.tpcInnerParam() < kCfgTpcClasses[1])) {
39833984
debugHistos.fill(HIST("debug/qa/h1TPCncrMidPPos"), track.tpcNClsCrossedRows());
39843985
}
3985-
if (track.tpcInnerParam() >= cfgTpcClasses[1]) {
3986+
if (track.tpcInnerParam() >= kCfgTpcClasses[1]) {
39863987
debugHistos.fill(HIST("debug/qa/h1TPCncrHighPPos"), track.tpcNClsCrossedRows());
39873988
}
39883989
} else {
39893990
debugHistos.fill(HIST("debug/qa/h2TPCncrVsPtNeg"), track.tpcInnerParam(), track.tpcNClsCrossedRows());
39903991
debugHistos.fill(HIST("debug/qa/h2TPCncrVsTPCsignalNeg"), track.tpcSignal(), track.tpcNClsCrossedRows());
39913992

3992-
if (track.tpcInnerParam() < cfgTpcClasses[0]) {
3993+
if (track.tpcInnerParam() < kCfgTpcClasses[0]) {
39933994
debugHistos.fill(HIST("debug/qa/h1TPCncrLowPNeg"), track.tpcNClsCrossedRows());
39943995
}
3995-
if ((track.tpcInnerParam() >= cfgTpcClasses[0]) && (track.tpcInnerParam() < cfgTpcClasses[1])) {
3996+
if ((track.tpcInnerParam() >= kCfgTpcClasses[0]) && (track.tpcInnerParam() < kCfgTpcClasses[1])) {
39963997
debugHistos.fill(HIST("debug/qa/h1TPCncrMidPNeg"), track.tpcNClsCrossedRows());
39973998
}
3998-
if (track.tpcInnerParam() >= cfgTpcClasses[1]) {
3999+
if (track.tpcInnerParam() >= kCfgTpcClasses[1]) {
39994000
debugHistos.fill(HIST("debug/qa/h1TPCncrHighPNeg"), track.tpcNClsCrossedRows());
40004001
}
40014002
}
@@ -4008,7 +4009,7 @@ struct LFNucleiBATask {
40084009
histos.fill(HIST("tracks/eff/h2pVsTPCmomentum"), track.tpcInnerParam(), track.p());
40094010

40104011
if (filterOptions.enableFiltering) {
4011-
if (track.tpcNSigmaKa() < 5)
4012+
if (track.tpcNSigmaKa() < kCfgKaonCut)
40124013
continue;
40134014
}
40144015

@@ -6398,7 +6399,7 @@ struct LFNucleiBATask {
63986399
for (const auto& mcPart : mcParticles) {
63996400
if (!mcPart.isPhysicalPrimary())
64006401
continue;
6401-
if (std::abs(mcPart.y()) >= 0.5)
6402+
if (std::abs(mcPart.y()) >= kCfgTpcClasses[0])
64026403
continue;
64036404
if (mcPart.pdgCode() == PDGDeuteron) {
64046405
evLossHistos.fill(HIST("evLoss/pt/hDeuteronGen"), mcPart.pt());

0 commit comments

Comments
 (0)