Message: Re: How to define multiple sources intensities from the beam power with GPS  Not Logged In (login) 
Dear Leonardo,
Please find below an example of how to use two sources in Geant4. Best regards Laurent # two beams in a generator # # beam #1 # default intensity is 1 now change to 5./gps/source/intensity 5. # /gps/particle proton /gps/pos/type Beam # # the incident surface is in the yz plane/gps/pos/rot1 0 1 0 /gps/pos/rot2 0 0 1 # # the beam spot is centered at the origin and is # of 1d gaussian shape with a 1 mm central plateau/gps/pos/shape Circle /gps/pos/centre 0. 0. 0. mm /gps/pos/radius 1. mm /gps/pos/sigma_r .2 mm # # the beam is travelling along the X_axis with 5 degrees dispersion/gps/ang/rot1 0 0 1 /gps/ang/rot2 0 1 0 /gps/ang/type beam1d /gps/ang/sigma_r 5. deg # # the beam energy is in gaussian profile centered at 400 MeV/gps/ene/type Gauss /gps/ene/mono 400 MeV /gps/ene/sigma 50. MeV
# beam #2 # 2x the instensity of beam #1/gps/source/add 10. # # this is a electron beam/gps/particle e /gps/pos/type Beam # it beam spot is of 2d gaussian profile with a 1x2 mm2 central plateau # it is in the xy plane centred at the orgin/gps/pos/centre 0. 0. 0. mm /gps/pos/halfx 0.5 mm /gps/pos/halfy 1. mm /gps/pos/sigma_x 0.1 mm # the spread in y direction is stronger /gps/pos/sigma_y 0.2 mm # #the beam is travelling along Z_axis/gps/ang/type beam2d /gps/ang/sigma_x 2. deg /gps/ang/sigma_y 1. deg # gaussian energy profile /gps/ene/type Gauss /gps/ene/mono 600 MeV /gps/ene/sigma 50. MeV Resulted Distribution Plots: On Fri, 19 Feb 2016 12:09:58 GMT, Leonardo Ghizoni wrote: > Hello, > > I'm simulating a photon shutter on a synchrotron accelerator, and I'm a > little stuck on defining multiple sources with GPS. > > I have two sources. One is an XRay source (100 keV photons) with 200W, > and the other source is of bremsstrahlung radiation (3 GeV photons) but > with 2.513E3 W. > > My question is how to define both sources in the same run. I tried the > following, but not sure on the intensities: > > # beam #1 > > /gps/source/intensity ?? > > /gps/particle gamma > > /gps/ene/type Mono > > /gps/ene/mono 100 keV > > /gps/pos/type Plane > > /gps/pos/shape Rectangle > > /gps/pos/rot1 0 0 1 > > /gps/pos/rot2 0 1 0 > > /gps/pos/halfx 0.015 mm > > /gps/pos/halfy 0.15 mm > > /gps/pos/halfz 0.15 mm > > /gps/pos/centre 0. 0. 50. cm > > /gps/direction 0 0 1 > > # beam #2 > > /gps/source/intensity ?? > > /gps/particle gamma > > /gps/ene/type Mono > > /gps/ene/mono 3 GeV > > /gps/pos/type Plane > > /gps/pos/shape Rectangle > > /gps/pos/rot1 0 0 1 > > /gps/pos/rot2 0 1 0 > > /gps/pos/halfx 0.015 mm > > /gps/pos/halfy 0.15 mm > > /gps/pos/halfz 0.15 mm > > /gps/pos/centre 0. 0. 50. cm > > /gps/direction 0 0 1 > > Thanks a lot! >

to: 