diff options
author | Kevin Wells <kwells@icloud.com> | 2017-08-07 04:39:33 +0000 |
---|---|---|
committer | tlatorre <tlatorre@uchicago.edu> | 2021-05-09 08:42:39 -0700 |
commit | f9b2c7648cf0b20553c23856d2fc079f9d4a7b2e (patch) | |
tree | 72ebc055950c5a608e9b4235fc3f1fd9f8f857c3 | |
parent | c784d169ae8f170e10ec2d5a61387e2b01643468 (diff) | |
download | chroma-f9b2c7648cf0b20553c23856d2fc079f9d4a7b2e.tar.gz chroma-f9b2c7648cf0b20553c23856d2fc079f9d4a7b2e.tar.bz2 chroma-f9b2c7648cf0b20553c23856d2fc079f9d4a7b2e.zip |
SetParticleEnergy() requires kinetic energy only, not total energy.
See; https://www-zeuthen.desy.de/geant4/geant4.9.3.b01/classG4ParticleGun.html
-rw-r--r-- | chroma/generator/g4gen.py | 5 |
1 files changed, 2 insertions, 3 deletions
diff --git a/chroma/generator/g4gen.py b/chroma/generator/g4gen.py index 27f1454..3e185a9 100644 --- a/chroma/generator/g4gen.py +++ b/chroma/generator/g4gen.py @@ -130,9 +130,8 @@ class G4Generator(object): try: for vertex in vertices: self.particle_gun.SetParticleByName(vertex.particle_name) - mass = G4ParticleTable.GetParticleTable().FindParticle(vertex.particle_name).GetPDGMass() - total_energy = vertex.ke*MeV + mass - self.particle_gun.SetParticleEnergy(total_energy) + #mass = G4ParticleTable.GetParticleTable().FindParticle(vertex.particle_name).GetPDGMass() + self.particle_gun.SetParticleEnergy(vertex.ke*MeV) # Must be float type to call GEANT4 code pos = np.asarray(vertex.pos, dtype=np.float64) |