[PWGLF] Corrected error in accessing McCollision indices

PWGLF/TableProducer/Strangeness/sigmaminustask.cxx

@@ -114,7 +114,7 @@
   }
   PROCESS_SWITCH(sigmaminustask, processData, "Data processing", true);
 
-  void processMC(CollisionsFullMC const& collisions, aod::KinkCands const& KinkCands, aod::McTrackLabels const& trackLabelsMC, aod::McParticles const& particlesMC, TracksFull const&)
+  void processMC(CollisionsFullMC const& collisions, aod::KinkCands const& KinkCands, aod::McTrackLabels const& trackLabelsMC, aod::McParticles const& particlesMC, aod::McCollisions const&, TracksFull const&)
   {
     for (const auto& collision : collisions) {
       if (std::abs(collision.posZ()) > cutzvertex || !collision.sel8()) {
@@ -149,26 +149,28 @@
           if (!mcTrackPiDau.has_mothers()) {
             continue;
           }
           for (auto& piMother : mcTrackPiDau.mothers_as<aod::McParticles>()) {
             if (piMother.globalIndex() != mcTrackSigma.globalIndex()) {
               continue;
             }
             if (std::abs(mcTrackSigma.pdgCode()) != 3112) {
               continue;
             }
             if (std::abs(mcTrackPiDau.pdgCode()) != 211 && std::abs(mcTrackPiDau.pdgCode()) != 2212) {
               continue;
             }
 
             float MotherMassMC = std::sqrt(piMother.e() * piMother.e() - piMother.p() * piMother.p());
             float MotherpTMC = piMother.pt();
             float deltaXMother = mcTrackPiDau.vx() - piMother.vx();
             float deltaYMother = mcTrackPiDau.vy() - piMother.vy();
             float decayRadiusMC = std::sqrt(deltaXMother * deltaXMother + deltaYMother * deltaYMother);
 
             // Check coherence of MCcollision Id for daughter MCparticle and reconstructed collision
-            auto mcCollision = mcTrackPiDau.template mcCollision_as<aod::McCollisions>();
-            bool mcCollisionIdCheck = collision.mcCollisionId() == mcCollision.globalIndex();
+            bool mcCollisionIdCheck = false;
+            if (collision.has_mcCollision()) {
+              mcCollisionIdCheck = collision.mcCollision().globalIndex() == mcTrackPiDau.mcCollisionId();
+            }
 
             rSigmaMinus.fill(HIST("h2MassPtMCRec"), kinkCand.mothSign() * kinkCand.ptMoth(), kinkCand.mSigmaMinus());
             if (mcTrackSigma.pdgCode() > 0) {
@@ -198,7 +200,7 @@
 
     // Loop over all generated particles to fill MC histograms
     for (const auto& mcPart : particlesMC) {
       if (std::abs(mcPart.pdgCode()) != 3112 || std::abs(mcPart.y()) > cutEtaMotherMC) { // only sigma mothers and rapidity cut
         continue;
       }
       if (!mcPart.has_daughters()) {
@@ -209,7 +211,7 @@
       std::array<float, 3> secVtx;
       std::array<float, 3> momDaug;
       for (const auto& daughter : mcPart.daughters_as<aod::McParticles>()) {
         if (std::abs(daughter.pdgCode()) == 211 || std::abs(daughter.pdgCode()) == 2212) { // Pi or proton daughter
           hasSigmaDaughter = true;
           secVtx = {daughter.vx(), daughter.vy(), daughter.vz()};
           momDaug = {daughter.px(), daughter.py(), daughter.pz()};