Message: GPS particle type Not Logged In (login)
 Next-in-Thread Next-in-Thread
 Next-in-Forum Next-in-Forum

Question GPS particle type 

Forum: Event and Track Management
Date: 01 Mar, 2011
From: Yuri Soloviev <Yuri Soloviev>

   Attachment:
      http://hypernews.slac.stanford.edu/HyperNews/geant4/get/AUX/2011/03/01/05.47-17415-gps1.log

  Hello everyone, 

 I am having the problem with using GPS as particle gun. 
 I have to simulate a spectrum of incident particle therefore 
 I changed G4ParticleGun to GPS. 
 Here is the macro : 

 /run/verbose 2
/event/verbose 2
/tracking/verbose 0
# 
# Set particle type, it's energy
# and number of events generated
#
/gps/verbose 2
/gps/particle neutron
/gps/position -135. 0. 0. cm
/gps/direction 1. 0. 0.
/gps/ene/type Arb
/gps/ene/emspec 0
/gps/hist/type arb
/gps/hist/point   0.   0.
/gps/hist/point   2.5 21609.
/gps/hist/point   7.5 56574.
/gps/hist/point  12.5 72274.
/gps/hist/point  17.5 76944.
/gps/hist/point  22.5 75408.
/gps/hist/point  27.5 72675.
/gps/hist/point  32.5 59494.
/gps/hist/point  37.5 43723.
/gps/hist/point  42.5 27558.
/gps/hist/point  47.5 18851.
/gps/hist/point  52.5 13470.
/gps/hist/point  57.5 9913.
/gps/hist/point  62.5 7274.
/gps/hist/point  67.5 5269.
/gps/hist/point  72.5 3603.
/gps/hist/point  77.5 2751.
/gps/hist/point  82.5 1914.
/gps/hist/point  87.5 1573.
/gps/hist/point  92.5 1612.
/gps/hist/point  97.5 1656.
/gps/hist/point 102.5 1863.
/gps/hist/point 107.5 1854.
/gps/hist/point 112.5 1091.
/gps/hist/point 117.5  501.
/gps/hist/point 122.5  334.
/gps/hist/point 127.5	43.
/gps/hist/point 132.5	17.
/gps/hist/point 137.5	21.
/gps/hist/point 142.5	17.
/gps/hist/point 147.5	11.
/gps/hist/point 152.5	 4.
/gps/hist/point 157.5	14.
/gps/hist/point 162.5	 8.
/gps/hist/point 167.5	 6.
/gps/hist/point 172.5	 3.
/gps/hist/point 177.5	 9.
/gps/hist/point 182.5	 8.
/gps/hist/point 187.5	13.
/gps/hist/point 192.5	 5.
/gps/hist/point 197.5	 4.
/gps/hist/point 202.5	 2.
/gps/hist/point 207.5	 3.
/gps/hist/point 212.5	 5.
/gps/hist/point 217.5	 0.
/gps/hist/point 222.5	 0.
/gps/hist/point 227.5	 3.
/gps/hist/point 232.5	 2.
/gps/hist/point 237.5	 3.
/gps/hist/point 242.5	 0.
/gps/hist/point 247.5	 0.
/gps/hist/point 252.5	 0.
/gps/hist/point 257.5	 0.
/gps/hist/point 262.5	 5.
/gps/hist/point 267.5	 0.
/gps/hist/point 272.5	 0.
/gps/hist/point 277.5	 0.
/gps/hist/point 282.5	 0.
/gps/hist/point 287.5	 0.
/gps/hist/point 292.5	 0.
/gps/hist/point 297.5	 0.
/gps/hist/point 302.5	 0.
/gps/hist/point 307.5	 0.
/gps/hist/point 312.5	 1.
/gps/hist/point 317.5	 0.
/gps/hist/point 322.5	 0.
/gps/hist/inter Lin
/gps/number 1  
/run/beamOn 10000

 after running job first time I discovered that the momentum spectra 
 of neutron differs very much from the arb distribution on input. 
 I set gps/verbose 2 and I got the message (see atteched file) 
 "Error: particle not defined". 
 Though in the message from geant4 inself I see that particle 
 is neutron. (see attached file) 

 It looks for me like gps didnt recognize particle type and took 
 geantino as default. Then it recalculated the energy from momentum distribution
 assuming zero mass. And finally I see crazy momentum distribution on output. 
 I checked the calculation for a few points.

 Any idea, comments and proposals are welcome. 

 Best regards 
                         Yuri Soloviev

Inline Depth:
 1 1
 All All
Outline Depth:
 1 1
 2 2
 All All
Add message: (add)

1 None: Re: GPS particle type   (michel maire - 02 Mar, 2011)
(_ None: Re: GPS particle type   (Yuri Soloviev - 03 Mar, 2011)
 Add Message Add Message
to: "GPS particle type"

 Subscribe Subscribe

This site runs SLAC HyperNews version 1.11-slac-98, derived from the original HyperNews


[ Geant 4 Home | Geant 4 HyperNews | Search | Request New Forum | Feedback ]