lfhcal: scale event count with momentum

benchmarks/lfhcal/Snakefile

+def get_n_events(wildcards):
+    energy = float(wildcards.ENERGY.replace("GeV", "").replace("MeV", "e-3"))
+    n_events = 1000
+    n_events = int(n_events // (energy ** 0.5))
+    return n_events
+
+
 rule lfhcal_sim:
     input:
         steering_file="EPIC/EVGEN/SINGLE/{PARTICLE}/{ENERGY}/{PHASE_SPACE}/{PARTICLE}_{ENERGY}_{PHASE_SPACE}.steer",
@@ -12,7 +19,7 @@ rule lfhcal_sim:
         PHASE_SPACE="3to50deg",
         INDEX="\d{4}",
     params:
-        N_EVENTS=1000
+        N_EVENTS=get_n_events,
     shell:
         """
 ddsim \