3       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
    9       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
  115       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
  116      +iibar(210),k1(210),k2(210)
 
  121       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
  122       COMMON /nuccc/   jt,jtz,jp,jpz,jjproj,jbproj,jjtarg,jbtarg
 
  124       COMMON /cmhico/ cmhis
 
  126       COMMON /resona/ ireso
 
  128       COMMON /trafop/ gamp,bgamp,betp
 
  130       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
  132       COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
 
  133      +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
 
  134      +irvs14, irvv11,irvv12,irvv13,irvv14
 
  136       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
  138       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
  139      +ipadis,ishmal,lpauli
 
  141       COMMON /dnun/   nn,np,
nt 
  143       COMMON /nshmak/ nnshma,npshma,ntshma,nshmac,nshma2
 
  145       COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
 
  146      *              bsite(0:1,200),nstatb,nsiteb
 
  148       COMMON /seaqxx/ seaqx,seaqxn 
 
  149       COMMON /diffra/ isingd,idiftp,ioudif,iflagd
 
  150       COMMON /evflag/ numev
 
  151       COMMON /neutyy/ neutyp,neudec
 
  152       COMMON /fluctu/ifluct
 
  153        COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
 
  154      *idzre(3),izdre(3),idiqrz(7)
 
  155        COMMON /intneu/ndzsu,nzdsu
 
  159        COMMON /bufues/ bnnvv,bnnss,bnnsv,bnnvs,bnncc,
 
  160      *                 bnndv,bnnvd,bnnds,bnnsd,
 
  162      *                 bptvv,bptss,bptsv,bptvs,bptcc,bptdv,
 
  165      *                 beevv,beess,beesv,beevs,beecc,beedv,
 
  169      *                ,bnnzd,bnndz,bptzd,bptdz,beezd,beedz
 
  170        COMMON /ncoucs/ bcouvv,bcouss,bcousv,bcouvs,
 
  171      *                 bcouzz,bcouhh,bcouds,bcousd,
 
  172      *                 bcoudz,bcouzd,bcoudi,
 
  173      *                 bcoudv,bcouvd,bcoucc
 
  174        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
  175      *                 anndv,annvd,annds,annsd,
 
  177      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
  179      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
  182      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
  183        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
  184      *                 acouzz,acouhh,acouds,acousd,
 
  185      *                 acoudz,acouzd,acoudi,
 
  186      *                 acoudv,acouvd,acoucc
 
  187        common/popcck/pdbck,pdbse,pdbseu,
 
  188      *  ijpock,irejck,ick4,ihad4,ick6,ihad6
 
  189      *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
 
  190      *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
 
  191      *isea43,isea63,irejao
 
  192        COMMON /nstari/nstart
 
  193        COMMON /ncshxx/ncouxh,ncouxt
 
  195       COMMON /nucros/dsigsu,dsigmc,ndsig
 
  196       COMMON /nncms/  gamcm,bgcm,umo,pcm,eproj,pproj
 
  197       COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
 
  198      *               ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
 
  199      *               ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss 
 
  200      *              ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
 
  201      *               nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
 
  202      *               nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss 
 
  203       COMMON /hdjase/nhse1,nhse2,nhse3,nhase1,nhase2,nhase3
 
  204       COMMON /casadi/casaxx,icasad
 
  205       COMMON /vxsvd/vxsp(50),vxst(50),vxsap(50),vxsat(50),
 
  206      *              vxvp(50),vxvt(50),vxdp(50),vxdt(50),
 
  207      *      nxsp,nxst,nxsap,nxsat,nxvp,nxvt,nxdp,nxdt
 
  208       dimension vxsss(50,6),vxvvv(50,6),xxxx(50,6)
 
  209       dimension xb(200),bimpp(200)
 
  214       CHARACTER*8 projty,targty
 
  215       COMMON /user1/titled,projty,targty
 
  216       COMMON /user2/cmener,sdfrac,ptlar,istruf,isingx,idubld
 
  217       COMMON /strufu/istrum,istrut
 
  218       COMMON /ptsamp/ isampt
 
  219       common/collis/
s,ijprox,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
 
  220       COMMON /dropjj/dropjt,dropva
 
  221       COMMON /pomtyp/ipim,icon,isig,lmax,mmax,
nmax,difel,difnu
 
  223       COMMON /zentra/ icentr
 
  227       COMMON /taufo/  taufor,ktauge,itauve,incmod
 
  228       common/popcor/pdb,ajsdef
 
  229       COMMON /diquax/amedd,idiqua,idiquu
 
  230       COMMON /colle/nevhad,nvers,ihadrz,nfile
 
  231       COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
 
  232      +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
 
  233      +prebin,taebin,fermod,etacou
 
  234        COMMON /cronin/cronco,mkcron
 
  235       COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
 
  237       COMMON /secint/isecin
 
  266       OPEN(47,
file=
'GLAUBTAR.DAT',
 
  268       OPEN(37,
file=
'GLAUBCROSSPB.DAT',
 
  306       CALL rluxgo(lux_level,iseed,0,0)
 
  386       CALL 
parpt(1,pt1,pt2,ipt,nevt)
 
  388       CALL 
dminit(ncases,epn,ppn,ncount,iglaub)
 
  406       IF( iglaub.EQ.1) 
THEN 
  414       WRITE(6,*)
' Printout of important Parameters before DPMJET run.' 
  415      *,
' Please note for DPMJET input all numbers are floating point!'  
  416       WRITE(6,*)
'PROJPAR  ',ip,ipz
 
  417       WRITE(6,*)
'TARPAR   ',it,itz
 
  418       WRITE(6,*)
'MOMENTUM ',ppn
 
  419       WRITE(6,*)
'ENERGY   ',epn
 
  420       WRITE(6,*)
'CMENERGY ',umo
 
  421       WRITE(6,*)
'NOFINALE ',ifinal
 
  422       WRITE(6,*)
'EVAPORAT ',ievap
 
  423       WRITE(6,*)
'OUTLEVEL ',ipri,ipev,ippa,ipco,
init,iphkk
 
  424       auauau=rd2out(iseed1,iseed2)
 
  425       WRITE(6,*)
'RANDOMIZ ',iseed1,iseed2, 
' Initial RNDM (RM48) seeds' 
  426       WRITE(6,*)
'STRUCFUN ',istruf+100*istrut
 
  427       WRITE(6,*)
'SAMPT    ',isampt
 
  428       WRITE(6,*)
'SELHARD  ',0,iophrd, 0,dropjt,ptthr,ptthr2 
 
  429       WRITE(6,*)
'SIGMAPOM ',0,isig,ipim+10*icon,imax,mmax,
nmax 
  430       WRITE(6,*)
'PSHOWER  ',ipshow
 
  431       WRITE(6,*)
'CENTRAL  ',icentr
 
  432       WRITE(6,*)
'CMHISTO  ',cmhis
 
  433       WRITE(6,*)
'SEASU3   ',seasq
 
  434       WRITE(6,*)
'RECOMBIN ',irecom
 
  435       WRITE(6,*)
'SINGDIFF ',isingd
 
  436       WRITE(6,*)
'TAUFOR   ',taufor,ktauge,itauve
 
  437       WRITE(6,*)
'POPCORN  ',pdb
 
  438       WRITE(6,*)
'POPCORCK ',ijpock,pdbck
 
  439       WRITE(6,*)
'POPCORSE ',pdbse,pdbseu
 
  440       WRITE(6,*)
'CASADIQU ',icasad,casaxx
 
  441       WRITE(6,*)
'DIQUARKS ',idiqua,idiquu,amedd
 
  442       WRITE(6,*)
'HADRONIZ ',ihadrz
 
  443       WRITE(6,*)
'INTPT    ',intpt
 
  444       WRITE(6,*)
'PAULI    ',lpauli
 
  445       WRITE(6,*)
'FERMI    ',fermp,fermod
 
  446       WRITE(6,*)
'CRONINPT ',mkcron,cronco
 
  447       WRITE(6,*)
'SEADISTR ',xseacu+0.95d0,unon,unom,unosea
 
  448       WRITE(6,*)
'SEAQUARK ',seaqx,seaqxn
 
  449       WRITE(6,*)
'SECINTER ',isecin
 
  450       WRITE(6,*)
'XCUTS    ',cvq,cdq,csea,ssmima
 
  451       WRITE(6,*)
' Printout of important Parameters before DPMJET run.' 
  452      *,
' Please note for DPMJET input all numbers are floating point!'  
  458         ndone=(iiii-1)*ncaset
 
  460         WRITE(6,
'(A,4I5)')
' KKINC call ',iit,iitz,iip,iipz
 
  462  1111   
FORMAT(
' NDONE= ',i10)
 
  465           numev = i+(iiii-1)*ncaset
 
  466           IF ((i.EQ.486).OR.(i.EQ.803).OR.(i.EQ.1368).OR.
 
  467      &    (i.EQ.1465).OR.(i.EQ.1693).OR.(i.EQ.1808)) 
THEN 
  535           CALL 
kkinc(epn,iit,iitz,iip,iipz,iiproj,kkmat,
 
  536      *    iitarg,nhkkh1,irej)
 
  538           IF(irej.EQ.1)go to 765
 
  615           WRITE(6,
'(A,4I5)')
' KKINC call ',iit,iitz,iip,iipz
 
  617       WRITE(6,
'(2I4,I6,4I4,5F10.2,2I3,I2,I4)')kk,isthkk(kk),idhkk(kk),
 
  618      *jmohkk(1,kk),jmohkk(2,kk),jdahkk(1,kk),jdahkk(2,kk), 
 
  619      *(phkk(ll,kk),ll=1,5),idres(kk),idxres(kk),nobam(kk),idbam(kk)
 
  625           WRITE(6,
'(A,4I5)')
' KKINC call ',iit,iitz,iip,iipz
 
  627       WRITE(6,
'(2I4,I6,4I4,5F10.2,2I3,I2,I4)')kk,isthkk(kk),idhkk(kk),
 
  628      *jmohkk(1,kk),jmohkk(2,kk),jdahkk(1,kk),jdahkk(2,kk), 
 
  629      *(phkk(ll,kk),ll=1,5),idres(kk),idxres(kk),nobam(kk),idbam(kk)
 
  637       WRITE(6,*)
' Printout of important Parameters after DPMJET run.' 
  638      *,
' Please note for DPMJET input all numbers are floating point!'  
  639       WRITE(6,*)
'PROJPAR  ',ip,ipz
 
  640       WRITE(6,*)
'TARPAR   ',it,itz
 
  641       WRITE(6,*)
'MOMENTUM ',ppn
 
  642       WRITE(6,*)
'ENERGY   ',epn
 
  643       WRITE(6,*)
'CMENERGY ',umo
 
  644       WRITE(6,*)
'NOFINALE ',ifinal
 
  645       WRITE(6,*)
'EVAPORAT ',ievap
 
  646       WRITE(6,*)
'OUTLEVEL ',ipri,ipev,ippa,ipco,
init,iphkk
 
  647       auauau=rd2out(iseed1,iseed2)
 
  648       WRITE(6,*)
'RANDOMIZ ',iseed1,iseed2, 
' Final RNDM (RM48) seeds' 
  649       WRITE(6,*)
'STRUCFUN ',istruf+100*istrut
 
  650       WRITE(6,*)
'SAMPT    ',isampt
 
  651       WRITE(6,*)
'SELHARD  ',0,iophrd, 0,dropjt,ptthr,ptthr2 
 
  652       WRITE(6,*)
'SIGMAPOM ',0,isig,ipim+10*icon,imax,mmax,
nmax 
  653       WRITE(6,*)
'PSHOWER  ',ipshow
 
  654       WRITE(6,*)
'CENTRAL  ',icentr
 
  655       WRITE(6,*)
'CMHISTO  ',cmhis
 
  656       WRITE(6,*)
'SEASU3   ',seasq
 
  657       WRITE(6,*)
'RECOMBIN ',irecom
 
  658       WRITE(6,*)
'SINGDIFF ',isingd
 
  659       WRITE(6,*)
'TAUFOR   ',taufor,ktauge,itauve
 
  660       WRITE(6,*)
'POPCORN  ',pdb
 
  661       WRITE(6,*)
'POPCORCK ',ijpock,pdbck
 
  662       WRITE(6,*)
'POPCORSE ',pdbse,pdbseu
 
  663       WRITE(6,*)
'CASADIQU ',icasad,casaxx
 
  664       WRITE(6,*)
'DIQUARKS ',idiqua,idiquu,amedd
 
  665       WRITE(6,*)
'HADRONIZ ',ihadrz
 
  666       WRITE(6,*)
'INTPT    ',intpt
 
  667       WRITE(6,*)
'PAULI    ',lpauli
 
  668       WRITE(6,*)
'FERMI    ',fermp,fermod
 
  669       WRITE(6,*)
'CRONINPT ',mkcron,cronco
 
  670       WRITE(6,*)
'SEADISTR ',xseacu+0.95d0,unon,unom,unosea
 
  671       WRITE(6,*)
'SEAQUARK ',seaqx,seaqxn
 
  672       WRITE(6,*)
'SECINTER ',isecin
 
  673       WRITE(6,*)
'XCUTS    ',cvq,cdq,csea,ssmima
 
  674       WRITE(6,*)
' Printout of important Parameters after DPMJET run.' 
  675      *,
' Please note for DPMJET input all numbers are floating point!'  
  681       auauau=rd2out(iseed1,iseed2)
 
  682       WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
 
  683       WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
 
  684       WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
 
  685       WRITE (6,*)
' Final RNDM seeds (RM48) ',iseed1,iseed2
 
  697         WRITE(6,1100) irvv11,irvv12,irvv13,irvv14, irsv11,irsv12,irsv13,
 
  698      +  irsv14, irvs11,irvs12,irvs13,irvs14, irss11,irss12,irss13,irss14
 
  699  1100 
FORMAT (
' REJECTION COUNTERS FROM KKEVT',/, 5
x,
' V-V CHAINS',4i6/
 
  700      +5
x,
' S-V CHAINS',4i6/ 5
x,
' V-S CHAINS',4i6/ 5
x,
' S-S CHAINS',4i6)
 
  701     WRITE(6,
'(A,4I10)')
' POPCCK/SE/S3/S0 rejections ',
 
  702      *  irejck,irejse,irejs3,irejs0
 
  703     WRITE(6,
'(A,4I10)')
' POPCCK/ASE/AS3/AS0 rejections ',
 
  704      *  irejsa,ireja3,ireja0
 
  705     WRITE(6,
'(2A,8I6)')
' POPCCK ICK4,ICK6,IHAD4,IHAD6,ISE4,ISE6 ',
 
  706      *  
'ISE43,ISE63 ', ick4,ick6,ihad4,ihad6,ise4,ise6,ise43,ise63
 
  707     WRITE(6,
'(2A,8I6)')
' POPCSAQ IHADA4,IHADA6,ISEA4,ISEA6 ',
 
  708      *  
'ISEA43,ISEA63 ', ihada4,ihada6,isea4,isea6,isea43,isea63
 
  709       WRITE(6,*)   
' NDVUU,NDVUS,NDVSS,NVDUU,NVDUS,NVDSS',
 
  710      *             
' NDSUU,NDSUS,NDSSS,NSDUU,NSDUS,NSDSS',
 
  711      *             
' NDZUU,NDZUS,NDZSS,NZDUU,NZDUS,NZDSS' ,
 
  712      *               ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
 
  713      *               ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
 
  714      *               ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss  
 
  715       WRITE(6,*)   
' NADVUU,NADVUS,NADVSS,NAVDUU,NAVDUS,NAVDSS',
 
  716      *             
' NADSUU,NADSUS,NADSSS,NASDUU,NASDUS,NASDSS',
 
  717      *             
' NADZUU,NADZUS,NADZSS,NAZDUU,NAZDUS,NAZDSS' ,
 
  718      *               nadvuu,nadvus,nadvss,navduu,navdus,navdss,
 
  719      *               nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
 
  720      *               nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss  
 
  721       WRITE(6,*)
' NHSE1,NHSE2,NHSE3,NHASE1,NHASE2,NHASE3 ',
 
  722      * nhse1,nhse2,nhse3,nhase1,nhase2,nhase3 
 
  729       WRITE(6,
'(A/7I8)')
' Diquark rejection IDIQRE(1-7),N,ss,su,ud',
 
  730      &                                   (idiqre(jj),jj=1,7)
 
  731       WRITE(6,
'(A/7I8)')
' Diquark rejection IDIQRZ(1-7),N,ss,su,ud',
 
  732      &                                   (idiqrz(jj),jj=1,7)
 
  733       WRITE(6,*)
' Diquark rej. IDVRE(1-3),ud,us,ss ',(idvre(jj),jj=1,3)
 
  734       WRITE(6,*)
' Diquark rej. IVDRE(1-3),ud,us,ss ',(ivdre(jj),jj=1,3)
 
  735       WRITE(6,*)
' Diquark rej. IDSRE(1-3),ud,us,ss ',(idsre(jj),jj=1,3)
 
  736       WRITE(6,*)
' Diquark rej. ISDRE(1-3),ud,us,ss ',(isdre(jj),jj=1,3)
 
  737       WRITE(6,*)
' Diquark rej. IDZRE(1-3),ud,us,ss ',(idzre(jj),jj=1,3)
 
  738       WRITE(6,*)
' Diquark rej. IZDRE(1-3),ud,us,ss ',(izdre(jj),jj=1,3)
 
  739       WRITE(6,*)
' NDZSU,NZDSU ',ndzsu,nzdsu
 
  741       IF ((cmhis.EQ.1.d0).AND.(ioudif.EQ.1))
 
  745       WRITE(6,*)
' Output of x-distribution survey',
 
  746      * 
'VXSP(II),VXST(II),VXSAP(II),VXSAT(II),',
 
  747      *
'VXVP(II),VXVT(II),VXDP(II),VXDT(II)' ,
 
  748      *nxsp,nxst,nxsap,nxsat,nxvp,nxvt,nxdp,nxdt
 
  750         IF(nxsp.GE.1)vxsp(ii)=50.d0*vxsp(ii)/nxsp
 
  751         IF(nxst.GE.1)vxst(ii)=50.d0*vxst(ii)/nxst
 
  752         IF(nxsap.GE.1)vxsap(ii)=50.d0*vxsap(ii)/nxsap
 
  753         IF(nxsat.GE.1)vxsat(ii)=50.d0*vxsat(ii)/nxsat
 
  754         IF(nxvp.GE.1)vxvp(ii)=50.d0*vxvp(ii)/nxvp
 
  755         IF(nxvt.GE.1)vxvt(ii)=50.d0*vxvt(ii)/nxvt
 
  756         IF(nxdp.GE.1)vxdp(ii)=50.d0*vxdp(ii)/nxdp
 
  757         IF(nxdt.GE.1)vxdt(ii)=50.d0*vxdt(ii)/nxdt
 
  758     xxxxx=ii*0.02d0-0.01d0
 
  765     fxvvv=(1.-xxxxx)**3/
sqrt(xxxxx)
 
  766     fxddd=2.d0*xxxxx**3.0d0/
sqrt(1.d0-xxxxx)
 
  767     vxsss(ii,1)=log10(vxsp(ii))
 
  768     vxsss(ii,2)=log10(vxst(ii))
 
  769     vxsss(ii,3)=log10(vxsap(ii))
 
  770     vxsss(ii,4)=log10(vxsat(ii))
 
  771     vxvvv(ii,1)=log10(vxvp(ii))
 
  772     vxvvv(ii,2)=log10(vxvt(ii))
 
  773     vxvvv(ii,3)=log10(vxdp(ii))
 
  774     vxvvv(ii,4)=log10(vxdt(ii))
 
  775     vxvvv(ii,5)=log10(fxvvv)
 
  776     vxvvv(ii,6)=log10(fxddd)
 
  777     axsp=axsp+0.02d0*vxsp(ii)*xxxxx
 
  778     axst=axst+0.02d0*vxst(ii)*xxxxx
 
  779     axsap=axsap+0.02d0*vxsap(ii)*xxxxx
 
  780     axsat=axsat+0.02d0*vxsat(ii)*xxxxx
 
  781     axvp=axvp+0.02d0*vxvp(ii)*xxxxx
 
  782     axvt=axvt+0.02d0*vxvt(ii)*xxxxx
 
  783     axdp=axdp+0.02d0*vxdp(ii)*xxxxx
 
  784     axdt=axdt+0.02d0*vxdt(ii)*xxxxx
 
  785     WRITE(6,*)vxsp(ii),vxst(ii),vxsap(ii),vxsat(ii),
 
  786      *        vxvp(ii),vxvt(ii),vxdp(ii),vxdt(ii)
 
  789      *axsp,axst,axsap,axsat,axvp,axvt,axdp,axdt
 
  790       CALL 
plot(xxxx,vxsss,200,4,50,0.d0,0.02d0,-3.d0,0.05d0)
 
  791       CALL 
plot(xxxx,vxvvv,300,6,50,0.d0,0.02d0,-3.d0,0.05d0)
 
  795       CALL 
parpt(3,pt1,pt2,ipt,ncases)
 
  799         fracxs=float(ncouxh)/(float(ncouxh)+float(ncouxt))
 
  801       WRITE(6,*)
' Fraction of x-sect: ',fracxs,ncouxh,ncouxt
 
  806     IF(ndsig.GE.1) dsigmc=dsigsu/ndsig
 
  807     WRITE(6,*)
' Neutrino-nucleon cross section DSIGMC,NDSIG ',
 
  808      &  dsigmc,
' *10**(-38) cm**2 ',ndsig,
' evts' 
  841       IF(ishmal) CALL 
shmak(3,nshmac,np,
nt,ip,it,umo,bimp)
 
  842       IF(ishmal) CALL 
shmak1(3,nshma2,np,
nt,ip,it,umo,bimp)
 
  844       IF (ireso.EQ.1)  CALL 
distrp(3,ncases,ppn)
 
  845       IF (cmhis.EQ.0.d0) CALL 
distr(3,ncases,ppn,idummy)
 
  846       IF (cmhis.EQ.1.d0) CALL 
distrc(3,ncases,ppn,idummy)
 
  847       IF (cmhis.EQ.2.d0) CALL 
distco(3,ncases,ppn,idummy)
 
  848       IF (ireso.EQ.1) CALL 
disres(3,ncases,ppn)
 
  851        IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN 
  852          CALL 
plomb(5,pp,char,xfxfxf,itif,ijproj)
 
  854        IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN 
  855          CALL 
plombc(5,pp,char,xfxfxf,itif,ijproj)
 
  866       SUBROUTINE dminit(NCASES,EPN,PPN,NCOUNT,IGLAUB)
 
  867       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
  979       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 1077       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 1078      +iibar(210),k1(210),k2(210)
 
 1087       COMMON /paname/ btype(30)
 
 1089       COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
 
 1090      +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
 
 1092       COMMON /factmo/ ifacto
 
 1094       COMMON /taufo/  taufor,ktauge,itauve,incmod
 
 1096       COMMON /rptshm/ rproj,rtarg,bimpac
 
 1098       COMMON /trafop/ gamp,bgamp,betp
 
 1100       COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
 
 1101      +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
 
 1102      +irvs14, irvv11,irvv12,irvv13,irvv14
 
 1104       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 1106       COMMON /dnun/   nn,np,
nt 
 1108       COMMON /nshmak/ nnshma,npshma,ntshma,nshmac,nshma2
 
 1110       COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
 
 1111      *              bsite(0:1,200),nstatb,nsiteb
 
 1113       COMMON /hadthr/ ehadth,inthad
 
 1117       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 1118       COMMON /nuccc/   jt,jtz,jp,jpz,jjproj,jbproj,jjtarg,jbtarg
 
 1120       COMMON /zentra/ icentr
 
 1122       COMMON /cmhico/ cmhis
 
 1124       COMMON /resona/ ireso
 
 1126       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
 1128       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
 1129      +ipadis,ishmal,lpauli
 
 1131       COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
 
 1134       COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
 
 1135      +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
 
 1136      +prebin,taebin,fermod,etacou
 
 1142       COMMON /projk/ iprojk
 
 1144       parameter(lunber=14)
 
 1147       COMMON /seaqxx/ seaqx,seaqxn 
 
 1148        COMMON /cronin/cronco,mkcron
 
 1150        COMMON /seadiq/lseadi
 
 1151        COMMON /final/ifinal
 
 1152        COMMON /recom/irecom
 
 1154       COMMON /neutyy/ neutyp,neudec
 
 1155        COMMON /nstari/nstart
 
 1156        common/popcor/pdb,ajsdef
 
 1157        common/popcck/pdbck,pdbse,pdbseu,
 
 1158      *  ijpock,irejck,ick4,ihad4,ick6,ihad6
 
 1159      *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
 
 1160      *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
 
 1161      *isea43,isea63,irejao
 
 1163        COMMON /bufues/ bnnvv,bnnss,bnnsv,bnnvs,bnncc,
 
 1164      *                 bnndv,bnnvd,bnnds,bnnsd,
 
 1166      *                 bptvv,bptss,bptsv,bptvs,bptcc,bptdv,
 
 1167      *                 bptvd,bptds,bptsd,
 
 1169      *                 beevv,beess,beesv,beevs,beecc,beedv,
 
 1170      *                 beevd,beeds,beesd,
 
 1172      *                ,bnndi,bptdi,beedi
 
 1173      *                ,bnnzd,bnndz,bptzd,bptdz,beezd,beedz
 
 1174        COMMON /ncoucs/ bcouvv,bcouss,bcousv,bcouvs,
 
 1175      *                 bcouzz,bcouhh,bcouds,bcousd,
 
 1176      *                 bcoudz,bcouzd,bcoudi,
 
 1177      *                 bcoudv,bcouvd,bcoucc
 
 1178        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 1179      *                 anndv,annvd,annds,annsd,
 
 1181      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 1183      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 1186      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 1187        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 1188      *                 acouzz,acouhh,acouds,acousd,
 
 1189      *                 acoudz,acouzd,acoudi,
 
 1190      *                 acoudv,acouvd,acoucc
 
 1191       COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
 
 1192      *               ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
 
 1193      *               ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss 
 
 1194      *              ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
 
 1195      *               nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
 
 1196      *               nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss 
 
 1197       COMMON /hdjase/nhse1,nhse2,nhse3,nhase1,nhase2,nhase3
 
 1199       COMMON /diffra/ isingd,idiftp,ioudif,iflagd
 
 1200       COMMON /seasu3/seasq
 
 1201       COMMON /ifragm/ifrag
 
 1202       COMMON /fluctu/ifluct
 
 1203       COMMON /diquax/amedd,idiqua,idiquu
 
 1204       COMMON /vxsvd/vxsp(50),vxst(50),vxsap(50),vxsat(50),
 
 1205      *              vxvp(50),vxvt(50),vxdp(50),vxdt(50),
 
 1206      *      nxsp,nxst,nxsap,nxsat,nxvp,nxvt,nxdp,nxdt
 
 1209       COMMON /nncms/  gamcm,bgcm,umo,pcm,eproj,pproj
 
 1212       COMMON /xsecpt/ ptcut,sigs,dsigh
 
 1213       COMMON /kglaub/jglaub
 
 1214       COMMON /xsecnu/ecmuu,ecmoo,ngritt,nevtt
 
 1258       CHARACTER*8 projty,targty
 
 1261       COMMON /user1/titled,projty,targty
 
 1262       COMMON /user2/cmener,sdfrac,ptlar,istruf,isingx,idubld
 
 1286       COMMON /colle/nevhad,nvers,ihadrz,nfile
 
 1298       common/collis/
s,ijprox,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
 
 1307       common/booklt/btypex(30),nbook(30)
 
 1325        COMMON /polmn0/pdifr,phard,psoft,alfah,betah,
 
 1326      *              sigtot,sigqel,sigel,sigine,sighin,sigd,sigdd
 
 1333       COMMON /pomtyp/ipim,icon,isig,lmax,mmax,
nmax,difel,difnu
 
 1340       COMMON /dropjj/dropjt,dropva
 
 1341       COMMON /gluspl/nugluu,nsgluu
 
 1344       COMMON /ptlarg/xsmax
 
 1346       COMMON /ptsamp/ isampt
 
 1347       COMMON /stars/istar2,istar3
 
 1350       COMMON /strufu/istrum,istrut
 
 1351       COMMON /cutofn/ncutox
 
 1354       COMMON /harlun/ qlun,iharlu
 
 1355       COMMON /pomtab/ipomta   
 
 1356       COMMON /sincha/isichaa
 
 1358       COMMON /evappp/ievap
 
 1360       parameter( frdiff = 0.2
d+00 )
 
 1361       parameter( ethsea = 1.0
d+00 )
 
 1363       LOGICAL ldiffr, linctv, levprt, lheavy, ldeexg, lgdhpr, lpreex,
 
 1364      &        lhlfix, lprfix, lparwv, lpower, lsngch, llvmod, lschdf
 
 1365       COMMON / parevt / dpower, fsprd0, fshpfn, rn1gsc, rn2gsc,
 
 1366      &                  ldiffr(39),lpower, linctv, levprt, lheavy,
 
 1367      &                  ldeexg, lgdhpr, lpreex, lhlfix, lprfix, lparwv,
 
 1368      &                  ilvmod, jlvmod, llvmod, lsngch, lschdf
 
 1378       parameter( 
mxpsst =   300 )
 
 1379       parameter( 
mxpsfb = 41000 )
 
 1380       LOGICAL lfrmbk, lncmss
 
 1381       COMMON / frbkcm /  amufbk, eexfbk(
mxpsst), amfrbk(
mxpsst),
 
 1383      &          exmxfb, r0frbk, r0cfbk, c1cfbk, c2cfbk,
 
 1388      &          ifbcha(5,
mxpsfb), iposst, iposfb, ifbstf,
 
 1389      &          ifbfrb, nbufbk, lfrmbk, lncmss
 
 1391       COMMON /inpflg/ iang,ifiss,ib0,igeom,istrag,keydk
 
 1392       COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
 
 1393       COMMON /secint/isecin
 
 1394        COMMON /nuclea/ pfermp(2),pfermn(2),fermdd,
 
 1395      & ebindp(2),ebindn(2),epot(2,210),
 
 1397       COMMON /ferfor/iferfo
 
 1398       COMMON /casadi/casaxx,icasad
 
 1399       COMMON /infore/ifrej         
 
 1409       CHARACTER*8 code,codewd,blank,sdum
 
 1412      1
'TITLE   ',
'PROJPAR ',
'TARPAR  ',
'ENERGY  ',
'HADRONIZ',
 
 1413      2
'RANDOMIZ',
'FERMI   ',
'EVENTAPE',
'START   ',
'PARTEV  ',
 
 1414      3
'INTPT   ',
'TECALBAM',
'RESONANC',
'VALVAL  ',
'COMMENT ',
 
 1415      4
'OUTLEVEL',
'LEPTOEVT',
'SEASEA  ',
'PARTICLE',
'ALLPART ',
 
 1416      5
'TAUFOR  ',
'SEAVAL  ',
'VALSEA  ',
'MOMENTUM',
'PAULI   ',
 
 1417      6
'PROJKASK',
'CENTRAL ',
'SEADISTR',
'CMHISTO ',
'SIGTEST ',
 
 1418      7
'XCUTS   ',
'HADRIN  ',
'FACTOMOM',
'COULOMB ',
'GLAUBERI',
 
 1419      8
'EDENSITY',
'CMENERGY',
'INFOREJE',
'RECOMBIN',
'SINGDIFF',
 
 1420      9
'NOFINALE',
'SEASU3  ',
'CRONINPT',
'POPCORN ',
'STOP    ',
 
 1421      9
'FLUCTUAT',
'DIQUARKS',
'HBOOKHIS',
'GLAUBERA',
'POMTABLE',
 
 1422      9
'SINGLECH',
'HADRINTH',
'EVAPORAT',
'SEAQUARK',
'SECINTER',
 
 1423      9
'POPCORCK',
'CASADIQU',
'POPCORSE',
'NEUTRINO',
'DIFFNUC ',
 
 1424      9
'XSECNUC ',
'        ',
'        ',
'        ',
'        '/
 
 1454       IF (ncount.EQ.1)
THEN 
 1462  1000 
FORMAT( 
'1    **************************************************',
 
 1463      +
'**************************************************', //
 
 1464      +
'     DPMJET VERSION II.5  (Sept. 1999)      ' /
 
 1465      +
'     DUAL PARTON MODEL FOR HADRON NUCLEUS COLLISIONS '/ /
 
 1466      +
'                AND NUCLEUS NUCLEUS COLLISIONS   '/
 
 1467      +
'      INCLUDING A FORMATION TIME INTRANUCLEAR CASCADE'/
 
 1468      4
'      Minijets and DTUJET like multiple soft jets    '/
 
 1469      4
'      Nuclear evaporation and residual target and    '/
 
 1470      4
'      projectile nuclei    '/
 
 1471      +
'     **************************************************',
 
 1472      +
'**************************************************',//)
 
 1574       IF (ihadrz.GE.2)
THEN 
 1616             ptthr=2.1+0.15*(log10(cmener/50.))**3
 
 1618           ELSEIF(istrut.EQ.2)
THEN 
 1619             ptthr=2.5+0.12*(log10(cmener/50.))**3
 
 1637       IF(ijproj.NE.0) nnpp=ijproj
 
 1638       epn=
sqrt(ppn**2+aam(nnpp)**2)
 
 1644       pproj = 
sqrt((epn-amproj)*(epn+amproj))
 
 1645       umo = 
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
 
 1647       gamcm = (eproj+amtar)/umo
 
 1650       pcm=gamcm*pproj - bgcm*eproj
 
 1652             ptthr=2.1+0.15*(log10(cmener/50.))**3
 
 1654           ELSEIF(istrut.EQ.2)
THEN 
 1655             ptthr=2.5+0.12*(log10(cmener/50.))**3
 
 1679         istruf=istruf-istrut*100
 
 1776       IF(ieof.EQ.1)                                             go to 40
 
 1778       READ(5,1010)codewd,(
what(i),i=1,6),sdum
 
 1779       WRITE(6,1020)codewd,(
what(i),i=1,6),sdum
 
 1782         IF(codewd.EQ.
code(isw))                                 go to 30
 
 1790      +  50      ,  60      ,  90      ,  120     ,  130     ,
 
 1794      +  140     ,  150     ,  160     ,  170     ,  210     ,
 
 1798      +  220     ,  230     ,  240     ,  250     ,  260     ,
 
 1802      +  280     ,  290     ,  300     ,  310     ,  320     ,
 
 1806      +  330     ,  340     ,  350     ,  360     ,  370     ,
 
 1810      +  380     ,  390     ,  400     ,  410     ,  420     ,
 
 1814      +  430     ,  440     , 450      , 460      , 470      ,
 
 1818      +  480     ,  490     , 500      , 510      , 520      ,
 
 1822      +  530     ,  535     ,  538,     539,      540  ,    
 
 1826      +  541     ,  542     ,  543,     544, 545,
 
 1831      +   551      , 552      , 553 ,  554      ,555 ,
 
 1835      +  556     , 557      ,558 ,    559   ,560,        
 
 1839      +  620     , 630      ,640 ,    650   ,660,610),isw
 
 1861  1010 
FORMAT(a8,2
x,6e10.0,a8 )
 
 1862  1020 
FORMAT(
' *****NEXT CONTROL CARD ***** ',a10,6(1
x,g11.4), 2
x,a10)
 
 1864  1030 
FORMAT(/,
' UNKNOWN CODEWORD - CONTROL CARD IGNORED')
 
 1865  1040 
FORMAT(/,
' UNEXPECTED END OF INPUT - STOP ASSUMED.')
 
 1866  1050 
FORMAT(/,
' UNEXPECTED END OF INPUT - START ASSUMED.')
 
 1906       IF (codewd.GT.
'-zzzzzzz')
 
 1907      1    
WRITE(6,91) codewd,(
what(i),i=1,6),sdum
 
 1908   91      
FORMAT(
' ---- control input card : ----' 
 1909      1        /1
x,a8,2
x,6(f10.3),a8)
 
 1926       IF(codewd.EQ.
'STRUCFUN') 
THEN 
 1932         istruf=istruf-istrut*100
 
 1934         WRITE(6,*)
' ISTRUF,ISTRUT ',istruf,istrut
 
 1948       ELSEIF(codewd.EQ.
'PSHOWER ') 
THEN 
 1976  1070 
FORMAT(//,5
x,a80,//)
 
 1991       IF(sdum.EQ.blank) 
THEN 
 1997         IF(ip.EQ.1) ijproj=1
 
 1998         IF(ip.EQ.1) ijprox=1
 
 2004         IF(ip.EQ.1) jjproj=1
 
 2005         IF(ip.EQ.1) jjprox=1
 
 2008           IF(sdum.EQ.btype(ii)) 
THEN 
 2011             ibproj=iibar(ijproj)
 
 2016             jbproj=iibar(ijproj)
 
 2022         WRITE(6,
'(A)') 
' WRONG STRUCTURE OF PROJPAR CARD' 
 2039       IF(sdum.EQ.blank) 
THEN 
 2046           IF(sdum.EQ.btype(ii)) 
THEN 
 2056         WRITE(6,
'(A)') 
' WRONG STRUCTURE OF TARPAR CARD' 
 2072       IF(ijproj.NE.0) nnpp=ijproj
 
 2073       ppn=
sqrt((epn-aam(nnpp))*(epn+aam(nnpp)))
 
 2080       pproj = 
sqrt((epn-amproj)*(epn+amproj))
 
 2081       umo = 
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
 
 2084             ptthr=2.1+0.15*(log10(cmener/50.))**3
 
 2086           ELSEIF(istrut.EQ.2)
THEN 
 2087             ptthr=2.5+0.12*(log10(cmener/50.))**3
 
 2090       gamcm = (eproj+amtar)/umo
 
 2093       pcm=gamcm*pproj - bgcm*eproj
 
 2095        print 1033, eproj,pproj,
 
 2096      +amproj,amtar,umo,gamcm,bgcm,pcm
 
 2097  1033 
FORMAT(
' CMS: ' , 
 
 2098      +
'    EPROJ,PPROJ,AMPROJ,AMTAR,UMO,GAMCM,BGCM,PCM'/8e22.13)
 
 2119       IF (ihadrz.GE.2)
THEN 
 2138       auauau=rd2in(iseed1,iseed2)
 
 2156       IF (
what(1).EQ.1.d0)
THEN 
 2163       IF(fermod.LT.0.0d0.OR.scafer.GT.2.0d0) scafer=1.0d0
 
 2181  1080 
FORMAT (
' THIS FILE CONTAINS EVENTS FROM KKEVT ')
 
 2200       IF (ireso.EQ.1)  CALL 
distrp(1,ijproj,ppn)
 
 2201       IF (cmhis.EQ.0.d0) CALL 
distr(1,ijproj,ppn,idummy)
 
 2202       IF (cmhis.EQ.1.d0) CALL 
distrc(1,ijproj,ppn,idummy)
 
 2203       IF (cmhis.EQ.2.d0) CALL 
distco(1,ijproj,ppn,idummy)
 
 2204       IF (ireso.EQ.1) CALL 
disres(1,nhkkh1,ppn)
 
 2205       IF (ipadis) CALL 
distpa(1)
 
 2206       IF (ioudif.EQ.1) CALL 
diadif(1,0)
 
 2207       CALL 
shmak(1,nn,np,
nt,ip,it,umo,bimp)
 
 2208       CALL 
shmak1(1,nn,np,
nt,ip,it,umo,bimp)
 
 2210       WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
 
 2213      *     ,form=
'UNFORMATTED')
 
 2217       WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
 
 2228        IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN 
 2229          CALL 
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
 
 2231        IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN 
 2232          CALL 
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
 
 2241       IF(lpauli .AND. (.NOT.fermp)) 
THEN 
 2243      +  
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
 
 2244      +  
' IF FERMI ACTIVE', 
' LPAULI CHANGED TO .FALSE.' 
 2251                 IF(nevnts.LE.0) nevnts=1000
 
 2255       IF(ncases.LE.0) ncases=100
 
 2258       IF(iglaub.NE.1) iglaub=0
 
 2261       IF(iglaub.EQ.1) 
THEN 
 2264         CALL 
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
 
 2267         CALL 
shmakf(ip,ipz,it,itz)
 
 2276                 IF(ipim.EQ.2)CALL 
prblm2(cmener)
 
 2345             ptthr=2.1+0.15*(log10(cmener/50.))**3
 
 2347           ELSEIF(istrut.EQ.2)
THEN 
 2348             ptthr=2.5+0.12*(log10(cmener/50.))**3
 
 2366       IF (
what(1).EQ.1.d0) 
THEN 
 2384       IF (
what(1).EQ.1.d0) 
THEN 
 2427       IF(
what(1).GT.0.5d0) ireso=1
 
 2440       IF (
what(1).EQ.1.d0) 
THEN 
 2462   270 
WRITE(6,1120)
title 
 2515       OPEN(29,
file=
'lepto.evt',
 
 2519       IF (ireso.EQ.1)  CALL 
distrp(1,ijproj,ppn)
 
 2520       IF (cmhis.EQ.0.d0) CALL 
distr(1,ijproj,ppn,idummy)
 
 2521       IF (cmhis.EQ.1.d0) CALL 
distrc(1,ijproj,ppn,idummy)
 
 2522       IF (cmhis.EQ.2.d0) CALL 
distco(1,ijproj,ppn,idummy)
 
 2523       IF (ireso.EQ.1) CALL 
disres(1,nhkkh1,ppn)
 
 2524       IF (ipadis) CALL 
distpa(1)
 
 2525       IF (ioudif.EQ.1) CALL 
diadif(1,0)
 
 2527       WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
 
 2530      *     ,form=
'UNFORMATTED')
 
 2534       WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
 
 2538     WRITE(6,*)
' NEUTRINO: after INCINI call' 
 2542        IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN 
 2543          CALL 
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
 
 2545        IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN 
 2546          CALL 
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
 
 2555       IF(lpauli .AND. (.NOT.fermp)) 
THEN 
 2557      +  
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
 
 2558      +  
' IF FERMI ACTIVE', 
' LPAULI CHANGED TO .FALSE.' 
 2564                 IF(nevnts.LE.0) nevnts=1000
 
 2568       IF(ncases.LE.0) ncases=100
 
 2574       CALL 
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
 
 2637       WRITE(6,*)
' NEUTRINO initialization finished' 
 2654       IF (
what(1).EQ.1.d0) 
THEN 
 2684       IF (
what(1).EQ.1.d0) 
THEN 
 2724       IF (
what(1).EQ.1.d0) 
THEN 
 2741       IF (
what(1).EQ.1.d0) 
THEN 
 2758       IF(ijproj.NE.0) nnpp=ijproj
 
 2759       epn=
sqrt(ppn**2+aam(nnpp)**2)
 
 2765       pproj = 
sqrt((epn-amproj)*(epn+amproj))
 
 2766       umo = 
sqrt(amproj**2 + amtar**2 + 2.*amtar*eproj)
 
 2769             ptthr=2.1+0.15*(log10(cmener/50.))**3
 
 2771           ELSEIF(istrut.EQ.2)
THEN 
 2772             ptthr=2.5+0.12*(log10(cmener/50.))**3
 
 2775       gamcm = (eproj+amtar)/umo
 
 2778       pcm=gamcm*pproj - bgcm*eproj
 
 2780        print 1033, eproj,pproj,
 
 2781      +amproj,amtar,umo,gamcm,bgcm,pcm
 
 2795       IF (
what(1).EQ.1.d0) 
THEN 
 2855       IF(unon.LT.0.1d0) unon=2.0
 
 2857       IF(unom.LT.0.1d0) unom=1.5
 
 2859       IF(unosea.LT.0.1d0) unosea=2.0
 
 2905       IF(cvq.LT.0.5d0) cvq=1.0
 
 2907       IF(cdq.LT.1.0d0) cdq=2.0
 
 2909       IF(csea.LT.0.1d0) csea =0.1
 
 2911       IF(ssmima.LT.0.0d0) ssmima=0.14
 
 2913       IF(
what(5).GT.2.0d0) vvmthr=
what(5)
 
 2931       IF(inthad.LT.0 .OR. inthad.GT.2) inthad=0
 
 2932       IF(inthad.EQ.1) 
WRITE(6,
'(/5X,A/)')
 
 2933      +
' FHAD: INELASTIC INTERACTION FORCED' 
 2934       IF(inthad.EQ.2) 
WRITE(6,
'(/5X,A/)')
 
 2935      +
' FHAD: ELASTIC INTERACTION FORCED' 
 2994       Write(47,473)ip,ipz,it,itz
 
 2995   473 
FORMAT(
' NUCLEUS  ',4i10)
 
 3008       CALL 
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
 
 3011         WRITE(47,
'(4F10.5)') bmax,bstep,rproj,rtarg
 
 3013         WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
 
 3021       CALL 
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
 
 3022         WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
 
 3056             ptthr=2.1+0.15*(log10(cmener/50.))**3
 
 3058           ELSEIF(istrut.EQ.2)
THEN 
 3059             ptthr=2.5+0.12*(log10(cmener/50.))**3
 
 3066       IF(ijproj.NE.0) nnpp=ijproj
 
 3067       epn=(cmener**2 + aam(nnpp)**2 - aam(1)**2)/(2.*aam(1))
 
 3068       ppn=
sqrt((epn-aam(nnpp))*(epn+aam(nnpp)))
 
 3069       ecm = 
sqrt(aam(idp)**2+aam(1)**2+2.0d0*aam(1)*epn)      
 
 3074       pproj = 
sqrt((epn-amproj)*(epn+amproj))
 
 3075       eproj=
sqrt(pproj**2+amproj**2)
 
 3079       gamcm = (eproj+amtar)/umo
 
 3081       gamcm = (eproj+aam(1))/umo
 
 3083       pcm=gamcm*pproj - bgcm*eproj
 
 3084        print 1033, eproj,pproj,
 
 3085      +amproj,amtar,umo,gamcm,bgcm,pcm
 
 3117       IF (
what(1).EQ.1.d0) irecom=1
 
 3118       IF (
what(1).NE.1.d0) irecom=0
 
 3119       IF (
what(1).EQ.1.d0) lseadi=.true.
 
 3153       IF (
what(1).EQ.1.d0) ifinal=1
 
 3154       IF (
what(1).EQ.0.d0) ifinal=0
 
 3218       IF(ifluct.EQ.1)CALL 
fluini 
 3231       IF(
what(3).GT.0.d0)
THEN 
 3269       Write(47,1473)ip,ipz,it,itz
 
 3270  1473 
FORMAT(
' NUCLEUS  ',4i10)
 
 3279       CALL 
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
 
 3282         WRITE(47,
'(4F10.5)') bmax,bstep,rproj,rtarg
 
 3284         WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
 
 3292       CALL 
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
 
 3293         WRITE(47,
'(5E16.8)') (bsite(1,ib),ib=1,200)
 
 3333       WRITE(6,
'(A,F10.2)')
' Threshold for HADRIN events = (GeV)',
 
 3417       IF ( 
nint(
what(1)) .GE. 10000 ) 
THEN 
 3430       ELSE IF ( 
nint(whtsav) .NE. 0 ) 
THEN 
 3436          IF ( abs(
nint(
what(2))) .GE. 10  ) 
THEN 
 3440          ELSE IF ( 
nint(
what(2)) .NE. 0 ) 
THEN 
 3461       IF(
what(1).EQ.0.)
THEN 
 3488       IF ( 
nint(
what(2)) .LT. 0 ) ldeexg = .false.
 
 3548       IF(
what(2).GE.0.1d0)
THEN 
 3589       OPEN(29,
file=
'qel.evt',
 
 3593       IF (ireso.EQ.1)  CALL 
distrp(1,ijproj,ppn)
 
 3594       IF (cmhis.EQ.0.d0) CALL 
distr(1,ijproj,ppn,idummy)
 
 3595       IF (cmhis.EQ.1.d0) CALL 
distrc(1,ijproj,ppn,idummy)
 
 3596       IF (cmhis.EQ.2.d0) CALL 
distco(1,ijproj,ppn,idummy)
 
 3597       IF (ireso.EQ.1) CALL 
disres(1,nhkkh1,ppn)
 
 3598       IF (ipadis) CALL 
distpa(1)
 
 3599       IF (ioudif.EQ.1) CALL 
diadif(1,0)
 
 3601       WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
 
 3604      *     ,form=
'UNFORMATTED')
 
 3608       WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
 
 3616     WRITE(6,*)
' NEUTRINO: after INCINI call' 
 3620        IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN 
 3621          CALL 
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
 
 3623        IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN 
 3624          CALL 
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
 
 3633       IF(lpauli .AND. (.NOT.fermp)) 
THEN 
 3635      +  
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
 
 3636      +  
' IF FERMI ACTIVE', 
' LPAULI CHANGED TO .FALSE.' 
 3643                 IF(nevnts.LE.0) nevnts=1000
 
 3647       IF(ncases.LE.0) ncases=100
 
 3652       CALL 
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
 
 3715       WRITE(6,*)
' NEUTRINO initialization finished' 
 3729       OPEN(29,
file=
'diffnuc.evt',
 
 3733       IF (ireso.EQ.1)  CALL 
distrp(1,ijproj,ppn)
 
 3734       IF (cmhis.EQ.0.d0) CALL 
distr(1,ijproj,ppn,idummy)
 
 3735       IF (cmhis.EQ.1.d0) CALL 
distrc(1,ijproj,ppn,idummy)
 
 3736       IF (cmhis.EQ.2.d0) CALL 
distco(1,ijproj,ppn,idummy)
 
 3737       IF (ireso.EQ.1) CALL 
disres(1,nhkkh1,ppn)
 
 3738       IF (ipadis) CALL 
distpa(1)
 
 3739       IF (ioudif.EQ.1) CALL 
diadif(1,0)
 
 3741       WRITE(6,*)
' before NUCLEAR.BIN opened LUNBER= ',lunber
 
 3744      *     ,form=
'UNFORMATTED')
 
 3748       WRITE(6,*)
'NUCLEAR.BIN opened LUNBER= ',lunber
 
 3759        IF(ihbook.EQ.1.AND.cmhis.EQ.0.d0)
THEN 
 3760          CALL 
plomb(1,ppnpn,char,xfxfxf,itif,ijproj)
 
 3762        IF(ihbook.EQ.1.AND.cmhis.EQ.1.d0)
THEN 
 3763          CALL 
plombc(1,ppnpn,char,xfxfxf,itif,ijproj)
 
 3772       IF(lpauli .AND. (.NOT.fermp)) 
THEN 
 3774      +  
' ACTIVATION OF PAULI PRINCIPLE ONLY REASONABLE',
 
 3775      +  
' IF FERMI ACTIVE', 
' LPAULI CHANGED TO .FALSE.' 
 3782                 IF(nevnts.LE.0) nevnts=1000
 
 3786       IF(ncases.LE.0) ncases=100
 
 3791       CALL 
shmaki(ip,ipz,it,itz,rproj,rtarg,ppn)
 
 3873       WRITE(6,*)
'call xsglau' 
 3874       CALL 
xsglau(ip,it,ijproj,1)       
 
 3943       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 3953       COMMON /paname/ btype(30)
 
 3955       COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
 
 3956      +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
 
 3957      +prebin,taebin,fermod,etacou
 
 3959       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 3961       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
 3963       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
 3964      +ipadis,ishmal,lpauli
 
 3966       COMMON /hadthr/ ehadth,inthad
 
 3968       COMMON /nshmak/ nnshma,npshma,ntshma,nshmac,nshma2
 
 3970       COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
 
 3971      +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
 
 3972      +irvs14, irvv11,irvv12,irvv13,irvv14
 
 3974       COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
 
 3975      *              bsite(0:1,200),nstatb,nsiteb
 
 3978       COMMON /damp/   ca,ci,ga
 
 3982       COMMON /diffra/ isingd,idiftp,ioudif,iflagd
 
 3983       COMMON /nomije/ ptmije(10),nnmije(10)
 
 3984       DATA ptmije /5.d0,7.d0,9.d0,11.d0,13.d0,15.d0,17.d0
 
 3985      +,19.d0,21.d0,23.d0 /
 
 3988       DATA irco1,irco2,irco3,irco4,irco5 /5*0/
 
 3989       DATA irss11,irss12,irss13,irss14,irsv11,irsv12,irsv13,irsv14 /8*0/
 
 3990       DATA irvs11,irvs12,irvs13,irvs14,irvv11,irvv12,irvv13,irvv14 /8*0/
 
 3994       DATA prepot /210*0.0/
 
 3995       DATA taepot /210*0.0/
 
 3996       DATA taebin,prebin,fermod /2*0.0d0,0.6d0/
 
 3998       DATA btype /
'PROTON  ' , 
'APROTON ' , 
'ELECTRON' , 
'POSITRON' ,
 
 3999      +
'NEUTRIE ' , 
'ANEUTRIE' , 
'PHOTON  ' , 
'NEUTRON ' , 
'ANEUTRON' ,
 
 4000      +
'MUON+   ' , 
'MUON-   ' , 
'KAONLONG' , 
'PION+   ' , 
'PION-   ' ,
 
 4001      +
'KAON+   ' , 
'KAON-   ' , 
'LAMBDA  ' , 
'ALAMBDA ' , 
'KAONSHRT' ,
 
 4002      +
'SIGMA-  ' , 
'SIGMA+  ' , 
'SIGMAZER' , 
'PIZERO  ' , 
'KAONZERO' ,
 
 4003      +
'AKAONZER' , 
'RESERVED' , 
'BLANK   ' , 
'BLANK   ' , 
'BLANK   ' ,
 
 4006       DATA ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr /0, 0, 0, -1, 0,
 
 4009       DATA intpt, fermp, ihadss,ihadsv,ihadvs,ihadvv, ihada /.true.,
 
 4010      +.true., 4*.false., .true./
 
 4011       DATA ipadis, ishmal, lpauli /.false., .false., .true./
 
 4016       DATA nstatb, nsiteb /2000, 200/
 
 4020       DATA isingd,idiftp,ioudif,iflagd /0,0,0,0/
 
 4041       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 4044       CHARACTER*8 projty,targty
 
 4047       COMMON /user1/
title,projty,targty
 
 4048       COMMON /user2/cmener,sdfrac,ptlar,istruf,isingd,idubld
 
 4052       COMMON /colle/nevhad,nvers,ihadrz,nfile
 
 4056       common/collis/
s,ijproj,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
 
 4063       common/booklt/btype(30),nbook(30)
 
 4068        COMMON /polmn0/pdifr,phard,psoft,alfah,betah,
 
 4069      *              sigtot,sigqel,sigel,sigine,sighin,sigd,sigdd
 
 4073       COMMON /pomtyp/ipim,icon,isig,lmax,mmax,
nmax,difel,difnu
 
 4078       COMMON /dropjj/dropjt,dropva
 
 4079       COMMON /gluspl/nugluu,nsgluu
 
 4080       COMMON /ptlarg/xsmax
 
 4081       COMMON /ptsamp/ isampt
 
 4082       COMMON /stars/istar2,istar3
 
 4083       COMMON /strufu/istrum,istrut
 
 4084       COMMON /popcor/pdb,ajsdef 
 
 4089       CHARACTER*8 codewd,sdum
 
 4099  9    
FORMAT( 
' special code word was used ')
 
 4177       IF(codewd.EQ.
'SIGMAPOM') 
THEN 
 4196                    IF (itest.EQ.1)CALL 
pomdi 
 4210       ELSEIF(codewd.EQ.
'GLUSPLIT') 
THEN 
 4233       ELSEIF(codewd.EQ.
'PARTEV  ') 
THEN 
 4236                 IF (
what(2).EQ.0.d0)npev=30
 
 4238                 IF (
what(3).EQ.0.d0)nvers=1
 
 4243                   IF(ipim.EQ.2)CALL 
prblm2(cmener)
 
 4265       ELSEIF(codewd.EQ.
'SELHARD ') 
THEN 
 4270         IF(
what(5).NE.0.d0)
THEN 
 4272           IF(cmener.LT.2000.0d0.AND.isig.EQ.3)ptthr=
what(5)
 
 4273           IF (cmener.GE.2000.0d0.AND.isig.EQ.3)
 
 4274      *                   ptthr=0.25*
log(cmener/2000.)+2.
 
 4275           IF(ptthr2.LT.ptthr)ptthr2=ptthr
 
 4277             ptthr=2.1+0.15*(log10(cmener/50.))**3
 
 4279           ELSEIF(istrut.EQ.2)
THEN 
 4280             ptthr=2.5+0.12*(log10(cmener/50.))**3
 
 4284  1244     
FORMAT (
' THRESHOLD PT FOR HARD SCATTERING PTTHR=',f12.2)
 
 4299       ELSEIF(codewd.EQ.
'XSLAPT  ') 
THEN 
 4314       ELSEIF(codewd.EQ.
'SAMPT   ') 
THEN 
 4317         IF( isampt.LT.0 .OR. isampt.GT.4 ) isampt=0
 
 4355       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 4358       common/booklt/btype(30),nbook(30)
 
 4360       DATA btype /
'PROTON  ' , 
'APROTON ' , 
'ELECTRON' ,
 
 4361      1            
'POSITRON' , 
'NEUTRIE ' , 
'ANEUTRIE' ,
 
 4362      2            
'PHOTON  ' , 
'NEUTRON ' , 
'ANEUTRON' ,
 
 4363      3            
'MUON+   ' , 
'MUON-   ' , 
'KAONLONG' ,
 
 4364      4            
'PION+   ' , 
'PION-   ' , 
'KAON+   ' ,
 
 4365      5            
'KAON-   ' , 
'LAMBDA  ' , 
'ALAMBDA ' ,
 
 4366      6            
'KAONSHRT' , 
'SIGMA-  ' , 
'SIGMA+  ' ,
 
 4367      7            
'SIGMAZER' , 
'PIZERO  ' , 
'KAONZERO' ,
 
 4368      9            
'AKAONZER' , 
'        ' , 
'        ' ,
 
 4372       DATA nbook / 2212      , -2212      ,  11        ,
 
 4374      2             22        ,  2112      , -2112      ,
 
 4376      4             211       , -211       ,  321       ,
 
 4377      5            -321       ,  3122      , -3122      ,
 
 4378      6             310       ,  3114      ,  3224      ,
 
 4379      7             3214      ,  111       ,  311       ,
 
 4396       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 4398       parameter( 
zero=0.d0, 
one=1.d0)
 
 4399       parameter( alfa=0.56268
d-01, 
beta=0.17173
d+03 )
 
 4400       parameter( acc = 0.0001d0 )
 
 4401       COMMON /xsecpt/ ptcut,sigs,dsigh
 
 4402       COMMON /sigma / sigsof,bs,zsof,sighar,
fill(7)
 
 4403       COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
 
 4405       common/collis/
s,ijproj,ijtar,ptthr,ptthr2,iophrd,ijprlu,ijtalu
 
 4407       CHARACTER*8 projty,targty
 
 4410       COMMON /user1/
title,projty,targty
 
 4411       COMMON /user2/cmener,sdfrac,ptlar,istruf,isingd,idubld
 
 4413       common/ptsamp/ isampt
 
 4414       dimension pptt(50),dpptt(50)
 
 4417       IF ( mode.EQ.0 ) 
THEN 
 4419           pptt(ii)=ii*ptcut/50.
 
 4423         IF(ecm.LT.1000.)
THEN 
 4424           aacucu=0.85*(ecm-400.)/600.
 
 4425           sigs=(1.-aacucu)*sigsof
 
 4431     WRITE(6,
'(A,4E12.4)')
' SAMPPT:ECM,PTCUT,SIGS,DSIGH',
 
 4432      *  ecm,ptcut,sigs,dsigh
 
 4433         IF( isampt.EQ.0 ) 
THEN 
 4434           c  = dsigh/(2.*sigs*ptcut)
 
 4436         IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,
b 
 4437      * ,
c,sigsof,sighar,rmin
 
 4438         ELSEIF( isampt.EQ.1 ) 
THEN 
 4442         IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,
b 
 4443      * ,
c,sigsof,sighar,rmin
 
 4444         ELSEIF( isampt.EQ.2 ) 
THEN 
 4446         IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,
b 
 4447      * ,
c,sigsof,sighar,rmin
 
 4448         ELSEIF( isampt.EQ.3 ) 
THEN 
 4450         IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,
b 
 4451      * ,
c,sigsof,sighar,rmin
 
 4452         ELSEIF( isampt.EQ.4)
THEN 
 4453           aaaa=ptcut**2*(sigsof+sighar)
 
 4454           IF (aaaa.LE.0.00001d0) 
THEN 
 4455             aaaa=abs(aaaa)+0.0002
 
 4456             WRITE(6,5559)ptcut,sigsof,sighar
 
 4457  5559       
FORMAT(
' SAMPPT:PTCUT,SIDSOF.SIGHATD:',3e12.3)
 
 4461         IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,
b 
 4462      * ,
c,sigsof,sighar,rmin
 
 4465         rmin = 
exp(
b*ptcut**2)
 
 4469         IF( ioutpa.GE.-1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,
b 
 4470      * ,
c,sigsof,sighar,rmin
 
 4471 9010  
FORMAT(
' SAMPPT MODE,ISAMPT,PTCUT,SIGS,DSIGH,B,C,SIGSOF',
 
 4475       ELSEIF ( mode.EQ.1 ) 
THEN 
 4476         IF( ioutpa.GE.1 )
WRITE(6,9010)mode,isampt,ptcut,sigs,dsigh,
b 
 4477      * ,
c,sigsof,sighar,rmin
 
 4478         ptt   =
log(1.0-
rndm(v)*(1.0-rmin))/(
b+0.00001d0)
 
 4480         iipt=
pt*50./ptcut+1.
 
 4482         dpptt(iipt)=dpptt(iipt)+1./(
pt+0.000001d0)
 
 4485       ELSEIF(mode.EQ.2)
THEN 
 4487           dpptt(ii)=log10(1.
e-8+dpptt(ii))
 
 4489         IF(iouxev.GE.-1)
THEN 
 4491   203   
FORMAT(
' PT DISTRIBUTION OF SOFT PARTONS AS SAMPLED IN BSOFPT')
 
 4500      *
 FUNCTION bsofpt(ACC,CC,PPTCUT)
 
 4501       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 4504       COMMON /bsoff1/
c,ptcut
 
 4505       COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
 
 4507       dimension 
x(50),
y(50)
 
 4515       IF(
c.LT.1.
d-10) 
THEN 
 4534       IF (.NOT.succes)
THEN 
 4535         IF (kkkk.EQ.0)go to 400
 
 4538       IF(iouxev.GE.1)
WRITE(6,111)b1,b2
 
 4549       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 4551       COMMON /bsoff1/
c,ptcut
 
 4554       df=
c*ptcut**2*aaa-aaa
 
 4564       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 4566       COMMON /bsoff1/
c,ptcut
 
 4569       IF(qqq.GT.-60.) 
THEN 
 4581      *
  FUNCTION rtsafe(FUNCD,X1,X2,XACC)
 
 4582       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 4584       parameter(maxit=200,itepri=0)
 
 4585       COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
 
 4586       CALL funcd(
x1,fl,df)
 
 4588       CALL funcd(
x2,fh,df)
 
 4590       IF(fl*fh.GE.0.) pause 
'ROOT MUST BE BRACKETED' 
 4616      *      .OR. abs(2.*
f).GT.abs(dxold*df) ) 
THEN 
 4629         IF(abs(
dx).LT.xacc) 
RETURN 
 4639       pause 
'RTSAFE EXCEEDING MAXIMUM ITERATIONS' 
 4641 9995  
FORMAT(
'  VR1,VR2:',2e12.5)
 
 4642 9996  
FORMAT(
'  RTSAFE,XH,XL,DXOLD,F,DF IN LOOP 11 J=1,MAXIT')
 
 4643 9997  
FORMAT(3
x,6e10.3)
 
 4644 9998  
FORMAT(
'  RTSAFE: RTSAFE,F,DF =',3e12.5)
 
 4645 9999  
FORMAT(
'  RTSAFE: F,DF =',2e12.5)
 
 4651      *
  FUNCTION var(A,B,C,D)
 
 4652       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 4654       parameter( ambmax = 1.0
d+38, epsi = 1.2
d-38, 
one=1.d0 )
 
 4658       abl = log10( abl + epsi )
 
 4661       ccl = log10( ccl + epsi )
 
 4663       IF( rcheck .LE. 38.d0 ) 
THEN 
 4666         var = ambmax*siab*sicc - 
d 
 4668       IF( 
var .GT. 1.0
d+18 ) 
var = 1.0
e+18
 
 4669       IF( 
var .LT. -1.0
d+18 ) 
var = -1.0
e+18
 
 4674       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 4677       parameter(factor=1.6d0,ntry=50)
 
 4679       IF(
x1.EQ.
x2)pause 
'You have to guess an initial range' 
 4684         IF(
f1*
f2.LT.0.d0)
RETURN 
 4685         IF(abs(
f1).LT.abs(
f2))
THEN 
 4718       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 4722       parameter(amuamu=0.93149432d0)
 
 4725       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 4813       COMMON /delp/ delpx,delpy,delpz,delpe
 
 4815       COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
 
 4817       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 4819       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 4824       IF(help.GT.5.d0)phelp=help-5.
 
 4839         IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13) 
THEN 
 4844           pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
 
 4845           pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
 
 4848         IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14) 
THEN 
 4855         IF(isthkk(i).EQ.1) 
THEN 
 4862         IF(isthkk(i).EQ.13.AND.jdahkk(2,i).EQ.0) 
THEN 
 4867           pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
 
 4868           pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
 
 4870         IF(isthkk(i).EQ.14.AND.jdahkk(2,i).EQ.0) 
THEN 
 4877         IF(isthkk(i).EQ.16) 
THEN 
 4883           eext=eext + phkk(4,i) - phkk(4,imo)
 
 4885         IF(isthkk(i).EQ.15) 
THEN 
 4891           pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
 
 4892           pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
 
 4893           eext=eext + phkk(4,i) - phkk(4,imo)
 
 4895         IF(isthkk(i).EQ.1) 
THEN 
 4898         IF(isthkk(i).EQ.-1) 
THEN 
 4899       eeem1=eeem1+phkk(4,i)
 
 4901         IF(isthkk(i).EQ.1001) 
THEN 
 4902       ee1001=ee1001+phkk(4,i)
 
 4905       eee=eee1+eeem1+ee1001
 
 4916       aip=aip+(ait*amuamu+1.
d-3*
energy(ait,aitz))/epnto
 
 4922         IF(it.EQ.ip)tole=0.02
 
 4925       IF(delle.GE.tole)irej=1
 
 4928     IF(icheck.LE.100)
THEN 
 4929           WRITE(6,
'(A,I5,E10.3,5F10.4)')
 
 4930      *    
' IP,EPN,AEEE,AEEEE,AEEE1,AEEEM1,AEE101:',
 
 4931      *    ip,epn,aeee,aeeee,aeee1,aeeem1,aee101
 
 4932           WRITE(6,
'(A,I5,E10.3,7E12.4)')
 
 4933      *    
' IP,EPN,EEE,EEEE,EEE1,EEEM1,EE1001,DELLE,ELLE:',
 
 4934      *    ip,epn,eee,eeee,eee1,eeem1,ee1001,delle,elle
 
 4950  1000 
FORMAT(
' CHECKF: PX,PY,PZ,PE,EEXT,EEXP',2f7.3,2e12.3,2f7.3
 
 4951      * / 8
x,
' DELPX/Y/Z/E',4f7.3,i10,
' IORIG')
 
 4957             IF(isthkk(i).EQ.11) 
THEN 
 4958               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 4959      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 4960      +        (vhkk(khkk,i),khkk=1,4)
 
 4963             IF(isthkk(i).EQ.12) 
THEN 
 4964               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 4965      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 4966      +        (vhkk(khkk,i),khkk=1,4)
 
 4969             IF(isthkk(i).EQ.1) 
THEN 
 4970               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 4971      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 4972      +        (vhkk(khkk,i),khkk=1,4)
 
 4974  1010 
FORMAT (i6,i4,5i6,9(1pe10.2))
 
 4976             IF(isthkk(i).EQ.16) 
THEN 
 4978               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 4979      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 4980      +        (vhkk(khkk,i),khkk=1,4)
 
 4993       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 4997       parameter(amuamu=0.93149432d0)
 
 5000       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 5088       COMMON /delp/ delpx,delpy,delpz,delpe
 
 5090       COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
 
 5092       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 5094       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 5099       IF(help.GT.5.d0)phelp=help-5.d0
 
 5100       pthelp=12.d0+phelp*5.d0
 
 5120         IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13) 
THEN 
 5125           pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
 
 5126           pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
 
 5129         IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14) 
THEN 
 5136         IF(isthkk(i).EQ.1) 
THEN 
 5143         IF(isthkk(i).EQ.13.AND.jdahkk(2,i).EQ.0) 
THEN 
 5148           pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
 
 5149           pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
 
 5151         IF(isthkk(i).EQ.14.AND.jdahkk(2,i).EQ.0) 
THEN 
 5158         IF(isthkk(i).EQ.16) 
THEN 
 5164           eext=eext + phkk(4,i) - phkk(4,imo)
 
 5166         IF(isthkk(i).EQ.15) 
THEN 
 5172           pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
 
 5173           pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
 
 5174           eext=eext + phkk(4,i) - phkk(4,imo)
 
 5176         IF(isthkk(i).EQ.1) 
THEN 
 5181         IF(isthkk(i).EQ.-1) 
THEN 
 5182       eeem1=eeem1+phkk(4,i)
 
 5186         IF(isthkk(i).EQ.1001) 
THEN 
 5187       ee1001=ee1001+phkk(4,i)
 
 5188       pz1001=pz1001+phkk(3,i)
 
 5189       px1001=px1001+phkk(1,i)
 
 5192       eee=eee1+eeem1+ee1001
 
 5193       pzpz=pz1+pzm1+pz1001
 
 5194       pxpx=px1+pxm1+px1001
 
 5205         IF(isthkk(i).EQ.1001) 
THEN 
 5206       phkk(3,i)=phkk(3,i)+delpz
 
 5207       phkk(4,i)=
sqrt(phkk(1,i)**2+phkk(2,i)**2+phkk(3,i)**2
 
 5209       ee1001=ee1001+phkk(4,i)
 
 5212       eee=eee1+eeem1+ee1001
 
 5220       bip=epn+(ait*amuamu+1.
d-3*
energy(ait,aitz))
 
 5226       IF(delle.GE.tole)irej=1
 
 5229     IF(icheck.LE.20)
THEN 
 5230           WRITE(6,
'(A,I5,E10.3,4F10.4)')
 
 5231      *    
' IP,EPN,PXPX,PX1,PXM1,PX1001:',
 
 5232      *    ip,epn,pxpx,px1,pxm1,px1001
 
 5233           WRITE(6,
'(A,I5,E10.3,6F10.4)')
 
 5234      *    
' IP,PPN,PZPZ,PZ1,PZM1,PZ1001,BIP,BMI:',
 
 5235      *    ip,ppn,pzpz,pz1,pzm1,pz1001,bip,bmi
 
 5236           WRITE(6,
'(A,I5,E10.3,5E12.4)')
 
 5237      *    
' IP,EPN,EEE,EEE1,EEEM1,EE1001,DELLE:',
 
 5238      *    ip,epn,eee,eee1,eeem1,ee1001,delle
 
 5254  1000 
FORMAT(
' CHECKF: PX,PY,PZ,PE,EEXT,EEXP',2f7.3,2e12.3,2f7.3
 
 5255      * / 8
x,
' DELPX/Y/Z/E',4f7.3,i10,
' IORIG')
 
 5261             IF(isthkk(i).EQ.11) 
THEN 
 5262               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5263      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5264      +        (vhkk(khkk,i),khkk=1,4)
 
 5267             IF(isthkk(i).EQ.12) 
THEN 
 5268               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5269      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5270      +        (vhkk(khkk,i),khkk=1,4)
 
 5273             IF(isthkk(i).EQ.1) 
THEN 
 5274               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5275      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5276      +        (vhkk(khkk,i),khkk=1,4)
 
 5278  1010 
FORMAT (i6,i4,5i6,9(1pe10.2))
 
 5280             IF(isthkk(i).EQ.16) 
THEN 
 5282               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5283      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5284      +        (vhkk(khkk,i),khkk=1,4)
 
 5296       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 5300       parameter(amuamu=0.93149432d0)
 
 5303       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 5390       COMMON /zentra/ icentr
 
 5392       COMMON /delp/ delpx,delpy,delpz,delpe
 
 5394       COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
 
 5396       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 5398       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 5403       IF(help.GT.5.d0)phelp=help-5.
 
 5418         IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13) 
THEN 
 5423           pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
 
 5424           pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
 
 5427         IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14) 
THEN 
 5434         IF(isthkk(i).EQ.1) 
THEN 
 5441         IF(isthkk(i).EQ.13.AND.jdahkk(2,i).EQ.0) 
THEN 
 5446           pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
 
 5447           pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
 
 5449         IF(isthkk(i).EQ.14.AND.jdahkk(2,i).EQ.0) 
THEN 
 5456         IF(isthkk(i).EQ.16) 
THEN 
 5462           eext=eext + phkk(4,i) - phkk(4,imo)
 
 5464         IF(isthkk(i).EQ.15) 
THEN 
 5470           pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
 
 5471           pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
 
 5472           eext=eext + phkk(4,i) - phkk(4,imo)
 
 5474         IF(isthkk(i).EQ.1) 
THEN 
 5504       IF(ip.EQ.it.AND.it.GT.1)tole=0.05d0*ip
 
 5506       IF(epn.LE.5.d0)tole=3.d0*tole
 
 5508         IF (abs(
px).GT.pthelp.OR. abs(
py).GT.pthelp.OR. 
 
 5509      *  abs(
pz)/epn.GT.tole.
 
 5510      +  or. abs(pe)/epn.GT.tole) 
THEN 
 5513           IF(icheck.LE.50.AND.irej.EQ.1)
THEN  
 5514       WRITE(6,1000) 
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe,iorig
 
 5516  1000 
FORMAT(
' CHECKO: PX,PY,PZ,PE,EEXT,EEXP',2f7.3,2e12.3,2f7.3
 
 5517      * / 8
x,
' DELPX/Y/Z/E',4f7.3,i10,
' IORIG')
 
 5521       WRITE(6,1000) 
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe,iorig
 
 5525             IF(isthkk(i).EQ.11) 
THEN 
 5526               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5527      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5528      +        (vhkk(khkk,i),khkk=1,4)
 
 5531             IF(isthkk(i).EQ.12) 
THEN 
 5532               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5533      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5534      +        (vhkk(khkk,i),khkk=1,4)
 
 5537             IF(isthkk(i).EQ.1) 
THEN 
 5538               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5539      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5540      +        (vhkk(khkk,i),khkk=1,4)
 
 5542  1010 
FORMAT (i6,i4,5i6,9(1pe10.2))
 
 5544             IF(isthkk(i).EQ.16) 
THEN 
 5546               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5547      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5548      +        (vhkk(khkk,i),khkk=1,4)
 
 5555       WRITE(6,1000) 
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe,iorig
 
 5561       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 5567       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 5655       COMMON /delp/ delpx,delpy,delpz,delpe
 
 5657       COMMON /tanuin/ tasuma,tasubi,tabi,tamasu,tama,taimma
 
 5659       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 5669         IF(isthkk(i).EQ.11.OR.isthkk(i).EQ.13) 
THEN 
 5674           pz=
pz - gam*phkk(3,i) - bgam*phkk(4,i)
 
 5675           pe=pe - gam*phkk(4,i) - bgam*phkk(3,i)
 
 5678         IF(isthkk(i).EQ.12.OR.isthkk(i).EQ.14) 
THEN 
 5685         IF(isthkk(i).EQ.1) 
THEN 
 5692         IF(isthkk(i).EQ.13.AND.jdahkk(1,i).EQ.0) 
THEN 
 5697           pz=
pz + gam*phkk(3,i) + bgam*phkk(4,i)
 
 5698           pe=pe + gam*phkk(4,i) + bgam*phkk(3,i)
 
 5700         IF(isthkk(i).EQ.14.AND.jdahkk(1,i).EQ.0) 
THEN 
 5707         IF(isthkk(i).EQ.16) 
THEN 
 5713           eext=eext + phkk(4,i) - phkk(4,imo)
 
 5715         IF(isthkk(i).EQ.15) 
THEN 
 5721           eext=eext + phkk(4,i) - phkk(4,imo)
 
 5728       WRITE(6,1000) 
px,
py,
pz,pe,eext,eexp,delpx,delpy,delpz,delpe
 
 5729  1000 
FORMAT(
' CHECKE: PX,PY,PZ,PE,EEXT,EEXP',6f7.3/ 8
x,
' DELPX/Y/Z/E',4
 
 5731       WRITE(6,
'(8X,A,6F8.3)') 
' TASUMA,TASUBI,TABI,TAMASU,TAMA,TAIMMA',
 
 5732      +tasuma,tasubi,tabi,tamasu,tama,taimma
 
 5734         IF (abs(
px).GT.0.004.OR. abs(
py).GT.0.004.OR. abs(
pz).GT.0.004.
 
 5735      +  or. abs(pe).GT.0.004) 
THEN 
 5739             IF(isthkk(i).EQ.11) 
THEN 
 5740               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5741      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5742      +        (vhkk(khkk,i),khkk=1,4)
 
 5745             IF(isthkk(i).EQ.12) 
THEN 
 5746               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5747      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5748      +        (vhkk(khkk,i),khkk=1,4)
 
 5751             IF(isthkk(i).EQ.1) 
THEN 
 5752               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5753      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5754      +        (vhkk(khkk,i),khkk=1,4)
 
 5756  1010 
FORMAT (i6,i4,5i6,9(1pe10.2))
 
 5758             IF(isthkk(i).EQ.16) 
THEN 
 5760               WRITE(6,1010) i,isthkk(i),idhkk(i),jmohkk(1,i),jmohkk
 
 5761      +        (2,i), jdahkk(1,i),jdahkk(2,i),(phkk(khkk,i),khkk=1,5),
 
 5762      +        (vhkk(khkk,i),khkk=1,4)
 
 5790       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 5797       DATA a1,a2,a3,a4,a5 /0.01575, 0.0178, 0.000710, 0.0237, 0.034/
 
 5800       IF(ia.LE.1.OR.
iz.EQ.0)
THEN 
 5805       ebind = a1*aa - a2*aa**0. 666667- a3*
iz*
iz*aa**(-0.333333) - a4
 
 5807       IF (
mod(ia,2).EQ.1) 
THEN 
 5809       ELSEIF (
mod(
iz,2).EQ.1) 
THEN 
 5821       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 5828       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 5926       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 5927      +iibar(210),k1(210),k2(210)
 
 5930       COMMON /factmo/ ifacto
 
 5932       COMMON /taufo/  taufor,ktauge,itauve,incmod
 
 5934       COMMON /hadthr/ ehadth,inthad
 
 5938       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 5939       COMMON /nuccc/   jt,jtz,jp,jpz,jjproj,jbproj,jjtarg,jbtarg
 
 5941       COMMON /zentra/ icentr
 
 5943       COMMON /cmhico/ cmhis
 
 5945       COMMON /resona/ ireso
 
 5947       COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
 
 5954       COMMON /projk/ iprojk
 
 5956       COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
 
 5957      +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
 
 5958      +prebin,taebin,fermod,etacou
 
 5960        COMMON /recom/irecom
 
 5970       ppn=
sqrt((epn-aam(ijproj))*(epn+aam(ijproj)))
 
 5971       ibproj=iibar(ijproj)
 
 5974       jbproj=iibar(ijproj)
 
 6026       SUBROUTINE hadhad(EPN,PPN,NHKKH1,IHTAWW,ITTA,IREJFO)
 
 6027       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6037       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 6125       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 6137       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 6138      +iibar(210),k1(210),k2(210)
 
 6141       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 6143       parameter(maxfin=10)
 
 6144       COMMON /dfinls/ itrh(maxfin),cxrh(maxfin),cyrh(maxfin), czrh
 
 6145      +(maxfin),elrh(maxfin),plrh(maxfin),irh
 
 6147       COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
 
 6148      +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
 
 6149      +prebin,taebin,fermod,etacou
 
 6156       IF(ipri.GE.2) 
WRITE(6,1001) ijproj,ijproj,ppn,epn,cccxp,cccyp,
 
 6157      +ccczp,ihtaww,itta,ieline
 
 6158  1001 
FORMAT(
' HADHAD 1:',
 
 6159      +
' IJPROJ,IJPROJ,PPN,EPN,CCCXP,CCCYP,CCCZP,IHTAWW,ITTA,IELINE'/ 2
 
 6160      +i3,2e12.3,3f7.3,3i4)
 
 6165       CALL 
sihnin(ijproj,itta,ppn,sight)
 
 6166       CALL 
sihnel(ijproj,itta,ppn,sighte)
 
 6167       sigtot=sight + sighte
 
 6168       IF (sigtot*
rndm(bb).LE.sighte)ieline=1
 
 6169       IF(ipri.GE.2) 
WRITE(6,1000) ijproj,ijproj,ppn,epn,cccxp,cccyp,
 
 6170      +ccczp,ihtaww,itta,ieline
 
 6171  1000 
FORMAT(
' HADHAD 2 nach si...:',
 
 6172      +
' IJPROJ,IJPROJ,PPN,EPN,CCCXP,CCCYP,CCCZP,IHTAWW,ITTA,IELINE'/ 2
 
 6173      +i3,2e12.3,3f7.3,3i4)
 
 6177       IF(ipri.GE.2) 
WRITE(6,1012) ijproj,ijproj,ppn,epn,cccxp,cccyp,
 
 6178      +ccczp,ihtaww,itta,ieline
 
 6179  1012 
FORMAT(
' HADHAD 12 loop:',
 
 6180      +
' IJPROJ,IJPROJ,PPN,EPN,CCCXP,CCCYP,CCCZP,IHTAWW,ITTA,IELINE'/ 2
 
 6181      +i3,2e12.3,3f7.3,3i4)
 
 6184           WRITE(6,
'(I3,5(1PE12.4),I5/3X,5(1PE12.4))') ii,elrh(ii),plrh
 
 6185      +    (ii),cxrh(ii),cyrh(ii),czrh(ii),itrh(ii), (phkk(jjj,nhkk),jjj
 
 6190       CALL 
fhad(ijproj,ijproj,ppn,epn,cccxp,cccyp,ccczp, ihtaww,itta,
 
 6195      + 
WRITE(6,
'(A)')
'  exit from hadhad with irejfo=1 ' 
 6201       IF (ihadha.LT.3)
THEN 
 6204           IF(itsec.EQ.1.AND.elrh(ii).LE.taefep+aam(itsec))     goto 12
 
 6205           IF(itsec.EQ.8.AND.elrh(ii).LE.taefen+aam(itsec))     goto 12
 
 6206           IF(iibar(itsec).NE.1.AND.elrh(ii)-aam(itsec)
 
 6207      +                                   .LE.taepot(itsec))  goto 12
 
 6212       IF (ipri.GE.2) 
WRITE (6,1010)irh,nhkkh1,ihtaww,itta
 
 6213  1010 
FORMAT (
' HADHAD IRH,NHKKH1,IHTAWW,ITTA = ',4i5)
 
 6217      +  
' HADHAD - PARTICLE TRANSFER FROM /FINLSP/ INTO /HKKEVT/',
 
 6218      +  
' II, ELRH, PLRH, CXRH, CYRH, CZRH / PHKK(1-5)' 
 6227            WRITE (6,
'(A,2I5)') .EQ.
' HADHAD:NHKKNMXHKK ',nhkk,
nmxhkk 
 6231           idhkk(nhkk)=
mpdgha(itsec)
 
 6233           jmohkk(2,nhkk)=ihtaww
 
 6236         phkk(1,nhkk)=plrh(ii)*cxrh(ii)
 
 6237         phkk(2,nhkk)=plrh(ii)*cyrh(ii)
 
 6238         phkk(3,nhkk)=plrh(ii)*czrh(ii)
 
 6239         phkk(4,nhkk)=elrh(ii)
 
 6240         IF(phkk(4,nhkk)-aam(itsec).LE.taepot(itsec).
 
 6241      +                                    and.iibar(itsec).EQ.1)
THEN 
 6246         phkk(5,nhkk)=aam(itrh(ii))
 
 6249           WRITE(6,
'(I3,5(1PE12.4),I5/3X,5(1PE12.4),I5)') 
 
 6251      +    (ii),cxrh(ii),cyrh(ii),czrh(ii),itrh(ii), (phkk(jjj,nhkk),jjj
 
 6254         vhkk(1,nhkk)=vhkk(1,ihtaww)
 
 6255         vhkk(2,nhkk)=vhkk(2,ihtaww)
 
 6256         vhkk(3,nhkk)=vhkk(3,ihtaww)
 
 6257         vhkk(4,nhkk)=vhkk(4,1)
 
 6260       jdahkk(1,1)=nhkkh1+1
 
 6262       jdahkk(1,ihtaww)=nhkkh1+1
 
 6263       jdahkk(2,ihtaww)=nhkk
 
 6267      + 
WRITE(6,
'(A)')
'  exit from hadhad with irejfo=0 ' 
 6271       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6277       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 6377       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 6378      +iibar(210),k1(210),k2(210)
 
 6381       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 6383       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 6385       COMMON /chabai/chargi,barnui
 
 6386       COMMON /evappp/ievap
 
 6401         IF (isthkk(i).EQ.13)
THEN 
 6403       IF(idhkk(i).EQ.2212)chaeve=chaeve+1.d0
 
 6405         IF (isthkk(i).EQ.14)
THEN 
 6407       IF(idhkk(i).EQ.2212)chaeve=chaeve+1.d0
 
 6410       DO 521 i=nhkkh1,nhkk
 
 6411         IF (isthkk(i).EQ.1.OR.isthkk(i).EQ.15.OR.isthkk(i).EQ.16)
THEN 
 6414           IF (nrhkk.LE.0.OR.nrhkk.GT.410)
THEN 
 6415             WRITE(6,1389)nrhkk,i,idhkk(i),nhkkh1,nhkk
 
 6416  1389       
FORMAT (
' distr: NRHKK ERROR ',5i10)
 
 6421       chaeve=chaeve+ichhkk
 
 6425       ELSEIF(ievap.EQ.1)
THEN 
 6427         IF (isthkk(i).EQ.1)
THEN 
 6431       chaeve=chaeve+ichhkk
 
 6437         IF (isthkk(i).EQ.-1)
THEN 
 6438       IF(idhkk(i).EQ.2112)
THEN 
 6442       IF(idhkk(i).EQ.2212)
THEN 
 6448     IF((idhkk(i).EQ.80000).AND.(isthkk(i).NE.1000))
THEN 
 6449       chaeve=chaeve+idxres(i)
 
 6450       baeve=baeve+idres(i)
 
 6455       IF(ievl.LE.10)
WRITE(6,
'(2A,4F10.2)')
' Event charge and B-number',
 
 6456      * 
'=',chaeve,baeve,chargi,barnui
 
 6457       IF(chaeve-chargi.NE.0.d0.OR.baeve-barnui.NE.0.d0)
THEN 
 6459       IF(ievl.LE.1000)
WRITE(6,
'(2A,4F10.2)')
'Event charge and B-numb',
 
 6460      *
'(violated)  =',chaeve,baeve,chargi,barnui
 
 6471       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6473       dimension 
pt(50,10),ypt(50,10)
 
 6479           pt(j,i)=j*dpt-dpt/2.
 
 6486       IF(ipt1.GT.50)ipt1=50
 
 6487       IF(ipt2.GT.50)ipt2=50
 
 6488       ypt(ipt1,ipt)=ypt(ipt1,ipt)+1.
 
 6489       ypt(ipt2,ipt)=ypt(ipt2,ipt)+1.
 
 6490       ypt(ipt1,10)=ypt(ipt1,10)+1.
 
 6491       ypt(ipt2,10)=ypt(ipt2,10)+1.
 
 6496           ypt(j,i)=ypt(j,i)/nevt
 
 6497           ypt(j,i)=log10(ypt(j,i)+1.
d-18)
 
 6509       SUBROUTINE hkkfil(IST,ID,M1,M2,PX,PY,PZ,E,NHKKAU,KORMO,ICALL)
 
 6511       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6513       parameter(lout=6,llook=9)
 
 6514       parameter(tiny10=1.0
d-10,tiny4=1.0
d-3)
 
 6520       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 6607       COMMON /nncms/  gacms,bgcms,umo,pcm,eproj,pproj
 
 6608       COMMON /trafop/galab,bglab,blab
 
 6609       COMMON /projk/ iprojk
 
 6621          WRITE(lout,1000) nhkk
 
 6622  1000    
FORMAT(1
x,
'HKKFIL: NHKK exeeds NMXHKK = ',i7,
 
 6623      &             
'! program execution stopped..')
 
 6626       IF (m1.LT.0) mo1 = nhkk+m1
 
 6627       IF (m2.LT.0) mo2 = nhkk+m2
 
 6630       IF(kormo.EQ.999)
THEN 
 6631         jmohkk(1,nhkk) = mo1
 
 6632         jmohkk(2,nhkk) = mo2
 
 6634         jmohkk(1,nhkk)=nhkkau+kormo-1
 
 6640          IF (jdahkk(1,mo1).NE.0) 
THEN 
 6641             jdahkk(2,mo1) = nhkk
 
 6643             jdahkk(1,mo1) = nhkk
 
 6645      jdahkk(1,mo1)=nhkkau
 
 6648          IF (jdahkk(1,mo2).NE.0) 
THEN 
 6649             jdahkk(2,mo2) = nhkk
 
 6651             jdahkk(1,mo2) = nhkk
 
 6653          jdahkk(1,mo2) = nhkkau
 
 6659       phkk(5,nhkk) = phkk(4,nhkk)**2-phkk(1,nhkk)**2-
 
 6660      &               phkk(2,nhkk)**2-phkk(3,nhkk)**2
 
 6661       IF ((phkk(5,nhkk).LT.0.0d0).AND.(abs(phkk(5,nhkk)).GT.tiny4))
 
 6662      &   
WRITE(lout,
'(1X,A,G10.3)')
 
 6663      &     
'HKKFIL: negative mass**2 ',phkk(5,nhkk)
 
 6664       phkk(5,nhkk) = 
sqrt(abs(phkk(5,nhkk)))
 
 6665       IF (ist.EQ.88888.OR.ist.EQ.88887.OR.ist.EQ.88889) 
THEN 
 6671             vhkk(i,nhkk) = vhkk(i,mo2)
 
 6673          vhkk(4,nhkk) = vhkk(3,mo2)/blab-vhkk(3,mo1)/bglab
 
 6677             vhkk(i,nhkk) = vhkk(i,mo1)
 
 6678             IF (iprojk.EQ.1) 
THEN 
 6679               whkk(i,nhkk) = whkk(i,mo1)
 
 6692       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6698       common/
pyjets/
n,npad,k(4000,5),
p(4000,5),v(4000,5)
 
 6699       common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
 
 6700       common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
 
 6701       common/
pydat3/mdcy(500,3),mdme(4000,2),brat(4000),kfdp(4000,5)
 
 6703       CHARACTER chap*16,chan*16,chad(5)*16
 
 6707         WRITE(mstu(11),6800)
 
 6711         IF(mstu(2).NE.0) kfmax=mstu(2)
 
 6716         IF(kc.EQ.0) goto 220
 
 6717         IF(mstu(14).EQ.0.AND.kf.GT.100.AND.kc.LE.100) goto 220
 
 6718         IF(mstu(14).GT.0.AND.kf.GT.100.AND.
max(
mod(kf/1000,10),
 
 6719      &  
mod(kf/100,10)).GT.mstu(14)) goto 220
 
 6727         IF(kf.LE.100.AND.chap.EQ.
' '.AND.mdcy(kc,2).EQ.0) goto 220
 
 6732     idc2=mdcy(kc,2)+mdcy(kc,3)-1
 
 6733         WRITE(mstu(11),6900)kbam,
 
 6734      &  kf,kc,idc1,idc2,chap,chan,kchg(kc,1),kchg(kc,2),
 
 6735      &  kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
 
 6737      &  kf,kc,idc1,idc2,chap,chan,kchg(kc,1),kchg(kc,2),
 
 6738      &  kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
 
 6742         IF(kf.GT.100.AND.kc.LE.100) goto 220
 
 6743         DO 210 idc=mdcy(kc,2),mdcy(kc,2)+mdcy(kc,3)-1
 
 6746           CALL 
pyname(kfdp(idc,j),chad(j))
 
 6748       kbamdp(j)=
mcihad(kfdp(idc,j))
 
 6749       IF(kbamdp(j).EQ.26)kbamdp(j)=0
 
 6751         WRITE(26,7001) idc,mdme(idc,1),mdme(idc,2),brat(idc),
 
 6753   210   
WRITE(mstu(11),7000) idc,mdme(idc,1),mdme(idc,2),brat(idc),
 
 6756     IF(kabam.NE.410)
THEN 
 6757         WRITE(mstu(11),6900)kabam,
 
 6758      &  -kf,-kc,idc1,idc2,chan,chap,-kchg(kc,1),kchg(kc,2),
 
 6759      &  kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
 
 6760         WRITE(26,6900)kabam,
 
 6761      &  -kf,-kc,idc1,idc2,chan,chap,-kchg(kc,1),kchg(kc,2),
 
 6762      &  kchg(kc,3),pm,pmas(kc,2),pmas(kc,3),pmas(kc,4),mdcy(kc,1)
 
 6763         DO 211 idc=mdcy(kc,2),mdcy(kc,2)+mdcy(kc,3)-1
 
 6767     IF(kcdp.LE.0.OR.kcdp.GT.500)
THEN 
 6773       IF(kchg(kcdp,3).EQ.0)kfdpm=kfdp(idc,j)
 
 6775       IF(kbamdp(j).EQ.26)kbamdp(j)=0
 
 6777           CALL 
pyname(kfdpm,chad(j))
 
 6779         WRITE(26,7001) idc,mdme(idc,1),mdme(idc,2),brat(idc),
 
 6781   211   
WRITE(mstu(11),7000) idc,mdme(idc,1),mdme(idc,2),brat(idc),
 
 6789  6800 
FORMAT(///30
x,
'Particle/parton data table'//1
x,
'BAM',
 
 6790      &1
x,
'ABAM',1
x,
'KF',1
x,
'KC',1
x,
'DCF',1
x,
'DCL',1
x,
 
 6791      &
'particle',8
x,
'antiparticle',6
x,
'chg  col  anti',8
x,
'mass',7
x,
 
 6792      &
'width',7
x,
'w-cut',5
x,
'lifetime',1
x,
'decay'/11
x,
'IDC',1
x,
'on/off',
 
 6793      &1
x,
'ME',3
x,
'Br.rat.',4
x,
'decay products')
 
 6794  6900 
FORMAT(/1
x,i4,i6,i4,2i5,a16,a16,3i3,1
x,f12.5,2(1
x,f11.5),
 
 6796  7000 
FORMAT(10
x,i4,2
x,i3,2
x,i3,2
x,f8.5,4
x,5a16)
 
 6797  7001 
FORMAT(10
x,i4,2
x,i3,2
x,i3,2
x,f8.5,4
x,5i5)