diff --git a/scripts/momentum_scan.rb b/scripts/momentum_scan.rb
index 0a00cd7..2d47b0f 100755
--- a/scripts/momentum_scan.rb
+++ b/scripts/momentum_scan.rb
@@ -7,8 +7,8 @@
 require 'pycall/import'
 
 ## settings
-NumEvents      = 100                 # number of events per fixed momentum
-NumPoints      = 20                  # number of momenta to sample
+NumEvents      = 50                  # number of events per fixed momentum
+NumPoints      = 10                  # number of momenta to sample
 PoolSize       = 6                   # number of parallel threads to run
 OutputDir      = 'out/momentum_scan' # output directory ( ! will be overwritten ! )
 RunSimRec      = true                # if false, do not run simulation+reconstruction, only draw the result
diff --git a/simulate.py b/simulate.py
index 4cc119c..b3c0580 100755
--- a/simulate.py
+++ b/simulate.py
@@ -478,7 +478,8 @@
         f'--runType {runType}',
         f'--compactFile {compactFile}',
         f'--outputFile {outputFileName_npsim}',
-        "--part.userParticleHandler=''", # necessary for opticalphotons truth output
+        # "--part.userParticleHandler=''", # necessary for opticalphotons truth output (REMOVE THIS?)
+        '--part.minimalKineticEnergy "0*eV"', # allow opticalphotons in trackers (dRICH) to be included in output
         # '--random.seed 1',
         # '--part.keepAllParticles True',
         ]