@@ -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,13 @@ 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)
911+ fillTrainingTable (coll, casc, pdgCode); // I only store cascades that passed PID and track quality selections
909912 }
910913 }
911914
@@ -967,8 +970,12 @@ struct cascadeFlow {
967970 auto bachExtra = casc.bachTrackExtra_as <DauTracks>();
968971
969972 int counter = 0 ;
970- IsCascAccepted (casc, negExtra, posExtra, bachExtra, counter);
973+ bool isCascCandidate = 0 ;
974+ isCascCandidate = IsCascAccepted (casc, negExtra, posExtra, bachExtra, counter);
971975 histos.fill (HIST (" hCascade" ), counter);
976+ histos.fill (HIST (" hCascadeDauSel" ), (int )isCascCandidate);
977+ if (!isCascCandidate)
978+ continue ;
972979
973980 // ML selections
974981 bool isSelectedCasc[2 ]{false , false };
@@ -1242,8 +1249,12 @@ struct cascadeFlow {
12421249 auto bachExtra = casc.bachTrackExtra_as <DauTracks>();
12431250
12441251 int counter = 0 ;
1245- IsCascAccepted (casc, negExtra, posExtra, bachExtra, counter);
1252+ bool isCascCandidate = 0 ;
1253+ isCascCandidate = IsCascAccepted (casc, negExtra, posExtra, bachExtra, counter);
12461254 histos.fill (HIST (" hCascade" ), counter);
1255+ histos.fill (HIST (" hCascadeDauSel" ), (int )isCascCandidate);
1256+ if (!isCascCandidate)
1257+ continue ;
12471258
12481259 // ML selections
12491260 bool isSelectedCasc[nParticles]{false , false };
@@ -1682,8 +1693,12 @@ struct cascadeFlow {
16821693 auto bachExtra = casc.bachTrackExtra_as <DauTracks>();
16831694
16841695 int counter = 0 ;
1685- IsCascAccepted (casc, negExtra, posExtra, bachExtra, counter);
1696+ bool isCascCandidate = 0 ;
1697+ isCascCandidate = IsCascAccepted (casc, negExtra, posExtra, bachExtra, counter);
16861698 histos.fill (HIST (" hCascade" ), counter);
1699+ histos.fill (HIST (" hCascadeDauSel" ), (int )isCascCandidate);
1700+ if (!isCascCandidate)
1701+ continue ;
16871702
16881703 // ML selections
16891704 bool isSelectedCasc[nParticles]{false , false };
@@ -1828,8 +1843,12 @@ struct cascadeFlow {
18281843 auto bachExtra = casc.bachTrackExtra_as <DauTracks>();
18291844
18301845 int counter = 0 ;
1831- IsCascAccepted (casc, negExtra, posExtra, bachExtra, counter);
1846+ bool isCascCandidate = 0 ;
1847+ isCascCandidate = IsCascAccepted (casc, negExtra, posExtra, bachExtra, counter);
18321848 histos.fill (HIST (" hCascade" ), counter);
1849+ histos.fill (HIST (" hCascadeDauSel" ), (int )isCascCandidate);
1850+ if (!isCascCandidate)
1851+ continue ;
18331852
18341853 // ML selections
18351854 bool isSelectedCasc[nParticles]{false , false };
0 commit comments