[PWGLF] Update on mother's decay length. (#8949)

Author: mtorresc15

PWGLF/Tasks/Nuspex/spectraTOF.cxx

@@ -340,7 +340,7 @@ struct tofSpectra {
     const AxisSpec phiAxis{200, 0, 7, "#it{#varphi} (rad)"};
     const AxisSpec dcaZAxis{binsOptions.binsDca, "DCA_{z} (cm)"};
     const AxisSpec lengthAxis{100, 0, 600, "Track length (cm)"};
-    const AxisSpec decayLengthAxis{100, 0, 0.1, "Decay Length (cm)"};
+    const AxisSpec decayLengthAxis{100, 0, 0.5, "Decay Length (cm)"};
 
     if (enableTrackCutHistograms) {
       const AxisSpec chargeAxis{2, -2.f, 2.f, "Charge"};
@@ -1650,12 +1650,15 @@ struct tofSpectra {
           histos.fill(HIST(hdcaxystr[i]), track.pt(), track.dcaXY());
           histos.fill(HIST(hdcazstr[i]), track.pt(), track.dcaZ());
         }
-        if (mcParticle.has_daughters()) {
-          auto daughter0 = mcParticle.template daughters_as<aod::McParticles>().begin();
-          double vertexDau[3] = {daughter0.vx(), daughter0.vy(), daughter0.vz()};
-          double vertexPrimary[3] = {mcCollision.posX(), mcCollision.posY(), mcCollision.posZ()};
-          auto decayLength = RecoDecay::distance(vertexPrimary, vertexDau) / 10000;
-          hDecayLengthStr[i]->Fill(track.pt(), decayLength);
+
+        if (mcParticle.has_mothers()) {
+          for (const auto& mother : mcParticle.template mothers_as<aod::McParticles>()) {
+            auto daughter0 = mother.template daughters_as<aod::McParticles>().begin();
+            double vertexDau[3] = {daughter0.vx(), daughter0.vy(), daughter0.vz()};
+            double vertexMoth[3] = {mother.vx(), mother.vy(), mother.vz()};
+            auto decayLength = RecoDecay::distance(vertexMoth, vertexDau);
+            hDecayLengthStr[i]->Fill(track.pt(), decayLength);
+          }
         }
       } else {
         if (enableDCAxyzHistograms) {
@@ -1722,23 +1725,25 @@ struct tofSpectra {
           hDcaZMCNotHF[i]->Fill(track.pt(), track.dcaZ());
         }
 
-        if (mcParticle.has_daughters()) {
-          auto daughter0 = mcParticle.template daughters_as<aod::McParticles>().begin();
-          double vertexDau[3] = {daughter0.vx(), daughter0.vy(), daughter0.vz()};
-          double vertexPrimary[3] = {mcCollision.posX(), mcCollision.posY(), mcCollision.posZ()};
-          auto decayLength = RecoDecay::distance(vertexPrimary, vertexDau) / 10000;
+        if (mcParticle.has_mothers()) {
+          for (const auto& mother : mcParticle.template mothers_as<aod::McParticles>()) {
+            auto daughter0 = mother.template daughters_as<aod::McParticles>().begin();
+            double vertexDau[3] = {daughter0.vx(), daughter0.vy(), daughter0.vz()};
+            double vertexMoth[3] = {mother.vx(), mother.vy(), mother.vz()};
+            auto decayLength = RecoDecay::distance(vertexMoth, vertexDau);
 
-          if (IsD0Mother) {
-            hDecayLengthMCD0[i]->Fill(track.pt(), decayLength);
-          }
-          if (IsCharmMother) {
-            hDecayLengthMCCharm[i]->Fill(track.pt(), decayLength);
-          }
-          if (IsBeautyMother) {
-            hDecayLengthMCBeauty[i]->Fill(track.pt(), decayLength);
-          }
-          if (IsNotHFMother) {
-            hDecayLengthMCNotHF[i]->Fill(track.pt(), decayLength);
+            if (IsD0Mother) {
+              hDecayLengthMCD0[i]->Fill(track.pt(), decayLength);
+            }
+            if (IsCharmMother) {
+              hDecayLengthMCCharm[i]->Fill(track.pt(), decayLength);
+            }
+            if (IsBeautyMother) {
+              hDecayLengthMCBeauty[i]->Fill(track.pt(), decayLength);
+            }
+            if (IsNotHFMother) {
+              hDecayLengthMCNotHF[i]->Fill(track.pt(), decayLength);
+            }
           }
         }
       }