reconstruction.py: ionBeamEnergy = str(os.environ["PBEAM"])

options/reconstruction.py

@@ -23,6 +23,11 @@ if "JUGGLER_DETECTOR_VERSION" in os.environ:
 
 compact_path = os.path.join(detector_path, detector_name)
 
+if "PBEAM" in os.environ:
+    ionBeamEnergy = str(os.environ["PBEAM"])
+else:
+    ionBeamEnergy = 100
+
 # RICH reconstruction
 qe_data = [(1.0, 0.25), (7.5, 0.25),]
 
@@ -705,7 +710,7 @@ fast_ff = FFSmearedParticles("fast_ff",
         enableB0         = True,
         enableRP         = True,
         enableOMD        = True,
-        ionBeamEnergy    = 100,
+        ionBeamEnergy    = ionBeamEnergy,
         crossingAngle    = -0.025)
 algorithms.append(fast_ff)