| 9ef1c2d9 |
1 | |
| 2 | parameter (mmrye=1,mmryk=1) !memory saving factor |
| 3 | |
| 4 | parameter (kollmx=15000/mmrye,npommx=80) |
| 5 | parameter (mamx=210/mmryk,ngrmx=1,ntymi=0,ntymx=2) |
| 6 | |
| 7 | !---------------------------------------------------------------------- |
| 8 | ! POMERON |
| 9 | !---------------------------------------------------------------------- |
| 10 | |
| 11 | double precision om1int,om1intc |
| 12 | * ,atilde,btildep,btildepp |
| 13 | common/comega/om1int(kollmx),om1intc(kollmx) |
| 14 | common/comtilde/atilde(ntymi:ntymx,kollmx) |
| 15 | *,btildep(ntymi:ntymx,kollmx),btildepp(ntymi:ntymx,kollmx) |
| 16 | *,ntymin |
| 17 | |
| 18 | !---------------------------------------------------------------------- |
| 19 | ! POMERON CONFIGURATION |
| 20 | !---------------------------------------------------------------------- |
| 21 | |
| 22 | common/col/koll,iproj(kollmx),itarg(kollmx),coord(4,kollmx) |
| 23 | common/col1/kremn(mamx,mamx,2),lremn(mamx,2) |
| 24 | dimension kproj(mamx,mamx),ktarg(mamx,mamx) |
| 25 | *,lproj(mamx),ltarg(mamx) |
| 26 | equivalence (kproj(1,1),kremn(1,1,1)),(ktarg(1,1),kremn(1,1,2)) |
| 27 | equivalence (lproj(1),lremn(1,1)),(ltarg(1),lremn(1,2)) |
| 28 | double precision gammaV |
| 29 | common/col6/kproj3(mamx,mamx),ktarg3(mamx,mamx) |
| 30 | *,lproj3(mamx),ltarg3(mamx) |
| 31 | *,epsilongp(kollmx,0:1),epsilongt(kollmx,0:1) |
| 32 | *,epsilongs(kollmx,0:1),gammaV(kollmx) |
| 33 | *,epsilongf(kollmx),epsilongb(kollmx) |
| 34 | common/col4/bk(kollmx),bhpr(npommx,kollmx),bkx(kollmx) |
| 35 | *,bky(kollmx) |
| 36 | common/col7/zparremn(kollmx,2) |
| 37 | real zparpro(kollmx),zpartar(kollmx) |
| 38 | equivalence (zparpro(1),zparremn(1,1)),(zpartar(1),zparremn(1,2)) |
| 39 | |
| 40 | !---------------------------------------------------------------------- |
| 41 | ! koll: number of relevant pairs |
| 42 | ! iproj(k): projectile index for k-th pair |
| 43 | ! itarg(k): target index for k-th pair |
| 44 | ! lproj(ip): number of pairs linked to projectile ip |
| 45 | ! kproj(ip,l): corresponding pair index of the l-th pair |
| 46 | ! ltarg(it): number of pairs linked to target it |
| 47 | ! ktarg(it,l): corresponding pair index of the l-th pair |
| 48 | ! bk(k): b value of the k-th pair (=sqrt(bx**2.+by**2.)) |
| 49 | ! bkx(k): bx value of the k-th pair |
| 50 | ! bky(k): by value of the k-th pair |
| 51 | ! lproj3(ip): number of pairs linked to projectile ip (no diffraction) |
| 52 | ! kproj3(ip,l): corresponding pair index of the l-th pair (no diffraction) |
| 53 | ! ltarg3(it): number of pairs linked to target it (no diffraction) |
| 54 | ! ktarg3(it,l): corresponding pair index of the l-th pair (no diffraction) |
| 55 | !---------------------------------------------------------------------- |
| 56 | |
| 57 | common/cems3/nprt(kollmx),npr(0:3,kollmx) |
| 58 | common/cems4/nprmx(kollmx) |
| 59 | |
| 60 | !---------------------------------------------------------------------- |
| 61 | ! nprt(k): number of pomerons of collision k |
| 62 | ! npr(id,k): number of pomerons of type id of collision k |
| 63 | ! (id: 0=zero, 1=soft, 2=reggeon, 3=hard) |
| 64 | ! nprmx(k): maximum number of pomerons of collision k |
| 65 | !---------------------------------------------------------------------- |
| 66 | |
| 67 | common/cems2/itpr(kollmx),idpr(npommx,kollmx),ivpr(npommx,kollmx) |
| 68 | *,idfpr(npommx,kollmx) |
| 69 | *,nvpr(npommx,kollmx),nbkpr(npommx,kollmx) |
| 70 | |
| 71 | !---------------------------------------------------------------------- |
| 72 | ! itpr(k): interaction type of collision k |
| 73 | ! (0=nothing, 1=inelastic, 2=diffractive) |
| 74 | ! idpr(n,k): type of pomeron (n,k) |
| 75 | ! (0=zero, 1=soft, 2=reggeon, 3=hard) |
| 76 | ! ivpr(n,k): status of pomeron (n,k) (0=virtual, 1=real, 2=backup) |
| 77 | ! nbkpr(n,k): position of the backup pomeron |
| 78 | ! nvpr(n,k): position of the pomeron for this backup pomeron |
| 79 | ! idfpr(n,k): link of pomeron (n,k) (0=no, 1=pro&tar, 2=pro, 3=tar) |
| 80 | ! idfpr=-1 means backup pomeron |
| 81 | !---------------------------------------------------------------------- |
| 82 | |
| 83 | double precision xpr,ypr,xppr,xmpr,xpprbor,xmprbor |
| 84 | common/cems1/xpr(npommx,kollmx),ypr(npommx,kollmx) |
| 85 | common/cems1c/xppr(npommx,kollmx),xmpr(npommx,kollmx) |
| 86 | common/cems1d/xpprbor(npommx,kollmx),xmprbor(npommx,kollmx) |
| 87 | common/cems1e/nemispr(2,npommx,kollmx),ptprboo(2,npommx,kollmx) |
| 88 | |
| 89 | !---------------------------------------------------------------------- |
| 90 | ! xpr(n,k): x=x+x- of pomeron (n,k) |
| 91 | ! ypr(n,k): rapidity of pomeron (n,k) |
| 92 | ! xpprbor(n,k): x+ of in Born parton of pomeron (n,k) |
| 93 | ! xmprbor(n,k): x- of in Born parton of pomeron (n,k) |
| 94 | ! nemispr(i,n,k): number of up (i=1) and down (i=2) emissions |
| 95 | !---------------------------------------------------------------------- |
| 96 | |
| 97 | double precision xp1pr,xp2pr,xm1pr,xm2pr,xxp1pr,xyp1pr |
| 98 | * ,xxp2pr,xyp2pr,xxm1pr,xym1pr,xxm2pr,xym2pr |
| 99 | common/cems7/xp1pr(npommx,kollmx),xp2pr(npommx,kollmx) |
| 100 | * ,xm1pr(npommx,kollmx),xm2pr(npommx,kollmx) |
| 101 | * ,idp1pr(npommx,kollmx),idp2pr(npommx,kollmx) |
| 102 | * ,idm1pr(npommx,kollmx),idm2pr(npommx,kollmx) |
| 103 | common/cems8/xxp1pr(npommx,kollmx),xyp1pr(npommx,kollmx) |
| 104 | * ,xxp2pr(npommx,kollmx),xyp2pr(npommx,kollmx) |
| 105 | * ,xxm1pr(npommx,kollmx),xym1pr(npommx,kollmx) |
| 106 | * ,xxm2pr(npommx,kollmx),xym2pr(npommx,kollmx) |
| 107 | |
| 108 | !---------------------------------------------------------------------- |
| 109 | ! xp1pr(n,k): x of string end p1 p1 p2 |
| 110 | ! xp2pr(n,k): x of string end p2 \ / |
| 111 | ! xm1pr(n,k): x of string end m1 \ / |
| 112 | ! xm2pr(n,k): x of string end m2 \/ |
| 113 | ! idp1pr(n,k): id of string end p1 ( 0=0 ) /\ |
| 114 | ! idp2pr(n,k): id of string end p2 ( 1=s 2=v ) / \ |
| 115 | ! idm1pr(n,k): id of string end m1 ( 4=ss 5=vv ) / \ |
| 116 | ! idm2pr(n,k): id of string end m2 ( 8=f ) m1 m2 |
| 117 | ! xxp1pr(n,k),xyp1pr(n,k): px,py of string end p1 |
| 118 | ! xxp2pr(n,k),xyp2pr(n,k): px,py of string end p2 |
| 119 | ! xxm1pr(n,k),xym1pr(n,k): px,py of string end m1 |
| 120 | ! xxm2pr(n,k),xym2pr(n,k): px,py of string end m2 |
| 121 | !---------------------------------------------------------------------- |
| 122 | |
| 123 | double precision xme |
| 124 | common/cems11/xme(mamx) |
| 125 | |
| 126 | !---------------------------------------------------------------------- |
| 127 | ! xme(j): x- exchange for nucleon j |
| 128 | !---------------------------------------------------------------------- |
| 129 | |
| 130 | common/cems22/idhpr(npommx,kollmx),idrpr(npommx,kollmx) |
| 131 | * ,idsppr(npommx,kollmx),idstpr(npommx,kollmx) |
| 132 | |
| 133 | !---------------------------------------------------------------------- |
| 134 | ! idhpr(n,k): flav hard Pom (n,k) (0=gg,1=qg,2=gq,3=qq) |
| 135 | ! idrpr(n,k): flav Regg (n,k) (0=ss,1=0va0,2=v00a,3=0as0,4=s00a,5=0sv0,6=v00s) |
| 136 | ! idsppr(n,k): flav soft P (n,k) (proj) (0=ss,1=vs,2=sa,3=va,4=dd,5=svv,6=aas) |
| 137 | ! idstpr(n,k): flav soft P (n,k) (targ) (0=ss,1=vs,2=sa,3=va,4=dd,5=svv,6=aas) |
| 138 | !---------------------------------------------------------------------- |
| 139 | |
| 140 | |
| 141 | !---------------------------------------------------------------------- |
| 142 | ! REMNANT CONFIGURATION |
| 143 | !---------------------------------------------------------------------- |
| 144 | |
| 145 | common/remn2/zzremn(mamx,2) |
| 146 | common/remn/icremn(2,mamx,2) |
| 147 | dimension icproj(2,mamx),ictarg(2,mamx) |
| 148 | double precision xpp,xmp,xpt,xmt,xxp,xyp,xxt,xyt |
| 149 | double precision xpz,xmz,xxz,xyz |
| 150 | common/cemsr/xpz(mamx,2),xmz(mamx,2) |
| 151 | * ,xxz(mamx,2),xyz(mamx,2) |
| 152 | * ,idz(mamx,2),ivz(mamx,2) |
| 153 | * ,isz(mamx,2),iaz(mamx,2) |
| 154 | dimension xpp(mamx),xmp(mamx),xpt(mamx),xmt(mamx) |
| 155 | * ,xxp(mamx),xyp(mamx),xxt(mamx),xyt(mamx) |
| 156 | * ,idp(mamx),ivp(mamx),ivt(mamx),idt(mamx) |
| 157 | * ,isp(mamx),iap(mamx),ist(mamx),iat(mamx) |
| 158 | |
| 159 | equivalence (icproj(1,1),icremn(1,1,1)) |
| 160 | * ,(ictarg(1,1),icremn(1,1,2)) |
| 161 | equivalence (xpp(1),xpz(1,1)),(xmt(1),xpz(1,2)) |
| 162 | * ,(xmp(1),xmz(1,1)),(xpt(1),xmz(1,2)) |
| 163 | * ,(xxp(1),xxz(1,1)),(xxt(1),xxz(1,2)) |
| 164 | * ,(xyp(1),xyz(1,1)),(xyt(1),xyz(1,2)) |
| 165 | * ,(idp(1),idz(1,1)),(idt(1),idz(1,2)) |
| 166 | * ,(ivp(1),ivz(1,1)),(ivt(1),ivz(1,2)) |
| 167 | * ,(isp(1),isz(1,1)),(ist(1),isz(1,2)) |
| 168 | * ,(iap(1),iaz(1,1)),(iat(1),iaz(1,2)) |
| 169 | |
| 170 | !---------------------------------------------------------------------- |
| 171 | ! xpp(ip),xmp(ip),xpt(it),xmt(it): x+,x- of projectile ip / target it |
| 172 | ! xxp(ip),xyp(ip),xxt(it),xyt(it): px,py of projectile ip / target it |
| 173 | ! idp(ip),idt(it): diquark counter |
| 174 | ! ivp(ip),ivt(it): valenve quark counter |
| 175 | ! isp(ip),ist(it): sea quark counter |
| 176 | ! iap(ip),iat(it): anti valence quark counter |
| 177 | !---------------------------------------------------------------------- |
| 178 | |
| 179 | double precision xpmn,xtmn,xpmx,xtmx,xppmx,xptmx,xmpmx,xmtmx |
| 180 | * ,xppmn ,xptmn,xmpmn,xmtmn,xpos,xtos |
| 181 | double precision xzmn,xzmx,xpzmx,xmzmx |
| 182 | * ,xpzmn ,xmzmn,xzos |
| 183 | common/cemsr1/xzmn(mamx,2),xzmx(mamx,2) |
| 184 | * ,xpzmx(mamx,2) |
| 185 | * ,xmzmx(mamx,2) |
| 186 | * ,xpzmn(mamx,2) |
| 187 | * ,xmzmn(mamx,2) |
| 188 | * ,xzos(mamx,2) |
| 189 | dimension xpmn(mamx),xtmn(mamx),xpmx(mamx),xtmx(mamx) |
| 190 | * ,xppmx(mamx),xptmx(mamx) |
| 191 | * ,xmpmx(mamx),xmtmx(mamx) |
| 192 | * ,xppmn(mamx),xptmn(mamx) |
| 193 | * ,xmpmn(mamx),xmtmn(mamx) |
| 194 | * ,xpos(mamx),xtos(mamx) |
| 195 | |
| 196 | equivalence ( xpmn(1), xzmn(1,1)),( xtmn(1), xzmn(1,2)) |
| 197 | * ,( xpmx(1), xzmx(1,1)),( xtmx(1), xzmx(1,2)) |
| 198 | * ,(xppmx(1),xpzmx(1,1)),(xptmx(1),xpzmx(1,2)) |
| 199 | * ,(xmpmx(1),xmzmx(1,1)),(xmtmx(1),xmzmx(1,2)) |
| 200 | * ,(xppmn(1),xpzmn(1,1)),(xptmn(1),xpzmn(1,2)) |
| 201 | * ,(xmpmn(1),xmzmn(1,1)),(xmtmn(1),xmzmn(1,2)) |
| 202 | * ,( xpos(1), xzos(1,1)),( xtos(1), xzos(1,2)) |
| 203 | |
| 204 | !---------------------------------------------------------------------- |
| 205 | ! xpmn(ip),xtmn(it),xpmx(ip),xtmx(it): x proj/targ minimum/maximum |
| 206 | ! xppmx(ip),xptmx(it),xmpmx(ip),xmtmx(it): x+/x- proj/targ maximum |
| 207 | ! xppmn(ip),xptmn(it),xmpmn(ip),xmtmn(it): x+/x- proj/targ minimum |
| 208 | ! xpos(ip),xtos(it): x value for on shell condition proj/targ |
| 209 | !---------------------------------------------------------------------- |
| 210 | |
| 211 | common/cemsr3/kolp(mamx),kolt(mamx),npp(mamx),npt(mamx) |
| 212 | common/cemsr4/iez(mamx,2),ifz(mamx,2) |
| 213 | dimension iep(mamx),iet(mamx),ifp(mamx),ift(mamx) |
| 214 | equivalence (iep(1),iez(1,1)),(iet(1),iez(1,2)) |
| 215 | * ,(ifp(1),ifz(1,1)),(ift(1),ifz(1,2)) |
| 216 | |
| 217 | !---------------------------------------------------------------------- |
| 218 | ! kolp(ip),kolt(it): number of cut Pomerons of proj ip / targ it |
| 219 | ! iep(ip),iet(it): excitation type of proj ip / targ it |
| 220 | ! (0=none,1=inelastic,2=diffractive) |
| 221 | ! ifp(ip),ift(it): remnant flip flag proj/targ (0=no,1=yes) |
| 222 | ! if3dp(ip), if3dp(it): triple pomeron diffraction (if 1) |
| 223 | !---------------------------------------------------------------------- |
| 224 | |
| 225 | |
| 226 | !---------------------------------------------------------------------- |
| 227 | ! transfer from emsini |
| 228 | !---------------------------------------------------------------------- |
| 229 | |
| 230 | double precision amzmn,ampmn,amtmn,ammn,amemn,amemx,amprmn,xis |
| 231 | &,delxi,ddelx,r3pompi,r3pommi,r3pompi3,r3pommi3 |
| 232 | parameter (ntypmx=19) |
| 233 | common/cems12/ammn(0:ntypmx),amprmn(0:ntypmx) |
| 234 | common/cemsr6/amzmn(0:1,2),amemn(0:1,0:6),amemx(0:2) |
| 235 | dimension ampmn(0:1),amtmn(0:1) |
| 236 | equivalence (ampmn(0),amzmn(0,1)), (amtmn(0),amzmn(0,2)) |
| 237 | common/cemsr13/xis,delxi,ddelx,r3pompi,r3pommi |
| 238 | &,r3pompi3,r3pommi3,ucfpro,ucftar |
| 239 | |
| 240 | !---------------------------------------------------------------------- |
| 241 | ! STORAGE |
| 242 | !---------------------------------------------------------------------- |
| 243 | |
| 244 | double precision |
| 245 | *xxpr,xxpr0,yx,yx0,xxppr,xxmpr,xxppr0,xxmpr0 |
| 246 | *,xmtx,xmtx0,xmtlx,xmtlx0 |
| 247 | *,xppx,xppx0,xpplx,xpplx0 |
| 248 | *,xppst,xmpst,xposst,xpest,xppstx0,xmpstx0,xppstx,xmpstx |
| 249 | *,xmtst,xptst,xtosst,xmest,xmtstx0,xptstx0,xmtstx,xptstx |
| 250 | dimension nprx(0:3),nprx0(0:3) |
| 251 | *,xppst(mamx),xmpst(mamx),xposst(mamx),xpest(mamx) |
| 252 | *,xmtst(mamx),xptst(mamx),xtosst(mamx),xmest(mamx) |
| 253 | |
| 254 | |
| 255 | common/emsstore/ |
| 256 | *xxpr,xxpr0,yx,yx0,xxppr,xxmpr,xxppr0,xxmpr0 |
| 257 | *,xmtx,xmtx0,xmtlx,xmtlx0 |
| 258 | *,xppx,xppx0,xpplx,xpplx0 |
| 259 | *,itx,itx0,idx,idx0 |
| 260 | *,npplx0,nptlx0,npplx,nptlx |
| 261 | *,nppx,nppx0,nptx,nptx0,nprx,nprx0,nprtx,nprtx0 |
| 262 | *,xppst,xmpst,xposst,xpest,xppstx0,xmpstx0,xppstx,xmpstx |
| 263 | *,xmtst,xptst,xtosst,xmest,xmtstx0,xptstx0,xmtstx,xptstx |
| 264 | |
| 265 | |
| 266 | !---------------------------------------------------------------------- |
| 267 | ! utilities |
| 268 | !---------------------------------------------------------------------- |
| 269 | |
| 270 | double precision facto,fctrl |
| 271 | common/factori/facto(0:npommx),fctrl(0:npommx), nfctrl |
| 272 | |
| 273 | |
| 274 | |
| 275 | |
| 276 | |
git://git.uio.no / u / mrichter / AliRoot.git / blame