]>
Commit | Line | Data |
---|---|---|
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 |