63 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
69 common/
pyjets/
n,npad,k(4000,5),
p(4000,5),v(4000,5)
70 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
71 common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
72 common/
pydat3/mdcy(500,3),mdme(4000,2),brat(4000),kfdp(4000,5)
73 common/popcor/pdb,ajsdef
77 COMMON /harlun/ qlun,iharlu
111 IF(ipromu.NE.0)ipromm=ipromu
115 IF(ajsdef.EQ.1.d0) go to 100
119 IF(pdb.EQ.0.d0)mstj(12)=1
126 IF(pdb.GT.0.d0)parj(5)=pdb
194 ELSEIF(ifrag.EQ.10)
THEN
249 WRITE (6,2355)pdb,ajsdef,mstj(12),parj(42),parj(21)
250 2355
FORMAT(
' LUNDIN initialization PDB,AJSDEF= ',2f10.3/
251 +
' MSTJ(12) popcorn default=2 : ',i10/
252 +
' PARJ(42 )Lund b,default=0.9 : ',f10.3/
253 +
' PARJ(21) sigma in pt distr,default=0.35 : ',f10.3)
310 IF(ipromm.EQ.1)go to 199
371 SUBROUTINE bamlun(IHAD,KFA1,KFA2,KFA3,KFA4,AEO,IOPT,IREJ)
378 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
383 common/
pyjets/
n,npad,k(4000,5),
p(4000,5),v(4000,5)
384 common/
pydat1/mstu(200),paru(200),mstj(200),parj(200)
385 common/pydat2/kchg(500,4),pmas(500,4),parf(2000),vckm(4,4)
387 parameter(nfimax=249)
388 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
389 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
390 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
393 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
395 common/dpar/aname(210),am(210),ga(210),
tau(210),ich(210),
396 *ibar(210),k1(210),k2(210)
397 common/diffra/isingd,idiftp,ioudif,iflagd
398 common/capkop/xxx1,xxx3
400 dimension koriii(4000)
402 common/popcor/pdb,ajsdef
405 COMMON /harlun/ qlun,iharlu
406 COMMON /nncms/ gamcm,bgcm,umo,pcm,eproj,pproj
407 COMMON /jspart/pxp(1000),
pyp(1000),pzp(1000),hep(1000),nnnp
421 DATA nbamlu /2,1,3,4,5,6,-2,-1,-3,-4/
425 IF(ipromu.NE.0)ipromm=ipromu
463 ELSEIF(aeo.GE.30.d0)
THEN
476 ainter=(aeo-15.d0)/15.d0
481 parj(42)=1.20d0-ainter*0.10d0
486 parj(41)=0.30d0+ainter*0.10d0
489 parj(21)=0.35d0+ainter*0.03d0
509 ELSEIF(ifrag.EQ.10)
THEN
514 IF(pproj.LE.30.d0)
THEN
516 ELSEIF(pproj.GE.100.d0)
THEN
519 dupar=(pproj-30.d0)/70.d0
520 parj(42)=0.6d0+dupar*5.4d0
531 ELSEIF(aeo.LT.7.d0)
THEN
533 ainter=-(aeo- 7.d0)/3.d0
535 IF(pproj.LE.30.d0)
THEN
537 ELSEIF(pproj.GE.100.d0)
THEN
540 dupar=(pproj-30.d0)/70.d0
541 dd42=0.25d0+dupar*5.4d0
543 parj(42)=0.35d0+ainter*dd42
553 ELSEIF(aeo.LT.10.d0)
THEN
570 ELSEIF(aeo.GE.30.d0)
THEN
585 ainter=-(aeo-30.d0)/20.d0
589 parj(42)=0.5d0-ainter*0.15d0
590 parj(42)=0.35d0-ainter*0.00d0
593 parj(21)=0.42d0+ainter*0.13d0
594 parj(21)=0.52d0+ainter*0.13d0
672 IF((kfa1.LE.0.OR.kfa1.GT.10).OR.
673 * (kfa2.LE.0.OR.kfa2.GT.10))
THEN
679 IF(iwarn.LE.20)
WRITE(6,*)
680 *
' BAMLUN KFA1,KFA2:',kfa11,kfa22,
681 *
' changed into :',kfa1,kfa2
690 IF(iwarn.LE.20)
WRITE(6,*)
691 *
' BAMLUN KFA1,KFA2 at energy AEO:',kfa11,kfa22,aeo,
692 *
' changed into:',kfa1,kfa2
693 ELSEIF((kfa1.EQ.10).AND.
699 IF(iwarn.LE.20)
WRITE(6,*)
700 *
' BAMLUN KFA1,KFA2 at energy AEO:',kfa11,kfa22,aeo,
701 *
' changed into:',kfa1,kfa2
711 IF(iwarn.LE.20)
WRITE(6,*)
712 *
' BAMLUN KFA1,KFA2 at energy AEO:',kfa11,kfa22,aeo,
713 *
' changed into:',kfa1,kfa2
720 IF(iwarn.LE.20)
WRITE(6,*)
721 *
' BAMLUN KFA1,KFA2 at energy AEO:',kfa11,kfa22,aeo,
722 *
' changed into:',kfa1,kfa2
729 IF(iwarn.LE.20)
WRITE(6,*)
730 *
' BAMLUN KFA1,KFA2 at energy AEO:',kfa11,kfa22,aeo,
731 *
' changed into:',kfa1,kfa2
738 IF(iwarn.LE.20)
WRITE(6,*)
739 *
' BAMLUN KFA1,KFA2 at energy AEO:',kfa11,kfa22,aeo,
740 *
' changed into:',kfa1,kfa2
743 IF(((kfa1.EQ.3.OR.kfa1.EQ.9).AND.
744 * (kfa2.EQ.3.OR.kfa2.EQ.9)).AND.aeo.LE.1.5d0)
THEN
749 IF(iwarn.LE.20)
WRITE(6,*)
750 *
' BAMLUN KFA1,KFA2 at energy AEO:',kfa11,kfa22,aeo,
751 *
' changed into:',kfa1,kfa2
753 IF(((kfa1.EQ.3.OR.kfa1.EQ.9).OR.
754 * (kfa2.EQ.3.OR.kfa2.EQ.9)).AND.aeo.LE.1.0d0)
THEN
760 IF(iwarn.LE.20)
WRITE(6,*)
761 *
' BAMLUN KFA1,KFA2 at energy AEO:',kfa11,kfa22,aeo,
762 *
' changed into:',kfa1,kfa2
767 IF(nbaml.LT.20)
WRITE(6,*)
' REJ. IN BAMLUN q-aq A < 0.8',aeo
779 ELSEIF(iopt.EQ.4.OR.iopt.EQ.6)
THEN
783 IF(nbaml.LT.20)
WRITE(6,*)
' REJ. IN BAMLUN q-qq E< 1.5 ',aeo
786 IF((kfa1.LE.0.OR.kfa1.GT.10).OR.
787 * (kfa2.LE.0.OR.kfa2.GT.10).OR.
788 * (kfa2.LE.0.OR.kfa3.GT.10))
THEN
796 IF(iwarn.LE.20)
WRITE(6,*)
797 *
' BAMLUN IOPT KFA1,KFA2,KFA3:',
799 * kfa33,
' changed into :',kfa1,kfa2,kfa3
809 IF(iwarn.LE.20)
WRITE(6,*)
810 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
811 *
' changed into:',kfa1,kfa2,kfa3
820 IF(iwarn.LE.20)
WRITE(6,*)
821 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
822 *
' changed into:',kfa1,kfa2,kfa3
831 IF(iwarn.LE.20)
WRITE(6,*)
832 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
833 *
' changed into:',kfa1,kfa2,kfa3
841 IF(iwarn.LE.20)
WRITE(6,*)
842 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
843 *
' changed into:',kfa1,kfa2,kfa3
852 IF(iwarn.LE.20)
WRITE(6,*)
853 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
854 *
' changed into:',kfa1,kfa2,kfa3
863 IF(iwarn.LE.20)
WRITE(6,*)
864 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
865 *
' changed into:',kfa1,kfa2,kfa3
876 IF(iwarn.LE.20)
WRITE(6,*)
877 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
878 *
' changed into:',kfa1,kfa2,kfa3
887 IF(iwarn.LE.20)
WRITE(6,*)
888 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
889 *
' changed into:',kfa1,kfa2,kfa3
898 IF(iwarn.LE.20)
WRITE(6,*)
899 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
900 *
' changed into:',kfa1,kfa2,kfa3
908 IF(iwarn.LE.20)
WRITE(6,*)
909 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
910 *
' changed into:',kfa1,kfa2,kfa3
919 IF(iwarn.LE.20)
WRITE(6,*)
920 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
921 *
' changed into:',kfa1,kfa2,kfa3
930 IF(iwarn.LE.20)
WRITE(6,*)
931 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
932 *
' changed into:',kfa1,kfa2,kfa3
936 IF(aeo.LT.1.25d0)
THEN
944 IF(iwarn.LE.20)
WRITE(6,*)
945 *
' BAMLUN KFA1,KFA2,KFA3 at energy AEO:',kfa11,kfa22,kfa33,aeo,
946 *
' changed into:',kfa1,kfa2,kfa3
954 IF(abs(ifl3).GT.abs(ifl2))
THEN
959 ifla2=1000*abs(ifl2)+100*abs(ifl3)
961 IF(abs(ifl3).LT.abs(ifl2).AND.
rndm(u).LE.0.25d0)iflz=1
963 IF(ifl2.LT.0)ifla2=-ifla2
972 ELSEIF(iopt.EQ.5)
THEN
975 IF(abs(ifl2).GT.abs(ifl1))
THEN
982 IF(abs(ifl4).GT.abs(ifl3))
THEN
987 ifla1=1000*abs(ifl1)+100*abs(ifl2)
988 ifla2=1000*abs(ifl3)+100*abs(ifl4)
990 IF(abs(ifl2).LT.abs(ifl1).AND.
rndm(u).LE.0.25d0)iflz=1
992 IF(ifl1.LT.0)ifla1=-ifla1
994 IF(abs(ifl4).LT.abs(ifl3).AND.
rndm(u).LE.0.25d0)iflz=1
996 IF(ifl3.LT.0)ifla2=-ifla2
1004 ELSEIF(iopt.EQ.10)
THEN
1009 IF(ipri.GE.1)
WRITE(6,*)ifla1,ifla2,aeo,ifla3
1011 103
FORMAT(
' BAMLUN',2i10,f10.2,i10)
1041 CALL
xseapa(aeo,xxxx,isq,isaq,xsq,xsaq,irej)
1042 WRITE(6,*)isq,isaq,xsq,xsaq,xx1
1044 IF(xsq.GE.xx1/2.d0.OR.xsaq.GE.xx3/2.d0)go to 1234
1053 IF(abs(iflasq).GT.abs(ifla2))
THEN
1058 iflad=1000*abs(ifla2)+100*abs(iflasq)
1060 IF(abs(iflasq).LT.abs(ifla2).AND.
rndm(u).LE.0.25d0)iflz=1
1062 WRITE(6,
'(4I10)')ifla1,iflad,iflasa,ifla3
1063 CALL
py4ent(1,ifla1,iflad,iflasa,ifla3,aeo,
1064 * xx1,xx2+xsq,xx3,0.d0,4.d0*xx1*xx3)
1074 ELSEIF(ifla1.LT.0)
THEN
1076 IF(abs(iflasa).GT.abs(ifla2))
THEN
1081 iflad=1000*abs(ifla2)+100*abs(iflasa)
1083 IF(abs(iflasa).LT.abs(ifla2).AND.
rndm(u).LE.0.25d0)iflz=1
1085 IF(ifla2.LT.0)iflad=-iflad
1086 WRITE(6,
'(4I10)')ifla1,iflad,iflasq,ifla3
1087 CALL
py4ent(1,ifla1,iflad,iflasq,ifla3,aeo,
1088 * xx1,xx2+xsaq,xx3,0.d0,xx1*xx3)
1099 CALL
py2ent(0,ifla1,ifla2,aeo)
1109 ELSEIF(iharlu.EQ.1)
THEN
1110 CALL
py2ent(-1,ifla1,ifla2,aeo)
1115 IF(nnnp.LT.1000)
THEN
1123 IF(ipri.GE.2)CALL
pylist(1)
1130 IF(k(iiii,1).EQ.4)k(iiii,1)=5
1134 IF(ipri.GE.2)
WRITE(6,*)
' After PYEXEC'
1135 IF(ipri.GE.2)CALL
pylist(1)
1147 IF(ipri.GE.2)
WRITE(6,*)
' After PYEDIT'
1149 IF(ipri.GE.2)CALL
pylist(1)
1150 IF(iharlu.EQ.1.AND.ndone.EQ.-107801)
THEN
1151 WRITE(6,*)
'NDONE ',ndone
1163 IF(ipri.GE.2)
WRITE(6,*)
' DPMJET COMMON particles'
1169 IF((k(i,1).EQ.11).AND.(k(i,2).EQ.92))
THEN
1177 IF(korjjj.EQ.0)korjjj=i
1180 IF((k(i,1).EQ.11).AND.(k(i,2).EQ.91))
THEN
1184 IF(korjjj.EQ.0)korjjj=i
1189 IF((k(i,1).EQ.1).OR.(k(i,1).EQ.4).OR.(k(i,1).EQ.11)
1190 * .OR.(k(i,1).EQ.15))
THEN
1192 IF((k(i,2).EQ.91).OR.(k(i,2).EQ.92).OR.(k(i,2).EQ.94))
THEN
1197 IF(abs(k(i,2)).LE.6) go to 101
1198 IF(abs(k(i,2)).GE.1000)
THEN
1203 IF(kll.EQ.1.OR.kll.EQ.3)go to 101
1206 IF(ihad.GT.nfimax)
THEN
1208 1112
FORMAT(.GT.
' BAMLUN: IHADNFIMAX INCREASE NFIMAX IHAD=',i10)
1211 IF(k(i,3).EQ.koriii(i))
THEN
1216 IF(korjjj.EQ.koriii(k(i,3)))
THEN
1219 iormo(ihad)=k(i,3)-korjjj
1223 iormo(ihad)=k(i,3)-korjjj-1-korkkk
1232 IF(ipri.GE.3)
WRITE(6,
'(4I10)')i,k(i,1),k(i,2),nrefb
1233 IF(nrefb.LT.1.OR.nrefb.GT.183)nrefb=4
1235 anf(ihad)=aname(nrefb)
1239 ichf(ihad)=ich(nrefb)
1240 ibarf(ihad)=ibar(nrefb)
1242 IF(k(i,1).EQ.11.OR.k(i,1).EQ.15)
THEN
1255 IF(ndone.EQ.-107801.AND.iharlu.EQ.1)
THEN
1256 WRITE(6,*)i,icoun,ihad,k(i,2),nrefb,k(i,3),
1257 * istath(ihad),ibarf(ihad),iormo(ihad),anf(ihad),
1258 * pxf(ihad),pyf(ihad),
1259 * pzf(ihad),hef(ihad),amf(ihad),koriii(i),
1262 102
FORMAT(
' BAMLUN',2i5,5i10,a8,5e12.3)
1268 IF(jni.NE.7)go to 777
1273 IF(k(i,1).EQ.11.AND.k(i,2).EQ.92)
THEN
1277 IF(k(i,1).EQ.11.OR.k(i,1).EQ.15)
THEN
1278 IF(k(i,2).EQ.92) go to 105
1279 IF(k(i,2).EQ.94) go to 105
1280 IF(abs(k(i,2)).LE.6) go to 105
1282 IF(ihad.GT.nfimax)
THEN
1286 IF(k(i,3).EQ.korii)
THEN
1290 iormo(ihad)=k(i,3)-korii
1298 IF(ipri.GE.3)
WRITE(6,
'(5I10)')i,k(i,1),k(i,2),k(i,3),nrefb
1299 IF(nrefb.LT.1.OR.nrefb.GT.183)nrefb=4
1301 anf(ihad)=aname(nrefb)
1304 ichf(ihad)=ich(nrefb)
1307 IF(ipri.GE.3)
WRITE(6,102)icoun,ihad,k(i,2),nrefb,k(i,3),
1308 * iormo(ihad),anf(ihad),
1309 * pxf(ihad),pyf(ihad),
1310 * pzf(ihad),hef(ihad),amf(ihad)
1320 IF(ifla1.EQ.-1.AND.ifla2.EQ.-3303)ipri=0
1323 SUBROUTINE xseapa(ECM,XXXX,IPSQ1,IPSAQ1,XPSQ1,XPSAQ1,IREJ)
1324 IMPLICIT DOUBLE PRECISION(
a-h,o-
z)
1326 COMMON /seasu3/seasq
1328 IF(ecm.LT.20.d0)
THEN
1329 xpthro=1.5d0*log10(ecm/200.d0)+3.5d0
1330 ELSEIF(ecm.GE.20.d0)
THEN
1333 xpthr=1.5d0*xpthro/(ecm**1.5d0*14.d0)
1338 xpthrx=xpthr-0.5d0*ai/ecm**2
1339 IF (xpthrx.LT.4.d0/ecm**2)xpthrx=4.d0/ecm**2
1347 ipsq1=1.+
rndm(ai)*(2.d0+seasqq)
1368 xglu1=
sampex(xpthrx,sox1)
1370 xpsq1=(0.2d0+(0.36d0*
rndm(ai))**0.50d0)*xglu1
1372 ELSEIF(ipsq1.EQ.3)
THEN
1374 IF (xpsq1.LE.0.3d0/ecm)go to 500
1375 IF (xpsaq1.LE.0.3d0/ecm)go to 500
1377 IF(xpsaq1.GE.xxxx)go to 500
1406 SUBROUTINE shmaki(NA,NCA,NB,NCB,RPROJ,RTARG,PPN)
1407 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1411 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
1412 * bsite(0:1,200),nstatb,nsiteb
1413 COMMON /dshms/ sigshs
1415 COMMON /rtar/ rtarnu
1417 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
1419 COMMON /zentra/ icentr
1420 COMMON /sigla/siglau
1421 COMMON /kglaub/jglaub
1425 rash=1.12d0*rna**0.33d0
1426 rbsh=1.12d0*rnb**0.33d0
1428 IF(na.EQ.9)rash=2.52d0
1429 IF(na.EQ.10)rash=2.45d0
1430 IF(na.EQ.11)rash=2.37d0
1431 IF(na.EQ.12)rash=2.45d0
1432 IF(na.EQ.13)rash=2.44d0
1433 IF(na.EQ.14)rash=2.55d0
1434 IF(na.EQ.15)rash=2.58d0
1435 IF(na.EQ.16)rash=2.71d0
1436 IF(na.EQ.17)rash=2.66d0
1437 IF(na.EQ.18)rash=2.71d0
1438 IF(nb.EQ.9)rbsh=2.52d0
1439 IF(nb.EQ.10)rbsh=2.45d0
1440 IF(nb.EQ.11)rbsh=2.37d0
1441 IF(nb.EQ.12)rbsh=2.45d0
1442 IF(nb.EQ.13)rbsh=2.44d0
1443 IF(nb.EQ.14)rbsh=2.55d0
1444 IF(nb.EQ.15)rbsh=2.58d0
1445 IF(nb.EQ.16)rbsh=2.71d0
1446 IF(nb.EQ.17)rbsh=2.66d0
1447 IF(nb.EQ.18)rbsh=2.71d0
1449 WRITE(6,*)
' SHMAKI: RASH, RBSH = ',rash,rbsh
1454 IF((icentr.EQ.1).AND.(na.GE.200).AND.(nb.GE.200))
THEN
1458 WRITE(6, 1010)nstatb,nsiteb
1460 1010
FORMAT(
' STATISTIC ON POINT ON B AND NUMBER OF POINTS',2i10)
1464 WRITE(6,*)
' SMAKI, PPN = ',ppn
1465 sigsh=
dshpto(ijproj,ppn)/10.d0
1470 ELSEIF(jglaub.EQ.1)
THEN
1471 sigshs=(
dshpto(ijproj,ppnn))/10.d0
1481 IF (ijproj.LE.12)bslope=8.5d0*(1.d0+0.065d0*
log(sss))
1482 IF (ijproj.GT.12)bslope=6.d0*(1.d0+0.065d0*
log(sss))
1483 gsh=1.d0/(2.d0*bslope*0.038938d0)
1485 IF (ijproj.LE.12)
THEN
1486 IF(ecm.GT.3.0d0.AND.ecm.LE.50.d0) rosh=-0.63d0+
1488 IF(ecm.GT.50.) rosh=0.1d0
1490 IF (ijproj.GT.12) rosh=0.01d0
1491 WRITE(6, 1030)sigsh,rosh,gsh,bslope,ecm
1493 1030
FORMAT(
' PARAMETERS OF THE NN AMPLITUDE SIG,RO,G,BSLOPE,ECM ' /5
1495 CALL
title(na,nb,nca,ncb)
1496 WRITE(6,*)
' vor PREVIO(RA,RB,NSTB,BMAX,BSTEP,SIG,RO,G)',
1497 &rash,rbsh,nsiteb,bmax,bstep,sigsh,rosh,gsh
1498 CALL
previo(rash,rbsh,nsiteb,bmax,bstep,sigsh,rosh,gsh)
1499 WRITE(6,*)
' SHMAKI: RASH, RBSH = ',rash,rbsh
1500 CALL
profb(bstep,nstatb,na,rash,nb,rbsh,bsite,nsiteb)
1501 WRITE(37,
'(4I5,F15.2,E15.5,F15.2)')
1502 * na,nca,nb,ncb,sigsh,ppn,siglau
1507 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1516 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1527 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
1528 * bsite(0:1,200),nstatb,nsiteb
1529 COMMON /dshms/ sigshs
1531 COMMON /rtar/ rtarnu
1533 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
1535 COMMON /dtumat/ bsiten(200,24,50),bsitem(200,24,50),
1536 + rprojj(50),rtargg(50),bstepp(50),bmaxx(50),
1537 + ntaxx(50),nztaxx(50),nprxx(50),nzprxx(50)
1540 COMMON /damp/ ca,ci,ga
1542 COMMON /kkmatu/kkmato,kkmata,ipoo,ipoa,ipzoo,ipzoa
1543 *,ibproo,ibproa,ireado
1554 IF(matnum.GT.50)
THEN
1555 WRITE(6,
'(2A,I3/A)')
1556 &
' Too large number of materials requested for Glauber',
1557 &
' initialization in SHMAKF / MATNUM=',matnum,
1558 &
' execution stopped in SHMAKF'
1561 WRITE(6,
'(A,I3,A)')
' Read Glauber data for material no.',matnum,
1565 READ(47,
'(A10)',
end=79) bnuc
1566 IF(bnuc.EQ.
' NUCLEUS ')
THEN
1568 READ(47,
'(A10,4I10)') bnuc,nprx,nzprx,ntax,nztax
1569 IF(nb.EQ.ntax.AND.ncb.EQ.nztax.AND.na.EQ.nprx.
1570 & and.nca.EQ.nzprx)
THEN
1572 nztaxx(matnum)=nztax
1574 nzprxx(matnum)=nzprx
1575 READ(47,
'(4F10.5)') bmaxx(matnum),bstepp(matnum),
1576 & rprojj(matnum),rtargg(matnum)
1578 READ(47,
'(5E16.8)') (bsiten(ida,ie,matnum),ida=1,200)
1582 READ(47,
'(5E16.8)') (bsitem(ida,ie,matnum),ida=1,200)
1590 WRITE(6,
'(A)')
' GLAUBER DATA NOT FOUND'
1595 help(i)=i*bstepp(matnum)
1596 WRITE (6,1040) help(i),(bsiten(i,ie,matnum),ie=1,24)
1597 1040
FORMAT (f10.4,10(1pe12.4)/10(1pe12.4))
1606 SUBROUTINE shmako(NA,NB,B,INTT,INTA,INTB,JS,JT,PPN,KKMAT)
1607 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1611 dimension js(namx),jt(namx)
1613 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
1618 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
1619 * bsite(0:1,200),nstatb,nsiteb
1620 COMMON /dshms/ sigshs
1622 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
1629 COMMON /dtumat/ bsiten(200,24,50),bsitem(200,24,50),
1630 + rprojj(50),rtargg(50),bstepp(50),bmaxx(50),
1631 + ntaxx(50),nztaxx(50),nprxx(50),nzprxx(50)
1633 COMMON /rtar/ rtarnu
1636 COMMON /damp/ ca,ci,ga
1638 COMMON /rptshm/ rproj,rtarg,bimpac
1640 COMMON /kkmatu/kkmato,kkmata,ipoo,ipoa,ipzoo,ipzoa
1641 *,ibproo,ibproa,ireado
1642 COMMON /kglaub/jglaub
1665 ELSEIF(jglaub.EQ.1)
THEN
1666 sigshs=(
dshpto(ijproj,ppn))/10.d0
1675 IF (ijproj.LE.12)bslope=8.5d0*(1.d0+0.065d0*
log(sss))
1676 IF (ijproj.GT.12)bslope=6.d0*(1.d0+0.065d0*
log(sss))
1678 gsh=1.d0/(2.d0*bslope*0.038938d0)
1680 IF (ibproj.LE.12)
THEN
1681 IF(ecm.GT.3.0d0.AND.ecm.LE.50.d0) rosh=-0.63d0+
1683 IF(ecm.GT.50.d0) rosh=0.1
1685 IF (ijproj.GT.12) rosh=0.01
1686 rca=gsh*sigshs/6.2831854d0
1687 fca=-gsh*sigshs*rosh/6.2831854d0
1693 ppo=wu10**(ipo+1) + 10.
1695 iread=iread + ipo - 1
1700 IF(iread.LE.0)iread=1
1701 IF(iread.GE.25)iread=24
1706 IF ((kkmat.NE.kkmato).OR.(ip.NE.ipoo).OR.
1707 * (ibproj.NE.ibproo).OR.(ipz.NE.ipzoo).OR.(iread.NE.ireado) )
THEN
1708 IF(ibproj.NE.0)
THEN
1710 bsite(1,ii)=bsiten(ii,iread,kkmat)
1718 bsite(1,ii)=bsitem(ii,iread,kkmat)
1728 WRITE(6,
'(A)')
' SHMAKO - BEFORE DIAGR'
1729 IF(ipri.GE.6)
WRITE(6, 1030)sigsh,rosh,gsh,bslope,ecm
1730 1030
FORMAT(
' PARAMETERS OF THE NN AMPLITUDE SIG,RO,G,BSLOPE,ECM ' /5
1733 CALL
diagr(na,nb,b,js,jt,intt,inta,intb)
1734 IF(ipev.GE.6)
WRITE(6,1000)na,nb,b,intt,inta,intb,js(1),jt(1),
1735 + pkoo(1,1),tkoo(1,1)
1736 1000
FORMAT(
' NA,NB,B,INTT,INTA,INTB,JS(1),JT(1),PKOO(1,1),TKOO(1,1)
1737 + IN SHMAKO '/2i5,f10.4,5i6,2f10.3)
1745 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1747 WRITE(6,1000)na,nca,nb,ncb
1748 1000
FORMAT(//10x,39hglauber
s approach is used to simulate ,
1749 +26hnucleus-nucleus collisions/ 24x,
1750 +40hthe calculation nas been carried out for/
1751 +27h projected nuclei with
a = ,i5,12h charge
a = ,i5/
1752 +24h
TARGET nuclei with b = ,i5,12h charge b = ,i5///)
1760 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1763 COMMON /kglaub/jglaub
1764 dimension x(3,
n),wd(4),rd(3)
1766 DATA sqr2/1.414216d0/
1767 DATA pdif/0.545d0/,r2min/0.16d0/
1768 DATA wd/0.d0,0.178d0,0.465d0,1.d0/
1769 DATA rd/2.09d0,0.935d0,0.697d0/
1776 ELSEIF (
n.EQ.4)
THEN
1778 ELSEIF (
n.GE.12)
THEN
1791 IF ((eps.GE.wd(i)).AND.(eps.LE.wd(i+1))) goto 40
1825 rad=rmax*(
rndm(v))**0.3333333
1827 fi=6.283185d0*
rndm(v)
1829 x(1,i)=rad*st*
cos(fi)
1830 x(2,i)=rad*st*
sin(fi)
1832 rr=
sqrt(x(1,i)**2+x(2,i)**2+x(3,i)**2)
1833 IF(jglaub.EQ.2)
f=1./(1.+
exp((rr-
r)/pdif))
1835 IF(
n.GE.11.OR.
n.LE.18)
THEN
1836 rr0=
r*
r/(2.5d0-4.d0/aan)
1837 f=(1.d0+(aan-4.d0)*rr**2/(6.d0*rr0))*6.d0/(aan-4.d0)
1838 * *
exp(-rr**2/rr0+(aan-10.d0)/(aan-4.d0))
1843 ELSEIF(
n.GE.9.AND.
n.LE.10)
THEN
1844 rr0=
r*
r/(2.5d0-4.d0/aan)
1845 f=(1.d0+(aan-4.d0)*rr**2/(6.d0*rr0))
1848 f=1./(1.+
exp((rr-
r)/pdif))
1851 IF (
rndm(v).GT.
f) goto 120
1853 IF (i.LT.2) goto 140
1857 dist2=(x(1,i)-x(1,i2))**2+(x(2,i)-x(2,i2))**2+(x(3,i)-x(3,i2))
1859 IF (dist2.LE.r2min) goto 120
1870 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1873 COMMON /zentra/ icentr
1874 COMMON /dshm/ rash,rbsh,bmax,bstap,sigsh,rosh,gsh,
1875 * bsita(0:1,200),nstatb,nsiteb
1876 COMMON /dshms/ sigshs
1878 dimension bsite(0:1,
n)
1886 IF (icentr.EQ.1)
THEN
1887 IF(rash.EQ.rbsh)
THEN
1888 IF(rash.LE.15.d0)
THEN
1892 ELSEIF(rash.GT.5.d0)
THEN
1897 ELSEIF(rash.LT.rbsh)
THEN
1898 bb=
rndm(v)*(rbsh-rash+3.d0)**2
1901 ELSEIF(rash.GT.rbsh)
THEN
1902 bb=
rndm(v)*(rash-rbsh+3.d0)**2
1913 left=((bsite(1,i0)-
y)*(bsite(1,i1)-
y)).LT.0.d0
1938 b=x0*(
y-y1)*(
y-y2)/((y0-y1)*(y0-y2)+1.
e-9)+ x1*(
y-y0)*(
y-y2)
1939 +/((y1-y0)*(y1-y2)+1.
e-9)+ x2*(
y-y0)*(
y-y1)/((y2-y0)*(y2-y1)+1.
e-9)
1942 IF (b.LT.0.0d0) b = x1
1943 IF (b.GT.bmax) b = bmax
1952 SUBROUTINE diagr(NA,NB,B,JS,JT,INT,INTA,INTB)
1953 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1967 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
1969 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
1970 * bsite(0:1,200),nstatb,nsiteb
1971 COMMON /dshms/ sigshs
1977 DOUBLE COMPLEX ca,ci
1978 COMMON /damp/ ca,ci,ga
1980 COMMON /kkmatu/kkmato,kkmata,ipoo,ipoa,ipzoo,ipzoa
1981 *,ibproo,ibproa,ireado
1982 COMMON /fluctu/ifluct
1983 COMMON /fluarr/flusi(1000),fluix(1000),fluixx(1000)
1985 dimension js(namx),jt(namx)
1992 CALL
modb(bsite,nsiteb,bstep,b)
2004 IF(intco.GE.500)
THEN
2006 CALL
modb(bsite,nsiteb,bstep,b)
2011 CALL
conucl(pkoo,na,rash)
2013 CALL
conucl(tkoo,nb,rbsh)
2030 IF(ipev.GE.6)
WRITE (6,1000)icnt,pkoo(1,1),tkoo(1,1)
2031 1000
FORMAT (
' 111 FORM IN DIAGR ICNT,PKOO(1,1),TKOO(1,1) ',i6,2f10.3)
2038 ELSEIF(ifluct.EQ.1)
THEN
2039 ifuk=(
rndm(v)+0.001)*1000.
2046 IF(
xy.GT.15.d0) go to 40
2052 IF(
rndm(v).LT.
p) go to 40
2061 IF (js(i).NE.0) inta=inta+1
2064 IF (jt(j).NE.0) intb=intb+1
2068 +
' DIAGR - AFTER 30 CONTINUE: ICNT, INT, B, NA,RA, NB,RB'
2069 WRITE(6,
'(I10,I5,1PE11.3,2(I5,1PE11.3))') icnt,
int, b, na,rash,
2072 IF(
int.EQ.0) go to 30
2077 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2081 COMMON /fluarr/ flusi(1000),fluix(1000),fluixx(1000)
2092 flus=((x-b)/(om*b))**
n
2094 flusi(i)=(x/b)*
exp(-((x-b)/(om*b))**
n)/(x/b+
a)
2098 flusu=flusu+flusi(i)
2101 flusuu=flusuu+flusi(i)/flusu
2105 3
FORMAT(
' FLUCTUATIONS')
2106 CALL
plot(fluix,flusi,1000,1,1000,0.d0,0.06d0,0.d0,0.01d0)
2110 IF(af.LE.flusi(j))
THEN
2118 fluixx(1000)=fluix(1000)
2138 SUBROUTINE calbam(NNCH,I1,I2,IFB11,IFB22,IFB33,IFB44,
2140 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2186 parameter(nfimax=249)
2187 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
2188 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
2189 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
2192 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
2202 IF(nobam.EQ.4.AND.isymm.EQ.1)
THEN
2204 IF (amch.LT.3.d0)
THEN
2206 IF (rr.LT.0.33333d0)
THEN
2210 ELSEIF (rr.GT.0.666666d0)
THEN
2219 ELSEIF(amch.GT.7.d0)
THEN
2225 ssss=(7.d0-amch)/4.d0
2227 IF(rrr.LT.1.d0-ssss)
THEN
2233 IF (rr.LT.0.33333d0)
THEN
2237 ELSEIF (rr.GT.0.666666d0)
THEN
2248 ELSEIF(nobam.EQ.6.AND.isymm.EQ.1)
THEN
2250 IF (amch.LT.3.d0)
THEN
2252 IF (rr.LT.0.33333d0)
THEN
2256 ELSEIF (rr.GT.0.666666d0)
THEN
2265 ELSEIF(amch.GT.7.d0)
THEN
2271 ssss=(7.d0-amch)/4.d0
2273 IF(rrr.LT.1.d0-ssss)
THEN
2279 IF (rr.LT.0.33333d0)
THEN
2283 ELSEIF (rr.GT.0.666666d0)
THEN
2302 WRITE (6,1000)nnch,i1,i2,ifb1,ifb2,ifb3,ifb4,amch,nobam,ihad
2303 1000
FORMAT(
' CALBAM:NNCH,I1,I2,IFB1,IFB2,IFB3,IFB4,AMCH,NOBAM,IHAD' /7
2307 CALL
dbamje(ihad,ifb1,ifb2,ifb3,ifb4,amch,nobam)
2312 IF(nnch.EQ.-99)
THEN
2327 IF (ifb1.LE.6) go to 40
2329 IF(nobam.EQ.3) CALL
dbamje(ihad,ifb2,ifb1,ifb3,ifb4,amch,nobam)
2332 CALL
dbamje(ihad,ifb1,ifb2,ifb3,ifb4,amch,nobam)
2334 IF(nobam.EQ.6) CALL
dbamje(ihad,ifb3,ifb1,ifb2,ifb4,amch,4)
2336 CALL
dbamje(ihad,ifb3,ifb4,ifb1,ifb2,amch,nobam)
2341 IF (nobam.EQ.3.OR.nobam.EQ.4.OR.nobam.EQ.5)
THEN
2342 CALL
dbamje(ihad,ifb1,ifb2,ifb3,ifb4,amch,nobam)
2343 ELSEIF(nobam.EQ.6)
THEN
2345 CALL
dbamje(ihad,ifb3,ifb1,ifb2,ifb4,amch,4)
2354 IF (iturn.EQ.0) go to 100
2363 WRITE(6,1245)i,pzf(i),hef(i),anf(i)
2364 1245
FORMAT(i5,2f10.4,a8)
2374 SUBROUTINE dbamje(IHAD,KFA1,KFA2,KFA3,KFA4,AE0,IOPT)
2375 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2395 parameter(nfimax=249)
2396 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
2397 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
2398 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
2401 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
2402 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
2405 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210),ich(210), ibar
2406 +(210),k1(210),k2(210)
2407 common/dremai/ rpxr,rpyr,rpzr,rer,kr1r,kr2r
2410 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
2412 COMMON /diffra/ isingd,idiftp,ioudif,iflagd
2414 dimension rpx(100),rpy(100),re(100)
2415 dimension kfr1(100),kfr2(100),iv(100)
2416 parameter(pimass=0.15d0)
2420 IF (lt.EQ.1)
WRITE(6, 1000)ihad,kfa1,kfa2,kfa3,kfa4,ae0,iopt
2421 1000
FORMAT (5i5,e12.4,i5,
2422 +
' BAMJET,IHAD,KFA1,KFA2,KFA3,KFA4,AE0,IOPT')
2429 CALL
bamlun(ihad,kfa1,kfa2,kfa3,kfa4,ae0,iopt,irej)
2435 IF(ifrag.EQ.1.OR.ifrag.EQ.2.OR.ifrag.GE.10)
THEN
2437 CALL
bamlun(ihad,kfa1,kfa2,kfa3,kfa4,ae0,iopt,irej)
2475 IF(itry.GT.10000)
THEN
2476 WRITE(6,
'(/1X,A)')
'DBAMJE:ERROR: FRAGMENTATION IMPOSSIBLE'
2477 WRITE(6,
'(1X,A,5I5,E12.3,I5)')
2478 &
'DBAMJE:IHAD,KFA1,KFA2,KFA3,KFA4,AE0,IOPT ',
2479 & ihad,kfa1,kfa2,kfa3,kfa4,ae0,iopt
2493 IF(itry.GT.100)
THEN
2495 IF(kfa1.GT.(kfa2-6))
THEN
2496 e0=ae0-max(ae0*0.1,pimass)
2498 e0=max(ae0*0.1,pimass)
2503 IF(iopt.EQ.1.OR.iopt.EQ.2) e0=ae0
2506 IF(kfa1.GT.6.AND.iopt.EQ.1) ll=1
2507 IF(kfa1.LE.6.AND.iopt.EQ.2) ll=1
2508 IF(kfa1.GT.6.AND.iopt.EQ.4) ll=1
2524 IF(iopt.EQ.3.AND.ll.EQ.0) goto 50
2525 IF(iopt.EQ.4.AND.kfa1.LE.6.AND.ll.EQ.0) goto 50
2526 IF(iopt.EQ.4.AND.kfa1.GT.6.AND.ll.EQ.1) goto 50
2527 IF(iopt.EQ.5.AND.ll.EQ.0) goto 50
2537 IF(iopt.EQ.4.AND.kfa1.GT.6) ll=0
2551 IF(amf(it).GT.rx) goto 10
2552 IF(amf(it).LE.rx) goto 70
2568 90 rpx(it)=rpx(j)-hpx
2571 IF (iopt.EQ.1.AND.ll.EQ.1)hpz=-hpz
2572 IF(iopt.EQ.2.AND.ll.EQ.1) hpz=-hpz
2573 IF(iopt.EQ.4.AND.kfa1.GT.6) hpz=-hpz
2574 IF(iopt.EQ.5) hpz=-hpz
2582 WRITE(6, 1010)pgx,pgy,pgz
2583 1010
FORMAT(1h0,12hpgx,pgy,pgz=,3f8.4)
2588 IF(iopt.EQ.1.OR.iopt.EQ.2) goto 150
2596 IF(iopt.EQ.3.OR.iopt.EQ.4.OR.iopt.EQ.5) goto 160
2597 IF(ll.EQ.0) goto 160
2605 IF(le.EQ.0) goto 180
2608 WRITE(6, 1020)nref(j),anf(j),amf(j),ichf(j),
2609 + ibarf(j),pxf(j),pyf(j),
2611 1020
FORMAT(2x,i3,a6,f6.3,2i4,4f8.4)
2612 1030
FORMAT(2x,
'NF,NAME,MASS,IQ,IB,PX,PY,PZ,E')
2615 1040
FORMAT(1h0,38hnumber of events with prest gt. erest=,i4, /,
2616 +21hnumber of all events=,i4)
2617 1050
FORMAT(1h0,
'NUMBER OF EVENTS WITH ONLY ONE PARTICLE=',i4)
2622 1060
FORMAT(1h0,
' MULTIPLIZITAET=',i3)
2623 1070
FORMAT(1h0,13hhadronanzahl=,i3)
2638 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2645 IF(kp.EQ.4.AND.ks.EQ.5)ind=4
2647 IF(kp.EQ.6.AND.ks.EQ.7)ind=6
2648 IF(kp.EQ.4.AND.ks.EQ.4)ind=7
2649 IF(kp.EQ.6.AND.ks.EQ.5)ind=8
2652 IF(kp.EQ.12.AND.ks.EQ.8)ind=11
2654 IF(kp.EQ.12.AND.ks.EQ.7)ind=13
2672 DOUBLE PRECISION FUNCTION dbeta(X1,X2,BET)
2673 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2677 IF(betx1.LT.70.) ax=-1./bet**2*(betx1+1.)*
exp(-betx1)
2678 ay=1./bet**2*(bet*x2+1.)*
exp(-bet*x2)
2687 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2689 x1=cops*x-sips*cote*
y+sips*site*
z
2690 x2=sips*x+cops*cote*
y-cops*site*
z
2748 SUBROUTINE dthrep(UMO,ECM1,ECM2,ECM3,PCM1,PCM2,PCM3,COD1,COF1,
2749 &sif1,cod2,cof2,sif2,cod3,cof3,sif3,am1,am2,am3)
2796 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2798 parameter( kalgnm = 2 )
2799 parameter( anglgb = 5.0
d-16 )
2800 parameter( anglsq = 2.5
d-31 )
2801 parameter( axcssv = 0.2
d+16 )
2802 parameter( andrfl = 1.0
d-38 )
2803 parameter( avrflw = 1.0
d+38 )
2804 parameter( ainfnt = 1.0
d+30 )
2805 parameter( azrzrz = 1.0
d-30 )
2806 parameter( einfnt = +69.07755278982137
d+00 )
2807 parameter( ezrzrz = -69.07755278982137
d+00 )
2808 parameter( onemns = 0.999999999999999
d+00 )
2809 parameter( onepls = 1.000000000000001
d+00 )
2810 parameter( csnnrm = 2.0
d-15 )
2811 parameter( dmxtrn = 1.0
d+08 )
2843 parameter( zerzer = 0.
d+00 )
2844 parameter( oneone = 1.
d+00 )
2845 parameter( twotwo = 2.
d+00 )
2846 parameter( thrthr = 3.
d+00 )
2847 parameter( foufou = 4.
d+00 )
2848 parameter( fivfiv = 5.
d+00 )
2849 parameter( sixsix = 6.
d+00 )
2850 parameter( sevsev = 7.
d+00 )
2851 parameter( eigeig = 8.
d+00 )
2852 parameter( aninen = 9.
d+00 )
2853 parameter( tenten = 10.
d+00 )
2854 parameter( hlfhlf = 0.5
d+00 )
2855 parameter( onethi = oneone / thrthr )
2856 parameter( twothi = twotwo / thrthr )
2857 parameter( pipipi = 3.1415926535897932270
d+00 )
2858 parameter( eneper = 2.7182818284590452354
d+00 )
2859 parameter( sqrent = 1.6487212707001281468
d+00 )
2901 parameter( clight = 2.99792458
d+10 )
2902 parameter( avogad = 6.0221367
d+23 )
2903 parameter( amelgr = 9.1093897
d-28 )
2904 parameter( plckbr = 1.05457266
d-27 )
2905 parameter( elccgs = 4.8032068
d-10 )
2906 parameter( elcmks = 1.60217733
d-19 )
2907 parameter( amugrm = 1.6605402
d-24 )
2908 parameter( ammumu = 0.113428913
d+00 )
2923 parameter( alpfsc = 7.2973530791728595
d-03 )
2924 parameter( fscto2 = 5.3251361962113614
d-05 )
2925 parameter( fscto3 = 3.8859399018437826
d-07 )
2926 parameter( fscto4 = 2.8357075508200407
d-09 )
2927 parameter( plabrc = 0.197327053
d+00 )
2928 parameter( amelct = 0.51099906
d-03 )
2929 parameter( amugev = 0.93149432
d+00 )
2930 parameter( ammuon = 0.105658389
d+00 )
2931 parameter( rclsel = 2.8179409183694872
d-13 )
2932 parameter( gevmev = 1.0
d+03 )
2933 parameter( emvgev = 1.0
d-03 )
2934 parameter( raddeg = 180.
d+00 / pipipi )
2935 parameter( degrad = pipipi / 180.
d+00 )
2962 parameter( lunin = 5 )
2963 parameter( lunout = 6 )
2964 parameter( lunerr = 66 )
2965 parameter( lunber = 14 )
2966 parameter( lunech = 8 )
2967 parameter( lunflu = 86 )
2968 parameter( lungeo = 16 )
2969 parameter( lunpgs = 12 )
2970 parameter( lunran = 2 )
2971 parameter( lunxsc = 81 )
2972 parameter( lunrdb = 1 )
3001 parameter( mxxrgn = 500 )
3002 parameter( mxxmdf = 56 )
3003 parameter( mxxmde = 50 )
3004 parameter( mfstck = 1000 )
3005 parameter( mestck = 100 )
3006 parameter( nallwp = 39 )
3007 parameter( mpdpdx = 8 )
3008 parameter( icomax = 180 )
3009 parameter( nstbis = 304 )
3010 parameter( idmaxp = 210 )
3016 dimension
f(5),
xx(5)
3018 COMMON /dgamre/ redu,amo,amm(15 )
3019 common/ddrei/uumo,aam1,aam2,aam3,s22,umo2,
3020 *am11,am22,am33,s2sup,s2sap(2)
3059 IF(rho2.LT.rho1) go to 125
3061 125 s2sup=(s22-s21)*.5d0+s21
3062 suprho=
dxlamb(s2sup,umo2,am11)*
dxlamb(s2sup,am22,am33)/
3064 suprho=suprho*1.05d0
3066 IF (gu.LT.go.AND.xo.LT.gu) xo=gu
3067 IF (gu.GT.go.AND.xo.GT.gu) xo=gu
3076 x4=(
xx(1)+
xx(2))*0.5d0
3077 x5=(
xx(2)+
xx(3))*0.5d0
3087 IF (
f(ii).GE.
f(iii)) go to 128
3100 IF (
xx(ii).GE.
xx(iii)) go to 129
3114 IF (ith.GT.200) redu=-9.d0
3115 IF (ith.GT.200) go to 400
3118 s2=am23+c*(umo-am1-am2-am3)*(umo-am1+am2+am3)
3122 IF(
y.GT.
rho) go to 1
3125 s1=s1*
rho+am11+am22-(s2-umo2+am11)*(s2+am22-am33)/(2.d0*s2)-
3127 s3=umo2+am11+am22+am33-s1-s2
3128 ecm1=(umo2+am11-s2)/umoo
3129 ecm2=(umo2+am22-s3)/umoo
3130 ecm3=(umo2+am33-s1)/umoo
3131 pcm1=
sqrt((ecm1+am1)*(ecm1-am1))
3132 pcm2=
sqrt((ecm2+am2)*(ecm2-am2))
3133 pcm3=
sqrt((ecm3+am3)*(ecm3-am3))
3138 IF ( pcm12 .LT. anglsq ) go to 200
3139 costh=(ecm1*ecm2+0.5
d+00*(am11+am22-s1))/pcm12
3143 costh=(uw-0.5
d+00)*2.
d+00
3147 IF(abs(costh).GT.oneone)
3148 &costh=sign(oneone,costh)
3149 IF (redu.LT.1.
d+00)
RETURN
3150 costh2=(pcm3*pcm3+pcm2*pcm2-pcm1*pcm1)/(2.
d+00*pcm2*pcm3)
3153 IF(abs(costh2).GT.oneone)
3154 &costh2=sign(oneone,costh2)
3155 sinth2=
sqrt((oneone-costh2)*(oneone+costh2))
3156 sinth =
sqrt((oneone-costh)*(oneone+costh))
3157 sinth1=costh2*sinth-costh*sinth2
3158 costh1=costh*costh2+sinth2*sinth
3170 cod3=twotwo*
rndm(cod3)-oneone
3171 sid3=
sqrt((1.
d+00-cod3)*(1.
d+00+cod3))
3173 cod1=cx11*cod3+cz11*sid3
3174 chlp=(oneone-cod1)*(oneone+cod1)
3175 IF(chlp.LT.1.
d-14)
WRITE(isys,2)cod1,cof3,sid3,
3178 cof1=(cx11*sid3*cof3-cy11*sif3-cz11*cod3*cof3)/sid1
3179 sif1=(cx11*sid3*sif3+cy11*cof3-cz11*cod3*sif3)/sid1
3180 cod2=cx22*cod3+cz22*sid3
3181 sid2=
sqrt((oneone-cod2)*(oneone+cod2))
3182 cof2=(cx22*sid3*cof3-cy22*sif3-cz22*cod3*cof3)/sid2
3183 sif2=(cx22*sid3*sif3+cy22*cof3-cz22*cod3*sif3)/sid2
3186 eochck = umo - ecm1 - ecm2 - ecm3
3190 pzchck = pcm1 * cod1 + pcm2 * cod2 + pcm3 * cod3
3191 pxchck = pcm1 * cof1 * sid1 + pcm2 * cof2 * sid2
3192 & + pcm3 * cof3 * sid3
3193 pychck = pcm1 * sif1 * sid1 + pcm2 * sif2 * sid2
3194 & + pcm3 * sif3 * sid3
3195 eocmpr = 1.
d-12 * umo
3196 IF ( abs(eochck) + abs(pxchck) + abs(pychck) + abs(pzchck)
3197 & .GT. eocmpr )
THEN
3199 &
' *** Threpd: energy/momentum conservation failure! ***',
3200 & eochck,pxchck,pychck,pzchck
3201 WRITE(lunerr,*)
' *** SID1,SID2,SID3',sid1,sid2,sid3
3208 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3211 COMMON /dgamre/ redu,amo,amm(15 )
3212 common/ddrei/
test(12)
3216 IF (idgb.LE.0) go to 11
3219 WRITE(6,10) xlam,x,
y,
z,
test
3221 12
FORMAT(/,10x,
' DXLAMB PRINT')
3222 13
FORMAT(10x,60(1h*))
3223 10
FORMAT(4e20.8,
'DXLAMB',/,12f10.5)
3244 SUBROUTINE strafo(GAM,BGAM,CX,CY,CZ,COD,COF,SIF,P,ECM,
3246 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3249 sid=
sqrt((1.-cod)*(1.+cod)+1.
e-22)
3250 sif=
sqrt((1.-cof)*(1.+cof)+1.
e-22)
3254 plz=gam*pcmz+bgam*ecm
3255 pl=
sqrt(plx*plx+ply*ply+plz*plz)
3256 el=gam*ecm+bgam*pcmz
3259 IF(coz.GE.1.)coz=0.999999999999
3260 siz=
sqrt((1.-coz)*(1.+coz))
3261 CALL
drtran(cx,cy,cz,coz,siz,sif,cof,cxl,cyl,czl)
3265 SUBROUTINE drtran(XO,YO,ZO,CDE,SDE,SFE,CFE,X,Y,Z)
3266 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3268 IF (abs(xo)-0.0001) 10,10,30
3269 10
IF (abs(yo)-0.0001) 20,20,30
3280 x=-yo*xi/
a-zo*xo*yi/
a+xo*zi
3281 y=xo*xi/
a-zo*yo*yi/
a+yo*zi
3290 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3295 parameter(idmax9=602)
3297 common/ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
3301 common/dpar/aname(210),am(210),ga(210),
tau(210),ich(210),ibar(210)
3326 SUBROUTINE dtwopd(UMO,ECM1,ECM2,PCM1,PCM2,COD1,COF1,SIF1, COD2,
3328 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3332 IF(umo.LT.(am1+am2))
THEN
3333 WRITE(6,
'(/,A/A,3(1PE12.4))')
3334 +
' INCONSISTENT CALL OF TWOPAD / EXECUTION STOPPED',
3335 +
' UMO, AM1, AM2 :', umo, am1, am2
3339 ecm1=((umo-am2)*(umo+am2) + am1*am1)/(2.*umo)
3341 pcm1=
sqrt((ecm1-am1)*(ecm1+am1))
3355 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3365 IF(
rndm(v).LT.0.5d0) goto20
3375 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3377 CHARACTER*8 zkname,
z
3380 parameter(idmax9=602)
3382 common/ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
3386 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210), ich(210),ibar
3387 +(210),k1(210),k2(210)
3389 dimension ichar(210)
3390 equivalence(ich(1),ichar(1))
3394 1000
FORMAT(1h1,
' ')
3396 1010
FORMAT(///
' TABLE OF USED PARTICLES AND RESONANCES (I)',//
3397 +
' I = NUMBER OF PARTICLE OR RESONANCE',/
3398 +
' IPDG = P D G NUMBER OF PARTICLE OR RESONANCE',/
3399 +
' ANAME = NAME OF I'/,
' AM = MASS OF I (GEV)',/
3400 +
' GA = WIDTH OF I (GEV)',/
' TAU = LIFE TIME OF I (SEC.)',/
3401 +
' ICH = ELECTRIC CHARGE OF I, IBAR = BARYONIC CHARGE OF I',/
' ',
'
3402 +K1 = FIRST DECAY CHANNEL NUMBER, K2 = LAST DECAY CHANNEL NUMBER OF
3408 +
' I ANAME AM GA TAU ICH IBAR K1 K2'/)
3412 WRITE(6, 1030)i,ipdg,aname(i),am(i),
3413 + ga(i),
tau(i),ich(i),ibar(i), k1
3415 1030
FORMAT (1i4,i6,2x,1a8,3e11.4,4i4)
3416 IF(i.EQ.43)
WRITE(6, 1000)
3417 IF(i.EQ.43)
WRITE(6, 1020)
3418 IF(i.EQ.99)
WRITE(6, 1000)
3419 IF(i.EQ.99)
WRITE(6, 1020)
3420 IF(i.EQ.155)
WRITE(6, 1000)
3421 IF(i.EQ.155)
WRITE(6, 1020)
3425 1040
FORMAT(///
' DECAY CHANNELS OF PARTICLES AND RESONANCES',//)
3427 1050
FORMAT(
' ANAME = PARTICLE AND RESONANCE NAME'/,
3428 +
' DNAME = DECAY CHANNEL NAME'/,
' J = DECAY CHANNEL NUMBER'/,
3429 +
' I = NUMBER OF DECAYING PARTICLE'/,
3430 +
' WT = SUM OF DECAY CHANNEL WEIGHTS FROM K1(I) UP TO J'/,
3431 +
' NZK = PROGRAM INTERNAL NUMBERS OF DECAY PRODUCTS')
3434 1060
FORMAT(///
' I J ANAME DNAME DECAY
3439 IF (ik1.LE.0) go to 60
3453 WRITE(6, 1070)i,ik,aname(i),zkname(ik),(
z(j),j=1,3),wt(ik),j1,j2,
3455 1070
FORMAT(2i5,
' DECAY OF ',1a8,
' (CHANNEL: ',1a6,
' ) TO ',3(1a6,2x),
3457 amtest=am(i)-am(j1)-am(j2)-am(j3)
3458 ibtest=ibar(i)-ibar(j1)-ibar(j2)-ibar(j3)
3459 ictest=ichar(i)-ichar(j1)-ichar(j2)-ichar(j3)
3460 IF (amtest) 20,30,30
3465 IF (mtest+ibtest**2+ictest**2.NE.0)
WRITE(6, 1080)amtest,
3468 1080
FORMAT (
' ***** ERROR IN MASS, BAR.CH. OR E.CH. ',f10.5,2i6)
3469 IF(ik.EQ.27)
WRITE(6, 1000)
3470 IF(ik.EQ.27)
WRITE(6, 1060)
3471 IF(ik.EQ.62)
WRITE(6, 1000)
3472 IF(ik.EQ.62)
WRITE(6, 1060)
3473 IF(ik.EQ.101)
WRITE(6, 1000)
3474 IF(ik.EQ.101)
WRITE(6, 1060)
3475 IF(ik.EQ.144)
WRITE(6, 1000)
3476 IF(ik.EQ.144)
WRITE(6, 1060)
3477 IF(ik.EQ.183)
WRITE(6, 1000)
3478 IF(ik.EQ.183)
WRITE(6, 1060)
3479 IF(ik.EQ.222)
WRITE(6, 1000)
3480 IF(ik.EQ.222)
WRITE(6, 1060)
3481 IF(ik.EQ.261)
WRITE(6, 1000)
3482 IF(ik.EQ.261)
WRITE(6, 1060)
3483 IF(ik.EQ.300)
WRITE(6, 1000)
3484 IF(ik.EQ.300)
WRITE(6, 1060)
3485 IF(ik.EQ.362)
WRITE(6, 1000)
3486 IF(ik.EQ.362)
WRITE(6, 1060)
3487 IF(ik.EQ.401)
WRITE(6, 1000)
3488 IF(ik.EQ.401)
WRITE(6, 1060)
3489 IF(ik.EQ.440)
WRITE(6, 1000)
3490 IF(ik.EQ.440)
WRITE(6, 1060)
3491 IF(ik.EQ.479)
WRITE(6, 1000)
3492 IF(ik.EQ.479)
WRITE(6, 1060)
3493 IF(ik.EQ.518)
WRITE(6, 1000)
3494 IF(ik.EQ.518)
WRITE(6, 1060)
3504 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3506 common/dinpda/imps(6,6),imve(6,6),ib08(6,21),ib10(6,21),
3507 *ia08(6,21),ia10(6,21),a1,b1,b2,b3,lt,lb,bet,as,b8,ame,diq,isu
3508 dimension iv(36),ip(36),ib(126),ibb(126),ia(126),iaa(126)
3515 *23,14,16,116,0,0,13,23,25,117,0,0,15,24,31,120,0,0,119,118,121,
3525 *33,34,38,123,0,0,32,33,39,124,0,0,36,37,96,127,0,0,126,125,128,
3535 *0,1,21,140,0,0,8,22,137,0,0,97,138,0,0,146,5*0,
3536 *1,8,22,137,0,0,0,20,142,0,0,98,139,0,0,147,5*0,
3537 *21,22,97,138,0,0,20,98,139,0,0,0,145,0,0,148,5*0,
3538 *140,137,138,146,0,0,142,139,147,0,0,145,148,50*0/
3547 *53,54,104,161,0,0,55,105,162,0,0,107,164,0,0,167,5*0,
3548 *54,55,105,162,0,0,56,106,163,0,0,108,165,0,0,168,5*0,
3549 *104,105,107,164,0,0,106,108,165,0,0,109,166,0,0,169,5*0,
3550 *161,162,164,167,0,0,163,165,168,0,0,166,169,0,0,170,47*0/
3559 *0,2,99,152,0,0,9,100,149,0,0,102,150,0,0,158,5*0,
3560 *2,9,100,149,0,0,0,101,154,0,0,103,151,0,0,159,5*0,
3561 *99,100,102,150,0,0,101,103,151,0,0,0,157,0,0,160,5*0,
3562 *152,149,150,158,0,0,154,151,159,0,0,157,160,50*0/
3571 *67,68,110,171,0,0,69,111,172,0,0,113,174,0,0,177,5*0,
3572 *68,69,111,172,0,0,70,112,173,0,0,114,175,0,0,178,5*0,
3573 *110,111,113,174,0,0,112,114,175,0,0,115,176,0,0,179,5*0,
3574 *171,172,174,177,0,0,173,175,178,0,0,176,179,0,0,180,47*0/
3636 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3656 parameter(nfimax=249)
3657 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
3658 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
3659 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
3662 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
3663 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
3665 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
3669 COMMON /dkpl/uplo,ipq
3670 COMMON /dinv/pnuc(3),inucvt
3671 common/capkop/xx1,xx3
3674 1000
FORMAT (
' ##############################################'/
3675 +
' PROGRAM TECABAPT'/
3676 +
' ######################################################')
3686 READ(5,1010)jni,nevt,ip,ncases,poo,aoo,znuc
3687 1010
FORMAT(4i10,3f10.2)
3689 IF (jni.LE.0) go to 120
3690 WRITE(6, 1010)jni,nevt,ip,ncases,poo,aoo,znuc,it
3694 go to(20,30,40,50,60,70,80,100),jni
3713 CALL
distcm(1,ipq,poo,ipq,ipq)
3716 IF (ip.EQ.103)CALL
calbam(0,1,1,1,7,1,1,poo,3,nhad)
3717 IF (ip.EQ.109)CALL
calbam(0,1,1,7,1,1,1,poo,3,nhad)
3718 IF (ip.EQ.104)CALL
calbam(0,1,1,1,2,2,1,poo,4,nhad)
3719 IF (ip.EQ.1010)CALL
calbam(0,1,1,1,2,3,1,poo,4,nhad)
3720 IF (ip.EQ.105)CALL
calbam(0,1,1,1,2,7,8,poo,5,nhad)
3721 IF (ip.EQ.1011)CALL
calbam(0,1,1,7,7,1,1,poo,5,nhad)
3722 IF (ip.EQ.106)CALL
calbam(0,1,1,1,1,1,1,poo,6,nhad)
3723 IF (ip.EQ.1012)CALL
calbam(0,1,1,7,7,7,1,poo,6,nhad)
3724 IF (ip.EQ.1050)CALL
calbam(0,1,1,1,1,2,1,poo,10,nhad)
3727 CALL
distcm(2,nhad,poo,ipq,ncases)
3730 WRITE(6, 1020)poo,ip,ncases
3731 1020
FORMAT (
' BAMJET (POO,IP,NCASES) = ',1f10.2,2i10)
3732 CALL
distcm(3,nevt,poo,ipq,ncases)
3745 SUBROUTINE distcm(IOP,NHAD,POLAB,KPROJ,KTARG)
3746 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3753 parameter(nfimax=249)
3754 COMMON /dfinpa/ anh(nfimax),
px(nfimax),
py(nfimax),
pz(nfimax),
3755 +he(nfimax),am(nfimax), ich(nfimax),ibar(nfimax),nr(nfimax)
3756 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
3760 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
3761 +iibar(210),k1(210),k2(210)
3763 COMMON /histo / xmult(100,10),ymult(100,10),xxfl(50,20), yxfl
3764 +(50,20),xyl(50,20),yyl(50,20), yylps(50,20),ptp(50,20),pty(50,20),
3766 dimension avmult(12,30),ave(12,30),indx(25),mu(12,30), akno
3767 +(100,2),xkno(100,2),ake(12,30),aaso(12,30)
3768 COMMON /dkpl/uplo,kpl
3772 DATA indx/1,8,10,10,10,10,7,2,7,10,10,7,3,4,5,6,
3773 *11,12,7,13,14,15,16,17,18/
3775 go to(10,60,100),iop
3784 IF(jni.EQ.7) umo=polab
3788 WRITE(6, 1000)eeo,po,nhad
3789 1000
FORMAT (
' EEO',f10.2,f10.2,i10)
3793 avmult(kpl,i)=1.
e-18
3801 WRITE(6, 1000)eeo,po,nhad
3804 xxfl(j,i)=j*dxfl -1.
3812 WRITE(6, 1000)eeo,po,nhad
3816 avmult(kpl,i)=1.
d-18
3820 WRITE(6, 1000)eeo,po,nhad
3825 avmult(kpl,30)=avmult(kpl,30)+nhad
3827 IF (nnhad.GT.100) nnhad=100
3828 ymult(nnhad,10)=ymult(nnhad,10)+1.
3834 IF(ibar(i).NE.500)
THEN
3839 IF(eetot.GT.polab+1.
d-6)
THEN
3840 WRITE(6,*).gt.
' eetotpolab ',eetot,polab
3844 IF(ibar(i).NE.500)
THEN
3846 IF (nre.GT.25) nre=3
3847 IF (nre.LT. 1) nre=3
3850 ave(kpl,nre)=ave(kpl,nre)+he(i)
3851 ave(kpl,30)=ave(kpl,30)+he(i)
3852 IF (ni.NE.6) ave(kpl,29)=ave(kpl,29)+he(i)
3853 avmult(kpl,nre)=avmult(kpl,nre)+1.
3854 IF (ni.NE.6) avmult(kpl,29)=avmult(kpl,29)+1.
3855 mu(kpl,ni)=mu(kpl,ni)+1
3856 IF (ich(i).NE.0)mu(kpl,9)=mu(kpl,9)+1
3859 IF (ich(i).NE.0)ave(kpl,27)=ave(kpl,27)+he(i)
3860 IF (ich(i).NE.0)avmult(kpl,27)=avmult(kpl,27)+1
3862 xfl=(
pz(i)/abs(
pz(i)))*he(i)/po
3864 IF (ixfl.LT. 1) ixfl=1
3865 IF (ixfl.GT.50) ixfl=50
3868 IF (ich(i).NE.0)yxfl(ixfl,9)=yxfl(ixfl,9)+xxxfl
3869 yxfl(ixfl,ni)=yxfl(ixfl,ni)+xxxfl
3870 yxfl(ixfl,10)=yxfl(ixfl,10)+xxxfl
3871 ptt=
px(i)**2+
py(i)**2
3872 yl=0.5*
log(abs((he(i)+
pz(i)+1.
e-10)/(he(i)-
pz(i)+1.
e-10)))
3875 IF (iylps.LT.1)iylps=1
3876 IF (iylps.GT.50)iylps=50
3877 yylps(iylps,ni)=yylps(iylps,ni)+1.
3878 yylps(iylps,10)=yylps(iylps,10)+1.
3879 IF (ich(i).NE.0)yylps(iylps,9)=yylps(iylps,9)+1.
3882 IF (iyl.GT.50) iyl=50
3883 IF (ich(i).NE.0)yyl(iyl,9)=yyl(iyl,9)+1.
3884 yyl(iyl,ni)=yyl(iyl,ni)+1.
3885 yyl(iyl,10)=yyl(iyl,10)+1.
3891 IF (ipt.GT.50) ipt=50
3892 IF (ich(i).NE.0)pty(ipt,9)=pty(ipt,9)+1./pt
3893 pty(ipt,ni)=pty(ipt,ni)+1./pt
3894 pty(ipt,10)=pty(ipt,10)+1./pt
3899 IF (im.GT.100)im=100
3900 ymult(im,i)=ymult(im,i)+1.
3905 WRITE(6, 1000)eeo,po,nhad
3908 avmult(kpl,i)=avmult(kpl,i)/nhad
3909 ave(kpl,i)=ave(kpl,i)/nhad
3912 WRITE (6,1030)avpt,navpt
3913 1030
FORMAT (
' AVERAGE PT= ',f12.4,i10)
3915 1040
FORMAT(
' PARTICLE REF,CHAR,IBAR, MASS AVERAGE',
3916 +
' ENERGY, MULTIPLICITY, INELASTICITY')
3918 ake(kpl,i)=ave(kpl,i)/eeo
3919 WRITE(6, 1050)aname(i),i,iich(i),iibar(i),
3920 + aam(i), ave(kpl,i),avmult
3921 + (kpl,i),ake(kpl,i)
3922 1050
FORMAT (
' ',a8,3i5,f10.3,3f18.6)
3926 ymult(j,i)=ymult(j,i)/nhad
3930 yxfl(j,i)=yxfl(j,i)/(nhad*dxfl)
3931 yy l(j,i)=
yy l(j,i)/(nhad*
dy)
3932 yylps(j,i)=yylps(j,i)/(nhad*
dy)
3933 pty(j,i)=pty(j,i)/(nhad*dpt)
3937 1060
FORMAT(
'1 RAPIDITY DISTRIBUTION')
3939 WRITE(6, 1070)xyl(j,1),(yyl(j,i),i=1,10)
3940 1070
FORMAT (f10.2,10e11.3)
3943 WRITE(6, 1070)xyl(j,1),(yyl(j,i),i=11,20)
3946 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3947 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3948 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
3949 CALL
plot(xyl,yyl,1000,20,50,-5.d0,
dy,0.d0,0.1d0)
3950 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3951 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3952 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
3953 CALL
plot(xyl,yylps,1000,20,50,-5.d0,
dy,0.d0,0.1d0)
3956 1080
FORMAT (
'1 LONG MOMENTUM (SCALED) DISTRIBUTION')
3958 WRITE(6, 1070)xxfl(j,1),(yxfl(j,i),i=1,10)
3961 WRITE(6, 1070)xxfl(j,1),(yxfl(j,i),i=11,20)
3965 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3966 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3967 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
3968 CALL
plot(xxfl,yxfl,1000,20,50,-1.d0,dxfl,0.d0,0.05d0)
3970 1090
FORMAT (
'1 MULTIPLICITY DISTRIBUTIONS')
3974 sumul=sumul+ymult(j,10)
3975 simul=simul+ymult(j,9)
3977 WRITE(6, 1100)(xmult(j,1),ymult(j,9),ymult(j,10),j=1,100)
3978 1100
FORMAT(f6.1,2e12.4,f6.1,2e12.4,f6.1,2e12.4,f6.1,2e12.4)
3980 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3981 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3982 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
3983 CALL
plot(xmult,ymult,1000,10,100,0.d0,1.d0,0.d0,0.01d0)
3985 xkno(i,1)=i/avmult(kpl,30)
3986 xkno(i,2)=i/avmult(kpl,27)
3987 akno(i,1)=ymult(i,10)*avmult(kpl,30)/sumul
3988 akno(i,2)=ymult(i,9)*avmult(kpl,27)/simul
3989 akno(i,1)=log10(akno(i,1)+1.
d-9)
3990 akno(i,2)=log10(akno(i,2)+1.
d-9)
3993 1110
FORMAT (
'1 KNO MULTIPLICITY DISTRIBUTIONS')
3994 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3995 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3996 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
3997 CALL
plot(xkno,akno,200,2,100,0.d0,0.08d0,-4.d0,0.05d0)
4000 ymult(j,i)=log10(ymult(j,i))
4004 yxfl(j,i)=log10(abs(yxfl(j,i)+1.
d-8))
4005 yyl(j,i)=log10(yyl(j,i)+1.
d-8)
4006 pty(j,i)=log10(pty(j,i)+1.
d-8)
4010 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4011 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4012 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4013 CALL
plot(xyl,yyl,1000,20,50,-5.d0,
dy,-3.5d0,0.05d0)
4015 WRITE(6, 1070)xxfl(j,1),(yxfl(j,i),i=1,10)
4018 WRITE(6, 1070)xxfl(j,1),(yxfl(j,i),i=11,20)
4021 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4022 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4023 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4024 CALL
plot(xxfl,yxfl,1000,20,50,-1.d0,dxfl,-4.5d0,0.05d0)
4026 1120
FORMAT (
'1 PT DISTRIBUTION DN/PTDPT')
4027 CALL
plot(ptp,pty,1000,20,50,0.d0,dpt,-2.0d0,0.05d0)
4028 IF (ipriop.EQ.1) go to 250
4030 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4031 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4032 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4033 CALL
plot(xmult,ymult,1000,10,100,0.d0,1.d0, -3.5d0,0.05d0)
4035 IF (kpl.NE.12) go to 270
4038 aaso(i,j)=log10(aaso(i,j)+1.
d-18)
4039 avmult(i,j)=log10(avmult(i,j)+1.
d-18)
4040 ake(i,j)=log10(ake(i,j)+1.
d-18)
4042 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4043 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4044 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4045 CALL
plot(aaso,avmult,360,30,12,0.d0,0.1d0,-3.d0,0.05d0)
4046 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4047 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4048 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4049 CALL
plot(aaso,ake,360,30,12,0.d0,0.1d0,-5.d0,0.05d0)
4058 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4063 parameter(nfimax=249)
4064 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
4065 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
4066 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
4073 parameter(idmax9=602)
4075 common/ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
4078 common/dpar/aname(210),am(210),ga(210),
tau(210),ich(210),ibar(210)
4080 common/dmetls/ cxs(149),cys(149),czs(149),els(149),
4081 *pls(149),is,its(149)
4082 common/ddre/
test(12)
4088 pls(i)=
sqrt(pxf(i)**2+pyf(i)**2+pzf(i)**2)
4089 IF(pls(i).NE.0.)cxs(i)=pxf(i)/pls(i)
4090 IF(pls(i).NE.0.)cys(i)=pyf(i)/pls(i)
4091 IF(pls(i).NE.0.)czs(i)=pzf(i)/pls(i)
4100 IF(istab.EQ.1) goto 30
4101 IF(istab.EQ.2) goto 50
4102 IF(istab.EQ.3) goto 40
4103 30
IF(its(ist).EQ.135.OR.its(ist).EQ.136) goto 60
4104 IF(its(ist).GE.1.AND.its(ist).LE.7) goto 60
4106 40
IF(its(ist).GE.1.AND.its(ist).LE.23) goto 60
4107 IF(its(ist).GE. 97.AND.its(ist).LE.103) goto 60
4109 IF(its(ist).EQ.109.OR.its(ist).EQ.115) goto 60
4110 IF(its(ist).GE.133.AND.its(ist).LE.136) goto 60
4112 50
IF(its(ist).GE. 1.AND.its(ist).LE. 30) goto 60
4113 IF(its(ist).GE. 97.AND.its(ist).LE.103) goto 60
4114 IF(its(ist).GE.115.AND.its(ist).LE.122) goto 60
4115 IF(its(ist).GE.131.AND.its(ist).LE.136) goto 60
4116 IF(its(ist).EQ.109) goto 60
4117 IF(its(ist).GE.137.AND.its(ist).LE.160) goto 60
4120 IF (ir.GT.nfimax)
THEN
4121 WRITE (6,1000)ir,nfimax
4122 1000
FORMAT(.GT.
' DECAY IRNFIMAX RETURN ',2i10)
4132 pxf(ir)=cxs(ist)*pls(ist)
4133 pyf(ir)=cys(ist)*pls(ist)
4134 pzf(ir)=czs(ist)*pls(ist)
4136 IF(ist.GE.1) goto 20
4140 bgam=pls(ist)/am(it)
4147 IF (vv.GT.wt(iik)) go to 90
4151 IF (it2-1.LT.0) go to 120
4154 IF(it3.EQ.0) go to 100
4155 CALL
dthrep(eco,ecm1,ecm2,ecm3,pcm1,pcm2,pcm3,cod1,cof1,sif1,
4156 *cod2,cof2,sif2,cod3,cof3,sif3,am(it1),am(it2),am(it3))
4158 100 CALL
dtwopd(eco,ecm1,ecm2,pcm1,pcm2,cod1,cof1,sif1,cod2,cof2,sif2,
4163 IF (it2-1.LT.0) go to 130
4169 CALL
dtrafo(gam,bgam,rx,ry,rz,cod1,cof1,sif1,pcm1,ecm1,
4170 *pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
4172 CALL
dtrafo(gam,bgam,rx,ry,rz,cod2,cof2,sif2,pcm2,ecm2,
4173 *pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
4174 IF (it3.LE.0) go to 130
4176 CALL
dtrafo(gam,bgam,rx,ry,rz,cod3,cof3,sif3,pcm3,ecm3,
4177 *pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
4181 IF(ir.GT.7998)
WRITE(isys,1010)
4182 1010
FORMAT(2x,
' NUMBER OF STAB. FINAL PART. IS GREATER THAN 7998')
4192 SUBROUTINE shmak(ICASE,NN,NNA,NNB,NA,NB,UMO,BIMP)
4193 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4200 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4203 dimension fnua(namx),fnub(namx),fnut(namx)
4205 dimension xb(200),bimpp(200)
4237 go to(10,30,40),icase
4261 bimpp(ib)=bimpp(ib)+1
4263 IF((nn.EQ.1).AND.(nna.EQ.1).AND.(nnb.EQ.1))
THEN
4264 CALL
sihndi(umo,1,1,singdif,sigdih)
4266 anusd=anusd + singdif/sigabs
4269 IF (intt.GT.namx)intt=namx
4271 IF (nua.GT.namx) nua=namx
4273 IF (nub.GT.namx) nub=namx
4274 fnua(nua)=fnua(nua)+1.
4275 fnut(intt)=fnut(intt)+1.
4276 fnub(nub)=fnub(nub)+1.
4298 WRITE(6,*)
' shmak(3,NN,... ) NN= ',nn
4313 WRITE(6,
'(1H1,50(1H*))')
4314 WRITE(6,
'(/10X,A/)')
' OUTPUT FROM SHMAK all events before',
4315 *
' diffraction modification'
4316 WRITE(6,
'(50(1H*))')
4317 WRITE(6,
'(A,I10)')
' NUMBER OF TOTALLY SAMPLED GLAUBER EVENTS',nn
4318 WRITE(6, 1000) bnut,bnua,bnub
4319 WRITE(6,*)
' Fraction of diffractive evnts: ',anusd
4320 1000
FORMAT(
' AVERAGE NO OF COLLISIONS BNUT,BNUA,BNUB=',3f10.3)
4321 WRITE(6,
'(/A)')
' AVERAGE NUMBERS OF DIFFERENT COLLISION TYPES'
4322 WRITE(6,
'(4(5X,A,F8.2/))')
' VAL-VAL:',bnvv,
' SEA-VAL:',bnsv,
4323 +
' VAL-SEA:',bnvs,
' SEA-SEA:',bnss
4329 1010
FORMAT (
' FNUA')
4332 fnu
a(i)=log10(fnu
a(i)+1.
d-5)
4334 CALL
plot(ann,fnu
a,namx,1,namx,0.d0,dnna,0.d0,0.05d0)
4337 1020
FORMAT (
' FNUB')
4339 fnu b(i)=log10(fnu b(i)+1.
d-5)
4341 CALL
plot(ann,fnu b,namx,1,namx,0.d0,dnnb,0.d0,0.05d0)
4344 1030
FORMAT (
' FNUT')
4346 fnu
t(i)=log10(fnu
t(i)+1.
e-5)
4348 CALL
plot(ann,fnu
t,namx,1,namx,0.d0,dnnt,0.d0,0.05d0)
4349 1040
FORMAT (10f12.2)
4358 SUBROUTINE shmak1(ICASE,NN,NNA,NNB,NA,NB,UMO,BIMP)
4359 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4366 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4369 dimension fnua(namx),fnub(namx),fnut(namx)
4371 dimension xb(200),bimpp(200)
4373 go to(10,30,40),icase
4393 bimpp(ib)=bimpp(ib)+1
4395 IF((nn.EQ.1).AND.(nna.EQ.1).AND.(nnb.EQ.1))
THEN
4396 CALL
sihndi(umo,1,1,singdif,sigdih)
4398 anusd=anusd + singdif/sigabs
4401 IF (intt.GT.namx)intt=namx
4403 IF (nua.GT.namx) nua=namx
4405 IF (nub.GT.namx) nub=namx
4406 fnua(nua)=fnua(nua)+1.
4407 fnut(intt)=fnut(intt)+1.
4408 fnub(nub)=fnub(nub)+1.
4415 WRITE(6,*)
' shmak1(3,NN,... ) NN= ',nn
4426 WRITE(6,
'(1H1,50(1H*))')
4427 WRITE(6,
'(/10X,A/)')
' OUTPUT FROM SHMAK1 after modification',
4428 *
' of Glauber events for diffractive cross section'
4429 WRITE(6,
'(50(1H*))')
4430 WRITE(6,
'(A,I10)')
' NUMBER OF TOTALLY SAMPLED GLAUBER EVENTS',nn
4431 WRITE(6, 1000) bnut,bnua,bnub
4432 1000
FORMAT(
' AVERAGE NO OF COLLISIONS BNUT,BNUA,BNUB=',3f10.3)
4433 WRITE(6,*)
' Fraction of diffractive evnts: ',anusd
4442 1010
FORMAT (
' FNUA')
4445 fnu
a(i)=log10(fnu
a(i)+1.
d-5)
4447 CALL
plot(ann,fnu
a,namx,1,namx,0.d0,dnna,0.d0,0.05d0)
4450 1020
FORMAT (
' FNUB')
4452 fnu b(i)=log10(fnu b(i)+1.
d-5)
4454 CALL
plot(ann,fnu b,namx,1,namx,0.d0,dnnb,0.d0,0.05d0)
4457 1030
FORMAT (
' FNUT')
4459 fnu
t(i)=log10(fnu
t(i)+1.
e-5)
4461 CALL
plot(ann,fnu
t,namx,1,namx,0.d0,dnnt,0.d0,0.05d0)
4462 1040
FORMAT (10f12.2)
4470 SUBROUTINE previo(RA,RB,NSTB,BMAX,BSTEP,SIG,RO,G)
4471 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4474 common/damp/ca,ci,ga
4499 rca=ga*sig/6.2831854
4504 fca=-ga*sig*ro/6.2831854
4515 WRITE(6,*)
' PREVIO(RA,RB,NSTB,BMAX,BSTEP,SIG,RO,G)',
4516 &ra,rb,nstb,bmax,bstep,sig,ro,g
4517 WRITE(6,*)
' /CA,CI,GA/ ',ca,ci,ga
4518 WRITE(6,*)
' PREVIO: RA, RB ,CI= ',ra,rb,ci
4525 SUBROUTINE profb(BSTEP,NSTAT,NA,RA,NB,RB,BSITE,NSITEB)
4526 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4534 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
4536 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
4537 +ipadis,ishmal,lpauli
4543 DOUBLE COMPLEX ca,ci
4544 COMMON /damp/ ca,ci,ga
4546 dimension bsite(0:1,nsiteb)
4548 dimension helpp(200)
4551 COMMON /sigla/siglau
4556 WRITE(6,*)
' PROFB: RA, RB = ',ra,rb
4557 WRITE(6, 1000)bstep,nstat,na,ra,nb,rb,irw,nsiteb
4558 1000
FORMAT (
' PROFB',e15.5,2i10,f15.5,i10,e15.5,2i10)
4575 IF(pi.LT.1.
d-100)go to 31
4582 IF(
xy.GT.15.) go to 20
4595 bs(i3+1)=bs(i3+1)+1.-pi
4601 bs(i)=bs(i)*(i-1)*bst/
ns
4606 bsite(1,i)=bs(i)/sumb+bsite(1,i-1)
4611 sumb=sumb*bst*6.2831854
4614 1020
FORMAT(/5x,7hsigma =,f7.3)
4618 WRITE (6,1030) help(i),helpp(i),bs(i),bsite(1,i)
4619 1030
FORMAT (f10.4,3e15.5)
4621 CALL
plot(help,bsite,50,1,50,0.d0,0.5d0,0.d0,0.01d0)
4622 CALL
plot(help,bs ,50,1,50,0.d0,0.5d0,0.d0,0.07d0)
4631 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4635 parameter(nfimax=249)
4636 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
4637 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
4638 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
4641 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
4642 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
4644 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4647 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210), ich(210),ibar
4648 +(210),k1(210),k2(210)
4660 WRITE(6,1000)ihad,i,pxf(1),pyf(1),pzf(1),hef(1),amf(1)
4661 1000
FORMAT(
' PARJET: IHAD,I,PXF(1),PYF(1),PZF(1),HEP(1),AMF(1)'/ 2i5,5
4671 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4681 IF (
a(3,i).LE.
a(3,j)) go to 20
4701 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4711 IF (
a(3,i).GE.
a(3,j)) go to 20
4730 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4777 parameter( kalgnm = 2 )
4778 parameter( anglgb = 5.0
d-16 )
4779 parameter( anglsq = 2.5
d-31 )
4780 parameter( axcssv = 0.2
d+16 )
4781 parameter( andrfl = 1.0
d-38 )
4782 parameter( avrflw = 1.0
d+38 )
4783 parameter( ainfnt = 1.0
d+30 )
4784 parameter( azrzrz = 1.0
d-30 )
4785 parameter( einfnt = +69.07755278982137
d+00 )
4786 parameter( ezrzrz = -69.07755278982137
d+00 )
4787 parameter( onemns = 0.999999999999999
d+00 )
4788 parameter( onepls = 1.000000000000001
d+00 )
4789 parameter( csnnrm = 2.0
d-15 )
4790 parameter( dmxtrn = 1.0
d+08 )
4822 parameter( zerzer = 0.
d+00 )
4823 parameter( oneone = 1.
d+00 )
4824 parameter( twotwo = 2.
d+00 )
4825 parameter( thrthr = 3.
d+00 )
4826 parameter( foufou = 4.
d+00 )
4827 parameter( fivfiv = 5.
d+00 )
4828 parameter( sixsix = 6.
d+00 )
4829 parameter( sevsev = 7.
d+00 )
4830 parameter( eigeig = 8.
d+00 )
4831 parameter( aninen = 9.
d+00 )
4832 parameter( tenten = 10.
d+00 )
4833 parameter( hlfhlf = 0.5
d+00 )
4834 parameter( onethi = oneone / thrthr )
4835 parameter( twothi = twotwo / thrthr )
4836 parameter( pipipi = 3.1415926535897932270
d+00 )
4837 parameter( eneper = 2.7182818284590452354
d+00 )
4838 parameter( sqrent = 1.6487212707001281468
d+00 )
4880 parameter( clight = 2.99792458
d+10 )
4881 parameter( avogad = 6.0221367
d+23 )
4882 parameter( amelgr = 9.1093897
d-28 )
4883 parameter( plckbr = 1.05457266
d-27 )
4884 parameter( elccgs = 4.8032068
d-10 )
4885 parameter( elcmks = 1.60217733
d-19 )
4886 parameter( amugrm = 1.6605402
d-24 )
4887 parameter( ammumu = 0.113428913
d+00 )
4902 parameter( alpfsc = 7.2973530791728595
d-03 )
4903 parameter( fscto2 = 5.3251361962113614
d-05 )
4904 parameter( fscto3 = 3.8859399018437826
d-07 )
4905 parameter( fscto4 = 2.8357075508200407
d-09 )
4906 parameter( plabrc = 0.197327053
d+00 )
4907 parameter( amelct = 0.51099906
d-03 )
4908 parameter( amugev = 0.93149432
d+00 )
4909 parameter( ammuon = 0.105658389
d+00 )
4910 parameter( rclsel = 2.8179409183694872
d-13 )
4911 parameter( gevmev = 1.0
d+03 )
4912 parameter( emvgev = 1.0
d-03 )
4913 parameter( raddeg = 180.
d+00 / pipipi )
4914 parameter( degrad = pipipi / 180.
d+00 )
4941 parameter( lunin = 5 )
4942 parameter( lunout = 6 )
4943 parameter( lunerr = 66 )
4944 parameter( lunber = 14 )
4945 parameter( lunech = 8 )
4946 parameter( lunflu = 86 )
4947 parameter( lungeo = 16 )
4948 parameter( lunpgs = 12 )
4949 parameter( lunran = 2 )
4950 parameter( lunxsc = 81 )
4951 parameter( lunrdb = 1 )
4980 parameter( mxxrgn = 500 )
4981 parameter( mxxmdf = 56 )
4982 parameter( mxxmde = 50 )
4983 parameter( mfstck = 1000 )
4984 parameter( mestck = 100 )
4985 parameter( nallwp = 39 )
4986 parameter( mpdpdx = 8 )
4987 parameter( icomax = 180 )
4988 parameter( nstbis = 304 )
4989 parameter( idmaxp = 210 )
4993 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210),
4994 + ich(210),ibar(210),k1(210),k2(210)
5002 DATA (am(k),k=1,85) /
5003 & .9383
d+00, .9383
d+00, amelct , amelct , .0000
d+00,
5004 & .0000
d+00, .0000
d+00, .9396
d+00, .9396
d+00, ammuon ,
5005 & ammuon , .4977
d+00, .1396
d+00, .1396
d+00, .4936
d+00,
5006 & .4936
d+00, .1116
d+01, .1116
d+01, .4977
d+00, .1197
d+01,
5007 & .1189
d+01, .1193
d+01, .1350
d+00, .4977
d+00, .4977
d+00,
5008 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5009 & .5488
d+00, .7669
d+00, .7700
d+00, .7669
d+00, .7820
d+00,
5010 & .8921
d+00, .8962
d+00, .8921
d+00, .8962
d+00, .1300
d+01,
5011 & .1300
d+01, .1300
d+01, .1300
d+01, .1421
d+01, .1421
d+01,
5012 & .1421
d+01, .1421
d+01, .1383
d+01, .1384
d+01, .1387
d+01,
5013 & .1820
d+01, .2030
d+01, .1231
d+01, .1232
d+01, .1233
d+01,
5014 & .1234
d+01, .1675
d+01, .1675
d+01, .1675
d+01, .1675
d+01,
5015 & .1500
d+01, .1500
d+01, .1515
d+01, .1515
d+01, .1775
d+01,
5016 & .1775
d+01, .1231
d+01, .1232
d+01, .1233
d+01, .1234
d+01,
5017 & .1675
d+01, .1675
d+01, .1675
d+01, .1675
d+01, .1515
d+01,
5018 & .1515
d+01, .2500
d+01, .4890
d+00, .4890
d+00, .4890
d+00,
5019 & .1300
d+01, .1300
d+01, .1300
d+01, .1300
d+01, .2200
d+01 /
5020 DATA (am(k),k=86,183) /
5021 & .2200
d+01, .2200
d+01, .2200
d+01, .1700
d+01, .1700
d+01,
5022 & .1700
d+01, .1700
d+01, .1820
d+01, .2030
d+01, .9575
d+00,
5023 & .1019
d+01, .1315
d+01, .1321
d+01, .1189
d+01, .1193
d+01,
5024 & .1197
d+01, .1315
d+01, .1321
d+01, .1383
d+01, .1384
d+01,
5025 & .1387
d+01, .1532
d+01, .1535
d+01, .1672
d+01, .1383
d+01,
5026 & .1384
d+01, .1387
d+01, .1532
d+01, .1535
d+01, .1672
d+01,
5027 & .1865
d+01, .1869
d+01, .1869
d+01, .1865
d+01, .1969
d+01,
5028 & .1969
d+01, .2980
d+01, .2007
d+01, .2010
d+01, .2010
d+01,
5029 & .2007
d+01, .2113
d+01, .2113
d+01, .3686
d+01, .3097
d+01,
5030 & .1777
d+01, .1777
d+01, .0000
d+00, .0000
d+00, .0000
d+00,
5031 & .0000
d+00, .2285
d+01, .2460
d+01, .2460
d+01, .2452
d+01,
5032 & .2453
d+01, .2454
d+01, .2560
d+01, .2560
d+01, .2730
d+01,
5033 & .3610
d+01, .3610
d+01, .3790
d+01, .2285
d+01, .2460
d+01,
5034 & .2460
d+01, .2452
d+01, .2453
d+01, .2454
d+01, .2560
d+01,
5035 & .2560
d+01, .2730
d+01, .3610
d+01, .3610
d+01, .3790
d+01,
5036 & .2490
d+01, .2490
d+01, .2490
d+01, .2610
d+01, .2610
d+01,
5037 & .2770
d+01, .3670
d+01, .3670
d+01, .3850
d+01, .4890
d+01,
5038 & .2490
d+01, .2490
d+01, .2490
d+01, .2610
d+01, .2610
d+01,
5039 & .2770
d+01, .3670
d+01, .3670
d+01, .3850
d+01, .4890
d+01,
5040 & .1250
d+01, .1250
d+01, .1250
d+01 /
5041 DATA ( am( i ), i = 184,210 ) /
5042 & 1.44000000000000
d+00, 1.44000000000000
d+00, 1.30000000000000
d+00,
5043 & 1.30000000000000
d+00, 1.30000000000000
d+00, 1.40000000000000
d+00,
5044 & 1.46000000000000
d+00, 1.46000000000000
d+00, 1.46000000000000
d+00,
5045 & 1.46000000000000
d+00, 1.60000000000000
d+00, 1.60000000000000
d+00,
5046 & 1.66000000000000
d+00, 1.66000000000000
d+00, 1.66000000000000
d+00,
5047 & 1.66000000000000
d+00, 1.66000000000000
d+00, 1.66000000000000
d+00,
5048 & 1.95000000000000
d+00, 1.95000000000000
d+00, 1.95000000000000
d+00,
5049 & 1.95000000000000
d+00, 2.25000000000000
d+00, 2.25000000000000
d+00,
5050 & 1.44000000000000
d+00, 1.44000000000000
d+00, 0.00000000000000
d+00/
5054 DATA (
tau(k),k=1,183) /
5055 & .1000
d+19, .1000
d+19, .1000
d+19, .1000
d+19, .1000
d+19,
5056 & .1000
d+19, .1000
d+19, .9180
d+03, .9180
d+03, .2200
d-05,
5057 & .2200
d-05, .5200
d-07, .2600
d-07, .2600
d-07, .1200
d-07,
5058 & .1200
d-07, .2600
d-09, .2600
d-09, .9000
d-10, .1500
d-09,
5059 & .8000
d-10, .5000
d-14, .8000
d-16, .0000
d+00, .0000
d+00,
5061 & .0000
d+00, .3000
d-09, .1700
d-09, .8000
d-10, .1000
d-13,
5062 & .1500
d-09, .3000
d-09, .1700
d-09, .0000
d+00, .0000
d+00,
5063 & .0000
d+00, .0000
d+00, .0000
d+00, .1000
d-09, .0000
d+00,
5064 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .1000
d-09,
5065 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5066 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5067 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5068 & .9000
d-11, .9000
d-11, .9000
d-11, .9000
d-11, .1000
d+19,
5069 & .1000
d+19, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5071 & .0000
d+00, .0000
d+00, .0000
d+00 /
5072 DATA (
tau( i ), i = 184,210 ) /
5073 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5074 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5075 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5076 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5077 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5078 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5079 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5080 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5081 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00/
5085 DATA (ga(k),k= 1,85) /
5087 & .8500
d-06, .1520
d+00, .1520
d+00, .1520
d+00, .1000
d-01,
5088 & .7900
d-01, .7900
d-01, .7900
d-01, .7900
d-01, .4500
d+00,
5089 & .4500
d+00, .4500
d+00, .4500
d+00, .1080
d+00, .1080
d+00,
5090 & .1080
d+00, .1080
d+00, .5000
d-01, .5000
d-01, .5000
d-01,
5091 & .8500
d-01, .1800
d+00, .1150
d+00, .1150
d+00, .1150
d+00,
5092 & .1150
d+00, .2000
d+00, .2000
d+00, .2000
d+00, .2000
d+00,
5093 & .2000
d+00, .2000
d+00, .1000
d+00, .1000
d+00, .2000
d+00,
5094 & .2000
d+00, .1150
d+00, .1150
d+00, .1150
d+00, .1150
d+00,
5095 & .2000
d+00, .2000
d+00, .2000
d+00, .2000
d+00, .1000
d+00,
5096 & .1000
d+00, .2000
d+00, .1000
d+00, .1000
d+00, .1000
d+00,
5097 & .1000
d+00, .1000
d+00, .1000
d+00, .1000
d+00, .2000
d+00 /
5098 DATA (ga(k),k= 86,183) /
5099 & .2000
d+00, .2000
d+00, .2000
d+00, .1500
d+00, .1500
d+00,
5100 & .1500
d+00, .1500
d+00, .8500
d-01, .1800
d+00, .2000
d-02,
5101 & .4000
d-02, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5102 & .0000
d+00, .0000
d+00, .0000
d+00, .3400
d-01, .3400
d-01,
5103 & .3600
d-01, .9000
d-02, .9000
d-02, .0000
d+00, .3400
d-01,
5104 & .3400
d-01, .3600
d-01, .9000
d-02, .9000
d-02, .0000
d+00,
5105 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5106 & .0000
d+00, .0000
d+00, .5000
d-02, .2000
d-02, .2000
d-02,
5107 & .5000
d-02, .2000
d-02, .2000
d-02, .2000
d-03, .7000
d-03,
5109 & .3000
d+00, .3000
d+00, .3000
d+00 /
5110 DATA ( ga( i ), i = 184,210 ) /
5111 & 2.00000000000000
d-01, 2.00000000000000
d-01, 3.00000000000000
d-01,
5112 & 3.00000000000000
d-01, 3.00000000000000
d-01, 2.70000000000000
d-01,
5113 & 2.50000000000000
d-01, 2.50000000000000
d-01, 2.50000000000000
d-01,
5114 & 2.50000000000000
d-01, 1.50000000000000
d-01, 1.50000000000000
d-01,
5115 & 1.00000000000000
d-01, 1.00000000000000
d-01, 1.00000000000000
d-01,
5116 & 1.00000000000000
d-01, 1.00000000000000
d-01, 1.00000000000000
d-01,
5117 & 6.00000000000000
d-02, 6.00000000000000
d-02, 6.00000000000000
d-02,
5118 & 6.00000000000000
d-02, 5.50000000000000
d-02, 5.50000000000000
d-02,
5119 & 2.00000000000000
d-01, 2.00000000000000
d-01, 0.00000000000000
d+00/
5128 DATA (aname(k),k=1,85) /
5129 &
'P ',
'AP ',
'E- ',
'E+ ',
'NUE ',
5130 &
'ANUE ',
'GAM ',
'NEU ',
'ANEU ',
'MUE+ ',
5131 &
'MUE- ',
'K0L ',
'PI+ ',
'PI- ',
'K+ ',
5132 &
'K- ',
'LAM ',
'ALAM ',
'K0S ',
'SIGM- ',
5133 &
'SIGM+ ',
'SIGM0 ',
'PI0 ',
'K0 ',
'AK0 ',
5134 &
'BLANK ',
'BLANK ',
'BLANK ',
'BLANK ',
'BLANK ',
5135 &
'ETA550 ',
'RHO+77 ',
'RHO077 ',
'RHO-77 ',
'OM0783 ',
5136 &
'K*+892 ',
'K*0892 ',
'K*-892 ',
'AK*089 ',
'KA+125 ',
5137 &
'KA0125 ',
'KA-125 ',
'AKA012 ',
'K*+142 ',
'K*0142 ',
5138 &
'K*-142 ',
'AK*014 ',
'S+1385 ',
'S01385 ',
'S-1385 ',
5139 &
'L01820 ',
'L02030 ',
'N*++12 ',
'N*+ 12 ',
'N*012 ',
5140 &
'N*-12 ',
'N*++16 ',
'N*+16 ',
'N*016 ',
'N*-16 ',
5141 &
'N*+14 ',
'N*014 ',
'N*+15 ',
'N*015 ',
'N*+18 ',
5142 &
'N*018 ',
'AN--12 ',
'AN*-12 ',
'AN*012 ',
'AN*+12 ',
5143 &
'AN--16 ',
'AN*-16 ',
'AN*016 ',
'AN*+16 ',
'AN*-15 ',
5144 &
'AN*015 ',
'DE*=24 ',
'RPI+49 ',
'RPI049 ',
'RPI-49 ',
5145 &
'PIN++ ',
'PIN+0 ',
'PIN+- ',
'PIN-0 ',
'PPPI ' /
5146 DATA (aname(k),k=86,183) /
5147 &
'PNPI ',
'APPPI ',
'APNPI ',
'K+PPI ',
'K-PPI ',
5148 &
'K+NPI ',
'K-NPI ',
'S+1820 ',
'S-2030 ',
'ETA* ',
5149 &
'PHI ',
'TETA0 ',
'TETA- ',
'ASIG- ',
'ASIG0 ',
5150 &
'ASIG+ ',
'ATETA0 ',
'ATETA+ ',
'SIG*+ ',
'SIG*0 ',
5151 &
'SIG*- ',
'TETA*0 ',
'TETA* ',
'OMEGA- ',
'ASIG*- ',
5152 &
'ASIG*0 ',
'ASIG*+ ',
'ATET*0 ',
'ATET*+ ',
'OMEGA+ ',
5153 &
'D0 ',
'D+ ',
'D- ',
'AD0 ',
'DS+ ',
5154 &
'DS- ',
'ETAC ',
'D*0 ',
'D*+ ',
'D*- ',
5155 &
'AD*0 ',
'DS*+ ',
'DS*- ',
'CHI1C ',
'JPSI ',
5156 &
'TAU+ ',
'TAU- ',
'NUET ',
'ANUET ',
'NUEM ',
5157 &
'ANUEM ',
'LAMC+ ',
'XIC+ ',
'XIC0 ',
'SIGC++ ',
5158 &
'SIGC+ ',
'SIGC0 ',
'S+ ',
'S0 ',
'T0 ',
5159 &
'XU++ ',
'XD+ ',
'XS+ ',
'ALAMC- ',
'AXIC- ',
5160 &
'AXIC0 ',
'ASIGC-- ',
'ASIGC- ',
'ASIGC0 ',
'AS- ',
5161 &
'AS0 ',
'AT0 ',
'AXU-- ',
'AXD- ',
'AXS ',
5162 &
'C1*++ ',
'C1*+ ',
'C1*0 ',
'S*+ ',
'S*0 ',
5163 &
'T*0 ',
'XU*++ ',
'XD*+ ',
'XS*+ ',
'TETA++ ',
5164 &
'AC1*-- ',
'AC1*- ',
'AC1*0 ',
'AS*- ',
'AS*0 ',
5165 &
'AT*0 ',
'AXU*-- ',
'AXD*- ',
'AXS*- ',
'ATET-- ',
5166 &
'RO ',
'R+ ',
'R- ' /
5167 DATA ( aname( i ), i = 184,210 ) /
5168 &
'AN*-14 ',
'AN*014 ',
'PI+130 ',
'PI0130 ',
'PI-130 ',
'F01400 ',
5169 &
'K*+146 ',
'K*-146 ',
'K*0146 ',
'AK0146 ',
'L01600 ',
'AL0160 ',
5170 &
'S+1660 ',
'S01660 ',
'S-1660 ',
'AS-166 ',
'AS0166 ',
'AS+166 ',
5171 &
'X01950 ',
'X-1950 ',
'AX0195 ',
'AX+195 ',
'OM-225 ',
'AOM+22 ',
5172 &
'N*+14 ',
'N*014 ',
'BLANK '/
5176 DATA ( ich( i ), i = 1,210 ) /
5177 & 1, -1, -1, 1, 0, 0, 0, 0, 0, 1, -1, 0, 1, -1, 1,
5178 & -1, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5179 & 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0,
5180 & -1, 0, 1, 0, -1, 0, 0, 2, 1, 0, -1, 2, 1, 0, -1,
5181 & 1, 0, 1, 0, 1, 0, -2, -1, 0, 1, -2, -1, 0, 1, -1,
5182 & 0, 1, 1, 0, -1, 2, 1, 0, -1, 2, 1, 0, -1, 2, 0,
5183 & 1, -1, 1, -1, 0, 0, 0, -1, -1, 0, 1, 0, 1, 1, 0,
5184 & -1, 0, -1, -1, -1, 0, 1, 0, 1, 1, 0, 1, -1, 0, 1,
5185 & -1, 0, 0, 1, -1, 0, 1, -1, 0, 0, 1, -1, 0, 0, 0,
5186 & 0, 1, 1, 0, 2, 1, 0, 1, 0, 0, 2, 1, 1, -1, -1,
5187 & 0, -2, -1, 0, -1, 0, 0, -2, -1, -1, 2, 1, 0, 1, 0,
5188 & 0, 2, 1, 1, 2, -2, -1, 0, -1, 0, 0, -2, -1, -1, -2,
5189 & 0, 1, -1, -1, 0, 1, 0, -1, 0, 1, -1, 0, 0, 0, 0,
5190 & 1, 0, -1, -1, 0, 1, 0, -1, 0, 1, -1, 1, 1, 0, 0/
5194 DATA ( ibar( i ), i = 1,210 ) /
5195 & 1, -1, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0,
5196 & 0, 1, -1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0,
5197 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5198 & 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
5199 & 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
5200 & -1, 2, 0, 0, 0, 1, 1, 1, 1, 2, 2, 0, 0, 1, 1,
5201 & 1, 1, 1, 1, 0, 0, 1, 1, -1, -1, -1, -1, -1, 1, 1,
5202 & 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0,
5203 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5204 & 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -1,
5205 & -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1,
5206 & 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
5207 & 0, 0, 0, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1,
5208 & 1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 0/
5213 DATA k1/ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 16, 17,
5214 & 18, 24, 30, 34, 38, 40, 41, 43, 44, 136, 138, 330, 327, 328,
5215 & 2*330, 46, 51, 52, 54, 55, 58,
5216 & 60, 62, 64, 66, 68, 70, 72, 74, 82, 90, 98, 106, 109, 112, 114,
5217 & 123, 140, 141, 143, 145, 146, 150, 157, 164, 168, 174, 180, 187,
5218 & 194, 202, 210, 211, 213, 215, 216, 220, 227, 234, 238, 245, 252,
5219 & 254, 255, 256, 257, 259, 262, 265, 267, 269, 272, 276, 279, 282,
5220 & 286, 290, 293, 299, 331, 335, 339, 340, 341, 343, 344, 345, 346,
5221 & 347, 350, 353, 356, 358, 360, 363, 366, 369, 372, 374, 376, 379,
5222 & 383, 385, 387, 391, 394, 397, 400, 402, 405, 408, 410, 412, 414,
5223 & 417, 420, 425, 430, 431, 432, 433, 434, 448, 452, 457, 458, 459,
5224 & 460, 461, 462, 466, 468, 470, 472, 486, 490, 495, 496, 497, 498,
5225 & 499, 500, 504, 506, 508, 510, 511, 512, 513, 514, 515, 516, 517,
5226 & 518, 519, 522, 523, 524, 525, 526, 527, 528, 529, 530, 531, 534,
5227 & 537, 539, 541, 547, 553, 558, 563, 568, 572, 573, 574, 575, 576,
5228 & 577, 578, 579, 580, 581, 582, 583, 584, 585, 586, 587, 588, 589,
5234 DATA k2/ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 15, 16, 17,
5235 & 23, 29, 31, 35, 39, 40, 42, 43, 45, 137, 139, 330, 327, 328,
5236 & 2* 330, 50, 51, 53, 54, 57,
5237 & 59, 61, 63, 65, 67, 69, 71, 73, 81, 89, 97, 105, 108, 111, 113,
5238 & 122, 135, 140, 142, 144, 145, 149, 156, 163, 167, 173, 179, 186,
5239 & 193, 201, 209, 210, 212, 214, 215, 219, 226, 233, 237, 244, 251,
5240 & 253, 254, 255, 256, 258, 261, 264, 266, 268, 271, 275, 278, 281,
5241 & 285, 289, 292, 298, 307, 334, 338, 339, 340, 342, 343, 344, 345,
5242 & 346, 349, 352, 355, 357, 359, 362, 365, 368, 371, 373, 375, 378,
5243 & 382, 384, 386, 390, 393, 396, 399, 401, 404, 407, 409, 411, 413,
5244 & 416, 419, 424, 429, 430, 431, 432, 433, 447, 451, 456, 457, 458,
5245 & 459, 460, 461, 465, 467, 469, 471, 485, 489, 494, 495, 496, 497,
5246 & 498, 499, 503, 505, 507, 509, 510, 511, 512, 513, 514, 515, 516,
5247 & 517, 518, 521, 522, 523, 524, 525, 526, 527, 528, 529, 530, 533,
5248 & 536, 538, 540, 546, 552, 557, 562, 567, 571, 572, 573, 574, 575,
5249 & 576, 577, 578, 579, 580, 581, 582, 583, 584, 585, 586, 587, 588,
5250 & 589, 595, 601, 602 /
5259 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5264 parameter(idmax9=602)
5266 common/ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
5277 DATA (zkname(k),k= 1, 85) /
5278 &
'P ',
'AP ',
'E- ',
'E+ ',
'NUE ',
5279 &
'ANUE ',
'GAM ',
'PE-NUE ',
'APEANU ',
'EANUNU ',
5280 &
'E-NUAN ',
'3PI0 ',
'PI+-0 ',
'PIMUNU ',
'PIE-NU ',
5281 &
'MU+NUE ',
'MU-NUE ',
'MU+NUE ',
'PI+PI0 ',
'PI++- ',
5282 &
'PI+00 ',
'M+P0NU ',
'E+P0NU ',
'MU-NU ',
'PI-0 ',
5283 &
'PI+-- ',
'PI-00 ',
'M-P0NU ',
'E-P0NU ',
'PPI- ',
5284 &
'NPI0 ',
'PD-NUE ',
'PM-NUE ',
'APPI+ ',
'ANPI0 ',
5285 &
'APE+NU ',
'APM+NU ',
'PI+PI- ',
'PI0PI0 ',
'NPI- ',
5286 &
'PPI0 ',
'NPI+ ',
'LAGA ',
'GAGA ',
'GAE+E- ',
5287 &
'GAGA ',
'GAGAP0 ',
'PI000 ',
'PI+-0 ',
'PI+-GA ',
5288 &
'PI+0 ',
'PI+- ',
'PI00 ',
'PI-0 ',
'PI+-0 ',
5289 &
'PI+- ',
'PI0GA ',
'K+PI0 ',
'K0PI+ ',
'KOPI0 ',
5290 &
'K+PI- ',
'K-PI0 ',
'AK0PI- ',
'AK0PI0 ',
'K-PI+ ',
5291 &
'K+PI0 ',
'K0PI+ ',
'K0PI0 ',
'K+PI- ',
'K-PI0 ',
5292 &
'K0PI- ',
'AK0PI0 ',
'K-PI+ ',
'K+PI0 ',
'K0PI+ ',
5293 &
'K+89P0 ',
'K08PI+ ',
'K+RO77 ',
'K0RO+7 ',
'K+OM07 ',
5294 &
'K+E055 ',
'K0PI0 ',
'K+PI+ ',
'K089P0 ',
'K+8PI- ' /
5295 DATA (zkname(k),k= 86,170) /
5296 &
'K0R077 ',
'K+R-77 ',
'K+R-77 ',
'K0OM07 ',
'K0E055 ',
5297 &
'K-PI0 ',
'K0PI- ',
'K-89P0 ',
'AK08P- ',
'K-R077 ',
5298 &
'AK0R-7 ',
'K-OM07 ',
'K-E055 ',
'AK0PI0 ',
'K-PI+ ',
5299 &
'AK08P0 ',
'K-8PI+ ',
'AK0R07 ',
'AK0OM7 ',
'AK0E05 ',
5300 &
'LA0PI+ ',
'SI0PI+ ',
'SI+PI0 ',
'LA0PI0 ',
'SI+PI- ',
5301 &
'SI-PI+ ',
'LA0PI- ',
'SI0PI- ',
'NEUAK0 ',
'PK- ',
5302 &
'SI+PI- ',
'SI0PI0 ',
'SI-PI+ ',
'LA0ET0 ',
'S+1PI- ',
5303 &
'S-1PI+ ',
'SO1PI0 ',
'NEUAK0 ',
'PK- ',
'LA0PI0 ',
5304 &
'LA0OM0 ',
'LA0RO0 ',
'SI+RO- ',
'SI-RO+ ',
'SI0RO0 ',
5305 &
'LA0ET0 ',
'SI0ET0 ',
'SI+PI- ',
'SI-PI+ ',
'SI0PI0 ',
5306 &
'K0S ',
'K0L ',
'K0S ',
'K0L ',
'P PI+ ',
5307 &
'P PI0 ',
'N PI+ ',
'P PI- ',
'N PI0 ',
'N PI- ',
5308 &
'P PI+ ',
'N*#PI0 ',
'N*+PI+ ',
'PRHO+ ',
'P PI0 ',
5309 &
'N PI+ ',
'N*#PI- ',
'N*+PI0 ',
'N*0PI+ ',
'PRHO0 ',
5310 &
'NRHO+ ',
'P PI- ',
'N PI0 ',
'N*+PI- ',
'N*0PI0 ',
5311 &
'N*-PI+ ',
'PRHO- ',
'NRHO0 ',
'N PI- ',
'N*0PI- ',
5312 &
'N*-PI0 ',
'NRHO- ',
'PETA0 ',
'N*#PI- ',
'N*+PI0 ' /
5313 DATA (zkname(k),k=171,255) /
5314 &
'N*0PI+ ',
'PRHO0 ',
'NRHO+ ',
'NETA0 ',
'N*+PI- ',
5315 &
'N*0PI0 ',
'N*-PI+ ',
'PRHO- ',
'NRHO0 ',
'P PI0 ',
5316 &
'N PI+ ',
'N*#PI- ',
'N*+PI0 ',
'N*0PI+ ',
'PRHO0 ',
5317 &
'NRHO+ ',
'P PI- ',
'N PI0 ',
'N*+PI- ',
'N*0PI0 ',
5318 &
'N*-PI+ ',
'PRHO- ',
'NRHO0 ',
'P PI0 ',
'N PI+ ',
5319 &
'PRHO0 ',
'NRHO+ ',
'LAMK+ ',
'S+ K0 ',
'S0 K+ ',
5320 &
'PETA0 ',
'P PI- ',
'N PI0 ',
'PRHO- ',
'NRHO0 ',
5321 &
'LAMK0 ',
'S0 K0 ',
'S- K+ ',
'NETA/ ',
'APPI- ',
5322 &
'APPI0 ',
'ANPI- ',
'APPI+ ',
'ANPI0 ',
'ANPI+ ',
5323 &
'APPI- ',
'AN*=P0 ',
'AN*-P- ',
'APRHO- ',
'APPI0 ',
5324 &
'ANPI- ',
'AN*=P+ ',
'AN*-P0 ',
'AN*0P- ',
'APRHO0 ',
5325 &
'ANRHO- ',
'APPI+ ',
'ANPI0 ',
'AN*-P+ ',
'AN*0P0 ',
5326 &
'AN*+P- ',
'APRHO+ ',
'ANRHO0 ',
'ANPI+ ',
'AN*0P+ ',
5327 &
'AN*+P0 ',
'ANRHO+ ',
'APPI0 ',
'ANPI- ',
'AN*=P+ ',
5328 &
'AN*-P0 ',
'AN*0P- ',
'APRHO0 ',
'ANRHO- ',
'APPI+, ',
5329 &
'ANPI0 ',
'AN*-P+ ',
'AN*0P0 ',
'AN*+P- ',
'APRHO+ ',
5330 &
'ANRHO0 ',
'PN*014 ',
'NN*=14 ',
'PI+0 ',
'PI+- ' /
5331 DATA (zkname(k),k=256,340) /
5332 &
'PI-0 ',
'P+0 ',
'N++ ',
'P+- ',
'P00 ',
5333 &
'N+0 ',
'N+- ',
'N00 ',
'P-0 ',
'N-0 ',
5334 &
'P-- ',
'PPPI0 ',
'PNPI+ ',
'PNPI0 ',
'PPPI- ',
5335 &
'NNPI+ ',
'APPPI0 ',
'APNPI+ ',
'ANNPI0 ',
'ANPPI- ',
5336 &
'APNPI0 ',
'APPPI- ',
'ANNPI- ',
'K+PPI0 ',
'K+NPI+ ',
5337 &
'K0PPI0 ',
'K-PPI0 ',
'K-NPI+ ',
'AKPPI- ',
'AKNPI0 ',
5338 &
'K+NPI0 ',
'K+PPI- ',
'K0PPI0 ',
'K0NPI+ ',
'K-NPI0 ',
5339 &
'K-PPI- ',
'AKNPI- ',
'PAK0 ',
'SI+PI0 ',
'SI0PI+ ',
5340 &
'SI+ETA ',
'S+1PI0 ',
'S01PI+ ',
'NEUK- ',
'LA0PI- ',
5341 &
'SI-OM0 ',
'LA0RO- ',
'SI0RO- ',
'SI-RO0 ',
'SI-ET0 ',
5342 &
'SI0PI- ',
'SI-0 ',
'BLANC ',
'BLANC ',
'BLANC ',
5343 &
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
5344 &
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
5345 &
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
5346 &
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
5347 &
'EPI+- ',
'EPI00 ',
'GAPI+- ',
'GAGA* ',
'K+- ',
5348 &
'KLKS ',
'PI+-0 ',
'EGA ',
'LPI0 ',
'LPI ' /
5349 DATA (zkname(k),k=341,425) /
5350 &
'APPI0 ',
'ANPI- ',
'ALAGA ',
'ANPI ',
'ALPI0 ',
5351 &
'ALPI+ ',
'LAPI+ ',
'SI+PI0 ',
'SI0PI+ ',
'LAPI0 ',
5352 &
'SI+PI- ',
'SI-PI+ ',
'LAPI- ',
'SI-PI0 ',
'SI0PI- ',
5353 &
'TE0PI0 ',
'TE-PI+ ',
'TE0PI- ',
'TE-PI0 ',
'TE0PI ',
5354 &
'TE-PI ',
'LAK- ',
'ALPI- ',
'AS-PI0 ',
'AS0PI- ',
5355 &
'ALPI0 ',
'AS+PI- ',
'AS-PI+ ',
'ALPI+ ',
'AS+PI0 ',
5356 &
'AS0PI+ ',
'AT0PI0 ',
'AT+PI- ',
'AT0PI+ ',
'AT+PI0 ',
5357 &
'AT0PI ',
'AT+PI ',
'ALK+ ',
'K-PI+ ',
'K-PI+0 ',
5358 &
'K0PI+- ',
'K0PI0 ',
'K-PI++ ',
'AK0PI+ ',
'K+PI-- ',
5359 &
'K0PI- ',
'K+PI- ',
'K+PI-0 ',
'AKPI-+ ',
'AK0PI0 ',
5360 &
'ETAPIF ',
'K++- ',
'K+AK0 ',
'ETAPI- ',
'K--+ ',
5361 &
'K-K0 ',
'PI00 ',
'PI+- ',
'GAGA ',
'D0PI0 ',
5362 &
'D0GA ',
'D0PI+ ',
'D+PI0 ',
'DFGA ',
'AD0PI- ',
5363 &
'D-PI0 ',
'D-GA ',
'AD0PI0 ',
'AD0GA ',
'F+GA ',
5364 &
'F+GA ',
'F-GA ',
'F-GA ',
'PSPI+- ',
'PSPI00 ',
5365 &
'PSETA ',
'E+E- ',
'MUE+- ',
'PI+-0 ',
'M+NN ',
5366 &
'E+NN ',
'RHO+NT ',
'PI+ANT ',
'K*+ANT ',
'M-NN ' /
5367 DATA (zkname(k),k=426,510) /
5368 &
'E-NN ',
'RHO-NT ',
'PI-NT ',
'K*-NT ',
'NUET ',
5369 &
'ANUET ',
'NUEM ',
'ANUEM ',
'SI+ETA ',
'SI+ET* ',
5370 &
'PAK0 ',
'TET0K+ ',
'SI*+ET ',
'N*+AK0 ',
'N*++K- ',
5371 &
'LAMRO+ ',
'SI0RO+ ',
'SI+RO0 ',
'SI+OME ',
'PAK*0 ',
5372 &
'N*+AK* ',
'N*++K* ',
'SI+AK0 ',
'TET0PI ',
'SI+AK* ',
5373 &
'TET0RO ',
'SI0AK* ',
'SI+K*- ',
'TET0OM ',
'TET-RO ',
5374 &
'SI*0AK ',
'C0+PI+ ',
'C0+PI0 ',
'C0+PI- ',
'A+GAM ',
5375 &
'A0GAM ',
'TET0AK ',
'TET0K* ',
'OM-RO+ ',
'OM-PI+ ',
5376 &
'C1++AK ',
'A+PI+ ',
'C0+AK0 ',
'A0PI+ ',
'A+AK0 ',
5377 &
'T0PI+ ',
'ASI-ET ',
'ASI-E* ',
'APK0 ',
'ATET0K ',
5378 &
'ASI*-E ',
'AN*-K0 ',
'AN*--K ',
'ALAMRO ',
'ASI0RO ',
5379 &
'ASI-RO ',
'ASI-OM ',
'APK*0 ',
'AN*-K* ',
'AN*--K ',
5380 &
'ASI-K0 ',
'ATETPI ',
'ASI-K* ',
'ATETRO ',
'ASI0K* ',
5381 &
'ASI-K* ',
'ATE0OM ',
'ATE+RO ',
'ASI*0K ',
'AC-PI- ',
5382 &
'AC-PI0 ',
'AC-PI+ ',
'AA-GAM ',
'AA0GAM ',
'ATET0K ',
5383 &
'ATE0K* ',
'AOM+RO ',
'AOM+PI ',
'AC1--K ',
'AA-PI- ',
5384 &
'AC0-K0 ',
'AA0PI- ',
'AA-K0 ',
'AT0PI- ',
'C1++GA ' /
5385 DATA (zkname(k),k=511,540) /
5386 &
'C1++GA ',
'C10GAM ',
'S+GAM ',
'S0GAM ',
'T0GAM ',
5387 &
'XU++GA ',
'XD+GAM ',
'XS+GAM ',
'A+AKPI ',
'T02PI+ ',
5388 &
'C1++2K ',
'AC1--G ',
'AC1-GA ',
'AC10GA ',
'AS-GAM ',
5389 &
'AS0GAM ',
'AT0GAM ',
'AXU--G ',
'AXD-GA ',
'AXS-GA ',
5390 &
'AA-KPI ',
'AT02PI ',
'AC1--K ',
'RH-PI+ ',
'RH+PI- ',
5391 &
'RH3PI0 ',
'RH0PI+ ',
'RH+PI0 ',
'RH0PI- ',
'RH-PI0 ' /
5392 DATA (zkname(i),i=541,602)/
5393 &
'APETA ',
'AN=P+ ',
'AN-PO ',
'ANOPO ',
'APRHO0',
'ANRHO-',
'ANETA ',
5394 &
'AN-P+ ',
'AN0PO ',
'AN+P- ',
'APRHO+',
'ANRHO0',
'RH0PI+',
'RH+PI0',
5395 &
'3PI+00',
'3PI-++',
'F0PI+ ',
'RH+PI-',
'RH0PI0',
'3PI000',
'3PI0+-',
5396 &
'F0PI0 ',
'RH0PI-',
'RH-PI0',
'3PI-00',
'3PI--+',
'F0PI- ',
'PI0PI0',
5397 &
'PI+PI-',
'K+K- ',
'K0AK0 ',
'L01600',
'AL0160',
'K*+146',
'K*-146',
5398 &
'K*0146',
'AK0146',
'S+1660',
'S01660',
'S-1660',
'AS-166',
'AS0166',
5399 &
'AS+166',
'X01690',
'X-1690',
'AX0169',
'AX+169',
'OM-225',
'AOM+22',
5400 &
'N*PPI0',
'N*NPI+',
'N*P2P0',
'N*PP+-',
'N*D+P0',
'N*D0P+',
'N*NPI0',
5401 &
'N*PPI-',
'N*N2P0',
'N*NP+-',
'N*D+P-',
'N*D0P0',
'BLANK '/
5405 DATA (wt(k),k= 1, 85) /
5406 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5407 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5408 & .1000
d+01, .2100
d+00, .1200
d+00, .2700
d+00, .4000
d+00,
5409 & .1000
d+01, .1000
d+01, .6400
d+00, .2100
d+00, .6000
d-01,
5410 & .2000
d-01, .3000
d-01, .4000
d-01, .6400
d+00, .2100
d+00,
5411 & .6000
d-01, .2000
d-01, .3000
d-01, .4000
d-01, .6400
d+00,
5412 & .3600
d+00, .0000
d+00, .0000
d+00, .6400
d+00, .3600
d+00,
5413 & .0000
d+00, .0000
d+00, .6900
d+00, .3100
d+00, .1000
d+01,
5414 & .5200
d+00, .4800
d+00, .1000
d+01, .9900
d+00, .1000
d-01,
5415 & .3800
d+00, .3000
d-01, .3000
d+00, .2400
d+00, .5000
d-01,
5416 & .1000
d+01, .1000
d+01, .0000
d+00, .1000
d+01, .9000
d+00,
5417 & .1000
d-01, .9000
d-01, .3300
d+00, .6700
d+00, .3300
d+00,
5418 & .6700
d+00, .3300
d+00, .6700
d+00, .3300
d+00, .6700
d+00,
5419 & .3300
d+00, .6700
d+00, .3300
d+00, .6700
d+00, .3300
d+00,
5420 & .6700
d+00, .3300
d+00, .6700
d+00, .1900
d+00, .3800
d+00,
5421 & .9000
d-01, .2000
d+00, .3000
d-01, .4000
d-01, .5000
d-01,
5422 & .2000
d-01, .1900
d+00, .3800
d+00, .9000
d-01, .2000
d+00 /
5423 DATA (wt(k),k= 86,170) /
5424 & .3000
d-01, .4000
d-01, .5000
d-01, .2000
d-01, .1900
d+00,
5425 & .3800
d+00, .9000
d-01, .2000
d+00, .3000
d-01, .4000
d-01,
5426 & .5000
d-01, .2000
d-01, .1900
d+00, .3800
d+00, .9000
d-01,
5427 & .2000
d+00, .3000
d-01, .4000
d-01, .5000
d-01, .2000
d-01,
5428 & .8800
d+00, .6000
d-01, .6000
d-01, .8800
d+00, .6000
d-01,
5429 & .6000
d-01, .8800
d+00, .1200
d+00, .1900
d+00, .1900
d+00,
5430 & .1600
d+00, .1600
d+00, .1700
d+00, .3000
d-01, .3000
d-01,
5431 & .3000
d-01, .4000
d-01, .1000
d+00, .1000
d+00, .2000
d+00,
5432 & .1200
d+00, .1000
d+00, .4000
d-01, .4000
d-01, .5000
d-01,
5433 & .7500
d-01, .7500
d-01, .3000
d-01, .3000
d-01, .4000
d-01,
5434 & .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00, .1000
d+01,
5435 & .6700
d+00, .3300
d+00, .3300
d+00, .6700
d+00, .1000
d+01,
5436 & .2500
d+00, .2700
d+00, .1800
d+00, .3000
d+00, .1700
d+00,
5437 & .8000
d-01, .1800
d+00, .3000
d-01, .2400
d+00, .2000
d+00,
5438 & .1000
d+00, .8000
d-01, .1700
d+00, .2400
d+00, .3000
d-01,
5439 & .1800
d+00, .1000
d+00, .2000
d+00, .2500
d+00, .1800
d+00,
5440 & .2700
d+00, .3000
d+00, .5000
d+00, .3000
d+00, .1250
d+00 /
5444 DATA (wt(k),k=171,255) /
5445 & .7500
d-01, .0000
d+00, .0000
d+00, .5000
d+00, .7500
d-01,
5446 & .1250
d+00, .3000
d+00, .0000
d+00, .0000
d+00, .1800
d+00,
5447 & .3700
d+00, .1300
d+00, .8000
d-01, .4000
d-01, .7000
d-01,
5448 & .1300
d+00, .3700
d+00, .1800
d+00, .4000
d-01, .8000
d-01,
5449 & .1300
d+00, .1300
d+00, .7000
d-01, .7000
d-01, .1300
d+00,
5450 & .2300
d+00, .4700
d+00, .5000
d-01, .2000
d-01, .1000
d-01,
5451 & .2000
d-01, .1300
d+00, .7000
d-01, .4700
d+00, .2300
d+00,
5452 & .5000
d-01, .1000
d-01, .2000
d-01, .2000
d-01, .1000
d+01,
5453 & .6700
d+00, .3300
d+00, .3300
d+00, .6700
d+00, .1000
d+01,
5454 & .2500
d+00, .2700
d+00, .1800
d+00, .3000
d+00, .1700
d+00,
5455 & .8000
d-01, .1800
d+00, .3000
d-01, .2400
d+00, .2000
d+00,
5456 & .1000
d+00, .8000
d-01, .1700
d+00, .2400
d+00, .3000
d-01,
5457 & .1800
d+00, .1000
d+00, .2000
d+00, .2500
d+00, .1800
d+00,
5458 & .2700
d+00, .3000
d+00, .1800
d+00, .3700
d+00, .1300
d+00,
5459 & .8000
d-01, .4000
d-01, .7000
d-01, .1300
d+00, .3700
d+00,
5460 & .1800
d+00, .4000
d-01, .8000
d-01, .1300
d+00, .1300
d+00,
5461 & .7000
d-01, .5000
d+00, .5000
d+00, .1000
d+01, .1000
d+01 /
5462 DATA (wt(k),k=256,340) /
5463 & .1000
d+01, .8000
d+00, .2000
d+00, .6000
d+00, .3000
d+00,
5464 & .1000
d+00, .6000
d+00, .3000
d+00, .1000
d+00, .8000
d+00,
5465 & .2000
d+00, .3300
d+00, .6700
d+00, .6600
d+00, .1700
d+00,
5466 & .1700
d+00, .3200
d+00, .1700
d+00, .3200
d+00, .1900
d+00,
5467 & .3300
d+00, .3300
d+00, .3400
d+00, .3000
d+00, .5000
d-01,
5468 & .6500
d+00, .3800
d+00, .1200
d+00, .3800
d+00, .1200
d+00,
5469 & .3800
d+00, .1200
d+00, .3800
d+00, .1200
d+00, .3000
d+00,
5470 & .5000
d-01, .6500
d+00, .3800
d+00, .2500
d+00, .2500
d+00,
5471 & .2000
d-01, .5000
d-01, .5000
d-01, .2000
d+00, .2000
d+00,
5472 & .1200
d+00, .1000
d+00, .7000
d-01, .7000
d-01, .1400
d+00,
5473 & .5000
d-01, .5000
d-01, .1000
d+01, .1000
d+01, .1000
d+01,
5474 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5475 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5476 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5477 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5478 & .4800
d+00, .2400
d+00, .2600
d+00, .2000
d-01, .4700
d+00,
5479 & .3500
d+00, .1500
d+00, .3000
d-01, .1000
d+01, .1000
d+01 /
5480 DATA (wt(k),k=341,425) /
5481 & .5200
d+00, .4800
d+00, .1000
d+01, .1000
d+01, .1000
d+01,
5482 & .1000
d+01, .9000
d+00, .5000
d-01, .5000
d-01, .9000
d+00,
5483 & .5000
d-01, .5000
d-01, .9000
d+00, .5000
d-01, .5000
d-01,
5484 & .3300
d+00, .6700
d+00, .6700
d+00, .3300
d+00, .2500
d+00,
5485 & .2500
d+00, .5000
d+00, .9000
d+00, .5000
d-01, .5000
d-01,
5486 & .9000
d+00, .5000
d-01, .5000
d-01, .9000
d+00, .5000
d-01,
5487 & .5000
d-01, .3300
d+00, .6700
d+00, .6700
d+00, .3300
d+00,
5488 & .2500
d+00, .2500
d+00, .5000
d+00, .1000
d+00, .5000
d+00,
5489 & .1600
d+00, .2400
d+00, .7000
d+00, .3000
d+00, .7000
d+00,
5490 & .3000
d+00, .1000
d+00, .5000
d+00, .1600
d+00, .2400
d+00,
5491 & .3000
d+00, .4000
d+00, .3000
d+00, .3000
d+00, .4000
d+00,
5492 & .3000
d+00, .4900
d+00, .4900
d+00, .2000
d-01, .5500
d+00,
5493 & .4500
d+00, .6800
d+00, .3000
d+00, .2000
d-01, .6800
d+00,
5494 & .3000
d+00, .2000
d-01, .5500
d+00, .4500
d+00, .9000
d+00,
5495 & .1000
d+00, .9000
d+00, .1000
d+00, .6000
d+00, .3000
d+00,
5496 & .1000
d+00, .1000
d+00, .1000
d+00, .8000
d+00, .2800
d+00,
5497 & .2800
d+00, .3500
d+00, .7000
d-01, .2000
d-01, .2800
d+00 /
5498 DATA (wt(k),k=426,510) /
5499 & .2800
d+00, .3500
d+00, .7000
d-01, .2000
d-01, .1000
d+01,
5500 & .1000
d+01, .1000
d+01, .1000
d+01, .2000
d-01, .3000
d-01,
5501 & .7000
d-01, .2000
d-01, .2000
d-01, .4000
d-01, .1300
d+00,
5502 & .7000
d-01, .6000
d-01, .6000
d-01, .2000
d+00, .1400
d+00,
5503 & .4000
d-01, .1000
d+00, .2500
d+00, .3000
d-01, .3000
d+00,
5504 & .4200
d+00, .2200
d+00, .3500
d+00, .1900
d+00, .1600
d+00,
5505 & .8000
d-01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5506 & .1000
d+01, .3700
d+00, .2000
d+00, .3600
d+00, .7000
d-01,
5507 & .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00,
5508 & .5000
d+00, .2000
d-01, .3000
d-01, .7000
d-01, .2000
d-01,
5509 & .2000
d-01, .4000
d-01, .1300
d+00, .7000
d-01, .6000
d-01,
5510 & .6000
d-01, .2000
d+00, .1400
d+00, .4000
d-01, .1000
d+00,
5511 & .2500
d+00, .3000
d-01, .3000
d+00, .4200
d+00, .2200
d+00,
5512 & .3500
d+00, .1900
d+00, .1600
d+00, .8000
d-01, .1000
d+01,
5513 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .3700
d+00,
5514 & .2000
d+00, .3600
d+00, .7000
d-01, .5000
d+00, .5000
d+00,
5515 & .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00, .1000
d+01 /
5516 DATA (wt(k),k=511,540) /
5517 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5518 & .1000
d+01, .1000
d+01, .1000
d+01, .3000
d+00, .3000
d+00,
5519 & .4000
d+00, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5520 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5521 & .3000
d+00, .3000
d+00, .4000
d+00, .3300
d+00, .3300
d+00,
5522 & .3400
d+00, .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00 /
5524 DATA (wt(i),i=541,602) / .0
d+00, .3334
d+00, .2083
d+00, 2*.125
d+00,
5525 & .2083
d+00, .0
d+00, .125
d+00, .2083
d+00, .3334
d+00, .2083
d+00,
5526 & .125
d+00, 0.2
d+00, 0.2
d+00, 0.3
d+00, 0.3
d+00, 0.0
d+00, 0.2
d+00,
5527 & 0.2
d+00, 0.3
d+00, 0.3
d+00, 0.0
d+00, 0.2
d+00, 0.2
d+00, 0.3
d+00,
5528 & 0.3
d+00, 0.0
d+00, 0.31
d+00, 0.62
d+00, 0.035
d+00, 0.035
d+00,
5529 & 18*1.
d+00, 0.5
d+00, 0.16
d+00, 2*0.12
d+00, 2*0.05
d+00, 0.5
d+00,
5530 & 0.16
d+00, 2*0.12
d+00, 2*0.05
d+00, 1.
d+00 /
5534 DATA (nzk(k,1),k= 1,170) /
5535 & 1, 2, 3, 4, 5, 6, 7, 1, 2, 4,
5536 & 3, 23, 13, 13, 13, 10, 11, 10, 13, 13,
5537 & 13, 10, 4, 11, 14, 14, 14, 11, 3, 1,
5538 & 8, 1, 1, 2, 9, 2, 2, 13, 23, 8,
5539 & 1, 8, 17, 7, 7, 7, 23, 23, 13, 13,
5540 & 13, 13, 23, 14, 13, 13, 23, 15, 24, 24,
5541 & 15, 16, 25, 25, 16, 15, 24, 24, 15, 16,
5542 & 24, 25, 16, 15, 24, 36, 37, 15, 24, 15,
5543 & 15, 24, 15, 37, 36, 24, 15, 24, 24, 16,
5544 & 24, 38, 39, 16, 25, 16, 16, 25, 16, 39,
5545 & 38, 25, 16, 25, 25, 17, 22, 21, 17, 21,
5546 & 20, 17, 22, 8, 1, 21, 22, 20, 17, 48,
5547 & 50, 49, 8, 1, 17, 17, 17, 21, 20, 22,
5548 & 17, 22, 21, 20, 22, 19, 12, 19, 12, 1,
5549 & 1, 8, 1, 8, 8, 1, 53, 54, 1, 1,
5550 & 8, 53, 54, 55, 1, 8, 1, 8, 54, 55,
5551 & 56, 1, 8, 8, 55, 56, 8, 1, 53, 54 /
5552 DATA (nzk(k,1),k=171,340) /
5553 & 55, 1, 8, 8, 54, 55, 56, 1, 8, 1,
5554 & 8, 53, 54, 55, 1, 8, 1, 8, 54, 55,
5555 & 56, 1, 8, 1, 8, 1, 8, 17, 21, 22,
5556 & 1, 1, 8, 1, 8, 17, 22, 20, 8, 2,
5557 & 2, 9, 2, 9, 9, 2, 67, 68, 2, 2,
5558 & 9, 67, 68, 69, 2, 9, 2, 9, 68, 69,
5559 & 70, 2, 9, 9, 69, 70, 9, 2, 9, 67,
5560 & 68, 69, 2, 9, 2, 9, 68, 69, 70, 2,
5561 & 9, 1, 8, 13, 13, 14, 1, 8, 1, 1,
5562 & 8, 8, 8, 1, 8, 1, 1, 1, 1, 1,
5563 & 8, 2, 2, 9, 9, 2, 2, 9, 15, 15,
5564 & 24, 16, 16, 25, 25, 15, 15, 24, 24, 16,
5565 & 16, 25, 1, 21, 22, 21, 48, 49, 8, 17,
5566 & 20, 17, 22, 20, 20, 22, 20, 0, 0, 0,
5567 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5568 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5569 & 31, 31, 13, 7, 15, 12, 13, 31, 17, 17 /
5570 DATA (nzk(k,1),k=341,510) /
5571 & 2, 9, 18, 9, 18, 18, 17, 21, 22, 17,
5572 & 21, 20, 17, 20, 22, 97, 98, 97, 98, 97,
5573 & 98, 17, 18, 99, 100, 18, 101, 99, 18, 101,
5574 & 100, 102, 103, 102, 103, 102, 103, 18, 16, 16,
5575 & 24, 24, 16, 25, 15, 24, 15, 15, 25, 25,
5576 & 31, 15, 15, 31, 16, 16, 23, 13, 7, 116,
5577 & 116, 116, 117, 117, 119, 118, 118, 119, 119, 120,
5578 & 120, 121, 121, 130, 130, 130, 4, 10, 13, 10,
5579 & 4, 32, 13, 36, 11, 3, 34, 14, 38, 133,
5580 & 134, 135, 136, 21, 21, 1, 97, 104, 54, 53,
5581 & 17, 22, 21, 21, 1, 54, 53, 21, 97, 21,
5582 & 97, 22, 21, 97, 98, 105, 137, 137, 137, 138,
5583 & 139, 97, 97, 109, 109, 140, 138, 137, 139, 138,
5584 & 145, 99, 99, 2, 102, 110, 68, 67, 18, 100,
5585 & 99, 99, 2, 68, 67, 99, 102, 99, 102, 100,
5586 & 99, 102, 103, 111, 149, 149, 149, 150, 151, 113,
5587 & 113, 115, 115, 152, 150, 149, 151, 150, 157, 140 /
5588 DATA (nzk(k,1),k=511,540) /
5589 & 141, 142, 143, 144, 145, 146, 147, 148, 138, 145,
5590 & 140, 152, 153, 154, 155, 156, 157, 158, 159, 160,
5591 & 150, 157, 152, 34, 32, 33, 33, 32, 33, 34 /
5592 DATA (nzk(i,1),i=541,602) / 2, 67, 68, 69, 2, 9, 9, 68, 69,
5593 & 70, 2, 9, 33, 32, 13, 14, 189, 32, 34, 23, 23, 189, 33, 34, 14,
5594 & 14, 189, 23, 13, 15, 24, 36, 38, 37, 39, 194, 195, 196, 197,
5595 & 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 1, 8, 1, 1, 54,
5596 & 55, 8, 1, 8, 8, 54, 55, 210/
5597 DATA (nzk(k,2),k= 1,170) /
5598 & 0, 0, 0, 0, 0, 0, 0, 3, 4, 6,
5599 & 5, 23, 14, 11, 3, 5, 5, 5, 23, 13,
5600 & 23, 23, 23, 5, 23, 13, 23, 23, 23, 14,
5601 & 23, 3, 11, 13, 23, 4, 10, 14, 23, 14,
5602 & 23, 13, 7, 7, 4, 7, 7, 23, 14, 14,
5603 & 23, 14, 23, 23, 14, 14, 7, 23, 13, 23,
5604 & 14, 23, 14, 23, 13, 23, 13, 23, 14, 23,
5605 & 14, 23, 13, 23, 13, 23, 13, 33, 32, 35,
5606 & 31, 23, 14, 23, 14, 33, 34, 35, 31, 23,
5607 & 14, 23, 14, 33, 34, 35, 31, 23, 13, 23,
5608 & 13, 33, 32, 35, 31, 13, 13, 23, 23, 14,
5609 & 13, 14, 14, 25, 16, 14, 23, 13, 31, 14,
5610 & 13, 23, 25, 16, 23, 35, 33, 34, 32, 33,
5611 & 31, 31, 14, 13, 23, 0, 0, 0, 0, 13,
5612 & 23, 13, 14, 23, 14, 13, 23, 13, 78, 23,
5613 & 13, 14, 23, 13, 79, 78, 14, 23, 14, 23,
5614 & 13, 80, 79, 14, 14, 23, 80, 31, 14, 23 /
5615 DATA (nzk(k,2),k=171,340) /
5616 & 13, 79, 78, 31, 14, 23, 13, 80, 79, 23,
5617 & 13, 14, 23, 13, 79, 78, 14, 23, 14, 23,
5618 & 13, 80, 79, 23, 13, 33, 32, 15, 24, 15,
5619 & 31, 14, 23, 34, 33, 24, 24, 15, 31, 14,
5620 & 23, 14, 13, 23, 13, 14, 23, 14, 80, 23,
5621 & 14, 13, 23, 14, 79, 80, 13, 23, 13, 23,
5622 & 14, 78, 79, 13, 13, 23, 78, 23, 14, 13,
5623 & 23, 14, 79, 80, 13, 23, 13, 23, 14, 78,
5624 & 79, 62, 61, 23, 14, 23, 13, 13, 13, 23,
5625 & 13, 13, 23, 14, 14, 14, 1, 8, 8, 1,
5626 & 8, 1, 8, 8, 1, 8, 1, 8, 1, 8,
5627 & 1, 1, 8, 1, 8, 8, 1, 1, 8, 8,
5628 & 1, 8, 25, 23, 13, 31, 23, 13, 16, 14,
5629 & 35, 34, 34, 33, 31, 14, 23, 0, 0, 0,
5630 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5631 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5632 & 13, 23, 14, 7, 16, 19, 14, 7, 23, 14 /
5633 DATA (nzk(k,2),k=341,510) /
5634 & 23, 14, 7, 13, 23, 13, 13, 23, 13, 23,
5635 & 14, 13, 14, 23, 14, 23, 13, 14, 23, 14,
5636 & 23, 16, 14, 23, 14, 23, 14, 13, 13, 23,
5637 & 13, 23, 14, 13, 23, 13, 23, 15, 13, 13,
5638 & 13, 23, 13, 13, 14, 14, 14, 14, 14, 23,
5639 & 13, 16, 25, 14, 15, 24, 23, 14, 7, 23,
5640 & 7, 13, 23, 7, 14, 23, 7, 23, 7, 7,
5641 & 7, 7, 7, 13, 23, 31, 3, 11, 14, 135,
5642 & 5, 134, 134, 134, 136, 6, 133, 133, 133, 0,
5643 & 0, 0, 0, 31, 95, 25, 15, 31, 95, 16,
5644 & 32, 32, 33, 35, 39, 39, 38, 25, 13, 39,
5645 & 32, 39, 38, 35, 32, 39, 13, 23, 14, 7,
5646 & 7, 25, 37, 32, 13, 25, 13, 25, 13, 25,
5647 & 13, 31, 95, 24, 16, 31, 24, 15, 34, 34,
5648 & 33, 35, 37, 37, 36, 24, 14, 37, 34, 37,
5649 & 36, 35, 34, 37, 14, 23, 13, 7, 7, 24,
5650 & 39, 34, 14, 24, 14, 24, 14, 24, 14, 7 /
5651 DATA (nzk(k,2),k=511,540) /
5652 & 7, 7, 7, 7, 7, 7, 7, 7, 25, 13,
5653 & 25, 7, 7, 7, 7, 7, 7, 7, 7, 7,
5654 & 24, 14, 24, 13, 14, 23, 13, 23, 14, 23 /
5655 DATA (nzk(i,2),i=541,602) / 31, 13, 23, 14, 79, 80, 31, 13, 23,
5656 & 14, 78, 79, 13, 23, 23, 13, 13, 14, 13, 23, 13, 23, 14, 23, 23,
5657 & 14, 14, 23, 14, 16, 25,
5658 & 4*23, 14*0, 23, 13, 23, 13, 23, 13, 23, 14,
5659 & 23, 13, 14, 23, 0 /
5660 DATA (nzk(k,3),k= 1,170) /
5661 & 0, 0, 0, 0, 0, 0, 0, 5, 6, 5,
5662 & 6, 23, 23, 5, 5, 0, 0, 0, 0, 14,
5663 & 23, 5, 5, 0, 0, 14, 23, 5, 5, 0,
5664 & 0, 5, 5, 0, 0, 5, 5, 0, 0, 0,
5665 & 0, 0, 0, 0, 3, 0, 7, 23, 23, 7,
5666 & 0, 0, 0, 0, 23, 0, 0, 0, 0, 0,
5668 DATA (nzk(k,3),k=171,340) /
5670 & 0, 0, 0, 0, 0, 0, 23, 13, 14, 23,
5671 & 23, 14, 23, 23, 23, 14, 23, 13, 23, 14,
5672 & 13, 23, 13, 23, 14, 23, 14, 14, 23, 13,
5673 & 13, 23, 13, 14, 23, 23, 14, 23, 13, 23,
5674 & 14, 14, 0, 0, 0, 0, 0, 0, 0, 0,
5676 & 14, 23, 7, 0, 0, 0, 23, 0, 0, 0 /
5677 DATA (nzk(k,3),k=341,510) /
5679 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 23,
5680 & 14, 0, 13, 0, 14, 0, 0, 23, 13, 0,
5681 & 0, 15, 0, 0, 16, 0, 0, 0, 0, 0,
5682 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5683 & 0, 0, 0, 14, 23, 0, 0, 0, 23, 134,
5684 & 134, 0, 0, 0, 133, 133, 0, 0, 0, 0,
5686 DATA (nzk(k,3),k=511,540) /
5687 & 0, 0, 0, 0, 0, 0, 0, 0, 13, 13,
5688 & 25, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5689 & 14, 14, 24, 0, 0, 0, 0, 0, 0, 0 /
5690 DATA (nzk(i,3),i=541,602) / 12*0, 2*0, 23, 13, 0, 2*0, 23, 14, 0,
5691 & 2*0, 23, 13, 0, 4*0, 18*0, 2*0, 23, 14, 2*0, 2*0, 23, 14, 2*0, 0/
5712 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5714 parameter(lout=6,llook=9)
5716 COMPLEX*16 czero,cone,ctwo
5717 parameter(
zero=0.0d0,
one=1.0d0,two=2.0d0,three=3.0d0,
5718 & onethi=
one/three,tiny25=1.0
d-25)
5719 parameter(twopi = 6.283185307179586454
d+00,
5721 & gev2mb = 0.38938d0,
5722 & gev2fm = 0.1972d0,
5723 & alphem =
one/137.0d0,
5733 COMMON /dpar/ aname(210),aam(210),gam(210),
tau(210),iich(210),
5734 & iibar(210),ka1(210),ka2(210)
5736 parameter(ncompx=1,
neb=50)
5737 COMMON /dshmm/ rash,rbsh(ncompx),bmax(ncompx),bstep(ncompx),
5738 & sigsh,rosh,gsh,bsite(0:
neb,ncompx,ksiteb),
5740 COMMON /glaber/ ecmnn(
neb),ecmnow,
5749 COMMON /vdmpar/ rl2,epspol,intrge(2),idpdf,modega,ishad(3)
5750 COMMON /glapar/ jstatb
5753 COMMON /damp/ ca,ci,ga
5754 COMMON /xsecnu/ecmuu,ecmoo,ngritt,nevtt
5755 COMMON /kglaub/jglaub
5757 parameter(maxncl = 210)
5758 COMPLEX*16 pp11,pp12,pp21,pp22,
5759 & ompp11,ompp12,ompp21,ompp22
5760 dimension coop1(3,maxncl),coot1(3,maxncl),
5761 & coop2(3,maxncl),coot2(3,maxncl),
5762 & bprod(ksiteb),sigshh(
neb),
5766 Write(6,*)
' XSGLAU(NA,NB,IJPROJ,NTARG)',
5768 WRITE(6,*)
'/XSECNU/ECMUU,ECMOO,NGRITT,NEVTT',
5769 &ecmuu,ecmoo,ngritt,nevtt
5773 ctwo = dcmplx(two,
zero)
5777 dellog=(log10(ecmoo)-log10(ecmuu))/(ngritt-1)
5778 deldel=10.d0**dellog
5780 DO 1123 ieee=1,ngritt
5785 WRITE(6,*)
'IE,EC111,S',ie,ec111,
s
5794 rash = rnucle*dble(na)**onethi
5795 rbsh(ntarg) = rnucle*dble(nb)**onethi
5797 IF(na.EQ.9)rash=2.52d0
5798 IF(na.EQ.10)rash=2.45d0
5799 IF(na.EQ.11)rash=2.37d0
5800 IF(na.EQ.12)rash=2.45d0
5801 IF(na.EQ.13)rash=2.44d0
5802 IF(na.EQ.14)rash=2.55d0
5803 IF(na.EQ.15)rash=2.58d0
5804 IF(na.EQ.16)rash=2.71d0
5805 IF(na.EQ.17)rash=2.66d0
5806 IF(na.EQ.18)rash=2.71d0
5807 IF(nb.EQ.9)rbsh(ntarg)=2.52d0
5808 IF(nb.EQ.10)rbsh(ntarg)=2.45d0
5809 IF(nb.EQ.11)rbsh(ntarg)=2.37d0
5810 IF(nb.EQ.12)rbsh(ntarg)=2.45d0
5811 IF(nb.EQ.13)rbsh(ntarg)=2.44d0
5812 IF(nb.EQ.14)rbsh(ntarg)=2.55d0
5813 IF(nb.EQ.15)rbsh(ntarg)=2.58d0
5814 IF(nb.EQ.16)rbsh(ntarg)=2.71d0
5815 IF(nb.EQ.17)rbsh(ntarg)=2.66d0
5816 IF(nb.EQ.18)rbsh(ntarg)=2.71d0
5819 bmax(ntarg) = 4.0d0*(rash+rbsh(ntarg))
5820 bstep(ntarg)= bmax(ntarg)/dble(nsiteb-1)
5823 IF (ijproj.LE.12)
THEN
5824 bslope = 8.5d0*(1.0d0+0.065d0*
log(
s))
5825 IF (ecmnn(ie).LE.3.0d0)
THEN
5827 ELSEIF ((ecmnn(ie).GT.3.0d0).AND.(ecmnn(ie).LE.50.d0))
THEN
5828 rosh = -0.63d0+0.175d0*
log(ecmnn(ie))
5829 ELSEIF (ecmnn(ie).GT.50.0d0)
THEN
5833 bslope = 6.0d0*(1.0d0+0.065d0*
log(
s))
5838 elab = (
s-aam(ijproj)**2-amp2)/(two*amp)
5839 elabb(ie)=elab/1000.
5840 plab =
sqrt( (elab-aam(ijproj))*(elab+aam(ijproj)) )
5842 sigsh =
dshpto(ijproj,plab)/10.d0
5843 sigshh(ie)=sigsh*10.d0
5844 WRITE(6,*)
' NSTATB,NSITEB,RASH,RBSH(NTARG),BMAX(NTARG),
5845 &BSLOPE,ROSH,SIGSH,ECM ELAB',
5846 & nstatb,nsiteb,rash,rbsh(ntarg),bmax(ntarg),
5847 &bslope,rosh,sigsh,ec111,elab
5850 bsite( 0,ntarg,i) =
zero
5851 bsite(ie,ntarg,i) =
zero
5866 facn =
one/dble(nstatb)
5882 CALL
conuclx(coot1,nb,rbsh(ntarg),1)
5884 CALL
conuclx(coot2,nb,rbsh(ntarg),1)
5895 b = dble(ib)*bstep(ntarg)
5896 facb = 10.0d0*twopi*b*bstep(ntarg)
5910 gsh = 10.0d0/(two*bslope*gev2mb)
5913 rca = ga*sigsh/twopi
5915 ca = dcmplx(rca,fca)
5923 x11 = b+coot1(1,inb)-coop1(1,ina)
5924 y11 = coot1(2,inb)-coop1(2,ina)
5925 xy11 = ga*(x11*x11+y11*y11)
5926 x12 = b+coot2(1,inb)-coop1(1,ina)
5927 y12 = coot2(2,inb)-coop1(2,ina)
5928 xy12 = ga*(x12*x12+y12*y12)
5929 x21 = b+coot1(1,inb)-coop2(1,ina)
5930 y21 = coot1(2,inb)-coop2(2,ina)
5931 xy21 = ga*(x21*x21+y21*y21)
5932 x22 = b+coot2(1,inb)-coop2(1,ina)
5933 y22 = coot2(2,inb)-coop2(2,ina)
5934 xy22 = ga*(x22*x22+y22*y22)
5935 IF (xy11.LE.15.0d0)
THEN
5936 c = cone-ca*
exp(-xy11)
5939 IF (abs(ar).LT.tiny25) ar =
zero
5940 IF (abs(ai).LT.tiny25) ai =
zero
5941 pp11 = dcmplx(ar,ai)
5945 shi = shi+
log(ar*ar+ai*ai)
5947 IF (xy12.LE.15.0d0)
THEN
5948 c = cone-ca*
exp(-xy12)
5951 IF (abs(ar).LT.tiny25) ar =
zero
5952 IF (abs(ai).LT.tiny25) ai =
zero
5953 pp12 = dcmplx(ar,ai)
5956 IF (xy21.LE.15.0d0)
THEN
5957 c = cone-ca*
exp(-xy21)
5960 IF (abs(ar).LT.tiny25) ar =
zero
5961 IF (abs(ai).LT.tiny25) ai =
zero
5962 pp21 = dcmplx(ar,ai)
5965 IF (xy22.LE.15.0d0)
THEN
5966 c = cone-ca*
exp(-xy22)
5969 IF (abs(ar).LT.tiny25) ar =
zero
5970 IF (abs(ai).LT.tiny25) ai =
zero
5971 pp22 = dcmplx(ar,ai)
5979 ompp11 = ompp11+(cone-pp11)
5980 ompp21 = ompp21+(cone-pp21)
5983 ompp12 = ompp12+(cone-pp12)
5984 ompp22 = ompp22+(cone-pp22)
5986 stotm = dble(ompp11+ompp22)
5987 selam = dble(ompp11*dconjg(ompp22))
5989 sqepm = dble(ompp11*dconjg(ompp21))-selam
5990 sqetm = dble(ompp11*dconjg(ompp12))-selam
5991 sqe2m = dble(ompp11*dconjg(ompp11))-selam-sqepm-sqetm
5993 stotb = stotb+facm*stotm
5994 selab = selab+facm*selam
5995 IF (nb.GT.1) sqepb = sqepb+facm*sqepm
5996 IF (na.GT.1) sqetb = sqetb+facm*sqetm
5997 IF ((na.GT.1).AND.(nb.GT.1)) sqe2b = sqe2b+facm*sqe2m
5998 sprob = sprob+facm*sprom
6002 stotn = stotn+facb*stotb
6003 selan = selan+facb*selab
6004 sqepn = sqepn+facb*sqepb
6005 sqetn = sqetn+facb*sqetb
6006 sqe2n = sqe2n+facb*sqe2b
6007 spron = spron+facb*sprob
6008 bprod(ib+1)= bprod(ib+1)+facn*facb*sprob
6012 stot = stot +facn*stotn
6013 stot2 = stot2+facn*stotn**2
6014 sela = sela +facn*selan
6015 sela2 = sela2+facn*selan**2
6016 sqep = sqep +facn*sqepn
6017 sqep2 = sqep2+facn*sqepn**2
6018 sqet = sqet +facn*sqetn
6019 sqet2 = sqet2+facn*sqetn**2
6020 sqe2 = sqe2 +facn*sqe2n
6021 sqe22 = sqe22+facn*sqe2n**2
6022 spro = spro +facn*spron
6023 spro2 = spro2+facn*spron**2
6040 WRITE(6,*)
' STOT,SELA ,SQEP,SQET,SQE2,SPRO ',
6041 & stot,sela ,sqep,sqet,sqe2,spro
6043 xetot(ie) =
sqrt(abs(stot2-stot**2)/dble(nstatb-1))
6044 xeela(ie) =
sqrt(abs(sela2-sela**2)/dble(nstatb-1))
6045 xeqep(ie) =
sqrt(abs(sqep2-sqep**2)/dble(nstatb-1))
6046 xeqet(ie) =
sqrt(abs(sqet2-sqet**2)/dble(nstatb-1))
6047 xeqe2(ie) =
sqrt(abs(sqe22-sqe2**2)/dble(nstatb-1))
6048 xepro(ie) =
sqrt(abs(spro2-spro**2)/dble(nstatb-1))
6049 WRITE(6,*)
' XETOT,XEELA,XEQEP,XEQET,XEQE2,XEPRO ',
6050 & xetot(ie),xeela(ie),xeqep(ie),
6051 &xeqet(ie),xeqe2(ie),xepro(ie)
6054 bsite(ie,ntarg,i) = bprod(i)/spro+bsite(ie,ntarg,i-1)
6056 & bsite(0,ntarg,i) = bprod(i)/spro+bsite(0,ntarg,i-1)
6058 WRITE (6,*)
' ECMNN,ELABB,SIGSHH,SIGTO,SIGEL,SIGIN,SIGSD'
6061 sigto(i)=
dshnto(1,1,ecmnn(i))
6062 sigel(i)=
dshnel(1,1,ecmnn(i))
6063 sigin(i)=
siinel(1,1,ecmnn(i))
6064 sigsd(i)=
sippsd(ecmnn(i))
6065 CALL
sihndi(ecmnn(i),1,1,sigdif(i),sigdih)
6066 WRITE (6,
'(2F18.4,6F11.3)')ecmnn(i),elabb(i),sigshh(i),
6067 & sigto(i),sigel(i),sigin(i),sigsd(i),sigdif(i)
6069 WRITE (6,*)
' ECMNN,ELABB,XSQEP,XEQEP,XSQET,XEQET,XSQE2,XEQE2'
6071 WRITE (6,
'(2F18.4,6F11.3)')ecmnn(i),elabb(i),xsqep(i),xeqep(i),
6072 * xsqet(i),xeqet(i),xsqe2(i),xeqe2(i)
6074 WRITE (6,*)
' ECMNN,ELABB,XSTOT,XETOT,XSELA,XEELA,XSPRO,XEPRO'
6076 WRITE (6,
'(2F18.4,6F11.3)')ecmnn(i),elabb(i),xstot(i),xetot(i),
6077 * xsela(i),xeela(i),xspro(i),xepro(i)
6085 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6087 parameter(maxncl = 210)
6088 dimension coop1(3,maxncl)
6089 CALL
conucl(coop1,na,rash)
6097 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6103 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
6104 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
6128 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6131 dimension ec(30),sd(30)
6132 DATA ec /0.d0, 5.d0, 20.d0, 50.d0, 100.d0,
6133 * 200.d0, 500.d0, 1000.d0, 1500.d0, 2000.d0,
6134 * 3000.d0, 4000.d0, 6000.d0, 8000.d0, 10000.d0,
6135 * 15000.d0, 20000.d0, 30000.d0, 40000.d0, 60000.d0,
6136 * 80000.d0, 100000.d0, 150000.d0, 200000.d0, 300000.d0,
6137 * 400000.d0, 600000.d0, 800000.d0, 1000000.d0, 2000000.d0/
6138 DATA sd /0.d0, 0.d0, 5.00d0, 6.14d0, 6.93d0,
6139 * 7.64d0, 8.43d0, 8.87d0, 9.07d0, 9.17d0,
6140 * 9.33d0, 9.40d0, 9.49d0, 9.56d0, 9.58d0,
6141 * 9.69d0, 9.72d0, 9.82d0, 9.85d0, 9.97d0,
6142 * 10.02d0, 10.03d0, 10.13d0, 10.16d0, 10.25d0,
6143 * 10.28d0, 10.39d0, 10.42d0, 10.43d0, 10.53d0/
6146 IF((ecm.GE.ec(i)).AND.(ecm.LT.ec(i+1)))
THEN
6148 del=(ecm-ec(i))*(sd(i+1)-sd(i))/(ec(i+1)-ec(i))
6156 DOUBLE PRECISION FUNCTION siinel(KPROJ,KTARG,UMO)
6157 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6170 DOUBLE PRECISION FUNCTION phnsch ( KP, KTARG, PLAB )
6247 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6248 parameter( kalgnm = 2 )
6249 parameter( anglgb = 5.0
d-16 )
6250 parameter( anglsq = 2.5
d-31 )
6251 parameter( axcssv = 0.2
d+16 )
6252 parameter( andrfl = 1.0
d-38 )
6253 parameter( avrflw = 1.0
d+38 )
6254 parameter( ainfnt = 1.0
d+30 )
6255 parameter( azrzrz = 1.0
d-30 )
6256 parameter( einfnt = +69.07755278982137
d+00 )
6257 parameter( ezrzrz = -69.07755278982137
d+00 )
6258 parameter( onemns = 0.999999999999999
d+00 )
6259 parameter( onepls = 1.000000000000001
d+00 )
6260 parameter( csnnrm = 2.0
d-15 )
6261 parameter( dmxtrn = 1.0
d+08 )
6293 parameter( zerzer = 0.
d+00 )
6294 parameter( oneone = 1.
d+00 )
6295 parameter( twotwo = 2.
d+00 )
6296 parameter( thrthr = 3.
d+00 )
6297 parameter( foufou = 4.
d+00 )
6298 parameter( fivfiv = 5.
d+00 )
6299 parameter( sixsix = 6.
d+00 )
6300 parameter( sevsev = 7.
d+00 )
6301 parameter( eigeig = 8.
d+00 )
6302 parameter( aninen = 9.
d+00 )
6303 parameter( tenten = 10.
d+00 )
6304 parameter( hlfhlf = 0.5
d+00 )
6305 parameter( onethi = oneone / thrthr )
6306 parameter( twothi = twotwo / thrthr )
6307 parameter( pipipi = 3.1415926535897932270
d+00 )
6308 parameter( eneper = 2.7182818284590452354
d+00 )
6309 parameter( sqrent = 1.6487212707001281468
d+00 )
6351 parameter( clight = 2.99792458
d+10 )
6352 parameter( avogad = 6.0221367
d+23 )
6353 parameter( amelgr = 9.1093897
d-28 )
6354 parameter( plckbr = 1.05457266
d-27 )
6355 parameter( elccgs = 4.8032068
d-10 )
6356 parameter( elcmks = 1.60217733
d-19 )
6357 parameter( amugrm = 1.6605402
d-24 )
6358 parameter( ammumu = 0.113428913
d+00 )
6373 parameter( alpfsc = 7.2973530791728595
d-03 )
6374 parameter( fscto2 = 5.3251361962113614
d-05 )
6375 parameter( fscto3 = 3.8859399018437826
d-07 )
6376 parameter( fscto4 = 2.8357075508200407
d-09 )
6377 parameter( plabrc = 0.197327053
d+00 )
6378 parameter( amelct = 0.51099906
d-03 )
6379 parameter( amugev = 0.93149432
d+00 )
6380 parameter( ammuon = 0.105658389
d+00 )
6381 parameter( rclsel = 2.8179409183694872
d-13 )
6382 parameter( gevmev = 1.0
d+03 )
6383 parameter( emvgev = 1.0
d-03 )
6384 parameter( raddeg = 180.
d+00 / pipipi )
6385 parameter( degrad = pipipi / 180.
d+00 )
6415 parameter( mxxrgn = 500 )
6416 parameter( mxxmdf = 56 )
6417 parameter( mxxmde = 50 )
6418 parameter( mfstck = 1000 )
6419 parameter( mestck = 100 )
6420 parameter( nallwp = 39 )
6421 parameter( mpdpdx = 8 )
6422 parameter( icomax = 180 )
6423 parameter( nstbis = 304 )
6424 parameter( idmaxp = 210 )
6453 parameter( lunin = 5 )
6454 parameter( lunout = 6 )
6455 parameter( lunerr = 66 )
6456 parameter( lunber = 14 )
6457 parameter( lunech = 8 )
6458 parameter( lunflu = 86 )
6459 parameter( lungeo = 16 )
6460 parameter( lunpgs = 12 )
6461 parameter( lunran = 2 )
6462 parameter( lunxsc = 81 )
6463 parameter( lunrdb = 1 )
6508 dimension ichrge(39),am(39)
6546 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),
6547 + iich(210),iibar(210),k1(210),k2(210)
6548 dimension kptoip(210),iptokp(39)
6564 DATA kptoip / 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
6565 & 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 66*0,
6566 & 34, 36, 31, 32, 33, 35, 37, 5*0, 38, 5*0, 39, 19*0, 27, 28, 74*0/
6570 DATA iptokp / 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
6571 & 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 99,
6572 & 100, 101, 97, 102, 98, 103, 109, 115 /
6622 COMMON / qquark / iqechr(-6:6), iqbchr(-6:6), iqichr(-6:6),
6623 & iqschr(-6:6), iqcchr(-6:6), iquchr(-6:6),
6624 & iqtchr(-6:6), mquark(3,39)
6626 dimension sieapp(11), sitapp(16), plaetb(16)
6627 dimension sgtcoe(5,33), plalim(2,33), ihlp(nallwp)
6628 dimension sgtco1(5,10),sgtco2(5,8),sgtco3(5,15)
6629 SAVE plaetb, sieapp, sitapp, sgtcoe, plalim, ihlp
6630 SAVE iqfsc1, iqfsc2, iqbsc1, iqbsc2
6631 equivalence(sgtco1(1,1),sgtcoe(1,1))
6632 equivalence(sgtco2(1,1),sgtcoe(1,11))
6633 equivalence(sgtco3(1,1),sgtcoe(1,19))
6635 DATA ihlp/1,4,5*0,1,4,2*0,3,2*2,2*3,1,4,3,3*1,2,
6636 & 2*3, 2, 4*0, 3*4, 1, 4, 1, 4, 1, 4 /
6640 &0.147
d+00, zerzer , zerzer , 0.0022
d+00, -0.0170
d+00,
6642 &0.300
d+00, zerzer , zerzer , 0.0095
d+00, -0.057
d+00,
6644 &16.4
d+00, 19.3
d+00, -0.42
d+00, 0.19
d+00, zerzer ,
6646 &33.0
d+00, 14.0
d+00, -1.36
d+00, 0.456
d+00, -4.03
d+00,
6648 &56.8
d+00, 42.2
d+00, -1.45
d+00, 0.65
d+00, -5.39
d+00,
6650 &18.1
d+00, zerzer , zerzer , 0.26
d+00, -1.0
d+00,
6652 &18.7
d+00, zerzer , zerzer , 0.21
d+00, -0.89
d+00,
6654 &34.2
d+00, 7.9
d+00, -2.1
d+00, 0.346
d+00, -0.99
d+00,
6656 &32.1
d+00, zerzer , zerzer , 0.66
d+00, -5.6
d+00,
6658 &25.2
d+00, zerzer , zerzer , 0.38
d+00, -2.9
d+00/
6662 &57.6
d+00, zerzer , zerzer , 1.17
d+00, -9.5
d+00,
6664 &48.0
d+00, zerzer , zerzer , 0.522
d+00, -4.51
d+00,
6666 &47.30
d+00, zerzer , zerzer , 0.513
d+00, -4.27
d+00,
6668 &91.3
d+00, zerzer , zerzer , 1.05
d+00, -8.8
d+00,
6670 &38.4
d+00, 77.6
d+00, -0.64
d+00, 0.26
d+00, -1.2
d+00,
6672 &zerzer ,133.6
d+00, -0.70
d+00, -1.22
d+00, 13.7
d+00,
6674 &112.
d+00, 125.
d+00, -1.08
d+00, 1.14
d+00, -12.4
d+00,
6676 &30.4
d+00, zerzer , zerzer , zerzer , 1.6
d+00/
6680 &zerzer , 11.4
d+00, -0.4
d+00, 0.079
d+00, zerzer ,
6682 &1.76
d+00, 11.2
d+00, -0.64
d+00, 0.043
d+00, zerzer ,
6684 &5.0
d+00, 8.1
d+00, -1.8
d+00, 0.16
d+00, -1.3
d+00,
6686 &7.3
d+00, zerzer , zerzer , 0.29
d+00, -2.40
d+00,
6688 &11.9
d+00, 26.9
d+00, -1.21
d+00, 0.169
d+00, -1.85
d+00,
6690 &16.1
d+00, zerzer , zerzer , 0.32
d+00, -3.4
d+00,
6692 &10.2
d+00, 52.7
d+00, -1.16
d+00, 0.125
d+00, -1.28
d+00,
6694 &10.6
d+00, 53.1
d+00, -1.19
d+00, 0.136
d+00, -1.41
d+00,
6696 &36.5
d+00, zerzer , zerzer , zerzer , -11.9
d+00,
6698 &12.3
d+00, zerzer , zerzer , zerzer , -2.4
d+00,
6700 &7.24
d+00, 46.0
d+00, -4.71
d+00, 0.279
d+00, -2.35
d+00,
6702 &zerzer ,0.912
d+00, -1.22
d+00, zerzer , zerzer ,
6704 &zerzer , 3.39
d+00, -1.75
d+00, zerzer , zerzer ,
6706 &zerzer , 7.18
d+00, -2.01
d+00, zerzer , zerzer ,
6708 &zerzer , 18.8
d+00, -2.01
d+00, zerzer , zerzer /
6712 & 3.0
d+00, 183.
d+00, 2.0
d+00, 17.8
d+00, 4.0
d+00, 340.
d+00,
6714 & 2.5
d+00, 370.
d+00, 2.5
d+00, 370.
d+00, 2.0
d+00, 310.
d+00,
6716 & 2.0
d+00, 310.
d+00, 2.0
d+00, 310.
d+00, 3.0
d+00, 310.
d+00,
6718 & 1.8
d+00, 310.
d+00, 3.0
d+00, 310.
d+00, 3.0
d+00, 2100.
d+00,
6720 & 3.0
d+00, 370.
d+00, 3.0
d+00, 370.
d+00, 5.0
d+00, 1.73
d+06,
6722 & 1.1
d+00, 280.
d+00, 2.0
d+00, 280.
d+00, 0.6
d+00, 21.
d+00,
6724 & 2.0
d+00, 200.
d+00, 2.0
d+00, 360.
d+00, 2.0
d+00, 175.
d+00,
6726 & 3.0
d+00, 175.
d+00, 3.0
d+00, 2100.
d+00, 2.0
d+00, 384.
d+00,
6728 & 5.0
d+00, 1.73
d+06, 2.0
d+00, 1.59
d+05, 1.1
d+00, 5.55
d+00,
6730 & 0.6
d+00, 24.
d+00, 2.0
d+00, 175.
d+00, 3.5
d+00, 200.
d+00,
6732 & 2.0
d+00, 40.
d+00, 2.0
d+00, 12.8
d+00, 3.0
d+00, 350.
d+00/
6734 DATA plaetb / 0.1
d+00, 0.2
d+00,
6735 & 0.3
d+00, 0.4
d+00, 0.5
d+00, 0.6
d+00, 0.8
d+00, 1.
d+00,
6736 & 1.2
d+00, 1.5
d+00, 2.
d+00, 2.5
d+00, 3.
d+00, 4.
d+00,
6737 & 4.5
d+00, 5.
d+00 /
6740 DATA sieapp / 142.
d+00, 95.1
d+00,
6741 & 75.0
d+00, 70.0
d+00, 62.0
d+00, 57.0
d+00, 48.0
d+00,
6742 & 44.5
d+00, 43.5
d+00, 38.0
d+00, 33.0
d+00 /
6744 DATA sitapp /1129.
d+00, 424.
d+00,
6745 & 239.
d+00, 195.
d+00, 172.
d+00, 150.
d+00, 124.
d+00,
6746 & 117.
d+00, 109.
d+00, 100.
d+00, 90.2
d+00, 81.5
d+00,
6747 & 78.0
d+00, 72.0
d+00, 67.0
d+00, 64.8
d+00 /
6750 ichrge(ktarg)=iich(ktarg)
6751 am(ktarg)=aam(ktarg)
6753 IF ( kp .NE. 26 )
THEN
6769 IF ( iibar(kp) .GT. 0 )
THEN
6775 ELSE IF ( ip .EQ. 15 )
THEN
6781 ELSE IF ( ip .EQ. 24 )
THEN
6787 ELSE IF ( ip .GE. 38 )
THEN
6803 IF ( plab .GT. 50.
d+00 )
THEN
6806 amtsq = am(ktarg)**2
6807 eproj =
sqrt( plab**2 + ampsq )
6808 umosq = ampsq + amtsq + twotwo * am(ktarg) * eproj
6809 eproj =
sqrt( pla**2 + ampsq )
6810 umo50 = ampsq + amtsq + twotwo * am(ktarg) * eproj
6811 umorat =
sqrt( umosq / umo50 )
6815 ELSE IF ( plab .LT. 3.
d+00 )
THEN
6818 amtsq = am(ktarg)**2
6819 eproj =
sqrt( plab**2 + ampsq )
6820 umosq = ampsq + amtsq + twotwo * am(ktarg) * eproj
6821 eproj =
sqrt( pla**2 + ampsq )
6822 umo50 = ampsq + amtsq + twotwo * am(ktarg) * eproj
6823 umorat =
sqrt( umosq / umo50 )
6836 IF ( ihlp(ip) .EQ. 2 )
THEN
6843 sppptt = acof + bcof * pla**enne + ccof * algpla**2
6851 spppel = acof + bcof * pla**enne + ccof * algpla**2
6854 spppin = sppptt - spppel
6861 spmptt = acof + bcof * pla**enne + ccof * algpla**2
6869 spmpel = acof + bcof * pla**enne + ccof * algpla**2
6872 spmpin = spmptt - spmpel
6873 sigdia = spmpin - spppin
6880 IF ( ichrge(ip) .NE. 0 )
THEN
6882 jreac = 3 + ip - 13 + ichrge(ip) * khelp
6883 acof = sgtcoe(1,jreac)
6884 bcof = sgtcoe(2,jreac)
6885 enne = sgtcoe(3,jreac)
6886 ccof = sgtcoe(4,jreac)
6887 dcof = sgtcoe(5,jreac)
6889 shnctt = acof + bcof * pla**enne + ccof * algpla**2
6891 jreac = 19 + ip - 13 + ichrge(ip) * khelp
6892 acof = sgtcoe(1,jreac)
6893 bcof = sgtcoe(2,jreac)
6894 enne = sgtcoe(3,jreac)
6895 ccof = sgtcoe(4,jreac)
6896 dcof = sgtcoe(5,jreac)
6898 shncel = acof + bcof * pla**enne + ccof * algpla**2
6901 shncin = shnctt - shncel
6903 ndiagr = 1 + ip - 13 + ichrge(ip) * khelp
6905 iqfsc1 = 1 + ip - 13
6908 iqbsc2 = 1 + ip - 13
6916 k2hlp = ( kp - 23 ) / 3
6922 ndiagr = 2 + khelp * ( 2 * k2hlp - 1 ) - k2hlp
6923 shncin = hlfhlf * ( spppin + spmpin )
6935 ELSE IF ( ihlp(ip) .EQ. 3 )
THEN
6942 skpptt = acof + bcof * pla**enne + ccof * algpla**2
6950 skppel = acof + bcof * pla**enne + ccof * algpla**2
6953 skppin = skpptt - skppel
6960 skmptt = acof + bcof * pla**enne + ccof * algpla**2
6968 skmpel = acof + bcof * pla**enne + ccof * algpla**2
6971 skmpin = skmptt - skmpel
6972 sigdia = hlfhlf * ( skmpin - skppin )
6975 IF ( ichrge(ip) .NE. 0 )
THEN
6979 IF ( khelp .EQ. 0 )
THEN
6995 shnctt = acof + bcof * pla**enne + ccof * algpla**2
7000 shncin = shnctt - shncel
7019 IF ( khelp .EQ. 0 )
THEN
7020 shncin = skmpin - sigdia
7037 shnctt = acof + bcof * pla**enne + ccof * algpla**2
7042 shncin = shnctt - shncel + sigdia
7059 ELSE IF ( ihlp(ip) .EQ. 4 .AND. ip .LE. 9 )
THEN
7068 sapptt = acof + bcof * pla**enne + ccof * algpla**2
7070 IF ( pla .LT. fivfiv )
THEN
7075 acof = sgtcoe(1,jreac)
7076 bcof = sgtcoe(2,jreac)
7077 enne = sgtcoe(3,jreac)
7078 ccof = sgtcoe(4,jreac)
7079 dcof = sgtcoe(5,jreac)
7081 sappel = acof + bcof * pla**enne + ccof * algpla**2
7084 sappin = sapptt - sappel
7091 spptot = acof + bcof * pla**enne + ccof * algpla**2
7099 sppela = acof + bcof * pla**enne + ccof * algpla**2
7102 sppine = spptot - sppela
7103 sigdia = ( sappin - sppine ) / fivfiv
7107 IF ( ichrge(ip) .NE. 0 )
THEN
7111 IF ( khelp .EQ. 0 )
THEN
7126 shnctt = acof + bcof * pla**enne + ccof * algpla**2
7131 shncin = shnctt - shncel
7139 rnchen =
rndm(rnchen)
7140 IF ( rnchen .LT. puubar )
THEN
7145 iqbsc1 = -iqfsc1 + khelp
7155 IF ( khelp .EQ. 0 )
THEN
7156 shncin = sappin - sigdia
7168 shncin = shnctt - shncel
7176 rnchen =
rndm(rnchen)
7177 IF ( rnchen .LT. pddbar )
THEN
7182 iqbsc1 = -iqfsc1 + khelp - 1
7199 phnsch = ndiagr * sigdia / shncin
7200 iqechc = iqechr(iqfsc1) + iqechr(iqfsc2) + iqechr(iqbsc1)
7202 iqbchc = iqbchr(iqfsc1) + iqbchr(iqfsc2) + iqbchr(iqbsc1)
7206 iqschc = iqschr(iqfsc1) + iqschr(iqfsc2) + iqschr(iqbsc1)
7208 iqspro = iqschr(mquark(1,ip)) + iqschr(mquark(2,ip))
7209 & + iqschr(mquark(3,ip))
7213 WRITE (lunout,*)
' *** Phnsch,kp,ktarg,pla',
7214 &
phnsch,kp,ktarg,pla,
' ****'
7215 WRITE (lunerr,*)
' *** Phnsch,kp,ktarg,pla',
7216 &
phnsch,kp,ktarg,pla,
' ****'
7224 IF ( iqspro .NE. iqschc .OR. ichrge(ip) + ichrge(ktarg)
7225 & .NE. iqechc .OR. iibar(kp) + iibar(ktarg) .NE. iqbchc)
THEN
7227 &
' *** Phnsch,iqspro,iqschc,ichrge,iqechc,ibar,iqbchc,ktarg',
7228 & iqspro,iqschc,ichrge(ip),iqechc,iibar(kp),iqbchc,ktarg
7230 &
' *** Phnsch,iqspro,iqschc,ichrge,iqechc,ibar,iqbchc,ktarg',
7231 & iqspro,iqschc,ichrge(ip),iqechc,iibar(kp),iqbchc,ktarg
7236 IF ( umorat .GT. onepls )
phnsch = oneone / ( ( oneone /
phnsch
7237 & - oneone ) * umorat + oneone )
7240 entry schqua( jqfsc1, jqfsc2, jqbsc1, jqbsc2 )
7295 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7297 COMMON / qquark / iqechr(-6:6), iqbchr(-6:6), iqichr(-6:6),
7298 & iqschr(-6:6), iqcchr(-6:6), iquchr(-6:6),
7299 & iqtchr(-6:6), mquark(3,39)
7303 DATA iqechr / -2, 1, -2, 1, 1, -2, 0, 2, -1, -1, 2, -1, 2 /
7304 DATA iqbchr / 6*-1, 0, 6*1 /
7305 DATA iqichr / 4*0, 1, -1, 0, 1, -1, 4*0 /
7306 DATA iqschr / 3*0, 1, 5*0, -1, 3*0 /
7307 DATA iqcchr / 2*0, -1, 7*0, 1, 2*0 /
7308 DATA iquchr / 0, 1, 9*0, -1, 0 /
7309 DATA iqtchr / -1, 11*0, 1 /
7310 DATA mquark / 1,1,2, -1,-1,-2,
7311 * 0,0,0, 0,0,0, 0,0,0, 0,0,0, 0,0,0,
7312 * 1,2,2, -1,-2,-2, 0,0,0, 0,0,0, 0,0,0,
7313 * 1,-2,0, 2,-1,0, 1,-3,0, 3,-1,0,
7314 * 1,2,3, -1,-2,-3, 0,0,0,
7315 * 2,2,3, 1,1,3, 1,2,3, 1,-1,0,
7316 * 2,-3,0, 3,-2,0, 2,-2,0, 0,0,0,
7317 * 0,0,0, 0,0,0, 0,0,0,
7318 * -1,-1,-3, -1,-2,-3, -2,-2,-3,
7319 * 1,3,3, -1,-3,-3, 2,3,3, -2,-3,-3,
7326 +ptysa2,plaq2,eaq2, amch1,irej,ikvala,pttq1)
7327 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7331 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
7333 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
7335 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
7336 +ipadis,ishmal,lpauli
7338 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
7354 IF (ikvala.EQ.1)b33=6.0
7359 IF (icount.EQ.10)
THEN
7364 IF (icount.GE.1)
THEN
7366 ptxsq1=qtxsq1+hps*cfe
7367 ptysq1=qtysq1+hps*sfe
7368 ptxsa2=qtxsa2-hps*cfe
7369 ptysa2=qtysa2-hps*sfe
7374 es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
7375 hps=
sqrt(es*es+2.*es*0.94)
7378 IF (.NOT.intpt) hps=0.0000001
7382 ptxsq1=qtxsq1+hps*cfe
7383 ptysq1=qtysq1+hps*sfe
7384 ptxsa2=qtxsa2-hps*cfe
7385 ptysa2=qtysa2-hps*sfe
7389 IF (ipev.GE.6)
WRITE(6,1000)ptxsq1,ptysq1,
7391 1000
FORMAT (
' PT S ',8f12.6)
7393 pttq1=ptxsq1**2+ptysq1**2
7394 IF((eq1**2.LE.pttq1)) go to 10
7399 IF (ikvala.EQ.1)b33=6.0
7404 IF (icoun2.EQ.12)
THEN
7411 IF (ipev.GE.6)
WRITE(6,1000)ptxsq1,ptysq1,
7414 ptta2=ptxsa2**2+ptysa2**2
7415 IF((eaq2**2.LE.ptta2)) go to 12
7418 IF(ip.GE.1)go to 1779
7419 plq1=
sqrt(eq1**2-pttq1)
7420 plaq2=-
sqrt(eaq2**2-ptta2)
7425 amch1q=(eq1+eaq2)**2-(ptxsq1+ptxsa2)** 2-(ptysq1+ptysa2)**2-(plq1
7427 IF (amch1q.LE.0.d0)
THEN
7429 301
FORMAT(
' inconsistent Kinematics in SELPT AMCH1Q=',e12.3)
7430 WRITE(6,305) qtxsq1,qtysq1,
7431 +qlq1,qeq1,qtxsa1,qtysa1,qlaq1,qeaq1, qtxsq2,qtysq2,qlq2,qeq2,
7433 +qtysa2,qlaq2,qeaq2, amch1,amch2
7434 305
FORMAT(
'PTXSQ1,PTYSQ1,
7435 +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2,PTYSQ2,PLQ2,EQ2,PTXSA2,
7436 +PTYSA2,PLAQ2,EAQ2, AMCH1,AMCH2',5(4e15.5/))
subroutine calbam(NNCH, I1, I2, IFB11, IFB22, IFB33, IFB44, AMCH, NOBAM, IHAD)
subroutine shmak1(ICASE, NN, NNA, NNB, NA, NB, UMO, BIMP)
subroutine distcm(IOP, NHAD, POLAB, KPROJ, KTARG)
subroutine bamlun(IHAD, KFA1, KFA2, KFA3, KFA4, AEO, IOPT, IREJ)
subroutine modb(BSITE, N, BSTEP, B)
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
subroutine dtwopd(UMO, ECM1, ECM2, PCM1, PCM2, COD1, COF1, SIF1, COD2, COF2, SIF2, AM1, AM2)
double precision function dshnel(KPROJ, KTARG, UMO)
subroutine dtrafo(GAM, BGAM, CX, CY, CZ, COD, COF, SIF, P, ECM, PL, CXL, CYL, CZL, EL)
subroutine sihndi(ECM, KPROJ, KTARG, SIGDIF, SIGDIH)
subroutine indexd(KA, KB, IND)
subroutine ddecay(IHAD, ISTAB)
double precision function sippsd(ECM)
void fill(G4double x, G4double weight=1.)
subroutine shmak(ICASE, NN, NNA, NNB, NA, NB, UMO, BIMP)
double precision function dshnto(KPROJ, KTARG, UMO)
subroutine shmaki(NA, NCA, NB, NCB, RPROJ, RTARG, PPN)
double precision function sampex(X1, X2)
subroutine conucl(X, N, R)
subroutine plot(X, Y, N, M, MM, XO, DX, YO, DY)
double precision function dshpto(IT, PO)
double precision function dbeta(X1, X2, BET)
subroutine py4ent(IP, KF1, KF2, KF3, KF4, PECM, X1, X2, X4, X12, X14)
subroutine py2ent(IP, KF1, KF2, PECM)
subroutine profb(BSTEP, NSTAT, NA, RA, NB, RB, BSITE, NSITEB)
double precision function phnsch(KP, KTARG, PLAB)
subroutine previo(RA, RB, NSTB, BMAX, BSTEP, SIG, RO, G)
double precision function dxlamb(X, Y, Z)
subroutine shmakf(NA, NCA, NB, NCB)
subroutine strafo(GAM, BGAM, CX, CY, CZ, COD, COF, SIF, P, ECM, PL, CXL, CYL, CZL, EL)
virtual G4bool diff(const G4IT &right) const =0
double precision function rndm(RDUMMY)
subroutine pyshow(IP1, IP2, QMAX)
subroutine ddrela(X, Y, Z, COTE, SITE, COPS, SIPS)
subroutine dparje(IHAD, I)
subroutine dbamje(IHAD, KFA1, KFA2, KFA3, KFA4, AE0, IOPT)
double precision function sigsds(S)
subroutine dsfecf(SFE, CFE)
subroutine conuclx(COOP1, NA, RASH, I)
subroutine xsglau(NA, NB, IJPROJ, NTARG)
subroutine title(NA, NB, NCA, NCB)
static c2_log_p< float_type > & log()
make a *new object
G4double tau(G4double bet, G4double homega, G4double ef, G4double t)
subroutine pyjoin(NJOIN, IJOIN)
static c2_sqrt_p< float_type > & sqrt()
make a *new object
subroutine drtran(XO, YO, ZO, CDE, SDE, SFE, CFE, X, Y, Z)
double precision function siinel(KPROJ, KTARG, UMO)
static c2_cos_p< float_type > & cos()
make a *new object
subroutine dthrep(UMO, ECM1, ECM2, ECM3, PCM1, PCM2, PCM3, COD1, COF1, SIF1, COD2, COF2, SIF2, COD3, COF3, SIF3, AM1, AM2, AM3)
subroutine selpts(PTXSQ1, PTYSQ1, PLQ1, EQ1, PTXSA2, PTYSA2, PLAQ2, EAQ2, AMCH1, IREJ, IKVALA, PTTQ1)
subroutine xseapa(ECM, XXXX, IPSQ1, IPSAQ1, XPSQ1, XPSAQ1, IREJ)
subroutine dbklas(I, J, K, I8, I10)
subroutine diagr(NA, NB, B, JS, JT, INT, INTA, INTB)
static c2_sin_p< float_type > & sin()
make a *new object
subroutine shmako(NA, NB, B, INTT, INTA, INTB, JS, JT, PPN, KKMAT)
static c2_exp_p< float_type > & exp()
make a *new object