diff --git a/PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx b/PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx
index b9f32b8600c..a1947c8b6b5 100644
--- a/PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx
+++ b/PWGLF/Tasks/Strangeness/cascadeanalysisMC.cxx
@@ -80,6 +80,7 @@ struct cascadeGenerated {
   Configurable<int> nPtBins{"nPtBins", 200, "number of pT bins"};
   Configurable<float> rapidityCut{"rapidityCut", 0.5, "max (absolute) rapidity of generated cascade"};
   Configurable<int> nRapidityBins{"nRapidityBins", 200, "number of pT bins"};
+  Configurable<bool> requirePhysicalPrimary{"requirePhysicalPrimary", false, "require the generated cascade to be a physical primary"};
 
   void init(InitContext const&)
   {
@@ -108,6 +109,8 @@ struct cascadeGenerated {
     for (auto& particle : mcparts) {
       if (TMath::Abs(particle.y()) > rapidityCut)
         continue;
+      if (requirePhysicalPrimary && !particle.isPhysicalPrimary())
+        continue;
       if (particle.pdgCode() == 3312) {
         registry.fill(HIST("hPtXiMinus"), particle.pt());
         registry.fill(HIST("h2DXiMinus"), particle.pt(), particle.y());