1       SUBROUTINE saptre(AM1,G1,BGX1,BGY1,BGZ1,
 
    2      &                  am2,g2,bgx2,bgy2,bgz2)
 
    3       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
   19       esmax=min(esmax1,esmax2)
 
   20       IF(esmax.LE.0.05d0) 
RETURN 
   22       IF (b3*esmax.GT.60.d0)
THEN 
   28       axexp=(1.d0-(b3*esmax-1.d0)*exeb)/b3**2
 
   35         es=-2./(b3**2)*
log(
x*
y+1.
e-7)
 
   38         es=abs(-
log(
x+1.
e-7)/b3)
 
   40       IF(es.GT.esmax)                                             goto10
 
   42       hps=
sqrt((es-hma)*(es+hma))
 
   49       pz1nsq=pz1**2-hps**2-2.*px1*hpx-2.*py1*hpy
 
   50       pz2nsq=pz2**2-hps**2+2.*px2*hpx+2.*py2*hpy
 
   51       IF(pz1nsq.LT.0.001d0.OR.pz2nsq.LT.0.001d0) 
RETURN 
   52       pz1=sign(
sqrt(pz1nsq),pz1)
 
   53       pz2=sign(
sqrt(pz2nsq),pz2)
 
   72       SUBROUTINE sltraf(GA,BGA,EIN,PZIN,EOUT,PZOUT)
 
   73       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
   75       pzout=ga*pzin - bga*ein
 
   76       eout=ga*ein - bga*pzin
 
   84       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
   96       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
   98      *                ,xpsu(248),xtsu(248)
 
   99      *                ,xpsut(248),xtsut(248)
 
  101       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
  102      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
  103      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
  105      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
  119       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
  125      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
  127       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
  128       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
  135       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
  138      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
  143       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
  231       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
  233       COMMON /nncms/  gamcm,bgcm,umo,pcm,eproj,pproj
 
  235       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
  247       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
  248      +iibar(210),k1(210),k2(210)
 
  251       COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
 
  252      +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
 
  253      +prebin,taebin,fermod,etacou
 
  255       COMMON /projk/ iprojk
 
  257        IF(ijproj.EQ.5)
RETURN 
  268           CALL 
sltraf(gaproj,-bgproj, phkk(4,j),phkk(3,j),prmom4,prmom3)
 
  270           CALL 
sltraf(gamcm,+bgcm, prmom4,prmom3,prmom(4,j),prmom(3,j))
 
  272           prmom(5,j)=
sqrt( abs((prmom(4,j)-aam(kk)) *(prmom(4,j)+aam(kk)
 
  284           tamom(1,j)=phkk(1,ihkk)
 
  285           tamom(2,j)=phkk(2,ihkk)
 
  286           CALL 
sltraf(gamcm,bgcm, phkk(4,ihkk),phkk(3,ihkk),tamom(4,j),
 
  288           tamom(5,j)=
sqrt(abs( (tamom(4,j)-aam(kk)) 
 
  289      +    *(tamom(4,j)+aam(kk))))
 
  295         WRITE(6,
'(/A,I5/5X,A)') 
' NUCMOM: IP=',ip,
 
  296      +  
' J,IPVQ(J),IPPV1(J),IPPV2(J),ISTHKK,KKPROJ,PRMOM' 
  298           WRITE(6,
'(I4,5I3,5(1PE11.3))') j,isthkk(j),kkproj(j), ipvq(j),
 
  299      +    ippv1(j),ippv2(j), (prmom(jj,j),jj=1,5)
 
  303         WRITE(6,
'(/A,I5/5X,A)') 
' NUCMOM: IT=',it,
 
  304      +  
' J,ITVQ(J),ITTV1(J),ITTV2(J),ISTHKK,KKTARG,TAMOM' 
  308           WRITE(6,
'(I4,5I3,5(1PE11.3))') j,isthkk(ihkk),kktarg(j), itvq
 
  309      +    (j),ittv1(j),ittv2(j), (tamom(jj,j),jj=1,5)
 
  320       SUBROUTINE fer4m(PFERM,PXT,PYT,PZT,ET,KT)
 
  321       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
  337       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
  338      +iibar(210),k1(210),k2(210)
 
  341       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
  343       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
  344      +ipadis,ishmal,lpauli
 
  351         CALL 
dpoli(polc,pols)
 
  357         et=
sqrt(pabs*pabs+aam(kt)**2)
 
  373       SUBROUTINE fer4mp(IP,PFERM,PXT,PYT,PZT,ET,KT)
 
  374       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
  376       COMMON /ferfor/iferfo
 
  391       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
  392      +iibar(210),k1(210),k2(210)
 
  395       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
  397       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
  398      +ipadis,ishmal,lpauli
 
  406     IF(iferfo.EQ.2)CALL 
dfatpr(ip,pabs)
 
  408         CALL 
dpoli(polc,pols)
 
  414         et=
sqrt(pabs*pabs+aam(kt)**2)
 
  430       SUBROUTINE fer4mt(IT,PFERM,PXT,PYT,PZT,ET,KT)
 
  431       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
  433       COMMON /ferfor/iferfo
 
  448       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
  449      +iibar(210),k1(210),k2(210)
 
  452       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
  454       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
  455      +ipadis,ishmal,lpauli
 
  466     IF(iferfo.EQ.2)CALL 
dfatta(it,pabs)
 
  469         CALL 
dpoli(polc,pols)
 
  475         et=
sqrt(pabs*pabs+aam(kt)**2)
 
  490       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
  493       dimension par10(6),par20(6),par30(6),par40(6),par50(6),
 
  494      *          par60(6),par11(6),par21(6),par31(6),par41(6),
 
  495      *          aia(6),att(101),catt(101),aka(101)
 
  496       common/fattad/daka(101),fatt(101)
 
  497       DATA par10/1.61d0,2.74d0,3.24d0,3.57d0,1.80d0,0.d0/
 
  498       DATA par20/2.66d0,3.33d0,3.72d0,4.97d0,4.77d0,0.d0/
 
  499       DATA par30/3.54d0,6.66d0,0.d0,0.d0,0.d0,0.d0/
 
  500       DATA par40/0.d0,0.d0,11.1d0,19.8d0,25.5d0,0.d0/
 
  501       DATA par50/0.d0,0.d0,0.d0,15.d0,0.d0,0.d0/
 
  502       DATA par60/0.d0,0.d0,0.d0,0.d0,40.3d0,0.d0/
 
  503       DATA par11/.426d0,.326d0,.419d0,.230d0,.275d0,0.d0/
 
  504       DATA par21/1.6d0,1.4d0,1.77d0,1.2d0,1.01d0,0.d0/
 
  505       DATA par31/.0237d0,.0263d0,.0282d0,.0286d0,.0304d0,0.d0/
 
  506       DATA par41/.22d0,.22d0,.22d0,.22d0,.22d0,0.d0/
 
  507       DATA aia/12.d0,16.d0,40.d0,56.d0,208.d0,209.d0/
 
  514     IF(ait.GE.aia(i).AND.ait.LT.aia(i+1))
THEN 
  515       dait=(ait-aia(i))/(aia(i+1)-aia(i))
 
  520       IF(ait.LT.aia(1))
THEN 
  525       IF(ait.GE.aia(5))
THEN 
  530       a0=dbit*par10(iii)+dait*par10(iii+1)
 
  531       b0=dbit*par20(iii)+dait*par20(iii+1)
 
  532       c0=dbit*par30(iii)+dait*par30(iii+1)
 
  533       d0=dbit*par40(iii)+dait*par40(iii+1)
 
  534       e0=dbit*par50(iii)+dait*par50(iii+1)
 
  535       f0=dbit*par60(iii)+dait*par60(iii+1)
 
  536       a1=dbit*par11(iii)+dait*par11(iii+1)
 
  537       b1=dbit*par21(iii)+dait*par21(iii+1)
 
  538       c1=dbit*par31(iii)+dait*par31(iii+1)
 
  539       d1=dbit*par41(iii)+dait*par41(iii+1)
 
  548     att(i)=ak**2*(a0*
exp(-b0*ak**2)*(1.d0+c0*ak**2+
 
  549      *         d0*ak**4+e0*ak**6+f0*ak**8)+
 
  550      *         a1*
exp(-b1*ak**2)+
c1*
exp(-d1*ak**2))
 
  551     IF(i.GT.1)catt(i)=catt(i-1)+att(i)
 
  554     catt(i)=catt(i)/catt(101)
 
  561     IF(rndfa.LT.catt(i))
THEN 
  567       pabs=aka(iatt)*0.197d0
 
  568       fatt(iatt)=fatt(iatt)+1.d0/pabs**2
 
  572       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
  575       dimension par10(6),par20(6),par30(6),par40(6),par50(6),
 
  576      *          par60(6),par11(6),par21(6),par31(6),par41(6),
 
  577      *          aia(6),att(101),catt(101),aka(101)
 
  578       DATA par10/1.61d0,2.74d0,3.24d0,3.57d0,1.80d0,0.d0/
 
  579       DATA par20/2.66d0,3.33d0,3.72d0,4.97d0,4.77d0,0.d0/
 
  580       DATA par30/3.54d0,6.66d0,0.d0,0.d0,0.d0,0.d0/
 
  581       DATA par40/0.d0,0.d0,11.1d0,19.8d0,25.5d0,0.d0/
 
  582       DATA par50/0.d0,0.d0,0.d0,15.d0,0.d0,0.d0/
 
  583       DATA par60/0.d0,0.d0,0.d0,0.d0,40.3d0,0.d0/
 
  584       DATA par11/.426d0,.326d0,.419d0,.230d0,.275d0,0.d0/
 
  585       DATA par21/1.6d0,1.4d0,1.77d0,1.2d0,1.01d0,0.d0/
 
  586       DATA par31/.0237d0,.0263d0,.0282d0,.0286d0,.0304d0,0.d0/
 
  587       DATA par41/.22d0,.22d0,.22d0,.22d0,.22d0,0.d0/
 
  588       DATA aia/12.d0,16.d0,40.d0,56.d0,208.d0,209.d0/
 
  595     IF(ait.GE.aia(i).AND.ait.LT.aia(i+1))
THEN 
  596       dait=(ait-aia(i))/(aia(i+1)-aia(i))
 
  601       IF(ait.LT.aia(1))
THEN 
  606       IF(ait.GE.aia(5))
THEN 
  611       a0=dbit*par10(iii)+dait*par10(iii+1)
 
  612       b0=dbit*par20(iii)+dait*par20(iii+1)
 
  613       c0=dbit*par30(iii)+dait*par30(iii+1)
 
  614       d0=dbit*par40(iii)+dait*par40(iii+1)
 
  615       e0=dbit*par50(iii)+dait*par50(iii+1)
 
  616       f0=dbit*par60(iii)+dait*par60(iii+1)
 
  617       a1=dbit*par11(iii)+dait*par11(iii+1)
 
  618       b1=dbit*par21(iii)+dait*par21(iii+1)
 
  619       c1=dbit*par31(iii)+dait*par31(iii+1)
 
  620       d1=dbit*par41(iii)+dait*par41(iii+1)
 
  628     att(i)=ak**2*(a0*
exp(-b0*ak**2)*(1.d0+c0*ak**2+
 
  629      *         d0*ak**4+e0*ak**6+f0*ak**8)+
 
  630      *         a1*
exp(-b1*ak**2)+
c1*
exp(-d1*ak**2))
 
  631     IF(i.GT.1)catt(i)=catt(i-1)+att(i)
 
  634     catt(i)=catt(i)/catt(101)
 
  640     IF(rndfa.LT.catt(i))
THEN 
  646       pabs=aka(iatt)*0.197d0
 
  654       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
  665       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
  667      *                ,xpsu(248),xtsu(248)
 
  668      *                ,xpsut(248),xtsut(248)
 
  670       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
  671      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
  672      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
  674      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
  688       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
  694      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
  696       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
  697       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
  704       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
  707      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
  712       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
  800       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
  812       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
  813      +iibar(210),k1(210),k2(210)
 
  816       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
  818       COMMON /projk/ iprojk
 
  821       dimension ihkkq(-6:6),ihkkqq(-3:3,-3:3)
 
  822       DATA ihkkq/-6,-5,-4,-3,-1,-2,0,2,1,3,4,5,6/
 
  823       DATA ihkkqq/-3303,-3103,-3203,0,0,0,0, -3103,-1103,-2103,0,0,0,0,
 
  824      +-3203,-2103,-2203,0,0,0,0, 0,0,0,0,0,0,0, 0,0,0,0,2203,2103,3203,
 
  825      +0,0,0,0,2103,1103,3103, 0,0,0,0,3203,3103,3303/
 
  832       IF(ipev.GE.3) 
WRITE(6,
'(A,6I4)')
 
  833      +
' FLKSAM-ENTRY: IT,ITZ, IP,IPZ, IJPROJ,IBPROJ', it,itz,ip,ipz,
 
  843         CALL 
flahad(kproj,ibproj,ipvq(jp),ippv1(jp),ippv2(jp))
 
  845         IF (ipev.GE.6) 
WRITE (6,1000)ipvq(jp),ippv1(jp),ippv2(jp)
 
  846  1000 
FORMAT (
' FLKSAM: IPVQ,IPPV1,IPPV2 ',3i4)
 
  850         idhkk(jhkkq)=ihkkq(ipvq(jp))
 
  852           idhkk(jhkk)=ihkkq(ippv1(jp))
 
  854           idhkk(jhkk)=ihkkqq(ippv1(jp),ippv2(jp))
 
  869         idhkk(jhkkq)=ihkkq(ipsq(
n))
 
  870         idhkk(jhkkaq)=ihkkq(ipsaq(
n))
 
  883         CALL 
flahad(ktarg,ibtarg,itvq(jt),ittv1(jt),ittv2(jt))
 
  887         idhkk(jhkkq)=ihkkq(itvq(jt))
 
  889           idhkk(jhkk)=ihkkq(ittv1(jt))
 
  891           idhkk(jhkk)=ihkkqq(ittv1(jt),ittv2(jt))
 
  893         IF (ipev.GE.8) 
WRITE (6,
'(A,8I4)')
 
  894      +  
' FLKSAM: KTARG,ITVQ(JT),ITTV1(JT),ITTV2(JT)', ktarg,itvq(jt),
 
  895      +  ittv1(jt),ittv2(jt), idhkk(jhkkq),jhkkq,idhkk(jhkk),jhkk
 
  911         idhkk(jhkkq)=ihkkq(itsq(
n))
 
  912         idhkk(jhkkaq)=ihkkq(itsaq(
n))
 
  923       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
  935       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
  937      *                ,xpsu(248),xtsu(248)
 
  938      *                ,xpsut(248),xtsut(248)
 
  940       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
  941      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
  942      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
  944      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
  958       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
  964      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
  966       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
  967       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
  974       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
  977      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
  982       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 1070       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 1082       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 1083      +iibar(210),k1(210),k2(210)
 
 1086       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 1088       COMMON /projk/ iprojk
 
 1089       COMMON /seasu3/seasq 
 
 1091       parameter(ummm=0.3d0)
 
 1092       parameter(smmm=0.5d0)
 
 1093       parameter(cmmm=1.3d0)
 
 1121       betcha=betoo+1.3-log10(ecm)
 
 1140         WRITE(6,4567)pc,betcha,pu1,ps1
 
 1141  4567   
FORMAT(
' Charm at chain ends FLKSAA: PC,BETCHA,PU,PS ',4f10.5)
 
 1145       IF(ipev.GE.3) 
WRITE(6,
'(A,6I4)')
 
 1146      +
' FLKSAA-ENTRY: IT,ITZ, IP,IPZ, IJPROJ,IBPROJ', it,itz,ip,ipz,
 
 1161         is=1.d0+
rndm(v1)*(2.d0+seasq)
 
 1165         IF (ipev.GE.8) 
WRITE (6,1010) 
n,ipsq(
n),ipsaq(
n)
 
 1166  1010 
FORMAT (
' FLKSAA: N,IPSQ(N),IPSAQ(N) ',3i4)
 
 1179         is=1.d0+
rndm(v1)*(2.d0+seasq) 
 
 1183         IF (ipev.GE.8) 
WRITE (6,1020) 
n,itsq(
n),itsaq(
n)
 
 1184  1020 
FORMAT (
' FLKSAA: N,ITSQ(N),ITSAQ(N) ',3i4)
 
 1195       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 1205       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 1207       dimension mquark(3,30)
 
 1208       DATA mquark/ 2,1,1, -2,-1,-1, 0,0,0, 0,0,0, 0,0,0, 0,0,0, 0,0,0,
 
 1209      +2,2,1, -2,-2,-1, 0,0,0, 0,0,0, 0,0,0, 1,-2,0, 2,-1,0, 1,-3,0, 3,
 
 1210      +-1,0, 1,2,3, -1,-2,-3, 0,0,0, 2,2,3, 1,1,3, 1,2,3, 1,-1,0, 2,-3,0,
 
 1211      +3,-2,0, 2,-2,0, 3,-3,0, 0,0,0, 0,0,0, 0,0,0/
 
 1214         ipq1 = mquark(1,ityp)
 
 1215         ipq2 = mquark(2,ityp)
 
 1216         ipq3 = mquark(3,ityp)
 
 1218         IF(ipev.GE.3) 
print 1000, ityp,ibar
 
 1219  1000 
FORMAT(
' FLAHAD: ITYP,IBAR',2i5)
 
 1221         isam5=1. + 6.*
rndm(v)
 
 1222         go to(10,20,30,40,50,60),isam5
 
 1253         IF (ipev.GE.3) 
WRITE (6,1010) if1,if2,if3
 
 1254  1010 
FORMAT (
' FLAHAD: IF1,IF2,IF3 ',3i4)
 
 1261           WRITE(6,
'(A,6I4)') 
' FLAHAD (MESON): IF1,IF2,IF3', if1,if2,if3
 
 1274       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 1280       COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
 
 1281      +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
 
 1284       parameter(amis=0.8d0,amas=2.6d0,amiu=0.5d0,amau=2.6d0)
 
 1288       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 1290      *                ,xpsu(248),xtsu(248)
 
 1291      *                ,xpsut(248),xtsut(248)
 
 1293       COMMON /intnez/ndz,nzd
 
 1294       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 1295      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 1296      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 1298      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 1312       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 1318      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 1320       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 1321       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 1328       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 1331      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 1336       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 1428       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 1440       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 1441      +iibar(210),k1(210),k2(210)
 
 1444       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 1449       COMMON /projk/ iprojk
 
 1451       COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
 
 1455       COMMON /abrvv/ amcvv1(248),amcvv2(248),gacvv1(248),gacvv2(248),
 
 1456      +bgxvv1(248),bgyvv1(248),bgzvv1(248), bgxvv2(248),bgyvv2(248),
 
 1457      +bgzvv2(248), nchvv1(248),nchvv2(248),ijcvv1(248),ijcvv2(248),
 
 1458      +pqvva1(248,4),pqvva2(248,4), pqvvb1(248,4),pqvvb2(248,4)
 
 1467       COMMON /abrvs/ amcvs1(248),amcvs2(248),gacvs1(248),gacvs2(248),
 
 1468      +bgxvs1(248),bgyvs1(248),bgzvs1(248), bgxvs2(248),bgyvs2(248),
 
 1469      +bgzvs2(248), nchvs1(248),nchvs2(248),ijcvs1(248),ijcvs2(248),
 
 1470      +pqvsa1(248,4),pqvsa2(248,4), pqvsb1(248,4),pqvsb2(248,4)
 
 1472       COMMON /abrsv/ amcsv1(248),amcsv2(248),gacsv1(248),gacsv2(248),
 
 1473      +bgxsv1(248),bgysv1(248),bgzsv1(248), bgxsv2(248),bgysv2(248),
 
 1474      +bgzsv2(248), nchsv1(248),nchsv2(248),ijcsv1(248),ijcsv2(248),
 
 1475      +pqsva1(248,4),pqsva2(248,4), pqsvb1(248,4),pqsvb2(248,4)
 
 1478       COMMON /seadiq/lseadi
 
 1480       common/diquax/amedd,idiqua,idiquu
 
 1481       COMMON /xsvthr/ xsthr,xvthr,xdthr,xssthr
 
 1482       COMMON /seaqxx/ seaqx,seaqxn 
 
 1483       dimension isxpvq(248),isxpvd(248),isxtvq(248),isxtvd(248)
 
 1484       parameter(sqma=0.7d0)
 
 1491       parameter(nsea=3,nval=10)
 
 1508         IF (i.GT.248)                                           go to 30
 
 1512       IF(ecm.LE.1.
d-3)
THEN 
 1513     WRITE(6,*)
' xksamp: ECM=0.D0 ' 
 1517       IF(xsthr.LE.1.
d-12)
THEN 
 1518     WRITE(6,*)
' xksamp 30 : XSTHR=0.D0 ',csea,ecm,xsthr
 
 1528       IF(ip.EQ.1) xsthr=4./ecm**2
 
 1531       IF(xsthr.LE.1.
d-12)
THEN 
 1540       IF(ip.GE.150.AND.it.GE.150) xsthr=2.5/(ecm*
sqrt(ecm))
 
 1546       IF (ecm.LT.10.d0.AND.ip.GT.1)xsthr=((12.-ecm)/5.+1.)*csea/ecm
 
 1549       IF (xvthr+xdthr.GT.0.90d0)
THEN 
 1551     IF(xvthr.LE.0.05d0)
THEN 
 1555       IF(ecm.LE.1.
d-3)
THEN 
 1556     WRITE(6,*)
' xksamp: ECM=0.D0 ' 
 1561       IF(jcoun.EQ.1)
WRITE(6,
'(A,4E15.5)')
 
 1562      *
' XKSAMP: XSTHR,XVTHR,XDTHR,XSSTHR ',
 
 1563      * xsthr,xvthr,xdthr,xssthr 
 
 1564       IF(ipev.GE.1)
WRITE(6,
'(A,4E15.5)')
 
 1565      *
' XKSAMP: XSTHR,XVTHR,XDTHR,XSSTHR ',
 
 1566      * xsthr,xvthr,xdthr,xssthr 
 
 1569       IF (xvthr+xdthr.GT.0.95d0)
THEN 
 1571  1000   
FORMAT (
' PROGRAMM STOPPED IN XSAMP1 ECM = ',f6.2,
' TOO SMALL')
 
 1578       xxseam=1.0 - xvthr*(1.d0+0.3d0*
rndm(v1)) 
 
 1579      *            - xdthr*(1.d0+0.3d0*
rndm(v2))
 
 1580      *             -0.01*(1.d0+1.5d0*
rndm(v3))
 
 1583       IF(seaqxn.GE.0.75d0)
THEN 
 1587     xxseam=1.d0-xvthr-xdthr
 
 1589         xxseam=1.0 - xvthr*(1.d0+
rndm(v1)) - xdthr*(1.d0+
rndm(v2))
 
 1590      *             -0.01*(1.d0+5.d0*
rndm(v3))
 
 1593       IF(xsthr.LE.1.
d-12)
THEN 
 1596           WRITE(6,*)
' xksamp: XSTHR=0.D0 ' 
 1597           WRITE(6,
'(A,2E20.5,I10)')
 
 1598      *            
' XXSEAM,XSTHR,NSMAX',xxseam,xsthr,nsmax
 
 1602       nsmax=0.50*xxseam / xsthr
 
 1603       IF(ipev.GE.1)
WRITE(6,
'(A,E15.5,I10)')
 
 1604      *            
' XXSEAM,NSMAX',xxseam,nsmax
 
 1609       IF (xdthr.GT.0.14d0)xdthr=0.14d0
 
 1610       IF (xvthr.GT.0.14d0)xvthr=0.14d0
 
 1617       IF(ijproj.NE.0.AND.ibproj.EQ.0) unoprv=unom
 
 1618       IF(jcoun.EQ.1)
WRITE(6,
'(A,4E15.5)')
 
 1619      *
' XKSAMP: XSTHR,XVTHR,XDTHR,XSSTHR ',
 
 1620      * xsthr,xvthr,xdthr,xssthr 
 
 1623         IF (jssh(ipp).NE.0) 
THEN 
 1634           jipp=min(jipp,nsmax)
 
 1639         xsmax=xxseam - 1.5*jipp*xsthr
 
 1641             IF(xsthr.GE.xsmax) 
THEN 
 1648             IF(ipev.GE.1)
WRITE(6,
'(A)') 
' XKSAMP-40' 
 1651             IF (nscoun.GT.nsea) 
THEN 
 1657               IF(ipsq(ixps+1).LE.2)
THEN 
 1660                 IF(seaqxn.LE.0.75d0)
THEN 
 1661                   xpsqi=
sampex(xsthr,xsmax)
 
 1663                 ELSEIF(seaqxn.GT.0.75d0)
THEN 
 1664                   xpsqi=
sampey(xsthr,xsmax)
 
 1667         IF(ipev.GE.1)
WRITE(6,
'(A,3E15.5)')
 
 1668      *          
'XPSQI 1:XPSQI,XSTHR,XSMAX',
 
 1671                 IF(xsmax.GT.xsthr+bsqma)
THEN 
 1672                   xpsqi=
sampxb(xsthr+bsqma,xsmax,bsqma)
 
 1673               IF(ipev.GE.1)
WRITE(6,
'(A,4E15.5)')
 
 1674      *            
'XPSQI 2:XPSQI,XSTHR,XSMAX,BSQMA',
 
 1675      *            xpsqi,xsthr,xsmax,bsqma
 
 1679                   IF(seaqxn.LE.0.75d0)
THEN 
 1680                     xpsqi=
sampex(xsthr,xsmax)
 
 1682                   ELSEIF(seaqxn.GT.0.75d0)
THEN 
 1683                     xpsqi=
sampey(xsthr,xsmax)
 
 1686               IF(ipev.GE.1)
WRITE(6,
'(A,3E15.5)')
 
 1687      *            
'XPSQI 3:XPSQI,XSTHR,XSMAX',
 
 1692               IF(ipsaq(ixps+1).GE.-2)
THEN 
 1695                 IF(seaqxn.LE.0.75d0)
THEN 
 1696                   xpsaqi=
sampex(xsthr,xsmax)
 
 1698                 ELSEIF(seaqxn.GT.0.75d0)
THEN 
 1699                   xpsaqi=
sampey(xsthr,xsmax)
 
 1702         IF(ipev.GE.1)
WRITE(6,
'(A,3E15.5)')
 
 1703      *          
'XPSAQI 1:XPSAQI,XSTHR,XSMAX',
 
 1704      *          xpsaqi,xsthr,xsmax
 
 1706                 IF(xsmax.GT.xsthr+bsqma)
THEN 
 1707                   xpsaqi=
sampxb(xsthr+bsqma,xsmax,bsqma)
 
 1708           IF(ipev.GE.1)
WRITE(6,
'(A,4E15.5)')
 
 1709      *            
'XPSAQI 2:XPSAQI,XSTHR,XSMAX,BSQMA',
 
 1710      *            xpsaqi,xsthr,xsmax,bsqma
 
 1714                   IF(seaqxn.LE.0.75d0)
THEN 
 1715                     xpsaqi=
sampex(xsthr,xsmax)
 
 1717                   ELSEIF(seaqxn.GT.0.75d0)
THEN 
 1718                     xpsaqi=
sampey(xsthr,xsmax)
 
 1721           IF(ipev.GE.1)
WRITE(6,
'(A,3E15.5)')
 
 1722      *            
'XPSAQI 3:XPSAQI,XSTHR,XSMAX',
 
 1723      *            xpsaqi,xsthr,xsmax
 
 1729      *        
WRITE(6,
'(A,3E15.4)') 
' XKSAMP-50: XPSQI,XSTHR,XSMAX',
 
 1732               IF(ipev.GE.1)
WRITE(6,
'(A)') 
' XKSAMP-60' 
 1734      *        
WRITE(6,
'(A,3E15.4)') 
' XKSAMP-60: XPSAQI,XSTHR,XSMAX',
 
 1735      &        xpsaqi,xsthr,xsmax
 
 1736               xxsea=xxsea + xpsqi + xpsaqi
 
 1737               IF(xxsea.GE.xxseam) 
THEN 
 1742               IF(ipev.GE.1)
WRITE(6,
'(A,I10)') 
' XKSAMP-60: IXPS',ixps
 
 1755           IF(xvthr.GT.0.05d0)
THEN 
 1756             IF(xvthr.GT.1.d0-xxsea-xdthr)
THEN 
 1757               WRITE(6,*)
' xvthr,xxsea,xdthr ', xvthr,xxsea,xdthr
 
 1761             xpvqi=
betrej(0.1d0,unoprv,xvthr,1.d0-xxsea-xdthr)
 
 1765             IF(ipev.GE.1)
WRITE(6,
'(A)') 
' XKSAMP-90' 
 1768             xpvqi=
dbetar(0.1d0,unoprv)
 
 1769             IF ((xpvqi.LT.xvthr).OR.(1.d0-xpvqi-xxsea.LT.xdthr))
 
 1772           xpvdi=1. - xpvqi - xxsea
 
 1775           IF(xpvdi.LT.xdthr) 
THEN 
 1776             WRITE(6,
'(A/A/E12.3,4I4,3E11.3)')
 
 1777      +      
' INCONSISTENT X-SAMPLING / XKSAMP / PROJECTILE',
 
 1778      +      
' ECM, IP, IPP, JSSH(IPP), JIPP, XPVQI, XPVDI, XXSEA', ecm,
 
 1779      +      ip, ipp, jssh(ipp), jipp, xpvqi, xpvdi, xxsea
 
 1806          IF (jtsh(itt).NE.0) 
THEN 
 1817            jitt=min(jitt,nsmax)
 
 1822          xsmax=xxseam -1.5*jitt*xsthr
 
 1824              IF(xsthr.GE.xsmax) 
THEN 
 1830              IF(ipev.GE.1)
WRITE(6,
'(A)') 
' XKSAMP-110' 
 1833              IF (nscoun.GT.nsea)
THEN 
 1841                IF(itsq(ixts+1).LE.2)
THEN 
 1844                  IF(seaqxn.LE.0.75d0)
THEN 
 1845                    xtsqi=
sampex(xsthr,xsmax)
 
 1847                  ELSEIF(seaqxn.GT.0.75d0)
THEN 
 1848                    xtsqi=
sampey(xsthr,xsmax)
 
 1852                IF(xsmax.GT.xsthr+bsqma)
THEN 
 1853                   xtsqi=
sampxb(xsthr+bsqma,xsmax,bsqma)
 
 1857       IF(seaqxn.LE.0.75d0)
THEN 
 1858                 xtsqi=
sampex(xsthr,xsmax)
 
 1860       ELSEIF(seaqxn.GT.0.75d0)
THEN 
 1861                 xtsqi=
sampey(xsthr,xsmax)
 
 1867               IF(itsaq(ixts+1).GE.-2)
THEN 
 1870       IF(seaqxn.LE.0.75d0)
THEN 
 1871                 xtsaqi=
sampex(xsthr,xsmax)
 
 1873       ELSEIF(seaqxn.GT.0.75d0)
THEN 
 1874                 xtsaqi=
sampey(xsthr,xsmax)
 
 1878                 IF(xsmax.GT.xsthr+bsqma)
THEN 
 1879                   xtsaqi=
sampxb(xsthr+bsqma,xsmax,bsqma)
 
 1883       IF(seaqxn.LE.0.75d0)
THEN 
 1884                 xtsaqi=
sampex(xsthr,xsmax)
 
 1886       ELSEIF(seaqxn.GT.0.75d0)
THEN 
 1887                 xtsaqi=
sampey(xsthr,xsmax)
 
 1898          IF(ipev.GE.1)
WRITE(6,
'(A)') 
' XKSAMP-120' 
 1902          IF(ipev.GE.1)
WRITE(6,
'(A)') 
' XKSAMP-130' 
 1906               xxsea=xxsea + xtsqi + xtsaqi
 
 1907               IF(xxsea.GE.xxseam) 
THEN 
 1912          IF(ipev.GE.1)
WRITE(6,
'(A,I10)')
' XKSAMP-130: IXTS',ixts
 
 1923           IF(xvthr.GT.0.05d0)
THEN 
 1924             IF(xvthr.GT.1.d0-xxsea-xdthr)
THEN 
 1925               WRITE(6,*)
' xvthr,xxsea,xdthr ', xvthr,xxsea,xdthr
 
 1929             xtvqi=
betrej(0.1d0,unon,xvthr,1.-xxsea-xdthr)
 
 1933          IF(ipev.GE.1)
WRITE(6,
'(A)') 
' XKSAMP-160' 
 1937         xmist=1.-xtvqi-xxsea
 
 1938         IF(ipev.GE.1)
WRITE(6,
'(A,5E15.5)')
 
 1939      *      
' XTVQI,XVTHR,XXSEA,XMIST,XDTHR',
 
 1940      *        xtvqi,xvthr,xxsea,xmist,xdthr
 
 1941           IF((xtvqi.LT.xvthr).OR.(1.d0-xtvqi-xxsea.LT.xdthr+0.0001d0))
 
 1944           xtvdi=1. - xtvqi - xxsea
 
 1947           IF(xtvdi.LT.xdthr) 
THEN 
 1948             WRITE(6,
'(A/A/E12.3,4I4,3E11.3)')
 
 1949      +      
' INCONSISTENT X-SAMPLING / XKSAMP / TARGET',
 
 1950      +      
' ECM, IT, ITT, JTSH(ITT), JITT, XTVQI, XTVDI, XXSEA', ecm,
 
 1951      +      it, itt, jtsh(itt), jitt, xtvqi, xtvdi, xxsea
 
 1977  1010 
FORMAT(
' XKSAMP:',
 
 1978      +
' I,XPVQ(I),XPVD(I),IFROVP(I),ITOVP(I),ZUOVP(I),KKPROJ(I)')
 
 1980           WRITE(6,1020) i,xpvq(i),xpvd(i),ifrovp(i),itovp(i),zuovp(i),
 
 1982  1020 
FORMAT(i5,2e15.5,2i5,l5,i5)
 
 1985  1030 
FORMAT(
' XKSAMP :  I,XPSQ(I),XPSAQ(I),IFROSP(I),ZUOSP(I)')
 
 1987           WRITE(6,1040) i,xpsq(i),xpsaq(i),ifrosp(i),zuosp(i)
 
 1988  1040 
FORMAT(i5,2e15.5,i5,l5)
 
 1992  1050 
FORMAT(
' XKSAMP:',
 
 1993      +
' I,XTVQ(I),XTVD(I),IFROVT(I),ITOVT(I),ZUOVT(I),KKTARG(I)')
 
 1995           WRITE(6,1020) i,xtvq(i),xtvd(i),ifrovt(i),itovt(i),zuovt(i),
 
 1999  1060 
FORMAT(
' XKSAMP :  I,XTSQ(I),XTSAQ(I),IFROST(I),ZUOST(I)')
 
 2001           WRITE(6,1040) i,xtsq(i),xtsaq(i),ifrost(i),zuost(i)
 
 2006      +  
' XKSAMP :  I,ITOVP(I),ITOVT(I),JSSHS(I),JTSHS(I)' 
 2009           WRITE(6,1070) i,itovp(i),itovt(i),jsshs(i),jtshs(i)
 
 2013       WRITE(6,*)
' I,IPVQ(I),IPPV1(I),IPPV2(I) ',
 
 2014      *    i,ipvq(i),ippv1(i),ippv2(i)
 
 2017       WRITE(6,*)
' I,ITVQ(I),ITTV1(I),ITTV2(I) ',
 
 2018      *    i,itvq(i),ittv1(i),ittv2(i)
 
 2021       WRITE(6,*)
' I,IPSQ(I),IPSAQ(I) ',
 
 2022      *    i,ipsq(i),ipsaq(i)
 
 2025       WRITE(6,*)
' I,ITSQ(I),ITSAQ(I) ',
 
 2026      *    i,itsq(i),itsaq(i)
 
 2033       IF(ipev.GE.4)
WRITE(6,*)
' collect v-v pairs NVV',nvv
 
 2042         IF(zuovp(iippv).AND.zuovt(iittv)) 
THEN 
 2044       IF(ipev.GE.6)
WRITE(6,
'(A,5I5)')
 
 2045      * 
' XKSAMP v-v loop IIPP,IITT,IIPPV,IITTV,NVV',iipp,iitt,iippv,iittv,nvv
 
 2046           zuovp(iippv)=.false.
 
 2047           zuovt(iittv)=.false.
 
 2049       IF(ipev.GE.4)
WRITE(6,*)
' collect v-v pairs NVV',nvv
 
 2086             IF(zuosp(j).AND.(ifrosp(j).EQ.iipp).AND.zuovt(iittv)) 
THEN 
 2088       IF(ipev.GE.6)
WRITE(6,
'(A,6I5)')
 
 2089      *
' XKSAMP s-v loop I(NN),J(IXPS),iitt,iittv,NSV,NDV',
 
 2090      +    i,j, iitt,iittv,nsv,ndv
 
 2091               zuovt(iittv)=.false.
 
 2094               IF(
rndm(v).GT.amedd.AND.idiqua.EQ.1)
THEN 
 2096                 CALL 
diqsv(ecm,iittv,j,irej)
 
 2097                 IF(irej.EQ.0)go to 260
 
 2106               amsvq1=xpsq(j)*xtvd(iittv)*ecm**2
 
 2107               amsvq2=xpsaq(j)*xtvq(iittv)*ecm**2
 
 2109               IF(ipsq(j).EQ.3)
THEN 
 2110                 IF(amsvq1.GT.amas)
THEN 
 2111                   xpsqxx=(xtvd(iittv)*ecm**2)
 
 2112           IF(xpsqxx.LE.1.
d-1)xpsqxx=1.
d-1
 
 2114                   xpsqxx=
sampex(xpsqth,xpsq(j))
 
 2115                   dxpsq=xpsq(j)-xpsqxx
 
 2116                   xpsq(j)=xpsq(j)-dxpsq
 
 2117                   xpvd(jxpv)=xpvd(jxpv)+dxpsq
 
 2118                 ELSEIF(amsvq1.LT.amas)
THEN 
 2119           IF(xtvd(iittv)*ecm**2.LE.1.
d-12)
THEN 
 2120             WRITE(6,*)
' xksamp: XTVD(IITTV)=0 ',iittv
 
 2123                   xpsqw=amas/(xtvd(iittv)*ecm**2)
 
 2126                   IF(xpvd(jxpv).GE.xdthr+dxpsq)
THEN  
 2127                     xpvd(jxpv)=xpvd(jxpv)-dxpsq
 
 2131                 IF(amsvq2.GT.amis)
THEN 
 2132                 ELSEIF(amsvq2.LT.amis)
THEN 
 2133           IF(xtvq(iittv)*ecm**2.LE.1.
d-12)
THEN 
 2134             WRITE(6,*)
' xksamp: XTVQ(IITTV)=0 ',iittv
 
 2137                   xpsqw=amis/(xtvq(iittv)*ecm**2)
 
 2138                   dxpsq=xpsqw-xpsaq(j)
 
 2140                   IF(xpvd(jxpv).GE.xdthr+dxpsq)
THEN  
 2141                     xpvd(jxpv)=xpvd(jxpv)-dxpsq
 
 2146                 IF(amsvq1.GT.amau)
THEN 
 2147           IF(xtvd(iittv)*ecm**2.LE.1.
d-12)
THEN 
 2148             WRITE(6,*)
' xksamp: XTVD(IITTV)=0 ',iittv
 
 2151                   xpsqth=amau/(xtvd(iittv)*ecm**2)
 
 2152                   xpsqxx=
sampex(xpsqth,xpsq(j))
 
 2153                   dxpsq=xpsq(j)-xpsqxx
 
 2154                   xpsq(j)=xpsq(j)-dxpsq
 
 2155                   xpvd(jxpv)=xpvd(jxpv)+dxpsq
 
 2156                 ELSEIF(amsvq1.LT.amau)
THEN 
 2157           IF(xtvd(iittv)*ecm**2.LE.1.
d-12)
THEN 
 2158             WRITE(6,*)
' xksamp: XTVD(IITTV)=0 ',iittv
 
 2161                   xpsqw=amau/(xtvd(iittv)*ecm**2)
 
 2164                   IF(xpvd(jxpv).GE.xdthr+dxpsq)
THEN  
 2165                     xpvd(jxpv)=xpvd(jxpv)-dxpsq
 
 2169                 IF(amsvq2.GT.amiu)
THEN 
 2170                 ELSEIF(amsvq2.LT.amiu)
THEN 
 2171           IF(xtvq(iittv)*ecm**2.LE.1.
d-12)
THEN 
 2172             WRITE(6,*)
' xksamp: XTVQ(IITTV)=0 ',iittv
 
 2175                   xpsqw=amiu/(xtvq(iittv)*ecm**2)
 
 2176                   dxpsq=xpsqw-xpsaq(j)
 
 2178                   IF(xpvd(jxpv).GE.xdthr+dxpsq)
THEN  
 2179                     xpvd(jxpv)=xpvd(jxpv)-dxpsq
 
 2188       IF(ipev.GE.6)
WRITE(6,
'(A,6I5)')
 
 2189      *
' XKSAMP s-v loop I(NN),J(IXPS),iitt,iittv,NSV,NDV',
 
 2190      +    i,j, iitt,iittv,nsv,ndv
 
 2205             IF(zuovp(iippv).AND.zuost(j).AND.(ifrost(j).EQ.iitt)) 
THEN 
 2207               IF(ipev.GE.6)
WRITE(6,*)
 
 2208      *        
' XKSAMP v-s loop IIPP,IITT,IIPPV,NVS,NVD,I,J,IXTS',
 
 2209      *        iipp,iitt,iippv,nvs,nvd,i,j,ixts
 
 2210               zuovp(iippv)=.false.
 
 2213               IF(
rndm(v).GT.amedd.AND.idiqua.EQ.1)
THEN 
 2215                 CALL 
diqvs(ecm,iippv,j,irej)
 
 2216                 IF(ipev.GE.6)
WRITE(6,*)
 
 2217      *          
' XKSAMP v-s loop IIPP,IITT,IIPPV,NVS,NVD,I,J,IXTS,JXTV' 
 2218      *          ,iipp,iitt,iippv,nvs,nvd,i,j,ixts,jxtv
 
 2219                 IF(irej.EQ.0)go to 290
 
 2227               amvsq1=xpvq(iippv)*xtsaq(j)*ecm**2
 
 2228               amvsq2=xpvd(iippv)*xtsq(j)*ecm**2
 
 2230               IF(itsq(j).EQ.3)
THEN 
 2232                 IF(amvsq1.LT.amis)
THEN 
 2233           IF(xpvq(iippv)*ecm**2.LE.1.
d-12)
THEN 
 2234             WRITE(6,*)
' xksamp: XPVQ(IIPPV)=0 ',iippv
 
 2237                   xtsqw=amis/(xpvq(iippv)*ecm**2)
 
 2238                   dxtsq=xtsqw-xtsaq(j)
 
 2240                   IF(xtvd(jxtv).GE.xdthr+dxtsq)
THEN  
 2241                     xtvd(jxtv)=xtvd(jxtv)-dxtsq
 
 2245                 IF(amvsq2.GT.amas)
THEN 
 2246           IF(xpvd(iippv)*ecm**2.LE.1.
d-12)
THEN 
 2247             WRITE(6,*)
' xksamp: XPVD(IIPPV)=0 ',iippv
 
 2250                   xtsqth=amas/(xpvd(iippv)*ecm**2)
 
 2251                   xtsqxx=
sampex(xtsqth,xtsq(j))
 
 2252                   dxtsq=xtsq(j)-xtsqxx
 
 2253                   xtsq(j)=xtsq(j)-dxtsq
 
 2254                   xtvd(jxtv)=xtvd(jxtv)+dxtsq
 
 2255                 ELSEIF(amvsq2.LT.amas)
THEN 
 2256           IF(xpvd(iippv)*ecm**2.LE.1.
d-12)
THEN 
 2257             WRITE(6,*)
' xksamp: XPVD(IIPPV)=0 ',iippv
 
 2260                   xtsqw=amas/(xpvd(iippv)*ecm**2)
 
 2263                   IF(xtvd(jxtv).GE.xdthr+dxtsq)
THEN  
 2264                     xtvd(jxtv)=xtvd(jxtv)-dxtsq
 
 2270                 IF(amvsq1.LT.amiu)
THEN 
 2271           IF(xpvq(iippv)*ecm**2.LE.1.
d-12)
THEN 
 2272             WRITE(6,*)
' xksamp: XPVQ(IIPPV)=0 ',iippv
 
 2275                   xtsqw=amiu/(xpvq(iippv)*ecm**2)
 
 2276                   dxtsq=xtsqw-xtsaq(j)
 
 2278                   IF(xtvd(jxtv).GE.xdthr+dxtsq)
THEN  
 2279                     xtvd(jxtv)=xtvd(jxtv)-dxtsq
 
 2283                 IF(amvsq2.GT.amau)
THEN 
 2284           IF(xpvd(iippv)*ecm**2.LE.1.
d-12)
THEN 
 2285             WRITE(6,*)
' xksamp: XPVD(IIPPV)=0 ',iippv
 
 2288                   xtsqth=amau/(xpvd(iippv)*ecm**2)
 
 2289                   xtsqxx=
sampex(xtsqth,xtsq(j))
 
 2290                   dxtsq=xtsq(j)-xtsqxx
 
 2291                   xtsq(j)=xtsq(j)-dxtsq
 
 2292                   xtvd(jxtv)=xtvd(jxtv)+dxtsq
 
 2293                 ELSEIF(amvsq2.LT.amau)
THEN 
 2294           IF(xpvd(iippv)*ecm**2.LE.1.
d-12)
THEN 
 2295             WRITE(6,*)
' xksamp: XPVD(IIPPV)=0 ',iippv
 
 2298                   xtsqw=amau/(xpvd(iippv)*ecm**2)
 
 2301                   IF(xtvd(jxtv).GE.xdthr+dxtsq)
THEN  
 2302                     xtvd(jxtv)=xtvd(jxtv)-dxtsq
 
 2328             IF(zuost(j).AND.(ifrost(j).EQ.iitt)) 
THEN 
 2330                 IF(zuosp(jj).AND.(ifrosp(jj).EQ.iipp)) 
THEN 
 2332                   IF(ipev.GE.6)
WRITE(6,
'(A,5I5)')
 
 2333      *            
' XKSAMP s-s loop IIPP,IITT,NSS',iipp,iitt,nss
 
 2336                   IF(ipev.GE.6)
WRITE(6,*)
 
 2337      *            
' XKSAMP s-s loop ,NCHSS1(NSS),NCHSS2(NSS),NSS ',
 
 2338      *            nchss1(nss),nchss2(nss),nss
 
 2345                   ssma1q=xpsq(jj)*xtsaq(j)*ecm**2
 
 2346                   ssma2q=xpsaq(jj)*xtsq(j)*ecm**2
 
 2347                   IF(ssma1q.LT.1.2d0.OR.ssma2q.LT.1.2d0) 
THEN 
 2359                   allket=(nvv+ixps+ixts)
 
 2360           IF(allket.LE.1.
d-5)
THEN 
 2361             WRITE(6,*)
' xksamp ALLKET=0' , allket
 
 2366           anvvo=min(ixpv,ixtv)
 
 2369           ansso=(ixpv+ixps)-anvvo-ansvo-anvso
 
 2370           IF(anvvo+ansso.LE.1.
d-5)
THEN 
 2371             WRITE(6,*)
' xksamp (...)=0' ,anvvo,ansso
 
 2375           IF(anvvo+ansso.GT.1.
d-5)valfra=anvvo/(anvvo+ansso)
 
 2382                       IF (nchvv1(ivv).NE.99.AND.nchvv2(ivv).NE.99)
THEN 
 2384                         inucpr=ifrovp(ixvpr)
 
 2386                         inucta=ifrovt(ixvta)
 
 2387                         IF(iipp.EQ.inucpr.OR.iitt.EQ.inucta)
THEN 
 2406                           IF(ipev.GE.6)
WRITE(6,*)
 
 2407      *                    
' XKSAMP before DIQSV ,NCHSS1(NSS),',
 
 2408      *                    
'NCHSS2(NSS),NSS ',
 
 2409      *                    nchss1(nss),nchss2(nss),nss
 
 2414                           IF(
rndm(v).GT.amedd.AND.idiqua.EQ.1)
THEN 
 2416                             CALL 
diqsv(ecm,ixvta,jj,irej)
 
 2417                             IF(irej.EQ.0)go to 4202
 
 2419                           IF(ipev.GE.6)
WRITE(6,*)
 
 2420      *            
' XKSAMP: NSS,NSV,NVS ',nss,nsv,nvs
 
 2424                           IF(ipev.GE.6)
WRITE(6,*)
 
 2425      *            
' XKSAMP: NSS,NSV,NVS ',nss,nsv,nvs
 
 2427                           amsvq1=xpsq(jj)*xtvd(ixvta)*ecm**2
 
 2428                           amsvq2=xpsaq(jj)*xtvq(ixvta)*ecm**2
 
 2430                           IF(ipev.GE.6)
WRITE(6,
'(A,5I5)')
 
 2431      *
' XKSAMP s-s loop rec sv,vs IXVTA,JXPV,JJ',ixvta,jxpv,jj
 
 2432                           IF(ipsq(jj).EQ.3)
THEN 
 2433                             IF(amsvq1.GT.amas)
THEN 
 2434                       IF(xtvd(ixvta)*ecm**2.LE.1.
d-12)
THEN 
 2436      *              
' xksamp: XTVD(IXVTA)=0 ',ixvta
 
 2439                               xpsqth=amas/(xtvd(ixvta)*ecm**2)
 
 2440                               xpsqxx=
sampex(xpsqth,xpsq(jj))
 
 2441                               dxpsq=xpsq(jj)-xpsqxx
 
 2442                               xpsq(jj)=xpsq(jj)-dxpsq
 
 2443                               xpvd(jxpv)=xpvd(jxpv)+dxpsq
 
 2444                             ELSEIF(amsvq1.LT.amas)
THEN 
 2445                       IF(xtvd(ixvta)*ecm**2.LE.1.
d-12)
THEN 
 2447      *                  
'xksamp: XTVD(IXVTA)=0 ',ixvta
 
 2450                               xpsqw=amas/(xtvd(ixvta)*ecm**2)
 
 2451                               dxpsq=xpsqw-xpsq(jj)
 
 2453                               IF(xpvd(jxpv).GE.xdthr+dxpsq)
THEN  
 2454                                 xpvd(jxpv)=xpvd(jxpv)-dxpsq
 
 2459                             IF(amsvq2.LT.amis)
THEN 
 2460                       IF(xtvq(ixvta)*ecm**2.LE.1.
d-12)
THEN 
 2462      *              
' xksamp: XTVQ(IXVTA)=0 ',ixvta
 
 2465                               xpsqw=amis/(xtvq(ixvta)*ecm**2)
 
 2466                               dxpsq=xpsqw-xpsaq(jj)
 
 2468                               IF(xpvd(jxpv).GE.xdthr+dxpsq)
THEN  
 2469                                 xpvd(jxpv)=xpvd(jxpv)-dxpsq
 
 2475                             IF(amsvq1.GT.amau)
THEN 
 2476                       IF(xtvd(ixvta)*ecm**2.LE.1.
d-12)
THEN 
 2478      *              
' xksamp: XTVD(IXVTA)=0 ',ixvta
 
 2481                               xpsqth=amau/(xtvd(ixvta)*ecm**2)
 
 2482                               xpsqxx=
sampex(xpsqth,xpsq(jj))
 
 2483                               dxpsq=xpsq(jj)-xpsqxx
 
 2484                               xpsq(jj)=xpsq(jj)-dxpsq
 
 2485                               xpvd(jxpv)=xpvd(jxpv)+dxpsq
 
 2486                             ELSEIF(amsvq1.LT.amau)
THEN 
 2487                       IF(xtvd(ixvta)*ecm**2.LE.1.
d-12)
THEN 
 2489      *              
' xksamp: XTVD(IXVTA)=0 ',ixvta
 
 2492                               xpsqw=amau/(xtvd(ixvta)*ecm**2)
 
 2493                               dxpsq=xpsqw-xpsq(jj)
 
 2495                               IF(xpvd(jxpv).GE.xdthr+dxpsq)
THEN  
 2496                                 xpvd(jxpv)=xpvd(jxpv)-dxpsq
 
 2501                             IF(amsvq2.LT.amiu)
THEN 
 2502                       IF(xtvq(ixvta)*ecm**2.LE.1.
d-12)
THEN 
 2504      *              
' xksamp: XTVQ(IXVTA)=0 ',ixvta
 
 2507                               xpsqw=amiu/(xtvq(ixvta)*ecm**2)
 
 2508                               dxpsq=xpsqw-xpsaq(jj)
 
 2510                               IF(xpvd(jxpv).GE.xdthr+dxpsq)
THEN  
 2511                                 xpvd(jxpv)=xpvd(jxpv)-dxpsq
 
 2523                           IF(
rndm(v).GT.amedd.AND.idiqua.EQ.1)
THEN 
 2525                             CALL 
diqvs(ecm,ixvpr,j,irej)
 
 2526                             IF(irej.EQ.0)go to 4203
 
 2528                           IF(ipev.GE.6)
WRITE(6,*)
 
 2529      *            
' XKSAMP: NSS,NSV,NVS ',nss,nsv,nvs
 
 2533                           IF(ipev.GE.6)
WRITE(6,*)
 
 2534      *            
' XKSAMP: NSS,NSV,NVS ',nss,nsv,nvs
 
 2536                           amvsq1=xpvq(ixvpr)*xtsaq(j)*ecm**2
 
 2537                           amvsq2=xpvd(ixvpr)*xtsq(j)*ecm**2
 
 2539                           IF(ipev.GE.6)
WRITE(6,
'(A,5I5)')
 
 2540      *
' XKSAMP s-s loop rec vs IXVTA,JXPV,JJ',ixvta,jxpv,jj
 
 2541                           IF(itsq(j).EQ.3)
THEN 
 2543                             IF(amvsq1.LT.amis)
THEN 
 2544                       IF(xpvq(ixvpr)*ecm**2.LE.1.
d-12)
THEN 
 2546      *              
' xksamp: XPVQ(IXVPR)=0 ',ixvpr
 
 2549                               xtsqw=amis/(xpvq(ixvpr)*ecm**2)
 
 2550                               dxtsq=xtsqw-xtsaq(j)
 
 2552                               IF(xtvd(jxtv).GE.xdthr+dxtsq)
THEN  
 2553                                 xtvd(jxtv)=xtvd(jxtv)-dxtsq
 
 2557                             IF(amvsq2.GT.amas)
THEN 
 2558                       IF(xpvd(ixvpr)*ecm**2.LE.1.
d-12)
THEN 
 2560      *              
' xksamp: XPVD(IXVPR)=0 ',ixvpr
 
 2563                               xtsqth=amas/(xpvd(ixvpr)*ecm**2)
 
 2564                               xtsqxx=
sampex(xtsqth,xtsq(j))
 
 2565                               dxtsq=xtsq(j)-xtsqxx
 
 2566                               xtsq(j)=xtsq(j)-dxtsq
 
 2567                               xtvd(jxtv)=xtvd(jxtv)+dxtsq
 
 2568                             ELSEIF(amvsq2.LT.amas)
THEN 
 2569                       IF(xpvd(ixvpr)*ecm**2.LE.1.
d-12)
THEN 
 2571      *              
' xksamp: XPVD(IXVPR)=0 ',ixvpr
 
 2574                               xtsqw=amas/(xpvd(ixvpr)*ecm**2)
 
 2577                               IF(xtvd(jxtv).GE.xdthr+dxtsq)
THEN  
 2578                                 xtvd(jxtv)=xtvd(jxtv)-dxtsq
 
 2584                             IF(amvsq1.LT.amiu)
THEN 
 2585                       IF(xpvq(ixvpr)*ecm**2.LE.1.
d-12)
THEN 
 2587      *              
' xksamp: XPVQ(IXVPR)=0 ',ixvpr
 
 2590                               xtsqw=amiu/(xpvq(ixvpr)*ecm**2)
 
 2591                               dxtsq=xtsqw-xtsaq(j)
 
 2593                               IF(xtvd(jxtv).GE.xdthr+dxtsq)
THEN  
 2594                                 xtvd(jxtv)=xtvd(jxtv)-dxtsq
 
 2598                             IF(amvsq2.GT.amau)
THEN 
 2599                       IF(xpvd(ixvpr)*ecm**2.LE.1.
d-12)
THEN 
 2601      *              
' xksamp: XPVD(IXVPR)=0 ',ixvpr
 
 2604                               xtsqth=amau/(xpvd(ixvpr)*ecm**2)
 
 2605                               xtsqxx=
sampex(xtsqth,xtsq(j))
 
 2606                               dxtsq=xtsq(j)-xtsqxx
 
 2607                               xtsq(j)=xtsq(j)-dxtsq
 
 2608                               xtvd(jxtv)=xtvd(jxtv)+dxtsq
 
 2609                             ELSEIF(amvsq2.LT.amau)
THEN 
 2610                       IF(xpvd(ixvpr)*ecm**2.LE.1.
d-12)
THEN 
 2612      *              
' xksamp: XPVD(IXVPR)=0 ',ixvpr
 
 2615                               xtsqw=amau/(xpvd(ixvpr)*ecm**2)
 
 2618                               IF(xtvd(jxtv).GE.xdthr+dxtsq)
THEN  
 2619                                 xtvd(jxtv)=xtvd(jxtv)-dxtsq
 
 2640                   IF(
rndm(v).GT.2.d0*amedd-1.d0.AND.idiqua.EQ.1)
THEN 
 2642                     CALL 
diqdss(ecm,j,jj,irej)
 
 2646                       IF(ipev.GE.6)
WRITE(6,*)
 
 2647      *                
' XKSAMP AFTER DIQDSS IREJ=0',
 
 2648      *                
',NCHSS1(NSS),NCHSS2(NSS),NSS ',
 
 2649      *                nchss1(nss),nchss2(nss),nss
 
 2654                   IF(
rndm(v).GT.2.d0*amedd-1.d0.AND.idiqua.EQ.1)
THEN 
 2656                     CALL 
diqssd(ecm,j,jj,irej)
 
 2660                       IF(ipev.GE.6)
WRITE(6,*)
 
 2661      *                
' XKSAMP AFTER DIQSSD IREJ=0',
 
 2662      *                
',NCHSS1(NSS),NCHSS2(NSS),NSS ',
 
 2663      *                nchss1(nss),nchss2(nss),nss
 
 2667                   ssma1q=xpsq(jj)*xtsaq(j)*ecm**2
 
 2668                   ssma2q=xpsaq(jj)*xtsq(j)*ecm**2
 
 2669                   IF(ssma1q.LT.ssmimq.OR.ssma2q.LT.ssmimq) 
THEN 
 2672                     IF((xtvd(jxtv).GT.xdthr+3.5d0*xssthr)
 
 2674      +              .GT.xdthr+3.5d0*xssthr)) 
THEN 
 2676                       xspmax=1.0 - xpvq(jxpv) - xdthr - 1.2*xssthr
 
 2677                       xstmax=1.0 - xtvq(jxtv) - xdthr - 1.2*xssthr
 
 2679                       IF((xspmax.LE.xssthr+0.05d0) .OR.(xstmax.LE.xssthr
 
 2680      +                +0.05d0))                                 goto 380
 
 2685                       IF(xssthr.GT.0.05d0) 
THEN 
 2686                         xpsqi=
betrej(xseacu,unosea,xssthr,xspmax)
 
 2687                         xpsaqi=
betrej(xseacu,unosea,xssthr,xspmax)
 
 2690                         xpsqi=
dbetar(xseacu,unosea)
 
 2691                         IF(xpsqi.LT.xssthr.OR.xpsqi.GT.xspmax)  goto 320
 
 2693                         xpsaqi=
dbetar(xseacu,unosea)
 
 2694                         IF(xpsaqi.LT.xssthr.OR.xpsaqi.GT.xspmax)
 
 2698                       xpvdco=xpvd(jxpv) - xpsqi - xpsaqi + xpsq(jj) +
 
 2700                       IF(xpvdco.GT.xdthr) 
THEN 
 2703                       ELSEIF(icous.LT.5) 
THEN 
 2714                       IF(xssthr.GT.0.05d0)
THEN 
 2715                         xtsqi=
betrej(xseacu,unosea,xssthr,xstmax)
 
 2716                         xtsaqi=
betrej(xseacu,unosea,xssthr,xstmax)
 
 2719                         xtsqi=
dbetar(xseacu,unosea)
 
 2720                         IF(xtsqi.LT.xssthr.OR.xtsqi.GT.xstmax)  goto 360
 
 2722                         xtsaqi=
dbetar(xseacu,unosea)
 
 2723                         IF(xtsaqi.LT.xssthr.OR.xtsaqi.GT.xstmax)
 
 2727                       xtvdco=xtvd(jxtv) - xtsqi - xtsaqi + xtsq(j) +
 
 2729                       IF(xtvdco.LT.xdthr) 
THEN 
 2731                         IF(icous.LT.5)                          goto 350
 
 2769           xpvq(iitop)=xpvq(iitop) + xpsq(i) + xpsaq(i)
 
 2777           xtvq(iitot)=xtvq(iitot) + xtsq(i) + xtsaq(i)
 
 2796         WRITE(6,
'(A)') 
' XKSAMP: I,INTVV1,INTVV2,IFROVP,IFROVT' 
 2800           WRITE(6,
'(5I5)') i,inup,inut,ifrovp(inup),ifrovt(inut)
 
 2802        WRITE(6,
'(A)')
'XKSAMP:I(NSV),INTSV1,INTSV2,IFROSP,IFROVT' 
 2806           WRITE(6,
'(5I5)') i,inup,inut,ifrosp(inup),ifrovt(inut)
 
 2808         WRITE(6,
'(A)') 
' XKSAMP: I,INTVS1,INTVS2,IFROVP,IFROST' 
 2812           WRITE(6,
'(5I5)') i,inup,inut,ifrovp(inup),ifrost(inut)
 
 2814         WRITE(6,
'(A)') 
' XKSAMP: I,INTSS1,INTSS2,IFROSP,IFROST' 
 2818           WRITE(6,
'(5I5)') i,inup,inut,ifrosp(inup),ifrost(inut)
 
 2822      +  
' XKSAMP :  FINAL X-VALUES AFTER POTENTIAL CORRECTION' 
 2825           WRITE(6,1020) i,xpvq(i),xpvd(i),ifrovp(i),itovp(i),zuovp(i)
 
 2826           WRITE(6,*)
' I(1-IXPV),IPVQ(I),IPPV1(I),IPPV2(I) ',
 
 2827      *    i,ipvq(i),ippv1(i),ippv2(i)
 
 2831           WRITE(6,1040) i,xpsq(i),xpsaq(i),ifrosp(i),zuosp(i)
 
 2832           WRITE(6,*)
' I(1-IXPS),IPSQ(I),IPSAQ(I) ',
 
 2833      *    i,ipsq(i),ipsaq(i)
 
 2837           WRITE(6,1020) i,xtvq(i),xtvd(i),ifrovt(i),itovt(i),zuovt(i)
 
 2838           WRITE(6,*)
' I(1-IXTV),ITVQ(I),ITTV1(I),ITTV2(I) ',
 
 2839      *    i,itvq(i),ittv1(i),ittv2(i)
 
 2843           WRITE(6,1040) i,xtsq(i),xtsaq(i),ifrost(i),zuost(i)
 
 2844           WRITE(6,*)
' I(1-IXTS),ITSQ(I),ITSAQ(I) ',
 
 2845      *    i,itsq(i),itsaq(i)
 
 2848       IF(ipev.GE.6)
WRITE(6,
'(A,6I5)')
 
 2849      *
' XKSAMP NSV,NDV,NVS,NVD',
 
 2860       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 2875       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 2877      *                ,xpsu(248),xtsu(248)
 
 2878      *                ,xpsut(248),xtsut(248)
 
 2880       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 2881      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 2882      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 2884      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 2898       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 2904      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 2906       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 2907       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 2914       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 2917      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 2922       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 3014       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 3020           WRITE (6,
'(A,2I5)') .EQ.
' XKSAMP: NHKKNMXHKK ',nhkk,
nmxhkk 
 3032         phkk(3,nhkk)=xpvq(i)
 
 3033         phkk(4,nhkk)=xpvq(i)
 
 3036         vhkk(1,nhkk)=vhkk(1,kkk)
 
 3037         vhkk(2,nhkk)=vhkk(2,kkk)
 
 3038         vhkk(3,nhkk)=vhkk(3,kkk)
 
 3041         IF (iphkk.GE.3) 
WRITE(6,1000) nhkk,isthkk(nhkk),idhkk(nhkk),
 
 3042      +  jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 3043      +  (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 3048           WRITE (6,
'(A,2I5)') .EQ.
' XKSAMP: NHKKNMXHKK ',nhkk,
nmxhkk 
 3060         phkk(3,nhkk)=xpvd(i)
 
 3061         phkk(4,nhkk)=xpvd(i)
 
 3064         vhkk(1,nhkk)=vhkk(1,kkk)
 
 3065         vhkk(2,nhkk)=vhkk(2,kkk)
 
 3066         vhkk(3,nhkk)=vhkk(3,kkk)
 
 3070         IF (iphkk.GE.7) 
WRITE(6,1000) nhkk,isthkk(nhkk),idhkk(nhkk),
 
 3071      +  jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 3072      +  (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 3074  1000 
FORMAT (i6,i4,5i6,9e10.2)
 
 3080           WRITE (6,
'(A,2I5)') .EQ.
' XKSAMP: NHKKNMXHKK ',nhkk,
nmxhkk 
 3092         phkk(3,nhkk)=xpsq(i)
 
 3093         phkk(4,nhkk)=xpsq(i)
 
 3096         vhkk(1,nhkk)=vhkk(1,kkk)
 
 3097         vhkk(2,nhkk)=vhkk(2,kkk)
 
 3098         vhkk(3,nhkk)=vhkk(3,kkk)
 
 3101         IF (iphkk.GE.7) 
WRITE(6,1000) nhkk,isthkk(nhkk),idhkk(nhkk),
 
 3102      +  jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 3103      +  (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 3108           WRITE (6,
'(A,2I5)') .EQ.
' XKSAMP: NHKKNMXHKK ',nhkk,
nmxhkk 
 3120         phkk(3,nhkk)=xpsaq(i)
 
 3121         phkk(4,nhkk)=xpsaq(i)
 
 3124         vhkk(1,nhkk)=vhkk(1,kkk)
 
 3125         vhkk(2,nhkk)=vhkk(2,kkk)
 
 3126         vhkk(3,nhkk)=vhkk(3,kkk)
 
 3130         IF (iphkk.GE.7) 
WRITE(6,1000) nhkk,isthkk(nhkk),idhkk(nhkk),
 
 3131      +  jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 3132      +  (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 3139           WRITE (6,
'(A,2I5)') .EQ.
' XKSAMP: NHKKNMXHKK ',nhkk,
nmxhkk 
 3151         phkk(3,nhkk)=xtvq(i)
 
 3152         phkk(4,nhkk)=xtvq(i)
 
 3155         vhkk(1,nhkk)=vhkk(1,kkk)
 
 3156         vhkk(2,nhkk)=vhkk(2,kkk)
 
 3157         vhkk(3,nhkk)=vhkk(3,kkk)
 
 3160         IF (iphkk.GE.7) 
WRITE(6,1000) nhkk,isthkk(nhkk),idhkk(nhkk),
 
 3161      +  jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 3162      +  (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 3167           WRITE (6,
'(A,2I5)') .EQ.
' XKSAMP: NHKKNMXHKK ',nhkk,
nmxhkk 
 3179         phkk(3,nhkk)=xtvd(i)
 
 3180         phkk(4,nhkk)=xtvd(i)
 
 3183         vhkk(1,nhkk)=vhkk(1,kkk)
 
 3184         vhkk(2,nhkk)=vhkk(2,kkk)
 
 3185         vhkk(3,nhkk)=vhkk(3,kkk)
 
 3189         IF (iphkk.GE.7) 
WRITE(6,1000) nhkk,isthkk(nhkk),idhkk(nhkk),
 
 3190      +  jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 3191      +  (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 3198           WRITE (6,
'(A,2I5)') .EQ.
' XKSAMP: NHKKNMXHKK ',nhkk,
nmxhkk 
 3210         phkk(3,nhkk)=xtsq(i)
 
 3211         phkk(4,nhkk)=xtsq(i)
 
 3214         vhkk(1,nhkk)=vhkk(1,kkk)
 
 3215         vhkk(2,nhkk)=vhkk(2,kkk)
 
 3216         vhkk(3,nhkk)=vhkk(3,kkk)
 
 3219         IF (iphkk.GE.7) 
WRITE(6,1000) nhkk,isthkk(nhkk),idhkk(nhkk),
 
 3220      +  jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 3221      +  (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 3226           WRITE (6,
'(A,2I5)') .EQ.
' XKSAMP: NHKKNMXHKK ',nhkk,
nmxhkk 
 3238         phkk(3,nhkk)=xtsaq(i)
 
 3239         phkk(4,nhkk)=xtsaq(i)
 
 3242         vhkk(1,nhkk)=vhkk(1,kkk)
 
 3243         vhkk(2,nhkk)=vhkk(2,kkk)
 
 3244         vhkk(3,nhkk)=vhkk(3,kkk)
 
 3248         IF (iphkk.GE.7) 
WRITE(6,1000) nhkk,isthkk(nhkk),idhkk(nhkk),
 
 3249      +  jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 3250      +  (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 3260       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 3273       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 3275      *                ,xpsu(248),xtsu(248)
 
 3276      *                ,xpsut(248),xtsut(248)
 
 3278       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 3279      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 3280      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 3282      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 3296       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 3302      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 3304       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 3305       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 3312       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 3315      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 3320       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 3412       COMMON /nncms/  gamcm,bgcm,umo,pcm,eproj,pproj
 
 3414       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
 3416       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
 3417      +ipadis,ishmal,lpauli
 
 3419       COMMON /cmhico/ cmhis
 
 3421       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 3433       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 3434      +iibar(210),k1(210),k2(210)
 
 3438        COMMON /bufues/ bnnvv,bnnss,bnnsv,bnnvs,bnncc,
 
 3439      *                 bnndv,bnnvd,bnnds,bnnsd,
 
 3441      *                 bptvv,bptss,bptsv,bptvs,bptcc,bptdv,
 
 3442      *                 bptvd,bptds,bptsd,
 
 3444      *                 beevv,beess,beesv,beevs,beecc,beedv,
 
 3445      *                 beevd,beeds,beesd,
 
 3447      *                ,bnndi,bptdi,beedi
 
 3448      *                ,bnnzd,bnndz,bptzd,bptdz,beezd,beedz
 
 3449        COMMON /ncoucs/ bcouvv,bcouss,bcousv,bcouvs,
 
 3450      *                 bcouzz,bcouhh,bcouds,bcousd,
 
 3451      *                 bcoudz,bcouzd,bcoudi,
 
 3452      *                 bcoudv,bcouvd,bcoucc
 
 3453        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 3454      *                 anndv,annvd,annds,annsd,
 
 3456      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 3458      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 3461      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 3462        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 3463      *                 acouzz,acouhh,acouds,acousd,
 
 3464      *                 acoudz,acouzd,acoudi,
 
 3465      *                 acoudv,acouvd,acoucc
 
 3467       COMMON /jspa/pxs(40000),pys(40000),pzs(40000),hes(40000),nnnps
 
 3471       COMMON /seadiq/ lseadi
 
 3472       COMMON /minij/iminij,nomje,nomjer,nrejev,nomjt,nomjtr
 
 3473       common/diquax/amedd,idiqua,idiquu
 
 3474       COMMON /secint/isecin
 
 3480       IF (ipco.GE.1) 
WRITE(6,1000) nvv,nvs,nsv,nss
 
 3481  1000 
FORMAT (
' ENTERING HADRKK NVV,NVS,NSV,NSS '/5i5)
 
 3546       IF(ihada.OR.ihadss) 
THEN 
 3550       IF(ihada.OR.ihadsv) 
THEN 
 3553       IF(ihada.OR.ihadvs) 
THEN 
 3556         IF (iminij.EQ.1) CALL 
hadrhh 
 3558       IF(idiquu.EQ.1)  CALL 
hadrdz 
 3559       IF(idiquu.EQ.1)  CALL 
hadrzd 
 3568       IF(idiqua.EQ.1)  CALL 
hadrds 
 3576       IF(idiqua.EQ.1)  CALL 
hadrsd 
 3583       IF(ihada.OR.ihadsv) 
THEN 
 3594       IF(idiqua.EQ.1)  CALL 
hadrdv 
 3600       IF(ihada.OR.ihadvs) 
THEN 
 3610       IF(idiqua.EQ.1)  CALL 
hadrvd 
 3616       IF(ihada.OR.ihadvv) 
THEN 
 3640         WRITE(6,
'(A)') 
' HADRONS FROM HADRKK / NUCLEON-NUCLEON CMS' 
 3641         DO 10 i=nhkkh1+1,nhkk
 
 3642     IF(isthkk(i).EQ.1)
THEN 
 3643           pxsu=pxsu + phkk(1,i)
 
 3644           pysu=pysu + phkk(2,i)
 
 3645           pzsu=pzsu + phkk(3,i)
 
 3646           esum=esum + phkk(4,i)
 
 3648           ichsu=ichsu + iich(nref)
 
 3649           ibasu=ibasu + iibar(nref)
 
 3651      *    
WRITE(6,1010)i,(phkk(j,i),j=1,5), iich(nref),iibar(nref),nref,
 
 3653  1010 
FORMAT(5
x,i4,5(1pe11.3),2i2,i5,a10)
 
 3656         WRITE(6,1020) pxsu,pysu,pzsu,esum,ichsu,ibasu
 
 3657  1020 
FORMAT(
' PXSU,PYSU,PZSU,ESUM,ICHSU,IBASU'/4f10.3,2i5)
 
 3666       DO 20 i=nhkkh1+1,nhkk
 
 3669         IF (cmhiss.EQ.0.d0)
THEN 
 3672           phkk(3,i) = gamcm*pznn + bgcm*enn
 
 3673           phkk(4,i) = gamcm*enn  + bgcm*pznn
 
 3675         ehecc=
sqrt(phkk(1,i)** 2+ phkk(2,i)** 2+ phkk(3,i)** 2+ phkk
 
 3677         IF (abs(ehecc-phkk(4,i)).GT.0.001d0) 
THEN 
 3686       IF(isecin.EQ.1)CALL 
sewew(1,nhkkh1)
 
 3700         WRITE(6,
'(A)') 
' HADRONS FROM HADRKK / CMS SYSTEM' 
 3701         DO 30 i=nhkkh1+1,nhkk
 
 3702       IF(isthkk(i).EQ.1)
THEN 
 3703           pxsu=pxsu + phkk(1,i)
 
 3704           pysu=pysu + phkk(2,i)
 
 3705           pzsu=pzsu + phkk(3,i)
 
 3706           esum=esum + phkk(4,i)
 
 3708           ichsu=ichsu + iich(nref)
 
 3709           ibasu=ibasu + iibar(nref)
 
 3711      *    
WRITE(6,1010) i, (phkk(j,i),j=1,5), iich(nref),iibar(nref),
 
 3715         WRITE(6,1020) pxsu,pysu,pzsu,esum,ichsu,ibasu
 
 3727       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 3744       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 3746      *                ,xpsu(248),xtsu(248)
 
 3747      *                ,xpsut(248),xtsut(248)
 
 3749       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 3750      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 3751      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 3753      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 3767       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 3773      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 3775       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 3776       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 3783       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 3786      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 3791       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 3879       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 3881       COMMON /abrvv/ amcvv1(248),amcvv2(248),gacvv1(248),gacvv2(248),
 
 3882      +bgxvv1(248),bgyvv1(248),bgzvv1(248), bgxvv2(248),bgyvv2(248),
 
 3883      +bgzvv2(248), nchvv1(248),nchvv2(248),ijcvv1(248),ijcvv2(248),
 
 3884      +pqvva1(248,4),pqvva2(248,4), pqvvb1(248,4),pqvvb2(248,4)
 
 3886       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 3889       parameter(nfimax=249)
 
 3890       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 3891      +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 3892       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 3895       COMMON /projk/ iprojk
 
 3898        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 3899      *                 anndv,annvd,annds,annsd,
 
 3901      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 3903      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 3906      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 3907        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 3908      *                 acouzz,acouhh,acouds,acousd,
 
 3909      *                 acoudz,acouzd,acoudi,
 
 3910      *                 acoudv,acouvd,acoucc
 
 3911        common/popcck/pdbck,pdbse,pdbseu,
 
 3912      *  ijpock,irejck,ick4,ihad4,ick6,ihad6
 
 3913      *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
 
 3914      *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
 
 3915      *isea43,isea63,irejao
 
 3917       COMMON /zsea/zseaav,zseasu,anzsea
 
 3919       dimension poj(4),pat(4)
 
 3921       IF(iphkk.GE.6)
WRITE (6,
'( A)') 
' hadrVV' 
 3926         IF(nchvv1(i).EQ.99.AND.nchvv2(i).EQ.99) go to 50
 
 3930         IF (ipco.GE.1) 
WRITE (6,1000) ipvq(is1),ippv1(is1),ippv2(is1),
 
 3931      +  itvq(is2),ittv1(is2),ittv2(is2), amcvv1(i),amcvv2(i),gacvv1(i),
 
 3932      +  gacvv2(i), bgxvv1(i),bgyvv1(i),bgzvv1(i), bgxvv2(i),bgyvv2(i),
 
 3933      +  bgzvv2(i), nchvv1(i),nchvv2(i),ijcvv1(i),ijcvv2(i), pqvva1(i,4),
 
 3934      +  pqvva2(i,4),pqvvb1(i,4),pqvvb2(i,4)
 
 3938  1000 
FORMAT(6i5,10f9.2/10
x,4i5,4f12.4)
 
 3943         IF(ibproj.GE.0) 
THEN 
 3959         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 3960         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 3961         CALL 
parpt(2,pt1,pt2,1,nevt)
 
 3972        ippp = ifrovp(intvv1(i))
 
 3973        ittt = ifrovt(intvv2(i))
 
 3984         WRITE(6,*)
' VV q-qq ,IFB1,IFB2,IFB3,',
 
 3985      *  
'INTVV1=IS1,INTVV2=IS2,JIPP,JITT',
 
 3986      *  ifb1,ifb2,ifb3,intvv1(i),intvv2(i),jipp,jitt
 
 3988     IF(nobam.EQ.3.OR.nchvv1(i).NE.0)
THEN 
 3990         CALL 
hadjet(nhad,amcvv1(i),poj,pat,gacvv1(i),bgxvv1(i), bgyvv1
 
 3991      +  (i),bgzvv1(i),ifb1,ifb2,ifb3,ifb4, ijcvv1(i),ijcvv1(i),nobam,
 
 3994     aack=float(ick4)/float(ick4+ihad4+1)
 
 3995     IF((nchvv1(i).EQ.0).AND.
 
 3997         zseawu=
rndm(bb)*2.d0*zseaav
 
 3998         rseack=float(jitt)*pdbse+ zseawu*pdbseu
 
 3999     IF(ipco.GE.1)
WRITE(6,*)
'HADJSE JITT,RSEACK,PDBSE 1 dpmnuc3',
 
 4002         IF(
rndm(v).LE.rseack)
THEN 
 4004     IF(amcvv1(i).GT.2.3d0)
THEN 
 4006           CALL 
hadjse(nhad,amcvv1(i),poj,pat,gacvv1(i),bgxvv1(i), bgyvv1
 
 4007      +    (i),bgzvv1(i),ifb1,ifb2,ifb3,ifb4, ijcvv1(i),ijcvv1(i),nobam,
 
 4008      +    nchvv1(i),7,irejss,iissqq)
 
 4009       IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JITT,',
 
 4010      *    
'RSEACK,IREJSS 1 dpmnuc3 ',
 
 4011      +    jitt,rseack,irejss
 
 4014     IF(irejss.EQ.1)irejse=irejse+1
 
 4015     IF(irejss.EQ.3)irejs3=irejs3+1
 
 4016     IF(irejss.EQ.2)irejs0=irejs0+1
 
 4017         CALL 
hadjet(nhad,amcvv1(i),poj,pat,gacvv1(i),bgxvv1(i), bgyvv1
 
 4018      +  (i),bgzvv1(i),ifb1,ifb2,ifb3,ifb4, ijcvv1(i),ijcvv1(i),nobam,
 
 4029         ELSEIF((ijpock.EQ.1).AND.
 
 4030      *      (aack.LE.pdbck))
THEN 
 4032         CALL 
hadjck(nhad,amcvv1(i),poj,pat,gacvv1(i),bgxvv1(i), bgyvv1
 
 4033      +  (i),bgzvv1(i),ifb1,ifb2,ifb3,ifb4, ijcvv1(i),ijcvv1(i),nobam,
 
 4037         CALL 
hadjet(nhad,amcvv1(i),poj,pat,gacvv1(i),bgxvv1(i), bgyvv1
 
 4038      +  (i),bgzvv1(i),ifb1,ifb2,ifb3,ifb4, ijcvv1(i),ijcvv1(i),nobam,
 
 4042     IF(irej.EQ.0)ick4=ick4+1
 
 4044         CALL 
hadjet(nhad,amcvv1(i),poj,pat,gacvv1(i),bgxvv1(i), bgyvv1
 
 4045      +  (i),bgzvv1(i),ifb1,ifb2,ifb3,ifb4, ijcvv1(i),ijcvv1(i),nobam,
 
 4060             WRITE (6,
'(A,2I5/A)') .EQ.
' HADRVV: NHKKNMXHKK ',nhkk,
nmxhkk 
 4064           ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 4065         IF (abs(ehecc-hef(j)).GT.0.001d0) 
THEN 
 4074           ptvv=ptvv+
sqrt(pxf(j)**2+pyf(j)**2)
 
 4077           IF(ibarf(j).EQ.500)istist=2
 
 4079      *                pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),1)
 
 4081           IF(idhkk(nhkk).EQ.99999) 
WRITE (6,1010) nhkk,nref(j),idhkk
 
 4083  1010 
FORMAT (
' NHKK,NREF(J),  ',3i10)
 
 4084           imohkk=jmohkk(1,nhkk)
 
 4085           IF(imohkk.LE.0.OR.imohkk.GT.
nmxhkk)
THEN 
 4086             WRITE(6,
'(A,I10)')
' HADRVV out of range IMOHKK= ',i10
 
 4090       WRITE(6,*)
' From HADRVV 1 first chain after HKKFIL' 
 4091           IF (iphkk.GE.0) 
WRITE(6,1020) nhkk, isthkk(nhkk),idhkk(nhkk),
 
 4092      +    jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 4093      +    (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 4095  1020 
FORMAT (i6,i4,5i6,9e10.2)
 
 4108         IF(ibproj.GT.0) 
THEN 
 4113         ELSEIF(ibproj.EQ.0) 
THEN 
 4132         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 4133         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 4134         CALL 
parpt(2,pt1,pt2,1,nevt)
 
 4145        ippp = ifrovp(intvv1(i))
 
 4146        ittt = ifrovt(intvv2(i))
 
 4157         WRITE(6,*)
' VV qq-q ,IFB1,IFB2,IFB3,',
 
 4158      *  
'INTVV1=IS1,INTVV2=IS2,JIPP,JITT',
 
 4159      *  ifb1,ifb2,ifb3,intvv1(i),intvv2(i),jipp,jitt
 
 4161     IF(nobam.EQ.5.OR.nobam.EQ.3.OR.nchvv2(i).NE.0)
THEN 
 4163         CALL 
hadjet(nhad,amcvv2(i),poj,pat,gacvv2(i),bgxvv2(i), bgyvv2
 
 4164      +  (i),bgzvv2(i),ifb1,ifb2,ifb3,ifb4, ijcvv2(i),ijcvv2(i),nobam,
 
 4167     aack=float(ick6)/float(ick6+ihad6+1)
 
 4168     IF((nchvv2(i).EQ.0).AND.
 
 4170         zseawu=
rndm(bb)*2.d0*zseaav
 
 4171         rseack=float(jipp)*pdbse+ zseawu*pdbseu
 
 4172     IF(ipco.GE.1)
WRITE(6,*)
'HADJSE JIPP,RSEACK,PDBSE 2 dpmnuc3',
 
 4175     IF(
rndm(v).LE.rseack)
THEN 
 4177     IF(amcvv2(i).GT.2.3d0)
THEN 
 4179           CALL 
hadjse(nhad,amcvv2(i),poj,pat,gacvv2(i),bgxvv2(i), bgyvv2
 
 4180      +    (i),bgzvv2(i),ifb1,ifb2,ifb3,ifb4, ijcvv2(i),ijcvv2(i),nobam,
 
 4181      +    nchvv2(i),8,irejss,iissqq)
 
 4182       IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JIPP,',
 
 4183      *    
'RSEACK,IREJSS 2 dpmnux3 ',
 
 4184      +    jipp,rseack,irejss
 
 4187     IF(irejss.EQ.1)irejse=irejse+1
 
 4188     IF(irejss.EQ.3)irejs3=irejs3+1
 
 4189     IF(irejss.EQ.2)irejs0=irejs0+1
 
 4190         CALL 
hadjet(nhad,amcvv2(i),poj,pat,gacvv2(i),bgxvv2(i), bgyvv2
 
 4191      +  (i),bgzvv2(i),ifb1,ifb2,ifb3,ifb4, ijcvv2(i),ijcvv2(i),nobam,
 
 4202         ELSEIF((ijpock.EQ.1).AND.
 
 4203      *      (aack.LE.pdbck))
THEN 
 4205         CALL 
hadjck(nhad,amcvv2(i),poj,pat,gacvv2(i),bgxvv2(i), bgyvv2
 
 4206      +  (i),bgzvv2(i),ifb1,ifb2,ifb3,ifb4, ijcvv2(i),ijcvv2(i),nobam,
 
 4210         CALL 
hadjet(nhad,amcvv2(i),poj,pat,gacvv2(i),bgxvv2(i), bgyvv2
 
 4211      +  (i),bgzvv2(i),ifb1,ifb2,ifb3,ifb4, ijcvv2(i),ijcvv2(i),nobam,
 
 4215     IF(irej.EQ.0)ick6=ick6+1
 
 4217         CALL 
hadjet(nhad,amcvv2(i),poj,pat,gacvv2(i),bgxvv2(i), bgyvv2
 
 4218      +  (i),bgzvv2(i),ifb1,ifb2,ifb3,ifb4, ijcvv2(i),ijcvv2(i),nobam,
 
 4229             WRITE (6,
'(A,2I5/A)') .EQ.
' HADRVV: NHKKNMXHKK ',nhkk,
nmxhkk 
 4233           ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 4234         IF (abs(ehecc-hef(j)).GT.0.001d0) 
THEN 
 4244           ptvv=ptvv+
sqrt(pxf(j)**2+pyf(j)**2)
 
 4246           IF(ibarf(j).EQ.500)istist=2
 
 4248      *                pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),2)
 
 4250           IF(idhkk(nhkk).EQ.99999) 
WRITE (6,1010)nhkk,nref(j), idhkk
 
 4252           imohkk=jmohkk(1,nhkk)
 
 4253           IF(imohkk.LE.0.OR.imohkk.GT.
nmxhkk)
THEN 
 4254             WRITE(6,
'(A,I10)')
' HADRVV out of range IMOHKK= ',i10
 
 4258       WRITE(6,*)
' From HADRVV second chain after HKKFIL' 
 4259           IF (iphkk.GE.0) 
WRITE(6,1020) nhkk, isthkk(nhkk),idhkk(nhkk),
 
 4260      +    jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 4261      +    (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 4280       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 4296       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 4298      *                ,xpsu(248),xtsu(248)
 
 4299      *                ,xpsut(248),xtsut(248)
 
 4300        common/popcck/pdbck,pdbse,pdbseu,
 
 4301      *  ijpock,irejck,ick4,ihad4,ick6,ihad6
 
 4302      *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
 
 4303      *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
 
 4304      *isea43,isea63,irejao
 
 4306       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 4307      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 4308      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 4310      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 4324       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 4330      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 4332       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 4333       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 4340       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 4343      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 4348       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 4436       COMMON /abrsv/ amcsv1(248),amcsv2(248),gacsv1(248),gacsv2(248),
 
 4437      +bgxsv1(248),bgysv1(248),bgzsv1(248), bgxsv2(248),bgysv2(248),
 
 4438      +bgzsv2(248), nchsv1(248),nchsv2(248),ijcsv1(248),ijcsv2(248),
 
 4439      +pqsva1(248,4),pqsva2(248,4), pqsvb1(248,4),pqsvb2(248,4)
 
 4442       parameter(nfimax=249)
 
 4443       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 4444      +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 4445       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 4448       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 4450       COMMON /projk/ iprojk
 
 4452       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 4455        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 4456      *                 anndv,annvd,annds,annsd,
 
 4458      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 4460      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 4463      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 4464        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 4465      *                 acouzz,acouhh,acouds,acousd,
 
 4466      *                 acoudz,acouzd,acoudi,
 
 4467      *                 acoudv,acouvd,acoucc
 
 4469       COMMON /zsea/zseaav,zseasu,anzsea
 
 4470       COMMON /casadi/casaxx,icasad
 
 4472       dimension poj(4),pat(4)
 
 4478         IF(nchsv1(i).EQ.99.AND.nchsv2(i).EQ.99) go to 50
 
 4482         IF (ipco.GE.6) 
WRITE (6,1000) ipsq(is1),ipsaq(is1),itvq(is2),
 
 4483      +  ittv1(is2),ittv2(is2), amcsv1(i),amcsv2(i),gacsv1(i),gacsv2(i),
 
 4484      +  bgxsv1(i),bgysv1(i),bgzsv1(i), bgxsv2(i),bgysv2(i),bgzsv2(i),
 
 4485      +  nchsv1(i),nchsv2(i),ijcsv1(i),ijcsv2(i), pqsva1(i,4),pqsva2
 
 4486      +  (i,4),pqsvb1(i,4),pqsvb2(i,4)
 
 4487  1000 
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
 
 4497         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 4498         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 4499         CALL 
parpt(2,pt1,pt2,3,nevt)
 
 4504         IF (ipco.GE.6)
WRITE (6,1244) poj,pat
 
 4505  1244   
FORMAT (
'  S-V QUARK-DIQUARK POJ,PAT ',8e12.3)
 
 4514        ittt = ifrovt(intsv2(i))
 
 4524         WRITE(6,*)
' SV q-qq ,IFB1,IFB2,IFB3,',
 
 4525      *  
'INTSV1=IS1,INTSV2=IS2,JIPPX,JITT',
 
 4526      *  ifb1,ifb2,ifb3,intsv1(i),intsv2(i),jippx,jitt
 
 4530     IF((nchsv1(i).NE.0))
THEN 
 4531         CALL 
hadjet(nhad,amcsv1(i),poj,pat,gacsv1(i),bgxsv1(i), bgysv1
 
 4532      +  (i),bgzsv1(i),ifb1,ifb2,ifb3,ifb4, ijcsv1(i),ijcsv1(i),4,nchsv1
 
 4535     aack=float(ick4)/float(ick4+ihad4+1)
 
 4536     IF((nchsv1(i).EQ.0))
THEN 
 4537         zseawu=
rndm(bb)*2.d0*zseaav
 
 4538         rseack=float(jitt)*pdbse+ zseawu*pdbseu
 
 4539     IF(ipco.GE.1)
WRITE(6,*)
'HADJSE JITT,RSEACK,PDBSE 3 dpmnuc3',
 
 4542         IF(
rndm(v).LE.rseack)
THEN 
 4544     IF(amcsv1(i).GT.2.3d0)
THEN 
 4546           CALL 
hadjse(nhad,amcsv1(i),poj,pat,gacsv1(i),bgxsv1(i),
 
 4548      +    (i),bgzsv1(i),ifb1,ifb2,ifb3,ifb4, ijcsv1(i),ijcsv1(i),4,
 
 4550      +    (i),3,irejss,iissqq)
 
 4551       IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JITT,',
 
 4552      *    
'RSEACK,IREJSS 3 dpmnuc3 ',
 
 4553      +    jitt,rseack,irejss
 
 4556     IF(irejss.EQ.1)irejse=irejse+1
 
 4557     IF(irejss.EQ.3)irejs3=irejs3+1
 
 4558     IF(irejss.EQ.2)irejs0=irejs0+1
 
 4559         CALL 
hadjet(nhad,amcsv1(i),poj,pat,gacsv1(i),bgxsv1(i), bgysv1
 
 4560      +  (i),bgzsv1(i),ifb1,ifb2,ifb3,ifb4, ijcsv1(i),ijcsv1(i),4,nchsv1
 
 4571         ELSEIF((ijpock.EQ.1).AND.
 
 4572      *      (aack.LE.pdbck))
THEN 
 4574         CALL 
hadjck(nhad,amcsv1(i),poj,pat,gacsv1(i),bgxsv1(i), bgysv1
 
 4575      +  (i),bgzsv1(i),ifb1,ifb2,ifb3,ifb4, ijcsv1(i),ijcsv1(i),4,nchsv1
 
 4579         CALL 
hadjet(nhad,amcsv1(i),poj,pat,gacsv1(i),bgxsv1(i), bgysv1
 
 4580      +  (i),bgzsv1(i),ifb1,ifb2,ifb3,ifb4, ijcsv1(i),ijcsv1(i),4,nchsv1
 
 4584         IF(irej.EQ.0)ick4=ick4+1
 
 4586         CALL 
hadjet(nhad,amcsv1(i),poj,pat,gacsv1(i),bgxsv1(i), bgysv1
 
 4587      +  (i),bgzsv1(i),ifb1,ifb2,ifb3,ifb4, ijcsv1(i),ijcsv1(i),4,nchsv1
 
 4605             WRITE (6,
'(A,2I5/A)') .EQ.
' HADRSV: NHKKNMXHKK ',nhkk,
nmxhkk 
 4609           ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 4610         IF (abs(ehecc-hef(j)).GT.0.001d0.AND.ibarf(j).NE.500) 
THEN 
 4611             WRITE(6,
'(2A/3I5,3E15.6)')
 
 4612      &            
' HADRSV / CHAIN 1 : CORRECT INCONSISTENT ENERGY ',
 
 4613      *            
'  NCALSV, NHKK,NREF(J), HEF(J),EHECC, AMF(J)',
 
 4614      *            ncalsv, nhkk,nref(j), hef(j),ehecc, amf(j)
 
 4619           ptsv=ptsv+
sqrt(pxf(j)**2+pyf(j)**2)
 
 4622           IF(ibarf(j).EQ.500)istist=2
 
 4624      *                pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),3)
 
 4626           IF(idhkk(nhkk).EQ.99999) 
WRITE (6,1030)nhkk,nref(j), idhkk
 
 4633       WRITE(6,*)
' HADRSV / CHAIN 1' 
 4634           IF (iphkk.GE.0) 
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
 
 4635      +    jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 4636      +    (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 4639         IF(ipco.GE.6)
WRITE(6,1644)pixu,piyu,pizu,pieu
 
 4640  1644   
FORMAT(
' HADRSV,ch1 PIXU,PIYU,PIZU,PIEU ',4f12.5)
 
 4653         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 4654         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 4655         CALL 
parpt(2,pt1,pt2,3,nevt)
 
 4658         WRITE(6,*)
' SV aq-q ,IFB1,IFB2,',
 
 4659      *  
'INTSV1=IS1,INTSV2=IS2,JIPPX,JITTX',
 
 4660      *  ifb1,ifb2,intsv1(i),intsv2(i),jippx,jittx
 
 4664         IF (ipco.GE.6)
WRITE (6,1244) poj,pat
 
 4665         CALL 
hadjet(nhad,amcsv2(i),poj,pat,gacsv2(i),bgxsv2(i), bgysv2
 
 4666      +  (i),bgzsv2(i),ifb1,ifb2,ifb3,ifb4, ijcsv2(i),ijcsv2(i),3,nchsv2
 
 4677             WRITE (6,
'(A,2I5/A)') .EQ.
' HADRSV: NHKKNMXHKK ', nhkk,
 
 4683           ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 4684         IF (abs(ehecc-hef(j)).GT.0.001d0.AND.ibarf(j).NE.500) 
THEN 
 4685               WRITE(6,
'(2A/3I5,3E15.6)')
 
 4686      &            
' HADRSV / CHAIN 2 : CORRECT INCONSISTENT ENERGY ',
 
 4687      *            
'  NCALSV, NHKK,NREF(J), HEF(J),EHECC, AMF(J)',
 
 4688      *            ncalsv, nhkk,nref(j), hef(j),ehecc, amf(j)
 
 4693           ptsv=ptsv+
sqrt(pxf(j)**2+pyf(j)**2)
 
 4696           IF(ibarf(j).EQ.500)istist=2
 
 4698      *                pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),4)
 
 4699           IF(idhkk(nhkk).EQ.99999) 
WRITE (6,1030)nhkk,nref(j), idhkk
 
 4705           IF (iphkk.GE.7) 
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
 
 4706      +    jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 4707      +    (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 4710         IF(ipco.GE.6)
WRITE(6,1644)pixu,piyu,pizu,pieu
 
 4719  1010 
FORMAT (i6,i4,5i6,9e10.2)
 
 4720  1020 
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
 
 4721  1030 
FORMAT (
' NHKK,IDHKK(NHKK)  ',3i10)
 
 4728       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 4735       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 4737      *                ,xpsu(248),xtsu(248)
 
 4738      *                ,xpsut(248),xtsut(248)
 
 4740       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 4741      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 4742      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 4744      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 4758       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 4764      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 4766       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 4767       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 4774       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 4777      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 4782       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 4877       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 4880       parameter(nfimax=249)
 
 4881       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 4882      +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 4883       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 4886       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 4888       COMMON /projk/ iprojk
 
 4890       dimension poj(4),pat(4)
 
 4892        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 4893      *                 anndv,annvd,annds,annsd,
 
 4895      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 4897      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 4900      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 4901        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 4902      *                 acouzz,acouhh,acouds,acousd,
 
 4903      *                 acoudz,acouzd,acoudi,
 
 4904      *                 acoudv,acouvd,acoucc
 
 4906       COMMON /pshow/ ipshow
 
 4908       COMMON /harlun/ qlun,iharlu
 
 4909       COMMON /jspar/pxj(1000),pyj(1000),pzj(1000),hej(1000),nnnpj
 
 4910       COMMON /jspa/pxs(40000),pys(40000),pzs(40000),hes(40000),nnnps
 
 4911       COMMON /nomije/ ptmije(10),nnmije(10)
 
 4912       COMMON /casadi/casaxx,icasad
 
 4916         IF(nchss1(i).EQ.99.AND.nchss2(i).EQ.99) go to 60
 
 4921           IF (ipco.GE.6) 
WRITE (6,1000) ipsq(is1),ipsaq(is1),itsq(is2),
 
 4922      +    itsaq(is2), amcss1(i),amcss2(i),gacss1(i),gacss2(i), bgxss1
 
 4923      +    (i),bgyss1(i),bgzss1(i), bgxss2(i),bgyss2(i),bgzss2(i), nchss1
 
 4924      +    (i),nchss2(i),ijcss1(i),ijcss2(i), pqssa1(i,4),pqssa2(i,4),
 
 4925      +    pqssb1(i,4),pqssb2(i,4)
 
 4926  1000 
FORMAT(10
x,4i5,10f9.2/10
x,4i5,4f12.4)
 
 4936         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 4937         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 4938         CALL 
parpt(2,pt1,pt2,4,nevt)
 
 4943               pojpt=
sqrt(poj(2)**2+poj(1)**2)
 
 4944               patpt=
sqrt(pat(1)**2+pat(2)**2)
 
 4946               IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
 
 4948               IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
 
 4951               qlun=min(pojpt,patpt)
 
 4952               IF((qlun.LT.2.5d0).OR.(amcss1(i).LT.5.d0))
THEN 
 4960         WRITE(6,*)
' SS q-aq ,IFB1,IFB2,',
 
 4961      *  
'INTSS1=IS1,INTSS2=IS2',
 
 4962      *  ifb1,ifb2,intss1(i),intss2(i)
 
 4963         WRITE (6,*)
' projectile sea quark IFB1=',ifb1,
 
 4964      *  
' from IS1=',intss1(i)
 
 4965         WRITE(6,*)
' with IPSQ(IS1),XPSQ(IS1),IFROSP(IS1)',
 
 4966      *  ipsq(is1),xpsq(is1),ifrosp(is1)
 
 4969       IF(ifrosp(is1).EQ.ifrovp(ii))iii=ii
 
 4972         WRITE (6,*)
' projectile III=',iii
 
 4973         WRITE(6,*)
' corresp. XPVQ(i),XPVD(i),IPVQ(I),IPPV1(I),IPPV2(I)',
 
 4974      *   xpvq(iii),xpvd(iii),ipvq(iii),ippv1(iii),ippv2(iii)
 
 4981          IF(
rndm(vv).LE.casaxx)
THEN 
 4982        IF(
rndm(vvv).LE.0.5d0)
THEN 
 4989         WRITE(6,*)
' Cas SS1 q-aq 1 ,IFB1,IFB2,',
 
 4990      *  
'INTSS1=IS1,INTSS2=IS2,III',
 
 4991      *  ifb1,ifb2,intss1(i),intss2(i),iii
 
 4992      *  ,
'-----------------------------------------------------' 
 5001         WRITE(6,*)
' Cas SS1 q-aq 2 ,IFB1,IFB2,',
 
 5002      *  
'INTSS1=IS1,INTSS2=IS2,III',
 
 5003      *  ifb1,ifb2,intss1(i),intss2(i),iii
 
 5004      *  ,
'-----------------------------------------------------' 
 5012           CALL 
hadjet(nhad,amcss1(i),poj,pat,gacss1(i),bgxss1(i), bgyss1
 
 5013      +    (i),bgzss1(i),ifb1,ifb2,ifb3,ifb4, ijcss1(i),ijcss1(i),3,
 
 5022             ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 5023           IF (abs(ehecc-hef(j)).GT.0.001d0.AND.ibarf(j).NE.500) 
THEN 
 5024               WRITE(6,
'(A,2I5,2E16.6)')
 
 5025      +        
' HADRSS: CORRECT INCONSISTENT PARTICLE ENERGY ', nhkk,
 
 5026      +        nref(j), hef(j),ehecc
 
 5031           ptss=ptss+
sqrt(pxf(j)**2+pyf(j)**2)
 
 5035               WRITE (6,
'(A,2I5)') 
' HADRSS: NHKK.EQ NMXHKK',nhkk,
nmxhkk 
 5039           IF(ibarf(j).EQ.500)istist=2
 
 5041      *                pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),5)
 
 5042             IF(idhkk(nhkk).EQ.99999) 
WRITE (6,1030) nhkk,nref(j), idhkk
 
 5045             IF (iphkk.GE.7) 
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk
 
 5046      +      (nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk
 
 5047      +      (2,nhkk),(phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk
 
 5059             DO 137 j=nnnps,nnnpsu
 
 5061               IF(j.GT.40000.OR.jj.GT.1000)
THEN 
 5081         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 5082         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 5083         CALL 
parpt(2,pt1,pt2,4,nevt)
 
 5087               pojpt=
sqrt(poj(2)**2+poj(1)**2)
 
 5088               patpt=
sqrt(pat(1)**2+pat(2)**2)
 
 5090               IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
 
 5092               IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
 
 5095               qlun=min(pojpt,patpt)
 
 5096               IF((qlun.LT.2.5d0).OR.(amcss2(i).LT.5.d0))
THEN 
 5105         WRITE(6,*)
' SS aq-q ,IFB1,IFB2,',
 
 5106      *  
'INTSS1=IS1,INTSS2=IS2',
 
 5107      *  ifb1,ifb2,intss1(i),intss2(i)
 
 5108         WRITE (6,*)
' target sea quark IFB2=',ifb2,
 
 5109      *  
' from IS2=',intss2(i)
 
 5110         WRITE(6,*)
' with ITSQ(IS2),XTSQ(IS2),IFROST(IS2)',
 
 5111      *  itsq(is2),xtsq(is2),ifrost(is2)
 
 5114       IF(ifrost(is2).EQ.ifrovt(ii))iii=ii
 
 5117         WRITE (6,*)
' projectile III=',iii
 
 5118         WRITE(6,*)
' corresp. XTVQ(i),XTVD(i),ITVQ(I),ITTV1(I),ITTV2(I)',
 
 5119      *   xtvq(iii),xtvd(iii),itvq(iii),ittv1(iii),ittv2(iii)
 
 5126          IF(
rndm(vv).LE.casaxx)
THEN 
 5127        IF(
rndm(vvv).LE.0.5d0)
THEN 
 5134         WRITE(6,*)
' Cas SS2 aq-q 1 ,IFB1,IFB2,',
 
 5135      *  
'INTSS1=IS1,INTSS2=IS2,III',
 
 5136      *  ifb1,ifb2,intss1(i),intss2(i),iii
 
 5137      *  ,
'-----------------------------------------------------' 
 5146         WRITE(6,*)
' Cas SS2 aq-q 2 ,IFB1,IFB2,',
 
 5147      *  
'INTSS1=IS1,INTSS2=IS2,III',
 
 5148      *  ifb1,ifb2,intss1(i),intss2(i),iii
 
 5149      *  ,
'-----------------------------------------------------' 
 5157           CALL 
hadjet(nhad,amcss2(i),poj,pat,gacss2(i),bgxss2(i), bgyss2
 
 5158      +    (i),bgzss2(i),ifb1,ifb2,ifb3,ifb4, ijcss2(i),ijcss2(i),3,
 
 5166             ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 5167           IF (abs(ehecc-hef(j)).GT.0.001d0.AND.ibarf(j).NE.500) 
THEN 
 5168               WRITE(6,
'(A,2I5,2E16.6)')
 
 5169      +        
' HADRSS: CORRECT INCONSISTENT PARTICLE ENERGY ', nhkk,
 
 5170      +        nref(j), hef(j),ehecc
 
 5175           ptss=ptss+
sqrt(pxf(j)**2+pyf(j)**2)
 
 5179               WRITE (6,
'(A,2I5)') 
' HADRSS: NHKK.EQ NMXHKK',nhkk,
nmxhkk 
 5183           IF(ibarf(j).EQ.500)istist=2
 
 5185      *                pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),6)
 
 5186             IF(idhkk(nhkk).EQ.99999) 
WRITE (6,1030) nhkk,nref(j), idhkk
 
 5189             IF (iphkk.GE.7) 
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk
 
 5190      +      (nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk
 
 5191      +      (2,nhkk),(phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk
 
 5202             DO 187 j=nnnps,nnnpsu
 
 5204               IF(j.GT.40000.OR.jj.GT.1000)
THEN 
 5221  1010 
FORMAT (i6,i4,5i6,9e10.2)
 
 5222  1020 
FORMAT (.GT.
' HADRVS JNAUMAX SKIP NEXT PARTICLES ',3i10)
 
 5223  1030 
FORMAT (
' NHKK,NREF(J),IDHKK(NHKK)  ',3i10)
 
 5224  1040 
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
 
 5231       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 5248       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 5250      *                ,xpsu(248),xtsu(248)
 
 5251      *                ,xpsut(248),xtsut(248)
 
 5252        common/popcck/pdbck,pdbse,pdbseu,
 
 5253      *  ijpock,irejck,ick4,ihad4,ick6,ihad6
 
 5254      *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
 
 5255      *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
 
 5256      *isea43,isea63,irejao
 
 5258       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 5259      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 5260      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 5262      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 5276       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 5282      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 5284       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 5285       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 5292       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 5295      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 5300       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 5388       COMMON /abrvs/ amcvs1(248),amcvs2(248),gacvs1(248),gacvs2(248),
 
 5389      +bgxvs1(248),bgyvs1(248),bgzvs1(248), bgxvs2(248),bgyvs2(248),
 
 5390      +bgzvs2(248), nchvs1(248),nchvs2(248),ijcvs1(248),ijcvs2(248),
 
 5391      +pqvsa1(248,4),pqvsa2(248,4), pqvsb1(248,4),pqvsb2(248,4)
 
 5394       parameter(nfimax=249)
 
 5395       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 5396      +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 5397       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 5400       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 5402       COMMON /projk/ iprojk
 
 5404       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 5407        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 5408      *                 anndv,annvd,annds,annsd,
 
 5410      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 5412      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 5415      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 5416        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 5417      *                 acouzz,acouhh,acouds,acousd,
 
 5418      *                 acoudz,acouzd,acoudi,
 
 5419      *                 acoudv,acouvd,acoucc
 
 5421       COMMON /zsea/zseaav,zseasu,anzsea
 
 5422       COMMON /casadi/casaxx,icasad
 
 5424       dimension poj(4),pat(4)
 
 5428         IF(nchvs1(i).EQ.99.AND.nchvs2(i).EQ.99) go to 50
 
 5432         IF (ipco.GE.6) 
WRITE (6,1010) ipvq(is1),ippv1(is1),ippv2(is1),
 
 5433      +  itsq(is2),itsaq(is2), amcvs1(i),amcvs2(i),gacvs1(i),gacvs2(i),
 
 5434      +  bgxvs1(i),bgyvs1(i),bgzvs1(i), bgxvs2(i),bgyvs2(i),bgzvs2(i),
 
 5435      +  nchvs1(i),nchvs2(i),ijcvs1(i),ijcvs2(i), pqvsa1(i,4),pqvsa2
 
 5436      +  (i,4),pqvsb1(i,4),pqvsb2(i,4)
 
 5446         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 5447         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 5448         CALL 
parpt(2,pt1,pt2,2,nevt)
 
 5456         WRITE(6,*)
' VS q-aq ,IFB1,IFB2,',
 
 5457      *  
'INTVS1=IS1,INTVS2=IS2,JIPPX,JITTX',
 
 5458      *  ifb1,ifb2,intvs1(i),intvs2(i),jippx,jittx
 
 5460         CALL 
hadjet(nhad,amcvs1(i),poj,pat,gacvs1(i),bgxvs1(i), bgyvs1
 
 5461      +  (i),bgzvs1(i),ifb1,ifb2,ifb3,ifb4, ijcvs1(i),ijcvs1(i),3,nchvs1
 
 5468           ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 5469         IF (abs(ehecc-hef(j)).GT.0.001d0.AND.ibarf(j).NE.500)
THEN 
 5470             WRITE(6,
'(A,2I5,2E16.6)')
 
 5471      +      
' HADRVS: CORRECT INCONSISTENT PARTICLE ENERGY ', nhkk,nref
 
 5477           ptvs=ptvs+
sqrt(pxf(j)**2+pyf(j)**2)
 
 5481             WRITE (6,
'(A,2I5)') .EQ.
' HADRVS: NHKKNMXHKK',nhkk,
nmxhkk 
 5485           IF(ibarf(j).EQ.500)istist=2
 
 5487      *                pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),7)
 
 5488           IF(idhkk(nhkk).EQ.99999) 
WRITE (6,1030) nhkk, idhkk(nhkk)
 
 5490           IF (iphkk.GE.7) 
WRITE(6,1000)nhkk, isthkk(nhkk),idhkk(nhkk),
 
 5491      +    jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 5492      +    (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 5508         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 5509         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 5510         CALL 
parpt(2,pt1,pt2,2,nevt)
 
 5519        ippp = ifrovp(intvs1(i))
 
 5529         WRITE(6,*)
' VS qq-q ,IFB1,IFB2,IFB3,',
 
 5530      *  
'INTVS1=IS1,INTVS2=IS2,JIPP,JITTX',
 
 5531      *  ifb1,ifb2,ifb3,intvs1(i),intvs2(i),jipp,jittx
 
 5535     IF((nchvs2(i).NE.0))
THEN 
 5536         CALL 
hadjet(nhad,amcvs2(i),poj,pat,gacvs2(i),bgxvs2(i), bgyvs2
 
 5537      +  (i),bgzvs2(i),ifb1,ifb2,ifb3,ifb4, ijcvs2(i),ijcvs2(i),6,nchvs2
 
 5540     aack=float(ick6)/float(ick6+ihad6+1)
 
 5541     IF((nchvs2(i).EQ.0))
THEN 
 5542           zseawu=
rndm(bb)*2.d0*zseaav
 
 5543         rseack=float(jipp)*pdbse+ zseawu*pdbseu
 
 5544       IF(ipco.GE.1)
WRITE(6,*)
'HADJSE JIPP,RSEACK,PDBSE 4 dpmnuc3',
 
 5547       IF(
rndm(v).LE.rseack)
THEN 
 5549       IF(amcvs2(i).GT.2.3d0)
THEN 
 5551             CALL 
hadjse(nhad,amcvs2(i),poj,pat,gacvs2(i),bgxvs2(i),
 
 5553      +      (i),bgzvs2(i),ifb1,ifb2,ifb3,ifb4, ijcvs2(i),ijcvs2(i),6,
 
 5555      +      (i),6,irejss,iissqq)
 
 5556         IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JIPP,',
 
 5557      *      
'RSEACK,IREJSS 4 dpmnuc3',
 
 5558      +      jipp,rseack,irejss
 
 5561         IF(irejss.EQ.1)irejse=irejse+1
 
 5562         IF(irejss.EQ.3)irejs3=irejs3+1
 
 5563         IF(irejss.EQ.2)irejs0=irejs0+1
 
 5564             CALL 
hadjet(nhad,amcvs2(i),poj,pat,gacvs2(i),
 
 5566      +      (i),bgzvs2(i),ifb1,ifb2,ifb3,ifb4, ijcvs2(i),
 
 5567      *      ijcvs2(i),6,nchvs2
 
 5578         ELSEIF((ijpock.EQ.1).AND.
 
 5579      *      (aack.LE.pdbck))
THEN 
 5581         CALL 
hadjck(nhad,amcvs2(i),poj,pat,gacvs2(i),bgxvs2(i), bgyvs2
 
 5582      +  (i),bgzvs2(i),ifb1,ifb2,ifb3,ifb4, ijcvs2(i),ijcvs2(i),6,nchvs2
 
 5586         CALL 
hadjet(nhad,amcvs2(i),poj,pat,gacvs2(i),bgxvs2(i), bgyvs2
 
 5587      +  (i),bgzvs2(i),ifb1,ifb2,ifb3,ifb4, ijcvs2(i),ijcvs2(i),6,nchvs2
 
 5591     IF(irej.EQ.0)ick6=ick6+1
 
 5593         CALL 
hadjet(nhad,amcvs2(i),poj,pat,gacvs2(i),bgxvs2(i), bgyvs2
 
 5594      +  (i),bgzvs2(i),ifb1,ifb2,ifb3,ifb4, ijcvs2(i),ijcvs2(i),6,nchvs2
 
 5605           ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 5606         IF (abs(ehecc-hef(j)).GT.0.001d0.AND.ibarf(j).NE.500) 
THEN 
 5607             WRITE(6,
'(A,2I5,2E16.6)')
 
 5608      +      
' HADRVS: CORRECT INCONSISTENT PARTICLE ENERGY ', nhkk,nref
 
 5614           ptvs=ptvs+
sqrt(pxf(j)**2+pyf(j)**2)
 
 5618             WRITE (6,
'(A,2I5)') .EQ.
' HADRVS: NHKKNMXHKK ',nhkk,
nmxhkk 
 5622           IF(ibarf(j).EQ.500)istist=2
 
 5624      *                pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),8)
 
 5626           IF(idhkk(nhkk).EQ.99999) 
WRITE (6,1030)nhkk,nref(j), idhkk
 
 5629       WRITE(6,*)
' Second chain HADRVS' 
 5630           IF (iphkk.GE.0) 
WRITE(6,1000) nhkk, isthkk(nhkk),idhkk(nhkk),
 
 5631      +    jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 5632      +    (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 5642  1000 
FORMAT (i6,i4,5i6,9e10.2)
 
 5643  1010 
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
 
 5644  1020 
FORMAT (.GT.
' HADRVS JNAUMAX SKIP NEXT PARTICLES ',3i10)
 
 5645  1030 
FORMAT (
' NHKK,IDHKK(NHKK)  ',3i10)
 
 5651       SUBROUTINE hadjet(NHAD,AMCH,PPR,PTA,GAM,BGX,BGY,BGZ, IFB1,IFB2,
 
 5652      +ifb3,ifb4,i1,i2,nobam,nnch,norig)
 
 5653       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 5709       parameter(nfimax=249)
 
 5710       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 5711      +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 5712       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 5715       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 5719       COMMON /jspart/pxp(1000),
pyp(1000),pzp(1000),hepp(1000),nnnp
 
 5720       COMMON /jspar/pxj(1000),pyj(1000),pzj(1000),hej(1000),nnnpj
 
 5724       dimension ppr(4),pta(4),pprj(4),ptaj(4)
 
 5726       parameter(tiny=1.
d-10)
 
 5736         WRITE(6,1010) gam,bgx,bgy,bgz,ppr,pprj
 
 5737         WRITE(6,1000) nhad,amch,pta, ifb1,ifb2,ifb3,ifb4,i1,i2,nobam,
 
 5739  1000 
FORMAT(10
x,i10,5f10.3/10
x,9i10)
 
 5741       IF(abs(nnch).EQ.99) 
THEN 
 5747       CALL 
daltra(gam,-bgx,-bgy,-bgz,ppr(1),ppr(2),ppr(3),ppr(4),pprtot,
 
 5748      +pprj(1),pprj(2),pprj(3),pprj(4))
 
 5749       CALL 
daltra(gam,-bgx,-bgy,-bgz,pta(1),pta(2),pta(3),pta(4),ptatot,
 
 5750      +ptaj(1),ptaj(2),ptaj(3),ptaj(4))
 
 5752       IF(ipco.GE.3) 
WRITE(6,1010)gam,bgx,bgy,bgz,ppr,pprj,pta,ptaj
 
 5753  1010 
FORMAT(
' HADJET: GAM,BGX,BGY,BGZ,PPR(4),PPRJ(4) ',4f15.5/8f15.5/ 8
 
 5757       IF(pprtot.LT.tiny)pprtot=tiny
 
 5759       IF(cod.GE.1.d0)cod=0.999999999999
 
 5760       IF(cod.LE.-1.d0)cod=-0.999999999999
 
 5761       sid=
sqrt(abs((1.d0-cod)*(1.d0+cod)))
 
 5764       IF((abs(pprj(1)).GT.0.d0).OR.(abs(pprj(2)).GT.0.d0))
THEN 
 5765       IF(pprtot*sid.GT.1.
d-9) 
THEN 
 5766         cof=pprj(1)/(sid*pprtot)
 
 5767         sif=pprj(2)/(sid*pprtot)
 
 5768         anorf=
sqrt(abs(cof*cof+sif*sif))
 
 5773       IF (ipco.GE.6) 
WRITE(6,1020)cod,sid,cof,sif
 
 5774  1020 
FORMAT(
' COD,SID,COF,SIF ',4f15.8)
 
 5776       CALL 
calbam(nnch,i1,i2,ifb1,ifb2,ifb3,ifb4,amch,nobam,ihad)
 
 5793         IF(ibarf(i).EQ.500)go to 1011
 
 5794         pxcal=pxcal + pxf(i)
 
 5795         pycal=pycal + pyf(i)
 
 5796         pzcal=pzcal + pzf(i)
 
 5798         ehecc=
sqrt(abs(pxf(i)**2+pyf(i)**2+pzf(i)**2+amf(i)**2))
 
 5799       IF (abs(ehecc-hef(i)).GT.0.0001d0) 
THEN 
 5800         IF(ipri.GE.1) 
WRITE(6,
'(2A/A/3I5,3E16.6)')
 
 5801      +     
' HADJET / AFTER CALBAM:',
 
 5802      +     
'  CORRECT INCONSISTENT ENERGY IN JET CMS',
 
 5803      +            
'  I, IHAD,NREF(I), HEF(I),EHECC, AMF(I)',
 
 5804      *            i,ihad,nref(i), hef(i),ehecc, amf(i)
 
 5812       IF (abs(ecal-amch).GT.0.005d0) ltesha=.true.
 
 5815       IF (abs(ecal-amch).GT.0.005d0) 
THEN 
 5816         IF(icheca.LE.10)
WRITE(6,
'(A/10I4)')
 
 5817      +  
' HADJET/1:ICHECA,IFB1,IFB2,IFB3,IFB4,I1,I2,NOBAM,NNCH,NORIG',
 
 5819      +  ifb2,ifb3,ifb4,i1,i2,nobam,nnch,norig
 
 5820         IF(icheca.LE.10)
WRITE(6,1030) amch,ecal,pxcal,pycal,pzcal
 
 5821  1030 
FORMAT(
' CALBAM E. CHECK (5 MeV) AMCH,ECAL,PXCAL,PYCAL,PZCAL=',
 
 5828         WRITE(6,1040)i,pxf(i),pyf(i),pzf(i),hef(i),amf(i), ichf(i),
 
 5829      +    ibarf(i),nref(i),anf(i)
 
 5830  1040 
FORMAT(
' JET SYSTEM ',i5,5f12.4,3i5,a10)
 
 5864       phec2=pxf(i)**2+pyf(i)**2+pzf(i)**2
 
 5865         CALL 
dtrans(pxf(i),pyf(i),pzf(i),cod,sid,cof,sif,
xx,
yy,
zz)
 
 5866       prota2=
xx**2 + 
yy**2 + 
zz**2
 
 5870       pxcal=pxcal + pxf(i)
 
 5871       pycal=pycal + pyf(i)
 
 5872       pzcal=pzcal + pzf(i)
 
 5874       IF(abs(phec2-prota2).GT.0.0001d0) 
THEN 
 5875         WRITE(6,
'(2A/3I5,3E16.6)')
 
 5876      &            
' HADJET: INCONSISTENT MOMENTUM AFTER TRANS',
 
 5877      *            
'  I, IHAD,NREF(I), PHEC2,PROTA2, AMF(I)',
 
 5878      *            i,ihad,nref(i), phec2,prota2, amf(i)
 
 5885      +  
' HADJET/2: IFB1,IFB2,IFB3,IFB4,I1,I2,NOBAM,NNCH,NORIG', ifb1,
 
 5886      +  ifb2,ifb3,ifb4,i1,i2,nobam,nnch,norig
 
 5887       WRITE(6,1031) pxcal,pycal,pzcal
 
 5888  1031   
FORMAT(
' CALBAM ENERGY CHECK/2:  PXCAL,PYCAL,PZCAL='/3e20.8)
 
 5893         WRITE(6,1050)i,pxf(i),pyf(i),pzf(i),hef(i),amf(i), ichf(i),
 
 5894      +    ibarf(i),nref(i),anf(i)
 
 5895  1050 
FORMAT(
' ROTATET JET SYSTEM ',i5,5f12.4,3i5,a10)
 
 5902       IF(nnnp.GT.1000)nnnnp=1000
 
 5915       CALL 
lortmo(ihad,gam,bgx,bgy,bgz)
 
 5920         WRITE(6,1060) i,pxf(i),pyf(i),pzf(i),hef(i),amf(i), ichf(i),
 
 5921      +    ibarf(i),nref(i),anf(i)
 
 5922  1060 
FORMAT(
' CMS SYSTEM ',i5,5f12.4,3i5,a10)
 
 5930       IF(ibarf(i).EQ.500)go to 6060
 
 5932         ehecc=
sqrt(abs(pxf(i)**2+pyf(i)**2+pzf(i)**2+amf(i)**2))
 
 5933       IF (abs(ehecc-hef(i)).GT.0.001d0) 
THEN 
 5934         IF(abs(ehecc-hef(i)).GT.0.1d0)
WRITE(6,
'(2A/4I5,3E16.6)')
 
 5935      &            
' HADJET: CORRECT INCONSISTENT ENERGY AFTER LORTRA',
 
 5936      *            
'  NORIG, I, IHAD,NREF(I), HEF(I),EHECC, AMF(I)',
 
 5937      *            norig, i,ihad,nref(i), hef(i),ehecc, amf(i)
 
 5947       IF(nnnp.GT.1000)nnnnp=1000
 
 5949       CALL 
lortrp(nnnnp,1,gam,bgx,bgy,bgz)
 
 5961       IF (abs(eein-ehad).GT.0.005d0) 
THEN 
 5963       IF (abs(eein-ehad).GT.0.005d0) 
THEN 
 5964     IF(icheco.LT.10)
THEN 
 5966      +  
' HADJET/3:ICHECO,IFB1,IFB2,IFB3,IFB4,I1,I2,NOBAM,NNCH,NORIG', 
 
 5968      +  ifb2,ifb3,ifb4,i1,i2,nobam,nnch,norig
 
 5969       WRITE (6,1070) eein,ehad,amch,gam,bgx,bgy,bgz
 
 5970  1070 
FORMAT(
' HADJET ENERGY CHECK (5 MeV) EEIN,EHAD,AMCH',3e20.8/
 
 5971      +  20
x,
' GAM,BGX,BGY,BGZ ',4e20.8)
 
 5990       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 5993       parameter(mxnupa=2500)
 
 5994       COMMON /jspart/pxp(1000),
pyp(1000),pzp(1000),hepp(1000),nnnp
 
 5995       COMMON /jspar/pxj(1000),pyj(1000),pzj(1000),hej(1000),nnnpj
 
 6017       CALL 
daltra(gam,bgx,bgy,bgz,pxp(i),
pyp(i),pzp(i),hepp(i),
 
 6018      *ppa,pxj(j),pyj(j),pzj(j),hej(j))
 
 6025       CALL 
daltra(gam,bgx,bgy,bgz,pxsm,pysm,pzsm,esum,
 
 6026      *ppa,pxsm,pysm,pzsm,esum)
 
 6033       diffl=pxdif+pydif+pzdif+edif
 
 6034       IF(diffl.GE.1.
d-2*
one)
 
 6035      1
WRITE(6,2)num,pxdif,pydif,pzdif,edif,pxsm,pxsc,
 
 6037      1pysm,pysc,pzsm,pzsc,esum,esmc
 
 6038  2    
FORMAT(
' ',2
x,
'LORTRA:NUM=',i5,2
x,
'PXDIF=',1pe15.6,2
x,
'PYDIF=',
 
 6039      21pe15.6,2
x,
'PZDIF=',1pe15.6,2
x,
'EDIF=',1pe15.6/2
x,
'PXSM=',1pe15.6,
 
 6040      32
x,
'PXSC=',1pe15.6,2
x,
'PYSM=',1pe15.6,2
x,
'PYSC=',1pe15.6/2
x,
'PZSM' 
 6041      4,1pe15.6,2
x,
'PZSC=',1pe15.6,2
x,
'ESUM=',1pe15.6,2
x,
'ESMC=',1pe15.6/
 
 6042      52
x,
'LORTRA DIFFERENCES DUE TO ALTRA'/)
 
 6049       SUBROUTINE cobcma(IF1,IF2,IF3,IJNCH,NNCH,IREJ,AMCH,AMCHN,IKET)
 
 6050       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6074       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 6075      +iibar(210),k1(210),k2(210)
 
 6078       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 6080       COMMON /ketmas/ am8(2),am10(2),ib88(2),ib1010(2),amch1n,amch2n
 
 6083       CALL 
dbklas(if1,if2,if3,ib8,ibb10)
 
 6085       IF (ipev.GE.2)
WRITE(6,1000)if1,if2,if3,ib8,ibb10
 
 6086  1000 
FORMAT (
' COBCMA: IPQ,ITTQ1,ITTQ2,IB8,IBB10 ',5i5)
 
 6105       IF(amch.LT.am81) 
THEN 
 6107       ELSEIF (amch.LT.am101)
THEN 
 6113       ELSEIF(amch.LT.amff1) 
THEN 
 6124         WRITE(6,1010) amch,amchn,am81,am101
 
 6125         WRITE(6,1020) if1,if2,if3,ib8,ibb10,ijnch,nnch,irej
 
 6126  1010 
FORMAT(
' COBCMA: AMCH,AMCHN,AM81,AM101', 4f13.4)
 
 6127  1020 
FORMAT(
' COBCMA: IF1,IF2,IF3,IB8,IBB10,IJNCH,NNCH,IREJ',8i4)
 
 6135       SUBROUTINE comcma(IFQ,IFAQ,IJNCH,NNCH,IREJ,AMCH,AMCHN)
 
 6136       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6156       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 6157      +iibar(210),k1(210),k2(210)
 
 6160       COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
 
 6161      +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
 
 6163       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 6167       ifps=imps(iifaq,ifq)
 
 6169       IF (ipev.GE.2)
WRITE (6,1000)iifaq,ifq,ifps,ifv
 
 6170  1000 
FORMAT (
' COMCMA',5
x,
' IIPPAQ,ITQ,IFPS,IFV ',4i5)
 
 6177       IF(ipev.GE.2) 
WRITE(6,1010) amch,amps,amv,ifps,ifv
 
 6178  1010 
FORMAT(
' AMCH,AMPS,AMV,IFPS,IFV ',3f12.4,2i10)
 
 6186       IF(amch.LT.amps) 
THEN 
 6191       IF (amch.LT.amv) 
THEN 
 6196       ELSEIF(amch.LT.amff) 
THEN 
 6205         WRITE(6,1030) amch,amchn,amps,amv
 
 6206         WRITE(6,1020) ifq,ifaq,ifps,ifv,ijnch,nnch,irej
 
 6207  1030 
FORMAT(
' COMCMA: AMCH,AMCHN,AMPS,AMV', 4f13.4)
 
 6208  1020 
FORMAT(
' COMCMA: IFQ,IFAQ,IFPS,IFV,IJNCH,NNCH,IREJ',8i4)
 
 6218       SUBROUTINE comcm2(IQ1,IQ2,IAQ1,IAQ2,NNCH,IREJ,AMCH)
 
 6219       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6240       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 6241      +iibar(210),k1(210),k2(210)
 
 6244       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 6246       COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
 
 6247      +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
 
 6253       IF (iiaq1.EQ.iq1)                                         go to 10
 
 6254       IF (iiaq1.EQ.iq2)                                         go to 20
 
 6255       IF (iiaq2.EQ.iq1)                                         go to 30
 
 6256       IF (iiaq2.EQ.iq2)                                         go to 40
 
 6261         WRITE(6,
'(A/5X,4I5,1PE13.5)')
 
 6262      +  
' KKEVVV/COMCM2 (QU. NUMBERS): IQ1, IQ2, IAQ1, IAQ2, AMCH', iq1,
 
 6263      +  iq2, iaq1, iaq2, amch
 
 6291       IF (amch.LT.amff) 
THEN 
 6294           WRITE(6,
'(A/5X,4I5,1PE13.5)')
 
 6295      +    
' KKEVVV/COMCM2 (MASS!): IQ1, IQ2, IAQ1, IAQ2, AMCH', iq1,
 
 6296      +    iq2, iaq1, iaq2, amch
 
 6306      +pq1x,pq1y,pq1z,pq1e,pa1x,pa1y,pa1z,pa1e, pq2x,pq2y,pq2z,pq2e,pa2x,
 
 6307      +pa2y,pa2z,pa2e, pxch1,pych1,pzch1,ech1, pxch2,pych2,pzch2,ech2,
 
 6309       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6321       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 6328       IF(amch1.EQ.0.d0)
THEN 
 6330         WRITE(6,*) 
'Error in CORMOM : AMCH1=0. Event rejected' 
 6375       pa1x=pa1xol+pq1xol-pq1x   
 
 6376       pa1y=pa1yol+pq1yol-pq1y   
 
 6377       pa1z=pa1zol+pq1zol-pq1z   
 
 6378       pa1e=pa1eol+pq1eol-pq1e   
 
 6379       pq2x=pq2xol+pa2xol-pa2x   
 
 6380       pq2y=pq2yol+pa2yol-pa2y   
 
 6381       pq2z=pq2zol+pa2zol-pa2z   
 
 6382       pq2e=pq2eol+pa2eol-pa2e   
 
 6393       root =(ech1-amch1)*(ech1+amch1)
 
 6394       IF(root.LT.0.d0)
THEN 
 6396         WRITE(6,*)
'Error in CORMOM : ROOT<0. Event rejected' 
 6397         WRITE(6,*)
'ECH1=',ech1,
'   AMCH1=',amch1,
'   ROOT=',root
 
 6400       pch1 = 
sqrt(root) + 0.000001
 
 6407       pch2 =
sqrt(pxch2**2+pych2**2+pzch2**2)
 
 6408       amch22=ech2**2-pxch2**2-pych2**2-pzch2**2
 
 6410       IF(amch22.LT.0.d0)
THEN 
 6412         WRITE(6,*)
'Error in CORMOM : AMCH22<0. Event rejected' 
 6420         pxsum=pq1x+pa1x+pq2x+pa2x
 
 6421         pysum=pq1y+pa1y+pq2y+pa2y
 
 6422         pzsum=pq1z+pa1z+pq2z+pa2z
 
 6423         pesum=pq1e+pa1e+pq2e+pa2e
 
 6424         WRITE(6,
'(A)') 
' CORMOM: KINEMATIC TEST FOR PARTONS' 
 6425         WRITE(6,
'(A,4(1PE12.5))') 
' PXSUM,PYSUM,PZSUM,PESUM', pxsum,
 
 6436      +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
 
 6437      +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,nselpt)
 
 6438       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6443       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 6445       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
 6447       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
 6448      +ipadis,ishmal,lpauli
 
 6450       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 6474       IF (ikvala.EQ.1)b33=16.0
 
 6475       es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
 
 6476       hps=
sqrt(es*es+2.*es*0.94)
 
 6478       IF (.NOT.intpt) hps=0.0000001
 
 6483       IF (icount.EQ.48)
THEN 
 6486       IF (icount.EQ.50)
THEN 
 6491       IF (icount.GE.1)
THEN 
 6494         ptxsq1=qtxsq1+hps*cfe
 
 6495         ptysq1=qtysq1+hps*sfe
 
 6496         ptxsa1=qtxsa1-hps*cfe
 
 6497         ptysa1=qtysa1-hps*sfe
 
 6501       es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
 
 6502       hps=
sqrt(es*es+2.*es*0.94)
 
 6505       IF (.NOT.intpt) hps=0.0000001
 
 6509       ptxsq1=qtxsq1+hps*cfe
 
 6510       ptysq1=qtysq1+hps*sfe
 
 6511       ptxsa1=qtxsa1-hps*cfe
 
 6512       ptysa1=qtysa1-hps*sfe
 
 6515       IF (ipev.GE.7)
WRITE(6,1000)ptxsq1,ptysq1,ptxsa1,ptysa1 ,ptxsq2,
 
 6516      +ptysq2,ptxsa2,ptysa2
 
 6517  1000 
FORMAT (
' PT S  ',8f12.6)
 
 6519       pttq1=ptxsq1**2+ptysq1**2
 
 6520       ptta1=ptxsa1**2+ptysa1**2
 
 6521       IF((eq1**2.LE.pttq1).OR. (eaq1**2.LE.ptta1))            go to 10
 
 6526       IF (ikvala.EQ.1)b33=16.0
 
 6527       es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
 
 6528       hps=
sqrt(es*es+2.*es*0.94)
 
 6530       IF (.NOT.intpt) hps=0.0000001
 
 6535       IF (icoun2.EQ.48)
THEN 
 6538       IF (icoun2.EQ.50)
THEN 
 6546         ptxsq2=qtxsq2+hps*cfe
 
 6547         ptysq2=qtysq2+hps*sfe
 
 6548         ptxsa2=qtxsa2-hps*cfe
 
 6549         ptysa2=qtysa2-hps*sfe
 
 6554       es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
 
 6555       hps=
sqrt(es*es+2.*es*0.94d0)
 
 6558       IF (.NOT.intpt) hps=0.0000001
 
 6562       ptxsq2=qtxsq2+hps*cfe
 
 6563       ptysq2=qtysq2+hps*sfe
 
 6564       ptxsa2=qtxsa2-hps*cfe
 
 6565       ptysa2=qtysa2-hps*sfe
 
 6569       IF (ipev.GE.7)
WRITE(6,1000)ptxsq1,ptysq1,ptxsa1,ptysa1 ,ptxsq2,
 
 6570      +ptysq2,ptxsa2,ptysa2
 
 6572       pttq2=ptxsq2**2+ptysq2**2
 
 6573       ptta2=ptxsa2**2+ptysa2**2
 
 6574       IF((eq2**2.LE.pttq2).OR. (eaq2**2.LE.ptta2))            go to 12
 
 6577         plq1=
sqrt(eq1**2-pttq1)
 
 6578         plaq1=
sqrt(eaq1**2-ptta1)
 
 6579         plq2=-
sqrt(eq2**2-pttq2)
 
 6580         plaq2=-
sqrt(eaq2**2-ptta2)
 
 6584       amch1q=(eq1+eaq2)**2-(ptxsq1+ptxsa2)** 2-(ptysq1+ptysa2)**2-(plq1
 
 6586       IF (amch1q.LE.0.d0)
THEN 
 6588   301   
FORMAT(
' inconsistent Kinematics in SELPT AMCH1Q=',e12.3)
 
 6589        WRITE(6,305) qtxsq1,qtysq1,
 
 6590      +qlq1,qeq1,qtxsa1,qtysa1,qlaq1,qeaq1, qtxsq2,qtysq2,qlq2,qeq2,
 
 6592      +qtysa2,qlaq2,qeaq2, amch1,amch2
 
 6593   305 
FORMAT( 
'PTXSQ1,PTYSQ1, 
 6594      +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2,PTYSQ2,PLQ2,EQ2,PTXSA2, 
 6595      +PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2',5(4e15.5/))
 
 6602       amch2q=(eq2+eaq1)**2-(ptxsq2+ptxsa1)** 2-(ptysq2+ptysa1)**2-(plq2
 
 6604       IF (amch2q.LE.0.d0)
THEN 
 6606   302   
FORMAT(
' inconsistent Kinematics in SELPT AMCH2Q=',e12.3)
 
 6607        WRITE(6,305) qtxsq1,qtysq1,
 
 6608      +qlq1,qeq1,qtxsa1,qtysa1,qlaq1,qeaq1, qtxsq2,qtysq2,qlq2,qeq2,
 
 6610      +qtysa2,qlaq2,qeaq2, amch1,amch2
 
 6621      +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
 
 6622      +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
 
 6623      * pttq2,ptta2, nselpt)
 
 6624       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 6630       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 6632       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
 6634       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
 6635      +ipadis,ishmal,lpauli
 
 6637       COMMON /nncms/  gamcm,bgcm,umo,pcm,eproj,pproj
 
 6638       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 6641       IF(ipev.GE.4)
WRITE(6,6633) ptxsq1,ptysq1,
 
 6642      +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
 
 6643      +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
 
 6645  6633 
FORMAT(
' selpt input: ',
 
 6646      + 
' PTXSQ1,PTYSQ1, +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2, 
 6647      + PTYSQ2,PLQ2,EQ2,PTXSA2, PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2,IREJ, 
 6648      +IKVALA, NSELPT ', 2(8e12.4/),2e12.4,2i5,2e12.4,i5)
 
 6676       IF ( nselpt.EQ.0 .OR.umo.LE.20.d0) 
THEN 
 6682         IF (ikvala.EQ.1)b33=3.0+6.0/log10(umo+10.)
 
 6683         IF (ikvala.EQ.2)b33=4.0+3.0/log10(umo+10.)
 
 6684         es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
 
 6685         hps=
sqrt(es*es+2.*es*0.94)
 
 6688         IF (.NOT.intpt) hps=0.0000001
 
 6692         IF (icount.EQ.48)
THEN 
 6695         IF (icount.EQ.50)
THEN 
 6700         IF (icount.GE.2)
THEN 
 6714         es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
 
 6715         hps=
sqrt(es*es+2.*es*0.94)
 
 6717         IF (.NOT.intpt) hps=0.0000001
 
 6719       ELSEIF(nselpt.EQ.1)
THEN 
 6722           IF(ipev.GE.4)
WRITE(6,6638)hps
 
 6723       ELSEIF(nselpt.EQ.2)
THEN 
 6724         IF (nusept.EQ.0)
THEN 
 6727           IF(ipev.GE.4)
WRITE(6,6638)hps
 
 6730         ELSEIF(nusept.EQ.1)
THEN 
 6736       ptxsq1=qtxsq1+hps*cfe
 
 6737       ptysq1=qtysq1+hps*sfe
 
 6738       ptxsa1=qtxsa1-hps*cfe
 
 6739       ptysa1=qtysa1-hps*sfe
 
 6747       IF (ipev.GE.7)
WRITE(6,1000)ptxsq1,ptysq1,ptxsa1,ptysa1 ,ptxsq2,
 
 6748      +ptysq2,ptxsa2,ptysa2
 
 6749  1000 
FORMAT (
' PT S  ',8f12.6)
 
 6751       pttq1=ptxsq1**2+ptysq1**2
 
 6752       ptta1=ptxsa1**2+ptysa1**2
 
 6754       IF ( nselpt.EQ.0.OR.umo.LE.20.d0 ) 
THEN 
 6757         IF (ikvala.EQ.1)b33=3.0+6.0/log10(umo+10.)
 
 6758         IF (ikvala.EQ.2)b33=4.0+3.0/log10(umo+10.)
 
 6762         IF (icoun2.EQ.48)
THEN 
 6765         IF (icoun2.EQ.50)
THEN 
 6785         es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
 
 6786         hps=
sqrt(es*es+2.*es*0.94)
 
 6789         IF (.NOT.intpt) hps=0.0000001
 
 6791       ELSEIF(nselpt.EQ.1)
THEN 
 6792         IF (musept.EQ.0)
THEN 
 6795           IF(ipev.GE.4)
WRITE(6,6638)hps
 
 6796  6638     
FORMAT (
' SELPT:SAMPPT: HPS= ',e12.4)
 
 6799         ELSEIF(musept.EQ.1)
THEN 
 6805       ptxsq2=qtxsq2+hps*cfe
 
 6806       ptysq2=qtysq2+hps*sfe
 
 6807       ptxsa2=qtxsa2-hps*cfe
 
 6808       ptysa2=qtysa2-hps*sfe
 
 6816       IF (ipev.GE.7)
WRITE(6,1000)ptxsq1,ptysq1,ptxsa1,ptysa1 ,ptxsq2,
 
 6817      +ptysq2,ptxsa2,ptysa2
 
 6819       pttq1=ptxsq1**2+ptysq1**2
 
 6820       ptta1=ptxsa1**2+ptysa1**2
 
 6821       pttq2=ptxsq2**2+ptysq2**2
 
 6822       ptta2=ptxsa2**2+ptysa2**2
 
 6827       IF(plq1.GT.ptwq1.AND.abs(plaq2).GT.ptwq1)
THEN 
 6830       ELSEIF(plq1.GT.ptwa2.AND.abs(plaq2).GT.ptwa2)
THEN 
 6834       IF(plaq1.GT.ptwa1.AND.abs(plq2).GT.ptwa1)
THEN 
 6837       ELSEIF(plaq1.GT.ptwq2.AND.abs(plq2).GT.ptwq2)
THEN 
 6846       ptta1=ptta1+plaq1**2
 
 6848       ptta2=ptta2+plaq2**2
 
 6851       IF(amte1.GE.eq1**2)amte1=eq1**2/2.
 
 6853       IF(amte2.GE.eq2**2)amte2=eq2**2/2.
 
 6855       IF(amte1.GE.eaq1**2)amte1=eaq1**2/2.
 
 6857       IF(amte2.GE.eaq2**2)amte2=eaq2**2/2.
 
 6858         IF((eq1**2-amte1.LE.pttq1).OR.
 
 6859      *  (eq2**2-amte1.LE.pttq2)
 
 6860      *  .OR.(eaq1**2-amte3.LE.ptta1).OR.
 
 6861      *  (eaq2**2-amte4.LE.ptta2))
THEN 
 6862           IF ( nselpt.EQ.0.OR.umo.LE.20.d0 ) 
THEN 
 6866             useptm = useptm * 0.7
 
 6867             IF( usept.GT.0.01d0 .OR. useptm.GT.0.01d0 ) 
THEN 
 6869                 WRITE(6,*)
'  SELPT: JUMP AFTER REDUCTION OF USEPT' 
 6870                 WRITE(6,*)
'  SELPT: USEPT,USEPTM,HPS',usept,useptm,hps
 
 6875       IF(ipev.GE.4)
WRITE(6,6634) ptxsq1,ptysq1,
 
 6876      +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
 
 6877      +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
 
 6879  6634 
FORMAT(
' selpt rejec: ',
 
 6880      + 
' PTXSQ1,PTYSQ1, +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2, 
 6881      + PTYSQ2,PLQ2,EQ2,PTXSA2, PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2,IREJ, 
 6882      +IKVALA, NSELPT ', 2(8e12.4/),2e12.4,3i5)
 
 6891       IF(ip.GE.1)go to 1779
 
 6892         qqq1=qtxsq1**2+qtysq1**2+qlq1**2-pttq1
 
 6893         IF(qqq1.GT.0.d0)
THEN 
 6896           plq1=
sqrt(eq1**2-pttq1)
 
 6898         qqa1=qtxsa1**2+qtysa1**2+qlaq1**2-ptta1
 
 6899         IF(qqa1.GT.0.d0)
THEN 
 6902           plaq1=
sqrt(eaq1**2-ptta1)
 
 6904         qqq2=qtxsq2**2+qtysq2**2+qlq2**2-pttq2
 
 6905         IF(qqq2.GT.0.d0)
THEN 
 6908           plq2=-
sqrt(eq2**2-pttq2)
 
 6910         qqa2=qtxsa2**2+qtysa2**2+qlaq2**2-ptta2
 
 6911         IF(qqa2.GT.0.d0)
THEN 
 6914           plaq2=-
sqrt(eaq2**2-ptta2)
 
 6920       amch1q=(eq1+eaq2)**2-(ptxsq1+ptxsa2)** 2-(ptysq1+ptysa2)**2-(plq1
 
 6922       IF (amch1q.LE.0.d0)
THEN 
 6930   301   
FORMAT(
' inconsistent Kinematics in SELPT AMCH1Q=',e12.3)
 
 6931        WRITE(6,305) qtxsq1,qtysq1,
 
 6932      +qlq1,qeq1,qtxsa1,qtysa1,qlaq1,qeaq1, qtxsq2,qtysq2,qlq2,qeq2,
 
 6934      +qtysa2,qlaq2,qeaq2, amch1,amch2
 
 6935   305 
FORMAT( 
'PTXSQ1,PTYSQ1, 
 6936      +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2,PTYSQ2,PLQ2,EQ2,PTXSA2, 
 6937      +PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2',5(4e15.5/))
 
 6940       IF(ipev.GE.4)
WRITE(6,6635) ptxsq1,ptysq1,
 
 6941      +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
 
 6942      +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
 
 6944  6635 
FORMAT(
' selpt rejec: ',
 
 6945      + 
' PTXSQ1,PTYSQ1, +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2, 
 6946      + PTYSQ2,PLQ2,EQ2,PTXSA2, PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2,IREJ, 
 6947      +IKVALA, NSELPT ', 2(8e12.4/),2e12.4,3i5)
 
 6953       amch2q=(eq2+eaq1)**2-(ptxsq2+ptxsa1)** 2-(ptysq2+ptysa1)**2-(plq2
 
 6955       IF (amch2q.LE.0.d0)
THEN 
 6963   302   
FORMAT(
' inconsistent Kinematics in SELPT AMCH2Q=',e12.3)
 
 6964        WRITE(6,305) qtxsq1,qtysq1,
 
 6965      +qlq1,qeq1,qtxsa1,qtysa1,qlaq1,qeaq1, qtxsq2,qtysq2,qlq2,qeq2,
 
 6967      +qtysa2,qlaq2,qeaq2, amch1,amch2
 
 6970       IF(ipev.GE.4)
WRITE(6,6636) ptxsq1,ptysq1,
 
 6971      +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
 
 6972      +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
 
 6974  6636 
FORMAT(
' selpt rejec: ',
 
 6975      + 
' PTXSQ1,PTYSQ1, +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2, 
 6976      + PTYSQ2,PLQ2,EQ2,PTXSA2, PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2,IREJ, 
 6977      +IKVALA, NSELPT ', 2(8e12.4/),2e12.4,3i5)
 
 6981       IF(ipev.GE.4)
WRITE(6,6637) ptxsq1,ptysq1,
 
 6982      +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
 
 6983      +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
 
 6985  6637 
FORMAT(
' selpt exit : ',
 
 6986      + 
' PTXSQ1,PTYSQ1, +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2, 
 6987      + PTYSQ2,PLQ2,EQ2,PTXSA2, PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2,IREJ, 
 6988      +IKVALA, NSELPT ', 2(8e12.4/),2e12.4,2i5,2e12.4,i5)
 
 6996       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 7002       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 7092       parameter(idmax9=602)
 
 7094       COMMON /ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
 
 7096       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 7108       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 7109      +iibar(210),k1(210),k2(210)
 
 7112       dimension ecmf(3),pcmf(3),codf(3),coff(3),siff(3),itf(3)
 
 7113       dimension ecmff(3),pcmff(3)
 
 7114       dimension cxf(3),cyf(3),czf(3)
 
 7119       IF (iphkk.GE.2) 
WRITE(6,1000) ihkk,nhkk
 
 7120  1000 
FORMAT(
' DECHKK IHKK,NHKK= ',2i5)
 
 7125       IF (ihkk.GT.nhkk)
THEN 
 7130       IF(abs(isthkk(ihkk)).NE.1)       goto 10
 
 7135       IF (it.LT.1.OR.it.GT.210) 
THEN 
 7137  1003   
FORMAT (
' DECHKK IT ',i10)
 
 7148         IF(it.EQ.135.OR.it.EQ.136)                               goto 10
 
 7149         IF(it.GE.1.AND.it.LE.7)                                  goto 10
 
 7150       ELSEIF(istab.EQ.2) 
THEN 
 7151         IF(it.GE. 1.AND.it.LE. 30)                               goto 10
 
 7152         IF(it.GE. 97.AND.it.LE.103)                              goto 10
 
 7153         IF(it.GE.115.AND.it.LE.122)                              goto 10
 
 7154         IF(it.GE.131.AND.it.LE.136)                              goto 10
 
 7155         IF(it.EQ.109)                                            goto 10
 
 7156         IF(it.GE.137.AND.it.LE.160)                              goto 10
 
 7157       ELSEIF(istab.EQ.3) 
THEN 
 7158         IF(it.GE.1.AND.it.LE.23)                                 goto 10
 
 7159         IF(it.GE. 97.AND.it.LE.103)                              goto 10
 
 7160         IF(it.EQ.109.AND.it.EQ.115)                              goto 10
 
 7161         IF(it.GE.133.AND.it.LE.136)                              goto 10
 
 7165       pls=
sqrt(abs(phkk(1,ihkk)**2+phkk(2,ihkk)**2+phkk(3,ihkk)**2))
 
 7169       amtest=
sqrt(abs(phkk(4,ihkk)**2-pls**2))
 
 7170       IF(abs(amtest-phkk(5,ihkk)).GE.1.
d-3)
THEN 
 7173     plss=(phkk(4,ihkk)**2-phkk(5,ihkk))
 
 7174     IF(plss.LE.0.d0)
THEN 
 7175       WRITE(6,
'(A)')
' negative momentum square!' 
 7180     phkk(1,ihkk)=phkk(1,ihkk)*amodp
 
 7181     phkk(2,ihkk)=phkk(2,ihkk)*amodp
 
 7182     phkk(3,ihkk)=phkk(3,ihkk)*amodp
 
 7186       IF(pls.NE.0.d0) 
THEN 
 7187         cxs=phkk(1,ihkk)/pls
 
 7188         cys=phkk(2,ihkk)/pls
 
 7189         czs=phkk(3,ihkk)/pls
 
 7200       IF (vv.GT.wt(iik))                                        go to 20
 
 7208       IF (itf(2).LT.1)                                          go to 10
 
 7212       IF(iphkk.GE.1) 
WRITE(6,1010)it,iik,itf(1),itf(2),itf(3)
 
 7213  1010 
FORMAT(
' DECHKK IT,IIK,IT1,IT2,IT3 ',5i5)
 
 7216       IF(itf(3).EQ.0) 
THEN 
 7218         CALL 
dtwopd(eco,ecmf(1),ecmf(2),pcmf(1),pcmf(2), codf(1),coff
 
 7219      +  (1),siff(1),codf(2),coff(2),siff(2), aam(itf(1)),aam(itf(2)))
 
 7220         sid1=
sqrt(abs((1.-codf(1))*(1.+codf(1))))
 
 7221         sid2=
sqrt(abs((1.-codf(2))*(1.+codf(2))))
 
 7222         pix1=pcmf(1)*sid1*coff(1)
 
 7223         piy1=pcmf(1)*sid1*siff(1)
 
 7224         piz1=pcmf(1)*codf(1)
 
 7225         pix2=pcmf(2)*sid2*coff(2)
 
 7226         piy2=pcmf(2)*sid2*siff(2)
 
 7227         piz2=pcmf(2)*codf(2)
 
 7231         ecm12=ecmf(1)+ecmf(2)-eco
 
 7232         IF((abs(pix12).GT.0.000001d0).OR.
 
 7233      +     (abs(piy12).GT.0.000001d0).OR.
 
 7234      +     (abs(piz12).GT.0.000001d0).OR. 
 
 7235      +     (abs(ecm12).GT.0.000001d0))
THEN 
 7236            WRITE(6,778)pix12,piy12,piz12,ecm12
 
 7237   778      
FORMAT(
' DWOPD px,py,pz,e',4f10.6)
 
 7242        CALL 
dthrep(eco,ecmf(1),ecmf(2),ecmf(3),pcmf(1),pcmf(2),pcmf(3),
 
 7243      +  codf(1),coff(1),siff(1),codf(2),coff(2),siff(2), codf(3),coff
 
 7244      +  (3),siff(3), aam(itf(1)),aam(itf(2)),aam(itf(3)))
 
 7245         sid1=
sqrt((1.-codf(1))*(1.+codf(1)))
 
 7246         sid2=
sqrt((1.-codf(2))*(1.+codf(2)))
 
 7247         sid3=
sqrt((1.-codf(3))*(1.+codf(3)))
 
 7248         pix1=pcmf(1)*sid1*coff(1)
 
 7249         piy1=pcmf(1)*sid1*siff(1)
 
 7250         piz1=pcmf(1)*codf(1)
 
 7251         pix2=pcmf(2)*sid2*coff(2)
 
 7252         piy2=pcmf(2)*sid2*siff(2)
 
 7253         piz2=pcmf(2)*codf(2)
 
 7254         pix3=pcmf(3)*sid3*coff(3)
 
 7255         piy3=pcmf(3)*sid3*siff(3)
 
 7256         piz3=pcmf(3)*codf(3)
 
 7257         pix12=pix1+pix2+pix3
 
 7258         piy12=piy1+piy2+piy3
 
 7259         piz12=piz1+piz2+piz3
 
 7260         ecm12=ecmf(1)+ecmf(2)+ecmf(3)-eco
 
 7261         IF((abs(pix12).GT.0.000001d0).OR.
 
 7262      +     (abs(piy12).GT.0.000001d0).OR.
 
 7263      +     (abs(piz12).GT.0.000001d0).OR. 
 
 7264      +     (abs(ecm12).GT.0.000001d0))
THEN 
 7265            WRITE(6,779)pix12,piy12,piz12,ecm12
 
 7266   779      
FORMAT(
' DTHEPD px,py,pz,e',4f10.6)
 
 7271       jdahkk(1,ihkk)=nhkk + 1
 
 7272       jdahkk(2,ihkk)=nhkk + ndecpr
 
 7274       ehecc=
sqrt(abs(pcmf(id)** 2+ aam(itf(id))**2))
 
 7275         IF (abs(ehecc-ecmf(id)).GT.0.0001d0) 
THEN 
 7276              WRITE(6,
'(2A/3I5,3E15.6)')
 
 7277      &            
' DECHKK: CORRECT INCONSISTENT ENERGY ',
 
 7278      *            
'  IHKK,NHKK,ITF(ID), ECMF(ID),EHECC, AAM(ITF(ID))',
 
 7279      *            ihkk,nhkk,itf(id), ecmf(id),ehecc, aam(itf(id))
 
 7283        CALL 
dtrafo(gam,bgam,cxs,cys,czs, codf(id),coff(id),
 
 7284      * siff(id),pcmf(id),ecmf(id), pcmff(id),cxf(id),
 
 7285      *cyf(id),czf(id),ecmff(id))
 
 7286         IF (iphkk.GE.2) 
WRITE(6,
'(A,7E15.5/8E15.5)')
' DTRAFO ',
 
 7287      * gam,bgam,cxs,cys,czs, codf(id),coff(id),
 
 7288      * siff(id),pcmf(id),ecmf(id), pcmff(id),cxf(id),
 
 7289      *cyf(id),czf(id),ecmff(id)
 
 7296           WRITE (6,1020)nhkk,
nmxhkk 
 7297  1020 
FORMAT (.GT.
' NHKKNMXHKK IN DECHKK RETURN ',2i10)
 
 7304           WRITE (6,
'(A,2I5)') .EQ.
' DECHKK: NHKKNMXHKK ',nhkk,
nmxhkk 
 7310         idhkk(nhkk)=
mpdgha(itf(id))
 
 7315         phkk(1,nhkk)=cxf(id)*pcmff(id)
 
 7316         phkk(2,nhkk)=cyf(id)*pcmff(id)
 
 7317         phkk(3,nhkk)=czf(id)*pcmff(id)
 
 7318       ehecc=
sqrt(abs(pcmff(id)** 2+ aam(itf(id))**2))
 
 7319         IF (abs(ehecc-ecmff(id)).GT.0.003d0) 
THEN 
 7320              WRITE(6,
'(2A/3I5,3E15.6)')
 
 7321      &            
' DECHKK: CORRECT INCONSISTENT ENERGY ',
 
 7322      *            
'  IHKK,NHKK,ITF(ID), ECMF(ID),EHECC, AAM(ITF(ID))',
 
 7323      *            ihkk,nhkk,itf(id), ecmff(id),ehecc, aam(itf(id))
 
 7326         phkk(4,nhkk)=ecmff(id)
 
 7327       phkk(5,nhkk)=aam(itf(id))
 
 7328         vhkk(1,nhkk)=vhkk(1,ihkk)
 
 7329         vhkk(2,nhkk)=vhkk(2,ihkk)
 
 7330         vhkk(3,nhkk)=vhkk(3,ihkk)
 
 7331         vhkk(4,nhkk)=vhkk(4,ihkk)
 
 7333         IF (iphkk.GE.7) 
WRITE(6,1030)nhkk, isthkk(nhkk),idhkk(nhkk),
 
 7334      +  jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
 
 7335      +  (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
 
 7337  1030 
FORMAT (i6,i4,5i6,9e10.2)
 
 7349       SUBROUTINE dtrafo(GAM,BGAM,CX,CY,CZ,COD,COF,SIF,P,ECM,                     
 
 7351       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 7354       sid=
sqrt(1.d0-cod*cod)                                                    
 
 7357       pppt=
sqrt(plx**2+ply**2)
 
 7359       plz=gam*pcmz+bgam*ecm                                                     
 
 7360       pl=
sqrt(plx*plx+ply*ply+plz*plz)                                          
 
 7361       el=gam*ecm+bgam*pcmz                                                      
 
 7367       CALL 
sttran(cx,cy,cz,coz,siz,sif,cof,cxl,cyl,czl)                         
 
 7372       SUBROUTINE sttran(XO,YO,ZO,CDE,SDE,SFE,CFE,X,Y,Z)                         
 
 7373       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 7403       IF ( 
a .LT. anglsq ) 
THEN                                                  
 7413          x=-yo*xi/
a-zo*xo*yi/
a+xo*zi                                            
 
 7414          y=xo*xi/
a-zo*yo*yi/
a+yo*zi                                             
 
 7424       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 7431       parameter(nfimax=249)
 
 7432       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 7433      +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 7434       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 7438       parameter(tiny=1.
d-10)
 
 7462         CALL 
daltra(gam,bgx,bgy,bgz,pxi,pyi,pzi,eei, ppa,pxf(i),pyf(i),
 
 7471       CALL 
daltra(gam,bgx,bgy,bgz,pxsm,pysm,pzsm,esum, ppa,pxsm,pysm,
 
 7480       diffl=pxdif+pydif+pzdif+edif
 
 7481       IF(esum.LT.tiny)esum=tiny
 
 7483       IF(diffl.GE.1.
d-4)
WRITE(6,1000)num,pxdif,pydif,pzdif,edif,pxsm,
 
 7484      +pxsc, pysm,pysc,pzsm,pzsc,esum,esmc
 
 7485  1000 
FORMAT(
' ',2
x,
'LORTRA:NUM=',i5,2
x,
'PXDIF=',1pe15.6,2
x,
'PYDIF=', 1
 
 7486      +pe15.6,2
x,
'PZDIF=',1pe15.6,2
x,
'EDIF=',1pe15.6/2
x,
'PXSM=',1pe15.6,2
 
 7487      +
x,
'PXSC=',1pe15.6,2
x,
'PYSM=',1pe15.6,2
x,
'PYSC=',1pe15.6/2
x,
'PZSM',
 
 7488      +1pe15.6,2
x,
'PZSC=',1pe15.6,2
x,
'ESUM=',1pe15.6,2
x,
'ESMC=',1pe15.6/2
 
 7489      +
x,
'LORTRA DIFFERENCES DUE TO ALTRA'/)
 
 7500       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 7510       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 7517       COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
 
 7518      +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
 
 7519      +prebin,taebin,fermod,etacou
 
 7531       COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
 
 7532      +iibar(210),k1(210),k2(210)
 
 7538       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 7540      *                ,xpsu(248),xtsu(248)
 
 7541      *                ,xpsut(248),xtsut(248)
 
 7542       COMMON /intnez/ ndz,nzd
 
 7544       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 7545      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 7546      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 7548      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 7562       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 7568      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 7570       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 7571       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 7578       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 7581      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 7590       COMMON /nncms/  gamcm,bgcm,umo,pcm,eproj,pproj
 
 7592       COMMON /abrsv/ amcsv1(248),amcsv2(248),gacsv1(248),gacsv2(248),
 
 7593      +bgxsv1(248),bgysv1(248),bgzsv1(248), bgxsv2(248),bgysv2(248),
 
 7594      +bgzsv2(248), nchsv1(248),nchsv2(248),ijcsv1(248),ijcsv2(248),
 
 7595      +pqsva1(248,4),pqsva2(248,4), pqsvb1(248,4),pqsvb2(248,4)
 
 7597       COMMON /abrvs/ amcvs1(248),amcvs2(248),gacvs1(248),gacvs2(248),
 
 7598      +bgxvs1(248),bgyvs1(248),bgzvs1(248), bgxvs2(248),bgyvs2(248),
 
 7599      +bgzvs2(248), nchvs1(248),nchvs2(248),ijcvs1(248),ijcvs2(248),
 
 7600      +pqvsa1(248,4),pqvsa2(248,4), pqvsb1(248,4),pqvsb2(248,4)
 
 7602       COMMON /abrvv/ amcvv1(248),amcvv2(248),gacvv1(248),gacvv2(248),
 
 7603      +bgxvv1(248),bgyvv1(248),bgzvv1(248), bgxvv2(248),bgyvv2(248),
 
 7604      +bgzvv2(248), nchvv1(248),nchvv2(248),ijcvv1(248),ijcvv2(248),
 
 7605      +pqvva1(248,4),pqvva2(248,4), pqvvb1(248,4),pqvvb2(248,4)
 
 7607       COMMON /abrdv/ amcdv1(248),amcdv2(248),gacdv1(248),gacdv2(248),
 
 7608      +bgxdv1(248),bgydv1(248),bgzdv1(248), bgxdv2(248),bgydv2(248),
 
 7609      +bgzdv2(248), nchdv1(248),nchdv2(248),ijcdv1(248),ijcdv2(248),
 
 7610      +pqdva1(248,4),pqdva2(248,4), pqdvb1(248,4),pqdvb2(248,4)
 
 7613       COMMON /abrvd/ amcvd1(248),amcvd2(248),gacvd1(248),gacvd2(248),
 
 7614      +bgxvd1(248),bgyvd1(248),bgzvd1(248), bgxvd2(248),bgyvd2(248),
 
 7615      +bgzvd2(248), nchvd1(248),nchvd2(248),ijcvd1(248),ijcvd2(248),
 
 7616      +pqvda1(248,4),pqvda2(248,4), pqvdb1(248,4),pqvdb2(248,4)
 
 7618       COMMON /abrds/ amcds1(248),amcds2(248),gacds1(248),gacds2(248),
 
 7619      +bgxds1(248),bgyds1(248),bgzds1(248), bgxds2(248),bgyds2(248),
 
 7620      +bgzds2(248), nchds1(248),nchds2(248),ijcds1(248),ijcds2(248),
 
 7621      +pqdsa1(248,4),pqdsa2(248,4), pqdsb1(248,4),pqdsb2(248,4)
 
 7637       COMMON /abrsd/ amcsd1(248),amcsd2(248),gacsd1(248),gacsd2(248),
 
 7638      +bgxsd1(248),bgysd1(248),bgzsd1(248), bgxsd2(248),bgysd2(248),
 
 7639      +bgzsd2(248), nchsd1(248),nchsd2(248),ijcsd1(248),ijcsd2(248),
 
 7640      +pqsda1(248,4),pqsda2(248,4), pqsdb1(248,4),pqsdb2(248,4)
 
 7656       LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
 
 7658       COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
 
 7659      +ipadis,ishmal,lpauli
 
 7661       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 7663       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 7680       COMMON /nucjtn/nonuj1,nonujt,nonus1,nonust
 
 7765             IF(abs(nchss1(
n)).NE.99) 
THEN 
 7766               pxss=pxss + pqssa1(
n,1) + pqssa2(
n,1)
 
 7767               pyss=pyss + pqssa1(
n,2) + pqssa2(
n,2)
 
 7768               pzss=pzss + pqssa1(
n,3) + pqssa2(
n,3)
 
 7769               pess=pess + pqssa1(
n,4) + pqssa2(
n,4)
 
 7771             IF(abs(nchss2(
n)).NE.99) 
THEN 
 7772               pxss=pxss + pqssb1(
n,1) + pqssb2(
n,1)
 
 7773               pyss=pyss + pqssb1(
n,2) + pqssb2(
n,2)
 
 7774               pzss=pzss + pqssb1(
n,3) + pqssb2(
n,3)
 
 7775               pess=pess + pqssb1(
n,4) + pqssb2(
n,4)
 
 7779         pzbss=gamcm*pzss + bgcm*pess
 
 7780         pebss=gamcm*pess + bgcm*pzss
 
 7788           IF(abs(nchsv1(
n)).NE.99) 
THEN 
 7789             pxsv=pxsv +pqsva1(
n,1)+pqsva2(
n,1)
 
 7790             pysv=pysv +pqsva1(
n,2)+pqsva2(
n,2)
 
 7791             pzsv=pzsv +pqsva1(
n,3)+pqsva2(
n,3)
 
 7792             pesv=pesv +pqsva1(
n,4)+pqsva2(
n,4)
 
 7795      +  pxsv,pysv,pzsv,pesv
 
 7801           IF(abs(nchsv2(
n)).NE.99) 
THEN 
 7802             pxsv=pxsv + pqsvb1(
n,1)+pqsvb2(
n,1)
 
 7803             pysv=pysv + pqsvb1(
n,2)+pqsvb2(
n,2)
 
 7804             pzsv=pzsv + pqsvb1(
n,3)+pqsvb2(
n,3)
 
 7805             pesv=pesv + pqsvb1(
n,4)+pqsvb2(
n,4)
 
 7808      +  pxsv,pysv,pzsv,pesv
 
 7812         pzbsv=gamcm*pzsv + bgcm*pesv
 
 7813         pebsv=gamcm*pesv + bgcm*pzsv
 
 7816      +  pxsv,pysv,pzbsv,pebsv
 
 7825           IF(abs(nchvs1(
n)).NE.99) 
THEN 
 7826             pxvs=pxvs + pqvsa1(
n,1) + pqvsa2(
n,1)
 
 7827             pyvs=pyvs + pqvsa1(
n,2) + pqvsa2(
n,2)
 
 7828             pzvs=pzvs + pqvsa1(
n,3) + pqvsa2(
n,3)
 
 7829             pevs=pevs + pqvsa1(
n,4) + pqvsa2(
n,4)
 
 7831           IF(abs(nchvs2(
n)).NE.99) 
THEN 
 7832             pxvs=pxvs + pqvsb1(
n,1) + pqvsb2(
n,1)
 
 7833             pyvs=pyvs + pqvsb1(
n,2) + pqvsb2(
n,2)
 
 7834             pzvs=pzvs + pqvsb1(
n,3) + pqvsb2(
n,3)
 
 7835             pevs=pevs + pqvsb1(
n,4) + pqvsb2(
n,4)
 
 7838         pzbvs=gamcm*pzvs + bgcm*pevs
 
 7839         pebvs=gamcm*pevs + bgcm*pzvs
 
 7845           IF(abs(nchds1(
n)).NE.99) 
THEN 
 7846             pxds=pxds + pqdsa1(
n,1) + pqdsa2(
n,1)
 
 7847             pyds=pyds + pqdsa1(
n,2) + pqdsa2(
n,2)
 
 7848             pzds=pzds + pqdsa1(
n,3) + pqdsa2(
n,3)
 
 7849             peds=peds + pqdsa1(
n,4) + pqdsa2(
n,4)
 
 7851           IF(abs(nchds2(
n)).NE.99) 
THEN 
 7852             pxds=pxds + pqdsb1(
n,1) + pqdsb2(
n,1)
 
 7853             pyds=pyds + pqdsb1(
n,2) + pqdsb2(
n,2)
 
 7854             pzds=pzds + pqdsb1(
n,3) + pqdsb2(
n,3)
 
 7855             peds=peds + pqdsb1(
n,4) + pqdsb2(
n,4)
 
 7858         pzbds=gamcm*pzds + bgcm*peds
 
 7859         pebds=gamcm*peds + bgcm*pzds
 
 7865           IF(abs(nchdz1(
n)).NE.99) 
THEN 
 7866             pxdz=pxdz + pqdza1(
n,1) + pqdza2(
n,1)
 
 7867             pydz=pydz + pqdza1(
n,2) + pqdza2(
n,2)
 
 7868             pzdz=pzdz + pqdza1(
n,3) + pqdza2(
n,3)
 
 7869             pedz=pedz + pqdza1(
n,4) + pqdza2(
n,4)
 
 7871           IF(abs(nchdz2(
n)).NE.99) 
THEN 
 7872             pxdz=pxdz + pqdzb1(
n,1) + pqdzb2(
n,1)
 
 7873             pydz=pydz + pqdzb1(
n,2) + pqdzb2(
n,2)
 
 7874             pzdz=pzdz + pqdzb1(
n,3) + pqdzb2(
n,3)
 
 7875             pedz=pedz + pqdzb1(
n,4) + pqdzb2(
n,4)
 
 7878         pzbdz=gamcm*pzdz + bgcm*pedz
 
 7879         pebdz=gamcm*pedz + bgcm*pzdz
 
 7885           IF(abs(nchsd1(
n)).NE.99) 
THEN 
 7886             pxsd=pxsd + pqsda1(
n,1) + pqsda2(
n,1)
 
 7887             pysd=pysd + pqsda1(
n,2) + pqsda2(
n,2)
 
 7888             pzsd=pzsd + pqsda1(
n,3) + pqsda2(
n,3)
 
 7889             pesd=pesd + pqsda1(
n,4) + pqsda2(
n,4)
 
 7891           IF(abs(nchsd2(
n)).NE.99) 
THEN 
 7892             pxsd=pxsd + pqsdb1(
n,1) + pqsdb2(
n,1)
 
 7893             pysd=pysd + pqsdb1(
n,2) + pqsdb2(
n,2)
 
 7894             pzsd=pzsd + pqsdb1(
n,3) + pqsdb2(
n,3)
 
 7895             pesd=pesd + pqsdb1(
n,4) + pqsdb2(
n,4)
 
 7898         pzbsd=gamcm*pzsd + bgcm*pesd
 
 7899         pebsd=gamcm*pesd + bgcm*pzsd
 
 7905           IF(abs(nchzd1(
n)).NE.99) 
THEN 
 7906             pxzd=pxzd + pqzda1(
n,1) + pqzda2(
n,1)
 
 7907             pyzd=pyzd + pqzda1(
n,2) + pqzda2(
n,2)
 
 7908             pzzd=pzzd + pqzda1(
n,3) + pqzda2(
n,3)
 
 7909             pezd=pezd + pqzda1(
n,4) + pqzda2(
n,4)
 
 7911           IF(abs(nchzd2(
n)).NE.99) 
THEN 
 7912             pxzd=pxzd + pqzdb1(
n,1) + pqzdb2(
n,1)
 
 7913             pyzd=pyzd + pqzdb1(
n,2) + pqzdb2(
n,2)
 
 7914             pzzd=pzzd + pqzdb1(
n,3) + pqzdb2(
n,3)
 
 7915             pezd=pezd + pqzdb1(
n,4) + pqzdb2(
n,4)
 
 7918         pzbzd=gamcm*pzzd + bgcm*pezd
 
 7919         pebzd=gamcm*pezd + bgcm*pzzd
 
 7925           IF(abs(nchdv1(
n)).NE.99) 
THEN 
 7926             pxdv=pxdv + pqdva1(
n,1) + pqdva2(
n,1)
 
 7927             pydv=pydv + pqdva1(
n,2) + pqdva2(
n,2)
 
 7928             pzdv=pzdv + pqdva1(
n,3) + pqdva2(
n,3)
 
 7929             pedv=pedv + pqdva1(
n,4) + pqdva2(
n,4)
 
 7931           IF(abs(nchdv2(
n)).NE.99) 
THEN 
 7932             pxdv=pxdv + pqdvb1(
n,1) + pqdvb2(
n,1)
 
 7933             pydv=pydv + pqdvb1(
n,2) + pqdvb2(
n,2)
 
 7934             pzdv=pzdv + pqdvb1(
n,3) + pqdvb2(
n,3)
 
 7935             pedv=pedv + pqdvb1(
n,4) + pqdvb2(
n,4)
 
 7938         pzbdv=gamcm*pzdv + bgcm*pedv
 
 7939         pebdv=gamcm*pedv + bgcm*pzdv
 
 7945           IF(abs(nchvd1(
n)).NE.99) 
THEN 
 7946             pxvd=pxvd + pqvda1(
n,1) + pqvda2(
n,1)
 
 7947             pyvd=pyvd + pqvda1(
n,2) + pqvda2(
n,2)
 
 7948             pzvd=pzvd + pqvda1(
n,3) + pqvda2(
n,3)
 
 7949             pevd=pevd + pqvda1(
n,4) + pqvda2(
n,4)
 
 7951           IF(abs(nchvd2(
n)).NE.99) 
THEN 
 7952             pxvd=pxvd + pqvdb1(
n,1) + pqvdb2(
n,1)
 
 7953             pyvd=pyvd + pqvdb1(
n,2) + pqvdb2(
n,2)
 
 7954             pzvd=pzvd + pqvdb1(
n,3) + pqvdb2(
n,3)
 
 7955             pevd=pevd + pqvdb1(
n,4) + pqvdb2(
n,4)
 
 7958         pzbvd=gamcm*pzvd + bgcm*pevd
 
 7959         pebvd=gamcm*pevd + bgcm*pzvd
 
 7967           IF((nchvv1(
n).NE.99).AND.(nchvv2(
n).NE.99)) 
THEN 
 7968           pxvv=pxvv+pqvva1(
n,1)+pqvva2(
n,1)+pqvvb1(
n,1)+pqvvb2(
n,1)
 
 7969           pyvv=pyvv+pqvva1(
n,2)+pqvva2(
n,2)+pqvvb1(
n,2)+pqvvb2(
n,2)
 
 7970           pzvv=pzvv+pqvva1(
n,3)+pqvva2(
n,3)+pqvvb1(
n,3)+pqvvb2(
n,3)
 
 7971           pevv=pevv+pqvva1(
n,4)+pqvva2(
n,4)+pqvvb1(
n,4)+pqvvb2(
n,4)
 
 7974         pzbvv=gamcm*pzvv + bgcm*pevv
 
 7975         pebvv=gamcm*pevv + bgcm*pzvv
 
 7998             IF(abs(nchzz1(
n)).NE.99.AND.jhkksx(
n).EQ.1) 
THEN 
 7999             IF(abs(nchzz1(
n)).NE.88) 
THEN 
 8000               pxzz=pxzz + pqzza1(
n,1) + pqzza2(
n,1)
 
 8001               pyzz=pyzz + pqzza1(
n,2) + pqzza2(
n,2)
 
 8002               pzzz=pzzz + pqzza1(
n,3) + pqzza2(
n,3)
 
 8003               pezz=pezz + pqzza1(
n,4) + pqzza2(
n,4)
 
 8006             IF(abs(nchzz2(
n)).NE.99.AND.jhkksx(
n).EQ.1) 
THEN 
 8007             IF(abs(nchzz2(
n)).NE.88) 
THEN 
 8008               pxzz=pxzz + pqzzb1(
n,1) + pqzzb2(
n,1)
 
 8009               pyzz=pyzz + pqzzb1(
n,2) + pqzzb2(
n,2)
 
 8010               pzzz=pzzz + pqzzb1(
n,3) + pqzzb2(
n,3)
 
 8011               pezz=pezz + pqzzb1(
n,4) + pqzzb2(
n,4)
 
 8015         pzbzz=gamcm*pzzz + bgcm*pezz
 
 8016         pebzz=gamcm*pezz + bgcm*pzzz
 
 8023             IF(abs(nchhh1(
n)).NE.99.AND.jhkkex(
n).EQ.1) 
THEN 
 8024               pxhh=pxhh + pqhha1(
n,1) + pqhha2(
n,1)
 
 8025               pyhh=pyhh + pqhha1(
n,2) + pqhha2(
n,2)
 
 8026               pzhh=pzhh + pqhha1(
n,3) + pqhha2(
n,3)
 
 8027               pehh=pehh + pqhha1(
n,4) + pqhha2(
n,4)
 
 8029             IF(abs(nchhh2(
n)).NE.99.AND.jhkkex(
n).EQ.1) 
THEN 
 8030               pxhh=pxhh + pqhhb1(
n,1) + pqhhb2(
n,1)
 
 8031               pyhh=pyhh + pqhhb1(
n,2) + pqhhb2(
n,2)
 
 8032               pzhh=pzhh + pqhhb1(
n,3) + pqhhb2(
n,3)
 
 8033               pehh=pehh + pqhhb1(
n,4) + pqhhb2(
n,4)
 
 8036         pzbhh=gamcm*pzhh + bgcm*pehh
 
 8037         pebhh=gamcm*pehh + bgcm*pzhh
 
 8047         IF(isthkk(i).EQ.11)e0000=e0000+prmom(4,i)
 
 8048         IF(isthkk(i).EQ.11)p0000=p0000+prmom(3,i)
 
 8053         IF(isthkk(i).EQ.12)e0000=e0000+tamom(4,ii)
 
 8054         IF(isthkk(i).EQ.12)p0000=p0000+tamom(3,ii)
 
 8056       p000=gamcm*p0000+bgcm*e0000
 
 8057       e000=gamcm*e0000+bgcm*p0000
 
 8058       iprojo=(pzbal*1.001)/pproj
 
 8059       residu=abs(e000-pebal)/(e000)
 
 8061     WRITE(6,
'(A,2E15.5)')
' E000,PEBAL', e000,pebal
 
 8062         WRITE(6,1000)pxbal,pybal,pzbal,pebal, pxss,pyss,pzbss,pebss,
 
 8063      +  pxsv,pysv,pzbsv,pebsv, pxvs,pyvs,pzbvs,pebvs, pxvv,pyvv,pzbvv,
 
 8064      +  pebvv,pxcc,pycc,pzbcc,pebcc,
 
 8065      +  pxzz,pyzz,pzbzz,pebzz,
 
 8066      +  pxhh,pyhh,pzbhh,pebhh,
 
 8067      +  pxds,pyds,pzbds,pebds,
 
 8068      +  pxsd,pysd,pzbsd,pebsd,
 
 8069      +  pxdz,pydz,pzbdz,pebdz,
 
 8070      +  pxzd,pyzd,pzbzd,pebzd,
 
 8071      +  pxdv,pydv,pzbdv,pebdv,
 
 8072      +  pxvd,pyvd,pzbvd,pebvd
 
 8074       IF (residu.GT.0.02d0)
THEN 
 8077       IF (residu.GT.0.02d0.AND.iphkk.GE.2)
THEN 
 8079     WRITE(6,
'(A,2E15.5)')
' E000,PEBAL', e000,pebal
 
 8080         WRITE(6,1000)pxbal,pybal,pzbal,pebal, pxss,pyss,pzbss,pebss,
 
 8081      +  pxsv,pysv,pzbsv,pebsv, pxvs,pyvs,pzbvs,pebvs, pxvv,pyvv,pzbvv,
 
 8082      +  pebvv,pxcc,pycc,pzbcc,pebcc,
 
 8083      +  pxzz,pyzz,pzbzz,pebzz,
 
 8084      +  pxhh,pyhh,pzbhh,pebhh,
 
 8085      +  pxds,pyds,pzbds,pebds,
 
 8086      +  pxsd,pysd,pzbsd,pebsd,
 
 8087      +  pxdz,pydz,pzbdz,pebdz,
 
 8088      +  pxzd,pyzd,pzbzd,pebzd,
 
 8089      +  pxdv,pydv,pzbdv,pebdv,
 
 8090      +  pxvd,pyvd,pzbvd,pebvd
 
 8092  1000 
FORMAT (
' 4 MOMENTUM CONS.IN EVENT LEVEL OF PARTONS',/ 
' ALL',4e15
 
 8093      +.5/,
' SS ',4e15.5/,
' SV ',4e15.5/ 
' VS ',4e15.5/,
' VV ',4e15.5/,
 
 8156           IF(abs(nchss1(
n)).NE.99) 
THEN 
 8157             pxss=pxss+bgxss1(
n)*amcss1(
n)
 
 8158             pyss=pyss+bgyss1(
n)*amcss1(
n)
 
 8159             pzss=pzss+bgzss1(
n)*amcss1(
n)
 
 8160             pess=pess+gacss1(
n)*amcss1(
n)
 
 8162           IF(abs(nchss2(
n)).NE.99) 
THEN 
 8163             pxss=pxss+bgxss2(
n)*amcss2(
n)
 
 8164             pyss=pyss+bgyss2(
n)*amcss2(
n)
 
 8165             pzss=pzss+bgzss2(
n)*amcss2(
n)
 
 8166             pess=pess+gacss2(
n)*amcss2(
n)
 
 8170         pzbss=gamcm*pzss + bgcm*pess
 
 8171         pebss=gamcm*pess + bgcm*pzss
 
 8191           IF(abs(nchsv1(
n)).NE.99) 
THEN 
 8192           pxsv=pxsv+bgxsv1(
n)*amcsv1(
n)
 
 8193           pysv=pysv+bgysv1(
n)*amcsv1(
n)
 
 8194           pzsv=pzsv+bgzsv1(
n)*amcsv1(
n)
 
 8195           pesv=pesv+gacsv1(
n)*amcsv1(
n)
 
 8197           IF(abs(nchsv2(
n)).NE.99) 
THEN 
 8198           pxsv=pxsv+bgxsv2(
n)*amcsv2(
n)
 
 8199           pysv=pysv+bgysv2(
n)*amcsv2(
n)
 
 8200           pzsv=pzsv+bgzsv2(
n)*amcsv2(
n)
 
 8201           pesv=pesv+gacsv2(
n)*amcsv2(
n)
 
 8204         pzbsv=gamcm*pzsv + bgcm*pesv
 
 8205         pebsv=gamcm*pesv + bgcm*pzsv
 
 8211           IF(abs(nchds1(
n)).NE.99) 
THEN 
 8212           pxds=pxds+bgxds1(
n)*amcds1(
n)
 
 8213           pyds=pyds+bgyds1(
n)*amcds1(
n)
 
 8214           pzds=pzds+bgzds1(
n)*amcds1(
n)
 
 8215           peds=peds+gacds1(
n)*amcds1(
n)
 
 8217           IF(abs(nchds2(
n)).NE.99) 
THEN 
 8218           pxds=pxds+bgxds2(
n)*amcds2(
n)
 
 8219           pyds=pyds+bgyds2(
n)*amcds2(
n)
 
 8220           pzds=pzds+bgzds2(
n)*amcds2(
n)
 
 8221           peds=peds+gacds2(
n)*amcds2(
n)
 
 8224         pzbds=gamcm*pzds + bgcm*peds
 
 8225         pebds=gamcm*peds + bgcm*pzds
 
 8231           IF(abs(nchsd1(
n)).NE.99) 
THEN 
 8232           pxsd=pxsd+bgxsd1(
n)*amcsd1(
n)
 
 8233           pysd=pysd+bgysd1(
n)*amcsd1(
n)
 
 8234           pzsd=pzsd+bgzsd1(
n)*amcsd1(
n)
 
 8235           pesd=pesd+gacsd1(
n)*amcsd1(
n)
 
 8237           IF(abs(nchsd2(
n)).NE.99) 
THEN 
 8238           pxsd=pxsd+bgxsd2(
n)*amcsd2(
n)
 
 8239           pysd=pysd+bgysd2(
n)*amcsd2(
n)
 
 8240           pzsd=pzsd+bgzsd2(
n)*amcsd2(
n)
 
 8241           pesd=pesd+gacsd2(
n)*amcsd2(
n)
 
 8244         pzbsd=gamcm*pzsd + bgcm*pesd
 
 8245         pebsd=gamcm*pesd + bgcm*pzsd
 
 8251           IF(abs(nchdv1(
n)).NE.99) 
THEN 
 8252           pxdv=pxdv+bgxdv1(
n)*amcdv1(
n)
 
 8253           pydv=pydv+bgydv1(
n)*amcdv1(
n)
 
 8254           pzdv=pzdv+bgzdv1(
n)*amcdv1(
n)
 
 8255           pedv=pedv+gacdv1(
n)*amcdv1(
n)
 
 8257           IF(abs(nchdv2(
n)).NE.99) 
THEN 
 8258           pxdv=pxdv+bgxdv2(
n)*amcdv2(
n)
 
 8259           pydv=pydv+bgydv2(
n)*amcdv2(
n)
 
 8260           pzdv=pzdv+bgzdv2(
n)*amcdv2(
n)
 
 8261           pedv=pedv+gacdv2(
n)*amcdv2(
n)
 
 8264         pzbdv=gamcm*pzdv + bgcm*pedv
 
 8265         pebdv=gamcm*pedv + bgcm*pzdv
 
 8271           IF(abs(nchvd1(
n)).NE.99) 
THEN 
 8272           pxvd=pxvd+bgxvd1(
n)*amcvd1(
n)
 
 8273           pyvd=pyvd+bgyvd1(
n)*amcvd1(
n)
 
 8274           pzvd=pzvd+bgzvd1(
n)*amcvd1(
n)
 
 8275           pevd=pevd+gacvd1(
n)*amcvd1(
n)
 
 8277           IF(abs(nchvd2(
n)).NE.99) 
THEN 
 8278           pxvd=pxvd+bgxvd2(
n)*amcvd2(
n)
 
 8279           pyvd=pyvd+bgyvd2(
n)*amcvd2(
n)
 
 8280           pzvd=pzvd+bgzvd2(
n)*amcvd2(
n)
 
 8281           pevd=pevd+gacvd2(
n)*amcvd2(
n)
 
 8284         pzbvd=gamcm*pzvd + bgcm*pevd
 
 8285         pebvd=gamcm*pevd + bgcm*pzvd
 
 8293           IF(abs(nchvs1(
n)).NE.99) 
THEN 
 8294           pxvs=pxvs+bgxvs1(
n)*amcvs1(
n)
 
 8295           pyvs=pyvs+bgyvs1(
n)*amcvs1(
n)
 
 8296           pzvs=pzvs+bgzvs1(
n)*amcvs1(
n)
 
 8297           pevs=pevs+gacvs1(
n)*amcvs1(
n)
 
 8299           IF(abs(nchvs2(
n)).NE.99) 
THEN 
 8300           pxvs=pxvs+bgxvs2(
n)*amcvs2(
n)
 
 8301           pyvs=pyvs+bgyvs2(
n)*amcvs2(
n)
 
 8302           pzvs=pzvs+bgzvs2(
n)*amcvs2(
n)
 
 8303           pevs=pevs+gacvs2(
n)*amcvs2(
n)
 
 8306         pzbvs=gamcm*pzvs + bgcm*pevs
 
 8307         pebvs=gamcm*pevs + bgcm*pzvs
 
 8314           IF(abs(nchzz1(
n)).NE.99.AND.jhkksx(
n).EQ.1) 
THEN 
 8315             pxzz=pxzz+bgxzz1(
n)*amczz1(
n)
 
 8316             pyzz=pyzz+bgyzz1(
n)*amczz1(
n)
 
 8317             pzzz=pzzz+bgzzz1(
n)*amczz1(
n)
 
 8318             pezz=pezz+gaczz1(
n)*amczz1(
n)
 
 8320           IF(abs(nchzz2(
n)).NE.99.AND.jhkksx(
n).EQ.1) 
THEN 
 8321             pxzz=pxzz+bgxzz2(
n)*amczz2(
n)
 
 8322             pyzz=pyzz+bgyzz2(
n)*amczz2(
n)
 
 8323             pzzz=pzzz+bgzzz2(
n)*amczz2(
n)
 
 8324             pezz=pezz+gaczz2(
n)*amczz2(
n)
 
 8327         pzbzz=gamcm*pzzz + bgcm*pezz
 
 8328         pebzz=gamcm*pezz + bgcm*pzzz
 
 8334           IF(abs(nchhh1(
n)).NE.99.AND.jhkkex(
n).EQ.1) 
THEN 
 8335             pxhh=pxhh+bgxhh1(
n)*amchh1(
n)
 
 8336             pyhh=pyhh+bgyhh1(
n)*amchh1(
n)
 
 8337             pzhh=pzhh+bgzhh1(
n)*amchh1(
n)
 
 8338             pehh=pehh+gachh1(
n)*amchh1(
n)
 
 8340           IF(abs(nchhh2(
n)).NE.99.AND.jhkkex(
n).EQ.1) 
THEN 
 8341             pxhh=pxhh+bgxhh2(
n)*amchh2(
n)
 
 8342             pyhh=pyhh+bgyhh2(
n)*amchh2(
n)
 
 8343             pzhh=pzhh+bgzhh2(
n)*amchh2(
n)
 
 8344             pehh=pehh+gachh2(
n)*amchh2(
n)
 
 8347         pzbhh=gamcm*pzhh + bgcm*pehh
 
 8348         pebhh=gamcm*pehh + bgcm*pzhh
 
 8355           IF((nchvv1(
n).NE.99).AND.(nchvv2(
n).NE.99)) 
THEN 
 8356           pxvv=pxvv+bgxvv1(
n)*amcvv1(
n)+bgxvv2(
n)*amcvv2(
n)
 
 8357           pyvv=pyvv+bgyvv1(
n)*amcvv1(
n)+bgyvv2(
n)*amcvv2(
n)
 
 8358           pzvv=pzvv+bgzvv1(
n)*amcvv1(
n)+bgzvv2(
n)*amcvv2(
n)
 
 8359           pevv=pevv+gacvv1(
n)*amcvv1(
n)+gacvv2(
n)*amcvv2(
n)
 
 8362         pzbvv=gamcm*pzvv + bgcm*pevv
 
 8363         pebvv=gamcm*pevv + bgcm*pzvv
 
 8370       IF (ipev.GE.1) 
WRITE(6,1010)pxbal,pybal,pzbal,
 
 8371      +pebal, pxss,pyss,pzbss,pebss, pxsv,pysv,pzbsv,pebsv, pxvs,pyvs,
 
 8372      +pzbvs,pebvs, pxvv,pyvv,pzbvv,pebvv, pxcc,pycc,pzbcc,pebcc,
 
 8373      +  pxds,pyds,pzbds,pebds,
 
 8374      +  pxzz,pyzz,pzbzz,pebzz,
 
 8375      +  pxhh,pyhh,pzbhh,pebhh,
 
 8376      +  pxsd,pysd,pzbsd,pebsd,
 
 8377      +  pxdv,pydv,pzbdv,pebdv,
 
 8378      +  pxvd,pyvd,pzbvd,pebvd
 
 8379  1010 
FORMAT (
' 4 MOMENTUM CONS.IN EVENT LEVEL OF CHAINS',/ 
' ALL',4e15.
 
 8380      +5/,
' SS ',4e15.5/,
' SV ',4e15.5/ 
' VS ',4e15.5/,
' VV ',4e15.5/,
 
 8393       SUBROUTINE corval(AMMM,IREJ,AMCH1,AMCH2, QTX1,QTY1,QZ1,QE1,QTX2,
 
 8394      +qty2,qz2,qe2,norig)
 
 8395       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 8404       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 8408       IF(ammm.LE.amch1+amch2+0.4d0) 
THEN 
 8413       ek1=(ammm**2-amch2**2 + amch1**2)/(2.*ammm)
 
 8415       pzk1=
sqrt(ek1**2 - amch1**2)
 
 8417       pzk2=
sqrt(ek2**2 - amch2**2)
 
 8463       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 8484       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 8485      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 8486      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 8488      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 8501       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx),
 
 8502      *                ifrovt(248),itovt(248),ifrost(
intmx),
 
 8503      *             jsshs(
intmx),jtshs(
intmx),jhkknp(248),jhkknt(248),
 
 8509      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 8512       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 8515      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 8517       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 8518       COMMON /lozuo/  zuovp(248),zuosp(
intmx),zuovt(248),
 
 8533       COMMON /hardha/nhard1,nhkkha
 
 8536        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 8537      *                 anndv,annvd,annds,annsd,
 
 8539      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 8541      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 8544      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 8545        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 8546      *                 acouzz,acouhh,acouds,acousd,
 
 8547      *                 acoudz,acouzd,acoudi,
 
 8548      *                 acoudv,acouvd,acoucc
 
 8550       COMMON /pshow/ ipshow
 
 8552       COMMON /harlun/ qlun,iharlu
 
 8553       COMMON /jspar/pxj(1000),pyj(1000),pzj(1000),hej(1000),nnnpj
 
 8554       COMMON /jspa/pxs(40000),pys(40000),pzs(40000),hes(40000),nnnps
 
 8558       parameter(nfimax=249)
 
 8559       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 8560      +hep(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 8561       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 8565       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk),
 
 8569       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 8570       COMMON /projk/ iprojk
 
 8571       COMMON /nucjtn/nonuj1,nonujt,nonus1,nonust
 
 8572       COMMON /gluspl/nugluu,nsgluu
 
 8573       COMMON /nomije/ ptmije(10),nnmije(10)
 
 8576       dimension poj(4),pat(4)
 
 8583           IF (iphkk.GE.2)
WRITE(6,7789)nonujt,ncalhh
 
 8584  7789     
FORMAT (
' HADRHH NONUJT,NCALHH ',2i10)
 
 8585           IF (jhkkex(i).EQ.1)
THEN 
 8587             WRITE (6,7744)i,
intmx 
 8588  7744       
FORMAT (.GT.
'  HADRHH IINTMX  ',2i10)
 
 8600         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 8601         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 8602         CALL 
parpt(2,pt1,pt2,6,nevt)
 
 8606               pojpt=
sqrt(poj(2)**2+poj(1)**2)
 
 8607               patpt=
sqrt(pat(1)**2+pat(2)**2)
 
 8609               IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
 
 8611               IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
 
 8614               qlun=min(pojpt,patpt)
 
 8615               IF((qlun.LT.2.5d0).OR.(amjch1(i).LT.5.d0))
THEN 
 8623             IF (gamjh1(i).LT.0.001d0.OR.amjch1(i).LT.2.d0)
THEN 
 8625      *      i,nhad,amjch1(i),poj,pat,gamjh1(i),bgxjh1(i),
 
 8626      *      bgyjh1(i),bgzjh1(i),ifb1,ifb2,ifb3,ifb4,jhkkex(i)
 
 8627  7788       
FORMAT (
' HADRHH ',2i5,8e12.2/5e12.2,5i5)
 
 8630             CALL 
hadjet(nhad,amjch1(i),poj,pat,gamjh1(i),bgxjh1(i),
 
 8631      *              bgyjh1(i),bgzjh1(i),ifb1,ifb2,ifb3,ifb4,
 
 8637         IF(iphkk.GE.3)
WRITE(6,*)
' HADRHH:NHKK,NHKKAU ',nhkk,nhkkau
 
 8638             IF (nhad.GT.nfimax) 
THEN 
 8639               WRITE (6,7755)nhad,nfimax
 
 8640  7755         
FORMAT (.GT.
' NHADNFIMAX ',2i10)
 
 8643      IF(ndone.EQ.-107801) 
WRITE(6,*)
' First chain HADRHH' 
 8647                 WRITE (*,
'(A,2I5/A)') .EQ.
' HADRHH: NHKKNMXHKK ',
 
 8652               ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 8653               IF (abs(ehecc-hep(j)).GT.0.001d0) 
THEN 
 8662               pthh=
sqrt(pxf(j)**2+pyf(j)**2)+pthh
 
 8665           IF(ibarf(j).EQ.500)istist=2
 
 8667      *                pxf(j),pyf(j),pzf(j),hep(j),nhkkau,iormo(j),9)
 
 8668               IF(idhkk(nhkk).EQ.99999) 
WRITE (6,5009)nhkk,nref(j),
 
 8671               IF (ndone.EQ.-107801) 
WRITE(6,5001) j,nhkk,
 
 8672      *        isthkk(nhkk),idhkk(nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk),
 
 8673      &        jdahkk(1,nhkk),jdahkk(2,nhkk),(phkk(khkk,nhkk),khkk=1,5),
 
 8674      &        (vhkk(khkk,nhkk),khkk=1,4)
 
 8675      IF(iphkk.GE.3) 
WRITE(6,*)
' First chain HADRHH' 
 8676               IF (iphkk.GE.3) 
WRITE(6,5001) nhkk,
 
 8677      *        isthkk(nhkk),idhkk(nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk),
 
 8678      &        jdahkk(1,nhkk),jdahkk(2,nhkk),(phkk(khkk,nhkk),khkk=1,5),
 
 8679      &        (vhkk(khkk,nhkk),khkk=1,4)
 
 8687             DO 137 j=nnnps,nnnpsu
 
 8689               IF(j.GT.40000.OR.jj.GT.1000)
THEN 
 8702             IF (nugluu.EQ.1) go to 5111
 
 8710         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 8711         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 8712         CALL 
parpt(2,pt1,pt2,6,nevt)
 
 8716               pojpt=
sqrt(poj(2)**2+poj(1)**2)
 
 8717               patpt=
sqrt(pat(1)**2+pat(2)**2)
 
 8719               IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
 
 8721               IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
 
 8724               qlun=min(pojpt,patpt)
 
 8725               IF((qlun.LT.2.5d0).OR.(amjch2(i).LT.5.d0))
THEN 
 8733             CALL 
hadjet(nhad,amjch2(i),poj,pat,gamjh2(i),bgxjh2(i),
 
 8734      *                bgyjh2(i),bgzjh2(i),ifb1,ifb2,ifb3,ifb4,
 
 8744               WRITE (*,
'(A,2I5/A)') .EQ.
' HADRHH: NHKKNMXHKK ',
 
 8749             ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 8750             IF (abs(ehecc-hep(j)).GT.0.001d0) 
THEN 
 8759             pthh=
sqrt(pxf(j)**2+pyf(j)**2)+pthh
 
 8762           IF(ibarf(j).EQ.500)istist=2
 
 8764      *                pxf(j),pyf(j),pzf(j),hep(j),nhkkau,iormo(j),10)
 
 8765             IF(idhkk(nhkk).EQ.99999) 
WRITE (6,5009)nhkk,nref(j),
 
 8768             IF (iphkk.GE.7) 
WRITE(6,5001) nhkk,
 
 8769      *      isthkk(nhkk),idhkk(nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk),
 
 8770      &      jdahkk(1,nhkk),jdahkk(2,nhkk),(phkk(khkk,nhkk),khkk=1,5),
 
 8771      &      (vhkk(khkk,nhkk),khkk=1,4)
 
 8779             DO 187 j=nnnps,nnnpsu
 
 8781               IF(j.GT.40000.OR.jj.GT.1000)
THEN 
 8800  5001 
FORMAT (i6,i6,5i6,9e10.2)
 
 8801  5003 
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
 
 8802  5009 
FORMAT (
' NHKK,IDHKK(NHKK)  ',3i10)
 
 8809       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 8838       COMMON /nucjtn/nonuj1,nonujt,nonus1,nonust
 
 8840       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 8841      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 8842      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 8844      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 8857       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx),
 
 8858      *                ifrovt(248),itovt(248),ifrost(
intmx),
 
 8859      *            jsshs(
intmx),jtshs(
intmx),jhkknp(248),jhkknt(248),
 
 8865      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 8868       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 8871      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 8873       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 8874       COMMON /lozuo/  zuovp(248),zuosp(
intmx),zuovt(248),
 
 8883        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 8884      *                 anndv,annvd,annds,annsd,
 
 8886      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 8888      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 8891      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 8892        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 8893      *                 acouzz,acouhh,acouds,acousd,
 
 8894      *                 acoudz,acouzd,acoudi,
 
 8895      *                 acoudv,acouvd,acoucc
 
 8897       COMMON /pshow/ ipshow
 
 8899       COMMON /harlun/ qlun,iharlu
 
 8900       COMMON /jspar/pxj(1000),pyj(1000),pzj(1000),hej(1000),nnnpj
 
 8901       COMMON /jspa/pxs(40000),pys(40000),pzs(40000),hes(40000),nnnps
 
 8905       parameter(nfimax=249)
 
 8906       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 8907      +hep(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 8908       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 8912       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk),
 
 8916       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 8917       COMMON /projk/ iprojk
 
 8918       COMMON /gluspl/nugluu,nsgluu
 
 8919       COMMON /nomije/ ptmije(10),nnmije(10)
 
 8921       dimension poj(4),pat(4)
 
 8925           IF(nch1(i).EQ.99.OR.nch1(i).EQ.88)go to 20
 
 8926           IF(nch2(i).EQ.99.OR.nch2(i).EQ.88)go to 20
 
 8929           IF (iphkk.GE.7)
WRITE(6,7789)nonust,ncalzz,jhkksx(i)
 
 8930  7789     
FORMAT (
' HADRZZ NONUST,NCALZZ,Jhkksx(i) ',3i10)
 
 8931           IF (jhkksx(i).EQ.1)
THEN 
 8933             WRITE (6,7744)i,
intmx 
 8934  7744       
FORMAT (.GT.
'  HADRZZ IINTMX  ',2i10)
 
 8946         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 8947         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 8948         CALL 
parpt(2,pt1,pt2,5,nevt)
 
 8952               pojpt=
sqrt(poj(2)**2+poj(1)**2)
 
 8953               patpt=
sqrt(pat(1)**2+pat(2)**2)
 
 8955               IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
 
 8957               IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
 
 8960               qlun=min(pojpt,patpt)
 
 8961               IF((qlun.LT.2.5d0).OR.(amcch1(i).LT.5.d0))
THEN 
 8969             IF (gamch1(i).LT.0.001d0)
WRITE (6,7788)
 
 8970      *      i,nhad,amcch1(i),poj,pat,gamch1(i),bgxch1(i),
 
 8971      *      bgych1(i),bgzch1(i),ifb1,ifb2,ifb3,ifb4,jhkksx(i)
 
 8972  7788       
FORMAT (
' HADRZZ ',2i5,10e12.2/3e12.2,5i5)
 
 8973             CALL 
hadjet(nhad,amcch1(i),poj,pat,gamch1(i),bgxch1(i),
 
 8974      *              bgych1(i),bgzch1(i),ifb1,ifb2,ifb3,ifb4,
 
 8975      *              ijczz1(i),ijczz1(i),3,nchzz1(i),23)
 
 8980             IF (nhad.GT.nfimax) 
THEN 
 8981               WRITE (6,7755)nhad,nfimax
 
 8982  7755         
FORMAT (.GT.
' NHADNFIMAX ',2i10)
 
 8988                 WRITE (*,
'(A,2I5/A)') .EQ.
' HADRZZ: NHKKNMXHKK ',
 
 8993               ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 8994               IF (abs(ehecc-hep(j)).GT.0.001d0) 
THEN 
 9003               ptzz=
sqrt(pxf(j)**2+pyf(j)**2)+ptzz
 
 9006           IF(ibarf(j).EQ.500)istist=2
 
 9008      *                pxf(j),pyf(j),pzf(j),hep(j),nhkkau,iormo(j),11)
 
 9009               IF(idhkk(nhkk).EQ.99999) 
WRITE (6,5009)nhkk,nref(j),
 
 9012               IF (iphkk.GE.7) 
WRITE(6,5001) nhkk,
 
 9013      *        isthkk(nhkk),idhkk(nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk),
 
 9014      &        jdahkk(1,nhkk),jdahkk(2,nhkk),(phkk(khkk,nhkk),khkk=1,5),
 
 9015      &        (vhkk(khkk,nhkk),khkk=1,4)
 
 9023             DO 137 j=nnnps,nnnpsu
 
 9025               IF(j.GT.40000.OR.jj.GT.1000)
THEN 
 9044         pt1=
sqrt(poj(1)**2+poj(2)**2)
 
 9045         pt2=
sqrt(pat(1)**2+pat(2)**2)
 
 9046         CALL 
parpt(2,pt1,pt2,5,nevt)
 
 9050               pojpt=
sqrt(poj(2)**2+poj(1)**2)
 
 9051               patpt=
sqrt(pat(1)**2+pat(2)**2)
 
 9053               IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
 
 9055               IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
 
 9058               qlun=min(pojpt,patpt)
 
 9059               IF((qlun.LT.2.5d0).OR.(amcch2(i).LT.5.d0))
THEN 
 9067             CALL 
hadjet(nhad,amcch2(i),pat,poj,gamch2(i),bgxch2(i),
 
 9068      *                bgych2(i),bgzch2(i),ifb1,ifb2,ifb3,ifb4,
 
 9069      *              ijczz2(i),ijczz2(i),3,nchzz2(i),24)
 
 9078               WRITE (*,
'(A,2I5/A)') .EQ.
' HADRZZ: NHKKNMXHKK ',
 
 9083             ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
 
 9084             IF (abs(ehecc-hep(j)).GT.0.001d0) 
THEN 
 9093             ptzz=
sqrt(pxf(j)**2+pyf(j)**2)+ptzz
 
 9096           IF(ibarf(j).EQ.500)istist=2
 
 9098      *                pxf(j),pyf(j),pzf(j),hep(j),nhkkau,iormo(j),12)
 
 9099             IF(idhkk(nhkk).EQ.99999) 
WRITE (6,5009)nhkk,nref(j),
 
 9102             IF (iphkk.GE.7) 
WRITE(6,5001) nhkk,
 
 9103      *      isthkk(nhkk),idhkk(nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk),
 
 9104      &      jdahkk(1,nhkk),jdahkk(2,nhkk),(phkk(khkk,nhkk),khkk=1,5),
 
 9105      &      (vhkk(khkk,nhkk),khkk=1,4)
 
 9113             DO 187 j=nnnps,nnnpsu
 
 9115               IF(j.GT.40000.OR.jj.GT.1000)
THEN 
 9132  5001 
FORMAT (i6,i4,5i6,9e10.2)
 
 9133  5003 
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
 
 9134  5009 
FORMAT (
' NHKK,IDHKK(NHKK)  ',3i10)
 
 9146       IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
 
 9150       CHARACTER*8 projty,targty
 
 9153       COMMON /user1/
title,projty,targty
 
 9154       COMMON /user2/cmener,sdfrac,ptlar,istruf,isingd,idubld
 
 9156       c=4.*(0.15
d-24+0.01
d-24*
log(cmener))
 
 9158       IF ((
p). eq .(1.d0)) 
THEN 
 9163       r=dsqrt(-c*dlog(1.d00-
p))
 
 9174       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 9191       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 9193      *                ,xpsu(248),xtsu(248)
 
 9194      *                ,xpsut(248),xtsut(248)
 
 9195        common/popcck/pdbck,pdbse,pdbseu,
 
 9196      *  ijpock,irejck,ick4,ihad4,ick6,ihad6
 
 9197      *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
 
 9198      *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
 
 9199      *isea43,isea63,irejao
 
 9201       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 9202      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 9203      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 9205      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 9219       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 9225      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 9227       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 9228       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 9235       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 9238      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 9243       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 9331       COMMON /abrvs/ amcvs1(248),amcvs2(248),gacvs1(248),gacvs2(248),
 
 9332      +bgxvs1(248),bgyvs1(248),bgzvs1(248), bgxvs2(248),bgyvs2(248),
 
 9333      +bgzvs2(248), nchvs1(248),nchvs2(248),ijcvs1(248),ijcvs2(248),
 
 9334      +pqvsa1(248,4),pqvsa2(248,4), pqvsb1(248,4),pqvsb2(248,4)
 
 9337       parameter(nfimax=249)
 
 9338       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 9339      +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 9340       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 9343       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 9345       COMMON /projk/ iprojk
 
 9347       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 9350        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 9351      *                 anndv,annvd,annds,annsd,
 
 9353      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 9355      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 9358      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 9359        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 9360      *                 acouzz,acouhh,acouds,acousd,
 
 9361      *                 acoudz,acouzd,acoudi,
 
 9362      *                 acoudv,acouvd,acoucc
 
 9364       COMMON /zsea/zseaav,zseasu,anzsea
 
 9365       COMMON /casadi/casaxx,icasad
 
 9370         IF(nchvs1(i).EQ.99.AND.nchvs2(i).EQ.99) go to 50
 
 9374         IF (ipco.GE.6) 
WRITE (6,1010) ipvq(is1),ippv1(is1),ippv2(is1),
 
 9375      +  itsq(is2),itsaq(is2), amcvs1(i),amcvs2(i),gacvs1(i),gacvs2(i),
 
 9376      +  bgxvs1(i),bgyvs1(i),bgzvs1(i), bgxvs2(i),bgyvs2(i),bgzvs2(i),
 
 9377      +  nchvs1(i),nchvs2(i),ijcvs1(i),ijcvs2(i), pqvsa1(i,4),pqvsa2
 
 9378      +  (i,4),pqvsb1(i,4),pqvsb2(i,4)
 
 9391        ippp = ifrovp(intvs1(i))
 
 9397         WRITE(6,*)
' VS qq-q ,IFB1,IFB2,IFB3,',
 
 9398      *  
'INTVS1=IS1,INTVS2=IS2,JIPP,JITTX',
 
 9399      *  ifb1,ifb2,ifb3,intvs1(i),intvs2(i),jipp,jittx
 
 9400         WRITE (6,*)
' target sea quark IFB3=',ifb3,
 
 9401      *  
' from IS2=',intvs2(i)
 
 9402         WRITE(6,*)
' with ITSQ(IS2),XTSQ(IS2),IFROST(IS2)',
 
 9403      *  itsq(is2),xtsq(is2),ifrost(is2)
 
 9406       IF(ifrost(is2).EQ.ifrovt(ii))iii=ii
 
 9409         WRITE (6,*)
' projectile III=',iii
 
 9410         WRITE(6,*)
' corresp. XTVQ(i),XTVD(i),ITVQ(I),ITTV1(I),ITTV2(I)',
 
 9411      *   xtvq(iii),xtvd(iii),itvq(iii),ittv1(iii),ittv2(iii)
 
 9418          IF(
rndm(vv).LE.casaxx)
THEN 
 9419        IF(
rndm(vvv).LE.0.5d0)
THEN 
 9426         WRITE(6,*)
' Cas VS2 qq-q 1 ,IFB1,IFB2,IFB3,',
 
 9427      *  
'INTVS1=IS1,INTVS2=IS2,III',
 
 9428      *  ifb1,ifb2,ifb3,intvs1(i),intvs2(i),iii
 
 9429      *  ,
'-----------------------------------------------------' 
 9438         WRITE(6,*)
' Cas VS2 qq-q 2 ,IFB1,IFB2,IFB3,',
 
 9439      *  
'INTVS1=IS1,INTVS2=IS2,III',
 
 9440      *  ifb1,ifb2,ifb3,intvs1(i),intvs2(i),iii
 
 9441      *  ,
'-----------------------------------------------------' 
 9452  1010 
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
 
 9458       IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
 
 9474       COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
 
 9476      *                ,xpsu(248),xtsu(248)
 
 9477      *                ,xpsut(248),xtsut(248)
 
 9478        common/popcck/pdbck,pdbse,pdbseu,
 
 9479      *  ijpock,irejck,ick4,ihad4,ick6,ihad6
 
 9480      *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
 
 9481      *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
 
 9482      *isea43,isea63,irejao
 
 9484       COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
 
 9485      +ixpv,ixps,ixtv,ixts, intvv1(248),
 
 9486      +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
 
 9488      +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
 
 9502       COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
 
 9508      +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
 
 9510       LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
 
 9511       COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
 
 9518       COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
 
 9521      +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
 
 9526       COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
 
 9614       COMMON /abrsv/ amcsv1(248),amcsv2(248),gacsv1(248),gacsv2(248),
 
 9615      +bgxsv1(248),bgysv1(248),bgzsv1(248), bgxsv2(248),bgysv2(248),
 
 9616      +bgzsv2(248), nchsv1(248),nchsv2(248),ijcsv1(248),ijcsv2(248),
 
 9617      +pqsva1(248,4),pqsva2(248,4), pqsvb1(248,4),pqsvb2(248,4)
 
 9620       parameter(nfimax=249)
 
 9621       COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
 
 9622      +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
 
 9623       COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
 
 9626       COMMON /dprin/  ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
 
 9628       COMMON /projk/ iprojk
 
 9630       COMMON /nucc/   it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
 
 9633        COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
 
 9634      *                 anndv,annvd,annds,annsd,
 
 9636      *                 ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
 
 9638      *                 eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
 
 9641      *                ,annzd,anndz,ptzd,ptdz,eezd,eedz
 
 9642        COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
 
 9643      *                 acouzz,acouhh,acouds,acousd,
 
 9644      *                 acoudz,acouzd,acoudi,
 
 9645      *                 acoudv,acouvd,acoucc
 
 9647       COMMON /zsea/zseaav,zseasu,anzsea
 
 9648       COMMON /casadi/casaxx,icasad
 
 9655         IF(nchsv1(i).EQ.99.AND.nchsv2(i).EQ.99) go to 50
 
 9659         IF (ipco.GE.6) 
WRITE (6,1000) ipsq(is1),ipsaq(is1),itvq(is2),
 
 9660      +  ittv1(is2),ittv2(is2), amcsv1(i),amcsv2(i),gacsv1(i),gacsv2(i),
 
 9661      +  bgxsv1(i),bgysv1(i),bgzsv1(i), bgxsv2(i),bgysv2(i),bgzsv2(i),
 
 9662      +  nchsv1(i),nchsv2(i),ijcsv1(i),ijcsv2(i), pqsva1(i,4),pqsva2
 
 9663      +  (i,4),pqsvb1(i,4),pqsvb2(i,4)
 
 9664  1000 
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
 
 9677        ittt = ifrovt(intsv2(i))
 
 9683         WRITE(6,*)
' SV q-qq ,IFB1,IFB2,IFB3,',
 
 9684      *  
'INTSV1=IS1,INTSV2=IS2,JIPPX,JITT',
 
 9685      *  ifb1,ifb2,ifb3,intsv1(i),intsv2(i),jippx,jitt
 
 9686         WRITE (6,*)
' projectile sea quark IFB1=',ifb1,
 
 9687      *  
' from IS1=',intsv1(i)
 
 9688         WRITE(6,*)
' with IPSQ(IS1),XPSQ(IS1),IFROSP(IS1)',
 
 9689      *  ipsq(is1),xpsq(is1),ifrosp(is1)
 
 9692       IF(ifrosp(is1).EQ.ifrovp(ii))iii=ii
 
 9695         WRITE (6,*)
' projectile III=',iii
 
 9696         WRITE(6,*)
' corresp. XPVQ(i),XPVD(i),IPVQ(I),IPPV1(I),IPPV2(I)',
 
 9697      *   xpvq(iii),xpvd(iii),ipvq(iii),ippv1(iii),ippv2(iii)
 
 9704          IF(
rndm(vv).LE.casaxx)
THEN 
 9705        IF(
rndm(vvv).LE.0.5d0)
THEN 
 9712         WRITE(6,*)
' Cas SV1 q-qq 1 ,IFB1,IFB2,IFB3,',
 
 9713      *  
'INTSV1=IS1,INTSV2=IS2,JIPPX,JITT,III',
 
 9714      *  ifb1,ifb2,ifb3,intsv1(i),intsv2(i),jippx,jitt,iii
 
 9715      *  ,
'-----------------------------------------------------' 
 9724         WRITE(6,*)
' Cas SV1 q-qq 2 ,IFB1,IFB2,IFB3,',
 
 9725      *  
'INTSV1=IS1,INTSV2=IS2,JIPPX,JITT,III',
 
 9726      *  ifb1,ifb2,ifb3,intsv1(i),intsv2(i),jippx,jitt,iii
 
 9727      *  ,
'-----------------------------------------------------' 
subroutine calbam(NNCH, I1, I2, IFB11, IFB22, IFB33, IFB44, AMCH, NOBAM, IHAD)
 
subroutine xksamp(NN, ECM)
 
subroutine hadjse(NHAD, AMCH, PPR, PTA, GAM, BGX, BGY, BGZ, IFB1, IFB2, IFB3, IFB4, I1, I2, NOBAM, NNCH, NORIG, IREJ, IISSQQ)
 
subroutine flksaa(NN, ECM)
 
subroutine hadrkk(NHKKH1, PPN)
 
subroutine dtrans(XO, YO, ZO, CDE, SDE, SFE, CFE, X, Y, Z)
 
subroutine sltraf(GA, BGA, EIN, PZIN, EOUT, PZOUT)
 
subroutine saptre(AM1, G1, BGX1, BGY1, BGZ1, AM2, G2, BGX2, BGY2, BGZ2)
 
DOUBLE PRECISION function rndm(RDUMMY)
 
subroutine dtwopd(UMO, ECM1, ECM2, PCM1, PCM2, COD1, COF1, SIF1, COD2, COF2, SIF2, AM1, AM2)
 
DOUBLE PRECISION function sampex(X1, X2)
 
subroutine dtrafo(GAM, BGAM, CX, CY, CZ, COD, COF, SIF, P, ECM, PL, CXL, CYL, CZL, EL)
 
subroutine lortmo(N, GAM, BGX, BGY, BGZ)
 
subroutine cormom(AMCH1, AMCH2, AMCH1N, AMCH2N, PQ1X, PQ1Y, PQ1Z, PQ1E, PA1X, PA1Y, PA1Z, PA1E, PQ2X, PQ2Y, PQ2Z, PQ2E, PA2X, PA2Y, PA2Z, PA2E, PXCH1, PYCH1, PZCH1, ECH1, PXCH2, PYCH2, PZCH2, ECH2, IREJ)
 
subroutine diqvs(ECM, IPV, J, IREJ)
 
subroutine flahad(ITYP, IBAR, IF1, IF2, IF3)
 
subroutine fer4mt(IT, PFERM, PXT, PYT, PZT, ET, KT)
 
DOUBLE PRECISION function sampxb(X1, X2, B)
 
DOUBLE PRECISION function dbetar(GAM, ETA)
 
subroutine selpt(PTXSQ1, PTYSQ1, PLQ1, EQ1, PTXSA1, PTYSA1, PLAQ1, EAQ1, PTXSQ2, PTYSQ2, PLQ2, EQ2, PTXSA2, PTYSA2, PLAQ2, EAQ2, AMCH1, AMCH2, IREJ, IKVALA, PTTQ1, PTTA1, PTTQ2, PTTA2, NSELPT)
 
subroutine selpt4(PTXSQ1, PTYSQ1, PLQ1, EQ1, PTXSA1, PTYSA1, PLAQ1, EAQ1, PTXSQ2, PTYSQ2, PLQ2, EQ2, PTXSA2, PTYSA2, PLAQ2, EAQ2, AMCH1, AMCH2, IREJ, IKVALA, PTTQ1, PTTA1, NSELPT)
 
subroutine parpt(IFL, PT1, PT2, IPT, NEVT)
 
subroutine fer4m(PFERM, PXT, PYT, PZT, ET, KT)
 
subroutine sewew(IOP, NHKKH1)
 
subroutine sttran(XO, YO, ZO, CDE, SDE, SFE, CFE, X, Y, Z)
 
subroutine corval(AMMM, IREJ, AMCH1, AMCH2, QTX1, QTY1, QZ1, QE1, QTX2, QTY2, QZ2, QE2, NORIG)
 
subroutine diqssd(ECM, ITS, IPS, IREJ)
 
subroutine daltra(GA, BGX, BGY, BGZ, PCX, PCY, PCZ, EC, P, PX, PY, PZ, E)
 
subroutine hadjck(NHAD, AMCH, PPR, PTA, GAM, BGX, BGY, BGZ, IFB1, IFB2, IFB3, IFB4, I1, I2, NOBAM, NNCH, NORIG, IREJ)
 
subroutine diqdss(ECM, ITS, IPS, IREJ)
 
subroutine diqsv(ECM, ITV, J, IREJ)
 
subroutine dechkk(NHKKH1)
 
subroutine samppt(MODE, PT)
 
subroutine fer4mp(IP, PFERM, PXT, PYT, PZT, ET, KT)
 
subroutine comcm2(IQ1, IQ2, IAQ1, IAQ2, NNCH, IREJ, AMCH)
 
DOUBLE PRECISION function sampey(X1, X2)
 
subroutine dsfecf(SFE, CFE)
 
subroutine title(NA, NB, NCA, NCB)
 
static c2_log_p< float_type > & log()
make a *new object 
 
subroutine cobcma(IF1, IF2, IF3, IJNCH, NNCH, IREJ, AMCH, AMCHN, IKET)
 
subroutine dfatpr(IP, PABS)
 
G4double tau(G4double bet, G4double homega, G4double ef, G4double t)
 
static c2_sqrt_p< float_type > & sqrt()
make a *new object 
 
subroutine hkkfil(IST, ID, M1, M2, PX, PY, PZ, E, NHKKAU, KORMO, ICALL)
 
subroutine comcma(IFQ, IFAQ, IJNCH, NNCH, IREJ, AMCH, AMCHN)
 
subroutine lortrp(N, NAUX, GAM, BGX, BGY, BGZ)
 
void print(const std::vector< T > &data)
 
subroutine dthrep(UMO, ECM1, ECM2, ECM3, PCM1, PCM2, PCM3, COD1, COF1, SIF1, COD2, COF2, SIF2, COD3, COF3, SIF3, AM1, AM2, AM3)
 
DOUBLE PRECISION function dbeta(X1, X2, BET)
 
DOUBLE PRECISION function betrej(GAM, ETA, XMIN, XMAX)
 
subroutine dbklas(I, J, K, I8, I10)
 
subroutine dfatta(IT, PABS)
 
subroutine hadjet(NHAD, AMCH, PPR, PTA, GAM, BGX, BGY, BGZ, IFB1, IFB2, IFB3, IFB4, I1, I2, NOBAM, NNCH, NORIG)
 
static c2_exp_p< float_type > & exp()
make a *new object