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