Add histograms for generated daughters (#5572)

Co-authored-by: Andrea Barra <andrea.barra@to.infn.it>

Author: brr-andrea

PWGLF/Tasks/Resonances/deltaanalysis.cxx

@@ -159,6 +159,15 @@ struct deltaAnalysis {
       histos.add("hRecPionDeltaPlusPlus", "Pion from #Delta^{++}", kTH1F, {ptAxis});
       histos.add("hRecPionAntiDeltaZero", "Pion from #bar{#Delta^{0}}", kTH1F, {ptAxis});
       histos.add("hRecPionDeltaZero", "Pion from #Delta^{0}", kTH1F, {ptAxis});
+
+      histos.add("hGenProtonAntiDeltaPlusPlus", "Proton from #bar{#Delta^{++}}", kTH1F, {ptAxis});
+      histos.add("hGenProtonDeltaPlusPlus", "Proton from #Delta^{++}", kTH1F, {ptAxis});
+      histos.add("hGenProtonAntiDeltaZero", "Proton from #bar{#Delta^{0}}", kTH1F, {ptAxis});
+      histos.add("hGenProtonDeltaZero", "Proton from #Delta^{0}", kTH1F, {ptAxis});
+      histos.add("hGenPionAntiDeltaPlusPlus", "Pion from #bar{#Delta^{++}}", kTH1F, {ptAxis});
+      histos.add("hGenPionDeltaPlusPlus", "Pion from #Delta^{++}", kTH1F, {ptAxis});
+      histos.add("hGenPionAntiDeltaZero", "Pion from #bar{#Delta^{0}}", kTH1F, {ptAxis});
+      histos.add("hGenPionDeltaZero", "Pion from #Delta^{0}", kTH1F, {ptAxis});
     }
   }
 
@@ -564,13 +573,17 @@ struct deltaAnalysis {
       auto daughters = mcParticle.daughters_as<aod::McParticles>();
       bool daughtPr = false;
       bool daughtPi = false;
+      float ptPr = -999.;
+      float ptPi = -999.;
       double eDelta = 0.;
       for (auto daught : daughters) {
         if (std::abs(daught.pdgCode()) == protonPDG) {
           daughtPr = true;
+          ptPr = daught.pt();
           eDelta += daught.e();
         } else if (std::abs(daught.pdgCode()) == pionPDG) {
           daughtPi = true;
+          ptPi = daught.pt();
           eDelta += daught.e();
         }
       }
@@ -584,14 +597,22 @@ struct deltaAnalysis {
       if (pdg > 0) {
         if (std::abs(pdg) == deltaPlusPlusPDG) {
           histos.fill(HIST("hDeltaPlusPlusInvMassGen"), mcParticle.pt(), gen_mass);
+          histos.fill(HIST("hGenProtonDeltaPlusPlus"), ptPr);
+          histos.fill(HIST("hGenPionDeltaPlusPlus"), ptPi);
         } else if (std::abs(pdg) == deltaZeroPDG) {
           histos.fill(HIST("hDeltaZeroInvMassGen"), mcParticle.pt(), gen_mass);
+          histos.fill(HIST("hGenProtonDeltaZero"), ptPr);
+          histos.fill(HIST("hGenPionDeltaZero"), ptPi);
         }
       } else {
         if (std::abs(pdg) == deltaPlusPlusPDG) {
           histos.fill(HIST("hAntiDeltaPlusPlusInvMassGen"), mcParticle.pt(), gen_mass);
+          histos.fill(HIST("hGenProtonAntiDeltaPlusPlus"), ptPr);
+          histos.fill(HIST("hGenPionAntiDeltaPlusPlus"), ptPi);
         } else if (std::abs(pdg) == deltaZeroPDG) {
           histos.fill(HIST("hAntiDeltaZeroInvMassGen"), mcParticle.pt(), gen_mass);
+          histos.fill(HIST("hGenProtonAntiDeltaZero"), ptPr);
+          histos.fill(HIST("hGenPionAntiDeltaZero"), ptPi);
         }
       }
     }