Rev | Line | |
---|
[807] | 1 | |
---|
| 2 | SUBROUTINE UGINIT |
---|
| 3 | * |
---|
| 4 | * To initialise GEANT/USER program and read data cards |
---|
| 5 | * |
---|
| 6 | #include "geant321/gcphys.inc" |
---|
| 7 | #include "geant321/gccuts.inc" |
---|
| 8 | #include "ugmate.inc" |
---|
| 9 | #include "ggvalues.inc" |
---|
| 10 | * |
---|
| 11 | CHARACTER*20 filnam |
---|
| 12 | * |
---|
| 13 | * *** Define the GEANT parameters |
---|
| 14 | CALL GINIT |
---|
| 15 | |
---|
| 16 | * *** read data cards |
---|
| 17 | PRINT *, 'G3 > gives the filename of the data cards to be read:' |
---|
| 18 | READ (*,'(A)') filnam |
---|
| 19 | IF (filnam.EQ.' ') filnam = 'run01.dat' |
---|
| 20 | OPEN (unit=5,file=filnam,status='unknown',form='formatted') |
---|
| 21 | * |
---|
| 22 | * *** define data cards |
---|
| 23 | call FFKEY ('MATE' ,Imat,1,'INTEGER') |
---|
| 24 | call FFKEY ('G4VAL',g4value(1),200,'REAL') |
---|
| 25 | * |
---|
| 26 | * *** read data cards |
---|
| 27 | CALL GFFGO |
---|
| 28 | * |
---|
| 29 | * *** overwrite some initialization |
---|
| 30 | iloss = 4 |
---|
| 31 | ppcutm = cutele |
---|
| 32 | * |
---|
| 33 | CALL GZINIT |
---|
| 34 | CALL GPART |
---|
| 35 | * |
---|
| 36 | * *** materials definition |
---|
| 37 | CALL UGIMAT |
---|
| 38 | * |
---|
| 39 | * *** Energy loss and cross-sections initialisations |
---|
| 40 | CALL GPHYSI |
---|
| 41 | * |
---|
| 42 | CALL GPRINT('MATE',0) |
---|
| 43 | * |
---|
| 44 | * *** ready for geant3/geant4 comparison |
---|
| 45 | call UG3G4 |
---|
| 46 | * |
---|
| 47 | END |
---|
Note: See
TracBrowser
for help on using the repository browser.