Merge

Generators/FaserParticleGun/python/RadialPosSampler.py

@@ -89,3 +89,4 @@ if __name__ == "__main__":
     plt.hist(np.sqrt(xarr**2 + yarr**2))    
     plt.tight_layout()
     plt.show()
+

Generators/ForeseeGenerator/python/Validate.py

@@ -163,7 +163,8 @@ class EvgenValidation(EvgenAnalysisAlg):
         llp_vtx = None
         
         for i, p in enumerate(evt.particles):
-            #p.print()
+            print("--- ", i)
+            p.print()
             self.fillDaughter(p)
 
             if self.mother_stored:
@@ -184,6 +185,9 @@ class EvgenValidation(EvgenAnalysisAlg):
             if p.production_vertex():
                 vertices.append(p.production_vertex().point3d())
 
+            if p.production_vertex():
+                vertices.append(p.production_vertex().point3d())
+
         # Fill daughter plots
         for i in range(self.ndaughters):
             if i >= len(momenta): continue