| 1 |
|
|---|
| 2 | SUBROUTINE GUSTEP
|
|---|
| 3 | *
|
|---|
| 4 | * User routine called at the end of each tracking step
|
|---|
| 5 | *
|
|---|
| 6 | #include "geant321/gcflag.inc"
|
|---|
| 7 | #include "geant321/gckine.inc"
|
|---|
| 8 | #include "geant321/gcking.inc"
|
|---|
| 9 | #include "geant321/gctrak.inc"
|
|---|
| 10 | *
|
|---|
| 11 | #include "process.inc"
|
|---|
| 12 | #include "histo.inc"
|
|---|
| 13 | *
|
|---|
| 14 | character*20 cdum
|
|---|
| 15 | *
|
|---|
| 16 | * *** Debug event and store tracks for drawing
|
|---|
| 17 | IF (IDEBUG.NE.0) CALL GPCXYZ
|
|---|
| 18 | IF (IDEBUG.NE.0) CALL GPGKIN
|
|---|
| 19 | IF ((ISWIT(1).EQ.1).AND.(CHARGE.NE.0.)) CALL GSXYZ
|
|---|
| 20 | IF (ISWIT(1).EQ.2) CALL GSXYZ
|
|---|
| 21 | *
|
|---|
| 22 | * *** if no process: return
|
|---|
| 23 | IF (NMEC.EQ.0) return
|
|---|
| 24 | *
|
|---|
| 25 | * *** count nb of invoked processes
|
|---|
| 26 | DO IM = 1,NMEC
|
|---|
| 27 | IPROC = LMEC(IM)
|
|---|
| 28 | IF (IPROC.EQ.21) IPROC = 12
|
|---|
| 29 | IF (IPROC.LE.12) NBCALL(IPROC) = NBCALL(IPROC)+1
|
|---|
| 30 | ENDDO
|
|---|
| 31 | *
|
|---|
| 32 | * *** sum track length for discrete processes
|
|---|
| 33 | if ((iproc.gt.5)) then
|
|---|
| 34 | nbTot = nbTot + 1
|
|---|
| 35 | sumTrak = sumTrak + sleng
|
|---|
| 36 | sumTrak2 = sumTrak2 + sleng*sleng
|
|---|
| 37 | endif
|
|---|
| 38 | *
|
|---|
| 39 | * *** plot final state
|
|---|
| 40 | *
|
|---|
| 41 | * scattered primary particle (if still alive)
|
|---|
| 42 | if (istop.eq.0) then
|
|---|
| 43 | id = 1
|
|---|
| 44 | if (histo(id)) call hfill (id,gekin/histUnit(id),0.,1.)
|
|---|
| 45 | id = 2
|
|---|
| 46 | if (histo(id)) call hfill (id,vect(4),0.,1.)
|
|---|
| 47 | endif
|
|---|
| 48 | *
|
|---|
| 49 | * *** secondaries
|
|---|
| 50 | if (ngkine.gt.0) then
|
|---|
| 51 | do lp = 1,ngkine
|
|---|
| 52 | ipar = gkin(5,lp) + 0.1
|
|---|
| 53 | call gfpart(ipar,cdum,ndum,gmass,gcharg,dum,dum,ndum)
|
|---|
| 54 | ekin = gkin(4,lp) - gmass
|
|---|
| 55 | pc = sqrt(gkin(1,lp)**2 + gkin(2,lp)**2 + gkin(3,lp)**2)
|
|---|
| 56 | cost = gkin(1,lp)/pc
|
|---|
| 57 | if (gcharg.ne.0.) id = 3
|
|---|
| 58 | if (gcharg.eq.0.) id = 5
|
|---|
| 59 | if (histo(id)) call hfill (id,ekin/histUnit(id),0.,1.)
|
|---|
| 60 | id = id + 1
|
|---|
| 61 | if (histo(id)) call hfill (id,cost,0.,1.)
|
|---|
| 62 | enddo
|
|---|
| 63 | endif
|
|---|
| 64 | *
|
|---|
| 65 | * *** stop the tracking
|
|---|
| 66 | istop = 1
|
|---|
| 67 | *
|
|---|
| 68 | END
|
|---|
