source: trunk/MagicSoft/Simulation/Corsika/Mmcs/avept.f

Last change on this file was 286, checked in by harald, 25 years ago
This is the start point for further developments of the Magic Monte Carlo Simulation written by Jose Carlos Gonzales. Now it is under control of one CVS repository for the whole collaboration. Everyone should use this CVS repository for further developments.
File size: 2.8 KB
Line 
1 SUBROUTINE AVEPT( ECM,SLOG )
2
3C-----------------------------------------------------------------------
4C AVE(RAGE) PT (TRANSVERSE MOMENTUM)
5C
6C CALCULATES AVERAGE RATIO PT(PARTICLE)/PT(PION) DEPENDING ON ENERGY
7C THE DEPENDENCE OF PT ON ENERGY IS DONE IN ROUTINE PTRAM/PTRAN
8C THIS SUBROUTINE IS CALLED FROM HDPM
9C ARGUMENTS:
10C ECM = ENERGY IN THE CM SYSTEM
11C SLOG = LOG(S) ( = LOG(ECM**2) )
12C-----------------------------------------------------------------------
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
44C-----------------------------------------------------------------------
45
46CC IF ( DEBUG ) WRITE(MDEBUG,*) 'AVEPT : ECM =',SNGL(ECM)
47
48C AVERAGE TRANSVERSE MOMENTUM
49
50C ... 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
57C ... 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
64C ... 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
71C SET AVERAGE PT RELATED TO AVERAGE PT FOR PIONS
72C ... FOR STRANGE BARYONS (=AVPH)
73 AVPH = 1.3D0 * (1.45D0 * AVPN - 0.45D0 * AVPK) / AVPT
74C ... 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
Note: See TracBrowser for help on using the repository browser.