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'