| 9ef1c2d9 |
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 |
git://git.uio.no / u / mrichter / AliRoot.git / blame