ENH: Cleaner output for pctweplfit

Author: acoussat

applications/pctweplfit/pctweplfit.py

@@ -106,7 +106,7 @@ def add_detector(name, translation, attach_to_phantom=False):
         phase_space = sim.add_actor("PhaseSpaceActor", "PhaseSpace" + name)
         phase_space.attached_to = plane.name
         phase_space.output_filename = (
-            f"{output}/l{int(phantom_length_mm)}_ps{name}.root"
+            f"{output}/output/l{int(phantom_length_mm)}_ps{name}.root"
         )
         phase_space.attributes = [
             "EventID",
@@ -131,10 +131,6 @@ def add_detector(name, translation, attach_to_phantom=False):
         ):
             add_detector(str(x), [0 * mm, 0 * mm, x * mm], True)
 
-    # Particle stats
-    stat = sim.add_actor("SimulationStatisticsActor", "stat")
-    stat.output_filename = f"{output}/wepl{int(phantom_length_mm)}_stats.txt"
-
     sim.run()
 
 
@@ -154,7 +150,9 @@ def process_phantom_length(
     wepls_phantom_length = []
     elosses_phantom_length = []
 
-    data = uproot.concatenate(f"{output}/l{int(phantom_length)}_*.root", library="np")
+    data = uproot.concatenate(
+        f"{output}/output/l{int(phantom_length)}_*.root", library="np"
+    )
     pv(
         verbose,
         "Loaded",