| 1 | SUBROUTINE AVEPT( ECM,SLOG )
|
|---|
| 2 |
|
|---|
| 3 | C-----------------------------------------------------------------------
|
|---|
| 4 | C AVE(RAGE) PT (TRANSVERSE MOMENTUM)
|
|---|
| 5 | C
|
|---|
| 6 | C CALCULATES AVERAGE RATIO PT(PARTICLE)/PT(PION) DEPENDING ON ENERGY
|
|---|
| 7 | C THE DEPENDENCE OF PT ON ENERGY IS DONE IN ROUTINE PTRAM/PTRAN
|
|---|
| 8 | C THIS SUBROUTINE IS CALLED FROM HDPM
|
|---|
| 9 | C ARGUMENTS:
|
|---|
| 10 | C ECM = ENERGY IN THE CM SYSTEM
|
|---|
| 11 | C SLOG = LOG(S) ( = LOG(ECM**2) )
|
|---|
| 12 | C-----------------------------------------------------------------------
|
|---|
| 13 |
|
|---|
| 14 | IMPLICIT NONE
|
|---|
| 15 | *KEEP,AVPT.
|
|---|
| 16 | COMMON /AVPT/ AVPT,AVPK,AVPN,AVPH,AVPE
|
|---|
| 17 | DOUBLE PRECISION AVPT,AVPK,AVPN,AVPH,AVPE
|
|---|
| 18 | *KEEP,RUNPAR.
|
|---|
| 19 | COMMON /RUNPAR/ FIXHEI,THICK0,HILOECM,HILOELB,
|
|---|
| 20 | * STEPFC,NRRUN,NSHOW,PATAPE,MONIIN,
|
|---|
| 21 | * MONIOU,MDEBUG,NUCNUC,
|
|---|
| 22 | * CETAPE,
|
|---|
| 23 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
|---|
| 24 | * N1STTR,MDBASE,
|
|---|
| 25 | * DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
|---|
| 26 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
|---|
| 27 | * ,GHEISH,GHESIG
|
|---|
| 28 | COMMON /RUNPAC/ DSN,HOST,USER
|
|---|
| 29 | DOUBLE PRECISION FIXHEI,THICK0,HILOECM,HILOELB
|
|---|
| 30 | REAL STEPFC
|
|---|
| 31 | INTEGER NRRUN,NSHOW,PATAPE,MONIIN,MONIOU,MDEBUG,NUCNUC,
|
|---|
| 32 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
|---|
| 33 | * N1STTR,MDBASE
|
|---|
| 34 | INTEGER CETAPE
|
|---|
| 35 | CHARACTER*79 DSN
|
|---|
| 36 | CHARACTER*20 HOST,USER
|
|---|
| 37 |
|
|---|
| 38 | LOGICAL DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
|---|
| 39 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
|---|
| 40 | * ,GHEISH,GHESIG
|
|---|
| 41 | *KEND.
|
|---|
| 42 |
|
|---|
| 43 | DOUBLE PRECISION ECM,SLOG
|
|---|
| 44 | C-----------------------------------------------------------------------
|
|---|
| 45 |
|
|---|
| 46 | CC IF ( DEBUG ) WRITE(MDEBUG,*) 'AVEPT : ECM =',SNGL(ECM)
|
|---|
| 47 |
|
|---|
| 48 | C AVERAGE TRANSVERSE MOMENTUM
|
|---|
| 49 |
|
|---|
| 50 | C ... FOR PIONS (=AVPT)
|
|---|
| 51 | IF ( ECM .LT. 132.D0 ) THEN
|
|---|
| 52 | AVPT = 0.3D0 + 6.272D-3 * SLOG
|
|---|
| 53 | ELSE
|
|---|
| 54 | AVPT = ( 0.442105D0 + 0.016276D0 * SLOG )**2
|
|---|
| 55 | ENDIF
|
|---|
| 56 |
|
|---|
| 57 | C ... FOR KAONS (=AVPK)
|
|---|
| 58 | IF ( ECM .LT. 131.D0 ) THEN
|
|---|
| 59 | AVPK = 1.27D0 * AVPT
|
|---|
| 60 | ELSE
|
|---|
| 61 | AVPK = (0.403146D0 + 0.0281D0 * SLOG)**2
|
|---|
| 62 | ENDIF
|
|---|
| 63 |
|
|---|
| 64 | C ... FOR NUCLEONS (=AVPN)
|
|---|
| 65 | IF ( ECM .LT. 102.D0 ) THEN
|
|---|
| 66 | AVPN = 1.39D0 * AVPT
|
|---|
| 67 | ELSE
|
|---|
| 68 | AVPN = (0.389873D0 + 0.034127D0 * SLOG)**2
|
|---|
| 69 | ENDIF
|
|---|
| 70 |
|
|---|
| 71 | C SET AVERAGE PT RELATED TO AVERAGE PT FOR PIONS
|
|---|
| 72 | C ... FOR STRANGE BARYONS (=AVPH)
|
|---|
| 73 | AVPH = 1.3D0 * (1.45D0 * AVPN - 0.45D0 * AVPK) / AVPT
|
|---|
| 74 | C ... FOR ETA MESONS (=AVPE)
|
|---|
| 75 | AVPE = 1.3D0 * (0.88D0 * AVPK + 0.12D0 * AVPN) / AVPT
|
|---|
| 76 | AVPK = 1.3D0 * AVPK / AVPT
|
|---|
| 77 | AVPN = 1.3D0 * AVPN / AVPT
|
|---|
| 78 | AVPT = 1.3D0
|
|---|
| 79 |
|
|---|
| 80 | IF ( DEBUG ) WRITE(MDEBUG,100)
|
|---|
| 81 | * SNGL(AVPT),SNGL(AVPK),SNGL(AVPN),SNGL(AVPH),SNGL(AVPE)
|
|---|
| 82 | 100 FORMAT(' AVEPT : AVPT,AVPK,AVPN,AVPH,AVPE=',5F12.5)
|
|---|
| 83 |
|
|---|
| 84 | RETURN
|
|---|
| 85 | END
|
|---|
