Skip to content

Commit ef07197

Browse files
author
Chiara De Martin
committed
Fix daughter track selection
1 parent 1c2dc8b commit ef07197

File tree

1 file changed

+21
-6
lines changed

1 file changed

+21
-6
lines changed

PWGLF/TableProducer/Strangeness/cascadeflow.cxx

Lines changed: 21 additions & 6 deletions
Original file line numberDiff line numberDiff line change
@@ -662,6 +662,7 @@ struct cascadeFlow {
662662
histos.add("hLambdaCandidate", "hLambdaCandidate", HistType::kTH1F, {{5, -0.5, 4.5}});
663663
histos.add("hCascadeSignal", "hCascadeSignal", HistType::kTH1F, {{6, -0.5, 5.5}});
664664
histos.add("hCascade", "hCascade", HistType::kTH1F, {{6, -0.5, 5.5}});
665+
histos.add("hCascadeDauSel", "hCascadeDauSel", HistType::kTH1F, {{2, -0.5, 1.5}});
665666
histos.add("hLambdaDauSel", "hLambdaDauSel", HistType::kTH1F, {{3, -0.5, 2.5}});
666667
histos.add("hALambdaDauSel", "hALambdaDauSel", HistType::kTH1F, {{3, -0.5, 2.5}});
667668
histos.add("hXiPtvsCent", "hXiPtvsCent", HistType::kTH2F, {{100, 0, 100}, {400, 0, 20}});
@@ -901,11 +902,12 @@ struct cascadeFlow {
901902
auto bachExtra = casc.bachTrackExtra_as<DauTracks>();
902903

903904
int counter = 0;
904-
IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
905+
bool isCascCandidate = 0;
906+
isCascCandidate = IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
905907
histos.fill(HIST("hCascadeSignal"), counter);
906908

907909
// PDG cascades
908-
fillTrainingTable(coll, casc, pdgCode);
910+
if (isCascCandidate) fillTrainingTable(coll, casc, pdgCode); //I only store cascades that passed PID and track quality selections
909911
}
910912
}
911913

@@ -967,8 +969,11 @@ struct cascadeFlow {
967969
auto bachExtra = casc.bachTrackExtra_as<DauTracks>();
968970

969971
int counter = 0;
970-
IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
972+
bool isCascCandidate = 0;
973+
isCascCandidate = IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
971974
histos.fill(HIST("hCascade"), counter);
975+
histos.fill(HIST("hCascadeDauSel"), (int)isCascCandidate);
976+
if (!isCascCandidate) continue;
972977

973978
// ML selections
974979
bool isSelectedCasc[2]{false, false};
@@ -1242,8 +1247,11 @@ struct cascadeFlow {
12421247
auto bachExtra = casc.bachTrackExtra_as<DauTracks>();
12431248

12441249
int counter = 0;
1245-
IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
1250+
bool isCascCandidate = 0;
1251+
isCascCandidate = IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
12461252
histos.fill(HIST("hCascade"), counter);
1253+
histos.fill(HIST("hCascadeDauSel"), (int)isCascCandidate);
1254+
if (!isCascCandidate) continue;
12471255

12481256
// ML selections
12491257
bool isSelectedCasc[nParticles]{false, false};
@@ -1682,8 +1690,12 @@ struct cascadeFlow {
16821690
auto bachExtra = casc.bachTrackExtra_as<DauTracks>();
16831691

16841692
int counter = 0;
1685-
IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
1693+
bool isCascCandidate = 0;
1694+
isCascCandidate = IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
16861695
histos.fill(HIST("hCascade"), counter);
1696+
histos.fill(HIST("hCascadeDauSel"), (int)isCascCandidate);
1697+
if (!isCascCandidate) continue;
1698+
16871699

16881700
// ML selections
16891701
bool isSelectedCasc[nParticles]{false, false};
@@ -1828,8 +1840,11 @@ struct cascadeFlow {
18281840
auto bachExtra = casc.bachTrackExtra_as<DauTracks>();
18291841

18301842
int counter = 0;
1831-
IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
1843+
bool isCascCandidate = 0;
1844+
isCascCandidate = IsCascAccepted(casc, negExtra, posExtra, bachExtra, counter);
18321845
histos.fill(HIST("hCascade"), counter);
1846+
histos.fill(HIST("hCascadeDauSel"), (int)isCascCandidate);
1847+
if (!isCascCandidate) continue;
18331848

18341849
// ML selections
18351850
bool isSelectedCasc[nParticles]{false, false};

0 commit comments

Comments
 (0)