EPOS/epos167/epos.inc

   1 c---------------------------------------------------------------------------
   2 c                         dimensions
   3 c---------------------------------------------------------------------------
   4
   5       parameter (mmry=1)   !memory saving factor
   6
   7       parameter (mxptl=200000/mmry) !max nr of particles in epos ptl list
   8       parameter (nmxhep=9990)       !max nr of particles in hep ptl list
   9       parameter (myptl=1000)        !max nr of droplets in epos ptl list
  10       parameter (nzeta=60)          !max nr of zeta bins for droplets
  11       parameter (nflav=6)           !max nr of flavors
  12       parameter (mxstr=20000/mmry)  !max nr of strings in epos string list
  13       parameter (mystr=20000/mmry)
  14       parameter (mxtau=4,mxvol=10,mxeps=16)
  15       parameter (mxtrig=99,mxidco=99)
  16       parameter (mxpri=200)
  17       parameter (mxbins=10000)
  18       parameter (matau=10,mxcoox=40,mxcooy=10)
  19       parameter (mxnucl=20)
  20       parameter (mxhisarg=100)
  21       parameter (idxD0=0,idxD1=2,idxD=1,nclha=4,nclegy=100)
  22       parameter (mamxx=250)
  23       parameter (mxjerr=10)
  24
  25 c---------------------------------------------------------------------------
  26 c                   epos event common block
  27 c---------------------------------------------------------------------------
  28
  29       common/cevt/phievt,nevt,bimevt,kolevt,koievt,pmxevt,egyevt,npjevt
  30      *,ntgevt,npnevt,nppevt,ntnevt,ntpevt,jpnevt,jppevt,jtnevt,jtpevt
  31      *,xbjevt,qsqevt,nglevt,zppevt,zptevt,minfra,maxfra
  32
  33 c     nevt .......... error code. 1=valid event, 0=invalid event
  34 c     bimevt ........ absolute value of impact parameter
  35 c     phievt ........ angle of impact parameter
  36 c     kolevt ........ number of collisions
  37 c     koievt ........ number of inelastic collisions
  38 c     pmxevt ........ reference momentum
  39 c     egyevt ........ pp cm energy (hadron) or string energy (lepton)
  40 c     npjevt ........ number of primary projectile participants
  41 c     ntgevt ........ number of primary target participants
  42 c     npnevt ........ number of primary projectile neutron spectators
  43 c     nppevt ........ number of primary projectile proton spectators
  44 c     ntnevt ........ number of primary target neutron spectators
  45 c     ntpevt ........ number of primary target proton spectators
  46 c     jpnevt ........ number of absolute projectile neutron spectators
  47 c     jppevt ........ number of absolute projectile proton spectators
  48 c     jtnevt ........ number of absolute target neutron spectators
  49 c     jtpevt ........ number of absolute target proton spectators
  50 c     xbjevt ........ bjorken x for dis
  51 c     qsqevt ........ q**2 for dis
  52 c     sigtot ........ total cross section
  53 c     nglevt ........ number of collisions acc to  Glauber
  54 c     zppevt ........ average Z-parton-proj
  55 c     zptevt ........ average Z-parton-targ
  56
  57 c---------------------------------------------------------------------------
  58 c                   epos particle list common block
  59 c---------------------------------------------------------------------------
  60
  61       common/cptl/nptl,pptl(5,mxptl),iorptl(mxptl),idptl(mxptl)
  62      *,istptl(mxptl),tivptl(2,mxptl),ifrptl(2,mxptl),jorptl(mxptl)
  63      *,xorptl(4,mxptl),ibptl(4,mxptl),ityptl(mxptl)
  64       common/c1ptl/itsptl(mxptl)
  65
  66 c     nptl .......... current particle index (=number of ptls stored)
  67 c     idptl(i) ...... particle id
  68 c     pptl(1,i) ..... x-component of particle momentum
  69 c     pptl(2,i) ..... y-component of particle momentum
  70 c     pptl(3,i) ..... z-component of particle momentum
  71 c     pptl(4,i) ..... particle energy
  72 c     pptl(5,i) ..... particle mass
  73 c     iorptl(i) ..... particle number of father (if .le. 0 : no father)
  74 c     jorptl(i) ..... particle number of mother (if .le. 0 : no mother)
  75 c     istptl(i) ..... status: 40 and 41 : Remnant
  76 c                             30 and 31 : Pomeron
  77 c                             20 and 21 : Parton
  78 c                             10 and 11 : Droplet
  79 c                             00 and 01 : Particle
  80 c                            last digit = 0 : last generation
  81 c                            last digit = 1 : not last generation
  82 c     xorptl(1,i) ... x-component of formation point
  83 c     xorptl(2,i) ... y-component of formation point
  84 c     xorptl(3,i) ... z-component of formation point
  85 c     xorptl(4,i) ... formation time
  86 c     tivptl(1,i) ... formation time (always in the pp-cms!)
  87 c     tivptl(2,i) ... destruction time (always in the pp-cms!)
  88 c     ityptl(i)  .... type of particles origin:
  89 c                         10-19: target
  90 c                         20-29: soft Pom
  91 c                         30-39: hard Pom
  92 c                         40-49: projectile
  93 c                         50: string, droplet
  94 c     itsptl(i) ..... string type of particles origin (if string)
  95
  96       common/c2ptl/iaaptl(mxptl),radptl(mxptl)
  97       common/c3ptl/desptl(mxptl),dezptl(mxptl)
  98       common/c4ptl/nptlbd
  99       common/c6ptl/rinptl(mxptl),vrad,inbxxx
 100       common/c8ptl/qsqptl(mxptl),zpaptl(mxptl)
 101
 102
 103 c---------------------------------------------------------------------------
 104 c                   hep standard event commonblock.
 105 c---------------------------------------------------------------------------
 106
 107       double precision phep,vhep
 108
 109       common/hepevt/nevhep,nhep,isthep(nmxhep),idhep(nmxhep),
 110      &jmohep(2,nmxhep),jdahep(2,nmxhep),phep(5,nmxhep),vhep(4,nmxhep)
 111
 112 c---------------------------------------------------------------------------
 113 c
 114 c         nevhep      -   event number
 115 c         nhep        -   number of entries in the event record
 116 c
 117 c         isthep(i)   -   status code
 118 c         idhep(i)    -   particle id (particle data group standard)
 119 c
 120 c         jmohep(1,i) -   position of mother particle in list
 121 c         jmohep(2,i) -   position of second mother particle in list
 122 c         jdahep(1,i) -   position of first daughter in list
 123 c         jdahep(2,i) -   position of first daughter in list
 124 c
 125 c         phep(1,i)   -   p_x momentum in gev/c
 126 c         phep(2,i)   -   p_y momentum in gev/c
 127 c         phep(3,i)   -   p_z momentum in gev/c
 128 c         phep(4,i)   -   energy in gev
 129 c         phep(5,i)   -   mass in gev/c**2
 130 c
 131 c         vhep(1,i)   -   x position of production vertex in mm
 132 c         vhep(2,i)   -   y position of production vertex in mm
 133 c         vhep(3,i)   -   z position of production vertex in mm
 134 c         vhep(4,i)   -   time of production  in mm/c
 135 c
 136 c          (note:  1 mm = 10^-12 fm = 5.07 10^-12 1/gev)
 137
 138 c------------------------------------------------------------------------
 139 c  Parameters set in sr aaset and variables to communicate between moduls
 140 c------------------------------------------------------------------------
 141
 142       common/files/ifop,ifmt,ifch,ifcx,ifhi,ifdt,ifcp,ifdr
 143       character*80 fnch,fnhi,fndt,fnii,fnid,fnie,fnrj,fnmt
 144      * ,fngrv,fncp,fnnx,fncs,fndr,fnhy
 145       common/fname/  fnch, fnhi, fndt, fnii, fnid, fnie, fnrj, fnmt
 146      * ,fngrv,fncp,fnnx,fncs,fndr,fnhy
 147       common/nfname/nfnch,nfnhi,nfndt,nfnii,nfnid,nfnie,nfnrj,nfnmt
 148      *,nfngrv,nfncp,nfnnx,nfncs,nfndr,nfnhy
 149       common/resc2/delvol,deleps,dlzeta,etafac,facnuc,taurea,epscri(3)
 150       common/resc4/epsdfm
 151       common/frag1/ndecay,maxres,pud,pudk,pudr,strcut,diqcut
 152       common/frag2/pdiqua,pdiquak,ioptf,delrex,ptfraqq,ptfrasr,ptfrakd
 153       common/frag3/aouni,pbreak,pbreakk,ptfra,ptfrair,ptfrak,itflav
 154       common/frag4/pude,pdiquae,pdiquar,pbreake,pbreakr,ptfrae,ptfrab
 155       common/frag5/pbreakd,pudd,pdiquad,ptfradr
 156       common/frag6/qmass(0:6),puds,isospin(0:6)
 157       real ptq
 158       common/hadr1/pnll,ptq,exmass,cutmss,wproj,wtarg
 159       common/hadr10/rstrau(4),rstrad(4),rstras(4),rstrasi
 160       common/wgtqrk/wgtqqq(4),wgtval,wgtsea,wgtdiq
 161       double precision timeini,timefin
 162       common/time1/timeini,timefin
 163       common/ciotst/iotst1,iotst2,iotst3,iotst4
 164       common/resc1/taumin,deltau,factau,numtau,amsiac,amprif
 165       common/resc3/dscale,cepara,iceopt,delamf,deuamf
 166       common/sprio/ispherio,icotabm,icotabr,icocore
 167       common/urqm1/iurqmd
 168       common/incon/cutico,dssico
 169       common/othe1/istore,istmax,gaumx,irescl,ntrymx,nclean,iopdg,ioidch
 170       common/othe2/ifrade,iframe,idecay,jdecay
 171       common/othe3/jframe,kframe
 172       common/jpsif/jpsi,jpsifi,taumx,nsttau,sigj,ijphis,ijtauan
 173       common/strlt/iopenu,themas
 174       common/appli/iappl,model
 175       common/events/nevent,nfull,nfreeze,ninicon
 176       common/enrgy/egymin,egymax,elab,ecms,ekin
 177       common/prnt1/iprmpt,ish,ishsub,irandm,irewch,iecho,modsho,idensi
 178       common/lept1/engy,elepti,elepto,angmue,icinpu
 179       common/nucl1/laproj,maproj,latarg,matarg,core,fctrmx
 180       common/nucl2/bmaxim,bminim,phimax,phimin
 181       common/wana1/ymximi,imihis,iclhis,iwtime,wtimet,wtimei,wtimea
 182       common/wana2/isphis,ispall,wtmini,wtstep,iwcent,iana,nbdky
 183       common/drop4/asuhax(7),asuhay(7)
 184       common/gribo/grigam,grirsq,gridel,grislo,gricel,sigppi,sigppd
 185       common/drop3/bag4rt,dezzer,amuseg,taunll,yradmx,facts,factb,factq
 186       common/drop2/rcoll,ylongmx,nsegsu,nsegce,facecc,yradpp,yradmi
 187       common/drop7/ptclu,yradpi,yradpx,ioclude,iocluin,ioquen
 188       common/metr1/iospec,iocova,iopair,iozero,ioflac,iomom
 189       common/metr2/nadd,iograc,iocite,ioceau,iociau
 190       common/hadr2/iomodl,idproj,idtarg,wexcit
 191       common/hadr25/idprojin,idtargin,rexdifi(4),rexndii(4),irdmpr
 192       common/metr3/iostat,ioinco,ionlat,ioobsv,iosngl,iorejz,iompar
 193       common/metr4/ioinfl,ioinct,iowidn,epsgc
 194       common/lept2/nstmax,prob(99),icbac(99,2),icfor(99,2)
 195       common/lept3/iolept,igampr,idisco
 196       common/ebin/noebin,engmin,engmax,nrebin,iologe,iologl
 197       common/cnsta/pi,pii,hquer,prom,piom,ainfin
 198       common/versn/iversn,iverso
 199       common/accum/imsg,jerr(mxjerr),ntevt,nrevt,naevt,nrstr,nrptl
 200       common/accum2/nglacc
 201       common/cptlu/nptlu /cnrclu/nrclu
 202       common/drop6/tecm,volu
 203       common/metr5/iterma,iternc,iterpr,iterpl,iozinc,iozevt
 204       common/metr6/epsr,keepr
 205       common/drop5/keu,ked,kes,kec,keb,ket
 206       double precision seedi,seedj,seedc
 207       common/cseed/seedi,seedj,seedc
 208       common/cjintc/clust(mxtau,mxvol,mxeps)
 209       common/cjintd/volsum(mxtau),vo2sum(mxtau),nclsum(mxtau)
 210       common/ciutot/iutotc,iutote
 211       common/copen/nopen,nopenr
 212       common/kopen/kchopen,khiopen,kdtopen,kcpopen,klgopen,knxopen
 213       character*6 xvaria,yvaria
 214       common/vana1/xvaria,yvaria,normal,xminim,xmaxim,nrbins,hisfac
 215       common/vana3/iologb,iocnxb
 216       parameter(mxnody=200)
 217       common/nodcy/nrnody,nody(mxnody)
 218       character*20 subpri
 219       common/prnt2/nrpri,subpri(mxpri),ishpri(mxpri)
 220       common/prnt3/ishevt,ixtau,iwseed,ixgeometry
 221       parameter (mxcnt=20)
 222       common/vana4/ar(mxbins,5),ary(mxbins,mxcnt),ardy(mxbins,mxcnt)
 223      &,ionoerr
 224       common/xpars/xpar1,xpar2,xpar3,xpar4,xpar5,xpar6,xpar7,xpar8
 225       common/khist/khisto
 226       common/ctcor/nctcor/ccttim/ncttim
 227       common/densi/kdensi(matau,nzeta,mxcoox,mxcooy),tauv(matau)
 228       common/cjinti/iorsce,iorsdf,iorshh
 229       common/camim/amimfs,amimel
 230       common/craddf/scr,scs,hacore
 231       common/ckoll/iokoll
 232       common/cncnt/ncnt  /cicnt/inicnt /cnemsi/nemsi
 233       common/ems1/iemspl,iemsct,gfactor,gwidth
 234       common/chadron/amproj,amtarg,ypjtl,yhaha,pnullx
 235       common/vana5/xshift,etacut
 236       double precision rnucl
 237       common/nucl5/rnucl(mxnucl,2)
 238       common/nucl4/nrnucl(2),drnucl(2),rnuclo(mxnucl,2),ionudi
 239       common/sig/xsig(7),xpom(7)
 240       common/metr7/ktnbod
 241       common/hadr3/iregge,isopom,ishpom,iscreen,nprmax,inueff,irmdrop
 242       common/hadr5/sigtot,sigcut,sigela,sloela,sigsd,sigine,sigdif
 243      &,sigineaa
 244       common/hadr6/intpol,isigma,iomega,isetcs
 245       common/hadr4/alppom,slopom,gamhad(4),r2had(4),chad(4),wdiff(4)
 246      &            ,gamtil,facdif,r2hads(4),gamhads(4),slopoms,isplit
 247       common/hadr42/gamhadsi(4)
 248       common/hadr7/alpreg,sloreg,gamreg,r2reg,ptdiff,ptsend,xmaxremn
 249      &,xmaxdiff
 250       common/hadr8/alpqua,alppar,alpsea,alpval,alpdiq,alplea(4),alpdif
 251       common/hadr14/alpndi,alpdi,ptsendi,zrminc,zrmmax,ptsemx,irzptn
 252       common/hadr15/zbcut,zopinc,zopmax,zipinc,zipmax,zoeinc,zoemax
 253       common/hadr16/zidinc,zidmax,zodinc,zodmax,zosinc,zosmax
 254       common/hadr17/edmaxi,epmaxi
 255       common/hadr9/ammsqq,ammsqd,ammsdd,cumpom,rexndi(4),rexdif(4)
 256      &            ,reminv,rexres,zbarfl,rexndf
 257       common/had10/iclpro,icltar,iclegy
 258       common/had11/iclpro1,iclpro2,icltar1,icltar2,iclegy1,iclegy2
 259       common/had12/egylow,egyfac
 260       common/had13/amdrmax,alpdro(3)
 261       common/had14/alpcoso,alpcose,betcoso,betcose
 262       common/emsx1/accept,reject
 263       common/xems1/iemspr,iemspm,iemspx,iemsrx,iemspu,iemsi2,iemspbx
 264       common/xems2/iemsse,iemsi1,iemsb,iemsbg,ioems,iemsdr
 265       common/xspatim/ispacetime
 266       character*80 hisarg
 267       common/chisarg/ihisarg,hisarg(2,mxhisarg)
 268       common /psar10/difnuc(mamxx),radnuc(mamxx)
 269       parameter (mxbarray=100)
 270       common/cbarray/barray(mxbarray),nbarray
 271       common/nxsair/airznxs(3),airanxs(3),airwnxs(3)
 272      &             ,airavznxs,airavanxs
 273       common/mod2incs/qgsincs
 274       common/mod3incs/gheincs
 275       common/mod4incs/pytincs
 276       common/mod5incs/hijincs
 277       common/mod6incs/sibincs
 278       common/mod7incs/qgsIIincs
 279       common/testpom/antot,ansh,ansf,pp4max,pp4ini,andropl,anstrg0
 280      *,anshf,ansff,antotf
 281      *,anstrg1,anreso0,anreso1,anghadr,antotre
 282       common/cdiff/anintdiff,anintsdif,anintine
 283      *,sigineex,sigdifex,sigsdex
 284       common/cepszer/epszero,alpff(nclha),betff(2)
 285       common/cgss/tgss(7,7),wgss(7,7)
 286
 287       common/Dparams/alpDs(  idxD0:idxD, nclegy, nclha,nclha),
 288      &              alpDps( idxD0:idxD, nclegy, nclha,nclha),
 289      &              alpDpps(idxD0:idxD, nclegy, nclha,nclha),
 290      &              betDs(  idxD0:idxD, nclegy, nclha,nclha),
 291      &              betDps( idxD0:idxD, nclegy, nclha,nclha),
 292      &              betDpps(idxD0:idxD, nclegy, nclha,nclha),
 293      &              gamDs(  idxD0:idxD, nclegy, nclha,nclha),
 294      &              delDs(  idxD0:idxD, nclegy, nclha,nclha)
 295       common/Dparam/alpD(  idxD0:idxD1, nclha, nclha),
 296      &              alpDp( idxD0:idxD1, nclha, nclha),
 297      &              alpDpp(idxD0:idxD1, nclha, nclha),
 298      &              betD(  idxD0:idxD1, nclha, nclha),
 299      &              betDp( idxD0:idxD1, nclha, nclha),
 300      &              betDpp(idxD0:idxD1, nclha, nclha),
 301      &              gamD(  idxD0:idxD1, nclha, nclha),
 302      &              delD(  idxD0:idxD1, nclha, nclha),
 303      &     idxDmin,bmxdif(nclha, nclha),bkmxndif
 304       double precision alpUni,betUni,betpUni,zpUni,ztUni,betfom
 305       common/DparUni/alpUni(  idxD0:idxD1,2),
 306      &               betUni(  idxD0:idxD1,2),
 307      &               betpUni(idxD0:idxD1,2),zpUni,ztUni,
 308      &               betfom,alpfom,alpfomi,gamfom