Run a short dummy GEANT4 event to force it to load the tables before

the timer starts on the real events.

1 files changed, 8 insertions, 0 deletions

diff --git a/sim.py b/sim.py
index 6432e7d..b8deef6 100755
--- a/sim.py
+++ b/sim.py
@@ -74,6 +74,14 @@ if __name__ == '__main__':
     gpu_worker.setup_propagate()
     gpu_worker.setup_daq(max(detector.pmtids))
 
+    print >>sys.stderr, 'Loading tables in GEANT4...'
+    # Do easy event to force tables to load, throw away photons
+    generator.generate_photons(particle_name='e-', 
+                               total_energy=1.5,
+                               position=(0,0,0),
+                               direction=(1,0,0))
+
+
     # Create output file
     f = ROOT.TFile(output_filename, 'RECREATE')
     ev, T = root.make_tree('T')