| 29 |
* |
* |
| 30 |
* Let's trasnform energy and momentum in MeV. Remember, GEANT is GeV |
* Let's trasnform energy and momentum in MeV. Remember, GEANT is GeV |
| 31 |
* |
* |
| 32 |
open(77,file='testgig.dat',access='append') |
* open(77,file='testgig.dat',access='append') |
| 33 |
En_gamma=En_g/1000. |
En_gamma=En_g/1000. |
| 34 |
px_gamma=px_g/1000. |
px_gamma=px_g/1000. |
| 35 |
py_gamma=py_g/1000. |
py_gamma=py_g/1000. |
| 67 |
enddo |
enddo |
| 68 |
|
|
| 69 |
endif |
endif |
| 70 |
write(77,*) Number_N, step, en_gamma |
* write(77,*) Number_N, step, en_gamma |
| 71 |
close(77) |
* close(77) |
| 72 |
RETURN |
RETURN |
| 73 |
end |
end |
| 74 |
|
|