LMU/musrsim
LMU/musrsim
5
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Workaround on program precision.

Browse Source
This commit is contained in:
prokscha 2007-02-08 10:02:17 +00:00
parent ddff2b0694
commit 9494ea7b29
1 changed files with 7 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
geant3/src/lemsr/get_direction.F
View File

@ -32,12 +32,18 @@ c
real*4 x,y,z,r real*4 x,y,z,r
real*4 t_in_sample,tof,beta,gamma_sqr ! to add relaxation in sample real*4 t_in_sample,tof,beta,gamma_sqr ! to add relaxation in sample
real*4 mass, mass_sqr, pz_sqr, z0 real*4 mass, mass_sqr, pz_sqr, z0
real*4 arg
c c
c----------------------------------------------------------------------------- c-----------------------------------------------------------------------------
c c
if ( lsets(2) .eq. 0 ) then ! positron created in MCP2 or sample if ( lsets(2) .eq. 0 ) then ! positron created in MCP2 or sample
t_in_sample = tofg*1.e9 ! time in nsec t_in_sample = tofg*1.e9 ! time in nsec
p = bfield(5) / 100. * exp(-beam_parameter(4)*t_in_sample) ! expo relaxation arg = beam_parameter(4)*t_in_sample
if (arg > 10. ) then
p = 0.
else
p = bfield(5) / 100. * dexp(-dble(arg)) ! expo relaxation
endif
else else
x = vect(1) x = vect(1)
y = vect(2) y = vect(2)