| 1 | SUBROUTINE PARRAP
|
|---|
| 2 |
|
|---|
| 3 | C-----------------------------------------------------------------------
|
|---|
| 4 | C PAR(TICLE) RAP(IDITY)
|
|---|
| 5 | C
|
|---|
| 6 | C ROUTINE GIVES THE NEW PARTICLES OF HDPM THEIR RAPIDITIES
|
|---|
| 7 | C THIS SUBROUTINE IS CALLED FROM HDPM
|
|---|
| 8 | C-----------------------------------------------------------------------
|
|---|
| 9 |
|
|---|
| 10 | IMPLICIT DOUBLE PRECISION (A-H,O-Z)
|
|---|
| 11 | *KEEP,INTER.
|
|---|
| 12 | COMMON /INTER/ AVCH,AVCH3,DC0,DLOG,DMLOG,ECMDIF,ECMDPM,ELAB,
|
|---|
| 13 | * FNEUT,FNEUT2,GNU,PLAB,POSC2,POSC3,POSN2,POSN3,
|
|---|
| 14 | * RC3TO2,S,SEUGF,SEUGP,SLOG,SLOGSQ,SMLOG,
|
|---|
| 15 | * WIDC2,WIDC3,WIDN2,WIDN3,YCM,YY0,ZN,
|
|---|
| 16 | * IDIF,ITAR
|
|---|
| 17 | DOUBLE PRECISION AVCH,AVCH3,DC0,DLOG,DMLOG,ECMDIF,ECMDPM,ELAB,
|
|---|
| 18 | * FNEUT,FNEUT2,GNU,PLAB,POSC2,POSC3,POSN2,POSN3,
|
|---|
| 19 | * RC3TO2,S,SEUGF,SEUGP,SLOG,SLOGSQ,SMLOG,
|
|---|
| 20 | * WIDC2,WIDC3,WIDN2,WIDN3,YCM,YY0,ZN
|
|---|
| 21 | INTEGER IDIF,ITAR
|
|---|
| 22 | *KEEP,NEWPAR.
|
|---|
| 23 | COMMON /NEWPAR/ EA,PT2,PX,PY,TMAS,YR,ITYP,
|
|---|
| 24 | * IA1,IA2,IB1,IB2,IC1,IC2,ID1,ID2,IE1,IE2,IF1,IF2,
|
|---|
| 25 | * IG1,IG2,IH1,IH2,II1,II2,IJ1,NTOT
|
|---|
| 26 | DOUBLE PRECISION EA(3000),PT2(3000),PX(3000),PY(3000),TMAS(3000),
|
|---|
| 27 | * YR(3000)
|
|---|
| 28 | INTEGER ITYP(3000),
|
|---|
| 29 | * IA1,IA2,IB1,IB2,IC1,IC2,ID1,ID2,IE1,IE2,IF1,IF2,
|
|---|
| 30 | * IG1,IG2,IH1,IH2,II1,II2,IJ1,NTOT
|
|---|
| 31 | *KEEP,RANDPA.
|
|---|
| 32 | COMMON /RANDPA/ FAC,U1,U2,RD,NSEQ,ISEED,KNOR
|
|---|
| 33 | DOUBLE PRECISION FAC,U1,U2
|
|---|
| 34 | REAL RD(3000)
|
|---|
| 35 | INTEGER ISEED(103,10),NSEQ
|
|---|
| 36 | LOGICAL KNOR
|
|---|
| 37 | *KEEP,RUNPAR.
|
|---|
| 38 | COMMON /RUNPAR/ FIXHEI,THICK0,HILOECM,HILOELB,
|
|---|
| 39 | * STEPFC,NRRUN,NSHOW,PATAPE,MONIIN,
|
|---|
| 40 | * MONIOU,MDEBUG,NUCNUC,
|
|---|
| 41 | * CETAPE,
|
|---|
| 42 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
|---|
| 43 | * N1STTR,MDBASE,
|
|---|
| 44 | * DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
|---|
| 45 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
|---|
| 46 | * ,GHEISH,GHESIG
|
|---|
| 47 | COMMON /RUNPAC/ DSN,HOST,USER
|
|---|
| 48 | DOUBLE PRECISION FIXHEI,THICK0,HILOECM,HILOELB
|
|---|
| 49 | REAL STEPFC
|
|---|
| 50 | INTEGER NRRUN,NSHOW,PATAPE,MONIIN,MONIOU,MDEBUG,NUCNUC,
|
|---|
| 51 | * SHOWNO,ISHW,NOPART,NRECS,NBLKS,MAXPRT,NDEBDL,
|
|---|
| 52 | * N1STTR,MDBASE
|
|---|
| 53 | INTEGER CETAPE
|
|---|
| 54 | CHARACTER*79 DSN
|
|---|
| 55 | CHARACTER*20 HOST,USER
|
|---|
| 56 |
|
|---|
| 57 | LOGICAL DEBDEL,DEBUG,FDECAY,FEGS,FIRSTI,FIXINC,FIXTAR,
|
|---|
| 58 | * FIX1I,FMUADD,FNKG,FPRINT,FDBASE
|
|---|
| 59 | * ,GHEISH,GHESIG
|
|---|
| 60 | *KEND.
|
|---|
| 61 |
|
|---|
| 62 | REAL RAND(3000)
|
|---|
| 63 | C-----------------------------------------------------------------------
|
|---|
| 64 |
|
|---|
| 65 | IF ( DEBUG ) WRITE(MDEBUG,*) 'PARRAP: NTOT=',NTOT
|
|---|
| 66 |
|
|---|
| 67 | C PROTON ANTIPROTON PAIRS
|
|---|
| 68 | CALL RMMAR( RAND(3),IJ1-2,1 )
|
|---|
| 69 | DO 1013 K = 3,IB1
|
|---|
| 70 | C GENERATION OF RAPIDITY FOR EXTRA PARTICLES FROM TARGET.
|
|---|
| 71 | IF ( K .LE. IA2 ) THEN
|
|---|
| 72 | YR(K) = RANNOR(-POSC3,WIDC3)
|
|---|
| 73 | ELSE
|
|---|
| 74 | C GENERATION OF RAPIDITY FOR PARTICLES FROM PP-COLLISION AND PROJECTILE
|
|---|
| 75 | IF ( RAND(K) .LE. 0.5 ) THEN
|
|---|
| 76 | YR(K) = RANNOR(POSC2,WIDC2)
|
|---|
| 77 | ELSE
|
|---|
| 78 | YR(K) = RANNOR(-POSC2,WIDC2)
|
|---|
| 79 | ENDIF
|
|---|
| 80 | ENDIF
|
|---|
| 81 | 1013 CONTINUE
|
|---|
| 82 |
|
|---|
| 83 | C K+ K- PAIRS
|
|---|
| 84 | DO 1014 K = IB1+1,IC1
|
|---|
| 85 | IF ( K .LE. IB2 ) THEN
|
|---|
| 86 | C GENERATION OF RAPIDITY FOR EXTRA PARTICLES FROM TARGET.
|
|---|
| 87 | YR(K) = RANNOR(-POSC3,WIDC3)
|
|---|
| 88 | ELSE
|
|---|
| 89 | C GENERATION OF RAPIDITY FOR PARTICLES FROM PP-COLLISION AND PROJECTILE
|
|---|
| 90 | IF ( RAND(K) .LE. 0.5 ) THEN
|
|---|
| 91 | YR(K) = RANNOR(POSC2,WIDC2)
|
|---|
| 92 | ELSE
|
|---|
| 93 | YR(K) = RANNOR(-POSC2,WIDC2)
|
|---|
| 94 | ENDIF
|
|---|
| 95 | ENDIF
|
|---|
| 96 | 1014 CONTINUE
|
|---|
| 97 |
|
|---|
| 98 | C CHARGED HYPERON PAIRS
|
|---|
| 99 | DO 1015 K = IC1+1,ID1
|
|---|
| 100 | IF ( K .LE. IC2 ) THEN
|
|---|
| 101 | C GENERATION OF RAPIDITY FOR EXTRA PARTICLES FROM TARGET.
|
|---|
| 102 | YR(K) = RANNOR(-POSC3,WIDC3)
|
|---|
| 103 | ELSE
|
|---|
| 104 | C GENERATION OF RAPIDITY FOR PARTICLES FROM PP-COLLISION AND PROJECTILE
|
|---|
| 105 | IF ( RAND(K) .LE. 0.5 ) THEN
|
|---|
| 106 | YR(K) = RANNOR(POSC2,WIDC2)
|
|---|
| 107 | ELSE
|
|---|
| 108 | YR(K) = RANNOR(-POSC2,WIDC2)
|
|---|
| 109 | ENDIF
|
|---|
| 110 | ENDIF
|
|---|
| 111 | 1015 CONTINUE
|
|---|
| 112 |
|
|---|
| 113 | C PI +-
|
|---|
| 114 | DO 1017 K = ID1+1,IE1
|
|---|
| 115 | IF ( K .LE. ID2 ) THEN
|
|---|
| 116 | C GENERATION OF RAPIDITY FOR EXTRA PARTICLES FROM TARGET.
|
|---|
| 117 | YR(K) = RANNOR(-POSC3,WIDC3)
|
|---|
| 118 | ELSE
|
|---|
| 119 | C GENERATION OF RAPIDITY FOR PARTICLES FROM PP-COLLISION AND PROJECTILE
|
|---|
| 120 | IF ( RAND(K) .LE. 0.5 ) THEN
|
|---|
| 121 | YR(K) = RANNOR(POSC2,WIDC2)
|
|---|
| 122 | ELSE
|
|---|
| 123 | YR(K) = RANNOR(-POSC2,WIDC2)
|
|---|
| 124 | ENDIF
|
|---|
| 125 | ENDIF
|
|---|
| 126 | 1017 CONTINUE
|
|---|
| 127 |
|
|---|
| 128 | C NEUTRON ANTINEUTRON PAIRS
|
|---|
| 129 | DO 1021 K = IE1+1,IF1
|
|---|
| 130 | IF ( K .LE. IE2 ) THEN
|
|---|
| 131 | C GENERATION OF RAPIDITY FOR EXTRA PARTICLES FROM TARGET.
|
|---|
| 132 | YR(K) = RANNOR(-POSC3,WIDC3)
|
|---|
| 133 | ELSE
|
|---|
| 134 | C GENERATION OF RAPIDITY FOR PARTICLES FROM PP-COLLISION AND PROJECTILE
|
|---|
| 135 | IF ( RAND(K) .LE. 0.5 ) THEN
|
|---|
| 136 | YR(K) = RANNOR(POSC2,WIDC2)
|
|---|
| 137 | ELSE
|
|---|
| 138 | YR(K) = RANNOR(-POSC2,WIDC2)
|
|---|
| 139 | ENDIF
|
|---|
| 140 | ENDIF
|
|---|
| 141 | 1021 CONTINUE
|
|---|
| 142 |
|
|---|
| 143 | C K0L K0S PAIRS
|
|---|
| 144 | DO 1022 K = IF1+1,IG1
|
|---|
| 145 | IF ( K .LE. IF2 ) THEN
|
|---|
| 146 | C GENERATION OF RAPIDITY FOR EXTRA PARTICLES FROM TARGET.
|
|---|
| 147 | YR(K) = RANNOR(-POSC3,WIDC3)
|
|---|
| 148 | ELSE
|
|---|
| 149 | C GENERATION OF RAPIDITY FOR PARTICLES FROM PP-COLLISION AND PROJECTILE
|
|---|
| 150 | IF ( RAND(K) .LE. 0.5 ) THEN
|
|---|
| 151 | YR(K) = RANNOR(POSC2,WIDC2)
|
|---|
| 152 | ELSE
|
|---|
| 153 | YR(K) = RANNOR(-POSC2,WIDC2)
|
|---|
| 154 | ENDIF
|
|---|
| 155 | ENDIF
|
|---|
| 156 | 1022 CONTINUE
|
|---|
| 157 |
|
|---|
| 158 | C NEUTRAL HYPERON PAIRS
|
|---|
| 159 | DO 1023 K = IG1+1,IH1
|
|---|
| 160 | IF ( K .LE. IG2 ) THEN
|
|---|
| 161 | C GENERATION OF RAPIDITY FOR EXTRA PARTICLES FROM TARGET.
|
|---|
| 162 | YR(K) = RANNOR(-POSC3,WIDC3)
|
|---|
| 163 | ELSE
|
|---|
| 164 | C GENERATION OF RAPIDITY FOR PARTICLES FROM PP-COLLISION AND PROJECTILE
|
|---|
| 165 | IF ( RAND(K) .LE. 0.5 ) THEN
|
|---|
| 166 | YR(K) = RANNOR(POSC2,WIDC2)
|
|---|
| 167 | ELSE
|
|---|
| 168 | YR(K) = RANNOR(-POSC2,WIDC2)
|
|---|
| 169 | ENDIF
|
|---|
| 170 | ENDIF
|
|---|
| 171 | 1023 CONTINUE
|
|---|
| 172 |
|
|---|
| 173 | C ETA
|
|---|
| 174 | DO 1025 K = IH1+1,II1
|
|---|
| 175 | IF ( K .LE. IH2 ) THEN
|
|---|
| 176 | C GENERATION OF RAPIDITY FOR EXTRA PARTICLES FROM TARGET.
|
|---|
| 177 | YR(K) = RANNOR(-POSN3,WIDN3)
|
|---|
| 178 | ELSE
|
|---|
| 179 | C GENERATION OF RAPIDITY FOR PARTICLES FROM PP-COLLISION AND PROJECTILE
|
|---|
| 180 | IF ( RAND(K) .LE. 0.5 ) THEN
|
|---|
| 181 | YR(K) = RANNOR(POSN2,WIDN2)
|
|---|
| 182 | ELSE
|
|---|
| 183 | YR(K) = RANNOR(-POSN2,WIDN2)
|
|---|
| 184 | ENDIF
|
|---|
| 185 | ENDIF
|
|---|
| 186 | 1025 CONTINUE
|
|---|
| 187 |
|
|---|
| 188 | C PI(0)
|
|---|
| 189 | DO 1026 K = II1+1,IJ1
|
|---|
| 190 | IF ( K .LE. II2 ) THEN
|
|---|
| 191 | C GENERATION OF RAPIDITY FOR EXTRA PARTICLES FROM TARGET.
|
|---|
| 192 | YR(K) = RANNOR(-POSN3,WIDN3)
|
|---|
| 193 | ELSE
|
|---|
| 194 | C GENERATION OF RAPIDITY FOR PARTICLES FROM PP-COLLISION AND PROJECTILE
|
|---|
| 195 | IF ( RAND(K) .LE. 0.5 ) THEN
|
|---|
| 196 | YR(K) = RANNOR(POSN2,WIDN2)
|
|---|
| 197 | ELSE
|
|---|
| 198 | YR(K) = RANNOR(-POSN2,WIDN2)
|
|---|
| 199 | ENDIF
|
|---|
| 200 | ENDIF
|
|---|
| 201 | 1026 CONTINUE
|
|---|
| 202 |
|
|---|
| 203 | RETURN
|
|---|
| 204 | END
|
|---|
