diff --git a/analyzers/dataframe/MCParticle.cc b/analyzers/dataframe/MCParticle.cc
index adf88d713a5..08598180227 100644
--- a/analyzers/dataframe/MCParticle.cc
+++ b/analyzers/dataframe/MCParticle.cc
@@ -426,7 +426,6 @@ ROOT::VecOps::RVec<int> MCParticle::get_parentid(ROOT::VecOps::RVec<int> mcind,
   return result;
 }
 
-
 // ----------------------------------------------------------------------------------------------------------------------------------
 
 // returns one MCParticle selected by its index in the particle block
@@ -676,9 +675,20 @@ ROOT::VecOps::RVec<float> MCParticle::AngleBetweenTwoMCParticles( ROOT::VecOps::
   }
 
   return result;
+}
 
+ROOT::VecOps::RVec<edm4hep::MCParticleData> MCParticle::get_subMC(ROOT::VecOps::RVec<int> index, ROOT::VecOps::RVec<edm4hep::MCParticleData> in){
+  ROOT::VecOps::RVec<edm4hep::MCParticleData> result;
+  for (size_t i = 0; i < index.size(); ++i) {
+    if (index[i]>-1)
+      result.push_back(in.at(index[i]));
+    //else
+    //  std::cout << "electron index negative " << index[i]<<std::endl;
+  }  
+  return result;
 }
 
 
 
 
+
diff --git a/analyzers/dataframe/MCParticle.h b/analyzers/dataframe/MCParticle.h
index fd4059d2e7f..a813eae9d5b 100644
--- a/analyzers/dataframe/MCParticle.h
+++ b/analyzers/dataframe/MCParticle.h
@@ -49,6 +49,8 @@ namespace MCParticle{
     bool m_chargeconjugate = true;
     ROOT::VecOps::RVec<edm4hep::MCParticleData>  operator() (ROOT::VecOps::RVec<edm4hep::MCParticleData> in);
   };
+
+
   
   /// get MC history tree for a given MCParticle index
   struct get_tree{
@@ -91,7 +93,6 @@ namespace MCParticle{
                                         			     ROOT::VecOps::RVec<edm4hep::MCParticleData> in , 
 								     ROOT::VecOps::RVec<int> ind);
 
-
   /// return the parent index of a given list of MC particles
   ROOT::VecOps::RVec<int> get_parentid(ROOT::VecOps::RVec<int> mcind, ROOT::VecOps::RVec<edm4hep::MCParticleData> mc, ROOT::VecOps::RVec<int> parents);
 
@@ -200,7 +201,8 @@ namespace MCParticle{
                                                  ROOT::VecOps::RVec<edm4hep::MCParticleData> in,
                                                  ROOT::VecOps::RVec<int> ind) ;
 
-
+/// return a sub list of MCParticles from a given list of indices
+ROOT::VecOps::RVec<edm4hep::MCParticleData> get_subMC(ROOT::VecOps::RVec<int> index, ROOT::VecOps::RVec<edm4hep::MCParticleData> in);
 
 
 }