summaryrefslogtreecommitdiff
path: root/generator
diff options
context:
space:
mode:
authorStan Seibert <stan@mtrr.org>2011-08-16 18:54:15 -0400
committerStan Seibert <stan@mtrr.org>2011-08-16 18:54:15 -0400
commite5f51bde2015b12301062ab7903038ba2fcf3e34 (patch)
tree6a491e5badcf1ac995e7261aef492fdf430731e7 /generator
parentdfd4da145cb2c980f68785b33d1596945f67a672 (diff)
downloadchroma-e5f51bde2015b12301062ab7903038ba2fcf3e34.tar.gz
chroma-e5f51bde2015b12301062ab7903038ba2fcf3e34.tar.bz2
chroma-e5f51bde2015b12301062ab7903038ba2fcf3e34.zip
Preheat GEANT4 process by running a short event at startup
Diffstat (limited to 'generator')
-rw-r--r--generator/g4gen.py10
1 files changed, 7 insertions, 3 deletions
diff --git a/generator/g4gen.py b/generator/g4gen.py
index 4416c45..cec53f7 100644
--- a/generator/g4gen.py
+++ b/generator/g4gen.py
@@ -49,7 +49,13 @@ class G4Generator(object):
self.tracking_action = G4chroma.PhotonTrackingAction()
gRunManager.SetUserAction(self.tracking_action)
gRunManager.Initialize()
-
+
+ #preinitialize the process by running a simple event
+ self.generate_photons(event.Event(event_id=0, particle_name='e-',
+ gen_position=(0,0,0),
+ gen_direction=(1,0,0),
+ gen_total_energy=1.0))
+
def create_g4material(self, material):
g4material = G4Material('world_material', material.density * g / cm3,
len(material.composition))
@@ -133,8 +139,6 @@ if __name__ == '__main__':
import time
import optics
gen = G4Generator(optics.water)
- # prime things
- gen.generate_photons(event.Event('e-', (0,0,0), (1,0,0), 1.0))
start = time.time()
n = 0