@@ -79,17 +79,17 @@ struct nucleiFromHypertritonMap {
7979
8080 void init (InitContext const &)
8181 {
82- registryMC.add (" hypertritonPtgen " , " hypertritonPtGen" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
82+ registryMC.add (" hypertritonPtGen " , " hypertritonPtGen" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
8383 if (saveHelium) {
8484 registryMC.add (" he3SecPtRec_from_hypertriton" , " he3SecPtRec_from_hypertriton" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
85- registryMC.add (" hyperHe4Ptgen " , " hyperHe4PtGen" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
85+ registryMC.add (" hyperHe4PtGen " , " hyperHe4PtGen" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
8686 registryMC.add (" he3SecPtRec_from_hyperHe4" , " he3SecPtRec_from_hyperHe4" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
8787 registryMC.add (" he3PtRec" , " he3PtRec" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
8888 registryMC.add (" he3PtGen" , " he3PtGen" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
8989 } else {
9090 registryMC.add (" deutSecPtRec_from_hypertriton" , " deutSecPtRec_from_hypertriton" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
9191 registryMC.add (" deutPtRec" , " deutPtRec" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
92- registryMC.add (" deutPtgen " , " deutPtGen" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
92+ registryMC.add (" deutPtGen " , " deutPtGen" , HistType::kTH1F , {{nbin_pt, min_pt, max_pt, " p_{T} (GeV/c)" }});
9393 }
9494 }
9595
@@ -107,21 +107,21 @@ struct nucleiFromHypertritonMap {
107107 if (mcparticle.pdgCode () == AntiHypertritonPDG) {
108108 for (auto & daughter : mcparticle.daughters_as <aod::McParticles>()) {
109109 if (daughter.pdgCode () == selectedPDG) {
110- registryMC.fill (HIST (" hypertritonPtgen " ), mcparticle.pt ());
110+ registryMC.fill (HIST (" hypertritonPtGen " ), mcparticle.pt ());
111111 }
112112 }
113113 }
114114 if (mcparticle.pdgCode () == AntiHyperHelium4PDG) {
115115 for (auto & daughter : mcparticle.daughters_as <aod::McParticles>()) {
116116 if (daughter.pdgCode () == selectedPDG) {
117- registryMC.fill (HIST (" hyperHe4Ptgen " ), mcparticle.pt ());
117+ registryMC.fill (HIST (" hyperHe4PtGen " ), mcparticle.pt ());
118118 }
119119 }
120120 }
121- if (mcparticle.pdgCode () == AntihePDG) {
121+ if (mcparticle.pdgCode () == AntihePDG && mcparticle. isPhysicalPrimary () ) {
122122 registryMC.fill (HIST (" he3PtGen" ), mcparticle.pt ());
123123 }
124- if (mcparticle.pdgCode () == AntideuteronPDG) {
124+ if (mcparticle.pdgCode () == AntideuteronPDG && mcparticle. isPhysicalPrimary () ) {
125125 registryMC.fill (HIST (" deutPtGen" ), mcparticle.pt ());
126126 }
127127 }
@@ -148,10 +148,10 @@ struct nucleiFromHypertritonMap {
148148 track.itsChi2NCl () > 36 .f ) {
149149 continue ;
150150 }
151- if (mcparticle.pdgCode () == AntideuteronPDG) {
151+ if (mcparticle.pdgCode () == AntideuteronPDG && mcparticle. isPhysicalPrimary () ) {
152152 registryMC.fill (HIST (" deutPtRec" ), track.pt ());
153153 }
154- if (mcparticle.pdgCode () == AntihePDG) {
154+ if (mcparticle.pdgCode () == AntihePDG && mcparticle. isPhysicalPrimary () ) {
155155 registryMC.fill (HIST (" he3PtRec" ), 2 * track.pt ());
156156 }
157157
0 commit comments