1 SUBROUTINE diqsv(ECM,ITV,J,IREJ)
2 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6 COMMON /zsea/zseaav,zseasu,anzsea
7 common/popcck/pdbck,pdbse,pdbseu,
8 * ijpock,irejck,ick4,ihad4,ick6,ihad6
9 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
10 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
15 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
21 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
23 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
26 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
30 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
32 * ,xpsu(248),xtsu(248)
33 * ,xpsut(248),xtsut(248)
35 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
36 +ixpv,ixps,ixtv,ixts, intvv1(248),
37 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
39 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
53 COMMON /abrdv/ amcdv1(248),amcdv2(248),gacdv1(248),gacdv2(248),
54 +bgxdv1(248),bgydv1(248),bgzdv1(248), bgxdv2(248),bgydv2(248),
55 +bgzdv2(248), nchdv1(248),nchdv2(248),ijcdv1(248),ijcdv2(248),
56 +pqdva1(248,4),pqdva2(248,4), pqdvb1(248,4),pqdvb2(248,4)
58 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
60 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
62 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
63 *idzre(3),izdre(3),idiqrz(7)
64 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
65 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
66 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
67 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
68 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
69 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
92 betcha=betoo+1.3d0-log10(ecm)
111 WRITE(6,4567)pc,betcha,pu1,ps1,seasq
112 4567
FORMAT(
' Charm chain ends DIQSV: PC,BETCHA,PU,PS,SEASQ ',5f10.5)
118 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' diqsv'
119 ipsq2(j)=1.d0+
rndm(v1)*(2.d0+2.d0*seasq)
121 IF(rr.LT.pc)ipsq2(j)=4
136 xall=xdfree+xpsq(j)+xpsaq(j)-2.*xdthr
148 xpvd(iitop)=xdthr+dx1
153 amdvq1=xpsq(j)*xtvq(itv)*ecm**2
154 amdvq2=xpsaq(j)*xtvd(itv)*ecm**2
155 idiqre(1)=idiqre(1)+1
156 IF(ipsq(j).GE.3.AND.ipsq2(j).GE.3)
THEN
157 idiqre(2)=idiqre(2)+1
159 IF(amdvq2.LE.17.0d0.OR.amdvq1.LE.6.60d0)
THEN
161 idiqre(3)=idiqre(3)+1
162 idiqre(2)=idiqre(2)-1
163 idiqre(1)=idiqre(1)-1
169 ELSEIF(ipsq(j).GE.3.OR.ipsq2(j).GE.3)
THEN
170 idiqre(4)=idiqre(4)+1
172 IF(amdvq2.LE.13.6d0.OR.amdvq1.LE.5.80d0)
THEN
174 idiqre(5)=idiqre(5)+1
175 idiqre(4)=idiqre(4)-1
176 idiqre(1)=idiqre(1)-1
183 idiqre(6)=idiqre(6)+1
185 IF(amdvq2.LE.12.40d0.OR.amdvq1.LE.3.9d0)
THEN
187 idiqre(7)=idiqre(7)+1
188 idiqre(6)=idiqre(6)-1
189 idiqre(1)=idiqre(1)-1
209 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
213 COMMON /zsea/zseaav,zseasu,anzsea
214 common/popcck/pdbck,pdbse,pdbseu,
215 * ijpock,irejck,ick4,ihad4,ick6,ihad6
216 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
217 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
218 *isea43,isea63,irejao
227 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
315 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
321 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
323 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
326 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
328 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
329 +ixpv,ixps,ixtv,ixts, intvv1(248),
330 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
332 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
346 COMMON /abrdv/ amcdv1(248),amcdv2(248),gacdv1(248),gacdv2(248),
347 +bgxdv1(248),bgydv1(248),bgzdv1(248), bgxdv2(248),bgydv2(248),
348 +bgzdv2(248), nchdv1(248),nchdv2(248),ijcdv1(248),ijcdv2(248),
349 +pqdva1(248,4),pqdva2(248,4), pqdvb1(248,4),pqdvb2(248,4)
353 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
355 * ,xpsu(248),xtsu(248)
356 * ,xpsut(248),xtsut(248)
358 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
359 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
361 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
362 *idzre(3),izdre(3),idiqrz(7)
368 COMMON /trafop/ gamp,bgamp,betp
370 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
371 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
372 +prebin,taebin,fermod,etacou
374 COMMON /fermi/ pquar(4,248),paquar(4,248), tquar(4,248),taquar
387 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
388 +iibar(210),k1(210),k2(210)
391 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
393 COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
394 +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
395 +irvs14, irvv11,irvv12,irvv13,irvv14
397 COMMON /projk/ iprojk
398 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
399 common/rptshm/rproj,rtarg,bimpac
400 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
401 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
402 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
403 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
404 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
405 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
408 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' kkevdv'
412 IF(nchdv1(
n).EQ.99.AND.nchdv2(
n).EQ.99)go to 10
419 psqpx=xpsq(ixspr)*prmom(1,inucpr)
420 psqpy=xpsq(ixspr)*prmom(2,inucpr)
421 psqpz=xpsq(ixspr)*prmom(3,inucpr)
422 psqe=xpsq(ixspr)*prmom(4,inucpr)
423 psaqpx=xpsaq(ixspr)*prmom(1,inucpr)
424 psaqpy=xpsaq(ixspr)*prmom(2,inucpr)
425 psaqpz=xpsaq(ixspr)*prmom(3,inucpr)
426 psaqe=xpsaq(ixspr)*prmom(4,inucpr)
433 tvqpx=xtvq(ixvta)*tamom(1,inucta)
434 tvqpy=xtvq(ixvta)*tamom(2,inucta)
435 tvqpz=xtvq(ixvta)*tamom(3,inucta)
436 tvqe=xtvq(ixvta)*tamom(4,inucta)
437 tvdqpx=xtvd(ixvta)*tamom(1,inucta)
438 tvdqpy=xtvd(ixvta)*tamom(2,inucta)
439 tvdqpz=xtvd(ixvta)*tamom(3,inucta)
440 tvdqe=xtvd(ixvta)*tamom(4,inucta)
449 rtix=(vhkk(1,itnu))*1.e12-bimpac*0.1
450 rtiy=vhkk(2,itnu)*1.e12
451 rtiz=vhkk(3,itnu)*1.e12
452 CALL
cromsc(psqpx,psqpy,psqpz,psqe,rtix,rtiy,rtiz,
453 * psqnx,psqny,psqnz,psqne,51)
458 CALL
cromsc(psaqpx,psaqpy,psaqpz,psaqe,rtix,rtiy,rtiz,
459 * psaqnx,psaqny,psaqnz,psaqne,52)
470 rtix=(vhkk(1,itnu))*1.e12-bimpac*0.1
471 rtiy=vhkk(2,itnu)*1.e12
472 rtiz=vhkk(3,itnu)*1.e12
473 CALL
cromsc(tvqpx,tvqpy,tvqpz,tvqe,rtix,rtiy,rtiz,
474 * tvqnx,tvqny,tvqnz,tvqne,53)
479 CALL
cromsc(tvdqpx,tvdqpy,tvdqpz,tvdqe,rtix,rtiy,rtiz,
480 * tvdqnx,tvdqny,tvdqnz,tvdqne,54)
490 IF(ip.GE.0)go to 1779
491 psqpz2=psqe**2-psqpx**2-psqpy**2
500 psaqp2=psaqe**2-psaqpx**2-psaqpy**2
509 tvqpz2=tvqe**2-tvqpx**2-tvqpy**2
518 tdqpz2=tvdqe**2-tvdqpx**2-tvdqpy**2
564 WRITE(6,
'(A,I5)')
' KKEVDV - IRDV13=',irdv13
565 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
566 +
' DV: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
567 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
568 + amch1,amch2,irej,ikvala,pttq1,ptta1
574 IF(nselpt.EQ.1)CALL
selpt( ptxsq1,ptysq1,plq1,
575 + eq1,ptxsa1,ptysa1,plaq1,eaq1,
576 + ptxsa2,ptysa2,plaq2,eaq2,
577 + ptxsq2,ptysq2,plq2,eq2,
578 + amch1,amch2,irej,ikvala,pttq1,ptta1,
581 IF(nselpt.EQ.0)CALL
selpt4( ptxsq1,ptysq1,plq1,
582 + eq1,ptxsa1,ptysa1,plaq1,eaq1,
583 + ptxsa2,ptysa2,plaq2,eaq2,
584 + ptxsq2,ptysq2,plq2,eq2,
585 + amch1,amch2,irej,ikvala,pttq1,ptta1,nselpt)
587 WRITE(6,
'(A,I5)')
' KKEVDV - IRDV13=',irdv13
588 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
589 +
' DV: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
590 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
591 + amch1,amch2,irej,ikvala,pttq1,ptta1
593 IF (ipev.GE.7)
WRITE(6,
'(A,I10)')
598 WRITE(6,
'(A,I5)')
' KKEVDV - IRDV13=',irdv13
599 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
600 +
' DV: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
601 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
602 + amch1,amch2,irej,ikvala,pttq1,ptta1
609 ptxch1=ptxsq1 + ptxsq2
610 ptych1=ptysq1 + ptysq2
613 ptxch2=ptxsa2 + ptxsa1
614 ptych2=ptysa2 + ptysa1
617 ammm=
sqrt((ech1+ech2)**2-(ptxch1+ptxch2)**2
618 + -(ptych1+ptych2)**2-(ptzch1+ptzch2)**2)
621 IF (ipev.GE.6)
WRITE(6,
'(A,I10/A,5F12.5/A,5F12.5)')
622 +
' DV: IREJ ',irej,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
623 + amch1,ptxch1,ptych1,ptzch1,ech1,
624 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
632 CALL
cobcma(ipsq(ixspr),ipsq2(ixspr),itvq(ixvta), ijnch1,nnch1,
633 + irej,amch1,amch1n,1)
638 WRITE(6,
'(A,I5)')
' KKEVDV - IRDV11=',irdv11
639 WRITE(6,
'(A,6I5/6E12.4/2E12.4)')
' DV:', ipsq(ixspr),ittv1
640 + (ixvta),ittv2(ixvta),ijnch1,nnch1,irej, xpsq(ixspr),xpsaq
641 + (ixspr),xpsqcm,xpsacm, xtvq(ixvta),xtvd(ixvta),amch1,amch1n
649 CALL
cormom(amch1,amch2,amch1n,amch2n,
650 + ptxsq1,ptysq1,plq1,eq1,
651 + ptxsa1,ptysa1,plaq1,eaq1,
652 + ptxsa2,ptysa2,plaq2,eaq2,
653 + ptxsq2,ptysq2,plq2,eq2,
654 + ptxch1,ptych1,ptzch1,ech1, ptxch2,ptych2,ptzch2,ech2,
659 IF (ipev.GE.2)
WRITE(6,
'(A,I10/A,5F12.5/A,5F12.5)')
660 +
' DV(2): IREJ ',irej,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
661 + amch1,ptxch1,ptych1,ptzch1,ech1,
662 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
665 IF(ipev.GE.1)
WRITE(6,
'(A)')
' dv cormom rej.'
668 IF (amch2 .LT.3.)
THEN
669 IF(ipev.GE.1)
WRITE(6,
'(A,F10.2)')
' dv amch2',amch2
714 ihkkpd=jhkkps(ixspr )
715 ihkkpo=jhkkps(ixspr )-1
716 ihkktd=jhkktv(ixvta )
717 ihkkto=jhkktv(ixvta )-1
718 IF (ipev.GT.3)
WRITE(6,1000)ixspr,inucpr,jnucpr,ihkkpo,ihkkpd
719 1000
FORMAT (
' IXSPR,INUCPR,JNUCPR,IHKKPO,IHKKPD ',5i5)
720 IF (ipev.GT.3)
WRITE(6,1010)ixvta,inucta,jnucta,ihkkto,ihkktd
721 1010
FORMAT (
' IXVTA,INUCTA,JNUCTA,IHKKTO,IHKKTD ',5i5)
725 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
730 idhkk(ihkk)=idhkk(ihkkpo)
731 jmohkk(1,ihkk)=ihkkpo
732 jmohkk(2,ihkk)=jmohkk(1,ihkkpo)
733 jdahkk(1,ihkk)=ihkk+2
734 jdahkk(2,ihkk)=ihkk+2
735 phkk(1,ihkk)=pqdva1(
n,1)
736 phkk(2,ihkk)=pqdva1(
n,2)
737 phkk(3,ihkk)=pqdva1(
n,3)
738 phkk(4,ihkk)=pqdva1(
n,4)
742 vhkk(1,ihkk)=vhkk(1,ihkkpo)+xxpp
743 vhkk(2,ihkk)=vhkk(2,ihkkpo)+yypp
744 vhkk(3,ihkk)=vhkk(3,ihkkpo)
745 vhkk(4,ihkk)=vhkk(4,ihkkpo)
746 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
747 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
748 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
750 1020
FORMAT (i6,i4,5i6,9e10.2)
754 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
759 idhkk(ihkk)=idhkk(ihkktd)
760 jmohkk(1,ihkk)=ihkktd
761 jmohkk(2,ihkk)=jmohkk(1,ihkktd)
762 jdahkk(1,ihkk)=ihkk+1
763 jdahkk(2,ihkk)=ihkk+1
764 phkk(1,ihkk)=pqdva2(
n,1)
765 phkk(2,ihkk)=pqdva2(
n,2)
766 phkk(3,ihkk)=pqdva2(
n,3)
767 phkk(4,ihkk)=pqdva2(
n,4)
771 vhkk(1,ihkk)=vhkk(1,ihkktd)+xxpp
772 vhkk(2,ihkk)=vhkk(2,ihkktd)+yypp
773 vhkk(3,ihkk)=vhkk(3,ihkktd)
774 vhkk(4,ihkk)=vhkk(4,ihkktd)
775 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
776 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
777 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
783 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
788 idhkk(ihkk)=88888+nnch1
789 jmohkk(1,ihkk)=ihkk-2
790 jmohkk(2,ihkk)=ihkk-1
801 vhkk(1,nhkk)= vhkk(1,nhkk-1)
802 vhkk(2,nhkk)= vhkk(2,nhkk-1)
803 vhkk(3,nhkk)= vhkk(3,nhkk-1)
804 vhkk(4,nhkk)=vhkk(3,nhkk)/betp-vhkk(3,nhkk-2)/bgamp
807 whkk(1,nhkk)= vhkk(1,nhkk-2)
808 whkk(2,nhkk)= vhkk(2,nhkk-2)
809 whkk(3,nhkk)= vhkk(3,nhkk-2)
810 whkk(4,nhkk)=vhkk(4,nhkk)*gamp-vhkk(3,nhkk)*bgamp
811 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
812 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
813 + (phkk(khkk,ihkk),khkk=1,5), (whkk(khkk,ihkk),khkk=1,4)
816 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
817 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
818 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
825 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
830 idhkk(ihkk)=idhkk(ihkkpd)
831 jmohkk(1,ihkk)=ihkkpd
832 jmohkk(2,ihkk)=jmohkk(1,ihkkpd)
833 jdahkk(1,ihkk)=ihkk+2
834 jdahkk(2,ihkk)=ihkk+2
835 phkk(1,ihkk)=pqdvb1(
n,1)
836 phkk(2,ihkk)=pqdvb1(
n,2)
837 phkk(3,ihkk)=pqdvb1(
n,3)
838 phkk(4,ihkk)=pqdvb1(
n,4)
842 vhkk(1,ihkk)=vhkk(1,ihkkpd)+xxpp
843 vhkk(2,ihkk)=vhkk(2,ihkkpd)+yypp
844 vhkk(3,ihkk)=vhkk(3,ihkkpd)
845 vhkk(4,ihkk)=vhkk(4,ihkkpd)
846 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
847 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
848 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
853 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
858 idhkk(ihkk)=idhkk(ihkkto)
859 jmohkk(1,ihkk)=ihkkto
860 jmohkk(2,ihkk)=jmohkk(1,ihkkto)
861 jdahkk(1,ihkk)=ihkk+1
862 jdahkk(2,ihkk)=ihkk+1
863 phkk(1,ihkk)=pqdvb2(
n,1)
864 phkk(2,ihkk)=pqdvb2(
n,2)
865 phkk(3,ihkk)=pqdvb2(
n,3)
866 phkk(4,ihkk)=pqdvb2(
n,4)
870 vhkk(1,ihkk)=vhkk(1,ihkkto)+xxpp
871 vhkk(2,ihkk)=vhkk(2,ihkkto)+yypp
872 vhkk(3,ihkk)=vhkk(3,ihkkto)
873 vhkk(4,ihkk)=vhkk(4,ihkkto)
874 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
875 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
876 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
882 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
887 idhkk(ihkk)=88888+nnch2
888 jmohkk(1,ihkk)=ihkk-2
889 jmohkk(2,ihkk)=ihkk-1
900 vhkk(1,nhkk)= vhkk(1,nhkk-1)
901 vhkk(2,nhkk)= vhkk(2,nhkk-1)
902 vhkk(3,nhkk)= vhkk(3,nhkk-1)
903 vhkk(4,nhkk)=vhkk(3,nhkk)/betp-vhkk(3,nhkk-2)/bgamp
906 whkk(1,nhkk)= vhkk(1,nhkk-2)
907 whkk(2,nhkk)= vhkk(2,nhkk-2)
908 whkk(3,nhkk)= vhkk(3,nhkk-2)
909 whkk(4,nhkk)=vhkk(4,nhkk)*gamp-vhkk(3,nhkk)*bgamp
910 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
911 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
912 + (phkk(khkk,ihkk),khkk=1,5), (whkk(khkk,ihkk),khkk=1,4)
915 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
916 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
917 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
927 gacdv1(
n)=qech1/amch1
928 bgxdv1(
n)=qtxch1/amch1
929 bgydv1(
n)=qtych1/amch1
930 bgzdv1(
n)=qtzch1/amch1
931 gacdv2(
n)=qech2/amch2
932 bgxdv2(
n)=qtxch2/amch2
933 bgydv2(
n)=qtych2/amch2
934 bgzdv2(
n)=qtzch2/amch2
939 IF (ipev.GE.6)
WRITE(6,
'(A/I10,4F12.7,5I5/10X,4F12.6/10X,6F12.6,
940 +4I5/8F15.5/ 8F15.5)')
' DV / FINAL PRINT',
n
957 IF(issqq.EQ.3.AND.jssqq.EQ.3)
THEN
959 ELSEIF(issqq.EQ.3.OR.jssqq.EQ.3)
THEN
972 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
983 COMMON /zsea/zseaav,zseasu,anzsea
984 common/popcck/pdbck,pdbse,pdbseu,
985 * ijpock,irejck,ick4,ihad4,ick6,ihad6
986 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
987 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
988 *isea43,isea63,irejao
992 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
998 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
1000 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
1003 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
1007 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
1009 * ,xpsu(248),xtsu(248)
1010 * ,xpsut(248),xtsut(248)
1012 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
1013 +ixpv,ixps,ixtv,ixts, intvv1(248),
1014 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
1016 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
1030 COMMON /abrdv/ amcdv1(248),amcdv2(248),gacdv1(248),gacdv2(248),
1031 +bgxdv1(248),bgydv1(248),bgzdv1(248), bgxdv2(248),bgydv2(248),
1032 +bgzdv2(248), nchdv1(248),nchdv2(248),ijcdv1(248),ijcdv2(248),
1033 +pqdva1(248,4),pqdva2(248,4), pqdvb1(248,4),pqdvb2(248,4)
1035 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
1036 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
1046 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
1135 parameter(nfimax=249)
1136 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
1137 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
1138 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
1141 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
1143 COMMON /projk/ iprojk
1145 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
1147 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
1148 * anndv,annvd,annds,annsd,
1150 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
1152 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
1155 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
1156 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
1157 * acouzz,acouhh,acouds,acousd,
1158 * acoudz,acouzd,acoudi,
1159 * acoudv,acouvd,acoucc
1160 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
1161 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
1162 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
1163 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
1164 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
1165 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
1167 dimension poj(4),pat(4)
1169 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' hadrdv'
1174 IF(nchdv1(i).EQ.99.AND.nchdv2(i).EQ.99) go to 50
1177 IF(ipco.GE.3)
WRITE(6,*)
' hadrdv I IS1,IS2 ',i,is1,is2
1179 IF (ipco.GE.6)
WRITE (6,1000) ipsq(is1),ipsaq(is1),itvq(is2),
1180 + ittv1(is2),ittv2(is2), amcdv1(i),amcdv2(i),gacdv1(i),gacdv2(i),
1181 + bgxdv1(i),bgydv1(i),bgzdv1(i), bgxdv2(i),bgydv2(i),bgzdv2(i),
1182 + nchdv1(i),nchdv2(i),ijcdv1(i),ijcdv2(i), pqdva1(i,4),pqdva2
1183 + (i,4),pqdvb1(i,4),pqdvb2(i,4)
1184 1000
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
1194 IF((nchdv1(i).NE.0.OR.nchdv2(i).NE.0).AND.ip.NE.1)
1195 & CALL
saptre(amcdv1(i),gacdv1(i),bgxdv1(i),bgydv1(i),bgzdv1(i),
1196 & amcdv2(i),gacdv2(i),bgxdv2(i),bgydv2(i),bgzdv2(i))
1206 IF(ipco.GE.3)
WRITE(6,*)
' INTDV1(I) ',intdv1(i)
1209 IF(intdv1(i).GE.1)
THEN
1210 ipppp = ifrosp(intdv1(i))
1214 IF(ipco.GE.3)
WRITE(6,*)
' IPPP,IPPPP ',ippp,ipppp
1221 IF(ipco.GE.3)
WRITE(6,*)
' JIPP ',jipp
1223 IF(ipco.GE.3)
WRITE(6,
'(A,3I5)')
'HADRVS: I,IPPP,JIPP ',
1229 IF(ifb1.LE.2.AND.ifb2.LE.2)
THEN
1231 ELSEIF((ifb1.EQ.3.AND.ifb2.LE.2).OR.
1232 * (ifb2.EQ.3.AND.ifb1.LE.2))
THEN
1234 ELSEIF(ifb1.EQ.3.AND.ifb2.EQ.3)
THEN
1238 * CALL
hadjet(nhad,amcdv1(i),poj,pat,gacdv1(i),
1240 + (i),bgzdv1(i),ifb1,ifb2,ifb3,ifb4, ijcdv1(i),
1241 * ijcdv1(i),6,nchdv1
1244 zseawu=
rndm(bb)*2.d0*zseaav
1245 rseack=float(jipp)*pdbse +zseawu*pdbseu
1246 IF(ipco.GE.1)
WRITE(6,*)
'HADJSE JIPP,RSEACK,PDBSE 1 dpmnuc5',
1249 IF(nchdv1(i).EQ.0)
THEN
1252 IF(
rndm(v).LE.rseack)
THEN
1254 IF(amcdv1(i).GT.2.3d0)
THEN
1256 CALL
hadjse(nhad,amcdv1(i),poj,pat,gacdv1(i),bgxdv1(i),
1258 + (i),bgzdv1(i),ifb1,ifb2,ifb3,ifb4, ijcdv1(i),ijcdv1(i),6,
1260 + (i),6,irejss,iissqq)
1261 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JIPP,',
1262 *
'RSEACK,IREJSS 1 dpmnuc5 ',
1263 + jipp,rseack,irejss
1266 IF(irejss.EQ.1)irejse=irejse+1
1267 IF(irejss.EQ.3)irejs3=irejs3+1
1268 IF(irejss.EQ.2)irejs0=irejs0+1
1269 CALL
hadjet(nhad,amcdv1(i),poj,pat,gacdv1(i),
1271 + (i),bgzdv1(i),ifb1,ifb2,ifb3,ifb4, ijcdv1(i),
1272 * ijcdv1(i),6,nchdv1
1284 CALL
hadjet(nhad,amcdv1(i),poj,pat,gacdv1(i),
1286 + (i),bgzdv1(i),ifb1,ifb2,ifb3,ifb4, ijcdv1(i),
1287 * ijcdv1(i),6,nchdv1
1297 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRDV: NHKKNMXHKK ',nhkk,
nmxhkk
1302 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
1306 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
1307 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
1316 ptdv=ptdv+
sqrt(pxf(j)**2+pyf(j)**2)
1319 IF(ibarf(j).EQ.500)istist=2
1321 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),13)
1322 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
1324 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
1325 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
1326 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
1345 IF(ifb1.LE.8.AND.ifb2.LE.8)
THEN
1347 ELSEIF((ifb1.EQ.9.AND.ifb2.LE.8).OR.
1348 * (ifb2.EQ.9.AND.ifb1.LE.8))
THEN
1350 ELSEIF(ifb1.EQ.9.AND.ifb2.EQ.9)
THEN
1353 CALL
hadjet(nhad,amcdv2(i),poj,pat,gacdv2(i),
1355 + (i),bgzdv2(i),ifb1,ifb2,ifb3,ifb4, ijcdv2(i),
1356 * ijcdv2(i),5,nchdv2
1363 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRDV: NHKKNMXHKK ', nhkk,
1369 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
1373 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
1374 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
1383 ptdv=ptdv+
sqrt(pxf(j)**2+pyf(j)**2)
1386 IF(ibarf(j).EQ.500)istist=2
1388 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),14)
1389 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
1391 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
1392 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
1393 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
1404 1010
FORMAT (i6,i4,5i6,9e10.2)
1405 1020
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
1406 1030
FORMAT (
' NHKK,IDHKK(NHKK) ',3i10)
1412 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1416 COMMON /zsea/zseaav,zseasu,anzsea
1417 common/popcck/pdbck,pdbse,pdbseu,
1418 * ijpock,irejck,ick4,ihad4,ick6,ihad6
1419 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
1420 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
1421 *isea43,isea63,irejao
1425 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
1431 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
1433 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
1436 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
1440 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
1442 * ,xpsu(248),xtsu(248)
1443 * ,xpsut(248),xtsut(248)
1445 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
1446 +ixpv,ixps,ixtv,ixts, intvv1(248),
1447 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
1449 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
1463 COMMON /abrvd/ amcvd1(248),amcvd2(248),gacvd1(248),gacvd2(248),
1464 +bgxvd1(248),bgyvd1(248),bgzvd1(248), bgxvd2(248),bgyvd2(248),
1465 +bgzvd2(248), nchvd1(248),nchvd2(248),ijcvd1(248),ijcvd2(248),
1466 +pqvda1(248,4),pqvda2(248,4), pqvdb1(248,4),pqvdb2(248,4)
1468 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
1470 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
1472 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
1473 *idzre(3),izdre(3),idiqrz(7)
1474 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
1475 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
1476 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
1477 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
1478 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
1479 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
1482 parameter(ummm=0.3d0)
1483 parameter(smmm=0.5d0)
1484 parameter(cmmm=1.3d0)
1502 betcha=betoo+1.3d0-log10(ecm)
1521 WRITE(6,4567)pc,betcha,pu1,ps1,seasq
1522 4567
FORMAT(
' Charm chain ends DIQVS: PC,BETCHA,PU,PS,SEASQ ',5f10.5)
1525 IF(iphkk.GE.6)
WRITE (6,
'( A,2I10)')
' diqvs IPV,J ',ipv,j
1529 itsq2(j)=1.d0+
rndm(v1)*(2.d0+2.d0*seasq)
1531 IF(rr.LT.pc)itsq2(j)=4
1546 xall=xdfree+xtsq(j)+xtsaq(j)-2.*xdthr
1558 xtvd(iitot)=xdthr+dx1
1564 amvdq1=xtsq(j)*xpvq(ipv)*ecm**2
1565 amvdq2=xtsaq(j)*xpvd(ipv)*ecm**2
1566 idiqre(1)=idiqre(1)+1
1567 IF(itsq(j).GE.3.AND.itsq2(j).GE.3)
THEN
1568 idiqre(2)=idiqre(2)+1
1570 IF(amvdq2.LE.17.0d0.OR.amvdq1.LE.6.60d0)
THEN
1572 idiqre(3)=idiqre(3)+1
1573 idiqre(2)=idiqre(2)-1
1574 idiqre(1)=idiqre(1)-1
1580 ELSEIF(itsq(j).GE.3.OR.itsq2(j).GE.3)
THEN
1581 idiqre(4)=idiqre(4)+1
1583 IF(amvdq2.LE.13.6.OR.amvdq1.LE.5.80)
THEN
1585 idiqre(5)=idiqre(5)+1
1586 idiqre(4)=idiqre(4)-1
1587 idiqre(1)=idiqre(1)-1
1594 idiqre(6)=idiqre(6)+1
1596 IF(amvdq2.LE.12.40d0.OR.amvdq1.LE.3.9d0)
THEN
1598 idiqre(7)=idiqre(7)+1
1599 idiqre(6)=idiqre(6)-1
1600 idiqre(1)=idiqre(1)-1
1623 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1628 COMMON /zsea/zseaav,zseasu,anzsea
1629 common/popcck/pdbck,pdbse,pdbseu,
1630 * ijpock,irejck,ick4,ihad4,ick6,ihad6
1631 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
1632 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
1633 *isea43,isea63,irejao
1640 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
1728 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
1734 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
1736 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
1739 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
1741 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
1742 +ixpv,ixps,ixtv,ixts, intvv1(248),
1743 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
1745 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
1759 COMMON /abrvd/ amcvd1(248),amcvd2(248),gacvd1(248),gacvd2(248),
1760 +bgxvd1(248),bgyvd1(248),bgzvd1(248), bgxvd2(248),bgyvd2(248),
1761 +bgzvd2(248), nchvd1(248),nchvd2(248),ijcvd1(248),ijcvd2(248),
1762 +pqvda1(248,4),pqvda2(248,4), pqvdb1(248,4),pqvdb2(248,4)
1766 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
1768 * ,xpsu(248),xtsu(248)
1769 * ,xpsut(248),xtsut(248)
1770 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
1771 *idzre(3),izdre(3),idiqrz(7)
1773 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
1774 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
1781 COMMON /trafop/ gamp,bgamp,betp
1783 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
1784 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
1785 +prebin,taebin,fermod,etacou
1787 COMMON /fermi/ pquar(4,248),paquar(4,248), tquar(4,248),taquar
1800 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
1801 +iibar(210),k1(210),k2(210)
1804 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
1806 COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
1807 +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
1808 +irvs14, irvv11,irvv12,irvv13,irvv14
1810 COMMON /projk/ iprojk
1811 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
1812 common/rptshm/rproj,rtarg,bimpac
1813 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
1814 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
1815 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
1816 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
1817 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
1818 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
1821 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' kkevvd'
1825 IF(nchvd1(
n).EQ.99.AND.nchvd2(
n).EQ.99)go to 10
1829 inucpr=ifrovp(ixvpr)
1830 jnucpr=itovp(inucpr)
1832 pvqpx=xpvq(ixvpr)*prmom(1,inucpr)
1833 pvqpy=xpvq(ixvpr)*prmom(2,inucpr)
1834 pvqpz=xpvq(ixvpr)*prmom(3,inucpr)
1835 pvqe=xpvq(ixvpr)*prmom(4,inucpr)
1836 pvdqpx=xpvd(ixvpr)*prmom(1,inucpr)
1837 pvdqpy=xpvd(ixvpr)*prmom(2,inucpr)
1838 pvdqpz=xpvd(ixvpr)*prmom(3,inucpr)
1839 pvdqe=xpvd(ixvpr)*prmom(4,inucpr)
1843 inucta=ifrost(ixsta)
1844 jnucta=itovt(inucta)
1846 tsqpx=xtsq(ixsta)*tamom(1,inucta)
1847 tsqpy=xtsq(ixsta)*tamom(2,inucta)
1848 tsqpz=xtsq(ixsta)*tamom(3,inucta)
1849 tsqe=xtsq(ixsta)*tamom(4,inucta)
1850 tsaqpx=xtsaq(ixsta)*tamom(1,inucta)
1851 tsaqpy=xtsaq(ixsta)*tamom(2,inucta)
1852 tsaqpz=xtsaq(ixsta)*tamom(3,inucta)
1853 tsaqe=xtsaq(ixsta)*tamom(4,inucta)
1862 rtix=(vhkk(1,itnu))*1.e12-bimpac*0.1
1863 rtiy=vhkk(2,itnu)*1.e12
1864 rtiz=vhkk(3,itnu)*1.e12
1865 CALL
cromsc(pvqpx,pvqpy,pvqpz,pvqe,rtix,rtiy,rtiz,
1866 * pvqnx,pvqny,pvqnz,pvqne,55)
1871 CALL
cromsc(pvdqpx,pvdqpy,pvdqpz,pvdqe,rtix,rtiy,rtiz,
1872 * pvdqnx,pvdqny,pvdqnz,pvdqne,56)
1883 rtix=(vhkk(1,itnu))*1.e12-bimpac*0.1
1884 rtiy=vhkk(2,itnu)*1.e12
1885 rtiz=vhkk(3,itnu)*1.e12
1886 CALL
cromsc(tsqpx,tsqpy,tsqpz,tsqe,rtix,rtiy,rtiz,
1887 * tsqnx,tsqny,tsqnz,tsqne,57)
1892 CALL
cromsc(tsaqpx,tsaqpy,tsaqpz,tsaqe,rtix,rtiy,rtiz,
1893 * tsaqnx,tsaqny,tsaqnz,tsaqne,58)
1902 IF(ip.GE.0)go to 1779
1903 pvqpz2=pvqe**2-pvqpx**2-pvqpy**2
1904 IF(pvqpz2.GE.0.)
THEN
1912 pdqpz2=pvdqe**2-pvdqpx**2-pvdqpy**2
1913 IF(pdqpz2.GE.0.)
THEN
1921 tsqpz2=tsqe**2-tsqpx**2-tsqpy**2
1922 IF(tsqpz2.GE.0.)
THEN
1930 taqpz2=tsaqe**2-tsaqpx**2-tsaqpy**2
1931 IF(taqpz2.GE.0.)
THEN
1932 tsaqpz=-
sqrt(taqpz2)
1972 WRITE(6,
'(A,I5)')
' KKEVVD - IRVD13=',irvd13
1973 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
1974 +
' VD: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
1975 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
1976 + amch1,amch2,irej,ikvala,pttq1,ptta1
1982 IF(nselpt.EQ.1)CALL
selpt(ptxsq1,ptysq1,plq1,
1983 + eq1,ptxsa1,ptysa1,plaq1,eaq1,
1984 + ptxsa2,ptysa2,plaq2,eaq2,
1985 + ptxsq2,ptysq2,plq2,eq2,
1986 + amch1,amch2,irej,ikvala,pttq1,ptta1,
1989 IF(nselpt.EQ.0)CALL
selpt4(ptxsq1,ptysq1,plq1,
1990 + eq1,ptxsa1,ptysa1,plaq1,eaq1,
1991 + ptxsa2,ptysa2,plaq2,eaq2,
1992 + ptxsq2,ptysq2,plq2,eq2,
1993 + amch1,amch2,irej,ikvala,pttq1,ptta1,nselpt)
1995 WRITE(6,
'(A,I5)')
' KKEVVD - IRVD13=',irvd13
1996 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
1997 +
' VD: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
1998 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
1999 + amch1,amch2,irej,ikvala,pttq1,ptta1
2000 WRITE(6,
'(A,I5)')
' KKEVVD - IRVD13=',irvd13
2001 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
2002 +
' VD: amch1,amch2 ',
2006 IF (ipev.GE.7)
WRITE(6,
'(A/5X,I10)')
2011 WRITE(6,
'(A,I5)')
' KKEVVD - IRVD13=',irvd13
2012 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
2013 +
' VD: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
2014 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
2015 + amch1,amch2,irej,ikvala,pttq1,ptta1
2023 ptxch1=ptxsq1 + ptxsq2
2024 ptych1=ptysq1 + ptysq2
2027 ptxch2=ptxsa2 + ptxsa1
2028 ptych2=ptysa2 + ptysa1
2029 ptzch2=plaq2 + plaq1
2031 ammm=
sqrt((ech1+ech2)**2-(ptxch1+ptxch2)**2
2032 + -(ptych1+ptych2)**2-(ptzch1+ptzch2)**2)
2034 IF (ipev.GE.6)
WRITE(6,
'(A,I10/A,5F12.5/A,5F12.5)')
2035 +
' VD: IREJ ', irej,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
2036 + amch1,ptxch1,ptych1,ptzch1,ech1,
2037 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
2045 CALL
cobcma(itsq(ixsta),itsq2(ixsta),ipvq(ixvpr), ijnch1,nnch1,
2046 + irej,amch1,amch1n,1)
2051 WRITE(6,
'(A,I5)')
' KKEVVD - IRVD11=',irvd11
2052 WRITE(6,
'(A,6I5/6E12.4/2E12.4)')
' VD:', ipvq(ixvpr),itsq
2053 + (ixsta),itsq2(ixsta),ijnch1,nnch1,irej, xpvq(ixvpr),xpvd
2054 + (ixvpr),xpvqcm,xpvdcm,
2055 + xtsq(ixsta),xtsaq(ixsta),amch1,amch1n
2063 CALL
cormom(amch1,amch2,amch1n,amch2n,
2064 + ptxsq1,ptysq1,plq1,eq1,
2065 + ptxsa1,ptysa1,plaq1,eaq1,
2066 + ptxsa2,ptysa2,plaq2,eaq2,
2067 + ptxsq2,ptysq2,plq2,eq2,
2068 + ptxch1,ptych1,ptzch1,ech1, ptxch2,ptych2,ptzch2,ech2,
2073 IF(ipev.GE.1)
WRITE(6,
'(A)')
' vd CORMOM rej.'
2077 IF (ipev.GE.6)
WRITE(6,
'(A,5F12.5,I10/A,5F12.5/A,5F12.5)')
2078 +
' VD(2): AMMM,GAMMM,BGGGX,BGGGY,BGGGZ,IREJ ', ammm,gammm,bgggx,
2079 + bgggy,bgggz,irej,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
2080 + amch1,ptxch1,ptych1,ptzch1,ech1,
2081 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
2124 ihkkpd=jhkkpv(ixvpr )
2125 ihkkpo=jhkkpv(ixvpr )-1
2126 ihkktd=jhkkts(ixsta )
2127 ihkkto=jhkkts(ixsta )-1
2128 IF (ipev.GT.3)
WRITE(6,1000)ixvpr,inucpr,jnucpr,ihkkpo,ihkkpd
2129 1000
FORMAT (
' IXVPR,INUCPR,JNUCPR,IHKKPO,IHKKPD ',5i5)
2130 IF (ipev.GT.3)
WRITE(6,1010)ixsta,inucta,jnucta,ihkkto,ihkktd
2131 1010
FORMAT (
' IXSTA,INUCTA,JNUCTA,IHKKTO,IHKKTD ',5i5)
2135 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
2140 idhkk(ihkk)=idhkk(ihkkpo)
2141 jmohkk(1,ihkk)=ihkkpo
2142 jmohkk(2,ihkk)=jmohkk(1,ihkkpo)
2143 jdahkk(1,ihkk)=ihkk+2
2144 jdahkk(2,ihkk)=ihkk+2
2145 phkk(1,ihkk)=pqvda1(
n,1)
2146 phkk(2,ihkk)=pqvda1(
n,2)
2147 phkk(3,ihkk)=pqvda1(
n,3)
2148 phkk(4,ihkk)=pqvda1(
n,4)
2152 vhkk(1,ihkk)=vhkk(1,ihkkpo)+xxpp
2153 vhkk(2,ihkk)=vhkk(2,ihkkpo)+yypp
2154 vhkk(3,ihkk)=vhkk(3,ihkkpo)
2155 vhkk(4,ihkk)=vhkk(4,ihkkpo)
2156 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
2157 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
2158 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
2160 1020
FORMAT (i6,i4,5i6,9e10.2)
2164 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
2169 idhkk(ihkk)=idhkk(ihkktd)
2170 jmohkk(1,ihkk)=ihkktd
2171 jmohkk(2,ihkk)=jmohkk(1,ihkktd)
2172 jdahkk(1,ihkk)=ihkk+1
2173 jdahkk(2,ihkk)=ihkk+1
2174 phkk(1,ihkk)=pqvda2(
n,1)
2175 phkk(2,ihkk)=pqvda2(
n,2)
2176 phkk(3,ihkk)=pqvda2(
n,3)
2177 phkk(4,ihkk)=pqvda2(
n,4)
2181 vhkk(1,ihkk)=vhkk(1,ihkktd)+xxpp
2182 vhkk(2,ihkk)=vhkk(2,ihkktd)+yypp
2183 vhkk(3,ihkk)=vhkk(3,ihkktd)
2184 vhkk(4,ihkk)=vhkk(4,ihkktd)
2185 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
2186 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
2187 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
2193 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
2198 idhkk(ihkk)=88888+nnch1
2199 jmohkk(1,ihkk)=ihkk-2
2200 jmohkk(2,ihkk)=ihkk-1
2211 vhkk(1,nhkk)= vhkk(1,nhkk-1)
2212 vhkk(2,nhkk)= vhkk(2,nhkk-1)
2213 vhkk(3,nhkk)= vhkk(3,nhkk-1)
2214 vhkk(4,nhkk)=vhkk(3,nhkk)/betp-vhkk(3,nhkk-2)/bgamp
2216 IF (iprojk.EQ.1)
THEN
2217 whkk(1,nhkk)= vhkk(1,nhkk-2)
2218 whkk(2,nhkk)= vhkk(2,nhkk-2)
2219 whkk(3,nhkk)= vhkk(3,nhkk-2)
2220 whkk(4,nhkk)=vhkk(4,nhkk)*gamp-vhkk(3,nhkk)*bgamp
2221 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
2222 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
2223 + (phkk(khkk,ihkk),khkk=1,5), (whkk(khkk,ihkk),khkk=1,4)
2226 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
2227 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
2228 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
2235 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
2240 idhkk(ihkk)=idhkk(ihkkpd)
2241 jmohkk(1,ihkk)=ihkkpd
2242 jmohkk(2,ihkk)=jmohkk(1,ihkkpd)
2243 jdahkk(1,ihkk)=ihkk+2
2244 jdahkk(2,ihkk)=ihkk+2
2245 phkk(1,ihkk)=pqvdb1(
n,1)
2246 phkk(2,ihkk)=pqvdb1(
n,2)
2247 phkk(3,ihkk)=pqvdb1(
n,3)
2248 phkk(4,ihkk)=pqvdb1(
n,4)
2252 vhkk(1,ihkk)=vhkk(1,ihkkpd)+xxpp
2253 vhkk(2,ihkk)=vhkk(2,ihkkpd)+yypp
2254 vhkk(3,ihkk)=vhkk(3,ihkkpd)
2255 vhkk(4,ihkk)=vhkk(4,ihkkpd)
2256 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
2257 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
2258 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
2263 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
2268 idhkk(ihkk)=idhkk(ihkkto)
2269 jmohkk(1,ihkk)=ihkkto
2270 jmohkk(2,ihkk)=jmohkk(1,ihkkto)
2271 jdahkk(1,ihkk)=ihkk+1
2272 jdahkk(2,ihkk)=ihkk+1
2273 phkk(1,ihkk)=pqvdb2(
n,1)
2274 phkk(2,ihkk)=pqvdb2(
n,2)
2275 phkk(3,ihkk)=pqvdb2(
n,3)
2276 phkk(4,ihkk)=pqvdb2(
n,4)
2280 vhkk(1,ihkk)=vhkk(1,ihkkto)+xxpp
2281 vhkk(2,ihkk)=vhkk(2,ihkkto)+yypp
2282 vhkk(3,ihkk)=vhkk(3,ihkkto)
2283 vhkk(4,ihkk)=vhkk(4,ihkkto)
2284 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
2285 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
2286 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
2292 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
2297 idhkk(ihkk)=88888+nnch2
2298 jmohkk(1,ihkk)=ihkk-2
2299 jmohkk(2,ihkk)=ihkk-1
2310 vhkk(1,nhkk)= vhkk(1,nhkk-1)
2311 vhkk(2,nhkk)= vhkk(2,nhkk-1)
2312 vhkk(3,nhkk)= vhkk(3,nhkk-1)
2313 vhkk(4,nhkk)=vhkk(3,nhkk)/betp-vhkk(3,nhkk-2)/bgamp
2315 IF (iprojk.EQ.1)
THEN
2316 whkk(1,nhkk)= vhkk(1,nhkk-2)
2317 whkk(2,nhkk)= vhkk(2,nhkk-2)
2318 whkk(3,nhkk)= vhkk(3,nhkk-2)
2319 whkk(4,nhkk)=vhkk(4,nhkk)*gamp-vhkk(3,nhkk)*bgamp
2320 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
2321 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
2322 + (phkk(khkk,ihkk),khkk=1,5), (whkk(khkk,ihkk),khkk=1,4)
2325 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
2326 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
2327 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
2337 gacvd1(
n)=qech1/amch1
2338 bgxvd1(
n)=qtxch1/amch1
2339 bgyvd1(
n)=qtych1/amch1
2340 bgzvd1(
n)=qtzch1/amch1
2341 gacvd2(
n)=qech2/amch2
2342 bgxvd2(
n)=qtxch2/amch2
2343 bgyvd2(
n)=qtych2/amch2
2344 bgzvd2(
n)=qtzch2/amch2
2349 IF (ipev.GE.2)
WRITE(6,
'(A/I10,4F12.7,5I5/10X,4F12.6/10X,6F12.6,
2350 +4I5/8F15.5/8F15.5/2I5)')
' VD / FINAL PRINT',
n
2368 IF(issqq.EQ.3.AND.jssqq.EQ.3)
THEN
2370 ELSEIF(issqq.EQ.3.OR.jssqq.EQ.3)
THEN
2383 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2394 COMMON /zsea/zseaav,zseasu,anzsea
2395 common/popcck/pdbck,pdbse,pdbseu,
2396 * ijpock,irejck,ick4,ihad4,ick6,ihad6
2397 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
2398 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
2399 *isea43,isea63,irejao
2403 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
2409 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
2411 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
2414 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
2418 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
2420 * ,xpsu(248),xtsu(248)
2421 * ,xpsut(248),xtsut(248)
2423 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
2424 +ixpv,ixps,ixtv,ixts, intvv1(248),
2425 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
2427 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
2441 COMMON /abrvd/ amcvd1(248),amcvd2(248),gacvd1(248),gacvd2(248),
2442 +bgxvd1(248),bgyvd1(248),bgzvd1(248), bgxvd2(248),bgyvd2(248),
2443 +bgzvd2(248), nchvd1(248),nchvd2(248),ijcvd1(248),ijcvd2(248),
2444 +pqvda1(248,4),pqvda2(248,4), pqvdb1(248,4),pqvdb2(248,4)
2446 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
2447 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
2457 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
2546 parameter(nfimax=249)
2547 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
2548 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
2549 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
2552 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
2554 COMMON /projk/ iprojk
2556 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
2558 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
2559 * anndv,annvd,annds,annsd,
2561 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
2563 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
2566 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
2567 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
2568 * acouzz,acouhh,acouds,acousd,
2569 * acoudz,acouzd,acoudi,
2570 * acoudv,acouvd,acoucc
2571 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
2572 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
2573 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
2574 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
2575 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
2576 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
2578 dimension poj(4),pat(4)
2581 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' hadrvd'
2585 IF(nchvd1(i).EQ.99.AND.nchvd2(i).EQ.99) go to 50
2589 IF (ipco.GE.6)
WRITE (6,1000) ipsq(is1),ipsaq(is1),itvq(is2),
2590 + ittv1(is2),ittv2(is2), amcvd1(i),amcvd2(i),gacvd1(i),gacvd2(i),
2591 + bgxvd1(i),bgyvd1(i),bgzvd1(i), bgxvd2(i),bgyvd2(i),bgzvd2(i),
2592 + nchvd1(i),nchvd2(i),ijcvd1(i),ijcvd2(i), pqvda1(i,4),pqvda2
2593 + (i,4),pqvdb1(i,4),pqvdb2(i,4)
2594 1000
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
2604 IF((nchvd1(i).NE.0.OR.nchvd2(i).NE.0).AND.ip.NE.1)
2605 & CALL
saptre(amcvd1(i),gacvd1(i),bgxvd1(i),
2606 * bgyvd1(i),bgzvd1(i),
2607 & amcvd2(i),gacvd2(i),bgxvd2(i),
2608 * bgyvd2(i),bgzvd2(i))
2637 IF(ifb2.LE.2.AND.ifb3.LE.2)
THEN
2639 ELSEIF((ifb2.EQ.3.AND.ifb3.LE.2).OR.
2640 * (ifb3.EQ.3.AND.ifb2.LE.2))
THEN
2642 ELSEIF(ifb2.EQ.3.AND.ifb3.EQ.3)
THEN
2645 IF((nchvd1(i).NE.0))
THEN
2646 CALL
hadjet(nhad,amcvd1(i),poj,pat,gacvd1(i),
2648 + (i),bgzvd1(i),ifb1,ifb2,ifb3,ifb4, ijcvd1(i),
2649 * ijcvd1(i),4,nchvd1
2653 aack=float(ick4)/float(ick4+ihad4+1)
2654 IF((nchvd1(i).EQ.0))
THEN
2655 zseawu=
rndm(bb)*2.d0*zseaav
2656 rseack=float(jitt)*pdbse +zseawu*pdbseu
2657 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,2F10.3)')
'HADJSE JITT,',
2658 *
'RSEACK,PDBSE 2 dpmnuc5 ',
2661 IF(
rndm(v).LE.rseack)
THEN
2663 IF(amcvd1(i).GT.2.3d0)
THEN
2665 CALL
hadjse(nhad,amcvd1(i),poj,pat,gacvd1(i),
2668 + (i),bgzvd1(i),ifb1,ifb2,ifb3,ifb4, ijcvd1(i),
2671 + (i),3,irejss,iissqq)
2672 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JITT,',
2673 *
'RSEACK,IREJSS 2 dpmnuc5 ',
2674 + jitt,rseack,irejss
2677 IF(irejss.EQ.1)irejse=irejse+1
2678 IF(irejss.EQ.3)irejs3=irejs3+1
2679 IF(irejss.EQ.2)irejs0=irejs0+1
2680 CALL
hadjet(nhad,amcvd1(i),poj,pat,gacvd1(i),
2682 + (i),bgzvd1(i),ifb1,ifb2,ifb3,ifb4,
2683 * ijcvd1(i),ijcvd1(i),4,nchvd1
2695 CALL
hadjet(nhad,amcvd1(i),poj,pat,gacvd1(i),
2697 + (i),bgzvd1(i),ifb1,ifb2,ifb3,ifb4,
2698 * ijcvd1(i),ijcvd1(i),4,nchvd1
2709 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRVD: NHKKNMXHKK ',nhkk,
nmxhkk
2714 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
2718 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
2719 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
2728 ptvd=ptvd+
sqrt(pxf(j)**2+pyf(j)**2)
2731 IF(ibarf(j).EQ.500)istist=2
2733 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),15)
2734 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
2736 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
2737 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
2738 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
2757 IF(amcvd2(i).LT.2.3)
THEN
2758 WRITE(6,
'(A,F10.2,I5)')
' HADRVD AMCVD2(I), I ',
2762 IF(ifb3.LE.8.AND.ifb4.LE.8)
THEN
2764 ELSEIF((ifb3.EQ.9.AND.ifb4.LE.8).OR.
2765 * (ifb4.EQ.9.AND.ifb3.LE.8))
THEN
2767 ELSEIF(ifb3.EQ.9.AND.ifb4.EQ.9)
THEN
2770 CALL
hadjet(nhad,amcvd2(i),poj,pat,gacvd2(i),
2772 + (i),bgzvd2(i),ifb1,ifb2,ifb3,ifb4, ijcvd2(i),
2773 * ijcvd2(i),5,nchvd2
2780 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRVD: NHKKNMXHKK ',
2787 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
2791 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
2792 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
2801 ptvd=ptvd+
sqrt(pxf(j)**2+pyf(j)**2)
2804 IF(ibarf(j).EQ.500)istist=2
2806 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),16)
2807 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
2809 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
2810 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
2811 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
2822 1010
FORMAT (i6,i4,5i6,9e10.2)
2823 1020
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
2824 1030
FORMAT (
' NHKK,IDHKK(NHKK) ',3i10)
2830 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2834 COMMON /zsea/zseaav,zseasu,anzsea
2835 common/popcck/pdbck,pdbse,pdbseu,
2836 * ijpock,irejck,ick4,ihad4,ick6,ihad6
2837 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
2838 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
2839 *isea43,isea63,irejao
2843 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
2849 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
2851 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
2854 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
2858 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
2860 * ,xpsu(248),xtsu(248)
2861 * ,xpsut(248),xtsut(248)
2863 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
2864 +ixpv,ixps,ixtv,ixts, intvv1(248),
2865 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
2867 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
2881 COMMON /abrds/ amcds1(248),amcds2(248),gacds1(248),gacds2(248),
2882 +bgxds1(248),bgyds1(248),bgzds1(248), bgxds2(248),bgyds2(248),
2883 +bgzds2(248), nchds1(248),nchds2(248),ijcds1(248),ijcds2(248),
2884 +pqdsa1(248,4),pqdsa2(248,4), pqdsb1(248,4),pqdsb2(248,4)
2886 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
2888 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
2890 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
2891 *idzre(3),izdre(3),idiqrz(7)
2892 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
2893 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
2894 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
2895 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
2896 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
2897 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
2900 parameter(ummm=0.3d0)
2901 parameter(smmm=0.5d0)
2902 parameter(cmmm=1.3d0)
2920 betcha=betoo+1.3d0-log10(ecm)
2939 WRITE(6,4567)pc,betcha,pu1,ps1,seasq
2940 4567
FORMAT(
' Charm chain ends DIQDSS: PC,BETCHA,PU,PS,SEASQ',5f10.5)
2943 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' diqdss'
2947 ipsq2(ips)=1.d0+
rndm(v1)*(2.d0+2.d0*seasq)
2949 IF(rr.LT.pc)ipsq2(ips)=4
2951 ipsaq2(ips)=-ipsq2(ips)
2964 xall=xdfree+xpsq(ips)+xpsaq(ips)-2.*xdthr
2976 xpvd(iitop)=xdthr+dx1
2978 xpsaq(ips)=xdthr+dx3
2981 amdsq1=xpsq(ips)*xtsq(its)*ecm**2
2982 amdsq2=xpsaq(ips)*xtsaq(its)*ecm**2
2983 idiqre(1)=idiqre(1)+1
2984 IF(ipsq(ips).GE.3.AND.ipsq2(ips).GE.3)
THEN
2985 idiqre(2)=idiqre(2)+1
2987 IF(amdsq2.LE.6.6d0.OR.amdsq1.LE.6.60d0)
THEN
2989 idiqre(3)=idiqre(3)+1
2990 idiqre(2)=idiqre(2)-1
2991 idiqre(1)=idiqre(1)-1
2997 ELSEIF(ipsq(ips).GE.3.OR.ipsq2(ips).GE.3)
THEN
2998 idiqre(4)=idiqre(4)+1
3000 IF(amdsq2.LE.5.8d0.OR.amdsq1.LE.5.80d0)
THEN
3002 idiqre(5)=idiqre(5)+1
3003 idiqre(4)=idiqre(4)-1
3004 idiqre(1)=idiqre(1)-1
3011 idiqre(6)=idiqre(6)+1
3013 IF(amdsq2.LE.3.9d0.OR.amdsq1.LE.3.9d0)
THEN
3015 idiqre(7)=idiqre(7)+1
3016 idiqre(6)=idiqre(6)-1
3017 idiqre(1)=idiqre(1)-1
3041 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3046 COMMON /zsea/zseaav,zseasu,anzsea
3047 common/popcck/pdbck,pdbse,pdbseu,
3048 * ijpock,irejck,ick4,ihad4,ick6,ihad6
3049 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
3050 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
3051 *isea43,isea63,irejao
3058 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
3146 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
3152 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
3154 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
3157 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
3159 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
3160 +ixpv,ixps,ixtv,ixts, intvv1(248),
3161 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
3163 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
3177 COMMON /abrds/ amcds1(248),amcds2(248),gacds1(248),gacds2(248),
3178 +bgxds1(248),bgyds1(248),bgzds1(248), bgxds2(248),bgyds2(248),
3179 +bgzds2(248), nchds1(248),nchds2(248),ijcds1(248),ijcds2(248),
3180 +pqdsa1(248,4),pqdsa2(248,4), pqdsb1(248,4),pqdsb2(248,4)
3181 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
3182 *idzre(3),izdre(3),idiqrz(7)
3186 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
3188 * ,xpsu(248),xtsu(248)
3189 * ,xpsut(248),xtsut(248)
3191 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
3192 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
3199 COMMON /trafop/ gamp,bgamp,betp
3201 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
3202 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
3203 +prebin,taebin,fermod,etacou
3205 COMMON /fermi/ pquar(4,248),paquar(4,248), tquar(4,248),taquar
3218 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
3219 +iibar(210),k1(210),k2(210)
3222 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
3224 COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
3225 +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
3226 +irvs14, irvv11,irvv12,irvv13,irvv14
3228 COMMON /projk/ iprojk
3229 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
3230 common/rptshm/rproj,rtarg,bimpac
3231 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
3232 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
3233 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
3234 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
3235 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
3236 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
3239 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' kkevds'
3243 IF(nchds1(
n).EQ.99.AND.nchds2(
n).EQ.99)go to 10
3246 IF(iphkk.GE.7)
WRITE(6,
'(A,2I10)')
' KKEVDS N,NDS',
n,nds
3248 IF(iphkk.GE.7)
WRITE(6,
'(A,2I10)')
' KKEVDS N,IXSPR',
n,ixspr
3249 inucpr=ifrosp(ixspr)
3250 jnucpr=itovp(inucpr)
3251 IF(iphkk.GE.7)
WRITE(6,
'(A,2I10)')
' KKEVDS INUCPR,JNUCPR',
3254 psqpx=xpsq(ixspr)*prmom(1,inucpr)
3255 psqpy=xpsq(ixspr)*prmom(2,inucpr)
3256 psqpz=xpsq(ixspr)*prmom(3,inucpr)
3257 psqe=xpsq(ixspr)*prmom(4,inucpr)
3258 psaqpx=xpsaq(ixspr)*prmom(1,inucpr)
3259 psaqpy=xpsaq(ixspr)*prmom(2,inucpr)
3260 psaqpz=xpsaq(ixspr)*prmom(3,inucpr)
3261 psaqe=xpsaq(ixspr)*prmom(4,inucpr)
3265 IF(iphkk.GE.7)
WRITE(6,
'(A,2I10)')
' KKEVDS N,IXSTA',
n,ixsta
3266 inucta=ifrost(ixsta)
3267 jnucta=itovt(inucta)
3268 IF(iphkk.GE.7)
WRITE(6,
'(A,2I10)')
' KKEVDS INUCTA,JNUCTA',
3271 tsqpx=xtsq(ixsta)*tamom(1,inucta)
3272 tsqpy=xtsq(ixsta)*tamom(2,inucta)
3273 tsqpz=xtsq(ixsta)*tamom(3,inucta)
3274 tsqe=xtsq(ixsta)*tamom(4,inucta)
3275 tsdqpx=xtsaq(ixsta)*tamom(1,inucta)
3276 tsdqpy=xtsaq(ixsta)*tamom(2,inucta)
3277 tsdqpz=xtsaq(ixsta)*tamom(3,inucta)
3278 tsdqe=xtsaq(ixsta)*tamom(4,inucta)
3287 rtix=(vhkk(1,itnu))*1.e12-bimpac*0.1
3288 rtiy=vhkk(2,itnu)*1.e12
3289 rtiz=vhkk(3,itnu)*1.e12
3290 CALL
cromsc(psqpx,psqpy,psqpz,psqe,rtix,rtiy,rtiz,
3291 * psqnx,psqny,psqnz,psqne,59)
3296 CALL
cromsc(psaqpx,psaqpy,psaqpz,psaqe,rtix,rtiy,rtiz,
3297 * psaqnx,psaqny,psaqnz,psaqne,60)
3308 rtix=(vhkk(1,itnu))*1.e12-bimpac*0.1
3309 rtiy=vhkk(2,itnu)*1.e12
3310 rtiz=vhkk(3,itnu)*1.e12
3311 CALL
cromsc(tsqpx,tsqpy,tsqpz,tsqe,rtix,rtiy,rtiz,
3312 * tsqnx,tsqny,tsqnz,tsqne,61)
3317 CALL
cromsc(tsdqpx,tsdqpy,tsdqpz,tsdqe,rtix,rtiy,rtiz,
3318 * tsdqnx,tsdqny,tsdqnz,tsdqne,62)
3326 IF(ip.GE.0)go to 1779
3327 psqpz2=psqe**2-psqpx**2-psqpy**2
3328 IF(psqpz2.GE.0.)
THEN
3336 paqpz2=psaqe**2-psaqpx**2-psaqpy**2
3337 IF(paqpz2.GE.0.)
THEN
3345 tsqpz2=tsqe**2-tsqpx**2-tsqpy**2
3346 IF(tsqpz2.GE.0.)
THEN
3354 tdqpz2=tsdqe**2-tsdqpx**2-tsdqpy**2
3355 IF(tdqpz2.GE.0.)
THEN
3356 tsdqpz=-
sqrt(tdqpz2)
3393 WRITE(6,
'(A,I5)')
' KKEVDS - IRDS13=',irds13
3394 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
3395 +
' DS: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
3396 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
3397 + amch1,amch2,irej,ikvala,pttq1,ptta1
3401 CALL
selpt( ptxsq1,ptysq1,plq1,
3402 + eq1,ptxsa1,ptysa1,plaq1,eaq1,
3403 + ptxsa2,ptysa2,plaq2,eaq2,
3404 + ptxsq2,ptysq2,plq2,eq2,
3405 + amch1,amch2,irej,ikvala,pttq1,ptta1,
3409 WRITE(6,
'(A,I5)')
' KKEVDS - IRDS13=',irds13
3410 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
3411 +
' DS: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
3412 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
3413 + amch1,amch2,irej,ikvala,pttq1,ptta1
3416 IF (ipev.GE.7)
WRITE(6,
'(A/5X,I10)')
3421 WRITE(6,
'(A,I5)')
' KKEVDS - IRDS13=',irds13
3422 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
3423 +
' DS: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
3424 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
3425 + amch1,amch2,irej,ikvala,pttq1,ptta1
3432 ptxch1=ptxsq1 + ptxsq2
3433 ptych1=ptysq1 + ptysq2
3436 ptxch2=ptxsa2 + ptxsa1
3437 ptych2=ptysa2 + ptysa1
3438 ptzch2=plaq2 + plaq1
3440 ammm=
sqrt((ech1+ech2)**2-(ptxch1+ptxch2)**2
3441 + -(ptych1+ptych2)**2-(ptzch1+ptzch2)**2)
3444 IF (ipev.GE.6)
WRITE(6,
'(A,I10/A,5F12.5/A,5F12.5)')
3445 +
' DS: IREJ ',irej,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
3446 + amch1,ptxch1,ptych1,ptzch1,ech1,
3447 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
3454 CALL
cobcma(ipsq(ixspr),ipsq2(ixspr),itsq(ixsta), ijnch1,nnch1,
3455 + irej,amch1,amch1n,1)
3460 WRITE(6,
'(A,I5)')
' KKEVDS - IRDS11=',irds11
3461 WRITE(6,
'(A,6I5/6E12.4/2E12.4)')
' DS:', ipsq(ixspr),ittv1
3462 + (ixsta),ittv2(ixsta),ijnch1,nnch1,irej, xpsq(ixspr),xpsaq
3463 + (ixspr),xpsqcm,xpsacm, xtvq(ixsta),xtvd(ixsta),amch1,amch1n
3470 CALL
cormom(amch1,amch2,amch1n,amch2n,
3471 + ptxsq1,ptysq1,plq1,eq1,
3472 + ptxsa1,ptysa1,plaq1,eaq1,
3473 + ptxsa2,ptysa2,plaq2,eaq2,
3474 + ptxsq2,ptysq2,plq2,eq2,
3475 + ptxch1,ptych1,ptzch1,ech1, ptxch2,ptych2,ptzch2,ech2,
3479 IF(irej.EQ.1)go to 11
3481 IF (ipev.GE.6)
WRITE(6,
'(A,I10/A,5F12.5/A,5F12.5)')
3482 +
' DS(2): IREJ ',irej,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
3483 + amch1,ptxch1,ptych1,ptzch1,ech1,
3484 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
3490 CALL
cobcma(ipsaq(ixspr),ipsaq2(ixspr),itsaq(ixsta),
3491 + ijnch2,nnch2,irej,amch2,amch2n,2)
3496 WRITE(6,1090) irds12
3497 WRITE(6,1100) ipsaq(ixspr),ipsaq2(ixspr),itsaq(ixsta),
3498 + ijnch2,nnch2,irej,
3499 + xpsq(ixspr),xpsaq(ixspr),xpsqcm,xpsacm, xtsq(ixsta),xtsaq
3500 + (ixsta),xtsqcm,xtsacm, amch2,amch2n
3501 1090
FORMAT(
' KKEVDS - IRDS12=',i5)
3502 1100
FORMAT(
' DS - 1100', 6i5/2(4e12.4/),2e12.4)
3512 ammm=
sqrt((ech1+ech2)**2-(ptxch1+ptxch2)**2
3513 + -(ptych1+ptych2)**2-(ptzch1+ptzch2)**2)
3518 gammm=eee/(ammm+1.
e-4)
3519 bgggx=pxxx/(ammm+1.
e-4)
3520 bgggy=pyyy/(ammm+1.
e-4)
3521 bgggz=pzzz/(ammm+1.
e-4)
3525 CALL
daltra(gammm,-bgggx,-bgggy,-bgggz,
3526 + ptxch1,ptych1,ptzch1,ech1,
3527 + pppch1, qtxch1,qtych1,qtzch1,qech1)
3529 CALL
daltra(gammm,-bgggx,-bgggy,-bgggz,
3530 + ptxch2,ptych2,ptzch2,ech2,
3531 + pppch2, qtxch2,qtych2,qtzch2,qech2)
3534 CALL
corval(ammm,irej,amch1,amch2, qtxch1,qtych1,qtzch1,qech1,
3535 + qtxch2,qtych2,qtzch2,qech2,norig)
3539 CALL
daltra(gammm,bgggx,bgggy,bgggz, qtxch1,qtych1,qtzch1,qech1,
3540 + pppch1, ptxch1,ptych1,ptzch1,ech1)
3542 CALL
daltra(gammm,bgggx,bgggy,bgggz, qtxch2,qtych2,qtzch2,qech2,
3543 + pppch2, ptxch2,ptych2,ptzch2,ech2)
3548 WRITE(6,
'(A/3(1PE15.4),3I5)')
3549 +
' DS - CALL CORVAL: AMMM, AMCH1, AMCH2, NNCH1, NNCH2, IREJ',
3550 + ammm, amch1, amch2, nnch1, nnch2, irej
3551 WRITE(6,1050) irej, amch1,
3552 + ptxch1,ptych1,ptzch1,ech1, amch2,ptxch2,ptych2,ptzch2,ech2
3553 1050
FORMAT (
' DS: IREJ || ',i10/
3554 +
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',5f12.5/
3555 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ',5f12.5)
3601 ihkkpd=jhkkps(ixspr )
3602 ihkkpo=jhkkps(ixspr )-1
3603 ihkktd=jhkkts(ixsta )
3604 ihkkto=jhkkts(ixsta )-1
3605 IF (ipev.GT.3)
WRITE(6,1000)ixspr,inucpr,jnucpr,ihkkpo,ihkkpd
3606 1000
FORMAT (
' IXSPR,INUCPR,JNUCPR,IHKKPO,IHKKPD ',5i5)
3607 IF (ipev.GT.3)
WRITE(6,1010)ixsta,inucta,jnucta,ihkkto,ihkktd
3608 1010
FORMAT (
' IXSTA,INUCTA,JNUCTA,IHKKTO,IHKKTD ',5i5)
3612 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
3617 idhkk(ihkk)=idhkk(ihkkpo)
3618 jmohkk(1,ihkk)=ihkkpo
3619 jmohkk(2,ihkk)=jmohkk(1,ihkkpo)
3620 jdahkk(1,ihkk)=ihkk+2
3621 jdahkk(2,ihkk)=ihkk+2
3622 phkk(1,ihkk)=pqdsa1(
n,1)
3623 phkk(2,ihkk)=pqdsa1(
n,2)
3624 phkk(3,ihkk)=pqdsa1(
n,3)
3625 phkk(4,ihkk)=pqdsa1(
n,4)
3629 vhkk(1,ihkk)=vhkk(1,ihkkpo)+xxpp
3630 vhkk(2,ihkk)=vhkk(2,ihkkpo)+yypp
3631 vhkk(3,ihkk)=vhkk(3,ihkkpo)
3632 vhkk(4,ihkk)=vhkk(4,ihkkpo)
3633 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
3634 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
3635 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
3637 1020
FORMAT (i6,i4,5i6,9e10.2)
3641 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
3646 idhkk(ihkk)=idhkk(ihkktd)
3647 jmohkk(1,ihkk)=ihkktd
3648 jmohkk(2,ihkk)=jmohkk(1,ihkktd)
3649 jdahkk(1,ihkk)=ihkk+1
3650 jdahkk(2,ihkk)=ihkk+1
3651 phkk(1,ihkk)=pqdsa2(
n,1)
3652 phkk(2,ihkk)=pqdsa2(
n,2)
3653 phkk(3,ihkk)=pqdsa2(
n,3)
3654 phkk(4,ihkk)=pqdsa2(
n,4)
3658 vhkk(1,ihkk)=vhkk(1,ihkktd)+xxpp
3659 vhkk(2,ihkk)=vhkk(2,ihkktd)+yypp
3660 vhkk(3,ihkk)=vhkk(3,ihkktd)
3661 vhkk(4,ihkk)=vhkk(4,ihkktd)
3662 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
3663 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
3664 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
3670 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
3675 idhkk(ihkk)=88888+nnch1
3676 jmohkk(1,ihkk)=ihkk-2
3677 jmohkk(2,ihkk)=ihkk-1
3688 vhkk(1,nhkk)= vhkk(1,nhkk-1)
3689 vhkk(2,nhkk)= vhkk(2,nhkk-1)
3690 vhkk(3,nhkk)= vhkk(3,nhkk-1)
3691 vhkk(4,nhkk)=vhkk(3,nhkk)/betp-vhkk(3,nhkk-2)/bgamp
3693 IF (iprojk.EQ.1)
THEN
3694 whkk(1,nhkk)= vhkk(1,nhkk-2)
3695 whkk(2,nhkk)= vhkk(2,nhkk-2)
3696 whkk(3,nhkk)= vhkk(3,nhkk-2)
3697 whkk(4,nhkk)=vhkk(4,nhkk)*gamp-vhkk(3,nhkk)*bgamp
3698 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
3699 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
3700 + (phkk(khkk,ihkk),khkk=1,5), (whkk(khkk,ihkk),khkk=1,4)
3703 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
3704 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
3705 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
3712 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
3717 idhkk(ihkk)=idhkk(ihkkpd)
3718 jmohkk(1,ihkk)=ihkkpd
3719 jmohkk(2,ihkk)=jmohkk(1,ihkkpd)
3720 jdahkk(1,ihkk)=ihkk+2
3721 jdahkk(2,ihkk)=ihkk+2
3722 phkk(1,ihkk)=pqdsb1(
n,1)
3723 phkk(2,ihkk)=pqdsb1(
n,2)
3724 phkk(3,ihkk)=pqdsb1(
n,3)
3725 phkk(4,ihkk)=pqdsb1(
n,4)
3729 vhkk(1,ihkk)=vhkk(1,ihkkpd)+xxpp
3730 vhkk(2,ihkk)=vhkk(2,ihkkpd)+yypp
3731 vhkk(3,ihkk)=vhkk(3,ihkkpd)
3732 vhkk(4,ihkk)=vhkk(4,ihkkpd)
3733 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
3734 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
3735 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
3740 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
3745 idhkk(ihkk)=idhkk(ihkkto)
3746 jmohkk(1,ihkk)=ihkkto
3747 jmohkk(2,ihkk)=jmohkk(1,ihkkto)
3748 jdahkk(1,ihkk)=ihkk+1
3749 jdahkk(2,ihkk)=ihkk+1
3750 phkk(1,ihkk)=pqdsb2(
n,1)
3751 phkk(2,ihkk)=pqdsb2(
n,2)
3752 phkk(3,ihkk)=pqdsb2(
n,3)
3753 phkk(4,ihkk)=pqdsb2(
n,4)
3757 vhkk(1,ihkk)=vhkk(1,ihkkto)+xxpp
3758 vhkk(2,ihkk)=vhkk(2,ihkkto)+yypp
3759 vhkk(3,ihkk)=vhkk(3,ihkkto)
3760 vhkk(4,ihkk)=vhkk(4,ihkkto)
3761 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
3762 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
3763 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
3769 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
3774 idhkk(ihkk)=88888+nnch2
3775 jmohkk(1,ihkk)=ihkk-2
3776 jmohkk(2,ihkk)=ihkk-1
3787 vhkk(1,nhkk)= vhkk(1,nhkk-1)
3788 vhkk(2,nhkk)= vhkk(2,nhkk-1)
3789 vhkk(3,nhkk)= vhkk(3,nhkk-1)
3790 vhkk(4,nhkk)=vhkk(3,nhkk)/betp-vhkk(3,nhkk-2)/bgamp
3792 IF (iprojk.EQ.1)
THEN
3793 whkk(1,nhkk)= vhkk(1,nhkk-2)
3794 whkk(2,nhkk)= vhkk(2,nhkk-2)
3795 whkk(3,nhkk)= vhkk(3,nhkk-2)
3796 whkk(4,nhkk)=vhkk(4,nhkk)*gamp-vhkk(3,nhkk)*bgamp
3797 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
3798 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
3799 + (phkk(khkk,ihkk),khkk=1,5), (whkk(khkk,ihkk),khkk=1,4)
3802 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
3803 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
3804 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
3812 gacds1(
n)=qech1/amch1
3813 bgxds1(
n)=qtxch1/amch1
3814 bgyds1(
n)=qtych1/amch1
3815 bgzds1(
n)=qtzch1/amch1
3816 gacds2(
n)=qech2/amch2
3817 bgxds2(
n)=qtxch2/amch2
3818 bgyds2(
n)=qtych2/amch2
3819 bgzds2(
n)=qtzch2/amch2
3824 IF (ipev.GE.6)
WRITE(6,
'(A/I10,4F12.7,5I5/10X,4F12.6/10X,6F12.6,
3825 +4I5/8F15.5/ 8F15.5)')
' DS / FINAL PRINT',
n
3838 xpvd(jnucpr)=xpvd(jnucpr) + xpsq(ixspr) + xpsaq(ixspr)
3839 xtvd(jnucta)=xtvd(jnucta) + xtsaq(ixsta) + xtsq(ixsta)
3840 issqq=abs(ipsaq(ixspr))
3841 jssqq=abs(ipsaq2(ixspr))
3842 IF(issqq.EQ.3.AND.jssqq.EQ.3)
THEN
3844 ELSEIF(issqq.EQ.3.OR.jssqq.EQ.3)
THEN
3859 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3870 COMMON /zsea/zseaav,zseasu,anzsea
3871 common/popcck/pdbck,pdbse,pdbseu,
3872 * ijpock,irejck,ick4,ihad4,ick6,ihad6
3873 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
3874 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
3875 *isea43,isea63,irejao
3879 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
3885 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
3887 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
3890 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
3894 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
3896 * ,xpsu(248),xtsu(248)
3897 * ,xpsut(248),xtsut(248)
3899 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
3900 +ixpv,ixps,ixtv,ixts, intvv1(248),
3901 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
3903 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
3917 COMMON /abrds/ amcds1(248),amcds2(248),gacds1(248),gacds2(248),
3918 +bgxds1(248),bgyds1(248),bgzds1(248), bgxds2(248),bgyds2(248),
3919 +bgzds2(248), nchds1(248),nchds2(248),ijcds1(248),ijcds2(248),
3920 +pqdsa1(248,4),pqdsa2(248,4), pqdsb1(248,4),pqdsb2(248,4)
3922 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
3923 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
3933 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
4022 parameter(nfimax=249)
4023 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
4024 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
4025 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
4028 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4030 COMMON /projk/ iprojk
4032 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
4034 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
4035 * anndv,annvd,annds,annsd,
4037 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
4039 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
4042 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
4043 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
4044 * acouzz,acouhh,acouds,acousd,
4045 * acoudz,acouzd,acoudi,
4046 * acoudv,acouvd,acoucc
4047 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
4048 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
4049 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
4050 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
4051 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
4052 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
4054 dimension poj(4),pat(4)
4057 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' hadrds'
4061 IF(nchds1(i).EQ.99.AND.nchds2(i).EQ.99) go to 50
4065 IF (ipco.GE.6)
WRITE (6,1000) ipsq(is1),ipsaq(is1),itsq(is2),
4066 + itsaq(is2),ittv2(is2), amcds1(i),amcds2(i),gacds1(i),gacds2(i),
4067 + bgxds1(i),bgyds1(i),bgzds1(i), bgxds2(i),bgyds2(i),bgzds2(i),
4068 + nchds1(i),nchds2(i),ijcds1(i),ijcds2(i), pqdsa1(i,4),pqdsa2
4069 + (i,4),pqdsb1(i,4),pqdsb2(i,4)
4070 1000
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
4080 IF((nchds1(i).NE.0.OR.nchds2(i).NE.0).AND.ip.NE.1)
4081 & CALL
saptre(amcds1(i),gacds1(i),bgxds1(i),bgyds1(i),bgzds1(i),
4082 & amcds2(i),gacds2(i),bgxds2(i),bgyds2(i),bgzds2(i))
4085 IF(ipco.GE.3)
WRITE (6,1244) poj,pat
4086 1244
FORMAT (
' D-S QUARK-DIQUARK POJ,PAT ',8e12.3)
4116 IF(ipco.GE.3)
WRITE(6,
'(A,3I5)')
'HADRDS: I,IPPP,JIPP ',
4122 IF(ifb1.LE.2.AND.ifb2.LE.2)
THEN
4124 ELSEIF((ifb1.EQ.3.AND.ifb2.LE.2).OR.
4125 * (ifb2.EQ.3.AND.ifb1.LE.2))
THEN
4127 ELSEIF(ifb1.EQ.3.AND.ifb2.EQ.3)
THEN
4130 IF((nchds1(i).NE.0))
4131 * CALL
hadjet(nhad,amcds1(i),poj,pat,gacds1(i),
4133 + (i),bgzds1(i),ifb1,ifb2,ifb3,ifb4, ijcds1(i),
4134 * ijcds1(i),6,nchds1
4137 aack=float(ick6)/float(ick6+ihad6+1)
4138 IF((nchds1(i).EQ.0))
THEN
4139 zseawu=
rndm(bb)*2.d0*zseaav
4140 rseack=float(jitt)*pdbse +zseawu*pdbseu
4141 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,2F10.3)')
'HADJSE JIPP,',
4142 *
'RSEACK,PDBSE 3 dpmnuc5',
4145 IF(
rndm(v).LE.rseack)
THEN
4147 IF(amcds1(i).GT.2.3d0)
THEN
4149 CALL
hadjse(nhad,amcds1(i),poj,pat,gacds1(i),bgxds1(i),
4151 + (i),bgzds1(i),ifb1,ifb2,ifb3,ifb4, ijcds1(i),ijcds1(i),6,
4153 + (i),6,irejss,iissqq)
4154 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JIPP,',
4155 *
'RSEACK,IREJSS 3 dpmnuc5 ',
4156 + jipp,rseack,irejss
4159 IF(irejss.EQ.1)irejse=irejse+1
4160 IF(irejss.EQ.3)irejs3=irejs3+1
4161 IF(irejss.EQ.2)irejs0=irejs0+1
4162 CALL
hadjet(nhad,amcds1(i),poj,pat,gacds1(i),
4164 + (i),bgzds1(i),ifb1,ifb2,ifb3,ifb4, ijcds1(i),
4165 * ijcds1(i),6,nchds1
4177 CALL
hadjet(nhad,amcds1(i),poj,pat,gacds1(i),
4179 + (i),bgzds1(i),ifb1,ifb2,ifb3,ifb4, ijcds1(i),
4180 * ijcds1(i),6,nchds1
4191 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRDS: NHKKNMXHKK ',nhkk,
nmxhkk
4196 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
4200 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
4201 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
4210 ptds=ptds+
sqrt(pxf(j)**2+pyf(j)**2)
4213 IF(ibarf(j).EQ.500)istist=2
4214 IF(ipco.GE.3)
WRITE(6,*)
' HADRDS before HKKFIL'
4216 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),17)
4217 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
4219 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
4220 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
4221 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
4273 IF(ifb1.LE.8.AND.ifb2.LE.8)
THEN
4275 ELSEIF((ifb1.EQ.9.AND.ifb2.LE.8).OR.
4276 * (ifb2.EQ.9.AND.ifb1.LE.8))
THEN
4278 ELSEIF(ifb1.EQ.9.AND.ifb2.EQ.9)
THEN
4284 IF((nchds2(i).NE.0))
4285 * CALL
hadjet(nhad,amcds2(i),poj,pat,gacds2(i),
4287 + (i),bgzds2(i),ifb1,ifb2,ifb3,ifb4, ijcds2(i),
4288 * ijcds2(i),6,nchds2
4292 aack=float(ick6)/float(ick6+ihad6+1)
4293 IF((nchds2(i).EQ.0))
THEN
4294 zseawu=
rndm(bb)*2.d0*zseaav
4295 rseack=float(jitt)*pdbse +zseawu*pdbseu
4296 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,2F10.3)')
'HADJSE JIPP,',
4300 IF(
rndm(v).LE.rseack)
THEN
4302 IF(amcds2(i).GT.2.3d0)
THEN
4305 CALL
hadjase(nhad,amcds2(i),poj,pat,gacds2(i),bgxds2(i),
4307 + (i),bgzds2(i),ifb1,ifb2,ifb3,ifb4, ijcds2(i),ijcds2(i),6,
4309 + (i),6,irejss,iissqq)
4310 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JIPP,',
4312 + jipp,rseack,irejss
4315 IF(irejss.EQ.1)irejsa=irejsa+1
4316 IF(irejss.EQ.3)ireja3=ireja3+1
4317 IF(irejss.EQ.2)ireja0=ireja0+1
4319 CALL
hadjet(nhad,amcds2(i),poj,pat,gacds2(i),
4321 + (i),bgzds2(i),ifb1,ifb2,ifb3,ifb4, ijcds2(i),
4322 * ijcds2(i),6,nchds2
4335 CALL
hadjet(nhad,amcds2(i),poj,pat,gacds2(i),
4337 + (i),bgzds2(i),ifb1,ifb2,ifb3,ifb4, ijcds2(i),
4338 * ijcds2(i),6,nchds2
4349 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRDS: NHKKNMXHKK ', nhkk,
4355 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
4359 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
4360 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
4369 ptds=ptds+
sqrt(pxf(j)**2+pyf(j)**2)
4372 IF(ibarf(j).EQ.500)istist=2
4373 IF(ipco.GE.3)
WRITE(6,*)
' HADRDS before 2 HKKFIL'
4375 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),18)
4376 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
4378 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
4379 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
4380 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
4391 1010
FORMAT (i6,i4,5i6,9e10.2)
4392 1020
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
4393 1030
FORMAT (
' NHKK,IDHKK(NHKK) ',3i10)
4399 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4403 COMMON /zsea/zseaav,zseasu,anzsea
4404 common/popcck/pdbck,pdbse,pdbseu,
4405 * ijpock,irejck,ick4,ihad4,ick6,ihad6
4406 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
4407 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
4408 *isea43,isea63,irejao
4412 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
4418 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
4420 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
4423 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
4427 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
4429 * ,xpsu(248),xtsu(248)
4430 * ,xpsut(248),xtsut(248)
4432 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
4433 +ixpv,ixps,ixtv,ixts, intvv1(248),
4434 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
4436 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
4450 COMMON /abrsd/ amcsd1(248),amcsd2(248),gacsd1(248),gacsd2(248),
4451 +bgxsd1(248),bgysd1(248),bgzsd1(248), bgxsd2(248),bgysd2(248),
4452 +bgzsd2(248), nchsd1(248),nchsd2(248),ijcsd1(248),ijcsd2(248),
4453 +pqsda1(248,4),pqsda2(248,4), pqsdb1(248,4),pqsdb2(248,4)
4455 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4457 COMMON /xseadi/ xseacu,unon,unom,unosea, cvq,cdq,csea,ssmima,
4459 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
4460 *idzre(3),izdre(3),idiqrz(7)
4461 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
4462 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
4463 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
4464 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
4465 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
4466 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
4469 parameter(ummm=0.3d0)
4470 parameter(smmm=0.5d0)
4471 parameter(cmmm=1.3d0)
4489 betcha=betoo+1.3d0-log10(ecm)
4508 WRITE(6,4567)pc,betcha,pu1,ps1,seasq
4509 4567
FORMAT(
' Charm chain ends DIQSSD: PC,BETCHA,PU,PS,SEASQ',4f10.5)
4513 is=1.d0+
rndm(v1)*(2.d0+2.d0*seasq)
4516 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' diqssd'
4520 itsq2(its)=1.d0+
rndm(v1)*(2.d0+2.d0*seasq)
4522 IF(rr.LT.pc)itsq2(its)=4
4524 itsaq2(its)=-itsq2(its)
4537 xall=xdfree+xtsq(its)+xtsaq(its)-2.*xdthr
4549 xtvd(iitot)=xdthr+dx1
4551 xtsaq(its)=xdthr+dx3
4554 amsdq1=xtsq(its)*xpsq(ips)*ecm**2
4555 amsdq2=xtsaq(its)*xpsaq(ips)*ecm**2
4556 idiqre(1)=idiqre(1)+1
4557 IF(itsq(its).GE.3.AND.itsq2(its).GE.3)
THEN
4558 idiqre(2)=idiqre(2)+1
4560 IF(amsdq2.LE.6.60d0.OR.amsdq1.LE.6.60d0)
THEN
4562 idiqre(3)=idiqre(3)+1
4563 idiqre(2)=idiqre(2)-1
4564 idiqre(1)=idiqre(1)-1
4570 ELSEIF(itsq(its).GE.3.OR.itsq2(its).GE.3)
THEN
4571 idiqre(4)=idiqre(4)+1
4573 IF(amsdq2.LE.5.8d0.OR.amsdq1.LE.5.80d0)
THEN
4575 idiqre(5)=idiqre(5)+1
4576 idiqre(4)=idiqre(4)-1
4577 idiqre(1)=idiqre(1)-1
4584 idiqre(6)=idiqre(6)+1
4586 IF(amsdq2.LE.3.9d0.OR.amsdq1.LE.3.9d0)
THEN
4588 idiqre(7)=idiqre(7)+1
4589 idiqre(6)=idiqre(6)-1
4590 idiqre(1)=idiqre(1)-1
4613 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4618 COMMON /zsea/zseaav,zseasu,anzsea
4619 common/popcck/pdbck,pdbse,pdbseu,
4620 * ijpock,irejck,ick4,ihad4,ick6,ihad6
4621 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
4622 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
4623 *isea43,isea63,irejao
4630 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
4718 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
4724 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
4726 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
4729 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
4731 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
4732 +ixpv,ixps,ixtv,ixts, intvv1(248),
4733 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
4735 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
4749 COMMON /abrsd/ amcsd1(248),amcsd2(248),gacsd1(248),gacsd2(248),
4750 +bgxsd1(248),bgysd1(248),bgzsd1(248), bgxsd2(248),bgysd2(248),
4751 +bgzsd2(248), nchsd1(248),nchsd2(248),ijcsd1(248),ijcsd2(248),
4752 +pqsda1(248,4),pqsda2(248,4), pqsdb1(248,4),pqsdb2(248,4)
4756 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
4758 * ,xpsu(248),xtsu(248)
4759 * ,xpsut(248),xtsut(248)
4760 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
4761 *idzre(3),izdre(3),idiqrz(7)
4763 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
4764 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
4771 COMMON /trafop/ gamp,bgamp,betp
4773 COMMON /nucimp/ prmom(5,248),tamom(5,248), prmfep,prmfen,tamfep,
4774 +tamfen, prefep,prefen,taefep,taefen, prepot(210),taepot(210),
4775 +prebin,taebin,fermod,etacou
4777 COMMON /fermi/ pquar(4,248),paquar(4,248), tquar(4,248),taquar
4790 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
4791 +iibar(210),k1(210),k2(210)
4794 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4796 COMMON /rejec/ irco1,irco2,irco3,irco4,irco5, irss11,irss12,
4797 +irss13,irss14, irsv11,irsv12,irsv13,irsv14, irvs11,irvs12,irvs13,
4798 +irvs14, irvv11,irvv12,irvv13,irvv14
4800 COMMON /projk/ iprojk
4801 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
4802 common/rptshm/rproj,rtarg,bimpac
4803 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
4804 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
4805 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
4806 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
4807 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
4808 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
4811 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' kkevsd'
4815 IF(nchsd1(
n).EQ.99.AND.nchsd2(
n).EQ.99)go to 10
4819 inucpr=ifrosp(ixspr)
4820 jnucpr=itovp(inucpr)
4822 psqpx=xpsq(ixspr)*prmom(1,inucpr)
4823 psqpy=xpsq(ixspr)*prmom(2,inucpr)
4824 psqpz=xpsq(ixspr)*prmom(3,inucpr)
4825 psqe=xpsq(ixspr)*prmom(4,inucpr)
4826 psdqpx=xpsaq(ixspr)*prmom(1,inucpr)
4827 psdqpy=xpsaq(ixspr)*prmom(2,inucpr)
4828 psdqpz=xpsaq(ixspr)*prmom(3,inucpr)
4829 psdqe=xpsaq(ixspr)*prmom(4,inucpr)
4833 inucta=ifrost(ixsta)
4834 jnucta=itovt(inucta)
4836 tsqpx=xtsq(ixsta)*tamom(1,inucta)
4837 tsqpy=xtsq(ixsta)*tamom(2,inucta)
4838 tsqpz=xtsq(ixsta)*tamom(3,inucta)
4839 tsqe=xtsq(ixsta)*tamom(4,inucta)
4840 tsaqpx=xtsaq(ixsta)*tamom(1,inucta)
4841 tsaqpy=xtsaq(ixsta)*tamom(2,inucta)
4842 tsaqpz=xtsaq(ixsta)*tamom(3,inucta)
4843 tsaqe=xtsaq(ixsta)*tamom(4,inucta)
4852 rtix=(vhkk(1,itnu))*1.e12-bimpac*0.1
4853 rtiy=vhkk(2,itnu)*1.e12
4854 rtiz=vhkk(3,itnu)*1.e12
4855 CALL
cromsc(psqpx,psqpy,psqpz,psqe,rtix,rtiy,rtiz,
4856 * psqnx,psqny,psqnz,psqne,63)
4861 CALL
cromsc(psdqpx,psdqpy,psdqpz,psdqe,rtix,rtiy,rtiz,
4862 * psdqnx,psdqny,psdqnz,psdqne,64)
4874 rtix=(vhkk(1,itnu))*1.e12-bimpac*0.1
4875 rtiy=vhkk(2,itnu)*1.e12
4876 rtiz=vhkk(3,itnu)*1.e12
4877 CALL
cromsc(tsqpx,tsqpy,tsqpz,tsqe,rtix,rtiy,rtiz,
4878 * tsqnx,tsqny,tsqnz,tsqne,65)
4883 CALL
cromsc(tsaqpx,tsaqpy,tsaqpz,tsaqe,rtix,rtiy,rtiz,
4884 * tsaqnx,tsaqny,tsaqnz,tsaqne,66)
4892 IF(ip.GE.0)go to 1779
4893 psqpz2=psqe**2-psqpx**2-psqpy**2
4894 IF(psqpz2.GE.0.)
THEN
4902 pdqpz2=psdqe**2-psdqpx**2-psdqpy**2
4903 IF(pdqpz2.GE.0.)
THEN
4911 tsqpz2=tsqe**2-tsqpx**2-tsqpy**2
4912 IF(tsqpz2.GE.0.)
THEN
4920 taqpz2=tsaqe**2-tsaqpx**2-tsaqpy**2
4921 IF(taqpz2.GE.0.)
THEN
4922 tsaqpz=-
sqrt(taqpz2)
4959 WRITE(6,
'(A,I5)')
' KKEVSD - IRSD13=',irsd13
4960 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
4961 +
' VD: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
4962 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
4963 + amch1,amch2,irej,ikvala,pttq1,ptta1
4967 CALL
selpt( ptxsq1,ptysq1,plq1,
4968 + eq1,ptxsa1,ptysa1,plaq1,eaq1,
4969 + ptxsa2,ptysa2,plaq2,eaq2,
4970 + ptxsq2,ptysq2,plq2,eq2,
4971 + amch1,amch2,irej,ikvala,pttq1,ptta1,
4975 WRITE(6,
'(A,I5)')
' KKEVSD - IRSD13=',irsd13
4976 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
4977 +
' VD: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
4978 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
4979 + amch1,amch2,irej,ikvala,pttq1,ptta1
4981 IF (ipev.GE.7)
WRITE(6,
'(A/5X,I10)')
4986 WRITE(6,
'(A,I5)')
' KKEVSD - IRSD13=',irsd13
4987 WRITE(6,
'(A/4(4E12.4/),2E12.4/2I5/4E12.4)')
4988 +
' VD: ...', ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
4989 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
4990 + amch1,amch2,irej,ikvala,pttq1,ptta1
4998 ptxch1=ptxsq1 + ptxsq2
4999 ptych1=ptysq1 + ptysq2
5002 ptxch2=ptxsa2 + ptxsa1
5003 ptych2=ptysa2 + ptysa1
5004 ptzch2=plaq2 + plaq1
5006 ammm=
sqrt((ech1+ech2)**2-(ptxch1+ptxch2)**2
5007 + -(ptych1+ptych2)**2-(ptzch1+ptzch2)**2)
5010 IF (ipev.GE.6)
WRITE(6,
'(A,I10/A,5F12.5/A,5F12.5)')
5011 +
' SD: IREJ ',irej,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
5012 + amch1,ptxch1,ptych1,ptzch1,ech1,
5013 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
5021 CALL
cobcma(itsq(ixsta),itsq2(ixsta),ipsq(ixspr), ijnch1,nnch1,
5022 + irej,amch1,amch1n,1)
5027 WRITE(6,
'(A,I5)')
' KKEVSD - IRSD11=',irsd11
5028 WRITE(6,
'(A,6I5/6E12.4/2E12.4)')
' SD:', ipvq(ixspr),itsq
5029 + (ixsta),itsq2(ixsta),ijnch1,nnch1,irej, xpvq(ixspr),xpvd
5030 + (ixspr),xpsqcm,xpsdcm, xtsq(ixsta),xtsaq(ixsta),amch1,amch1n
5037 CALL
cormom(amch1,amch2,amch1n,amch2n,
5038 + ptxsq1,ptysq1,plq1,eq1,
5039 + ptxsa1,ptysa1,plaq1,eaq1,
5040 + ptxsa2,ptysa2,plaq2,eaq2,
5041 + ptxsq2,ptysq2,plq2,eq2,
5042 + ptxch1,ptych1,ptzch1,ech1, ptxch2,ptych2,ptzch2,ech2,
5046 IF (irej.EQ.1)go to 11
5048 IF (ipev.GE.6)
WRITE(6,
'(A,I10/A,5F12.5/A,5F12.5)')
5049 +
' SD(2): IREJ ',irej,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
5050 + amch1,ptxch1,ptych1,ptzch1,ech1,
5051 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
5057 CALL
cobcma(ipsaq(ixspr),itsaq(ixsta),itsaq2(ixsta),
5058 + ijnch2,nnch2,irej,amch2,amch2n,2)
5063 WRITE(6,1090) irsd12
5064 WRITE(6,1100) ipsaq(ixspr),itsaq(ixsta),itsaq2(ixsta),
5065 + ijnch2,nnch2,irej,
5066 + xpsq(ixspr),xpsaq(ixspr),xpsqcm,xtsacm, xtsq(ixsta),xtsaq
5067 + (ixsta),xtsqcm,xtsacm, amch2,amch2n
5068 1090
FORMAT(
' KKEVSD - IRSD12=',i5)
5069 1100
FORMAT(
' SD - 1100', 6i5/2(4e12.4/),2e12.4)
5080 ammm=
sqrt((ech1+ech2)**2-(ptxch1+ptxch2)**2
5081 + -(ptych1+ptych2)**2-(ptzch1+ptzch2)**2)
5086 gammm=eee/(ammm+1.
e-4)
5087 bgggx=pxxx/(ammm+1.
e-4)
5088 bgggy=pyyy/(ammm+1.
e-4)
5089 bgggz=pzzz/(ammm+1.
e-4)
5092 CALL
daltra(gammm,-bgggx,-bgggy,-bgggz,
5093 + ptxch1,ptych1,ptzch1,ech1,
5094 + pppch1, qtxch1,qtych1,qtzch1,qech1)
5096 CALL
daltra(gammm,-bgggx,-bgggy,-bgggz,
5097 + ptxch2,ptych2,ptzch2,ech2,
5098 + pppch2, qtxch2,qtych2,qtzch2,qech2)
5101 CALL
corval(ammm,irej,amch1,amch2, qtxch1,qtych1,qtzch1,qech1,
5102 + qtxch2,qtych2,qtzch2,qech2,norig)
5107 CALL
daltra(gammm,bgggx,bgggy,bgggz, qtxch1,qtych1,qtzch1,qech1,
5108 + pppch1, ptxch1,ptych1,ptzch1,ech1)
5110 CALL
daltra(gammm,bgggx,bgggy,bgggz, qtxch2,qtych2,qtzch2,qech2,
5111 + pppch2, ptxch2,ptych2,ptzch2,ech2)
5116 WRITE(6,
'(A/3(1PE15.4),3I5)')
5117 +
' SD - CALL CORVAL: AMMM, AMCH1, AMCH2, NNCH1, NNCH2, IREJ',
5118 + ammm, amch1, amch2, nnch1, nnch2, irej
5119 WRITE(6,1050) irej, amch1,
5120 + ptxch1,ptych1,ptzch1,ech1, amch2,ptxch2,ptych2,ptzch2,ech2
5121 1050
FORMAT (
' SD: IREJ || ',i10/
5122 +
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',5f12.5/
5123 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ',5f12.5)
5169 ihkkpd=jhkkps(ixspr )
5170 ihkkpo=jhkkps(ixspr )-1
5171 ihkktd=jhkkts(ixsta )
5172 ihkkto=jhkkts(ixsta )-1
5173 IF (ipev.GT.3)
WRITE(6,1000)ixspr,inucpr,jnucpr,ihkkpo,ihkkpd
5174 1000
FORMAT (
' IXSPR,INUCPR,JNUCPR,IHKKPO,IHKKPD ',5i5)
5175 IF (ipev.GT.3)
WRITE(6,1010)ixsta,inucta,jnucta,ihkkto,ihkktd
5176 1010
FORMAT (
' IXSTA,INUCTA,JNUCTA,IHKKTO,IHKKTD ',5i5)
5180 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
5185 idhkk(ihkk)=idhkk(ihkkpo)
5186 jmohkk(1,ihkk)=ihkkpo
5187 jmohkk(2,ihkk)=jmohkk(1,ihkkpo)
5188 jdahkk(1,ihkk)=ihkk+2
5189 jdahkk(2,ihkk)=ihkk+2
5190 phkk(1,ihkk)=pqsda1(
n,1)
5191 phkk(2,ihkk)=pqsda1(
n,2)
5192 phkk(3,ihkk)=pqsda1(
n,3)
5193 phkk(4,ihkk)=pqsda1(
n,4)
5197 vhkk(1,ihkk)=vhkk(1,ihkkpo)+xxpp
5198 vhkk(2,ihkk)=vhkk(2,ihkkpo)+yypp
5199 vhkk(3,ihkk)=vhkk(3,ihkkpo)
5200 vhkk(4,ihkk)=vhkk(4,ihkkpo)
5201 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
5202 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
5203 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
5205 1020
FORMAT (i6,i4,5i6,9e10.2)
5209 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
5214 idhkk(ihkk)=idhkk(ihkktd)
5215 jmohkk(1,ihkk)=ihkktd
5216 jmohkk(2,ihkk)=jmohkk(1,ihkktd)
5217 jdahkk(1,ihkk)=ihkk+1
5218 jdahkk(2,ihkk)=ihkk+1
5219 phkk(1,ihkk)=pqsda2(
n,1)
5220 phkk(2,ihkk)=pqsda2(
n,2)
5221 phkk(3,ihkk)=pqsda2(
n,3)
5222 phkk(4,ihkk)=pqsda2(
n,4)
5226 vhkk(1,ihkk)=vhkk(1,ihkktd)+xxpp
5227 vhkk(2,ihkk)=vhkk(2,ihkktd)+yypp
5228 vhkk(3,ihkk)=vhkk(3,ihkktd)
5229 vhkk(4,ihkk)=vhkk(4,ihkktd)
5230 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
5231 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
5232 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
5238 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
5243 idhkk(ihkk)=88888+nnch1
5244 jmohkk(1,ihkk)=ihkk-2
5245 jmohkk(2,ihkk)=ihkk-1
5256 vhkk(1,nhkk)= vhkk(1,nhkk-1)
5257 vhkk(2,nhkk)= vhkk(2,nhkk-1)
5258 vhkk(3,nhkk)= vhkk(3,nhkk-1)
5259 vhkk(4,nhkk)=vhkk(3,nhkk)/betp-vhkk(3,nhkk-2)/bgamp
5261 IF (iprojk.EQ.1)
THEN
5262 whkk(1,nhkk)= vhkk(1,nhkk-2)
5263 whkk(2,nhkk)= vhkk(2,nhkk-2)
5264 whkk(3,nhkk)= vhkk(3,nhkk-2)
5265 whkk(4,nhkk)=vhkk(4,nhkk)*gamp-vhkk(3,nhkk)*bgamp
5266 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
5267 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
5268 + (phkk(khkk,ihkk),khkk=1,5), (whkk(khkk,ihkk),khkk=1,4)
5271 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
5272 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
5273 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
5280 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
5285 idhkk(ihkk)=idhkk(ihkkpd)
5286 jmohkk(1,ihkk)=ihkkpd
5287 jmohkk(2,ihkk)=jmohkk(1,ihkkpd)
5288 jdahkk(1,ihkk)=ihkk+2
5289 jdahkk(2,ihkk)=ihkk+2
5290 phkk(1,ihkk)=pqsdb1(
n,1)
5291 phkk(2,ihkk)=pqsdb1(
n,2)
5292 phkk(3,ihkk)=pqsdb1(
n,3)
5293 phkk(4,ihkk)=pqsdb1(
n,4)
5297 vhkk(1,ihkk)=vhkk(1,ihkkpd)+xxpp
5298 vhkk(2,ihkk)=vhkk(2,ihkkpd)+yypp
5299 vhkk(3,ihkk)=vhkk(3,ihkkpd)
5300 vhkk(4,ihkk)=vhkk(4,ihkkpd)
5301 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
5302 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
5303 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
5308 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
5313 idhkk(ihkk)=idhkk(ihkkto)
5314 jmohkk(1,ihkk)=ihkkto
5315 jmohkk(2,ihkk)=jmohkk(1,ihkkto)
5316 jdahkk(1,ihkk)=ihkk+1
5317 jdahkk(2,ihkk)=ihkk+1
5318 phkk(1,ihkk)=pqsdb2(
n,1)
5319 phkk(2,ihkk)=pqsdb2(
n,2)
5320 phkk(3,ihkk)=pqsdb2(
n,3)
5321 phkk(4,ihkk)=pqsdb2(
n,4)
5325 vhkk(1,ihkk)=vhkk(1,ihkkto)+xxpp
5326 vhkk(2,ihkk)=vhkk(2,ihkkto)+yypp
5327 vhkk(3,ihkk)=vhkk(3,ihkkto)
5328 vhkk(4,ihkk)=vhkk(4,ihkkto)
5329 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
5330 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
5331 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
5337 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
5342 idhkk(ihkk)=88888+nnch2
5343 jmohkk(1,ihkk)=ihkk-2
5344 jmohkk(2,ihkk)=ihkk-1
5355 vhkk(1,nhkk)= vhkk(1,nhkk-1)
5356 vhkk(2,nhkk)= vhkk(2,nhkk-1)
5357 vhkk(3,nhkk)= vhkk(3,nhkk-1)
5358 vhkk(4,nhkk)=vhkk(3,nhkk)/betp-vhkk(3,nhkk-2)/bgamp
5360 IF (iprojk.EQ.1)
THEN
5361 whkk(1,nhkk)= vhkk(1,nhkk-2)
5362 whkk(2,nhkk)= vhkk(2,nhkk-2)
5363 whkk(3,nhkk)= vhkk(3,nhkk-2)
5364 whkk(4,nhkk)=vhkk(4,nhkk)*gamp-vhkk(3,nhkk)*bgamp
5365 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
5366 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
5367 + (phkk(khkk,ihkk),khkk=1,5), (whkk(khkk,ihkk),khkk=1,4)
5370 IF (iphkk.GE.2)
WRITE(6,1020) ihkk,isthkk(ihkk),idhkk(ihkk),
5371 + jmohkk(1,ihkk),jmohkk(2,ihkk), jdahkk(1,ihkk),jdahkk(2,ihkk),
5372 + (phkk(khkk,ihkk),khkk=1,5), (vhkk(khkk,ihkk),khkk=1,4)
5382 gacsd1(
n)=qech1/amch1
5383 bgxsd1(
n)=qtxch1/amch1
5384 bgysd1(
n)=qtych1/amch1
5385 bgzsd1(
n)=qtzch1/amch1
5386 gacsd2(
n)=qech2/amch2
5387 bgxsd2(
n)=qtxch2/amch2
5388 bgysd2(
n)=qtych2/amch2
5389 bgzsd2(
n)=qtzch2/amch2
5394 IF (ipev.GE.2)
WRITE(6,
'(A/I10,4F12.7,5I5/10X,4F12.6/10X,6F12.6,
5395 +4I5/8F15.5/8F15.5/2I5)')
' SD / FINAL PRINT',
n
5409 xpvd(jnucpr)=xpvd(jnucpr) + xpsq(ixspr) + xpsaq(ixspr)
5410 xtvd(jnucta)=xtvd(jnucta) + xtsaq(ixsta) + xtsq(ixsta)
5411 issqq=abs(itsaq(ixsta))
5412 jssqq=abs(itsaq2(ixsta))
5413 IF(issqq.EQ.3.AND.jssqq.EQ.3)
THEN
5415 ELSEIF(issqq.EQ.3.OR.jssqq.EQ.3)
THEN
5430 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5441 COMMON /zsea/zseaav,zseasu,anzsea
5442 common/popcck/pdbck,pdbse,pdbseu,
5443 * ijpock,irejck,ick4,ihad4,ick6,ihad6
5444 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
5445 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
5446 *isea43,isea63,irejao
5450 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
5456 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
5458 COMMON /diqi/ ipvq(248),ippv1(248),ippv2(248), itvq(248),ittv1
5461 +itsaq(
intmx),itsaq2(
intmx),kkproj(248),kktarg(248)
5465 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
5467 * ,xpsu(248),xtsu(248)
5468 * ,xpsut(248),xtsut(248)
5470 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
5471 +ixpv,ixps,ixtv,ixts, intvv1(248),
5472 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
5474 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
5488 COMMON /abrsd/ amcsd1(248),amcsd2(248),gacsd1(248),gacsd2(248),
5489 +bgxsd1(248),bgysd1(248),bgzsd1(248), bgxsd2(248),bgysd2(248),
5490 +bgzsd2(248), nchsd1(248),nchsd2(248),ijcsd1(248),ijcsd2(248),
5491 +pqsda1(248,4),pqsda2(248,4), pqsdb1(248,4),pqsdb2(248,4)
5493 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
5494 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
5504 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
5593 parameter(nfimax=249)
5594 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
5595 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
5596 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
5599 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
5601 COMMON /projk/ iprojk
5603 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
5605 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
5606 * anndv,annvd,annds,annsd,
5608 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
5610 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
5613 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
5614 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
5615 * acouzz,acouhh,acouds,acousd,
5616 * acoudz,acouzd,acoudi,
5617 * acoudv,acouvd,acoucc
5618 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
5619 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
5620 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
5621 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
5622 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
5623 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
5625 dimension poj(4),pat(4)
5628 IF(iphkk.GE.6)
WRITE (6,
'( A)')
' hadrsd'
5632 IF(nchsd1(i).EQ.99.AND.nchsd2(i).EQ.99) go to 50
5636 IF (ipco.GE.6)
WRITE (6,1000) ipsq(is1),ipsaq(is1),itvq(is2),
5637 + ittv1(is2),ittv2(is2), amcsd1(i),amcsd2(i),gacsd1(i),gacsd2(i),
5638 + bgxsd1(i),bgysd1(i),bgzsd1(i), bgxsd2(i),bgysd2(i),bgzsd2(i),
5639 + nchsd1(i),nchsd2(i),ijcsd1(i),ijcsd2(i), pqsda1(i,4),pqsda2
5640 + (i,4),pqsdb1(i,4),pqsdb2(i,4)
5641 1000
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
5651 IF((nchsd1(i).NE.0.OR.nchsd2(i).NE.0).AND.ip.NE.1)
5652 & CALL
saptre(amcsd1(i),gacsd1(i),bgxsd1(i),bgysd1(i),bgzsd1(i),
5653 & amcsd2(i),gacsd2(i),bgxsd2(i),bgysd2(i),bgzsd2(i))
5689 IF(ifb2.LE.2.AND.ifb3.LE.2)
THEN
5691 ELSEIF((ifb2.EQ.3.AND.ifb3.LE.2).OR.
5692 * (ifb3.EQ.3.AND.ifb2.LE.2))
THEN
5694 ELSEIF(ifb2.EQ.3.AND.ifb3.EQ.3)
THEN
5697 IF((nchsd1(i).NE.0))
5698 * CALL
hadjet(nhad,amcsd1(i),poj,pat,gacsd1(i),
5700 + (i),bgzsd1(i),ifb1,ifb2,ifb3,ifb4, ijcsd1(i),
5701 * ijcsd1(i),4,nchsd1
5704 aack=float(ick4)/float(ick4+ihad4+1)
5705 IF((nchsd1(i).EQ.0))
THEN
5706 zseawu=
rndm(bb)*2.d0*zseaav
5707 rseack=float(jitt)*pdbse +zseawu*pdbseu
5708 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,2F10.3)')
'HADJSE JITT,',
5709 *
'RSEACK,PDBSE 4 dpmnuc5 ',
5712 IF(
rndm(v).LE.rseack)
THEN
5714 IF(amcsd1(i).GT.2.3d0)
THEN
5716 CALL
hadjse(nhad,amcsd1(i),poj,pat,gacsd1(i),bgxsd1(i),
5718 + (i),bgzsd1(i),ifb1,ifb2,ifb3,ifb4, ijcsd1(i),ijcsd1(i),4,
5720 + (i),3,irejss,iissqq)
5721 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JITT,',
5722 *
'RSEACK,IREJSS 4 dpmnuc5 NHAD',
5723 + jitt,rseack,irejss,nhad
5726 IF(irejss.EQ.1)irejse=irejse+1
5727 IF(irejss.EQ.3)irejs3=irejs3+1
5728 IF(irejss.EQ.2)irejs0=irejs0+1
5729 CALL
hadjet(nhad,amcsd1(i),poj,pat,gacsd1(i),
5731 + (i),bgzsd1(i),ifb1,ifb2,ifb3,ifb4, ijcsd1(i),
5732 * ijcsd1(i),4,nchsd1
5744 CALL
hadjet(nhad,amcsd1(i),poj,pat,gacsd1(i),
5746 + (i),bgzsd1(i),ifb1,ifb2,ifb3,ifb4, ijcsd1(i),
5747 * ijcsd1(i),4,nchsd1
5757 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRSD: NHKKNMXHKK ',nhkk,
nmxhkk
5762 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
5766 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
5767 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
5776 ptsd=ptsd+
sqrt(pxf(j)**2+pyf(j)**2)
5779 IF(ibarf(j).EQ.500)istist=2
5780 IF(ipco.GE.3)
WRITE(6,*)
' HADRSD before HKKFIL J,NHAD',j,nhad
5782 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),19)
5783 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
5785 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
5786 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
5787 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
5841 IF(ifb2.LE.8.AND.ifb3.LE.8)
THEN
5843 ELSEIF((ifb2.EQ.9.AND.ifb3.LE.8).OR.
5844 * (ifb3.EQ.9.AND.ifb2.LE.8))
THEN
5846 ELSEIF(ifb2.EQ.9.AND.ifb3.EQ.9)
THEN
5849 IF((nchsd2(i).NE.0))
5850 * CALL
hadjet(nhad,amcsd2(i),poj,pat,gacsd2(i),
5852 + (i),bgzsd2(i),ifb1,ifb2,ifb3,ifb4, ijcsd2(i),
5853 * ijcsd2(i),4,nchsd2
5856 aack=float(ick4)/float(ick4+ihad4+1)
5857 IF((nchsd2(i).EQ.0))
THEN
5858 zseawu=
rndm(bb)*2.d0*zseaav
5859 rseack=float(jitt)*pdbse +zseawu*pdbseu
5860 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,2F10.3)')
'HADJASE JITT,',
5864 IF(
rndm(v).LE.rseack)
THEN
5866 IF(amcsd2(i).GT.2.3d0)
THEN
5868 CALL
hadjase(nhad,amcsd2(i),poj,pat,gacsd2(i),bgxsd2(i),
5870 + (i),bgzsd2(i),ifb1,ifb2,ifb3,ifb4, ijcsd2(i),ijcsd2(i),4,
5872 + (i),3,irejss,iissqq)
5873 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JITT,',
5875 + jitt,rseack,irejss
5878 IF(irejss.EQ.1)irejsa=irejsa+1
5879 IF(irejss.EQ.3)ireja3=ireja3+1
5880 IF(irejss.EQ.2)ireja0=ireja0+1
5881 CALL
hadjet(nhad,amcsd2(i),poj,pat,gacsd2(i),
5883 + (i),bgzsd2(i),ifb1,ifb2,ifb3,ifb4, ijcsd2(i),
5884 * ijcsd2(i),4,nchsd2
5896 CALL
hadjet(nhad,amcsd2(i),poj,pat,gacsd2(i),
5898 + (i),bgzsd2(i),ifb1,ifb2,ifb3,ifb4, ijcsd2(i),
5899 * ijcsd2(i),4,nchsd2
5910 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRSD: NHKKNMXHKK ', nhkk,
5916 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
5920 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
5921 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
5930 ptsd=ptsd+
sqrt(pxf(j)**2+pyf(j)**2)
5933 IF(ibarf(j).EQ.500)istist=2
5935 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),20)
5936 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
5938 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
5939 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
5940 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
5951 1010
FORMAT (i6,i4,5i6,9e10.2)
5952 1020
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
5953 1030
FORMAT (
' NHKK,IDHKK(NHKK) ',3i10)
5958 SUBROUTINE diqdzz(ECM,XPSQ1,XPSAQ1,XPSQ2,XPSAQ2,
5959 * ipsq1,ipsaq1,ipsq2,ipsaq2,irejds)
5960 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5964 COMMON /zsea/zseaav,zseasu,anzsea
5965 common/popcck/pdbck,pdbse,pdbseu,
5966 * ijpock,irejck,ick4,ihad4,ick6,ihad6
5967 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
5968 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
5969 *isea43,isea63,irejao
5971 parameter(
intmd=252)
5972 COMMON /intnez/ndz,nzd
5973 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
5989 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
6001 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210),ich(210),
6002 +ibar(210),k1(210),k2(210)
6005 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
6006 *idzre(3),izdre(3),idiqrz(7)
6007 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
6008 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
6009 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
6010 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
6011 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
6012 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
6013 parameter(ummm=0.3d0)
6014 parameter(smmm=0.5d0)
6015 parameter(cmmm=1.3d0)
6033 betcha=betoo+1.3-log10(ecm)
6048 WRITE(6,4567)pc,betcha,pu1,ps1,seasq
6049 4567
FORMAT(
' Charm chain ends DIQDZZ: PC,BETCHA,PU,PS,SEASQ',5f10.5)
6052 IF(iphkk.GE.3)
WRITE (6,
'( A)')
' diqdss'
6056 ipsqq1=1.d0+
rndm(v1)*(2.d0+2.d0*seasq)
6058 IF(rr.LT.pc)ipsqq1=4
6061 amdsq1=xpsq1*xpsq2*ecm**2
6062 amdsq2=xpsaq1*xpsaq2*ecm**2
6063 idiqrz(1)=idiqrz(1)+1
6064 IF(ipsq1.EQ.3.AND.ipsqq1.EQ.3)
THEN
6065 idiqrz(2)=idiqrz(2)+1
6067 IF(amdsq2.LE.6.6d0.OR.amdsq1.LE.6.60d0)
THEN
6068 idiqrz(3)=idiqrz(3)+1
6069 idiqrz(2)=idiqrz(2)-1
6070 idiqrz(1)=idiqrz(1)-1
6074 ELSEIF(ipsq1.EQ.3.OR.ipsqq1.EQ.3)
THEN
6075 idiqrz(4)=idiqrz(4)+1
6077 IF(amdsq2.LE.5.8d0.OR.amdsq1.LE.5.80d0)
THEN
6078 idiqrz(5)=idiqrz(5)+1
6079 idiqrz(4)=idiqrz(4)-1
6080 idiqrz(1)=idiqrz(1)-1
6084 ELSEIF(((ipsq1.EQ.4).OR.(ipsqq1.EQ.4)).AND.
6085 * ((ipsq1.EQ.3).OR.(ipsqq1.EQ.3)))
THEN
6087 IF(amdsq2.LE.30.8d0.OR.amdsq1.LE.30.80d0)
THEN
6091 ELSEIF(ipsq1.EQ.4.OR.ipsqq1.EQ.4)
THEN
6093 IF(amdsq2.LE.25.8.OR.amdsq1.LE.25.80)
THEN
6098 idiqrz(6)=idiqrz(6)+1
6100 IF(amdsq2.LE.3.9.OR.amdsq1.LE.3.9)
THEN
6101 idiqrz(7)=idiqrz(7)+1
6102 idiqrz(6)=idiqrz(6)-1
6103 idiqrz(1)=idiqrz(1)-1
6109 IF(ndz.GE.
intmd)
THEN
6123 IF(iphkk.GE.3)
WRITE (6,
'( A,I10)')
' kkevdz',ndz
6128 IF(iphkk.GE.7)
WRITE(6,
'(A,2I10)')
' KKEVDZ N,NDZ',
n,ndz
6130 prmomz=
sqrt(ecm**2/4.-am(1)**2)
6137 psaqpz=xpsaq1*prmomz
6148 tsdqpz=-xpsaq2*prmomz
6155 pxxx=psqpx + psaqpx + tsqpx + tsdqpx
6156 pyyy=psqpy + psaqpy + tsqpy + tsdqpy
6157 pzzz=psqpz + psaqpz + tsqpz + tsdqpz
6158 eee =psqe + psaqe + tsqe + tsdqe
6159 pptoto=
sqrt(pxxx**2+pyyy**2+pzzz**2)
6160 ammm=
sqrt(abs((eee+pptoto)*(eee-pptoto)))
6161 gammm=eee/(ammm+1.
e-4)
6162 bgggx=pxxx/(ammm+1.
e-4)
6163 bgggy=pyyy/(ammm+1.
e-4)
6164 bgggz=pzzz/(ammm+1.
e-4)
6189 plq1 = xpsq1 *ecm/2.
6191 plaq1= xpsaq1*ecm/2.
6192 eaq1 = xpsaq1*ecm/2.
6193 plq2 =-xpsq2 *ecm/2.
6195 plaq2=-xpsaq2*ecm/2.
6196 eaq2 = xpsaq2*ecm/2.
6199 IF(ipev.GE.6)
WRITE(6,
'(A)')
' KKEVDZ call selpt'
6201 + ptxsq1,ptysq1,plq1,eq1,
6202 + ptxsa1,ptysa1,plaq1,eaq1,
6203 + ptxsa2,ptysa2,plaq2,eaq2,
6204 + ptxsq2,ptysq2,plq2,eq2,
6205 + amch1,amch2,irejds,ikvala,pttq1,ptta1,
6206 + pttq2,ptta2,nselpt)
6208 IF (ipev.GE.7)
WRITE(6,
'(A/5X,5F12.5,I10)')
6209 +
'DS AMMM,GAMMM,BGGGX,BGGGY,BGGGZ,IREJ ', ammm,gammm,bgggx,
6211 IF (irejds.EQ.1)
THEN
6215 WRITE(6,
'(A,I5)')
' KKEVDZ - IRDS13=',irds13
6216 WRITE(6,
'(A/5E12.4/4(4E12.4/),2E12.4/2I5/4E12.4)')
6217 +
' DS: XPSQCM,XPSACM,XTSQCM,XTSACM,AMMM ...', xpsqcm,xpsacm,
6218 + xtsqcm,xtsacm,ammm, ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
6219 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
6220 + amch1,amch2,irejds,ikvala,pttq1,ptta1
6227 ptxch1=ptxsq1 + ptxsq2
6228 ptych1=ptysq1 + ptysq2
6231 ptxch2=ptxsa2 + ptxsa1
6232 ptych2=ptysa2 + ptysa1
6233 ptzch2=plaq2 + plaq1
6238 IF (ipev.GE.6)
WRITE(6,
'(A,5F12.5,I10/A,5F12.5/A,5F12.5)')
6239 +
' DS: AMMM,GAMMM,BGGGX,BGGGY,BGGGZ,IREJDS ', ammm,gammm,bgggx,
6240 + bgggy,bgggz,irejds,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
6241 + amch1,ptxch1,ptych1,ptzch1,ech1,
6242 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
6250 CALL
zobcma(ipsq1,ipsqq1,ipsq2, ijnch1,nnch1,
6251 + irejds,amch1,amch1n,1)
6253 IF(irejds.EQ.1)
THEN
6261 + CALL
zormom(ammm,amch1,amch1n,amch2,
6262 + xpsq1,xpsaq1,xpsaq2,xpsq2,
6263 + ptxsq1,ptysq1,plq1,eq1,
6264 + ptxsa1,ptysa1,plaq1,eaq1,
6265 + ptxsa2,ptysa2,plaq2,eaq2,
6266 + ptxsq2,ptysq2,plq2,eq2,
6267 + ptxch1,ptych1,ptzch1,ech1, ptxch2,ptych2,ptzch2,ech2,
6273 IF (ipev.GE.6)
WRITE(6,
'(A,5F12.5,I10/A,5F12.5/A,5F12.5)')
6274 +
' DS(2): AMMM,GAMMM,BGGGX,BGGGY,BGGGZ,IREJDS',ammm,gammm,bgggx,
6275 + bgggy,bgggz,irejds,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
6276 + amch1,ptxch1,ptych1,ptzch1,ech1,
6277 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
6283 CALL
zobcma(ipsaq1,isaqq1,ipsaq2,
6284 + ijnch2,nnch2,irejds,amch2,amch2n,2)
6286 IF(irejds.EQ.1)
THEN
6290 WRITE(6,1090) irds12
6291 1090
FORMAT(
' KKEVDZ - IRDS12=',i5)
6292 1100
FORMAT(
' DS - 1100', 6i5/2(4e12.4/),2e12.4)
6303 CALL
zorval(ammm,irejds,amch1,amch2, ptxch1,ptych1,ptzch1,ech1,
6304 + ptxch2,ptych2,ptzch2,ech2,iori)
6305 IF(irejds.EQ.1)
THEN
6314 WRITE(6,
'(A/3(1PE15.4),3I5)')
6315 +
' DS - CALL ZORVAL: AMMM,AMCH1,AMCH2,NNCH1,NNCH2,IREJDS',
6316 + ammm, amch1, amch2, nnch1, nnch2, irejds
6317 WRITE(6,1050) ammm,gammm,bgggx,bgggy,bgggz,irejds, amch1,
6318 + ptxch1,ptych1,ptzch1,ech1, amch2,ptxch2,ptych2,ptzch2,ech2
6319 1050
FORMAT (
' DS: AMMM,GAMMM,BGGGX,BGGGY,BGGGZ,IREJDS ',5f12.5,i10/
6320 +
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',5f12.5/
6321 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ',5f12.5)
6323 IF(irejds.EQ.1)
THEN
6344 887
FORMAT(
' DS: QECH1,QECH2,QTZCH1,QTZCH2,AMCH1,AMCH2 ',6f10.2)
6384 gacds1(
n)=qech1/amch1
6385 bgxds1(
n)=qtxch1/amch1
6386 bgyds1(
n)=qtych1/amch1
6387 bgzds1(
n)=qtzch1/amch1
6388 gacds2(
n)=qech2/amch2
6389 bgxds2(
n)=qtxch2/amch2
6390 bgyds2(
n)=qtych2/amch2
6391 bgzds2(
n)=qtzch2/amch2
6396 IF (ipev.GE.3)
WRITE(6,
'(A/I10,4F12.7,5I5/10X,4F12.6/10X,6F12.6,
6397 +4I5/8F15.5/ 8F15.5)')
' DS / FINAL PRINT',
n
6410 IF(issqq.EQ.3.AND.jssqq.EQ.3)
THEN
6412 ELSEIF(issqq.EQ.3.OR.jssqq.EQ.3)
THEN
6433 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6444 COMMON /zsea/zseaav,zseasu,anzsea
6445 common/popcck/pdbck,pdbse,pdbseu,
6446 * ijpock,irejck,ick4,ihad4,ick6,ihad6
6447 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
6448 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
6449 *isea43,isea63,irejao
6451 parameter(
intmd=252)
6452 COMMON /intnez/ ndz,nzd
6467 parameter(
intmx=2488)
6469 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
6475 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
6484 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
6486 * ,xpsu(248),xtsu(248)
6487 * ,xpsut(248),xtsut(248)
6491 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
6492 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
6501 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
6590 parameter(nfimax=249)
6591 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
6592 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
6593 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
6596 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
6598 COMMON /projk/ iprojk
6600 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
6602 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
6603 * anndv,annvd,annds,annsd,
6605 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
6607 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
6610 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
6611 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
6612 * acouzz,acouhh,acouds,acousd,
6613 * acoudz,acouzd,acoudi,
6614 * acoudv,acouvd,acoucc
6615 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
6616 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
6617 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
6618 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
6619 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
6620 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
6622 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
6623 +ixpv,ixps,ixtv,ixts, intvv1(248),
6624 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
6626 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
6629 dimension poj(4),pat(4)
6633 IF(iphkk.GE.3)
WRITE (6,
'( A,4I10)')
' hadrdz',ndz,nzd,
6634 * nchds1(1),nchds2(1)
6638 IF(nchds1(i).EQ.99.AND.nchds2(i).EQ.99) go to 50
6642 IF (ipco.GE.6)
WRITE (6,*)
'IPSQ(IS1),IPSAQ(IS1),ITSQ(IS2)',
6643 +
'ITSAQ(IS2), AMCDS1(I),AMCDS2(I),GACDS1(I),GACDS2(I)',
6644 +
'BGXDS1(I),BGYDS1(I),BGZDS1(I), BGXDS2(I),BGYDS2(I),BGZDS2(I)',
6645 +
'NCHDS1(I),NCHDS2(I),IJCDS1(I),IJCDS2(I), PQDSA1(I,4),PQDSA2',
6646 +
'(I,4),PQDSB1(I,4),PQDSB2(I,4)',
6647 * ipsq(is1),ipsaq(is1),itsq(is2),
6648 + itsaq(is2), amcds1(i),amcds2(i),gacds1(i),gacds2(i),
6649 + bgxds1(i),bgyds1(i),bgzds1(i), bgxds2(i),bgyds2(i),bgzds2(i),
6650 + nchds1(i),nchds2(i),ijcds1(i),ijcds2(i), pqdsa1(i,4),pqdsa2
6651 + (i,4),pqdsb1(i,4),pqdsb2(i,4)
6652 1000
FORMAT(10
x,4i5,10f9.2/10
x,4i5,4f12.4)
6663 IF((nchds1(i).NE.0.OR.nchds2(i).NE.0).AND.ip.NE.1)
6664 & CALL
saptre(amcds1(i),gacds1(i),bgxds1(i),bgyds1(i),bgzds1(i),
6665 & amcds2(i),gacds2(i),bgxds2(i),bgyds2(i),bgzds2(i))
6668 IF (ipco.GE.6)
WRITE (6,1244) poj,pat
6669 1244
FORMAT (
' D-V QUARK-DIQUARK POJ,PAT ',8e12.3)
6683 IF(ipco.GE.4)
WRITE(6,*)
' IPPP,INTVS1(I)',ippp,intvs1(i)
6684 IF(intvs1(i).GT.0)
THEN
6685 ippp = ifrovp(intvs1(i))
6686 IF(ipco.GE.4)
WRITE(6,*)
' IPPP,INTVS1(I)',ippp,intvs1(i)
6688 ELSEIF(intvs1(i).EQ.0)
THEN
6691 IF(ipco.GE.4)
WRITE(6,*)
' JIPP,INTVS1(I)',jipp,intvs1(i)
6699 IF(ifb1.LE.2.AND.ifb2.LE.2)
THEN
6701 ELSEIF((ifb1.EQ.3.AND.ifb2.LE.2).OR.
6702 * (ifb2.EQ.3.AND.ifb1.LE.2))
THEN
6704 ELSEIF(ifb1.EQ.3.AND.ifb2.EQ.3)
THEN
6707 IF((nchds1(i).NE.0))
6708 * CALL
hadjet(nhad,amcds1(i),poj,pat,gacds1(i),bgxds1(i), bgyds1
6709 + (i),bgzds1(i),ifb1,ifb2,ifb3,ifb4, ijcds1(i),ijcds1(i),6,nchds1
6712 aack=float(ick6)/float(ick6+ihad6+1)
6713 IF((nchds1(i).EQ.0))
THEN
6714 zseawu=
rndm(bb)*2.d0*zseaav
6715 rseack=float(jitt)*pdbse +zseawu*pdbseu
6716 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,2F10.3)')
'HADJSE JIPP,',
6717 *
'RSEACK,PDBSE 1 dpmdiqqq',
6720 IF(
rndm(v).LE.rseack)
THEN
6722 IF(amcds1(i).GT.2.3d0)
THEN
6724 CALL
hadjse(nhad,amcds1(i),poj,pat,gacds1(i),bgxds1(i),
6726 + (i),bgzds1(i),ifb1,ifb2,ifb3,ifb4, ijcds1(i),ijcds1(i),6,
6728 + (i),6,irejss,iissqq)
6729 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JIPP,',
6730 *
'RSEACK,IREJSS 1 dpmdiqqq ',
6731 + jipp,rseack,irejss
6734 IF(irejss.EQ.1)irejse=irejse+1
6735 IF(irejss.EQ.3)irejs3=irejs3+1
6736 IF(irejss.EQ.2)irejs0=irejs0+1
6737 CALL
hadjet(nhad,amcds1(i),poj,pat,gacds1(i),bgxds1(i), bgyds1
6738 + (i),bgzds1(i),ifb1,ifb2,ifb3,ifb4, ijcds1(i),ijcds1(i),6,nchds1
6750 CALL
hadjet(nhad,amcds1(i),poj,pat,gacds1(i),bgxds1(i), bgyds1
6751 + (i),bgzds1(i),ifb1,ifb2,ifb3,ifb4, ijcds1(i),ijcds1(i),6,nchds1
6761 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRDZ: NHKKNMXHKK ',nhkk,
nmxhkk
6766 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
6770 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
6771 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
6780 ptdz=ptdz+
sqrt(pxf(j)**2+pyf(j)**2)
6783 IF(ibarf(j).EQ.500)istist=2
6785 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),21)
6786 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
6789 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
6790 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
6791 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
6823 IF(intvs1(i).GT.0)
THEN
6824 ippp = ifrovp(intvs1(i))
6826 ELSEIF(intvs1(i).EQ.0)
THEN
6836 IF(ifb1.LE.8.AND.ifb2.LE.8)
THEN
6838 ELSEIF((ifb1.EQ.9.AND.ifb2.LE.8).OR.
6839 * (ifb2.EQ.9.AND.ifb1.LE.8))
THEN
6841 ELSEIF(ifb1.EQ.9.AND.ifb2.EQ.9)
THEN
6844 IF((nchds2(i).NE.0))
6845 * CALL
hadjet(nhad,amcds2(i),poj,pat,gacds2(i),bgxds2(i), bgyds2
6846 + (i),bgzds2(i),ifb1,ifb2,ifb3,ifb4, ijcds2(i),ijcds2(i),6,nchds2
6849 IF((nchds2(i).EQ.0))
THEN
6850 zseawu=
rndm(bb)*2.d0*zseaav
6851 rseack=float(jitt)*pdbse +zseawu*pdbseu
6852 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,2F10.3)')
'HADJSE JIPP,',
6856 IF(
rndm(v).LE.rseack)
THEN
6858 IF(amcds2(i).GT.2.3d0)
THEN
6860 CALL
hadjase(nhad,amcds2(i),poj,pat,gacds2(i),bgxds2(i),
6862 + (i),bgzds2(i),ifb1,ifb2,ifb3,ifb4, ijcds2(i),ijcds2(i),6,
6864 + (i),6,irejss,iissqq)
6865 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JIPP,',
6867 + jipp,rseack,irejss
6870 IF(irejss.EQ.1)irejsa=irejsa+1
6871 IF(irejss.EQ.3)ireja3=ireja3+1
6872 IF(irejss.EQ.2)ireja0=ireja0+1
6873 CALL
hadjet(nhad,amcds2(i),poj,pat,gacds2(i),bgxds2(i), bgyds2
6874 + (i),bgzds2(i),ifb1,ifb2,ifb3,ifb4, ijcds2(i),ijcds2(i),6,nchds2
6886 CALL
hadjet(nhad,amcds2(i),poj,pat,gacds2(i),bgxds2(i), bgyds2
6887 + (i),bgzds2(i),ifb1,ifb2,ifb3,ifb4, ijcds2(i),ijcds2(i),6,nchds2
6898 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRDZ: NHKKNMXHKK ', nhkk,
6904 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
6908 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
6909 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
6919 ptdz=ptdz+
sqrt(pxf(j)**2+pyf(j)**2)
6922 IF(ibarf(j).EQ.500)istist=2
6924 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),22)
6925 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
6928 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
6929 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
6930 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
6942 1010
FORMAT (i6,i4,5i6,9e10.2)
6943 1020
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
6944 1030
FORMAT (
' NHKK,IDHKK(NHKK) ',3i10)
6951 SUBROUTINE diqzzd(ECM,XPSQ1,XPSAQ1,XPSQ2,XPSAQ2,
6952 * ipsq1,ipsaq1,ipsq2,ipsaq2,irejsd)
6953 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6957 COMMON /zsea/zseaav,zseasu,anzsea
6958 common/popcck/pdbck,pdbse,pdbseu,
6959 * ijpock,irejck,ick4,ihad4,ick6,ihad6
6960 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
6961 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
6962 *isea43,isea63,irejao
6964 parameter(
intmd=252)
6965 COMMON /intnez/ ndz,nzd
6966 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
6979 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210),ich(210),
6980 +ibar(210),k1(210),k2(210)
6995 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
6997 COMMON /diqrej/idiqre(7),idvre(3),ivdre(3),idsre(3),isdre(3),
6998 *idzre(3),izdre(3),idiqrz(7)
6999 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
7000 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
7001 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
7002 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
7003 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
7004 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
7006 parameter(ummm=0.3d0)
7007 parameter(smmm=0.5d0)
7008 parameter(cmmm=1.3d0)
7026 betcha=betoo+1.3-log10(ecm)
7041 WRITE(6,4567)pc,betcha,pu1,ps1,seasq
7042 4567
FORMAT(
' Charm chain ends DIQZZD: PC,BETCHA,PU,PS,SEASQ',5f10.5)
7045 IF(iphkk.GE.3)
WRITE (6,
'( A)')
' diqssd'
7053 ipsqq2=1.d0+
rndm(v1)*(2.d0+2.d0*seasq)
7055 IF(rr.LT.pc)ipsqq2=4
7057 amsdq1=xpsq2*xpsq1*ecm**2
7058 amsdq2=xpsaq2*xpsaq1*ecm**2
7059 idiqrz(1)=idiqrz(1)+1
7061 IF(ipsq2.EQ.3.AND.ipsqq2.EQ.3)
THEN
7062 idiqrz(2)=idiqrz(2)+1
7064 IF(amsdq2.LE.6.6d0.OR.amsdq1.LE.6.6d0)
THEN
7065 idiqrz(3)=idiqrz(3)+1
7066 idiqrz(2)=idiqrz(2)-1
7067 idiqrz(1)=idiqrz(1)-1
7071 ELSEIF(ipsq2.EQ.3.OR.ipsqq2.EQ.3)
THEN
7072 idiqrz(4)=idiqrz(4)+1
7074 IF(amsdq2.LE.5.8d0.OR.amsdq1.LE.5.80d0)
THEN
7075 idiqrz(5)=idiqrz(5)+1
7076 idiqrz(4)=idiqrz(4)-1
7077 idiqrz(1)=idiqrz(1)-1
7081 ELSEIF(((ipsq2.EQ.4).OR.(ipsqq2.EQ.4)).AND.
7082 * ((ipsq2.EQ.3).OR.(ipsqq2.EQ.3)))
THEN
7084 IF(amsdq2.LE.30.8d0.OR.amsdq1.LE.30.80d0)
THEN
7088 ELSEIF(ipsq2.EQ.4.OR.ipsqq2.EQ.4)
THEN
7090 IF(amsdq2.LE.25.8d0.OR.amsdq1.LE.25.80d0)
THEN
7095 idiqrz(6)=idiqrz(6)+1
7097 IF(amsdq2.LE.3.9d0.OR.amsdq1.LE.3.9d0)
THEN
7098 idiqrz(7)=idiqrz(7)+1
7099 idiqrz(6)=idiqrz(6)-1
7100 idiqrz(1)=idiqrz(1)-1
7106 IF(nzd.GE.
intmd)
THEN
7118 IF(iphkk.GE.3)
WRITE (6,
'( A,3I10)')
' kkevzd',nzd,
7119 * nchsd1(1),nchsd2(1)
7123 IF(nchsd1(
n).EQ.99.AND.nchsd2(
n).EQ.99)go to 10
7127 prmomz=
sqrt(ecm**2/4.-am(1)**2)
7134 psdqpz=xpsaq1*prmomz
7145 tsaqpz=-xpsaq2*prmomz
7152 pxxx=tsqpx + tsaqpx + psqpx + psdqpx
7153 pyyy=tsqpy + tsaqpy + psqpy + psdqpy
7154 pzzz=tsqpz + tsaqpz + psqpz + psdqpz
7155 eee =tsqe + tsaqe + psqe + psdqe
7156 pptoto=
sqrt(pxxx**2+pyyy**2+pzzz**2)
7157 ammm=
sqrt(abs((eee+pptoto)*(eee-pptoto)))
7158 gammm=eee/(ammm+1.
e-4)
7159 bgggx=pxxx/(ammm+1.
e-4)
7160 bgggy=pyyy/(ammm+1.
e-4)
7161 bgggz=pzzz/(ammm+1.
e-4)
7182 plq1 = xpsq1 *ecm/2.
7184 plaq1= xpsaq1*ecm/2.
7185 eaq1 = xpsaq1*ecm/2.
7186 plq2 =-xpsq2 *ecm/2.
7188 plaq2=-xpsaq2*ecm/2.
7189 eaq2 = xpsaq2*ecm/2.
7192 IF(ipev.GE.6)
WRITE(6,
'(A)')
' KKEVZD call selpt'
7193 CALL
selpt( ptxsq1,ptysq1,plq1,
7194 + eq1,ptxsa1,ptysa1,plaq1,eaq1,
7195 + ptxsa2,ptysa2,plaq2,eaq2,
7196 + ptxsq2,ptysq2,plq2,eq2,
7197 + amch1,amch2,irejsd,ikvala,pttq1,ptta1,
7198 + pttq2,ptta2,nselpt)
7204 IF (ipev.GE.6)
WRITE(6,
'(A/5X,5F12.5,I10)')
7205 +
'SD AMMM,GAMMM,BGGGX,BGGGY,BGGGZ,IREJSD ', ammm,gammm,bgggx,
7206 + bgggy,bgggz,irejsd
7207 IF (irejsd.EQ.1)
THEN
7210 WRITE(6,
'(A,I5)')
' KKEVZD - IRSD13=',irsd13
7211 WRITE(6,
'(A/5E12.4/4(4E12.4/),2E12.4/2I5/4E12.4)')
7212 +
' VD: XPVQCM,XPVDCM,XTSQCM,XTSACM,AMMM ...', xpsqcm,xpsdcm,
7213 + xtsqcm,xtsacm,ammm, ptxsq1,ptysq1,plq1,eq1,ptxsa1,ptysa1,
7214 + plaq1,eaq1, ptxsq2,ptysq2,plq2,eq2,ptxsa2,ptysa2,plaq2,eaq2,
7215 + amch1,amch2,irejsd,ikvala,pttq1,ptta1
7223 ptxch1=ptxsq1 + ptxsq2
7224 ptych1=ptysq1 + ptysq2
7227 ptxch2=ptxsa2 + ptxsa1
7228 ptych2=ptysa2 + ptysa1
7229 ptzch2=plaq2 + plaq1
7234 IF (ipev.GE.6)
WRITE(6,
'(A,5F12.5,I10/A,5F12.5/A,5F12.5)')
7235 +
' SD: AMMM,GAMMM,BGGGX,BGGGY,BGGGZ,IREJSD ', ammm,gammm,bgggx,
7236 + bgggy,bgggz,irejsd,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
7237 + amch1,ptxch1,ptych1,ptzch1,ech1,
7238 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
7246 CALL
zobcma(ipsq2,ipsqq2,ipsq1, ijnch1,nnch1,
7247 + irejsd,amch1,amch1n,1)
7249 IF(irejsd.EQ.1)
THEN
7256 + CALL
zormom(ammm,amch1,amch1n,amch2,
7257 + xpsq1,xpsaq1,xpsaq2,xpsq2,
7258 + ptxsq1,ptysq1,plq1,eq1,
7259 + ptxsa1,ptysa1,plaq1,eaq1,
7260 + ptxsa2,ptysa2,plaq2,eaq2,
7261 + ptxsq2,ptysq2,plq2,eq2,
7262 + ptxch1,ptych1,ptzch1,ech1, ptxch2,ptych2,ptzch2,ech2,
7268 IF (ipev.GE.6)
WRITE(6,
'(A,5F12.5,I10/A,5F12.5/A,5F12.5)')
7269 +
' SD(2): AMMM,GAMMM,BGGGX,BGGGY,BGGGZ,IREJSD',ammm,gammm,bgggx,
7270 + bgggy,bgggz,irejsd,
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',
7271 + amch1,ptxch1,ptych1,ptzch1,ech1,
7272 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ', amch2,ptxch2,ptych2,
7278 CALL
zobcma(ipsaq1,ipsaq2,isaqq2,
7279 + ijnch2,nnch2,irejsd,amch2,amch2n,2)
7281 IF(irejsd.EQ.1)
THEN
7284 WRITE(6,1090) irsd12
7285 WRITE(6,1100) ipsaq1,ipsaq2,isaqq2,
7286 + ijnch2,nnch2,irejsd,
7287 + xpsq1,xpsaq1,xpsqcm,xtsacm, xpsq2,xpsaq2
7288 + ,xtsqcm,xtsacm, amch2,amch2n
7289 1090
FORMAT(
' KKEVZD - IRSD12=',i5)
7290 1100
FORMAT(
' SD - 1100', 6i5/2(4e12.4/),2e12.4)
7301 CALL
zorval(ammm,irejsd,amch1,amch2, ptxch1,ptych1,ptzch1,ech1,
7302 + ptxch2,ptych2,ptzch2,ech2,iori)
7305 WRITE(6,
'(A/3(1PE15.4),3I5)')
7306 +
' SD - CALL ZORVAL: AMMM,AMCH1,AMCH2,NNCH1,NNCH2,IREJSD',
7307 + ammm, amch1, amch2, nnch1, nnch2, irejsd
7308 WRITE(6,1050) ammm,gammm,bgggx,bgggy,bgggz,irejsd, amch1,
7309 + ptxch1,ptych1,ptzch1,ech1, amch2,ptxch2,ptych2,ptzch2,ech2
7310 1050
FORMAT (
' SD: AMMM,GAMMM,BGGGX,BGGGY,BGGGZ,IREJSD ',5f12.5,i10/
7311 +
' AMCH1,PTXCH1,PTYCH1,PTZCH1,ECH1 ',5f12.5/
7312 +
' AMCH2,PTXCH2,PTYCH2,PTZCH2,ECH2 ',5f12.5)
7314 IF(irejsd.EQ.1)
THEN
7375 gacsd1(
n)=qech1/amch1
7376 bgxsd1(
n)=qtxch1/amch1
7377 bgysd1(
n)=qtych1/amch1
7378 bgzsd1(
n)=qtzch1/amch1
7379 gacsd2(
n)=qech2/amch2
7380 bgxsd2(
n)=qtxch2/amch2
7381 bgysd2(
n)=qtych2/amch2
7382 bgzsd2(
n)=qtzch2/amch2
7387 IF (ipev.GE.2)
WRITE(6,
'(A/I10,4F12.7,5I5/10X,4F12.6/10X,6F12.6,
7388 +4I5/8F15.5/8F15.5/2I5)')
' SD / FINAL PRINT',
n
7401 IF(issqq.EQ.3.AND.jssqq.EQ.3)
THEN
7403 ELSEIF(issqq.EQ.3.OR.jssqq.EQ.3)
THEN
7423 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7434 COMMON /zsea/zseaav,zseasu,anzsea
7435 common/popcck/pdbck,pdbse,pdbseu,
7436 * ijpock,irejck,ick4,ihad4,ick6,ihad6
7437 *,irejse,ise4,ise6,irejs3,ise43,ise63,irejs0
7438 *,ihada4,ihada6,irejsa,isea4,isea6,ireja3,
7439 *isea43,isea63,irejao
7440 parameter(
intmd=252)
7441 COMMON /intnez/ ndz,nzd
7456 parameter(
intmx=2488)
7458 COMMON /ifroto/ ifrovp(248),itovp(248),ifrosp(
intmx), ifrovt(248),
7464 +mhkkdv(248),mhkkvd(248), mhkkds(
intmd),mhkksd(
intmd)
7473 COMMON /dxqx/ xpvq(248),xpvd(248),xtvq(248),xtvd(248), xpsq
7475 * ,xpsu(248),xtsu(248)
7476 * ,xpsut(248),xtsut(248)
7496 LOGICAL zuovp,zuosp,zuovt,zuost,intlo,inloss
7497 COMMON /lozuo/ zuovp(248),zuosp(
intmx),zuovt(248),zuost(
intmx),
7507 COMMON /hkkevt/ nhkk,nevhkk,isthkk(
nmxhkk),idhkk(
nmxhkk), jmohkk
7596 parameter(nfimax=249)
7597 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
7598 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
7599 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
7602 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
7604 COMMON /projk/ iprojk
7606 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
7608 COMMON /bufueh/ annvv,annss,annsv,annvs,anncc,
7609 * anndv,annvd,annds,annsd,
7611 * ptvv,ptss,ptsv,ptvs,ptcc,ptdv,ptvd,ptds,ptsd,
7613 * eevv,eess,eesv,eevs,eecc,eedv,eevd,eeds,eesd,
7616 * ,annzd,anndz,ptzd,ptdz,eezd,eedz
7617 COMMON /ncouch/ acouvv,acouss,acousv,acouvs,
7618 * acouzz,acouhh,acouds,acousd,
7619 * acoudz,acouzd,acoudi,
7620 * acoudv,acouvd,acoucc
7621 COMMON /diqsum/ndvuu,ndvus,ndvss,nvduu,nvdus,nvdss,
7622 * ndsuu,ndsus,ndsss,nsduu,nsdus,nsdss,
7623 * ndzuu,ndzus,ndzss,nzduu,nzdus,nzdss
7624 * ,nadvuu,nadvus,nadvss,navduu,navdus,navdss,
7625 * nadsuu,nadsus,nadsss,nasduu,nasdus,nasdss,
7626 * nadzuu,nadzus,nadzss,nazduu,nazdus,nazdss
7629 COMMON /intnew/ nvv,nsv,nvs,nss,ndv,nvd,nds,nsd,
7630 +ixpv,ixps,ixtv,ixts, intvv1(248),
7631 +intvv2(248),intsv1(248),intsv2(248), intvs1(248),intvs2(248),
7633 +intdv1(248),intdv2(248),intvd1(248),intvd2(248),
7636 dimension poj(4),pat(4)
7640 IF(iphkk.GE.3)
WRITE (6,
'( A,4I10)')
' hadrzd',ndz,nzd,
7641 * nchsd1(1),nchsd2(1)
7645 IF(nchsd1(i).EQ.99.AND.nchsd2(i).EQ.99) go to 50
7654 1000
FORMAT(10
x,5i5,10f9.2/10
x,4i5,4f12.4)
7664 IF((nchsd1(i).NE.0.OR.nchsd2(i).NE.0).AND.ip.NE.1)
7665 & CALL
saptre(amcsd1(i),gacsd1(i),bgxsd1(i),bgysd1(i),bgzsd1(i),
7666 & amcsd2(i),gacsd2(i),bgxsd2(i),bgysd2(i),bgzsd2(i))
7684 IF(intsv2(i).GT.0)
THEN
7685 ittt = ifrovt(intsv2(i))
7687 ELSEIF(intsv2(i).EQ.0)
THEN
7697 IF(ifb2.LE.2.AND.ifb3.LE.2)
THEN
7699 ELSEIF((ifb2.EQ.3.AND.ifb3.LE.2).OR.
7700 * (ifb3.EQ.3.AND.ifb2.LE.2))
THEN
7702 ELSEIF(ifb2.EQ.3.AND.ifb3.EQ.3)
THEN
7705 IF((nchsd1(i).NE.0))
7706 * CALL
hadjet(nhad,amcsd1(i),poj,pat,gacsd1(i),bgxsd1(i), bgysd1
7707 + (i),bgzsd1(i),ifb1,ifb2,ifb3,ifb4, ijcsd1(i),ijcsd1(i),4,nchsd1
7710 aack=float(ick4)/float(ick4+ihad4+1)
7711 IF((nchsd1(i).EQ.0))
THEN
7712 zseawu=
rndm(bb)*2.d0*zseaav
7713 rseack=float(jitt)*pdbse +zseawu*pdbseu
7714 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,2F10.3)')
'HADJSE JITT,',
7715 *
'RSEACK,PDBSE 2 dpmdiqqq ',
7718 IF(
rndm(v).LE.rseack)
THEN
7720 IF(amcsd1(i).GT.2.3d0)
THEN
7722 CALL
hadjse(nhad,amcsd1(i),poj,pat,gacsd1(i),bgxsd1(i),
7724 + (i),bgzsd1(i),ifb1,ifb2,ifb3,ifb4, ijcsd1(i),ijcsd1(i),4,
7726 + (i),3,irejss,iissqq)
7727 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JITT,',
7728 *
'RSEACK,IREJSS 2 dpmdiqqq ',
7729 + jitt,rseack,irejss
7732 IF(irejss.EQ.1)irejse=irejse+1
7733 IF(irejss.EQ.3)irejs3=irejs3+1
7734 IF(irejss.EQ.2)irejs0=irejs0+1
7735 CALL
hadjet(nhad,amcsd1(i),poj,pat,gacsd1(i),bgxsd1(i), bgysd1
7736 + (i),bgzsd1(i),ifb1,ifb2,ifb3,ifb4, ijcsd1(i),ijcsd1(i),4,nchsd1
7748 CALL
hadjet(nhad,amcsd1(i),poj,pat,gacsd1(i),bgxsd1(i), bgysd1
7749 + (i),bgzsd1(i),ifb1,ifb2,ifb3,ifb4, ijcsd1(i),ijcsd1(i),4,nchsd1
7759 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRZD: NHKKNMXHKK ',nhkk,
nmxhkk
7764 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
7768 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
7769 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
7778 ptzd=ptzd+
sqrt(pxf(j)**2+pyf(j)**2)
7781 IF(ibarf(j).EQ.500)istist=2
7783 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),23)
7784 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
7787 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
7788 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
7789 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
7821 IF(intsv2(i).GT.0)
THEN
7822 ittt = ifrovt(intsv2(i))
7824 ELSEIF(intsv2(i).EQ.0)
THEN
7834 IF(ifb2.LE.8.AND.ifb3.LE.8)
THEN
7836 ELSEIF((ifb2.EQ.9.AND.ifb3.LE.8).OR.
7837 * (ifb3.EQ.9.AND.ifb2.LE.8))
THEN
7839 ELSEIF(ifb2.EQ.9.AND.ifb3.EQ.9)
THEN
7842 IF((nchsd2(i).NE.0))
7843 * CALL
hadjet(nhad,amcsd2(i),poj,pat,gacsd2(i),bgxsd2(i), bgysd2
7844 + (i),bgzsd2(i),ifb1,ifb2,ifb3,ifb4, ijcsd2(i),ijcsd2(i),4,nchsd2
7847 IF((nchsd1(i).EQ.0))
THEN
7848 zseawu=
rndm(bb)*2.d0*zseaav
7849 rseack=float(jitt)*pdbse +zseawu*pdbseu
7850 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,2F10.3)')
'HADJSE JITT,',
7854 IF(
rndm(v).LE.rseack)
THEN
7856 IF(amcsd2(i).GT.2.3d0)
THEN
7858 CALL
hadjase(nhad,amcsd2(i),poj,pat,gacsd2(i),bgxsd2(i),
7860 + (i),bgzsd2(i),ifb1,ifb2,ifb3,ifb4, ijcsd2(i),ijcsd2(i),4,
7862 + (i),3,irejss,iissqq)
7863 IF(ipco.GE.1)
WRITE(6,
'(2A,I5,F10.3,I5)')
'HADJSE JITT,',
7865 + jitt,rseack,irejss
7868 IF(irejss.EQ.1)irejsa=irejsa+1
7869 IF(irejss.EQ.3)ireja3=ireja3+1
7870 IF(irejss.EQ.2)ireja0=ireja0+1
7871 CALL
hadjet(nhad,amcsd2(i),poj,pat,gacsd2(i),bgxsd2(i), bgysd2
7872 + (i),bgzsd2(i),ifb1,ifb2,ifb3,ifb4, ijcsd2(i),ijcsd2(i),4,nchsd2
7884 CALL
hadjet(nhad,amcsd2(i),poj,pat,gacsd2(i),bgxsd2(i), bgysd2
7885 + (i),bgzsd2(i),ifb1,ifb2,ifb3,ifb4, ijcsd2(i),ijcsd2(i),4,nchsd2
7896 WRITE (6,
'(A,2I5/A)') .EQ.
' HADRZD: NHKKNMXHKK ', nhkk,
7902 WRITE (6,
'(A,2I5)').EQ.
' XKSAMP:NHKKNMXHKK ',nhkk,
nmxhkk
7906 ehecc=
sqrt(pxf(j)**2+pyf(j)**2+pzf(j)**2+amf(j)**2)
7907 IF (abs(ehecc-hef(j)).GT.0.001)
THEN
7916 ptzd=ptzd+
sqrt(pxf(j)**2+pyf(j)**2)
7919 IF(ibarf(j).EQ.500)istist=2
7921 * pxf(j),pyf(j),pzf(j),hef(j),nhkkau,iormo(j),24)
7922 IF(idhkk(nhkk).EQ.99999)
WRITE (6,1030)nhkk,nref(j), idhkk
7925 IF (iphkk.GE.2)
WRITE(6,1010) nhkk, isthkk(nhkk),idhkk(nhkk),
7926 + jmohkk(1,nhkk),jmohkk(2,nhkk), jdahkk(1,nhkk),jdahkk(2,nhkk),
7927 + (phkk(khkk,nhkk),khkk=1,5), (vhkk(khkk,nhkk),khkk=1,4)
7939 1010
FORMAT (i6,i4,5i6,9e10.2)
7940 1020
FORMAT (.GT.
' HADRKK JNAUMAX SKIP NEXT PARTICLES ',3i10)
7941 1030
FORMAT (
' NHKK,IDHKK(NHKK) ',3i10)
7947 SUBROUTINE zobcma(IF1,IF2,IF3,IJNCH,NNCH,IREJ,AMCH,AMCHN,IKET)
7948 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7970 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
7971 +iibar(210),k1(210),k2(210)
7973 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
7976 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
7978 COMMON /ketmas/ am8(2),am10(2),ib88(2),ib1010(2),amch1n,amch2n
7981 CALL
dbklas(if1,if2,if3,ib8,ibb10)
7983 IF (ipev.GE.6)
WRITE(6,1000)if1,if2,if3,ib8,ibb10
7984 1000
FORMAT (
' COBCMA: IPQ,ITTQ1,ITTQ2,IB8,IBB10 ',5i5)
7997 IF(amch.LT.am81)
THEN
7999 ELSEIF (amch.LT.am101)
THEN
8005 ELSEIF(amch.LT.amff1)
THEN
8016 WRITE(6,1010) amch,amchn,am81,am101
8017 WRITE(6,1020) if1,if2,if3,ib8,ibb10,ijnch,nnch,irej
8018 1010
FORMAT(
' COBCMA: AMCH,AMCHN,AM81,AM101', 4f13.4)
8019 1020
FORMAT(
' COBCMA: IF1,IF2,IF3,IB8,IBB10,IJNCH,NNCH,IREJ',8i4)
8028 SUBROUTINE zomcma(IFQ,IFAQ,IJNCH,NNCH,IREJ,AMCH,AMCHN)
8029 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
8046 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
8047 +iibar(210),k1(210),k2(210)
8051 COMMON /inpdat/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
8052 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
8054 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
8056 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
8059 ifps=imps(iifaq,ifq)
8061 IF (ipev.GE.6)
WRITE (6,1000)iifaq,ifq,ifps,ifv
8062 1000
FORMAT (
' COMCMA',5
x,
' IIPPAQ,ITQ,IFPS,IFV ',4i5)
8069 IF(ipev.GE.6)
WRITE(6,1010) amch,amps,amv,ifps,ifv
8070 1010
FORMAT(
' AMCH,AMPS,AMV,IFPS,IFV ',3f12.4,2i10)
8072 IF(amch.LT.amps)
THEN
8077 IF (amch.LT.amv)
THEN
8082 ELSEIF(amch.LT.amff)
THEN
8099 SUBROUTINE zomcm2(IQ1,IQ2,IAQ1,IAQ2,NNCH,IREJ,AMCH)
8100 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
8122 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
8123 +iibar(210),k1(210),k2(210)
8127 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
8129 COMMON /inpdat/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
8130 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
8132 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
8137 IF (iiaq1.EQ.iq1) go to 10
8138 IF (iiaq1.EQ.iq2) go to 20
8139 IF (iiaq2.EQ.iq1) go to 30
8140 IF (iiaq2.EQ.iq2) go to 40
8145 WRITE(6,
'(A/5X,4I5,1PE13.5)')
8146 +
' KKEVVV/COMCM2 (QU. NUMBERS): IQ1, IQ2, IAQ1, IAQ2, AMCH', iq1,
8147 + iq2, iaq1, iaq2, amch
8175 IF (amch.LT.amff)
THEN
8178 WRITE(6,
'(A/5X,4I5,1PE13.5)')
8179 +
' KKEVVV/COMCM2 (MASS!): IQ1, IQ2, IAQ1, IAQ2, AMCH', iq1,
8180 + iq2, iaq1, iaq2, amch
8190 SUBROUTINE zormom(AMMM,AMCH1,AMCH1N,AMCH2N, XP,XPP,XTVQ,XTVD,
8191 +pq1x,pq1y,pq1z,pq1e,pa1x,pa1y,pa1z,pa1e, pq2x,pq2y,pq2z,pq2e,pa2x,
8192 +pa2y,pa2z,pa2e, pxch1,pych1,pzch1,ech1, pxch2,pych2,pzch2,ech2,
8194 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
8206 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
8208 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
8215 xpp=xpp + xpsqol - xp
8218 xtvq=xtvq + xtvdol - xtvd
8220 xtvqcm=xtvq/(xtvq+xtvd)
8237 IF(ech1.LE.amch1)
THEN
8242 pch1 =
sqrt(abs((ech1-amch1)*(ech1+amch1)))+0.000001
8251 IF(ech2.LE.pch1)
THEN
8256 amch2n=
sqrt(abs((ech2-pch1)*(ech2+pch1)))
8260 IF(pch1.GT.(pa1e+pq2e))
THEN
8267 ct1=-(pch1**2 + (pa1e-pq2e)*(pa1e+pq2e))/(2.0*pch1*pa1e)
8268 if(abs(ct1).gt.1.0)
then
8271 ct1=sign(0.999999999,ct1)
8276 st1=
sqrt(abs((1.0+ct1)*(1.0-ct1)))
8278 CALL
drtran(cxch1,cych1,czch1,ct1,st1,sfe,cfe,cxa1,cya1,cza1)
8287 pxsum=pq1x+pa1x+pq2x+pa2x
8288 pysum=pq1y+pa1y+pq2y+pa2y
8289 pzsum=pq1z+pa1z+pq2z+pa2z
8290 pesum=pq1e+pa1e+pq2e+pa2e
8291 WRITE(6,
'(A)')
' ZORMOM: KINEMATIC TEST FOR PARTONS'
8292 WRITE(6,
'(A,1PE12.5)')
' AMMM',ammm
8293 WRITE(6,
'(A,4(1PE12.5))')
' PXSUM,PYSUM,PZSUM,PESUM', pxsum,
8304 SUBROUTINE zorval(AMMM,IREJ,AMCH1,AMCH2, QTX1,QTY1,QZ1,QE1,QTX2,
8306 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
8315 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
8318 COMMON /outlev/ioutpo,ioutpa,iouxev,ioucol
8320 IF(ammm.LE.amch1+amch2+0.4)
THEN
8325 ek1=(ammm**2-amch2**2 + amch1**2)/(2.*ammm)
8327 pzk1=
sqrt(ek1**2 - amch1**2)
8329 pzk2=
sqrt(ek2**2 - amch2**2)
8338 bgx=(qtx1+qtx2)/ammm
8339 bgy=(qty1+qty2)/ammm
8342 IF(abs(gam-1.).GT.1
e-3)
THEN
8351 CALL
daltra(gam,-bgx,-bgy,-bgz,pxk1,pyk1,pzk1,ek1,pppch1, qtx1,
8353 CALL
daltra(gam,-bgx,-bgy,-bgz,pxk2,pyk2,pzk2,ek2,pppch2, qtx2,
8356 WRITE(6,
'(2A)')
' ZORVAL - CORRECTION OF CHAIN MOMENTA',
8357 +
' IF MASS OF CHAIN 2 HAD TO BE CHANGED'
subroutine kkevdv(IREJDV)
subroutine hadjse(NHAD, AMCH, PPR, PTA, GAM, BGX, BGY, BGZ, IFB1, IFB2, IFB3, IFB4, I1, I2, NOBAM, NNCH, NORIG, IREJ, IISSQQ)
subroutine kkevvd(IREJVD)
subroutine cromsc(PX, PY, PZ, E, RX, RY, RZ, PXN, PYN, PZN, EN, IORIG)
subroutine saptre(AM1, G1, BGX1, BGY1, BGZ1, AM2, G2, BGX2, BGY2, BGZ2)
DOUBLE PRECISION function rndm(RDUMMY)
subroutine diqzzd(ECM, XPSQ1, XPSAQ1, XPSQ2, XPSAQ2, IPSQ1, IPSAQ1, IPSQ2, IPSAQ2, IREJSD)
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 zomcma(IFQ, IFAQ, IJNCH, NNCH, IREJ, AMCH, AMCHN)
subroutine zormom(AMMM, AMCH1, AMCH1N, AMCH2N, XP, XPP, XTVQ, XTVD, 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 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 kkevds(IREJDS)
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 diqdss(ECM, ITS, IPS, IREJ)
subroutine diqsv(ECM, ITV, J, IREJ)
subroutine zobcma(IF1, IF2, IF3, IJNCH, NNCH, IREJ, AMCH, AMCHN, IKET)
subroutine dsfecf(SFE, CFE)
subroutine cobcma(IF1, IF2, IF3, IJNCH, NNCH, IREJ, AMCH, AMCHN, IKET)
subroutine zomcm2(IQ1, IQ2, IAQ1, IAQ2, NNCH, IREJ, AMCH)
G4double tau(G4double bet, G4double homega, G4double ef, G4double t)
static c2_sqrt_p< float_type > & sqrt()
make a *new object
subroutine drtran(XO, YO, ZO, CDE, SDE, SFE, CFE, X, Y, Z)
subroutine hkkfil(IST, ID, M1, M2, PX, PY, PZ, E, NHKKAU, KORMO, ICALL)
subroutine zorval(AMMM, IREJ, AMCH1, AMCH2, QTX1, QTY1, QZ1, QE1, QTX2, QTY2, QZ2, QE2, IORI)
subroutine kkevsd(IREJSD)
DOUBLE PRECISION function dbeta(X1, X2, BET)
subroutine dbklas(I, J, K, I8, I10)
subroutine diqdzz(ECM, XPSQ1, XPSAQ1, XPSQ2, XPSAQ2, IPSQ1, IPSAQ1, IPSQ2, IPSAQ2, IREJDS)
subroutine hadjase(NHAD, AMCH, PPR, PTA, GAM, BGX, BGY, BGZ, IFB1, IFB2, IFB3, IFB4, I1, I2, NOBAM, NNCH, NORIG, IREJ, IISSQQ)
subroutine hadjet(NHAD, AMCH, PPR, PTA, GAM, BGX, BGY, BGZ, IFB1, IFB2, IFB3, IFB4, I1, I2, NOBAM, NNCH, NORIG)