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-9)
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
6721 IF(ipev.GE.4)
WRITE(6,6638)hps
6722 ELSEIF(nselpt.EQ.2)
THEN
6723 IF (nusept.EQ.0)
THEN
6725 IF(ipev.GE.4)
WRITE(6,6638)hps
6728 ELSEIF(nusept.EQ.1)
THEN
6734 ptxsq1=qtxsq1+hps*cfe
6735 ptysq1=qtysq1+hps*sfe
6736 ptxsa1=qtxsa1-hps*cfe
6737 ptysa1=qtysa1-hps*sfe
6745 IF (ipev.GE.7)
WRITE(6,1000)ptxsq1,ptysq1,ptxsa1,ptysa1 ,ptxsq2,
6746 +ptysq2,ptxsa2,ptysa2
6747 1000
FORMAT (
' PT S ',8f12.6)
6749 pttq1=ptxsq1**2+ptysq1**2
6750 ptta1=ptxsa1**2+ptysa1**2
6752 IF ( nselpt.EQ.0.OR.umo.LE.20.d0 )
THEN
6755 IF (ikvala.EQ.1)b33=3.0+6.0/log10(umo+10.)
6756 IF (ikvala.EQ.2)b33=4.0+3.0/log10(umo+10.)
6760 IF (icoun2.EQ.48)
THEN
6763 IF (icoun2.EQ.50)
THEN
6783 es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
6784 hps=
sqrt(es*es+2.*es*0.94)
6787 IF (.NOT.intpt) hps=0.0000001
6789 ELSEIF(nselpt.EQ.1)
THEN
6790 IF (musept.EQ.0)
THEN
6792 IF(ipev.GE.4)
WRITE(6,6638)hps
6793 6638
FORMAT (
' SELPT:SAMPPT: HPS= ',e12.4)
6796 ELSEIF(musept.EQ.1)
THEN
6802 ptxsq2=qtxsq2+hps*cfe
6803 ptysq2=qtysq2+hps*sfe
6804 ptxsa2=qtxsa2-hps*cfe
6805 ptysa2=qtysa2-hps*sfe
6813 IF (ipev.GE.7)
WRITE(6,1000)ptxsq1,ptysq1,ptxsa1,ptysa1 ,ptxsq2,
6814 +ptysq2,ptxsa2,ptysa2
6816 pttq1=ptxsq1**2+ptysq1**2
6817 ptta1=ptxsa1**2+ptysa1**2
6818 pttq2=ptxsq2**2+ptysq2**2
6819 ptta2=ptxsa2**2+ptysa2**2
6824 IF(plq1.GT.ptwq1.AND.abs(plaq2).GT.ptwq1)
THEN
6827 ELSEIF(plq1.GT.ptwa2.AND.abs(plaq2).GT.ptwa2)
THEN
6831 IF(plaq1.GT.ptwa1.AND.abs(plq2).GT.ptwa1)
THEN
6834 ELSEIF(plaq1.GT.ptwq2.AND.abs(plq2).GT.ptwq2)
THEN
6843 ptta1=ptta1+plaq1**2
6845 ptta2=ptta2+plaq2**2
6848 IF(amte1.GE.eq1**2)amte1=eq1**2/2.
6850 IF(amte2.GE.eq2**2)amte2=eq2**2/2.
6852 IF(amte1.GE.eaq1**2)amte1=eaq1**2/2.
6854 IF(amte2.GE.eaq2**2)amte2=eaq2**2/2.
6855 IF((eq1**2-amte1.LE.pttq1).OR.
6856 * (eq2**2-amte1.LE.pttq2)
6857 * .OR.(eaq1**2-amte3.LE.ptta1).OR.
6858 * (eaq2**2-amte4.LE.ptta2))
THEN
6859 IF ( nselpt.EQ.0.OR.umo.LE.20.d0 )
THEN
6863 useptm = useptm * 0.7
6864 IF( usept.GT.0.01d0 .OR. useptm.GT.0.01d0 )
THEN
6866 WRITE(6,*)
' SELPT: JUMP AFTER REDUCTION OF USEPT'
6867 WRITE(6,*)
' SELPT: USEPT,USEPTM,HPS',usept,useptm,hps
6872 IF(ipev.GE.4)
WRITE(6,6634) ptxsq1,ptysq1,
6873 +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
6874 +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
6876 6634
FORMAT(
' selpt rejec: ',
6877 +
' PTXSQ1,PTYSQ1, +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2,
6878 + PTYSQ2,PLQ2,EQ2,PTXSA2, PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2,IREJ,
6879 +IKVALA, NSELPT ', 2(8e12.4/),2e12.4,3i5)
6888 IF(ip.GE.1)go to 1779
6889 qqq1=qtxsq1**2+qtysq1**2+qlq1**2-pttq1
6890 IF(qqq1.GT.0.d0)
THEN
6893 plq1=
sqrt(eq1**2-pttq1)
6895 qqa1=qtxsa1**2+qtysa1**2+qlaq1**2-ptta1
6896 IF(qqa1.GT.0.d0)
THEN
6899 plaq1=
sqrt(eaq1**2-ptta1)
6901 qqq2=qtxsq2**2+qtysq2**2+qlq2**2-pttq2
6902 IF(qqq2.GT.0.d0)
THEN
6905 plq2=-
sqrt(eq2**2-pttq2)
6907 qqa2=qtxsa2**2+qtysa2**2+qlaq2**2-ptta2
6908 IF(qqa2.GT.0.d0)
THEN
6911 plaq2=-
sqrt(eaq2**2-ptta2)
6917 amch1q=(eq1+eaq2)**2-(ptxsq1+ptxsa2)** 2-(ptysq1+ptysa2)**2-(plq1
6919 IF (amch1q.LE.0.d0)
THEN
6927 301
FORMAT(
' inconsistent Kinematics in SELPT AMCH1Q=',e12.3)
6928 WRITE(6,305) qtxsq1,qtysq1,
6929 +qlq1,qeq1,qtxsa1,qtysa1,qlaq1,qeaq1, qtxsq2,qtysq2,qlq2,qeq2,
6931 +qtysa2,qlaq2,qeaq2, amch1,amch2
6932 305
FORMAT(
'PTXSQ1,PTYSQ1,
6933 +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2,PTYSQ2,PLQ2,EQ2,PTXSA2,
6934 +PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2',5(4e15.5/))
6937 IF(ipev.GE.4)
WRITE(6,6635) ptxsq1,ptysq1,
6938 +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
6939 +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
6941 6635
FORMAT(
' selpt rejec: ',
6942 +
' PTXSQ1,PTYSQ1, +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2,
6943 + PTYSQ2,PLQ2,EQ2,PTXSA2, PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2,IREJ,
6944 +IKVALA, NSELPT ', 2(8e12.4/),2e12.4,3i5)
6950 amch2q=(eq2+eaq1)**2-(ptxsq2+ptxsa1)** 2-(ptysq2+ptysa1)**2-(plq2
6952 IF (amch2q.LE.0.d0)
THEN
6960 302
FORMAT(
' inconsistent Kinematics in SELPT AMCH2Q=',e12.3)
6961 WRITE(6,305) qtxsq1,qtysq1,
6962 +qlq1,qeq1,qtxsa1,qtysa1,qlaq1,qeaq1, qtxsq2,qtysq2,qlq2,qeq2,
6964 +qtysa2,qlaq2,qeaq2, amch1,amch2
6967 IF(ipev.GE.4)
WRITE(6,6636) ptxsq1,ptysq1,
6968 +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
6969 +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
6971 6636
FORMAT(
' selpt rejec: ',
6972 +
' PTXSQ1,PTYSQ1, +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2,
6973 + PTYSQ2,PLQ2,EQ2,PTXSA2, PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2,IREJ,
6974 +IKVALA, NSELPT ', 2(8e12.4/),2e12.4,3i5)
6978 IF(ipev.GE.4)
WRITE(6,6637) ptxsq1,ptysq1,
6979 +plq1,eq1,ptxsa1,ptysa1,plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,
6980 +ptysa2,plaq2,eaq2, amch1,amch2,irej,ikvala,pttq1,ptta1,
6982 6637
FORMAT(
' selpt exit : ',
6983 +
' PTXSQ1,PTYSQ1, +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2,
6984 + PTYSQ2,PLQ2,EQ2,PTXSA2, PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2,IREJ,
6985 +IKVALA, NSELPT ', 2(8e12.4/),2e12.4,2i5,2e12.4,i5)
6993 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6999 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
7089 parameter(idmax9=602)
7091 COMMON /ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
7093 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
7105 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
7106 +iibar(210),k1(210),k2(210)
7109 dimension ecmf(3),pcmf(3),codf(3),coff(3),siff(3),itf(3)
7110 dimension ecmff(3),pcmff(3)
7111 dimension cxf(3),cyf(3),czf(3)
7116 IF (iphkk.GE.2)
WRITE(6,1000) ihkk,nhkk
7117 1000
FORMAT(
' DECHKK IHKK,NHKK= ',2i5)
7122 IF (ihkk.GT.nhkk)
THEN
7127 IF(abs(isthkk(ihkk)).NE.1) goto 10
7132 IF (it.LT.1.OR.it.GT.210)
THEN
7134 1003
FORMAT (
' DECHKK IT ',i10)
7145 IF(it.EQ.135.OR.it.EQ.136) goto 10
7146 IF(it.GE.1.AND.it.LE.7) goto 10
7147 ELSEIF(istab.EQ.2)
THEN
7148 IF(it.GE. 1.AND.it.LE. 30) goto 10
7149 IF(it.GE. 97.AND.it.LE.103) goto 10
7150 IF(it.GE.115.AND.it.LE.122) goto 10
7151 IF(it.GE.131.AND.it.LE.136) goto 10
7152 IF(it.EQ.109) goto 10
7153 IF(it.GE.137.AND.it.LE.160) goto 10
7154 ELSEIF(istab.EQ.3)
THEN
7155 IF(it.GE.1.AND.it.LE.23) goto 10
7156 IF(it.GE. 97.AND.it.LE.103) goto 10
7157 IF(it.EQ.109.AND.it.EQ.115) goto 10
7158 IF(it.GE.133.AND.it.LE.136) goto 10
7162 pls=
sqrt(abs(phkk(1,ihkk)**2+phkk(2,ihkk)**2+phkk(3,ihkk)**2))
7166 amtest=
sqrt(abs(phkk(4,ihkk)**2-pls**2))
7167 IF(abs(amtest-phkk(5,ihkk)).GE.1.
d-3)
THEN
7170 plss=(phkk(4,ihkk)**2-phkk(5,ihkk))
7171 IF(plss.LE.0.d0)
THEN
7172 WRITE(6,
'(A)')
' negative momentum square!'
7177 phkk(1,ihkk)=phkk(1,ihkk)*amodp
7178 phkk(2,ihkk)=phkk(2,ihkk)*amodp
7179 phkk(3,ihkk)=phkk(3,ihkk)*amodp
7183 IF(pls.NE.0.d0)
THEN
7184 cxs=phkk(1,ihkk)/pls
7185 cys=phkk(2,ihkk)/pls
7186 czs=phkk(3,ihkk)/pls
7197 IF (vv.GT.wt(iik)) go to 20
7205 IF (itf(2).LT.1) go to 10
7209 IF(iphkk.GE.1)
WRITE(6,1010)it,iik,itf(1),itf(2),itf(3)
7210 1010
FORMAT(
' DECHKK IT,IIK,IT1,IT2,IT3 ',5i5)
7213 IF(itf(3).EQ.0)
THEN
7215 CALL
dtwopd(eco,ecmf(1),ecmf(2),pcmf(1),pcmf(2), codf(1),coff
7216 + (1),siff(1),codf(2),coff(2),siff(2), aam(itf(1)),aam(itf(2)))
7217 sid1=
sqrt(abs((1.-codf(1))*(1.+codf(1))))
7218 sid2=
sqrt(abs((1.-codf(2))*(1.+codf(2))))
7219 pix1=pcmf(1)*sid1*coff(1)
7220 piy1=pcmf(1)*sid1*siff(1)
7221 piz1=pcmf(1)*codf(1)
7222 pix2=pcmf(2)*sid2*coff(2)
7223 piy2=pcmf(2)*sid2*siff(2)
7224 piz2=pcmf(2)*codf(2)
7228 ecm12=ecmf(1)+ecmf(2)-eco
7229 IF((abs(pix12).GT.0.000001d0).OR.
7230 + (abs(piy12).GT.0.000001d0).OR.
7231 + (abs(piz12).GT.0.000001d0).OR.
7232 + (abs(ecm12).GT.0.000001d0))
THEN
7233 WRITE(6,778)pix12,piy12,piz12,ecm12
7234 778
FORMAT(
' DWOPD px,py,pz,e',4f10.6)
7239 CALL
dthrep(eco,ecmf(1),ecmf(2),ecmf(3),pcmf(1),pcmf(2),pcmf(3),
7240 + codf(1),coff(1),siff(1),codf(2),coff(2),siff(2), codf(3),coff
7241 + (3),siff(3), aam(itf(1)),aam(itf(2)),aam(itf(3)))
7242 sid1=
sqrt((1.-codf(1))*(1.+codf(1)))
7243 sid2=
sqrt((1.-codf(2))*(1.+codf(2)))
7244 sid3=
sqrt((1.-codf(3))*(1.+codf(3)))
7245 pix1=pcmf(1)*sid1*coff(1)
7246 piy1=pcmf(1)*sid1*siff(1)
7247 piz1=pcmf(1)*codf(1)
7248 pix2=pcmf(2)*sid2*coff(2)
7249 piy2=pcmf(2)*sid2*siff(2)
7250 piz2=pcmf(2)*codf(2)
7251 pix3=pcmf(3)*sid3*coff(3)
7252 piy3=pcmf(3)*sid3*siff(3)
7253 piz3=pcmf(3)*codf(3)
7254 pix12=pix1+pix2+pix3
7255 piy12=piy1+piy2+piy3
7256 piz12=piz1+piz2+piz3
7257 ecm12=ecmf(1)+ecmf(2)+ecmf(3)-eco
7258 IF((abs(pix12).GT.0.000001d0).OR.
7259 + (abs(piy12).GT.0.000001d0).OR.
7260 + (abs(piz12).GT.0.000001d0).OR.
7261 + (abs(ecm12).GT.0.000001d0))
THEN
7262 WRITE(6,779)pix12,piy12,piz12,ecm12
7263 779
FORMAT(
' DTHEPD px,py,pz,e',4f10.6)
7268 jdahkk(1,ihkk)=nhkk + 1
7269 jdahkk(2,ihkk)=nhkk + ndecpr
7271 ehecc=
sqrt(abs(pcmf(id)** 2+ aam(itf(id))**2))
7272 IF (abs(ehecc-ecmf(id)).GT.0.0001d0)
THEN
7273 WRITE(6,
'(2A/3I5,3E15.6)')
7274 &
' DECHKK: CORRECT INCONSISTENT ENERGY ',
7275 *
' IHKK,NHKK,ITF(ID), ECMF(ID),EHECC, AAM(ITF(ID))',
7276 * ihkk,nhkk,itf(id), ecmf(id),ehecc, aam(itf(id))
7280 CALL
dtrafo(gam,bgam,cxs,cys,czs, codf(id),coff(id),
7281 * siff(id),pcmf(id),ecmf(id), pcmff(id),cxf(id),
7282 *cyf(id),czf(id),ecmff(id))
7283 IF (iphkk.GE.2)
WRITE(6,
'(A,7E15.5/8E15.5)')
' DTRAFO ',
7284 * gam,bgam,cxs,cys,czs, codf(id),coff(id),
7285 * siff(id),pcmf(id),ecmf(id), pcmff(id),cxf(id),
7286 *cyf(id),czf(id),ecmff(id)
7293 WRITE (6,1020)nhkk,
nmxhkk
7294 1020
FORMAT (.GT.
' NHKKNMXHKK IN DECHKK RETURN ',2i10)
7301 WRITE (6,
'(A,2I5)') .EQ.
' DECHKK: NHKKNMXHKK ',nhkk,
nmxhkk
7307 idhkk(nhkk)=
mpdgha(itf(id))
7312 phkk(1,nhkk)=cxf(id)*pcmff(id)
7313 phkk(2,nhkk)=cyf(id)*pcmff(id)
7314 phkk(3,nhkk)=czf(id)*pcmff(id)
7315 ehecc=
sqrt(abs(pcmff(id)** 2+ aam(itf(id))**2))
7316 IF (abs(ehecc-ecmff(id)).GT.0.003d0)
THEN
7317 WRITE(6,
'(2A/3I5,3E15.6)')
7318 &
' DECHKK: CORRECT INCONSISTENT ENERGY ',
7319 *
' IHKK,NHKK,ITF(ID), ECMF(ID),EHECC, AAM(ITF(ID))',
7320 * ihkk,nhkk,itf(id), ecmff(id),ehecc, aam(itf(id))
7323 phkk(4,nhkk)=ecmff(id)
7324 phkk(5,nhkk)=aam(itf(id))
7325 vhkk(1,nhkk)=vhkk(1,ihkk)
7326 vhkk(2,nhkk)=vhkk(2,ihkk)
7327 vhkk(3,nhkk)=vhkk(3,ihkk)
7328 vhkk(4,nhkk)=vhkk(4,ihkk)
7330 IF (iphkk.GE.7)
WRITE(6,1030)nhkk, isthkk(nhkk),idhkk(nhkk),
7331 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
7332 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
7334 1030
FORMAT (i6,i4,5i6,9e10.2)
7346 SUBROUTINE dtrafo(GAM,BGAM,CX,CY,CZ,COD,COF,SIF,P,ECM,
7348 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7351 sid=
sqrt(1.d0-cod*cod)
7354 pppt=
sqrt(plx**2+ply**2)
7356 plz=gam*pcmz+bgam*ecm
7357 pl=
sqrt(plx*plx+ply*ply+plz*plz)
7358 el=gam*ecm+bgam*pcmz
7364 CALL
sttran(cx,cy,cz,coz,siz,sif,cof,cxl,cyl,czl)
7369 SUBROUTINE sttran(XO,YO,ZO,CDE,SDE,SFE,CFE,X,Y,Z)
7370 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7400 IF (
a .LT. anglsq )
THEN
7410 x=-yo*xi/
a-zo*xo*yi/
a+xo*zi
7411 y=xo*xi/
a-zo*yo*yi/
a+yo*zi
7421 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7428 parameter(nfimax=249)
7429 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
7430 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
7431 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
7435 parameter(tiny=1.
d-10)
7459 CALL
daltra(gam,bgx,bgy,bgz,pxi,pyi,pzi,eei, ppa,pxf(i),pyf(i),
7468 CALL
daltra(gam,bgx,bgy,bgz,pxsm,pysm,pzsm,esum, ppa,pxsm,pysm,
7477 diffl=pxdif+pydif+pzdif+edif
7478 IF(esum.LT.tiny)esum=tiny
7480 IF(diffl.GE.1.
d-4)
WRITE(6,1000)num,pxdif,pydif,pzdif,edif,pxsm,
7481 +pxsc, pysm,pysc,pzsm,pzsc,esum,esmc
7482 1000
FORMAT(
' ',2
x,
'LORTRA:NUM=',i5,2
x,
'PXDIF=',1pe15.6,2
x,
'PYDIF=', 1
7483 +pe15.6,2
x,
'PZDIF=',1pe15.6,2
x,
'EDIF=',1pe15.6/2
x,
'PXSM=',1pe15.6,2
7484 +
x,
'PXSC=',1pe15.6,2
x,
'PYSM=',1pe15.6,2
x,
'PYSC=',1pe15.6/2
x,
'PZSM',
7485 +1pe15.6,2
x,
'PZSC=',1pe15.6,2
x,
'ESUM=',1pe15.6,2
x,
'ESMC=',1pe15.6/2
7486 +
x,
'LORTRA DIFFERENCES DUE TO ALTRA'/)
7497 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7507 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
7514 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
7515 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
7516 +prebin,taebin,fermod,etacou
7528 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
7529 +iibar(210),k1(210),k2(210)
7535 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
7537 * ,xpsu(248),xtsu(248)
7538 * ,xpsut(248),xtsut(248)
7539 COMMON /intnez/ ndz,nzd
7541 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
7542 +ixpv,ixps,ixtv,ixts, intvv1(248),
7543 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
7545 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
7559 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
7565 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
7567 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
7568 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
7575 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
7578 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
7587 COMMON /nncms/ gamcm,bgcm,umo,pcm,eproj,pproj
7589 COMMON /abrsv/ amcsv1(248),amcsv2(248),gacsv1(248),gacsv2(248),
7590 +bgxsv1(248),bgysv1(248),bgzsv1(248), bgxsv2(248),bgysv2(248),
7591 +bgzsv2(248), nchsv1(248),nchsv2(248),ijcsv1(248),ijcsv2(248),
7592 +pqsva1(248,4),pqsva2(248,4), pqsvb1(248,4),pqsvb2(248,4)
7594 COMMON /abrvs/ amcvs1(248),amcvs2(248),gacvs1(248),gacvs2(248),
7595 +bgxvs1(248),bgyvs1(248),bgzvs1(248), bgxvs2(248),bgyvs2(248),
7596 +bgzvs2(248), nchvs1(248),nchvs2(248),ijcvs1(248),ijcvs2(248),
7597 +pqvsa1(248,4),pqvsa2(248,4), pqvsb1(248,4),pqvsb2(248,4)
7599 COMMON /abrvv/ amcvv1(248),amcvv2(248),gacvv1(248),gacvv2(248),
7600 +bgxvv1(248),bgyvv1(248),bgzvv1(248), bgxvv2(248),bgyvv2(248),
7601 +bgzvv2(248), nchvv1(248),nchvv2(248),ijcvv1(248),ijcvv2(248),
7602 +pqvva1(248,4),pqvva2(248,4), pqvvb1(248,4),pqvvb2(248,4)
7604 COMMON /abrdv/ amcdv1(248),amcdv2(248),gacdv1(248),gacdv2(248),
7605 +bgxdv1(248),bgydv1(248),bgzdv1(248), bgxdv2(248),bgydv2(248),
7606 +bgzdv2(248), nchdv1(248),nchdv2(248),ijcdv1(248),ijcdv2(248),
7607 +pqdva1(248,4),pqdva2(248,4), pqdvb1(248,4),pqdvb2(248,4)
7610 COMMON /abrvd/ amcvd1(248),amcvd2(248),gacvd1(248),gacvd2(248),
7611 +bgxvd1(248),bgyvd1(248),bgzvd1(248), bgxvd2(248),bgyvd2(248),
7612 +bgzvd2(248), nchvd1(248),nchvd2(248),ijcvd1(248),ijcvd2(248),
7613 +pqvda1(248,4),pqvda2(248,4), pqvdb1(248,4),pqvdb2(248,4)
7615 COMMON /abrds/ amcds1(248),amcds2(248),gacds1(248),gacds2(248),
7616 +bgxds1(248),bgyds1(248),bgzds1(248), bgxds2(248),bgyds2(248),
7617 +bgzds2(248), nchds1(248),nchds2(248),ijcds1(248),ijcds2(248),
7618 +pqdsa1(248,4),pqdsa2(248,4), pqdsb1(248,4),pqdsb2(248,4)
7634 COMMON /abrsd/ amcsd1(248),amcsd2(248),gacsd1(248),gacsd2(248),
7635 +bgxsd1(248),bgysd1(248),bgzsd1(248), bgxsd2(248),bgysd2(248),
7636 +bgzsd2(248), nchsd1(248),nchsd2(248),ijcsd1(248),ijcsd2(248),
7637 +pqsda1(248,4),pqsda2(248,4), pqsdb1(248,4),pqsdb2(248,4)
7653 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
7655 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
7656 +ipadis,ishmal,lpauli
7658 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
7660 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
7677 COMMON /nucjtn/nonuj1,nonujt,nonus1,nonust
7762 IF(abs(nchss1(
n)).NE.99)
THEN
7763 pxss=pxss + pqssa1(
n,1) + pqssa2(
n,1)
7764 pyss=pyss + pqssa1(
n,2) + pqssa2(
n,2)
7765 pzss=pzss + pqssa1(
n,3) + pqssa2(
n,3)
7766 pess=pess + pqssa1(
n,4) + pqssa2(
n,4)
7768 IF(abs(nchss2(
n)).NE.99)
THEN
7769 pxss=pxss + pqssb1(
n,1) + pqssb2(
n,1)
7770 pyss=pyss + pqssb1(
n,2) + pqssb2(
n,2)
7771 pzss=pzss + pqssb1(
n,3) + pqssb2(
n,3)
7772 pess=pess + pqssb1(
n,4) + pqssb2(
n,4)
7776 pzbss=gamcm*pzss + bgcm*pess
7777 pebss=gamcm*pess + bgcm*pzss
7785 IF(abs(nchsv1(
n)).NE.99)
THEN
7786 pxsv=pxsv +pqsva1(
n,1)+pqsva2(
n,1)
7787 pysv=pysv +pqsva1(
n,2)+pqsva2(
n,2)
7788 pzsv=pzsv +pqsva1(
n,3)+pqsva2(
n,3)
7789 pesv=pesv +pqsva1(
n,4)+pqsva2(
n,4)
7792 + pxsv,pysv,pzsv,pesv
7798 IF(abs(nchsv2(
n)).NE.99)
THEN
7799 pxsv=pxsv + pqsvb1(
n,1)+pqsvb2(
n,1)
7800 pysv=pysv + pqsvb1(
n,2)+pqsvb2(
n,2)
7801 pzsv=pzsv + pqsvb1(
n,3)+pqsvb2(
n,3)
7802 pesv=pesv + pqsvb1(
n,4)+pqsvb2(
n,4)
7805 + pxsv,pysv,pzsv,pesv
7809 pzbsv=gamcm*pzsv + bgcm*pesv
7810 pebsv=gamcm*pesv + bgcm*pzsv
7813 + pxsv,pysv,pzbsv,pebsv
7822 IF(abs(nchvs1(
n)).NE.99)
THEN
7823 pxvs=pxvs + pqvsa1(
n,1) + pqvsa2(
n,1)
7824 pyvs=pyvs + pqvsa1(
n,2) + pqvsa2(
n,2)
7825 pzvs=pzvs + pqvsa1(
n,3) + pqvsa2(
n,3)
7826 pevs=pevs + pqvsa1(
n,4) + pqvsa2(
n,4)
7828 IF(abs(nchvs2(
n)).NE.99)
THEN
7829 pxvs=pxvs + pqvsb1(
n,1) + pqvsb2(
n,1)
7830 pyvs=pyvs + pqvsb1(
n,2) + pqvsb2(
n,2)
7831 pzvs=pzvs + pqvsb1(
n,3) + pqvsb2(
n,3)
7832 pevs=pevs + pqvsb1(
n,4) + pqvsb2(
n,4)
7835 pzbvs=gamcm*pzvs + bgcm*pevs
7836 pebvs=gamcm*pevs + bgcm*pzvs
7842 IF(abs(nchds1(
n)).NE.99)
THEN
7843 pxds=pxds + pqdsa1(
n,1) + pqdsa2(
n,1)
7844 pyds=pyds + pqdsa1(
n,2) + pqdsa2(
n,2)
7845 pzds=pzds + pqdsa1(
n,3) + pqdsa2(
n,3)
7846 peds=peds + pqdsa1(
n,4) + pqdsa2(
n,4)
7848 IF(abs(nchds2(
n)).NE.99)
THEN
7849 pxds=pxds + pqdsb1(
n,1) + pqdsb2(
n,1)
7850 pyds=pyds + pqdsb1(
n,2) + pqdsb2(
n,2)
7851 pzds=pzds + pqdsb1(
n,3) + pqdsb2(
n,3)
7852 peds=peds + pqdsb1(
n,4) + pqdsb2(
n,4)
7855 pzbds=gamcm*pzds + bgcm*peds
7856 pebds=gamcm*peds + bgcm*pzds
7862 IF(abs(nchdz1(
n)).NE.99)
THEN
7863 pxdz=pxdz + pqdza1(
n,1) + pqdza2(
n,1)
7864 pydz=pydz + pqdza1(
n,2) + pqdza2(
n,2)
7865 pzdz=pzdz + pqdza1(
n,3) + pqdza2(
n,3)
7866 pedz=pedz + pqdza1(
n,4) + pqdza2(
n,4)
7868 IF(abs(nchdz2(
n)).NE.99)
THEN
7869 pxdz=pxdz + pqdzb1(
n,1) + pqdzb2(
n,1)
7870 pydz=pydz + pqdzb1(
n,2) + pqdzb2(
n,2)
7871 pzdz=pzdz + pqdzb1(
n,3) + pqdzb2(
n,3)
7872 pedz=pedz + pqdzb1(
n,4) + pqdzb2(
n,4)
7875 pzbdz=gamcm*pzdz + bgcm*pedz
7876 pebdz=gamcm*pedz + bgcm*pzdz
7882 IF(abs(nchsd1(
n)).NE.99)
THEN
7883 pxsd=pxsd + pqsda1(
n,1) + pqsda2(
n,1)
7884 pysd=pysd + pqsda1(
n,2) + pqsda2(
n,2)
7885 pzsd=pzsd + pqsda1(
n,3) + pqsda2(
n,3)
7886 pesd=pesd + pqsda1(
n,4) + pqsda2(
n,4)
7888 IF(abs(nchsd2(
n)).NE.99)
THEN
7889 pxsd=pxsd + pqsdb1(
n,1) + pqsdb2(
n,1)
7890 pysd=pysd + pqsdb1(
n,2) + pqsdb2(
n,2)
7891 pzsd=pzsd + pqsdb1(
n,3) + pqsdb2(
n,3)
7892 pesd=pesd + pqsdb1(
n,4) + pqsdb2(
n,4)
7895 pzbsd=gamcm*pzsd + bgcm*pesd
7896 pebsd=gamcm*pesd + bgcm*pzsd
7902 IF(abs(nchzd1(
n)).NE.99)
THEN
7903 pxzd=pxzd + pqzda1(
n,1) + pqzda2(
n,1)
7904 pyzd=pyzd + pqzda1(
n,2) + pqzda2(
n,2)
7905 pzzd=pzzd + pqzda1(
n,3) + pqzda2(
n,3)
7906 pezd=pezd + pqzda1(
n,4) + pqzda2(
n,4)
7908 IF(abs(nchzd2(
n)).NE.99)
THEN
7909 pxzd=pxzd + pqzdb1(
n,1) + pqzdb2(
n,1)
7910 pyzd=pyzd + pqzdb1(
n,2) + pqzdb2(
n,2)
7911 pzzd=pzzd + pqzdb1(
n,3) + pqzdb2(
n,3)
7912 pezd=pezd + pqzdb1(
n,4) + pqzdb2(
n,4)
7915 pzbzd=gamcm*pzzd + bgcm*pezd
7916 pebzd=gamcm*pezd + bgcm*pzzd
7922 IF(abs(nchdv1(
n)).NE.99)
THEN
7923 pxdv=pxdv + pqdva1(
n,1) + pqdva2(
n,1)
7924 pydv=pydv + pqdva1(
n,2) + pqdva2(
n,2)
7925 pzdv=pzdv + pqdva1(
n,3) + pqdva2(
n,3)
7926 pedv=pedv + pqdva1(
n,4) + pqdva2(
n,4)
7928 IF(abs(nchdv2(
n)).NE.99)
THEN
7929 pxdv=pxdv + pqdvb1(
n,1) + pqdvb2(
n,1)
7930 pydv=pydv + pqdvb1(
n,2) + pqdvb2(
n,2)
7931 pzdv=pzdv + pqdvb1(
n,3) + pqdvb2(
n,3)
7932 pedv=pedv + pqdvb1(
n,4) + pqdvb2(
n,4)
7935 pzbdv=gamcm*pzdv + bgcm*pedv
7936 pebdv=gamcm*pedv + bgcm*pzdv
7942 IF(abs(nchvd1(
n)).NE.99)
THEN
7943 pxvd=pxvd + pqvda1(
n,1) + pqvda2(
n,1)
7944 pyvd=pyvd + pqvda1(
n,2) + pqvda2(
n,2)
7945 pzvd=pzvd + pqvda1(
n,3) + pqvda2(
n,3)
7946 pevd=pevd + pqvda1(
n,4) + pqvda2(
n,4)
7948 IF(abs(nchvd2(
n)).NE.99)
THEN
7949 pxvd=pxvd + pqvdb1(
n,1) + pqvdb2(
n,1)
7950 pyvd=pyvd + pqvdb1(
n,2) + pqvdb2(
n,2)
7951 pzvd=pzvd + pqvdb1(
n,3) + pqvdb2(
n,3)
7952 pevd=pevd + pqvdb1(
n,4) + pqvdb2(
n,4)
7955 pzbvd=gamcm*pzvd + bgcm*pevd
7956 pebvd=gamcm*pevd + bgcm*pzvd
7964 IF((nchvv1(
n).NE.99).AND.(nchvv2(
n).NE.99))
THEN
7965 pxvv=pxvv+pqvva1(
n,1)+pqvva2(
n,1)+pqvvb1(
n,1)+pqvvb2(
n,1)
7966 pyvv=pyvv+pqvva1(
n,2)+pqvva2(
n,2)+pqvvb1(
n,2)+pqvvb2(
n,2)
7967 pzvv=pzvv+pqvva1(
n,3)+pqvva2(
n,3)+pqvvb1(
n,3)+pqvvb2(
n,3)
7968 pevv=pevv+pqvva1(
n,4)+pqvva2(
n,4)+pqvvb1(
n,4)+pqvvb2(
n,4)
7971 pzbvv=gamcm*pzvv + bgcm*pevv
7972 pebvv=gamcm*pevv + bgcm*pzvv
7995 IF(abs(nchzz1(
n)).NE.99.AND.jhkksx(
n).EQ.1)
THEN
7996 IF(abs(nchzz1(
n)).NE.88)
THEN
7997 pxzz=pxzz + pqzza1(
n,1) + pqzza2(
n,1)
7998 pyzz=pyzz + pqzza1(
n,2) + pqzza2(
n,2)
7999 pzzz=pzzz + pqzza1(
n,3) + pqzza2(
n,3)
8000 pezz=pezz + pqzza1(
n,4) + pqzza2(
n,4)
8003 IF(abs(nchzz2(
n)).NE.99.AND.jhkksx(
n).EQ.1)
THEN
8004 IF(abs(nchzz2(
n)).NE.88)
THEN
8005 pxzz=pxzz + pqzzb1(
n,1) + pqzzb2(
n,1)
8006 pyzz=pyzz + pqzzb1(
n,2) + pqzzb2(
n,2)
8007 pzzz=pzzz + pqzzb1(
n,3) + pqzzb2(
n,3)
8008 pezz=pezz + pqzzb1(
n,4) + pqzzb2(
n,4)
8012 pzbzz=gamcm*pzzz + bgcm*pezz
8013 pebzz=gamcm*pezz + bgcm*pzzz
8020 IF(abs(nchhh1(
n)).NE.99.AND.jhkkex(
n).EQ.1)
THEN
8021 pxhh=pxhh + pqhha1(
n,1) + pqhha2(
n,1)
8022 pyhh=pyhh + pqhha1(
n,2) + pqhha2(
n,2)
8023 pzhh=pzhh + pqhha1(
n,3) + pqhha2(
n,3)
8024 pehh=pehh + pqhha1(
n,4) + pqhha2(
n,4)
8026 IF(abs(nchhh2(
n)).NE.99.AND.jhkkex(
n).EQ.1)
THEN
8027 pxhh=pxhh + pqhhb1(
n,1) + pqhhb2(
n,1)
8028 pyhh=pyhh + pqhhb1(
n,2) + pqhhb2(
n,2)
8029 pzhh=pzhh + pqhhb1(
n,3) + pqhhb2(
n,3)
8030 pehh=pehh + pqhhb1(
n,4) + pqhhb2(
n,4)
8033 pzbhh=gamcm*pzhh + bgcm*pehh
8034 pebhh=gamcm*pehh + bgcm*pzhh
8044 IF(isthkk(i).EQ.11)e0000=e0000+prmom(4,i)
8045 IF(isthkk(i).EQ.11)p0000=p0000+prmom(3,i)
8050 IF(isthkk(i).EQ.12)e0000=e0000+tamom(4,ii)
8051 IF(isthkk(i).EQ.12)p0000=p0000+tamom(3,ii)
8053 p000=gamcm*p0000+bgcm*e0000
8054 e000=gamcm*e0000+bgcm*p0000
8055 iprojo=(pzbal*1.001)/pproj
8056 residu=abs(e000-pebal)/(e000)
8058 WRITE(6,
'(A,2E15.5)')
' E000,PEBAL', e000,pebal
8059 WRITE(6,1000)pxbal,pybal,pzbal,pebal, pxss,pyss,pzbss,pebss,
8060 + pxsv,pysv,pzbsv,pebsv, pxvs,pyvs,pzbvs,pebvs, pxvv,pyvv,pzbvv,
8061 + pebvv,pxcc,pycc,pzbcc,pebcc,
8062 + pxzz,pyzz,pzbzz,pebzz,
8063 + pxhh,pyhh,pzbhh,pebhh,
8064 + pxds,pyds,pzbds,pebds,
8065 + pxsd,pysd,pzbsd,pebsd,
8066 + pxdz,pydz,pzbdz,pebdz,
8067 + pxzd,pyzd,pzbzd,pebzd,
8068 + pxdv,pydv,pzbdv,pebdv,
8069 + pxvd,pyvd,pzbvd,pebvd
8071 IF (residu.GT.0.02d0)
THEN
8074 IF (residu.GT.0.02d0.AND.iphkk.GE.2)
THEN
8076 WRITE(6,
'(A,2E15.5)')
' E000,PEBAL', e000,pebal
8077 WRITE(6,1000)pxbal,pybal,pzbal,pebal, pxss,pyss,pzbss,pebss,
8078 + pxsv,pysv,pzbsv,pebsv, pxvs,pyvs,pzbvs,pebvs, pxvv,pyvv,pzbvv,
8079 + pebvv,pxcc,pycc,pzbcc,pebcc,
8080 + pxzz,pyzz,pzbzz,pebzz,
8081 + pxhh,pyhh,pzbhh,pebhh,
8082 + pxds,pyds,pzbds,pebds,
8083 + pxsd,pysd,pzbsd,pebsd,
8084 + pxdz,pydz,pzbdz,pebdz,
8085 + pxzd,pyzd,pzbzd,pebzd,
8086 + pxdv,pydv,pzbdv,pebdv,
8087 + pxvd,pyvd,pzbvd,pebvd
8089 1000
FORMAT (
' 4 MOMENTUM CONS.IN EVENT LEVEL OF PARTONS',/
' ALL',4e15
8090 +.5/,
' SS ',4e15.5/,
' SV ',4e15.5/
' VS ',4e15.5/,
' VV ',4e15.5/,
8153 IF(abs(nchss1(
n)).NE.99)
THEN
8154 pxss=pxss+bgxss1(
n)*amcss1(
n)
8155 pyss=pyss+bgyss1(
n)*amcss1(
n)
8156 pzss=pzss+bgzss1(
n)*amcss1(
n)
8157 pess=pess+gacss1(
n)*amcss1(
n)
8159 IF(abs(nchss2(
n)).NE.99)
THEN
8160 pxss=pxss+bgxss2(
n)*amcss2(
n)
8161 pyss=pyss+bgyss2(
n)*amcss2(
n)
8162 pzss=pzss+bgzss2(
n)*amcss2(
n)
8163 pess=pess+gacss2(
n)*amcss2(
n)
8167 pzbss=gamcm*pzss + bgcm*pess
8168 pebss=gamcm*pess + bgcm*pzss
8188 IF(abs(nchsv1(
n)).NE.99)
THEN
8189 pxsv=pxsv+bgxsv1(
n)*amcsv1(
n)
8190 pysv=pysv+bgysv1(
n)*amcsv1(
n)
8191 pzsv=pzsv+bgzsv1(
n)*amcsv1(
n)
8192 pesv=pesv+gacsv1(
n)*amcsv1(
n)
8194 IF(abs(nchsv2(
n)).NE.99)
THEN
8195 pxsv=pxsv+bgxsv2(
n)*amcsv2(
n)
8196 pysv=pysv+bgysv2(
n)*amcsv2(
n)
8197 pzsv=pzsv+bgzsv2(
n)*amcsv2(
n)
8198 pesv=pesv+gacsv2(
n)*amcsv2(
n)
8201 pzbsv=gamcm*pzsv + bgcm*pesv
8202 pebsv=gamcm*pesv + bgcm*pzsv
8208 IF(abs(nchds1(
n)).NE.99)
THEN
8209 pxds=pxds+bgxds1(
n)*amcds1(
n)
8210 pyds=pyds+bgyds1(
n)*amcds1(
n)
8211 pzds=pzds+bgzds1(
n)*amcds1(
n)
8212 peds=peds+gacds1(
n)*amcds1(
n)
8214 IF(abs(nchds2(
n)).NE.99)
THEN
8215 pxds=pxds+bgxds2(
n)*amcds2(
n)
8216 pyds=pyds+bgyds2(
n)*amcds2(
n)
8217 pzds=pzds+bgzds2(
n)*amcds2(
n)
8218 peds=peds+gacds2(
n)*amcds2(
n)
8221 pzbds=gamcm*pzds + bgcm*peds
8222 pebds=gamcm*peds + bgcm*pzds
8228 IF(abs(nchsd1(
n)).NE.99)
THEN
8229 pxsd=pxsd+bgxsd1(
n)*amcsd1(
n)
8230 pysd=pysd+bgysd1(
n)*amcsd1(
n)
8231 pzsd=pzsd+bgzsd1(
n)*amcsd1(
n)
8232 pesd=pesd+gacsd1(
n)*amcsd1(
n)
8234 IF(abs(nchsd2(
n)).NE.99)
THEN
8235 pxsd=pxsd+bgxsd2(
n)*amcsd2(
n)
8236 pysd=pysd+bgysd2(
n)*amcsd2(
n)
8237 pzsd=pzsd+bgzsd2(
n)*amcsd2(
n)
8238 pesd=pesd+gacsd2(
n)*amcsd2(
n)
8241 pzbsd=gamcm*pzsd + bgcm*pesd
8242 pebsd=gamcm*pesd + bgcm*pzsd
8248 IF(abs(nchdv1(
n)).NE.99)
THEN
8249 pxdv=pxdv+bgxdv1(
n)*amcdv1(
n)
8250 pydv=pydv+bgydv1(
n)*amcdv1(
n)
8251 pzdv=pzdv+bgzdv1(
n)*amcdv1(
n)
8252 pedv=pedv+gacdv1(
n)*amcdv1(
n)
8254 IF(abs(nchdv2(
n)).NE.99)
THEN
8255 pxdv=pxdv+bgxdv2(
n)*amcdv2(
n)
8256 pydv=pydv+bgydv2(
n)*amcdv2(
n)
8257 pzdv=pzdv+bgzdv2(
n)*amcdv2(
n)
8258 pedv=pedv+gacdv2(
n)*amcdv2(
n)
8261 pzbdv=gamcm*pzdv + bgcm*pedv
8262 pebdv=gamcm*pedv + bgcm*pzdv
8268 IF(abs(nchvd1(
n)).NE.99)
THEN
8269 pxvd=pxvd+bgxvd1(
n)*amcvd1(
n)
8270 pyvd=pyvd+bgyvd1(
n)*amcvd1(
n)
8271 pzvd=pzvd+bgzvd1(
n)*amcvd1(
n)
8272 pevd=pevd+gacvd1(
n)*amcvd1(
n)
8274 IF(abs(nchvd2(
n)).NE.99)
THEN
8275 pxvd=pxvd+bgxvd2(
n)*amcvd2(
n)
8276 pyvd=pyvd+bgyvd2(
n)*amcvd2(
n)
8277 pzvd=pzvd+bgzvd2(
n)*amcvd2(
n)
8278 pevd=pevd+gacvd2(
n)*amcvd2(
n)
8281 pzbvd=gamcm*pzvd + bgcm*pevd
8282 pebvd=gamcm*pevd + bgcm*pzvd
8290 IF(abs(nchvs1(
n)).NE.99)
THEN
8291 pxvs=pxvs+bgxvs1(
n)*amcvs1(
n)
8292 pyvs=pyvs+bgyvs1(
n)*amcvs1(
n)
8293 pzvs=pzvs+bgzvs1(
n)*amcvs1(
n)
8294 pevs=pevs+gacvs1(
n)*amcvs1(
n)
8296 IF(abs(nchvs2(
n)).NE.99)
THEN
8297 pxvs=pxvs+bgxvs2(
n)*amcvs2(
n)
8298 pyvs=pyvs+bgyvs2(
n)*amcvs2(
n)
8299 pzvs=pzvs+bgzvs2(
n)*amcvs2(
n)
8300 pevs=pevs+gacvs2(
n)*amcvs2(
n)
8303 pzbvs=gamcm*pzvs + bgcm*pevs
8304 pebvs=gamcm*pevs + bgcm*pzvs
8311 IF(abs(nchzz1(
n)).NE.99.AND.jhkksx(
n).EQ.1)
THEN
8312 pxzz=pxzz+bgxzz1(
n)*amczz1(
n)
8313 pyzz=pyzz+bgyzz1(
n)*amczz1(
n)
8314 pzzz=pzzz+bgzzz1(
n)*amczz1(
n)
8315 pezz=pezz+gaczz1(
n)*amczz1(
n)
8317 IF(abs(nchzz2(
n)).NE.99.AND.jhkksx(
n).EQ.1)
THEN
8318 pxzz=pxzz+bgxzz2(
n)*amczz2(
n)
8319 pyzz=pyzz+bgyzz2(
n)*amczz2(
n)
8320 pzzz=pzzz+bgzzz2(
n)*amczz2(
n)
8321 pezz=pezz+gaczz2(
n)*amczz2(
n)
8324 pzbzz=gamcm*pzzz + bgcm*pezz
8325 pebzz=gamcm*pezz + bgcm*pzzz
8331 IF(abs(nchhh1(
n)).NE.99.AND.jhkkex(
n).EQ.1)
THEN
8332 pxhh=pxhh+bgxhh1(
n)*amchh1(
n)
8333 pyhh=pyhh+bgyhh1(
n)*amchh1(
n)
8334 pzhh=pzhh+bgzhh1(
n)*amchh1(
n)
8335 pehh=pehh+gachh1(
n)*amchh1(
n)
8337 IF(abs(nchhh2(
n)).NE.99.AND.jhkkex(
n).EQ.1)
THEN
8338 pxhh=pxhh+bgxhh2(
n)*amchh2(
n)
8339 pyhh=pyhh+bgyhh2(
n)*amchh2(
n)
8340 pzhh=pzhh+bgzhh2(
n)*amchh2(
n)
8341 pehh=pehh+gachh2(
n)*amchh2(
n)
8344 pzbhh=gamcm*pzhh + bgcm*pehh
8345 pebhh=gamcm*pehh + bgcm*pzhh
8352 IF((nchvv1(
n).NE.99).AND.(nchvv2(
n).NE.99))
THEN
8353 pxvv=pxvv+bgxvv1(
n)*amcvv1(
n)+bgxvv2(
n)*amcvv2(
n)
8354 pyvv=pyvv+bgyvv1(
n)*amcvv1(
n)+bgyvv2(
n)*amcvv2(
n)
8355 pzvv=pzvv+bgzvv1(
n)*amcvv1(
n)+bgzvv2(
n)*amcvv2(
n)
8356 pevv=pevv+gacvv1(
n)*amcvv1(
n)+gacvv2(
n)*amcvv2(
n)
8359 pzbvv=gamcm*pzvv + bgcm*pevv
8360 pebvv=gamcm*pevv + bgcm*pzvv
8367 IF (ipev.GE.1)
WRITE(6,1010)pxbal,pybal,pzbal,
8368 +pebal, pxss,pyss,pzbss,pebss, pxsv,pysv,pzbsv,pebsv, pxvs,pyvs,
8369 +pzbvs,pebvs, pxvv,pyvv,pzbvv,pebvv, pxcc,pycc,pzbcc,pebcc,
8370 + pxds,pyds,pzbds,pebds,
8371 + pxzz,pyzz,pzbzz,pebzz,
8372 + pxhh,pyhh,pzbhh,pebhh,
8373 + pxsd,pysd,pzbsd,pebsd,
8374 + pxdv,pydv,pzbdv,pebdv,
8375 + pxvd,pyvd,pzbvd,pebvd
8376 1010
FORMAT (
' 4 MOMENTUM CONS.IN EVENT LEVEL OF CHAINS',/
' ALL',4e15.
8377 +5/,
' SS ',4e15.5/,
' SV ',4e15.5/
' VS ',4e15.5/,
' VV ',4e15.5/,
8390 SUBROUTINE corval(AMMM,IREJ,AMCH1,AMCH2, QTX1,QTY1,QZ1,QE1,QTX2,
8391 +qty2,qz2,qe2,norig)
8392 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
8401 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
8405 IF(ammm.LE.amch1+amch2+0.4d0)
THEN
8410 ek1=(ammm**2-amch2**2 + amch1**2)/(2.*ammm)
8412 pzk1=
sqrt(ek1**2 - amch1**2)
8414 pzk2=
sqrt(ek2**2 - amch2**2)
8460 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
8481 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
8482 +ixpv,ixps,ixtv,ixts, intvv1(248),
8483 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
8485 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
8498 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx),
8499 * ifrovt(248),itovt(248),ifrost(
intmx),
8500 * jsshs(
intmx),jtshs(
intmx),jhkknp(248),jhkknt(248),
8506 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
8509 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
8512 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
8514 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
8515 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),
8530 COMMON /hardha/nhard1,nhkkha
8533 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
8534 * anndv,annvd,annds,annsd,
8536 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
8538 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
8541 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
8542 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
8543 * acouzz,acouhh,acouds,acousd,
8544 * acoudz,acouzd,acoudi,
8545 * acoudv,acouvd,acoucc
8547 COMMON /pshow/ ipshow
8549 COMMON /harlun/ qlun,iharlu
8550 COMMON /jspar/pxj(1000),pyj(1000),pzj(1000),hej(1000),nnnpj
8551 COMMON /jspa/pxs(40000),pys(40000),pzs(40000),hes(40000),nnnps
8555 parameter(nfimax=249)
8556 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
8557 +hep(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
8558 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
8562 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk),
8566 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
8567 COMMON /projk/ iprojk
8568 COMMON /nucjtn/nonuj1,nonujt,nonus1,nonust
8569 COMMON /gluspl/nugluu,nsgluu
8570 COMMON /nomije/ ptmije(10),nnmije(10)
8572 dimension poj(4),pat(4)
8579 IF (iphkk.GE.2)
WRITE(6,7789)nonujt,ncalhh
8580 7789
FORMAT (
' HADRHH NONUJT,NCALHH ',2i10)
8581 IF (jhkkex(i).EQ.1)
THEN
8583 WRITE (6,7744)i,
intmx
8584 7744
FORMAT (.GT.
' HADRHH IINTMX ',2i10)
8596 pt1=
sqrt(poj(1)**2+poj(2)**2)
8597 pt2=
sqrt(pat(1)**2+pat(2)**2)
8598 CALL
parpt(2,pt1,pt2,6,nevt)
8602 pojpt=
sqrt(poj(2)**2+poj(1)**2)
8603 patpt=
sqrt(pat(1)**2+pat(2)**2)
8605 IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
8607 IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
8610 qlun=min(pojpt,patpt)
8611 IF((qlun.LT.2.5d0).OR.(amjch1(i).LT.5.d0))
THEN
8619 IF (gamjh1(i).LT.0.001d0.OR.amjch1(i).LT.2.d0)
THEN
8621 * i,nhad,amjch1(i),poj,pat,gamjh1(i),bgxjh1(i),
8622 * bgyjh1(i),bgzjh1(i),ifb1,ifb2,ifb3,ifb4,jhkkex(i)
8623 7788
FORMAT (
' HADRHH ',2i5,8e12.2/5e12.2,5i5)
8626 CALL
hadjet(nhad,amjch1(i),poj,pat,gamjh1(i),bgxjh1(i),
8627 * bgyjh1(i),bgzjh1(i),ifb1,ifb2,ifb3,ifb4,
8633 IF(iphkk.GE.3)
WRITE(6,*)
' HADRHH:NHKK,NHKKAU ',nhkk,nhkkau
8634 IF (nhad.GT.nfimax)
THEN
8635 WRITE (6,7755)nhad,nfimax
8636 7755
FORMAT (.GT.
' NHADNFIMAX ',2i10)
8642 WRITE (*,
'(A,2I5/A)') .EQ.
' HADRHH: NHKKNMXHKK ',
8647 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
8648 IF (abs(ehecc-hep(j)).GT.0.001d0)
THEN
8657 pthh=
sqrt(pxf(j)**2+pyf(j)**2)+pthh
8660 IF(ibarf(j).EQ.500)istist=2
8662 * pxf(j),pyf(j),pzf(j),hep(j),nhkkau,iormo(j),9)
8663 IF(idhkk(nhkk).EQ.99999)
WRITE (6,5009)nhkk,nref(j),
8665 IF(iphkk.GE.3)
WRITE(6,*)
' First chain HADRHH'
8666 IF (iphkk.GE.3)
WRITE(6,5001) nhkk,
8667 * isthkk(nhkk),idhkk(nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk),
8668 & jdahkk(1,nhkk),jdahkk(2,nhkk),(phkk(khkk,nhkk),khkk=1,5),
8669 & (vhkk(khkk,nhkk),khkk=1,4)
8677 DO 137 j=nnnps,nnnpsu
8679 IF(j.GT.40000.OR.jj.GT.1000)
THEN
8692 IF (nugluu.EQ.1) go to 5111
8700 pt1=
sqrt(poj(1)**2+poj(2)**2)
8701 pt2=
sqrt(pat(1)**2+pat(2)**2)
8702 CALL
parpt(2,pt1,pt2,6,nevt)
8706 pojpt=
sqrt(poj(2)**2+poj(1)**2)
8707 patpt=
sqrt(pat(1)**2+pat(2)**2)
8709 IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
8711 IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
8714 qlun=min(pojpt,patpt)
8715 IF((qlun.LT.2.5d0).OR.(amjch2(i).LT.5.d0))
THEN
8723 CALL
hadjet(nhad,amjch2(i),poj,pat,gamjh2(i),bgxjh2(i),
8724 * bgyjh2(i),bgzjh2(i),ifb1,ifb2,ifb3,ifb4,
8734 WRITE (*,
'(A,2I5/A)') .EQ.
' HADRHH: NHKKNMXHKK ',
8739 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
8740 IF (abs(ehecc-hep(j)).GT.0.001d0)
THEN
8749 pthh=
sqrt(pxf(j)**2+pyf(j)**2)+pthh
8752 IF(ibarf(j).EQ.500)istist=2
8754 * pxf(j),pyf(j),pzf(j),hep(j),nhkkau,iormo(j),10)
8755 IF(idhkk(nhkk).EQ.99999)
WRITE (6,5009)nhkk,nref(j),
8758 IF (iphkk.GE.7)
WRITE(6,5001) nhkk,
8759 * isthkk(nhkk),idhkk(nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk),
8760 & jdahkk(1,nhkk),jdahkk(2,nhkk),(phkk(khkk,nhkk),khkk=1,5),
8761 & (vhkk(khkk,nhkk),khkk=1,4)
8769 DO 187 j=nnnps,nnnpsu
8771 IF(j.GT.40000.OR.jj.GT.1000)
THEN
8790 5001
FORMAT (i6,i4,5i6,9e10.2)
8791 5003
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
8792 5009
FORMAT (
' NHKK,IDHKK(NHKK) ',3i10)
8799 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
8828 COMMON /nucjtn/nonuj1,nonujt,nonus1,nonust
8830 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
8831 +ixpv,ixps,ixtv,ixts, intvv1(248),
8832 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
8834 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
8847 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx),
8848 * ifrovt(248),itovt(248),ifrost(
intmx),
8849 * jsshs(
intmx),jtshs(
intmx),jhkknp(248),jhkknt(248),
8855 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
8858 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
8861 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
8863 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
8864 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),
8873 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
8874 * anndv,annvd,annds,annsd,
8876 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
8878 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
8881 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
8882 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
8883 * acouzz,acouhh,acouds,acousd,
8884 * acoudz,acouzd,acoudi,
8885 * acoudv,acouvd,acoucc
8887 COMMON /pshow/ ipshow
8889 COMMON /harlun/ qlun,iharlu
8890 COMMON /jspar/pxj(1000),pyj(1000),pzj(1000),hej(1000),nnnpj
8891 COMMON /jspa/pxs(40000),pys(40000),pzs(40000),hes(40000),nnnps
8895 parameter(nfimax=249)
8896 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
8897 +hep(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
8898 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
8902 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk),
8906 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
8907 COMMON /projk/ iprojk
8908 COMMON /gluspl/nugluu,nsgluu
8909 COMMON /nomije/ ptmije(10),nnmije(10)
8911 dimension poj(4),pat(4)
8915 IF(nch1(i).EQ.99.OR.nch1(i).EQ.88)go to 20
8916 IF(nch2(i).EQ.99.OR.nch2(i).EQ.88)go to 20
8919 IF (iphkk.GE.7)
WRITE(6,7789)nonust,ncalzz,jhkksx(i)
8920 7789
FORMAT (
' HADRZZ NONUST,NCALZZ,Jhkksx(i) ',3i10)
8921 IF (jhkksx(i).EQ.1)
THEN
8923 WRITE (6,7744)i,
intmx
8924 7744
FORMAT (.GT.
' HADRZZ IINTMX ',2i10)
8936 pt1=
sqrt(poj(1)**2+poj(2)**2)
8937 pt2=
sqrt(pat(1)**2+pat(2)**2)
8938 CALL
parpt(2,pt1,pt2,5,nevt)
8942 pojpt=
sqrt(poj(2)**2+poj(1)**2)
8943 patpt=
sqrt(pat(1)**2+pat(2)**2)
8945 IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
8947 IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
8950 qlun=min(pojpt,patpt)
8951 IF((qlun.LT.2.5d0).OR.(amcch1(i).LT.5.d0))
THEN
8959 IF (gamch1(i).LT.0.001d0)
WRITE (6,7788)
8960 * i,nhad,amcch1(i),poj,pat,gamch1(i),bgxch1(i),
8961 * bgych1(i),bgzch1(i),ifb1,ifb2,ifb3,ifb4,jhkksx(i)
8962 7788
FORMAT (
' HADRZZ ',2i5,10e12.2/3e12.2,5i5)
8963 CALL
hadjet(nhad,amcch1(i),poj,pat,gamch1(i),bgxch1(i),
8964 * bgych1(i),bgzch1(i),ifb1,ifb2,ifb3,ifb4,
8965 * ijczz1(i),ijczz1(i),3,nchzz1(i),23)
8970 IF (nhad.GT.nfimax)
THEN
8971 WRITE (6,7755)nhad,nfimax
8972 7755
FORMAT (.GT.
' NHADNFIMAX ',2i10)
8978 WRITE (*,
'(A,2I5/A)') .EQ.
' HADRZZ: NHKKNMXHKK ',
8983 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
8984 IF (abs(ehecc-hep(j)).GT.0.001d0)
THEN
8993 ptzz=
sqrt(pxf(j)**2+pyf(j)**2)+ptzz
8996 IF(ibarf(j).EQ.500)istist=2
8998 * pxf(j),pyf(j),pzf(j),hep(j),nhkkau,iormo(j),11)
8999 IF(idhkk(nhkk).EQ.99999)
WRITE (6,5009)nhkk,nref(j),
9002 IF (iphkk.GE.7)
WRITE(6,5001) nhkk,
9003 * isthkk(nhkk),idhkk(nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk),
9004 & jdahkk(1,nhkk),jdahkk(2,nhkk),(phkk(khkk,nhkk),khkk=1,5),
9005 & (vhkk(khkk,nhkk),khkk=1,4)
9013 DO 137 j=nnnps,nnnpsu
9015 IF(j.GT.40000.OR.jj.GT.1000)
THEN
9034 pt1=
sqrt(poj(1)**2+poj(2)**2)
9035 pt2=
sqrt(pat(1)**2+pat(2)**2)
9036 CALL
parpt(2,pt1,pt2,5,nevt)
9040 pojpt=
sqrt(poj(2)**2+poj(1)**2)
9041 patpt=
sqrt(pat(1)**2+pat(2)**2)
9043 IF(pojpt.GE.ptmije(iiii))nnmije(iiii)=
9045 IF(patpt.GE.ptmije(iiii))nnmije(iiii)=
9048 qlun=min(pojpt,patpt)
9049 IF((qlun.LT.2.5d0).OR.(amcch2(i).LT.5.d0))
THEN
9057 CALL
hadjet(nhad,amcch2(i),pat,poj,gamch2(i),bgxch2(i),
9058 * bgych2(i),bgzch2(i),ifb1,ifb2,ifb3,ifb4,
9059 * ijczz2(i),ijczz2(i),3,nchzz2(i),24)
9068 WRITE (*,
'(A,2I5/A)') .EQ.
' HADRZZ: NHKKNMXHKK ',
9073 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
9074 IF (abs(ehecc-hep(j)).GT.0.001d0)
THEN
9083 ptzz=
sqrt(pxf(j)**2+pyf(j)**2)+ptzz
9086 IF(ibarf(j).EQ.500)istist=2
9088 * pxf(j),pyf(j),pzf(j),hep(j),nhkkau,iormo(j),12)
9089 IF(idhkk(nhkk).EQ.99999)
WRITE (6,5009)nhkk,nref(j),
9092 IF (iphkk.GE.7)
WRITE(6,5001) nhkk,
9093 * isthkk(nhkk),idhkk(nhkk),jmohkk(1,nhkk),jmohkk(2,nhkk),
9094 & jdahkk(1,nhkk),jdahkk(2,nhkk),(phkk(khkk,nhkk),khkk=1,5),
9095 & (vhkk(khkk,nhkk),khkk=1,4)
9103 DO 187 j=nnnps,nnnpsu
9105 IF(j.GT.40000.OR.jj.GT.1000)
THEN
9122 5001
FORMAT (i6,i4,5i6,9e10.2)
9123 5003
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
9124 5009
FORMAT (
' NHKK,IDHKK(NHKK) ',3i10)
9136 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
9140 CHARACTER*8 projty,targty
9143 COMMON /user1/
title,projty,targty
9144 COMMON /user2/cmener,sdfrac,ptlar,istruf,isingd,idubld
9146 c=4.*(0.15
d-24+0.01
d-24*
log(cmener))
9148 IF ((
p). eq .(1.d0))
THEN
9153 r=dsqrt(-c*dlog(1.d00-
p))
9164 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
9181 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
9183 * ,xpsu(248),xtsu(248)
9184 * ,xpsut(248),xtsut(248)
9185 common/popcck/pdbck,pdbse,pdbseu,
9186 * ijpock,irejck,ick4,ihad4,ick6,ihad6
9187 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
9188 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
9189 *isea43,isea63,irejao
9191 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
9192 +ixpv,ixps,ixtv,ixts, intvv1(248),
9193 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
9195 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
9209 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
9215 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
9217 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
9218 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
9225 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
9228 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
9233 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
9321 COMMON /abrvs/ amcvs1(248),amcvs2(248),gacvs1(248),gacvs2(248),
9322 +bgxvs1(248),bgyvs1(248),bgzvs1(248), bgxvs2(248),bgyvs2(248),
9323 +bgzvs2(248), nchvs1(248),nchvs2(248),ijcvs1(248),ijcvs2(248),
9324 +pqvsa1(248,4),pqvsa2(248,4), pqvsb1(248,4),pqvsb2(248,4)
9327 parameter(nfimax=249)
9328 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
9329 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
9330 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
9333 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
9335 COMMON /projk/ iprojk
9337 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
9340 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
9341 * anndv,annvd,annds,annsd,
9343 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
9345 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
9348 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
9349 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
9350 * acouzz,acouhh,acouds,acousd,
9351 * acoudz,acouzd,acoudi,
9352 * acoudv,acouvd,acoucc
9354 COMMON /zsea/zseaav,zseasu,anzsea
9355 COMMON /casadi/casaxx,icasad
9360 IF(nchvs1(i).EQ.99.AND.nchvs2(i).EQ.99) go to 50
9364 IF (ipco.GE.6)
WRITE (6,1010) ipvq(is1),ippv1(is1),ippv2(is1),
9365 + itsq(is2),itsaq(is2), amcvs1(i),amcvs2(i),gacvs1(i),gacvs2(i),
9366 + bgxvs1(i),bgyvs1(i),bgzvs1(i), bgxvs2(i),bgyvs2(i),bgzvs2(i),
9367 + nchvs1(i),nchvs2(i),ijcvs1(i),ijcvs2(i), pqvsa1(i,4),pqvsa2
9368 + (i,4),pqvsb1(i,4),pqvsb2(i,4)
9381 ippp = ifrovp(intvs1(i))
9387 WRITE(6,*)
' VS qq-q ,IFB1,IFB2,IFB3,',
9388 *
'INTVS1=IS1,INTVS2=IS2,JIPP,JITTX',
9389 * ifb1,ifb2,ifb3,intvs1(i),intvs2(i),jipp,jittx
9390 WRITE (6,*)
' target sea quark IFB3=',ifb3,
9391 *
' from IS2=',intvs2(i)
9392 WRITE(6,*)
' with ITSQ(IS2),XTSQ(IS2),IFROST(IS2)',
9393 * itsq(is2),xtsq(is2),ifrost(is2)
9396 IF(ifrost(is2).EQ.ifrovt(ii))iii=ii
9399 WRITE (6,*)
' projectile III=',iii
9400 WRITE(6,*)
' corresp. XTVQ(i),XTVD(i),ITVQ(I),ITTV1(I),ITTV2(I)',
9401 * xtvq(iii),xtvd(iii),itvq(iii),ittv1(iii),ittv2(iii)
9408 IF(
rndm(vv).LE.casaxx)
THEN
9409 IF(
rndm(vvv).LE.0.5d0)
THEN
9416 WRITE(6,*)
' Cas VS2 qq-q 1 ,IFB1,IFB2,IFB3,',
9417 *
'INTVS1=IS1,INTVS2=IS2,III',
9418 * ifb1,ifb2,ifb3,intvs1(i),intvs2(i),iii
9419 * ,
'-----------------------------------------------------'
9428 WRITE(6,*)
' Cas VS2 qq-q 2 ,IFB1,IFB2,IFB3,',
9429 *
'INTVS1=IS1,INTVS2=IS2,III',
9430 * ifb1,ifb2,ifb3,intvs1(i),intvs2(i),iii
9431 * ,
'-----------------------------------------------------'
9442 1010
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
9448 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
9464 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
9466 * ,xpsu(248),xtsu(248)
9467 * ,xpsut(248),xtsut(248)
9468 common/popcck/pdbck,pdbse,pdbseu,
9469 * ijpock,irejck,ick4,ihad4,ick6,ihad6
9470 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
9471 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
9472 *isea43,isea63,irejao
9474 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
9475 +ixpv,ixps,ixtv,ixts, intvv1(248),
9476 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
9478 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
9492 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
9498 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
9500 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
9501 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
9508 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
9511 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
9516 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
9604 COMMON /abrsv/ amcsv1(248),amcsv2(248),gacsv1(248),gacsv2(248),
9605 +bgxsv1(248),bgysv1(248),bgzsv1(248), bgxsv2(248),bgysv2(248),
9606 +bgzsv2(248), nchsv1(248),nchsv2(248),ijcsv1(248),ijcsv2(248),
9607 +pqsva1(248,4),pqsva2(248,4), pqsvb1(248,4),pqsvb2(248,4)
9610 parameter(nfimax=249)
9611 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
9612 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
9613 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
9616 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
9618 COMMON /projk/ iprojk
9620 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
9623 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
9624 * anndv,annvd,annds,annsd,
9626 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
9628 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
9631 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
9632 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
9633 * acouzz,acouhh,acouds,acousd,
9634 * acoudz,acouzd,acoudi,
9635 * acoudv,acouvd,acoucc
9637 COMMON /zsea/zseaav,zseasu,anzsea
9638 COMMON /casadi/casaxx,icasad
9645 IF(nchsv1(i).EQ.99.AND.nchsv2(i).EQ.99) go to 50
9649 IF (ipco.GE.6)
WRITE (6,1000) ipsq(is1),ipsaq(is1),itvq(is2),
9650 + ittv1(is2),ittv2(is2), amcsv1(i),amcsv2(i),gacsv1(i),gacsv2(i),
9651 + bgxsv1(i),bgysv1(i),bgzsv1(i), bgxsv2(i),bgysv2(i),bgzsv2(i),
9652 + nchsv1(i),nchsv2(i),ijcsv1(i),ijcsv2(i), pqsva1(i,4),pqsva2
9653 + (i,4),pqsvb1(i,4),pqsvb2(i,4)
9654 1000
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
9667 ittt = ifrovt(intsv2(i))
9673 WRITE(6,*)
' SV q-qq ,IFB1,IFB2,IFB3,',
9674 *
'INTSV1=IS1,INTSV2=IS2,JIPPX,JITT',
9675 * ifb1,ifb2,ifb3,intsv1(i),intsv2(i),jippx,jitt
9676 WRITE (6,*)
' projectile sea quark IFB1=',ifb1,
9677 *
' from IS1=',intsv1(i)
9678 WRITE(6,*)
' with IPSQ(IS1),XPSQ(IS1),IFROSP(IS1)',
9679 * ipsq(is1),xpsq(is1),ifrosp(is1)
9682 IF(ifrosp(is1).EQ.ifrovp(ii))iii=ii
9685 WRITE (6,*)
' projectile III=',iii
9686 WRITE(6,*)
' corresp. XPVQ(i),XPVD(i),IPVQ(I),IPPV1(I),IPPV2(I)',
9687 * xpvq(iii),xpvd(iii),ipvq(iii),ippv1(iii),ippv2(iii)
9694 IF(
rndm(vv).LE.casaxx)
THEN
9695 IF(
rndm(vvv).LE.0.5d0)
THEN
9702 WRITE(6,*)
' Cas SV1 q-qq 1 ,IFB1,IFB2,IFB3,',
9703 *
'INTSV1=IS1,INTSV2=IS2,JIPPX,JITT,III',
9704 * ifb1,ifb2,ifb3,intsv1(i),intsv2(i),jippx,jitt,iii
9705 * ,
'-----------------------------------------------------'
9714 WRITE(6,*)
' Cas SV1 q-qq 2 ,IFB1,IFB2,IFB3,',
9715 *
'INTSV1=IS1,INTSV2=IS2,JIPPX,JITT,III',
9716 * ifb1,ifb2,ifb3,intsv1(i),intsv2(i),jippx,jitt,iii
9717 * ,
'-----------------------------------------------------'
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