Fix: Missing alpha propagation

Author: Henrik Fribert

ALICE3/Core/DelphesO2TrackSmearer.cxx

@@ -106,10 +106,15 @@ bool TrackSmearer::loadTable(int pdg, const char* filename, bool forceReload)
     return false;
   }
   if (mLUTHeader[ipdg]->pdg != pdg) {
-    LOG(info) << " --- LUT header PDG mismatch: expected/detected = " << pdg << "/" << mLUTHeader[ipdg]->pdg << std::endl;
-    delete mLUTHeader[ipdg];
-    mLUTHeader[ipdg] = nullptr;
-    return false;
+    // Special case: Allow Alpha particles to use He3 LUT
+    if (pdg == o2::constants::physics::kAlpha && mLUTHeader[ipdg]->pdg == o2::constants::physics::kHelium3) {
+      LOG(info) << " --- Alpha particles (PDG " << pdg << ") will use He3 LUT data (PDG " << mLUTHeader[ipdg]->pdg << ")" << std::endl;
+    } else {
+      LOG(info) << " --- LUT header PDG mismatch: expected/detected = " << pdg << "/" << mLUTHeader[ipdg]->pdg << std::endl;
+      delete mLUTHeader[ipdg];
+      mLUTHeader[ipdg] = nullptr;
+      return false;
+    }
   }
   const int nnch = mLUTHeader[ipdg]->nchmap.nbins;
   const int nrad = mLUTHeader[ipdg]->radmap.nbins;

ALICE3/Core/DelphesO2TrackSmearer.h

@@ -214,6 +214,8 @@ class TrackSmearer
         return 6; // Triton
       case 1000020030:
         return 7; // Helium3
+      case 1000020040:
+        return 8; // Alphas
       default:
         return 2; // Default: pion
     }
@@ -238,6 +240,8 @@ class TrackSmearer
         return "triton";
       case 1000020030:
         return "helium3";
+      case 1000020040:
+        return "alpha";
       default:
         return "pion"; // Default: pion
     }
@@ -246,7 +250,7 @@ class TrackSmearer
   void setCcdbManager(o2::ccdb::BasicCCDBManager* mgr) { mCcdbManager = mgr; } //;
 
  protected:
-  static constexpr unsigned int nLUTs = 8; // Number of LUT available
+  static constexpr unsigned int nLUTs = 9; // Number of LUT available
   lutHeader_t* mLUTHeader[nLUTs] = {nullptr};
   lutEntry_t***** mLUTEntry[nLUTs] = {nullptr};
   bool mUseEfficiency = true;

ALICE3/TableProducer/OTF/onTheFlyTracker.cxx

@@ -266,6 +266,7 @@ struct OnTheFlyTracker {
       loadLUT(1000010020, lutDe.value);
       loadLUT(1000010030, lutTr.value);
       loadLUT(1000020030, lutHe3.value);
+      loadLUT(1000020040, lutAl.value);
 
       // interpolate efficiencies if requested to do so
       mSmearer.interpolateEfficiency(static_cast<bool>(interpolateLutEfficiencyVsNch));