Warning, /GeneratorInterface/ReggeGribovPartonMCInterface/src/epos.inc is written in an unsupported language. File is not indexed.
0001 c---------------------------------------------------------------------------
0002 c dimensions
0003 c---------------------------------------------------------------------------
0004
0005 integer mmry,mxptl,nmxhep,myptl,nzeta,nflav,mxstr,mystr,mxtau
0006 * ,mxtrig,mxpri,mxbins,matau,mxnucl,mxhisarg,idxD0,idxD1
0007 * ,idxD,nclha,nclegy,mamxx,mxjerr,mxvol,mxeps,mxidco,mxcoox
0008 * ,mxcooy
0009 parameter (mmry=1) !memory saving factor
0010
0011 parameter (mxptl=200000/mmry) !max nr of particles in epos ptl list
0012 parameter (nmxhep=99900) !max nr of particles in hep ptl list
0013 parameter (myptl=1000) !max nr of droplets in epos ptl list
0014 parameter (nzeta=60) !max nr of zeta bins for droplets
0015 parameter (nflav=6) !max nr of flavors
0016 parameter (mxstr=20000/mmry) !max nr of strings in epos string list
0017 parameter (mystr=20000/mmry)
0018 parameter (mxtau=4,mxvol=10,mxeps=16)
0019 parameter (mxtrig=99,mxidco=99)
0020 parameter (mxpri=200)
0021 parameter (mxbins=10000)
0022 parameter (matau=10,mxcoox=40,mxcooy=10)
0023 parameter (mxnucl=20)
0024 parameter (mxhisarg=100)
0025 parameter (idxD0=0,idxD1=2,idxD=1,nclha=4,nclegy=100)
0026 parameter (mamxx=250)
0027 parameter (mxjerr=10)
0028
0029 c---------------------------------------------------------------------------
0030 c epos event common block
0031 c---------------------------------------------------------------------------
0032
0033 real phievt,bimevt,pmxevt,egyevt
0034 *,xbjevt,qsqevt,zppevt,zptevt
0035 integer nevt,kolevt,koievt,kohevt,npjevt
0036 *,ntgevt,npnevt,nppevt,ntnevt,ntpevt,jpnevt,jppevt,jtnevt,jtpevt
0037 *,nglevt,minfra,maxfra
0038 common/cevt/phievt,nevt,bimevt,kolevt,koievt,pmxevt,egyevt,npjevt
0039 *,ntgevt,npnevt,nppevt,ntnevt,ntpevt,jpnevt,jppevt,jtnevt,jtpevt
0040 *,xbjevt,qsqevt,nglevt,zppevt,zptevt,minfra,maxfra,kohevt
0041 real rglevt,sglevt,eglevt,fglevt,typevt
0042 integer ng1evt,ng2evt,ikoevt
0043 common/c2evt/ng1evt,ng2evt,rglevt,sglevt,eglevt,fglevt,ikoevt
0044 *,typevt !qgsjet03.f and qgsjetII.f and phojet.f
0045
0046 c nevt .......... error code. 1=valid event, 0=invalid event
0047 c bimevt ........ absolute value of impact parameter
0048 c phievt ........ angle of impact parameter
0049 c kolevt ........ number of collisions
0050 c koievt ........ number of inelastic collisions
0051 c kohevt ........ number of hard collisions
0052 c pmxevt ........ reference momentum
0053 c egyevt ........ pp cm energy (hadron) or string energy (lepton)
0054 c npjevt ........ number of primary projectile participants
0055 c ntgevt ........ number of primary target participants
0056 c npnevt ........ number of primary projectile neutron spectators
0057 c nppevt ........ number of primary projectile proton spectators
0058 c ntnevt ........ number of primary target neutron spectators
0059 c ntpevt ........ number of primary target proton spectators
0060 c jpnevt ........ number of absolute projectile neutron spectators
0061 c jppevt ........ number of absolute projectile proton spectators
0062 c jtnevt ........ number of absolute target neutron spectators
0063 c jtpevt ........ number of absolute target proton spectators
0064 c xbjevt ........ bjorken x for dis
0065 c qsqevt ........ q**2 for dis
0066 c sigtot ........ total cross section
0067 c nglevt ........ number of collisions acc to Glauber
0068 c zppevt ........ average Z-parton-proj
0069 c zptevt ........ average Z-parton-targ
0070 c ng1evt ........ number of Glauber participants with at least one IAs
0071 c ng2evt ........ number of Glauber participants with at least two IAs
0072 c ikoevt ........ number of elementary parton-parton scatterings
0073 c typevt ........ type of event (1=Non Diff, 2=Double Diff, 3=Central Diff, 4=AB->XB, -4=AB->AX)
0074
0075 c---------------------------------------------------------------------------
0076 c epos particle list common block
0077 c---------------------------------------------------------------------------
0078
0079 real pptl,tivptl,xorptl
0080 integer nptl,iorptl,idptl,istptl,ifrptl,jorptl,ibptl,ityptl
0081 common/cptl/nptl,pptl(5,mxptl),iorptl(mxptl),idptl(mxptl)
0082 *,istptl(mxptl),tivptl(2,mxptl),ifrptl(2,mxptl),jorptl(mxptl)
0083 *,xorptl(4,mxptl),ibptl(4,mxptl),ityptl(mxptl)
0084 real ekievt
0085 integer itsptl
0086 common/c1ptl/ekievt,itsptl(mxptl)
0087
0088 c nptl .......... current particle index (=number of ptls stored)
0089 c idptl(i) ...... particle id
0090 c pptl(1,i) ..... x-component of particle momentum
0091 c pptl(2,i) ..... y-component of particle momentum
0092 c pptl(3,i) ..... z-component of particle momentum
0093 c pptl(4,i) ..... particle energy
0094 c pptl(5,i) ..... particle mass
0095 c iorptl(i) ..... particle number of father (if .le. 0 : no father)
0096 c jorptl(i) ..... particle number of mother (if .le. 0 : no mother)
0097 c istptl(i) ..... status: 40 and 41 : Remnant
0098 c 30 and 31 : Pomeron
0099 c 20 and 21 : Parton
0100 c 10 and 11 : Droplet
0101 c 00 and 01 : Particle
0102 c last digit = 0 : last generation
0103 c last digit = 1 : not last generation
0104 c xorptl(1,i) ... x-component of formation point
0105 c xorptl(2,i) ... y-component of formation point
0106 c xorptl(3,i) ... z-component of formation point
0107 c xorptl(4,i) ... formation time
0108 c tivptl(1,i) ... formation time (always in the pp-cms!)
0109 c tivptl(2,i) ... destruction time (always in the pp-cms!)
0110 c ityptl(i) .... type of particles origin:
0111 c 10-19: target
0112 c 20-29: soft Pom
0113 c 30-39: hard Pom
0114 c 40-49: projectile
0115 c 50: string, droplet
0116 c itsptl(i) ..... string type of particles origin (if string)
0117
0118 real radptl
0119 integer iaaptl
0120 common/c2ptl/iaaptl(mxptl),radptl(mxptl)
0121 real desptl,dezptl
0122 common/c3ptl/desptl(mxptl),dezptl(mxptl)
0123 integer nptlbd
0124 common/c4ptl/nptlbd
0125 real rinptl,vrad
0126 integer inbxxx
0127 common/c6ptl/rinptl(mxptl),vrad,inbxxx
0128 real qsqptl,zpaptl
0129 common/c8ptl/qsqptl(mxptl),zpaptl(2,mxptl)
0130
0131
0132 c---------------------------------------------------------------------------
0133 c hep standard event commonblock.
0134 c---------------------------------------------------------------------------
0135
0136 double precision phep,vhep
0137 integer nevhep,nhep,isthep,idhep,jmohep,jdahep
0138
0139 common/hepcom/nevhep,nhep,isthep(nmxhep),idhep(nmxhep),
0140 &jmohep(2,nmxhep),jdahep(2,nmxhep),phep(5,nmxhep),vhep(4,nmxhep)
0141
0142 c---------------------------------------------------------------------------
0143 c
0144 c nevhep - event number
0145 c nhep - number of entries in the event record
0146 c
0147 c isthep(i) - status code
0148 c idhep(i) - particle id (particle data group standard)
0149 c
0150 c jmohep(1,i) - position of mother particle in list
0151 c jmohep(2,i) - position of second mother particle in list
0152 c jdahep(1,i) - position of first daughter in list
0153 c jdahep(2,i) - position of first daughter in list
0154 c
0155 c phep(1,i) - p_x momentum in gev/c
0156 c phep(2,i) - p_y momentum in gev/c
0157 c phep(3,i) - p_z momentum in gev/c
0158 c phep(4,i) - energy in gev
0159 c phep(5,i) - mass in gev/c**2
0160 c
0161 c vhep(1,i) - x position of production vertex in mm
0162 c vhep(2,i) - y position of production vertex in mm
0163 c vhep(3,i) - z position of production vertex in mm
0164 c vhep(4,i) - time of production in mm/c
0165 c
0166 c (note: 1 mm = 10^-12 fm = 5.07 10^-12 1/gev)
0167
0168 c------------------------------------------------------------------------
0169 c Parameters set in sr aaset and variables to communicate between moduls
0170 c------------------------------------------------------------------------
0171
0172 integer ifop,ifmt,ifch,ifcx,ifhi,ifdt,ifcp,ifdr
0173 common/files/ifop,ifmt,ifch,ifcx,ifhi,ifdt,ifcp,ifdr
0174 integer ifin
0175 common/files2/ifin
0176 character*500 fnch,fnhi,fndt,fnii,fnid,fnie,fnrj,fnmt
0177 * ,fngrv,fncp,fnnx,fncs,fndr,fnhpf
0178 common/fname/ fnch, fnhi, fndt, fnii, fnid, fnie, fnrj, fnmt
0179 * ,fngrv,fncp,fnnx,fncs,fndr,fnhpf
0180 integer nfnch,nfnhi,nfndt,nfnii,nfnid,nfnie,nfnrj,nfnmt
0181 *,nfngrv,nfncp,nfnnx,nfncs,nfndr,nfnhpf
0182 common/nfname/nfnch,nfnhi,nfndt,nfnii,nfnid,nfnie,nfnrj,nfnmt
0183 *,nfngrv,nfncp,nfnnx,nfncs,nfndr,nfnhpf
0184 character*500 fnin
0185 common/fname2/ fnin
0186 integer nfnin
0187 common/nfname2/nfnin
0188 real delvol,deleps,dlzeta,etafac,facnuc,taurea,epscri
0189 common/resc2/delvol,deleps,dlzeta,etafac,facnuc,taurea,epscri(3)
0190 real epsdfm
0191 common/resc4/epsdfm
0192 character*3 hydt
0193 common/hydr1/hydt
0194 real pud,pmqu,pmqd,pmqs,pmqc,pmqq
0195 integer ndecay,maxres
0196 common/frag1/ndecay,maxres,pud,pmqu,pmqd,pmqs,pmqc,pmqq
0197 real pdiqua,delrex,ptfraqq,ptfra,ptfrasr
0198 integer ioptf
0199 common/frag2/pdiqua,delrex,ptfraqq,ptfra,ptfrasr,ioptf
0200 real aouni,pbreak,fkappa,strcut,diqcut,fkappag,pbreakg
0201 *,zetacut
0202 integer itflav
0203 common/frag3/aouni,pbreak,fkappa,itflav,strcut,diqcut
0204 *,fkappag,pbreakg,zetacut
0205 real difud,difus,difuc,pudd,puds,pudc,difuuu,difuud
0206 *,difuus,difuuc,difudd,difuds,difudc,difuss,difusc,difucc
0207 integer nrflav
0208 common/frag4/difud,difus,difuc,pudd,puds,pudc,difuuu,difuud
0209 *,difuus,difuuc,difudd,difuds,difudc,difuss,difusc,difucc,nrflav
0210 real qmass
0211 integer isospin
0212 common/frag5/qmass(0:6),isospin(0:6)
0213 real pnll,ptq,exmass,cutmss,wproj,wtarg
0214 common/hadr1/pnll,ptq,exmass,cutmss,wproj,wtarg
0215 real rstrau,rstrad,rstras,rstrac,rstrasi
0216 common/hadr10/rstrau(4),rstrad(4),rstras(4),rstrac(4),rstrasi
0217 real wgtqqq,wgtval,wgtsea,wgtdiq
0218 common/wgtqrk/wgtqqq(4),wgtval,wgtsea,wgtdiq
0219 double precision timeini,timefin
0220 common/time1/timeini,timefin
0221 integer iotst1,iotst2,iotst3,iotst4
0222 common/ciotst/iotst1,iotst2,iotst3,iotst4
0223 real taumin,deltau,factau,numtau,amsiac,amprif
0224 common/resc1/taumin,deltau,factau,numtau,amsiac,amprif
0225 real dscale,cepara,delamf,deuamf
0226 integer iceopt
0227 common/resc3/dscale,cepara,iceopt,delamf,deuamf
0228 integer ispherio,icotabm,icotabr,icocore
0229 common/sprio/ispherio,icotabm,icotabr,icocore
0230 integer iurqmd
0231 common/urqm1/iurqmd
0232 real cutico,dssico
0233 common/incon/cutico,dssico
0234 real gaumx
0235 integer istore,istmax,irescl,ntrymx,nclean,iopdg,ioidch
0236 common/othe1/istore,istmax,gaumx,irescl,ntrymx,nclean,iopdg,ioidch
0237 integer ifrade,iframe,idecay,jdecay,iremn
0238 common/othe2/ifrade,iframe,idecay,jdecay,iremn
0239 integer jframe,kframe
0240 common/othe3/jframe,kframe
0241 integer iselect
0242 common/othe4/iselect
0243 real taumx,sigj
0244 integer jpsi,jpsifi,nsttau,ijphis,ijtauan
0245 common/jpsif/jpsi,jpsifi,taumx,nsttau,sigj,ijphis,ijtauan
0246 real themas
0247 integer iopenu
0248 common/strlt/iopenu,themas
0249 integer iappl,model
0250 common/appli/iappl,model
0251 integer nevent,nfull,nfreeze,ninicon
0252 common/events/nevent,nfull,nfreeze,ninicon
0253 real egymin,egymax,elab,ecms,ekin
0254 common/enrgy/egymin,egymax,elab,ecms,ekin
0255 integer iprmpt,ish,ishsub,irandm,irewch,iecho,modsho,idensi
0256 common/prnt1/iprmpt,ish,ishsub,irandm,irewch,iecho,modsho,idensi
0257 real engy,elepti,elepto,angmue
0258 integer icinpu
0259 common/lept1/engy,elepti,elepto,angmue,icinpu
0260 real core,fctrmx
0261 integer laproj,maproj,latarg,matarg
0262 common/nucl1/laproj,maproj,latarg,matarg,core,fctrmx
0263 real bmaxim,bminim,phimax,phimin
0264 common/nucl2/bmaxim,bminim,phimax,phimin
0265 real ymximi,wtimet,wtimei,wtimea
0266 integer imihis,iclhis,iwtime
0267 common/wana1/ymximi,imihis,iclhis,iwtime,wtimet,wtimei,wtimea
0268 real wtmini,wtstep
0269 integer isphis,ispall,iwcent,iana,nbdky
0270 common/wana2/isphis,ispall,wtmini,wtstep,iwcent,iana,nbdky
0271 real asuhax,asuhay
0272 common/drop4/asuhax(7),asuhay(7)
0273 real grigam,grirsq,gridel,grislo,gricel,sigppi,sigppd
0274 common/gribo/grigam,grirsq,gridel,grislo,gricel,sigppi,sigppd
0275 real bag4rt,dezzer,amuseg,taunll,yradmx,facts,factb,factq
0276 common/drop3/bag4rt,dezzer,amuseg,taunll,yradmx,facts,factb,factq
0277 real rcoll,ylongmx,facecc,yradpp,yradmi,yrmaxi,fradflii
0278 integer nsegsu,nsegce
0279 common/drop2/rcoll,ylongmx,nsegsu,nsegce,facecc,yradpp,yradmi
0280 * ,yrmaxi,fradflii
0281 real ptclu,yradpi,yradpx,fploss,fvisco,fplmin
0282 integer ioclude,iocluin,ioquen,iohole,kigrid
0283 common/drop7/ptclu,yradpi,yradpx,fploss,fvisco,fplmin,ioclude
0284 * ,iocluin,ioquen,iohole,kigrid
0285 real fsgrid
0286 common/drop8/fsgrid
0287 integer iospec,iocova,iopair,iozero,ioflac,iomom
0288 common/metr1/iospec,iocova,iopair,iozero,ioflac,iomom
0289 integer nadd,iograc,iocite,ioceau,iociau
0290 common/metr2/nadd,iograc,iocite,ioceau,iociau
0291 integer iomodl,idproj,idtarg
0292 real wexcit
0293 common/hadr2/iomodl,idproj,idtarg,wexcit
0294 real rexdifi,rexndii
0295 integer idprojin,idtargin,irdmpr,isoproj,isotarg
0296 common/hadr25/idprojin,idtargin,rexdifi(4),rexndii(4),irdmpr,
0297 * isoproj,isotarg
0298 integer iostat,ioinco,ionlat,ioobsv,iosngl,iorejz,iompar
0299 common/metr3/iostat,ioinco,ionlat,ioobsv,iosngl,iorejz,iompar
0300 integer ioinfl,ioinct,iowidn
0301 real epsgc
0302 common/metr4/ioinfl,ioinct,iowidn,epsgc
0303 real prob
0304 integer nstmax,icbac,icfor
0305 common/lept2/nstmax,prob(99),icbac(99,2),icfor(99,2)
0306 integer iolept,igampr,idisco
0307 common/lept3/iolept,igampr,idisco
0308 real engmin,engmax
0309 integer noebin,nrebin,iologe,iologl
0310 common/ebin/noebin,engmin,engmax,nrebin,iologe,iologl
0311 real pi,pii,hquer,prom,piom,ainfin
0312 common/cnsta/pi,pii,hquer,prom,piom,ainfin
0313 integer iversn,iverso
0314 common/versn/iversn,iverso
0315 integer imsg,jerr,ntevt,nrevt,naevt,nrstr,nrptl
0316 common/accum/imsg,jerr(mxjerr),ntevt,nrevt,naevt,nrstr,nrptl
0317 integer nglacc
0318 common/accum2/nglacc
0319 integer nptlu,nrclu
0320 common/cptlu/nptlu /cnrclu/nrclu
0321 real tecm,volu
0322 common/drop6/tecm,volu
0323 integer iterma,iternc,iterpr,iterpl,iozinc,iozevt
0324 common/metr5/iterma,iternc,iterpr,iterpl,iozinc,iozevt
0325 real epsr
0326 integer keepr
0327 common/metr6/epsr,keepr
0328 integer keu,ked,kes,kec,keb,ket
0329 common/drop5/keu,ked,kes,kec,keb,ket
0330 double precision seedi,seedj,seedj2,seedc
0331 integer iseqini,iseqsim
0332 common/cseed/seedi,seedj,seedj2,seedc,iseqini,iseqsim
0333 real clust
0334 common/cjintc/clust(mxtau,mxvol,mxeps)
0335 real volsum,vo2sum
0336 integer nclsum
0337 common/cjintd/volsum(mxtau),vo2sum(mxtau),nclsum(mxtau)
0338 integer iutotc,iutote
0339 common/ciutot/iutotc,iutote
0340 integer nopen,nopenr
0341 common/copen/nopen,nopenr
0342 integer kchopen,khiopen,kdtopen,kcpopen,klgopen,knxopen
0343 common/kopen/kchopen,khiopen,kdtopen,kcpopen,klgopen,knxopen
0344 character*6 xvaria,yvaria
0345 real xminim,xmaxim,hisfac
0346 integer normal,nrbins
0347 common/vana1/xvaria,yvaria,normal,xminim,xmaxim,nrbins,hisfac
0348 integer iologb,iocnxb
0349 common/vana3/iologb,iocnxb
0350 integer mxnody,nrnody,nody
0351 parameter(mxnody=200)
0352 common/nodcy/nrnody,nody(mxnody)
0353 real ctaumin
0354 common/ctdcy/ctaumin
0355 character*20 subpri
0356 integer nrpri,ishpri
0357 common/prnt2/nrpri,subpri(mxpri),ishpri(mxpri)
0358 integer ishevt,ixtau,iwseed,jwseed,ixgeometry
0359 common/prnt3/ishevt,ixtau,iwseed,jwseed,ixgeometry
0360 integer mxcnt,ionoerr
0361 real ar,ary,ardy
0362 parameter (mxcnt=20)
0363 common/vana4/ar(mxbins,5),ary(mxbins,mxcnt),ardy(mxbins,mxcnt)
0364 *,ionoerr
0365 real xpar1,xpar2,xpar3,xpar4,xpar5,xpar6,xpar7,xpar8
0366 *,xpar99
0367 common/xpars/xpar1,xpar2,xpar3,xpar4,xpar5,xpar6,xpar7,xpar8
0368 *,xpar99
0369 integer khisto
0370 common/khist/khisto
0371 integer nctcor,ncttim
0372 common/ctcor/nctcor/ccttim/ncttim
0373 integer kdensi
0374 real tauv
0375 common/densi/kdensi(matau,nzeta,mxcoox,mxcooy),tauv(matau)
0376 integer iorsce,iorsdf,iorshh,ionudi
0377 common/cjinti/iorsce,iorsdf,iorshh,ionudi
0378 real amimfs,amimel
0379 common/camim/amimfs,amimel
0380 real scr,scs,hacore
0381 common/craddf/scr,scs,hacore
0382 integer iokoll
0383 common/ckoll/iokoll
0384 integer ncnt,inicnt,nemsi
0385 common/cncnt/ncnt /cicnt/inicnt /cnemsi/nemsi
0386 real gfactor,gwidth
0387 integer iemspl,iemsct
0388 common/ems1/iemspl,iemsct,gfactor,gwidth
0389 real amproj,amtarg,ypjtl,yhaha,pnullx
0390 common/chadron/amproj,amtarg,ypjtl,yhaha,pnullx
0391 real xshift,etacut
0392 common/vana5/xshift,etacut
0393 double precision rnucl
0394 real bnucl,xbtot
0395 integer ixbDens
0396 common/nucl6/infragm
0397 integer infragm
0398 common/nucl5/rnucl(mxnucl,2),bnucl(mxnucl,4),xbtot(4),ixbDens
0399 real drnucl,rnuclo
0400 integer nrnucl
0401 common/nucl4/nrnucl(2),drnucl(2),rnuclo(mxnucl,2)
0402 real xsig,xpom
0403 common/sig/xsig(7),xpom(7)
0404 integer ktnbod
0405 common/metr7/ktnbod
0406 integer iregge,isopom,ishpom,iscreen,nprmax,inueff,irmdrop
0407 common/hadr3/iregge,isopom,ishpom,iscreen,nprmax,inueff,irmdrop
0408 real sigtot,sigcut,sigela,sloela,sigsd,sigine,sigdif
0409 *,sigineaa,sigtotaa,sigelaaa,sigcutaa,sigdd
0410 common/hadr5/sigtot,sigcut,sigela,sloela,sigsd,sigine,sigdif
0411 *,sigineaa,sigtotaa,sigelaaa,sigcutaa,sigdd
0412 integer intpol,isigma,iomega,isetcs
0413 common/hadr6/intpol,isigma,iomega,isetcs
0414 real alppom,slopom,gamhad,r2had,chad,wdiff
0415 *,gamtil,facdif,facmc,r2hads,gamhads,slopoms
0416 integer isplit
0417 common/hadr4/alppom,slopom,gamhad(4),r2had(4),chad(4),wdiff(4)
0418 & ,gamtil,facdif,facmc,r2hads(4),gamhads(4),slopoms,isplit
0419 real gamhadsi
0420 common/hadr42/gamhadsi(4)
0421 real alpreg,sloreg,gamreg,r2reg,ptdiff,ptsend,xminremn
0422 &,xmindiff,ptsecu
0423 common/hadr7/alpreg,sloreg,gamreg,r2reg,ptdiff,ptsend,xminremn
0424 &,xmindiff,ptsecu
0425 real alpqua,alppar,alpsea,alpval,alpdiq,alplea,alpdif
0426 common/hadr8/alpqua,alppar,alpsea,alpval,alpdiq,alplea(4),alpdif
0427 real alpndi,alpdi,ptsendi,zdrinc,zmsinc,ptsems
0428 integer irzptn
0429 common/hadr14/alpndi,alpdi,ptsendi,zdrinc,zmsinc,ptsems,irzptn
0430 real zbcut,zopinc,zipinc,zoeinc,xmxrem
0431 common/hadr15/zbcut,zopinc,zipinc,zoeinc,xmxrem
0432 real fkainc,fkamax,zodinc,zbrmax,zdfinc,xzcut,ptvpom
0433 common/hadr16/fkainc,fkamax,zodinc,zbrmax,zdfinc,xzcut,ptvpom
0434 real edmaxi,epmaxi
0435 common/hadr17/edmaxi,epmaxi
0436 real ammsqq,ammsqd,ammsdd,cumpom,rexndi,rexdif
0437 * ,reminv,rexpdif,rexres,zrminc,rexndf
0438 common/hadr9/ammsqq,ammsqd,ammsdd,cumpom,rexndi(4),rexdif(4)
0439 * ,reminv,rexpdif(4),rexres(4),zrminc,rexndf
0440 integer iclpro,icltar,iclegy
0441 common/had10/iclpro,icltar,iclegy
0442 integer iclpro1,iclpro2,icltar1,icltar2,iclegy1,iclegy2
0443 common/had11/iclpro1,iclpro2,icltar1,icltar2,iclegy1,iclegy2
0444 real egylow,egyfac
0445 common/had12/egylow,egyfac
0446 real amdrmax,amdrmin,alpdro
0447 common/had13/amdrmax,amdrmin,alpdro(3)
0448 real alpcoso,alpcose,betcoso,betcose
0449 common/had14/alpcoso,alpcose,betcoso,betcose
0450 real accept,reject
0451 common/emsx1/accept,reject
0452 integer iemspr,iemspm,iemspx,iemsrx,iemspu,iemsi2,iemspbx
0453 common/xems1/iemspr,iemspm,iemspx,iemsrx,iemspu,iemsi2,iemspbx
0454 integer iemsse,iemsi1,iemsb,iemsbg,ioems,iemsdr
0455 common/xems2/iemsse,iemsi1,iemsb,iemsbg,ioems,iemsdr
0456 integer ispacetime
0457 common/xspatim/ispacetime
0458 real difnuc,radnuc
0459 common /psar10/difnuc(mamxx),radnuc(mamxx)
0460 integer mxbarray,nbarray
0461 real barray
0462 parameter (mxbarray=100)
0463 common/cbarray/barray(mxbarray),nbarray
0464 real airznxs,airanxs,airwnxs
0465 * ,airavznxs,airavanxs
0466 common/nxsair/airznxs(3),airanxs(3),airwnxs(3)
0467 * ,airavznxs,airavanxs
0468 real qgsincs
0469 common/mod2incs/qgsincs
0470 real gheincs
0471 common/mod3incs/gheincs
0472 real pytincs
0473 common/mod4incs/pytincs
0474 real hijincs
0475 common/mod5incs/hijincs
0476 real sibincs
0477 common/mod6incs/sibincs
0478 real qgsIIincs
0479 common/mod7incs/qgsIIincs
0480 real phoincs
0481 common/mod8incs/phoincs
0482 real fluincs
0483 common/mod9incs/fluincs
0484 real urqincs
0485 common/mod10incs/urqincs
0486 real antot,ansh,ansf,pp4max,pp4ini,andropl,anstrg0
0487 *,anshf,ansff,antotf
0488 *,anstrg1,anreso0,anreso1,anghadr,antotre
0489 common/testpom/antot,ansh,ansf,pp4max,pp4ini,andropl,anstrg0
0490 *,anshf,ansff,antotf
0491 *,anstrg1,anreso0,anreso1,anghadr,antotre
0492 real anintdiff,anintsdif,anintine
0493 *,sigineex,sigdifex,sigsdex
0494 common/cdiff/anintdiff,anintsdif,anintine
0495 *,sigineex,sigdifex,sigsdex
0496 real epszero,alpff,betff
0497 common/cepszer/epszero,alpff(nclha),betff(2)
0498 real tgss,wgss
0499 common/cgss/tgss(7,7),wgss(7,7)
0500
0501 real alpDs,alpDps,alpDpps,betDs,betDps,betDpps,gamDs,delDs
0502 * ,alpD,alpDp,alpDpp,betD,betDp,betDpp,gamD,delD
0503 common/Dparams/alpDs( idxD0:idxD, nclegy, nclha,nclha),
0504 * alpDps( idxD0:idxD, nclegy, nclha,nclha),
0505 * alpDpps(idxD0:idxD, nclegy, nclha,nclha),
0506 * betDs( idxD0:idxD, nclegy, nclha,nclha),
0507 * betDps( idxD0:idxD, nclegy, nclha,nclha),
0508 * betDpps(idxD0:idxD, nclegy, nclha,nclha),
0509 * gamDs( idxD0:idxD, nclegy, nclha,nclha),
0510 * delDs( idxD0:idxD, nclegy, nclha,nclha)
0511 real bmxdif,bkmxndif
0512 integer idxDmin
0513 common/Dparam/alpD( idxD0:idxD1, nclha, nclha),
0514 * alpDp( idxD0:idxD1, nclha, nclha),
0515 * alpDpp(idxD0:idxD1, nclha, nclha),
0516 * betD( idxD0:idxD1, nclha, nclha),
0517 * betDp( idxD0:idxD1, nclha, nclha),
0518 * betDpp(idxD0:idxD1, nclha, nclha),
0519 * gamD( idxD0:idxD1, nclha, nclha),
0520 * delD( idxD0:idxD1, nclha, nclha),
0521 * idxDmin,bmxdif(nclha, nclha),bkmxndif
0522 double precision alpUni,betUni,betpUni,fpUni,ftUni,betfom
0523 * ,epspUni,epstUni,zzpUni,zztUni
0524 real alpfom,alpfomi,gamfom
0525 common/DparUni/alpUni( idxD0:idxD1,2),
0526 * betUni( idxD0:idxD1,2),
0527 * betpUni(idxD0:idxD1,2),fpUni,ftUni,
0528 * epspUni(idxD0:idxD1),zzpUni,
0529 * epstUni(idxD0:idxD1),zztUni,
0530 * betfom,alpfom,alpfomi,gamfom
0531 integer idlead,ilprtg
0532 common/crvar/idlead,ilprtg
0533 integer iLHC
0534 common/LHCtune/iLHC
