| Line | |
|---|
| 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.
