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(//10
x,39hglauber
s approach is used to simulate ,
1749 +26hnucleus-nucleus collisions/ 24
x,
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
1941 CALL g4dpmjet_modb(
y,b)
1945 IF (b.LT.0.0d0) b =
x1
1946 IF (b.GT.bmax) b = bmax
1955 SUBROUTINE diagr(NA,NB,B,JS,JT,INT,INTA,INTB)
1956 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
1970 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
1972 COMMON /dshm/ rash,rbsh,bmax,bstep,sigsh,rosh,gsh,
1973 * bsite(0:1,200),nstatb,nsiteb
1974 COMMON /dshms/ sigshs
1980 DOUBLE COMPLEX ca,ci
1981 COMMON /damp/ ca,ci,ga
1983 COMMON /kkmatu/kkmato,kkmata,ipoo,ipoa,ipzoo,ipzoa
1984 *,ibproo,ibproa,ireado
1985 COMMON /fluctu/ifluct
1986 COMMON /fluarr/flusi(1000),fluix(1000),fluixx(1000)
1988 dimension js(namx),jt(namx)
1995 CALL
modb(bsite,nsiteb,bstep,b)
2007 IF(intco.GE.500)
THEN
2009 CALL
modb(bsite,nsiteb,bstep,b)
2014 CALL
conucl(pkoo,na,rash)
2016 CALL
conucl(tkoo,nb,rbsh)
2033 IF(ipev.GE.6)
WRITE (6,1000)icnt,pkoo(1,1),tkoo(1,1)
2034 1000
FORMAT (
' 111 FORM IN DIAGR ICNT,PKOO(1,1),TKOO(1,1) ',i6,2f10.3)
2041 ELSEIF(ifluct.EQ.1)
THEN
2042 ifuk=(
rndm(v)+0.001)*1000.
2049 IF(
xy.GT.15.d0) go to 40
2055 IF(
rndm(v).LT.
p) go to 40
2064 IF (js(i).NE.0) inta=inta+1
2067 IF (jt(j).NE.0) intb=intb+1
2071 +
' DIAGR - AFTER 30 CONTINUE: ICNT, INT, B, NA,RA, NB,RB'
2072 WRITE(6,
'(I10,I5,1PE11.3,2(I5,1PE11.3))') icnt,
int, b, na,rash,
2075 IF(
int.EQ.0) go to 30
2080 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2084 COMMON /fluarr/ flusi(1000),fluix(1000),fluixx(1000)
2095 flus=((
x-b)/(om*b))**
n
2097 flusi(i)=(
x/b)*
exp(-((
x-b)/(om*b))**
n)/(
x/b+
a)
2101 flusu=flusu+flusi(i)
2104 flusuu=flusuu+flusi(i)/flusu
2108 3
FORMAT(
' FLUCTUATIONS')
2109 CALL
plot(fluix,flusi,1000,1,1000,0.d0,0.06d0,0.d0,0.01d0)
2113 IF(af.LE.flusi(j))
THEN
2121 fluixx(1000)=fluix(1000)
2141 SUBROUTINE calbam(NNCH,I1,I2,IFB11,IFB22,IFB33,IFB44,
2143 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2189 parameter(nfimax=249)
2190 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
2191 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
2192 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
2195 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
2205 IF(nobam.EQ.4.AND.isymm.EQ.1)
THEN
2207 IF (amch.LT.3.d0)
THEN
2209 IF (rr.LT.0.33333d0)
THEN
2213 ELSEIF (rr.GT.0.666666d0)
THEN
2222 ELSEIF(amch.GT.7.d0)
THEN
2228 ssss=(7.d0-amch)/4.d0
2230 IF(rrr.LT.1.d0-ssss)
THEN
2236 IF (rr.LT.0.33333d0)
THEN
2240 ELSEIF (rr.GT.0.666666d0)
THEN
2251 ELSEIF(nobam.EQ.6.AND.isymm.EQ.1)
THEN
2253 IF (amch.LT.3.d0)
THEN
2255 IF (rr.LT.0.33333d0)
THEN
2259 ELSEIF (rr.GT.0.666666d0)
THEN
2268 ELSEIF(amch.GT.7.d0)
THEN
2274 ssss=(7.d0-amch)/4.d0
2276 IF(rrr.LT.1.d0-ssss)
THEN
2282 IF (rr.LT.0.33333d0)
THEN
2286 ELSEIF (rr.GT.0.666666d0)
THEN
2305 WRITE (6,1000)nnch,i1,i2,ifb1,ifb2,ifb3,ifb4,amch,nobam,ihad
2306 1000
FORMAT(
' CALBAM:NNCH,I1,I2,IFB1,IFB2,IFB3,IFB4,AMCH,NOBAM,IHAD' /7
2310 CALL
dbamje(ihad,ifb1,ifb2,ifb3,ifb4,amch,nobam)
2315 IF(nnch.EQ.-99)
THEN
2330 IF (ifb1.LE.6) go to 40
2332 IF(nobam.EQ.3) CALL
dbamje(ihad,ifb2,ifb1,ifb3,ifb4,amch,nobam)
2335 CALL
dbamje(ihad,ifb1,ifb2,ifb3,ifb4,amch,nobam)
2337 IF(nobam.EQ.6) CALL
dbamje(ihad,ifb3,ifb1,ifb2,ifb4,amch,4)
2339 CALL
dbamje(ihad,ifb3,ifb4,ifb1,ifb2,amch,nobam)
2344 IF (nobam.EQ.3.OR.nobam.EQ.4.OR.nobam.EQ.5)
THEN
2345 CALL
dbamje(ihad,ifb1,ifb2,ifb3,ifb4,amch,nobam)
2346 ELSEIF(nobam.EQ.6)
THEN
2348 CALL
dbamje(ihad,ifb3,ifb1,ifb2,ifb4,amch,4)
2357 IF (iturn.EQ.0) go to 100
2366 WRITE(6,1245)i,pzf(i),hef(i),anf(i)
2367 1245
FORMAT(i5,2f10.4,a8)
2377 SUBROUTINE dbamje(IHAD,KFA1,KFA2,KFA3,KFA4,AE0,IOPT)
2378 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2398 parameter(nfimax=249)
2399 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
2400 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
2401 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
2404 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
2405 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
2408 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210),ich(210), ibar
2409 +(210),k1(210),k2(210)
2410 common/dremai/ rpxr,rpyr,rpzr,rer,kr1r,kr2r
2413 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
2415 COMMON /diffra/ isingd,idiftp,ioudif,iflagd
2417 dimension rpx(100),rpy(100),re(100)
2418 dimension kfr1(100),kfr2(100),iv(100)
2419 parameter(pimass=0.15d0)
2423 IF (lt.EQ.1)
WRITE(6, 1000)ihad,kfa1,kfa2,kfa3,kfa4,ae0,iopt
2424 1000
FORMAT (5i5,e12.4,i5,
2425 +
' BAMJET,IHAD,KFA1,KFA2,KFA3,KFA4,AE0,IOPT')
2432 CALL
bamlun(ihad,kfa1,kfa2,kfa3,kfa4,ae0,iopt,irej)
2438 IF(ifrag.EQ.1.OR.ifrag.EQ.2.OR.ifrag.GE.10)
THEN
2440 CALL
bamlun(ihad,kfa1,kfa2,kfa3,kfa4,ae0,iopt,irej)
2478 IF(itry.GT.10000)
THEN
2479 WRITE(6,
'(/1X,A)')
'DBAMJE:ERROR: FRAGMENTATION IMPOSSIBLE'
2480 WRITE(6,
'(1X,A,5I5,E12.3,I5)')
2481 &
'DBAMJE:IHAD,KFA1,KFA2,KFA3,KFA4,AE0,IOPT ',
2482 & ihad,kfa1,kfa2,kfa3,kfa4,ae0,iopt
2496 IF(itry.GT.100)
THEN
2498 IF(kfa1.GT.(kfa2-6))
THEN
2499 e0=ae0-max(ae0*0.1,pimass)
2501 e0=max(ae0*0.1,pimass)
2506 IF(iopt.EQ.1.OR.iopt.EQ.2) e0=ae0
2509 IF(kfa1.GT.6.AND.iopt.EQ.1) ll=1
2510 IF(kfa1.LE.6.AND.iopt.EQ.2) ll=1
2511 IF(kfa1.GT.6.AND.iopt.EQ.4) ll=1
2527 IF(iopt.EQ.3.AND.ll.EQ.0) goto 50
2528 IF(iopt.EQ.4.AND.kfa1.LE.6.AND.ll.EQ.0) goto 50
2529 IF(iopt.EQ.4.AND.kfa1.GT.6.AND.ll.EQ.1) goto 50
2530 IF(iopt.EQ.5.AND.ll.EQ.0) goto 50
2540 IF(iopt.EQ.4.AND.kfa1.GT.6) ll=0
2554 IF(amf(it).GT.rx) goto 10
2555 IF(amf(it).LE.rx) goto 70
2571 90 rpx(it)=rpx(j)-hpx
2574 IF (iopt.EQ.1.AND.ll.EQ.1)hpz=-hpz
2575 IF(iopt.EQ.2.AND.ll.EQ.1) hpz=-hpz
2576 IF(iopt.EQ.4.AND.kfa1.GT.6) hpz=-hpz
2577 IF(iopt.EQ.5) hpz=-hpz
2585 WRITE(6, 1010)pgx,pgy,pgz
2586 1010
FORMAT(1h0,12hpgx,pgy,pgz=,3f8.4)
2591 IF(iopt.EQ.1.OR.iopt.EQ.2) goto 150
2599 IF(iopt.EQ.3.OR.iopt.EQ.4.OR.iopt.EQ.5) goto 160
2600 IF(ll.EQ.0) goto 160
2608 IF(le.EQ.0) goto 180
2611 WRITE(6, 1020)nref(j),anf(j),amf(j),ichf(j),
2612 + ibarf(j),pxf(j),pyf(j),
2614 1020
FORMAT(2
x,i3,a6,f6.3,2i4,4f8.4)
2615 1030
FORMAT(2
x,
'NF,NAME,MASS,IQ,IB,PX,PY,PZ,E')
2618 1040
FORMAT(1h0,38hnumber of events with prest gt. erest=,i4, /,
2619 +21hnumber of all events=,i4)
2620 1050
FORMAT(1h0,
'NUMBER OF EVENTS WITH ONLY ONE PARTICLE=',i4)
2625 1060
FORMAT(1h0,
' MULTIPLIZITAET=',i3)
2626 1070
FORMAT(1h0,13hhadronanzahl=,i3)
2641 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2648 IF(kp.EQ.4.AND.ks.EQ.5)ind=4
2650 IF(kp.EQ.6.AND.ks.EQ.7)ind=6
2651 IF(kp.EQ.4.AND.ks.EQ.4)ind=7
2652 IF(kp.EQ.6.AND.ks.EQ.5)ind=8
2655 IF(kp.EQ.12.AND.ks.EQ.8)ind=11
2657 IF(kp.EQ.12.AND.ks.EQ.7)ind=13
2675 DOUBLE PRECISION FUNCTION dbeta(X1,X2,BET)
2676 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2680 IF(betx1.LT.70.) ax=-1./bet**2*(betx1+1.)*
exp(-betx1)
2681 ay=1./bet**2*(bet*
x2+1.)*
exp(-bet*
x2)
2690 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2692 x1=cops*
x-sips*cote*
y+sips*site*
z
2693 x2=sips*
x+cops*cote*
y-cops*site*
z
2751 SUBROUTINE dthrep(UMO,ECM1,ECM2,ECM3,PCM1,PCM2,PCM3,COD1,COF1,
2752 &sif1,cod2,cof2,sif2,cod3,cof3,sif3,am1,am2,am3)
2799 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
2801 parameter( kalgnm = 2 )
2802 parameter( anglgb = 5.0
d-16 )
2803 parameter( anglsq = 2.5
d-31 )
2804 parameter( axcssv = 0.2
d+16 )
2805 parameter( andrfl = 1.0
d-38 )
2806 parameter( avrflw = 1.0
d+38 )
2807 parameter( ainfnt = 1.0
d+30 )
2808 parameter( azrzrz = 1.0
d-30 )
2809 parameter( einfnt = +69.07755278982137
d+00 )
2810 parameter( ezrzrz = -69.07755278982137
d+00 )
2811 parameter( onemns = 0.999999999999999
d+00 )
2812 parameter( onepls = 1.000000000000001
d+00 )
2813 parameter( csnnrm = 2.0
d-15 )
2814 parameter( dmxtrn = 1.0
d+08 )
2846 parameter( zerzer = 0.
d+00 )
2847 parameter( oneone = 1.
d+00 )
2848 parameter( twotwo = 2.
d+00 )
2849 parameter( thrthr = 3.
d+00 )
2850 parameter( foufou = 4.
d+00 )
2851 parameter( fivfiv = 5.
d+00 )
2852 parameter( sixsix = 6.
d+00 )
2853 parameter( sevsev = 7.
d+00 )
2854 parameter( eigeig = 8.
d+00 )
2855 parameter( aninen = 9.
d+00 )
2856 parameter( tenten = 10.
d+00 )
2857 parameter( hlfhlf = 0.5
d+00 )
2858 parameter( onethi = oneone / thrthr )
2859 parameter( twothi = twotwo / thrthr )
2860 parameter( pipipi = 3.1415926535897932270
d+00 )
2861 parameter( eneper = 2.7182818284590452354
d+00 )
2862 parameter( sqrent = 1.6487212707001281468
d+00 )
2904 parameter( clight = 2.99792458
d+10 )
2905 parameter( avogad = 6.0221367
d+23 )
2906 parameter( amelgr = 9.1093897
d-28 )
2907 parameter( plckbr = 1.05457266
d-27 )
2908 parameter( elccgs = 4.8032068
d-10 )
2909 parameter( elcmks = 1.60217733
d-19 )
2910 parameter( amugrm = 1.6605402
d-24 )
2911 parameter( ammumu = 0.113428913
d+00 )
2926 parameter( alpfsc = 7.2973530791728595
d-03 )
2927 parameter( fscto2 = 5.3251361962113614
d-05 )
2928 parameter( fscto3 = 3.8859399018437826
d-07 )
2929 parameter( fscto4 = 2.8357075508200407
d-09 )
2930 parameter( plabrc = 0.197327053
d+00 )
2931 parameter( amelct = 0.51099906
d-03 )
2932 parameter( amugev = 0.93149432
d+00 )
2933 parameter( ammuon = 0.105658389
d+00 )
2934 parameter( rclsel = 2.8179409183694872
d-13 )
2935 parameter( gevmev = 1.0
d+03 )
2936 parameter( emvgev = 1.0
d-03 )
2937 parameter( raddeg = 180.
d+00 / pipipi )
2938 parameter( degrad = pipipi / 180.
d+00 )
2965 parameter( lunin = 5 )
2966 parameter( lunout = 6 )
2967 parameter( lunerr = 66 )
2968 parameter( lunber = 14 )
2969 parameter( lunech = 8 )
2970 parameter( lunflu = 86 )
2971 parameter( lungeo = 16 )
2972 parameter( lunpgs = 12 )
2973 parameter( lunran = 2 )
2974 parameter( lunxsc = 81 )
2975 parameter( lunrdb = 1 )
3004 parameter( mxxrgn = 500 )
3005 parameter( mxxmdf = 56 )
3006 parameter( mxxmde = 50 )
3007 parameter( mfstck = 1000 )
3008 parameter( mestck = 100 )
3009 parameter( nallwp = 39 )
3010 parameter( mpdpdx = 8 )
3011 parameter( icomax = 180 )
3012 parameter( nstbis = 304 )
3013 parameter( idmaxp = 210 )
3019 dimension
f(5),
xx(5)
3021 COMMON /dgamre/ redu,amo,amm(15 )
3022 common/ddrei/uumo,aam1,aam2,aam3,s22,umo2,
3023 *am11,am22,am33,s2sup,s2sap(2)
3062 IF(rho2.LT.rho1) go to 125
3064 125 s2sup=(s22-s21)*.5d0+s21
3065 suprho=
dxlamb(s2sup,umo2,am11)*
dxlamb(s2sup,am22,am33)/
3067 suprho=suprho*1.05d0
3069 IF (gu.LT.go.AND.xo.LT.gu) xo=gu
3070 IF (gu.GT.go.AND.xo.GT.gu) xo=gu
3079 x4=(
xx(1)+
xx(2))*0.5d0
3080 x5=(
xx(2)+
xx(3))*0.5d0
3090 IF (
f(ii).GE.
f(iii)) go to 128
3103 IF (
xx(ii).GE.
xx(iii)) go to 129
3117 IF (ith.GT.200) redu=-9.d0
3118 IF (ith.GT.200) go to 400
3121 s2=am23+c*(umo-am1-am2-am3)*(umo-am1+am2+am3)
3125 IF(
y.GT.
rho) go to 1
3128 s1=s1*
rho+am11+am22-(s2-umo2+am11)*(s2+am22-am33)/(2.d0*s2)-
3130 s3=umo2+am11+am22+am33-s1-s2
3131 ecm1=(umo2+am11-s2)/umoo
3132 ecm2=(umo2+am22-s3)/umoo
3133 ecm3=(umo2+am33-s1)/umoo
3134 pcm1=
sqrt((ecm1+am1)*(ecm1-am1))
3135 pcm2=
sqrt((ecm2+am2)*(ecm2-am2))
3136 pcm3=
sqrt((ecm3+am3)*(ecm3-am3))
3141 IF ( pcm12 .LT. anglsq ) go to 200
3142 costh=(ecm1*ecm2+0.5
d+00*(am11+am22-s1))/pcm12
3146 costh=(uw-0.5
d+00)*2.
d+00
3150 IF(abs(costh).GT.oneone)
3151 &costh=sign(oneone,costh)
3152 IF (redu.LT.1.
d+00)
RETURN
3153 costh2=(pcm3*pcm3+pcm2*pcm2-pcm1*pcm1)/(2.
d+00*pcm2*pcm3)
3156 IF(abs(costh2).GT.oneone)
3157 &costh2=sign(oneone,costh2)
3158 sinth2=
sqrt((oneone-costh2)*(oneone+costh2))
3159 sinth =
sqrt((oneone-costh)*(oneone+costh))
3160 sinth1=costh2*sinth-costh*sinth2
3161 costh1=costh*costh2+sinth2*sinth
3173 cod3=twotwo*
rndm(cod3)-oneone
3174 sid3=
sqrt((1.
d+00-cod3)*(1.
d+00+cod3))
3176 cod1=cx11*cod3+cz11*sid3
3177 chlp=(oneone-cod1)*(oneone+cod1)
3178 IF(chlp.LT.1.
d-14)
WRITE(isys,2)cod1,cof3,sid3,
3181 cof1=(cx11*sid3*cof3-cy11*sif3-cz11*cod3*cof3)/sid1
3182 sif1=(cx11*sid3*sif3+cy11*cof3-cz11*cod3*sif3)/sid1
3183 cod2=cx22*cod3+cz22*sid3
3184 sid2=
sqrt((oneone-cod2)*(oneone+cod2))
3185 cof2=(cx22*sid3*cof3-cy22*sif3-cz22*cod3*cof3)/sid2
3186 sif2=(cx22*sid3*sif3+cy22*cof3-cz22*cod3*sif3)/sid2
3189 eochck = umo - ecm1 - ecm2 - ecm3
3193 pzchck = pcm1 * cod1 + pcm2 * cod2 + pcm3 * cod3
3194 pxchck = pcm1 * cof1 * sid1 + pcm2 * cof2 * sid2
3195 & + pcm3 * cof3 * sid3
3196 pychck = pcm1 * sif1 * sid1 + pcm2 * sif2 * sid2
3197 & + pcm3 * sif3 * sid3
3198 eocmpr = 1.
d-12 * umo
3199 IF ( abs(eochck) + abs(pxchck) + abs(pychck) + abs(pzchck)
3200 & .GT. eocmpr )
THEN
3202 &
' *** Threpd: energy/momentum conservation failure! ***',
3203 & eochck,pxchck,pychck,pzchck
3204 WRITE(lunerr,*)
' *** SID1,SID2,SID3',sid1,sid2,sid3
3211 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3214 COMMON /dgamre/ redu,amo,amm(15 )
3215 common/ddrei/
test(12)
3219 IF (idgb.LE.0) go to 11
3224 12
FORMAT(/,10
x,
' DXLAMB PRINT')
3225 13
FORMAT(10
x,60(1h*))
3226 10
FORMAT(4e20.8,
'DXLAMB',/,12f10.5)
3247 SUBROUTINE strafo(GAM,BGAM,CX,CY,CZ,COD,COF,SIF,P,ECM,
3249 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3252 sid=
sqrt((1.-cod)*(1.+cod)+1.
e-22)
3253 sif=
sqrt((1.-cof)*(1.+cof)+1.
e-22)
3257 plz=gam*pcmz+bgam*ecm
3258 pl=
sqrt(plx*plx+ply*ply+plz*plz)
3259 el=gam*ecm+bgam*pcmz
3262 IF(coz.GE.1.)coz=0.999999999999
3263 siz=
sqrt((1.-coz)*(1.+coz))
3264 CALL
drtran(cx,cy,cz,coz,siz,sif,cof,cxl,cyl,czl)
3268 SUBROUTINE drtran(XO,YO,ZO,CDE,SDE,SFE,CFE,X,Y,Z)
3269 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3271 IF (abs(xo)-0.0001) 10,10,30
3272 10
IF (abs(yo)-0.0001) 20,20,30
3283 x=-yo*xi/
a-zo*xo*yi/
a+xo*zi
3284 y=xo*xi/
a-zo*yo*yi/
a+yo*zi
3293 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3298 parameter(idmax9=602)
3300 common/ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
3304 common/dpar/aname(210),am(210),ga(210),
tau(210),ich(210),ibar(210)
3329 SUBROUTINE dtwopd(UMO,ECM1,ECM2,PCM1,PCM2,COD1,COF1,SIF1, COD2,
3331 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3335 IF(umo.LT.(am1+am2))
THEN
3336 WRITE(6,
'(/,A/A,3(1PE12.4))')
3337 +
' INCONSISTENT CALL OF TWOPAD / EXECUTION STOPPED',
3338 +
' UMO, AM1, AM2 :', umo, am1, am2
3342 ecm1=((umo-am2)*(umo+am2) + am1*am1)/(2.*umo)
3344 pcm1=
sqrt((ecm1-am1)*(ecm1+am1))
3358 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3368 IF(
rndm(v).LT.0.5d0) goto20
3378 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3380 CHARACTER*8 zkname,
z
3383 parameter(idmax9=602)
3385 common/ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
3389 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210), ich(210),ibar
3390 +(210),k1(210),k2(210)
3392 dimension ichar(210)
3393 equivalence(ich(1),ichar(1))
3397 1000
FORMAT(1
h1,
' ')
3399 1010
FORMAT(///
' TABLE OF USED PARTICLES AND RESONANCES (I)',//
3400 +
' I = NUMBER OF PARTICLE OR RESONANCE',/
3401 +
' IPDG = P D G NUMBER OF PARTICLE OR RESONANCE',/
3402 +
' ANAME = NAME OF I'/,
' AM = MASS OF I (GEV)',/
3403 +
' GA = WIDTH OF I (GEV)',/
' TAU = LIFE TIME OF I (SEC.)',/
3404 +
' ICH = ELECTRIC CHARGE OF I, IBAR = BARYONIC CHARGE OF I',/
' ',
'
3405 +K1 = FIRST DECAY CHANNEL NUMBER, K2 = LAST DECAY CHANNEL NUMBER OF
3411 +
' I ANAME AM GA TAU ICH IBAR K1 K2'/)
3415 WRITE(6, 1030)i,ipdg,aname(i),am(i),
3416 + ga(i),
tau(i),ich(i),ibar(i), k1
3418 1030
FORMAT (1i4,i6,2
x,1a8,3e11.4,4i4)
3419 IF(i.EQ.43)
WRITE(6, 1000)
3420 IF(i.EQ.43)
WRITE(6, 1020)
3421 IF(i.EQ.99)
WRITE(6, 1000)
3422 IF(i.EQ.99)
WRITE(6, 1020)
3423 IF(i.EQ.155)
WRITE(6, 1000)
3424 IF(i.EQ.155)
WRITE(6, 1020)
3428 1040
FORMAT(///
' DECAY CHANNELS OF PARTICLES AND RESONANCES',//)
3430 1050
FORMAT(
' ANAME = PARTICLE AND RESONANCE NAME'/,
3431 +
' DNAME = DECAY CHANNEL NAME'/,
' J = DECAY CHANNEL NUMBER'/,
3432 +
' I = NUMBER OF DECAYING PARTICLE'/,
3433 +
' WT = SUM OF DECAY CHANNEL WEIGHTS FROM K1(I) UP TO J'/,
3434 +
' NZK = PROGRAM INTERNAL NUMBERS OF DECAY PRODUCTS')
3437 1060
FORMAT(///
' I J ANAME DNAME DECAY
3442 IF (ik1.LE.0) go to 60
3456 WRITE(6, 1070)i,ik,aname(i),zkname(ik),(
z(j),j=1,3),wt(ik),j1,j2,
3458 1070
FORMAT(2i5,
' DECAY OF ',1a8,
' (CHANNEL: ',1a6,
' ) TO ',3(1a6,2
x),
3460 amtest=am(i)-am(j1)-am(j2)-am(j3)
3461 ibtest=ibar(i)-ibar(j1)-ibar(j2)-ibar(j3)
3462 ictest=ichar(i)-ichar(j1)-ichar(j2)-ichar(j3)
3463 IF (amtest) 20,30,30
3468 IF (mtest+ibtest**2+ictest**2.NE.0)
WRITE(6, 1080)amtest,
3471 1080
FORMAT (
' ***** ERROR IN MASS, BAR.CH. OR E.CH. ',f10.5,2i6)
3472 IF(ik.EQ.27)
WRITE(6, 1000)
3473 IF(ik.EQ.27)
WRITE(6, 1060)
3474 IF(ik.EQ.62)
WRITE(6, 1000)
3475 IF(ik.EQ.62)
WRITE(6, 1060)
3476 IF(ik.EQ.101)
WRITE(6, 1000)
3477 IF(ik.EQ.101)
WRITE(6, 1060)
3478 IF(ik.EQ.144)
WRITE(6, 1000)
3479 IF(ik.EQ.144)
WRITE(6, 1060)
3480 IF(ik.EQ.183)
WRITE(6, 1000)
3481 IF(ik.EQ.183)
WRITE(6, 1060)
3482 IF(ik.EQ.222)
WRITE(6, 1000)
3483 IF(ik.EQ.222)
WRITE(6, 1060)
3484 IF(ik.EQ.261)
WRITE(6, 1000)
3485 IF(ik.EQ.261)
WRITE(6, 1060)
3486 IF(ik.EQ.300)
WRITE(6, 1000)
3487 IF(ik.EQ.300)
WRITE(6, 1060)
3488 IF(ik.EQ.362)
WRITE(6, 1000)
3489 IF(ik.EQ.362)
WRITE(6, 1060)
3490 IF(ik.EQ.401)
WRITE(6, 1000)
3491 IF(ik.EQ.401)
WRITE(6, 1060)
3492 IF(ik.EQ.440)
WRITE(6, 1000)
3493 IF(ik.EQ.440)
WRITE(6, 1060)
3494 IF(ik.EQ.479)
WRITE(6, 1000)
3495 IF(ik.EQ.479)
WRITE(6, 1060)
3496 IF(ik.EQ.518)
WRITE(6, 1000)
3497 IF(ik.EQ.518)
WRITE(6, 1060)
3507 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3509 common/dinpda/imps(6,6),imve(6,6),ib08(6,21),ib10(6,21),
3510 *ia08(6,21),ia10(6,21),a1,b1,b2,b3,lt,lb,bet,as,b8,ame,diq,isu
3511 dimension iv(36),ip(36),ib(126),ibb(126),ia(126),iaa(126)
3518 *23,14,16,116,0,0,13,23,25,117,0,0,15,24,31,120,0,0,119,118,121,
3528 *33,34,38,123,0,0,32,33,39,124,0,0,36,37,96,127,0,0,126,125,128,
3538 *0,1,21,140,0,0,8,22,137,0,0,97,138,0,0,146,5*0,
3539 *1,8,22,137,0,0,0,20,142,0,0,98,139,0,0,147,5*0,
3540 *21,22,97,138,0,0,20,98,139,0,0,0,145,0,0,148,5*0,
3541 *140,137,138,146,0,0,142,139,147,0,0,145,148,50*0/
3550 *53,54,104,161,0,0,55,105,162,0,0,107,164,0,0,167,5*0,
3551 *54,55,105,162,0,0,56,106,163,0,0,108,165,0,0,168,5*0,
3552 *104,105,107,164,0,0,106,108,165,0,0,109,166,0,0,169,5*0,
3553 *161,162,164,167,0,0,163,165,168,0,0,166,169,0,0,170,47*0/
3562 *0,2,99,152,0,0,9,100,149,0,0,102,150,0,0,158,5*0,
3563 *2,9,100,149,0,0,0,101,154,0,0,103,151,0,0,159,5*0,
3564 *99,100,102,150,0,0,101,103,151,0,0,0,157,0,0,160,5*0,
3565 *152,149,150,158,0,0,154,151,159,0,0,157,160,50*0/
3574 *67,68,110,171,0,0,69,111,172,0,0,113,174,0,0,177,5*0,
3575 *68,69,111,172,0,0,70,112,173,0,0,114,175,0,0,178,5*0,
3576 *110,111,113,174,0,0,112,114,175,0,0,115,176,0,0,179,5*0,
3577 *171,172,174,177,0,0,173,175,178,0,0,176,179,0,0,180,47*0/
3639 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3659 parameter(nfimax=249)
3660 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
3661 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
3662 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
3665 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
3666 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
3668 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
3672 COMMON /dkpl/uplo,ipq
3673 COMMON /dinv/pnuc(3),inucvt
3674 common/capkop/xx1,xx3
3677 1000
FORMAT (
' ##############################################'/
3678 +
' PROGRAM TECABAPT'/
3679 +
' ######################################################')
3689 READ(5,1010)jni,nevt,ip,ncases,poo,aoo,znuc
3690 1010
FORMAT(4i10,3f10.2)
3692 IF (jni.LE.0) go to 120
3693 WRITE(6, 1010)jni,nevt,ip,ncases,poo,aoo,znuc,it
3697 go to(20,30,40,50,60,70,80,100),jni
3716 CALL
distcm(1,ipq,poo,ipq,ipq)
3719 IF (ip.EQ.103)CALL
calbam(0,1,1,1,7,1,1,poo,3,nhad)
3720 IF (ip.EQ.109)CALL
calbam(0,1,1,7,1,1,1,poo,3,nhad)
3721 IF (ip.EQ.104)CALL
calbam(0,1,1,1,2,2,1,poo,4,nhad)
3722 IF (ip.EQ.1010)CALL
calbam(0,1,1,1,2,3,1,poo,4,nhad)
3723 IF (ip.EQ.105)CALL
calbam(0,1,1,1,2,7,8,poo,5,nhad)
3724 IF (ip.EQ.1011)CALL
calbam(0,1,1,7,7,1,1,poo,5,nhad)
3725 IF (ip.EQ.106)CALL
calbam(0,1,1,1,1,1,1,poo,6,nhad)
3726 IF (ip.EQ.1012)CALL
calbam(0,1,1,7,7,7,1,poo,6,nhad)
3727 IF (ip.EQ.1050)CALL
calbam(0,1,1,1,1,2,1,poo,10,nhad)
3730 CALL
distcm(2,nhad,poo,ipq,ncases)
3733 WRITE(6, 1020)poo,ip,ncases
3734 1020
FORMAT (
' BAMJET (POO,IP,NCASES) = ',1f10.2,2i10)
3735 CALL
distcm(3,nevt,poo,ipq,ncases)
3748 SUBROUTINE distcm(IOP,NHAD,POLAB,KPROJ,KTARG)
3749 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
3756 parameter(nfimax=249)
3757 COMMON /dfinpa/ anh(nfimax),
px(nfimax),
py(nfimax),
pz(nfimax),
3758 +he(nfimax),am(nfimax), ich(nfimax),ibar(nfimax),nr(nfimax)
3759 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
3763 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),iich(210),
3764 +iibar(210),k1(210),k2(210)
3766 COMMON /histo / xmult(100,10),ymult(100,10),xxfl(50,20), yxfl
3767 +(50,20),xyl(50,20),yyl(50,20), yylps(50,20),ptp(50,20),pty(50,20),
3769 dimension avmult(12,30),ave(12,30),indx(25),mu(12,30), akno
3770 +(100,2),xkno(100,2),ake(12,30),aaso(12,30)
3771 COMMON /dkpl/uplo,kpl
3775 DATA indx/1,8,10,10,10,10,7,2,7,10,10,7,3,4,5,6,
3776 *11,12,7,13,14,15,16,17,18/
3778 go to(10,60,100),iop
3787 IF(jni.EQ.7) umo=polab
3791 WRITE(6, 1000)eeo,po,nhad
3792 1000
FORMAT (
' EEO',f10.2,f10.2,i10)
3796 avmult(kpl,i)=1.
e-18
3804 WRITE(6, 1000)eeo,po,nhad
3807 xxfl(j,i)=j*dxfl -1.
3815 WRITE(6, 1000)eeo,po,nhad
3819 avmult(kpl,i)=1.
d-18
3823 WRITE(6, 1000)eeo,po,nhad
3828 avmult(kpl,30)=avmult(kpl,30)+nhad
3830 IF (nnhad.GT.100) nnhad=100
3831 ymult(nnhad,10)=ymult(nnhad,10)+1.
3837 IF(ibar(i).NE.500)
THEN
3842 IF(eetot.GT.polab+1.
d-6)
THEN
3843 WRITE(6,*).gt.
' eetotpolab ',eetot,polab
3847 IF(ibar(i).NE.500)
THEN
3849 IF (nre.GT.25) nre=3
3850 IF (nre.LT. 1) nre=3
3853 ave(kpl,nre)=ave(kpl,nre)+he(i)
3854 ave(kpl,30)=ave(kpl,30)+he(i)
3855 IF (ni.NE.6) ave(kpl,29)=ave(kpl,29)+he(i)
3856 avmult(kpl,nre)=avmult(kpl,nre)+1.
3857 IF (ni.NE.6) avmult(kpl,29)=avmult(kpl,29)+1.
3858 mu(kpl,ni)=mu(kpl,ni)+1
3859 IF (ich(i).NE.0)mu(kpl,9)=mu(kpl,9)+1
3862 IF (ich(i).NE.0)ave(kpl,27)=ave(kpl,27)+he(i)
3863 IF (ich(i).NE.0)avmult(kpl,27)=avmult(kpl,27)+1
3865 xfl=(
pz(i)/abs(
pz(i)))*he(i)/po
3867 IF (ixfl.LT. 1) ixfl=1
3868 IF (ixfl.GT.50) ixfl=50
3871 IF (ich(i).NE.0)yxfl(ixfl,9)=yxfl(ixfl,9)+xxxfl
3872 yxfl(ixfl,ni)=yxfl(ixfl,ni)+xxxfl
3873 yxfl(ixfl,10)=yxfl(ixfl,10)+xxxfl
3874 ptt=
px(i)**2+
py(i)**2
3875 yl=0.5*
log(abs((he(i)+
pz(i)+1.
e-10)/(he(i)-
pz(i)+1.
e-10)))
3878 IF (iylps.LT.1)iylps=1
3879 IF (iylps.GT.50)iylps=50
3880 yylps(iylps,ni)=yylps(iylps,ni)+1.
3881 yylps(iylps,10)=yylps(iylps,10)+1.
3882 IF (ich(i).NE.0)yylps(iylps,9)=yylps(iylps,9)+1.
3885 IF (iyl.GT.50) iyl=50
3886 IF (ich(i).NE.0)yyl(iyl,9)=yyl(iyl,9)+1.
3887 yyl(iyl,ni)=yyl(iyl,ni)+1.
3888 yyl(iyl,10)=yyl(iyl,10)+1.
3894 IF (ipt.GT.50) ipt=50
3895 IF (ich(i).NE.0)pty(ipt,9)=pty(ipt,9)+1./
pt
3896 pty(ipt,ni)=pty(ipt,ni)+1./
pt
3897 pty(ipt,10)=pty(ipt,10)+1./
pt
3902 IF (im.GT.100)im=100
3903 ymult(im,i)=ymult(im,i)+1.
3908 WRITE(6, 1000)eeo,po,nhad
3911 avmult(kpl,i)=avmult(kpl,i)/nhad
3912 ave(kpl,i)=ave(kpl,i)/nhad
3915 WRITE (6,1030)avpt,navpt
3916 1030
FORMAT (
' AVERAGE PT= ',f12.4,i10)
3918 1040
FORMAT(
' PARTICLE REF,CHAR,IBAR, MASS AVERAGE',
3919 +
' ENERGY, MULTIPLICITY, INELASTICITY')
3921 ake(kpl,i)=ave(kpl,i)/eeo
3922 WRITE(6, 1050)aname(i),i,iich(i),iibar(i),
3923 + aam(i), ave(kpl,i),avmult
3924 + (kpl,i),ake(kpl,i)
3925 1050
FORMAT (
' ',a8,3i5,f10.3,3f18.6)
3929 ymult(j,i)=ymult(j,i)/nhad
3933 yxfl(j,i)=yxfl(j,i)/(nhad*dxfl)
3934 yy l(j,i)=
yy l(j,i)/(nhad*
dy)
3935 yylps(j,i)=yylps(j,i)/(nhad*
dy)
3936 pty(j,i)=pty(j,i)/(nhad*dpt)
3940 1060
FORMAT(
'1 RAPIDITY DISTRIBUTION')
3942 WRITE(6, 1070)xyl(j,1),(yyl(j,i),i=1,10)
3943 1070
FORMAT (f10.2,10e11.3)
3946 WRITE(6, 1070)xyl(j,1),(yyl(j,i),i=11,20)
3949 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3950 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3951 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
3952 CALL
plot(xyl,yyl,1000,20,50,-5.d0,
dy,0.d0,0.1d0)
3953 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3954 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3955 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
3956 CALL
plot(xyl,yylps,1000,20,50,-5.d0,
dy,0.d0,0.1d0)
3959 1080
FORMAT (
'1 LONG MOMENTUM (SCALED) DISTRIBUTION')
3961 WRITE(6, 1070)xxfl(j,1),(yxfl(j,i),i=1,10)
3964 WRITE(6, 1070)xxfl(j,1),(yxfl(j,i),i=11,20)
3968 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3969 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3970 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
3971 CALL
plot(xxfl,yxfl,1000,20,50,-1.d0,dxfl,0.d0,0.05d0)
3973 1090
FORMAT (
'1 MULTIPLICITY DISTRIBUTIONS')
3977 sumul=sumul+ymult(j,10)
3978 simul=simul+ymult(j,9)
3980 WRITE(6, 1100)(xmult(j,1),ymult(j,9),ymult(j,10),j=1,100)
3981 1100
FORMAT(f6.1,2e12.4,f6.1,2e12.4,f6.1,2e12.4,f6.1,2e12.4)
3983 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3984 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3985 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
3986 CALL
plot(xmult,ymult,1000,10,100,0.d0,1.d0,0.d0,0.01d0)
3988 xkno(i,1)=i/avmult(kpl,30)
3989 xkno(i,2)=i/avmult(kpl,27)
3990 akno(i,1)=ymult(i,10)*avmult(kpl,30)/sumul
3991 akno(i,2)=ymult(i,9)*avmult(kpl,27)/simul
3992 akno(i,1)=log10(akno(i,1)+1.
d-9)
3993 akno(i,2)=log10(akno(i,2)+1.
d-9)
3996 1110
FORMAT (
'1 KNO MULTIPLICITY DISTRIBUTIONS')
3997 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
3998 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
3999 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4000 CALL
plot(xkno,akno,200,2,100,0.d0,0.08d0,-4.d0,0.05d0)
4003 ymult(j,i)=log10(ymult(j,i))
4007 yxfl(j,i)=log10(abs(yxfl(j,i)+1.
d-8))
4008 yyl(j,i)=log10(yyl(j,i)+1.
d-8)
4009 pty(j,i)=log10(pty(j,i)+1.
d-8)
4013 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4014 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4015 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4016 CALL
plot(xyl,yyl,1000,20,50,-5.d0,
dy,-3.5d0,0.05d0)
4018 WRITE(6, 1070)xxfl(j,1),(yxfl(j,i),i=1,10)
4021 WRITE(6, 1070)xxfl(j,1),(yxfl(j,i),i=11,20)
4024 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4025 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4026 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4027 CALL
plot(xxfl,yxfl,1000,20,50,-1.d0,dxfl,-4.5d0,0.05d0)
4029 1120
FORMAT (
'1 PT DISTRIBUTION DN/PTDPT')
4030 CALL
plot(ptp,pty,1000,20,50,0.d0,dpt,-2.0d0,0.05d0)
4031 IF (ipriop.EQ.1) go to 250
4033 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4034 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4035 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4036 CALL
plot(xmult,ymult,1000,10,100,0.d0,1.d0, -3.5d0,0.05d0)
4038 IF (kpl.NE.12) go to 270
4041 aaso(i,j)=log10(aaso(i,j)+1.
d-18)
4042 avmult(i,j)=log10(avmult(i,j)+1.
d-18)
4043 ake(i,j)=log10(ake(i,j)+1.
d-18)
4045 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4046 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4047 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4048 CALL
plot(aaso,avmult,360,30,12,0.d0,0.1d0,-3.d0,0.05d0)
4049 WRITE(6,*)
' 1=*=p, 2=n, 3=pi+, 4=pi-, 5=K+, 6=K-,',
4050 &
' 8=ap, 9=chargd., 10=Z=all, 11=+=Lamb, 12=A=aLam,',
4051 &
' 13=O=sig-, 14=B=Sig+, 15=C=Sig0, 16=D=pi0'
4052 CALL
plot(aaso,ake,360,30,12,0.d0,0.1d0,-5.d0,0.05d0)
4061 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4066 parameter(nfimax=249)
4067 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
4068 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
4069 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
4076 parameter(idmax9=602)
4078 common/ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
4081 common/dpar/aname(210),am(210),ga(210),
tau(210),ich(210),ibar(210)
4083 common/dmetls/ cxs(149),cys(149),czs(149),els(149),
4084 *pls(149),is,its(149)
4085 common/ddre/
test(12)
4091 pls(i)=
sqrt(pxf(i)**2+pyf(i)**2+pzf(i)**2)
4092 IF(pls(i).NE.0.)cxs(i)=pxf(i)/pls(i)
4093 IF(pls(i).NE.0.)cys(i)=pyf(i)/pls(i)
4094 IF(pls(i).NE.0.)czs(i)=pzf(i)/pls(i)
4103 IF(istab.EQ.1) goto 30
4104 IF(istab.EQ.2) goto 50
4105 IF(istab.EQ.3) goto 40
4106 30
IF(its(ist).EQ.135.OR.its(ist).EQ.136) goto 60
4107 IF(its(ist).GE.1.AND.its(ist).LE.7) goto 60
4109 40
IF(its(ist).GE.1.AND.its(ist).LE.23) goto 60
4110 IF(its(ist).GE. 97.AND.its(ist).LE.103) goto 60
4112 IF(its(ist).EQ.109.OR.its(ist).EQ.115) goto 60
4113 IF(its(ist).GE.133.AND.its(ist).LE.136) goto 60
4115 50
IF(its(ist).GE. 1.AND.its(ist).LE. 30) goto 60
4116 IF(its(ist).GE. 97.AND.its(ist).LE.103) goto 60
4117 IF(its(ist).GE.115.AND.its(ist).LE.122) goto 60
4118 IF(its(ist).GE.131.AND.its(ist).LE.136) goto 60
4119 IF(its(ist).EQ.109) goto 60
4120 IF(its(ist).GE.137.AND.its(ist).LE.160) goto 60
4123 IF (ir.GT.nfimax)
THEN
4124 WRITE (6,1000)ir,nfimax
4125 1000
FORMAT(.GT.
' DECAY IRNFIMAX RETURN ',2i10)
4135 pxf(ir)=cxs(ist)*pls(ist)
4136 pyf(ir)=cys(ist)*pls(ist)
4137 pzf(ir)=czs(ist)*pls(ist)
4139 IF(ist.GE.1) goto 20
4143 bgam=pls(ist)/am(it)
4150 IF (vv.GT.wt(iik)) go to 90
4154 IF (it2-1.LT.0) go to 120
4157 IF(it3.EQ.0) go to 100
4158 CALL
dthrep(eco,ecm1,ecm2,ecm3,pcm1,pcm2,pcm3,cod1,cof1,sif1,
4159 *cod2,cof2,sif2,cod3,cof3,sif3,am(it1),am(it2),am(it3))
4161 100 CALL
dtwopd(eco,ecm1,ecm2,pcm1,pcm2,cod1,cof1,sif1,cod2,cof2,sif2,
4166 IF (it2-1.LT.0) go to 130
4172 CALL
dtrafo(gam,bgam,rx,ry,rz,cod1,cof1,sif1,pcm1,ecm1,
4173 *pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
4175 CALL
dtrafo(gam,bgam,rx,ry,rz,cod2,cof2,sif2,pcm2,ecm2,
4176 *pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
4177 IF (it3.LE.0) go to 130
4179 CALL
dtrafo(gam,bgam,rx,ry,rz,cod3,cof3,sif3,pcm3,ecm3,
4180 *pls(ist),cxs(ist),cys(ist),czs(ist),els(ist))
4184 IF(ir.GT.7998)
WRITE(isys,1010)
4185 1010
FORMAT(2
x,
' NUMBER OF STAB. FINAL PART. IS GREATER THAN 7998')
4195 SUBROUTINE shmak(ICASE,NN,NNA,NNB,NA,NB,UMO,BIMP)
4196 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4203 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4206 dimension fnua(namx),fnub(namx),fnut(namx)
4208 dimension xb(200),bimpp(200)
4240 go to(10,30,40),icase
4264 bimpp(ib)=bimpp(ib)+1
4266 IF((nn.EQ.1).AND.(nna.EQ.1).AND.(nnb.EQ.1))
THEN
4267 CALL
sihndi(umo,1,1,singdif,sigdih)
4269 anusd=anusd + singdif/sigabs
4272 IF (intt.GT.namx)intt=namx
4274 IF (nua.GT.namx) nua=namx
4276 IF (nub.GT.namx) nub=namx
4277 fnua(nua)=fnua(nua)+1.
4278 fnut(intt)=fnut(intt)+1.
4279 fnub(nub)=fnub(nub)+1.
4301 WRITE(6,*)
' shmak(3,NN,... ) NN= ',nn
4316 WRITE(6,
'(1H1,50(1H*))')
4317 WRITE(6,
'(/10X,A/)')
' OUTPUT FROM SHMAK all events before',
4318 *
' diffraction modification'
4319 WRITE(6,
'(50(1H*))')
4320 WRITE(6,
'(A,I10)')
' NUMBER OF TOTALLY SAMPLED GLAUBER EVENTS',nn
4321 WRITE(6, 1000) bnut,bnua,bnub
4322 WRITE(6,*)
' Fraction of diffractive evnts: ',anusd
4323 1000
FORMAT(
' AVERAGE NO OF COLLISIONS BNUT,BNUA,BNUB=',3f10.3)
4324 WRITE(6,
'(/A)')
' AVERAGE NUMBERS OF DIFFERENT COLLISION TYPES'
4325 WRITE(6,
'(4(5X,A,F8.2/))')
' VAL-VAL:',bnvv,
' SEA-VAL:',bnsv,
4326 +
' VAL-SEA:',bnvs,
' SEA-SEA:',bnss
4332 1010
FORMAT (
' FNUA')
4335 fnu
a(i)=log10(fnu
a(i)+1.
d-5)
4337 CALL
plot(ann,fnu
a,namx,1,namx,0.d0,dnna,0.d0,0.05d0)
4340 1020
FORMAT (
' FNUB')
4342 fnu b(i)=log10(fnu b(i)+1.
d-5)
4344 CALL
plot(ann,fnu b,namx,1,namx,0.d0,dnnb,0.d0,0.05d0)
4347 1030
FORMAT (
' FNUT')
4349 fnu
t(i)=log10(fnu
t(i)+1.
e-5)
4351 CALL
plot(ann,fnu
t,namx,1,namx,0.d0,dnnt,0.d0,0.05d0)
4352 1040
FORMAT (10f12.2)
4361 SUBROUTINE shmak1(ICASE,NN,NNA,NNB,NA,NB,UMO,BIMP)
4362 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4369 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4372 dimension fnua(namx),fnub(namx),fnut(namx)
4374 dimension xb(200),bimpp(200)
4376 go to(10,30,40),icase
4396 bimpp(ib)=bimpp(ib)+1
4398 IF((nn.EQ.1).AND.(nna.EQ.1).AND.(nnb.EQ.1))
THEN
4399 CALL
sihndi(umo,1,1,singdif,sigdih)
4401 anusd=anusd + singdif/sigabs
4404 IF (intt.GT.namx)intt=namx
4406 IF (nua.GT.namx) nua=namx
4408 IF (nub.GT.namx) nub=namx
4409 fnua(nua)=fnua(nua)+1.
4410 fnut(intt)=fnut(intt)+1.
4411 fnub(nub)=fnub(nub)+1.
4418 WRITE(6,*)
' shmak1(3,NN,... ) NN= ',nn
4429 WRITE(6,
'(1H1,50(1H*))')
4430 WRITE(6,
'(/10X,A/)')
' OUTPUT FROM SHMAK1 after modification',
4431 *
' of Glauber events for diffractive cross section'
4432 WRITE(6,
'(50(1H*))')
4433 WRITE(6,
'(A,I10)')
' NUMBER OF TOTALLY SAMPLED GLAUBER EVENTS',nn
4434 WRITE(6, 1000) bnut,bnua,bnub
4435 1000
FORMAT(
' AVERAGE NO OF COLLISIONS BNUT,BNUA,BNUB=',3f10.3)
4436 WRITE(6,*)
' Fraction of diffractive evnts: ',anusd
4445 1010
FORMAT (
' FNUA')
4448 fnu
a(i)=log10(fnu
a(i)+1.
d-5)
4450 CALL
plot(ann,fnu
a,namx,1,namx,0.d0,dnna,0.d0,0.05d0)
4453 1020
FORMAT (
' FNUB')
4455 fnu b(i)=log10(fnu b(i)+1.
d-5)
4457 CALL
plot(ann,fnu b,namx,1,namx,0.d0,dnnb,0.d0,0.05d0)
4460 1030
FORMAT (
' FNUT')
4462 fnu
t(i)=log10(fnu
t(i)+1.
e-5)
4464 CALL
plot(ann,fnu
t,namx,1,namx,0.d0,dnnt,0.d0,0.05d0)
4465 1040
FORMAT (10f12.2)
4473 SUBROUTINE previo(RA,RB,NSTB,BMAX,BSTEP,SIG,RO,G)
4474 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4477 common/damp/ca,ci,ga
4502 rca=ga*sig/6.2831854
4507 fca=-ga*sig*ro/6.2831854
4518 WRITE(6,*)
' PREVIO(RA,RB,NSTB,BMAX,BSTEP,SIG,RO,G)',
4519 &ra,rb,nstb,bmax,bstep,sig,ro,g
4520 WRITE(6,*)
' /CA,CI,GA/ ',ca,ci,ga
4521 WRITE(6,*)
' PREVIO: RA, RB ,CI= ',ra,rb,ci
4528 SUBROUTINE profb(BSTEP,NSTAT,NA,RA,NB,RB,BSITE,NSITEB)
4529 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4537 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
4539 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
4540 +ipadis,ishmal,lpauli
4546 DOUBLE COMPLEX ca,ci
4547 COMMON /damp/ ca,ci,ga
4549 dimension bsite(0:1,nsiteb)
4551 dimension helpp(200)
4554 COMMON /sigla/siglau
4559 WRITE(6,*)
' PROFB: RA, RB = ',ra,rb
4560 WRITE(6, 1000)bstep,nstat,na,ra,nb,rb,irw,nsiteb
4561 1000
FORMAT (
' PROFB',e15.5,2i10,f15.5,i10,e15.5,2i10)
4578 IF(pi.LT.1.
d-100)go to 31
4585 IF(
xy.GT.15.) go to 20
4598 bs(i3+1)=bs(i3+1)+1.-pi
4604 bs(i)=bs(i)*(i-1)*bst/
ns
4609 bsite(1,i)=bs(i)/sumb+bsite(1,i-1)
4614 sumb=sumb*bst*6.2831854
4617 1020
FORMAT(/5
x,7hsigma =,f7.3)
4621 WRITE (6,1030) help(i),helpp(i),bs(i),bsite(1,i)
4622 1030
FORMAT (f10.4,3e15.5)
4624 CALL
plot(help,bsite,50,1,50,0.d0,0.5d0,0.d0,0.01d0)
4625 CALL
plot(help,bs ,50,1,50,0.d0,0.5d0,0.d0,0.07d0)
4634 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4638 parameter(nfimax=249)
4639 COMMON /dfinpa/ anf(nfimax),pxf(nfimax),pyf(nfimax),pzf(nfimax),
4640 +hef(nfimax),amf(nfimax), ichf(nfimax),ibarf(nfimax),nref(nfimax)
4641 COMMON /dfinpz/iormo(nfimax),idaug1(nfimax),idaug2(nfimax),
4644 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
4645 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
4647 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
4650 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210), ich(210),ibar
4651 +(210),k1(210),k2(210)
4663 WRITE(6,1000)ihad,i,pxf(1),pyf(1),pzf(1),hef(1),amf(1)
4664 1000
FORMAT(
' PARJET: IHAD,I,PXF(1),PYF(1),PZF(1),HEP(1),AMF(1)'/ 2i5,5
4674 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4684 IF (
a(3,i).LE.
a(3,j)) go to 20
4704 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4714 IF (
a(3,i).GE.
a(3,j)) go to 20
4733 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
4780 parameter( kalgnm = 2 )
4781 parameter( anglgb = 5.0
d-16 )
4782 parameter( anglsq = 2.5
d-31 )
4783 parameter( axcssv = 0.2
d+16 )
4784 parameter( andrfl = 1.0
d-38 )
4785 parameter( avrflw = 1.0
d+38 )
4786 parameter( ainfnt = 1.0
d+30 )
4787 parameter( azrzrz = 1.0
d-30 )
4788 parameter( einfnt = +69.07755278982137
d+00 )
4789 parameter( ezrzrz = -69.07755278982137
d+00 )
4790 parameter( onemns = 0.999999999999999
d+00 )
4791 parameter( onepls = 1.000000000000001
d+00 )
4792 parameter( csnnrm = 2.0
d-15 )
4793 parameter( dmxtrn = 1.0
d+08 )
4825 parameter( zerzer = 0.
d+00 )
4826 parameter( oneone = 1.
d+00 )
4827 parameter( twotwo = 2.
d+00 )
4828 parameter( thrthr = 3.
d+00 )
4829 parameter( foufou = 4.
d+00 )
4830 parameter( fivfiv = 5.
d+00 )
4831 parameter( sixsix = 6.
d+00 )
4832 parameter( sevsev = 7.
d+00 )
4833 parameter( eigeig = 8.
d+00 )
4834 parameter( aninen = 9.
d+00 )
4835 parameter( tenten = 10.
d+00 )
4836 parameter( hlfhlf = 0.5
d+00 )
4837 parameter( onethi = oneone / thrthr )
4838 parameter( twothi = twotwo / thrthr )
4839 parameter( pipipi = 3.1415926535897932270
d+00 )
4840 parameter( eneper = 2.7182818284590452354
d+00 )
4841 parameter( sqrent = 1.6487212707001281468
d+00 )
4883 parameter( clight = 2.99792458
d+10 )
4884 parameter( avogad = 6.0221367
d+23 )
4885 parameter( amelgr = 9.1093897
d-28 )
4886 parameter( plckbr = 1.05457266
d-27 )
4887 parameter( elccgs = 4.8032068
d-10 )
4888 parameter( elcmks = 1.60217733
d-19 )
4889 parameter( amugrm = 1.6605402
d-24 )
4890 parameter( ammumu = 0.113428913
d+00 )
4905 parameter( alpfsc = 7.2973530791728595
d-03 )
4906 parameter( fscto2 = 5.3251361962113614
d-05 )
4907 parameter( fscto3 = 3.8859399018437826
d-07 )
4908 parameter( fscto4 = 2.8357075508200407
d-09 )
4909 parameter( plabrc = 0.197327053
d+00 )
4910 parameter( amelct = 0.51099906
d-03 )
4911 parameter( amugev = 0.93149432
d+00 )
4912 parameter( ammuon = 0.105658389
d+00 )
4913 parameter( rclsel = 2.8179409183694872
d-13 )
4914 parameter( gevmev = 1.0
d+03 )
4915 parameter( emvgev = 1.0
d-03 )
4916 parameter( raddeg = 180.
d+00 / pipipi )
4917 parameter( degrad = pipipi / 180.
d+00 )
4944 parameter( lunin = 5 )
4945 parameter( lunout = 6 )
4946 parameter( lunerr = 66 )
4947 parameter( lunber = 14 )
4948 parameter( lunech = 8 )
4949 parameter( lunflu = 86 )
4950 parameter( lungeo = 16 )
4951 parameter( lunpgs = 12 )
4952 parameter( lunran = 2 )
4953 parameter( lunxsc = 81 )
4954 parameter( lunrdb = 1 )
4983 parameter( mxxrgn = 500 )
4984 parameter( mxxmdf = 56 )
4985 parameter( mxxmde = 50 )
4986 parameter( mfstck = 1000 )
4987 parameter( mestck = 100 )
4988 parameter( nallwp = 39 )
4989 parameter( mpdpdx = 8 )
4990 parameter( icomax = 180 )
4991 parameter( nstbis = 304 )
4992 parameter( idmaxp = 210 )
4996 COMMON /dpar/ aname(210),am(210),ga(210),
tau(210),
4997 + ich(210),ibar(210),k1(210),k2(210)
5005 DATA (am(k),k=1,85) /
5006 & .9383
d+00, .9383
d+00, amelct , amelct , .0000
d+00,
5007 & .0000
d+00, .0000
d+00, .9396
d+00, .9396
d+00, ammuon ,
5008 & ammuon , .4977
d+00, .1396
d+00, .1396
d+00, .4936
d+00,
5009 & .4936
d+00, .1116
d+01, .1116
d+01, .4977
d+00, .1197
d+01,
5010 & .1189
d+01, .1193
d+01, .1350
d+00, .4977
d+00, .4977
d+00,
5011 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5012 & .5488
d+00, .7669
d+00, .7700
d+00, .7669
d+00, .7820
d+00,
5013 & .8921
d+00, .8962
d+00, .8921
d+00, .8962
d+00, .1300
d+01,
5014 & .1300
d+01, .1300
d+01, .1300
d+01, .1421
d+01, .1421
d+01,
5015 & .1421
d+01, .1421
d+01, .1383
d+01, .1384
d+01, .1387
d+01,
5016 & .1820
d+01, .2030
d+01, .1231
d+01, .1232
d+01, .1233
d+01,
5017 & .1234
d+01, .1675
d+01, .1675
d+01, .1675
d+01, .1675
d+01,
5018 & .1500
d+01, .1500
d+01, .1515
d+01, .1515
d+01, .1775
d+01,
5019 & .1775
d+01, .1231
d+01, .1232
d+01, .1233
d+01, .1234
d+01,
5020 & .1675
d+01, .1675
d+01, .1675
d+01, .1675
d+01, .1515
d+01,
5021 & .1515
d+01, .2500
d+01, .4890
d+00, .4890
d+00, .4890
d+00,
5022 & .1300
d+01, .1300
d+01, .1300
d+01, .1300
d+01, .2200
d+01 /
5023 DATA (am(k),k=86,183) /
5024 & .2200
d+01, .2200
d+01, .2200
d+01, .1700
d+01, .1700
d+01,
5025 & .1700
d+01, .1700
d+01, .1820
d+01, .2030
d+01, .9575
d+00,
5026 & .1019
d+01, .1315
d+01, .1321
d+01, .1189
d+01, .1193
d+01,
5027 & .1197
d+01, .1315
d+01, .1321
d+01, .1383
d+01, .1384
d+01,
5028 & .1387
d+01, .1532
d+01, .1535
d+01, .1672
d+01, .1383
d+01,
5029 & .1384
d+01, .1387
d+01, .1532
d+01, .1535
d+01, .1672
d+01,
5030 & .1865
d+01, .1869
d+01, .1869
d+01, .1865
d+01, .1969
d+01,
5031 & .1969
d+01, .2980
d+01, .2007
d+01, .2010
d+01, .2010
d+01,
5032 & .2007
d+01, .2113
d+01, .2113
d+01, .3686
d+01, .3097
d+01,
5033 & .1777
d+01, .1777
d+01, .0000
d+00, .0000
d+00, .0000
d+00,
5034 & .0000
d+00, .2285
d+01, .2460
d+01, .2460
d+01, .2452
d+01,
5035 & .2453
d+01, .2454
d+01, .2560
d+01, .2560
d+01, .2730
d+01,
5036 & .3610
d+01, .3610
d+01, .3790
d+01, .2285
d+01, .2460
d+01,
5037 & .2460
d+01, .2452
d+01, .2453
d+01, .2454
d+01, .2560
d+01,
5038 & .2560
d+01, .2730
d+01, .3610
d+01, .3610
d+01, .3790
d+01,
5039 & .2490
d+01, .2490
d+01, .2490
d+01, .2610
d+01, .2610
d+01,
5040 & .2770
d+01, .3670
d+01, .3670
d+01, .3850
d+01, .4890
d+01,
5041 & .2490
d+01, .2490
d+01, .2490
d+01, .2610
d+01, .2610
d+01,
5042 & .2770
d+01, .3670
d+01, .3670
d+01, .3850
d+01, .4890
d+01,
5043 & .1250
d+01, .1250
d+01, .1250
d+01 /
5044 DATA ( am( i ), i = 184,210 ) /
5045 & 1.44000000000000
d+00, 1.44000000000000
d+00, 1.30000000000000
d+00,
5046 & 1.30000000000000
d+00, 1.30000000000000
d+00, 1.40000000000000
d+00,
5047 & 1.46000000000000
d+00, 1.46000000000000
d+00, 1.46000000000000
d+00,
5048 & 1.46000000000000
d+00, 1.60000000000000
d+00, 1.60000000000000
d+00,
5049 & 1.66000000000000
d+00, 1.66000000000000
d+00, 1.66000000000000
d+00,
5050 & 1.66000000000000
d+00, 1.66000000000000
d+00, 1.66000000000000
d+00,
5051 & 1.95000000000000
d+00, 1.95000000000000
d+00, 1.95000000000000
d+00,
5052 & 1.95000000000000
d+00, 2.25000000000000
d+00, 2.25000000000000
d+00,
5053 & 1.44000000000000
d+00, 1.44000000000000
d+00, 0.00000000000000
d+00/
5057 DATA (
tau(k),k=1,183) /
5058 & .1000
d+19, .1000
d+19, .1000
d+19, .1000
d+19, .1000
d+19,
5059 & .1000
d+19, .1000
d+19, .9180
d+03, .9180
d+03, .2200
d-05,
5060 & .2200
d-05, .5200
d-07, .2600
d-07, .2600
d-07, .1200
d-07,
5061 & .1200
d-07, .2600
d-09, .2600
d-09, .9000
d-10, .1500
d-09,
5062 & .8000
d-10, .5000
d-14, .8000
d-16, .0000
d+00, .0000
d+00,
5064 & .0000
d+00, .3000
d-09, .1700
d-09, .8000
d-10, .1000
d-13,
5065 & .1500
d-09, .3000
d-09, .1700
d-09, .0000
d+00, .0000
d+00,
5066 & .0000
d+00, .0000
d+00, .0000
d+00, .1000
d-09, .0000
d+00,
5067 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .1000
d-09,
5068 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5069 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5070 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5071 & .9000
d-11, .9000
d-11, .9000
d-11, .9000
d-11, .1000
d+19,
5072 & .1000
d+19, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5074 & .0000
d+00, .0000
d+00, .0000
d+00 /
5075 DATA (
tau( i ), i = 184,210 ) /
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,
5082 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5083 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00,
5084 & 0.00000000000000
d+00, 0.00000000000000
d+00, 0.00000000000000
d+00/
5088 DATA (ga(k),k= 1,85) /
5090 & .8500
d-06, .1520
d+00, .1520
d+00, .1520
d+00, .1000
d-01,
5091 & .7900
d-01, .7900
d-01, .7900
d-01, .7900
d-01, .4500
d+00,
5092 & .4500
d+00, .4500
d+00, .4500
d+00, .1080
d+00, .1080
d+00,
5093 & .1080
d+00, .1080
d+00, .5000
d-01, .5000
d-01, .5000
d-01,
5094 & .8500
d-01, .1800
d+00, .1150
d+00, .1150
d+00, .1150
d+00,
5095 & .1150
d+00, .2000
d+00, .2000
d+00, .2000
d+00, .2000
d+00,
5096 & .2000
d+00, .2000
d+00, .1000
d+00, .1000
d+00, .2000
d+00,
5097 & .2000
d+00, .1150
d+00, .1150
d+00, .1150
d+00, .1150
d+00,
5098 & .2000
d+00, .2000
d+00, .2000
d+00, .2000
d+00, .1000
d+00,
5099 & .1000
d+00, .2000
d+00, .1000
d+00, .1000
d+00, .1000
d+00,
5100 & .1000
d+00, .1000
d+00, .1000
d+00, .1000
d+00, .2000
d+00 /
5101 DATA (ga(k),k= 86,183) /
5102 & .2000
d+00, .2000
d+00, .2000
d+00, .1500
d+00, .1500
d+00,
5103 & .1500
d+00, .1500
d+00, .8500
d-01, .1800
d+00, .2000
d-02,
5104 & .4000
d-02, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5105 & .0000
d+00, .0000
d+00, .0000
d+00, .3400
d-01, .3400
d-01,
5106 & .3600
d-01, .9000
d-02, .9000
d-02, .0000
d+00, .3400
d-01,
5107 & .3400
d-01, .3600
d-01, .9000
d-02, .9000
d-02, .0000
d+00,
5108 & .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00, .0000
d+00,
5109 & .0000
d+00, .0000
d+00, .5000
d-02, .2000
d-02, .2000
d-02,
5110 & .5000
d-02, .2000
d-02, .2000
d-02, .2000
d-03, .7000
d-03,
5112 & .3000
d+00, .3000
d+00, .3000
d+00 /
5113 DATA ( ga( i ), i = 184,210 ) /
5114 & 2.00000000000000
d-01, 2.00000000000000
d-01, 3.00000000000000
d-01,
5115 & 3.00000000000000
d-01, 3.00000000000000
d-01, 2.70000000000000
d-01,
5116 & 2.50000000000000
d-01, 2.50000000000000
d-01, 2.50000000000000
d-01,
5117 & 2.50000000000000
d-01, 1.50000000000000
d-01, 1.50000000000000
d-01,
5118 & 1.00000000000000
d-01, 1.00000000000000
d-01, 1.00000000000000
d-01,
5119 & 1.00000000000000
d-01, 1.00000000000000
d-01, 1.00000000000000
d-01,
5120 & 6.00000000000000
d-02, 6.00000000000000
d-02, 6.00000000000000
d-02,
5121 & 6.00000000000000
d-02, 5.50000000000000
d-02, 5.50000000000000
d-02,
5122 & 2.00000000000000
d-01, 2.00000000000000
d-01, 0.00000000000000
d+00/
5131 DATA (aname(k),k=1,85) /
5132 &
'P ',
'AP ',
'E- ',
'E+ ',
'NUE ',
5133 &
'ANUE ',
'GAM ',
'NEU ',
'ANEU ',
'MUE+ ',
5134 &
'MUE- ',
'K0L ',
'PI+ ',
'PI- ',
'K+ ',
5135 &
'K- ',
'LAM ',
'ALAM ',
'K0S ',
'SIGM- ',
5136 &
'SIGM+ ',
'SIGM0 ',
'PI0 ',
'K0 ',
'AK0 ',
5137 &
'BLANK ',
'BLANK ',
'BLANK ',
'BLANK ',
'BLANK ',
5138 &
'ETA550 ',
'RHO+77 ',
'RHO077 ',
'RHO-77 ',
'OM0783 ',
5139 &
'K*+892 ',
'K*0892 ',
'K*-892 ',
'AK*089 ',
'KA+125 ',
5140 &
'KA0125 ',
'KA-125 ',
'AKA012 ',
'K*+142 ',
'K*0142 ',
5141 &
'K*-142 ',
'AK*014 ',
'S+1385 ',
'S01385 ',
'S-1385 ',
5142 &
'L01820 ',
'L02030 ',
'N*++12 ',
'N*+ 12 ',
'N*012 ',
5143 &
'N*-12 ',
'N*++16 ',
'N*+16 ',
'N*016 ',
'N*-16 ',
5144 &
'N*+14 ',
'N*014 ',
'N*+15 ',
'N*015 ',
'N*+18 ',
5145 &
'N*018 ',
'AN--12 ',
'AN*-12 ',
'AN*012 ',
'AN*+12 ',
5146 &
'AN--16 ',
'AN*-16 ',
'AN*016 ',
'AN*+16 ',
'AN*-15 ',
5147 &
'AN*015 ',
'DE*=24 ',
'RPI+49 ',
'RPI049 ',
'RPI-49 ',
5148 &
'PIN++ ',
'PIN+0 ',
'PIN+- ',
'PIN-0 ',
'PPPI ' /
5149 DATA (aname(k),k=86,183) /
5150 &
'PNPI ',
'APPPI ',
'APNPI ',
'K+PPI ',
'K-PPI ',
5151 &
'K+NPI ',
'K-NPI ',
'S+1820 ',
'S-2030 ',
'ETA* ',
5152 &
'PHI ',
'TETA0 ',
'TETA- ',
'ASIG- ',
'ASIG0 ',
5153 &
'ASIG+ ',
'ATETA0 ',
'ATETA+ ',
'SIG*+ ',
'SIG*0 ',
5154 &
'SIG*- ',
'TETA*0 ',
'TETA* ',
'OMEGA- ',
'ASIG*- ',
5155 &
'ASIG*0 ',
'ASIG*+ ',
'ATET*0 ',
'ATET*+ ',
'OMEGA+ ',
5156 &
'D0 ',
'D+ ',
'D- ',
'AD0 ',
'DS+ ',
5157 &
'DS- ',
'ETAC ',
'D*0 ',
'D*+ ',
'D*- ',
5158 &
'AD*0 ',
'DS*+ ',
'DS*- ',
'CHI1C ',
'JPSI ',
5159 &
'TAU+ ',
'TAU- ',
'NUET ',
'ANUET ',
'NUEM ',
5160 &
'ANUEM ',
'LAMC+ ',
'XIC+ ',
'XIC0 ',
'SIGC++ ',
5161 &
'SIGC+ ',
'SIGC0 ',
'S+ ',
'S0 ',
'T0 ',
5162 &
'XU++ ',
'XD+ ',
'XS+ ',
'ALAMC- ',
'AXIC- ',
5163 &
'AXIC0 ',
'ASIGC-- ',
'ASIGC- ',
'ASIGC0 ',
'AS- ',
5164 &
'AS0 ',
'AT0 ',
'AXU-- ',
'AXD- ',
'AXS ',
5165 &
'C1*++ ',
'C1*+ ',
'C1*0 ',
'S*+ ',
'S*0 ',
5166 &
'T*0 ',
'XU*++ ',
'XD*+ ',
'XS*+ ',
'TETA++ ',
5167 &
'AC1*-- ',
'AC1*- ',
'AC1*0 ',
'AS*- ',
'AS*0 ',
5168 &
'AT*0 ',
'AXU*-- ',
'AXD*- ',
'AXS*- ',
'ATET-- ',
5169 &
'RO ',
'R+ ',
'R- ' /
5170 DATA ( aname( i ), i = 184,210 ) /
5171 &
'AN*-14 ',
'AN*014 ',
'PI+130 ',
'PI0130 ',
'PI-130 ',
'F01400 ',
5172 &
'K*+146 ',
'K*-146 ',
'K*0146 ',
'AK0146 ',
'L01600 ',
'AL0160 ',
5173 &
'S+1660 ',
'S01660 ',
'S-1660 ',
'AS-166 ',
'AS0166 ',
'AS+166 ',
5174 &
'X01950 ',
'X-1950 ',
'AX0195 ',
'AX+195 ',
'OM-225 ',
'AOM+22 ',
5175 &
'N*+14 ',
'N*014 ',
'BLANK '/
5179 DATA ( ich( i ), i = 1,210 ) /
5180 & 1, -1, -1, 1, 0, 0, 0, 0, 0, 1, -1, 0, 1, -1, 1,
5181 & -1, 0, 0, 0, -1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5182 & 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0, -1, 0, 1, 0,
5183 & -1, 0, 1, 0, -1, 0, 0, 2, 1, 0, -1, 2, 1, 0, -1,
5184 & 1, 0, 1, 0, 1, 0, -2, -1, 0, 1, -2, -1, 0, 1, -1,
5185 & 0, 1, 1, 0, -1, 2, 1, 0, -1, 2, 1, 0, -1, 2, 0,
5186 & 1, -1, 1, -1, 0, 0, 0, -1, -1, 0, 1, 0, 1, 1, 0,
5187 & -1, 0, -1, -1, -1, 0, 1, 0, 1, 1, 0, 1, -1, 0, 1,
5188 & -1, 0, 0, 1, -1, 0, 1, -1, 0, 0, 1, -1, 0, 0, 0,
5189 & 0, 1, 1, 0, 2, 1, 0, 1, 0, 0, 2, 1, 1, -1, -1,
5190 & 0, -2, -1, 0, -1, 0, 0, -2, -1, -1, 2, 1, 0, 1, 0,
5191 & 0, 2, 1, 1, 2, -2, -1, 0, -1, 0, 0, -2, -1, -1, -2,
5192 & 0, 1, -1, -1, 0, 1, 0, -1, 0, 1, -1, 0, 0, 0, 0,
5193 & 1, 0, -1, -1, 0, 1, 0, -1, 0, 1, -1, 1, 1, 0, 0/
5197 DATA ( ibar( i ), i = 1,210 ) /
5198 & 1, -1, 0, 0, 0, 0, 0, 1, -1, 0, 0, 0, 0, 0, 0,
5199 & 0, 1, -1, 0, 1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0,
5200 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5201 & 0, 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
5202 & 1, 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
5203 & -1, 2, 0, 0, 0, 1, 1, 1, 1, 2, 2, 0, 0, 1, 1,
5204 & 1, 1, 1, 1, 0, 0, 1, 1, -1, -1, -1, -1, -1, 1, 1,
5205 & 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, 0, 0, 0, 0, 0,
5206 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5207 & 0, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1, -1,
5208 & -1, -1, -1, -1, -1, -1, -1, -1, -1, -1, 1, 1, 1, 1, 1,
5209 & 1, 1, 1, 1, 1, -1, -1, -1, -1, -1, -1, -1, -1, -1, -1,
5210 & 0, 0, 0, -1, -1, 0, 0, 0, 0, 0, 0, 0, 0, 1, -1,
5211 & 1, 1, 1, -1, -1, -1, 1, 1, -1, -1, 1, -1, 1, 1, 0/
5216 DATA k1/ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 16, 17,
5217 & 18, 24, 30, 34, 38, 40, 41, 43, 44, 136, 138, 330, 327, 328,
5218 & 2*330, 46, 51, 52, 54, 55, 58,
5219 & 60, 62, 64, 66, 68, 70, 72, 74, 82, 90, 98, 106, 109, 112, 114,
5220 & 123, 140, 141, 143, 145, 146, 150, 157, 164, 168, 174, 180, 187,
5221 & 194, 202, 210, 211, 213, 215, 216, 220, 227, 234, 238, 245, 252,
5222 & 254, 255, 256, 257, 259, 262, 265, 267, 269, 272, 276, 279, 282,
5223 & 286, 290, 293, 299, 331, 335, 339, 340, 341, 343, 344, 345, 346,
5224 & 347, 350, 353, 356, 358, 360, 363, 366, 369, 372, 374, 376, 379,
5225 & 383, 385, 387, 391, 394, 397, 400, 402, 405, 408, 410, 412, 414,
5226 & 417, 420, 425, 430, 431, 432, 433, 434, 448, 452, 457, 458, 459,
5227 & 460, 461, 462, 466, 468, 470, 472, 486, 490, 495, 496, 497, 498,
5228 & 499, 500, 504, 506, 508, 510, 511, 512, 513, 514, 515, 516, 517,
5229 & 518, 519, 522, 523, 524, 525, 526, 527, 528, 529, 530, 531, 534,
5230 & 537, 539, 541, 547, 553, 558, 563, 568, 572, 573, 574, 575, 576,
5231 & 577, 578, 579, 580, 581, 582, 583, 584, 585, 586, 587, 588, 589,
5237 DATA k2/ 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 15, 16, 17,
5238 & 23, 29, 31, 35, 39, 40, 42, 43, 45, 137, 139, 330, 327, 328,
5239 & 2* 330, 50, 51, 53, 54, 57,
5240 & 59, 61, 63, 65, 67, 69, 71, 73, 81, 89, 97, 105, 108, 111, 113,
5241 & 122, 135, 140, 142, 144, 145, 149, 156, 163, 167, 173, 179, 186,
5242 & 193, 201, 209, 210, 212, 214, 215, 219, 226, 233, 237, 244, 251,
5243 & 253, 254, 255, 256, 258, 261, 264, 266, 268, 271, 275, 278, 281,
5244 & 285, 289, 292, 298, 307, 334, 338, 339, 340, 342, 343, 344, 345,
5245 & 346, 349, 352, 355, 357, 359, 362, 365, 368, 371, 373, 375, 378,
5246 & 382, 384, 386, 390, 393, 396, 399, 401, 404, 407, 409, 411, 413,
5247 & 416, 419, 424, 429, 430, 431, 432, 433, 447, 451, 456, 457, 458,
5248 & 459, 460, 461, 465, 467, 469, 471, 485, 489, 494, 495, 496, 497,
5249 & 498, 499, 503, 505, 507, 509, 510, 511, 512, 513, 514, 515, 516,
5250 & 517, 518, 521, 522, 523, 524, 525, 526, 527, 528, 529, 530, 533,
5251 & 536, 538, 540, 546, 552, 557, 562, 567, 571, 572, 573, 574, 575,
5252 & 576, 577, 578, 579, 580, 581, 582, 583, 584, 585, 586, 587, 588,
5253 & 589, 595, 601, 602 /
5262 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5267 parameter(idmax9=602)
5269 common/ddecac/ zkname(idmax9),wt(idmax9),nzk(idmax9,3)
5280 DATA (zkname(k),k= 1, 85) /
5281 &
'P ',
'AP ',
'E- ',
'E+ ',
'NUE ',
5282 &
'ANUE ',
'GAM ',
'PE-NUE ',
'APEANU ',
'EANUNU ',
5283 &
'E-NUAN ',
'3PI0 ',
'PI+-0 ',
'PIMUNU ',
'PIE-NU ',
5284 &
'MU+NUE ',
'MU-NUE ',
'MU+NUE ',
'PI+PI0 ',
'PI++- ',
5285 &
'PI+00 ',
'M+P0NU ',
'E+P0NU ',
'MU-NU ',
'PI-0 ',
5286 &
'PI+-- ',
'PI-00 ',
'M-P0NU ',
'E-P0NU ',
'PPI- ',
5287 &
'NPI0 ',
'PD-NUE ',
'PM-NUE ',
'APPI+ ',
'ANPI0 ',
5288 &
'APE+NU ',
'APM+NU ',
'PI+PI- ',
'PI0PI0 ',
'NPI- ',
5289 &
'PPI0 ',
'NPI+ ',
'LAGA ',
'GAGA ',
'GAE+E- ',
5290 &
'GAGA ',
'GAGAP0 ',
'PI000 ',
'PI+-0 ',
'PI+-GA ',
5291 &
'PI+0 ',
'PI+- ',
'PI00 ',
'PI-0 ',
'PI+-0 ',
5292 &
'PI+- ',
'PI0GA ',
'K+PI0 ',
'K0PI+ ',
'KOPI0 ',
5293 &
'K+PI- ',
'K-PI0 ',
'AK0PI- ',
'AK0PI0 ',
'K-PI+ ',
5294 &
'K+PI0 ',
'K0PI+ ',
'K0PI0 ',
'K+PI- ',
'K-PI0 ',
5295 &
'K0PI- ',
'AK0PI0 ',
'K-PI+ ',
'K+PI0 ',
'K0PI+ ',
5296 &
'K+89P0 ',
'K08PI+ ',
'K+RO77 ',
'K0RO+7 ',
'K+OM07 ',
5297 &
'K+E055 ',
'K0PI0 ',
'K+PI+ ',
'K089P0 ',
'K+8PI- ' /
5298 DATA (zkname(k),k= 86,170) /
5299 &
'K0R077 ',
'K+R-77 ',
'K+R-77 ',
'K0OM07 ',
'K0E055 ',
5300 &
'K-PI0 ',
'K0PI- ',
'K-89P0 ',
'AK08P- ',
'K-R077 ',
5301 &
'AK0R-7 ',
'K-OM07 ',
'K-E055 ',
'AK0PI0 ',
'K-PI+ ',
5302 &
'AK08P0 ',
'K-8PI+ ',
'AK0R07 ',
'AK0OM7 ',
'AK0E05 ',
5303 &
'LA0PI+ ',
'SI0PI+ ',
'SI+PI0 ',
'LA0PI0 ',
'SI+PI- ',
5304 &
'SI-PI+ ',
'LA0PI- ',
'SI0PI- ',
'NEUAK0 ',
'PK- ',
5305 &
'SI+PI- ',
'SI0PI0 ',
'SI-PI+ ',
'LA0ET0 ',
'S+1PI- ',
5306 &
'S-1PI+ ',
'SO1PI0 ',
'NEUAK0 ',
'PK- ',
'LA0PI0 ',
5307 &
'LA0OM0 ',
'LA0RO0 ',
'SI+RO- ',
'SI-RO+ ',
'SI0RO0 ',
5308 &
'LA0ET0 ',
'SI0ET0 ',
'SI+PI- ',
'SI-PI+ ',
'SI0PI0 ',
5309 &
'K0S ',
'K0L ',
'K0S ',
'K0L ',
'P PI+ ',
5310 &
'P PI0 ',
'N PI+ ',
'P PI- ',
'N PI0 ',
'N PI- ',
5311 &
'P PI+ ',
'N*#PI0 ',
'N*+PI+ ',
'PRHO+ ',
'P PI0 ',
5312 &
'N PI+ ',
'N*#PI- ',
'N*+PI0 ',
'N*0PI+ ',
'PRHO0 ',
5313 &
'NRHO+ ',
'P PI- ',
'N PI0 ',
'N*+PI- ',
'N*0PI0 ',
5314 &
'N*-PI+ ',
'PRHO- ',
'NRHO0 ',
'N PI- ',
'N*0PI- ',
5315 &
'N*-PI0 ',
'NRHO- ',
'PETA0 ',
'N*#PI- ',
'N*+PI0 ' /
5316 DATA (zkname(k),k=171,255) /
5317 &
'N*0PI+ ',
'PRHO0 ',
'NRHO+ ',
'NETA0 ',
'N*+PI- ',
5318 &
'N*0PI0 ',
'N*-PI+ ',
'PRHO- ',
'NRHO0 ',
'P PI0 ',
5319 &
'N PI+ ',
'N*#PI- ',
'N*+PI0 ',
'N*0PI+ ',
'PRHO0 ',
5320 &
'NRHO+ ',
'P PI- ',
'N PI0 ',
'N*+PI- ',
'N*0PI0 ',
5321 &
'N*-PI+ ',
'PRHO- ',
'NRHO0 ',
'P PI0 ',
'N PI+ ',
5322 &
'PRHO0 ',
'NRHO+ ',
'LAMK+ ',
'S+ K0 ',
'S0 K+ ',
5323 &
'PETA0 ',
'P PI- ',
'N PI0 ',
'PRHO- ',
'NRHO0 ',
5324 &
'LAMK0 ',
'S0 K0 ',
'S- K+ ',
'NETA/ ',
'APPI- ',
5325 &
'APPI0 ',
'ANPI- ',
'APPI+ ',
'ANPI0 ',
'ANPI+ ',
5326 &
'APPI- ',
'AN*=P0 ',
'AN*-P- ',
'APRHO- ',
'APPI0 ',
5327 &
'ANPI- ',
'AN*=P+ ',
'AN*-P0 ',
'AN*0P- ',
'APRHO0 ',
5328 &
'ANRHO- ',
'APPI+ ',
'ANPI0 ',
'AN*-P+ ',
'AN*0P0 ',
5329 &
'AN*+P- ',
'APRHO+ ',
'ANRHO0 ',
'ANPI+ ',
'AN*0P+ ',
5330 &
'AN*+P0 ',
'ANRHO+ ',
'APPI0 ',
'ANPI- ',
'AN*=P+ ',
5331 &
'AN*-P0 ',
'AN*0P- ',
'APRHO0 ',
'ANRHO- ',
'APPI+, ',
5332 &
'ANPI0 ',
'AN*-P+ ',
'AN*0P0 ',
'AN*+P- ',
'APRHO+ ',
5333 &
'ANRHO0 ',
'PN*014 ',
'NN*=14 ',
'PI+0 ',
'PI+- ' /
5334 DATA (zkname(k),k=256,340) /
5335 &
'PI-0 ',
'P+0 ',
'N++ ',
'P+- ',
'P00 ',
5336 &
'N+0 ',
'N+- ',
'N00 ',
'P-0 ',
'N-0 ',
5337 &
'P-- ',
'PPPI0 ',
'PNPI+ ',
'PNPI0 ',
'PPPI- ',
5338 &
'NNPI+ ',
'APPPI0 ',
'APNPI+ ',
'ANNPI0 ',
'ANPPI- ',
5339 &
'APNPI0 ',
'APPPI- ',
'ANNPI- ',
'K+PPI0 ',
'K+NPI+ ',
5340 &
'K0PPI0 ',
'K-PPI0 ',
'K-NPI+ ',
'AKPPI- ',
'AKNPI0 ',
5341 &
'K+NPI0 ',
'K+PPI- ',
'K0PPI0 ',
'K0NPI+ ',
'K-NPI0 ',
5342 &
'K-PPI- ',
'AKNPI- ',
'PAK0 ',
'SI+PI0 ',
'SI0PI+ ',
5343 &
'SI+ETA ',
'S+1PI0 ',
'S01PI+ ',
'NEUK- ',
'LA0PI- ',
5344 &
'SI-OM0 ',
'LA0RO- ',
'SI0RO- ',
'SI-RO0 ',
'SI-ET0 ',
5345 &
'SI0PI- ',
'SI-0 ',
'BLANC ',
'BLANC ',
'BLANC ',
5346 &
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
5347 &
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
5348 &
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
5349 &
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
'BLANC ',
5350 &
'EPI+- ',
'EPI00 ',
'GAPI+- ',
'GAGA* ',
'K+- ',
5351 &
'KLKS ',
'PI+-0 ',
'EGA ',
'LPI0 ',
'LPI ' /
5352 DATA (zkname(k),k=341,425) /
5353 &
'APPI0 ',
'ANPI- ',
'ALAGA ',
'ANPI ',
'ALPI0 ',
5354 &
'ALPI+ ',
'LAPI+ ',
'SI+PI0 ',
'SI0PI+ ',
'LAPI0 ',
5355 &
'SI+PI- ',
'SI-PI+ ',
'LAPI- ',
'SI-PI0 ',
'SI0PI- ',
5356 &
'TE0PI0 ',
'TE-PI+ ',
'TE0PI- ',
'TE-PI0 ',
'TE0PI ',
5357 &
'TE-PI ',
'LAK- ',
'ALPI- ',
'AS-PI0 ',
'AS0PI- ',
5358 &
'ALPI0 ',
'AS+PI- ',
'AS-PI+ ',
'ALPI+ ',
'AS+PI0 ',
5359 &
'AS0PI+ ',
'AT0PI0 ',
'AT+PI- ',
'AT0PI+ ',
'AT+PI0 ',
5360 &
'AT0PI ',
'AT+PI ',
'ALK+ ',
'K-PI+ ',
'K-PI+0 ',
5361 &
'K0PI+- ',
'K0PI0 ',
'K-PI++ ',
'AK0PI+ ',
'K+PI-- ',
5362 &
'K0PI- ',
'K+PI- ',
'K+PI-0 ',
'AKPI-+ ',
'AK0PI0 ',
5363 &
'ETAPIF ',
'K++- ',
'K+AK0 ',
'ETAPI- ',
'K--+ ',
5364 &
'K-K0 ',
'PI00 ',
'PI+- ',
'GAGA ',
'D0PI0 ',
5365 &
'D0GA ',
'D0PI+ ',
'D+PI0 ',
'DFGA ',
'AD0PI- ',
5366 &
'D-PI0 ',
'D-GA ',
'AD0PI0 ',
'AD0GA ',
'F+GA ',
5367 &
'F+GA ',
'F-GA ',
'F-GA ',
'PSPI+- ',
'PSPI00 ',
5368 &
'PSETA ',
'E+E- ',
'MUE+- ',
'PI+-0 ',
'M+NN ',
5369 &
'E+NN ',
'RHO+NT ',
'PI+ANT ',
'K*+ANT ',
'M-NN ' /
5370 DATA (zkname(k),k=426,510) /
5371 &
'E-NN ',
'RHO-NT ',
'PI-NT ',
'K*-NT ',
'NUET ',
5372 &
'ANUET ',
'NUEM ',
'ANUEM ',
'SI+ETA ',
'SI+ET* ',
5373 &
'PAK0 ',
'TET0K+ ',
'SI*+ET ',
'N*+AK0 ',
'N*++K- ',
5374 &
'LAMRO+ ',
'SI0RO+ ',
'SI+RO0 ',
'SI+OME ',
'PAK*0 ',
5375 &
'N*+AK* ',
'N*++K* ',
'SI+AK0 ',
'TET0PI ',
'SI+AK* ',
5376 &
'TET0RO ',
'SI0AK* ',
'SI+K*- ',
'TET0OM ',
'TET-RO ',
5377 &
'SI*0AK ',
'C0+PI+ ',
'C0+PI0 ',
'C0+PI- ',
'A+GAM ',
5378 &
'A0GAM ',
'TET0AK ',
'TET0K* ',
'OM-RO+ ',
'OM-PI+ ',
5379 &
'C1++AK ',
'A+PI+ ',
'C0+AK0 ',
'A0PI+ ',
'A+AK0 ',
5380 &
'T0PI+ ',
'ASI-ET ',
'ASI-E* ',
'APK0 ',
'ATET0K ',
5381 &
'ASI*-E ',
'AN*-K0 ',
'AN*--K ',
'ALAMRO ',
'ASI0RO ',
5382 &
'ASI-RO ',
'ASI-OM ',
'APK*0 ',
'AN*-K* ',
'AN*--K ',
5383 &
'ASI-K0 ',
'ATETPI ',
'ASI-K* ',
'ATETRO ',
'ASI0K* ',
5384 &
'ASI-K* ',
'ATE0OM ',
'ATE+RO ',
'ASI*0K ',
'AC-PI- ',
5385 &
'AC-PI0 ',
'AC-PI+ ',
'AA-GAM ',
'AA0GAM ',
'ATET0K ',
5386 &
'ATE0K* ',
'AOM+RO ',
'AOM+PI ',
'AC1--K ',
'AA-PI- ',
5387 &
'AC0-K0 ',
'AA0PI- ',
'AA-K0 ',
'AT0PI- ',
'C1++GA ' /
5388 DATA (zkname(k),k=511,540) /
5389 &
'C1++GA ',
'C10GAM ',
'S+GAM ',
'S0GAM ',
'T0GAM ',
5390 &
'XU++GA ',
'XD+GAM ',
'XS+GAM ',
'A+AKPI ',
'T02PI+ ',
5391 &
'C1++2K ',
'AC1--G ',
'AC1-GA ',
'AC10GA ',
'AS-GAM ',
5392 &
'AS0GAM ',
'AT0GAM ',
'AXU--G ',
'AXD-GA ',
'AXS-GA ',
5393 &
'AA-KPI ',
'AT02PI ',
'AC1--K ',
'RH-PI+ ',
'RH+PI- ',
5394 &
'RH3PI0 ',
'RH0PI+ ',
'RH+PI0 ',
'RH0PI- ',
'RH-PI0 ' /
5395 DATA (zkname(i),i=541,602)/
5396 &
'APETA ',
'AN=P+ ',
'AN-PO ',
'ANOPO ',
'APRHO0',
'ANRHO-',
'ANETA ',
5397 &
'AN-P+ ',
'AN0PO ',
'AN+P- ',
'APRHO+',
'ANRHO0',
'RH0PI+',
'RH+PI0',
5398 &
'3PI+00',
'3PI-++',
'F0PI+ ',
'RH+PI-',
'RH0PI0',
'3PI000',
'3PI0+-',
5399 &
'F0PI0 ',
'RH0PI-',
'RH-PI0',
'3PI-00',
'3PI--+',
'F0PI- ',
'PI0PI0',
5400 &
'PI+PI-',
'K+K- ',
'K0AK0 ',
'L01600',
'AL0160',
'K*+146',
'K*-146',
5401 &
'K*0146',
'AK0146',
'S+1660',
'S01660',
'S-1660',
'AS-166',
'AS0166',
5402 &
'AS+166',
'X01690',
'X-1690',
'AX0169',
'AX+169',
'OM-225',
'AOM+22',
5403 &
'N*PPI0',
'N*NPI+',
'N*P2P0',
'N*PP+-',
'N*D+P0',
'N*D0P+',
'N*NPI0',
5404 &
'N*PPI-',
'N*N2P0',
'N*NP+-',
'N*D+P-',
'N*D0P0',
'BLANK '/
5408 DATA (wt(k),k= 1, 85) /
5409 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5410 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5411 & .1000
d+01, .2100
d+00, .1200
d+00, .2700
d+00, .4000
d+00,
5412 & .1000
d+01, .1000
d+01, .6400
d+00, .2100
d+00, .6000
d-01,
5413 & .2000
d-01, .3000
d-01, .4000
d-01, .6400
d+00, .2100
d+00,
5414 & .6000
d-01, .2000
d-01, .3000
d-01, .4000
d-01, .6400
d+00,
5415 & .3600
d+00, .0000
d+00, .0000
d+00, .6400
d+00, .3600
d+00,
5416 & .0000
d+00, .0000
d+00, .6900
d+00, .3100
d+00, .1000
d+01,
5417 & .5200
d+00, .4800
d+00, .1000
d+01, .9900
d+00, .1000
d-01,
5418 & .3800
d+00, .3000
d-01, .3000
d+00, .2400
d+00, .5000
d-01,
5419 & .1000
d+01, .1000
d+01, .0000
d+00, .1000
d+01, .9000
d+00,
5420 & .1000
d-01, .9000
d-01, .3300
d+00, .6700
d+00, .3300
d+00,
5421 & .6700
d+00, .3300
d+00, .6700
d+00, .3300
d+00, .6700
d+00,
5422 & .3300
d+00, .6700
d+00, .3300
d+00, .6700
d+00, .3300
d+00,
5423 & .6700
d+00, .3300
d+00, .6700
d+00, .1900
d+00, .3800
d+00,
5424 & .9000
d-01, .2000
d+00, .3000
d-01, .4000
d-01, .5000
d-01,
5425 & .2000
d-01, .1900
d+00, .3800
d+00, .9000
d-01, .2000
d+00 /
5426 DATA (wt(k),k= 86,170) /
5427 & .3000
d-01, .4000
d-01, .5000
d-01, .2000
d-01, .1900
d+00,
5428 & .3800
d+00, .9000
d-01, .2000
d+00, .3000
d-01, .4000
d-01,
5429 & .5000
d-01, .2000
d-01, .1900
d+00, .3800
d+00, .9000
d-01,
5430 & .2000
d+00, .3000
d-01, .4000
d-01, .5000
d-01, .2000
d-01,
5431 & .8800
d+00, .6000
d-01, .6000
d-01, .8800
d+00, .6000
d-01,
5432 & .6000
d-01, .8800
d+00, .1200
d+00, .1900
d+00, .1900
d+00,
5433 & .1600
d+00, .1600
d+00, .1700
d+00, .3000
d-01, .3000
d-01,
5434 & .3000
d-01, .4000
d-01, .1000
d+00, .1000
d+00, .2000
d+00,
5435 & .1200
d+00, .1000
d+00, .4000
d-01, .4000
d-01, .5000
d-01,
5436 & .7500
d-01, .7500
d-01, .3000
d-01, .3000
d-01, .4000
d-01,
5437 & .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00, .1000
d+01,
5438 & .6700
d+00, .3300
d+00, .3300
d+00, .6700
d+00, .1000
d+01,
5439 & .2500
d+00, .2700
d+00, .1800
d+00, .3000
d+00, .1700
d+00,
5440 & .8000
d-01, .1800
d+00, .3000
d-01, .2400
d+00, .2000
d+00,
5441 & .1000
d+00, .8000
d-01, .1700
d+00, .2400
d+00, .3000
d-01,
5442 & .1800
d+00, .1000
d+00, .2000
d+00, .2500
d+00, .1800
d+00,
5443 & .2700
d+00, .3000
d+00, .5000
d+00, .3000
d+00, .1250
d+00 /
5447 DATA (wt(k),k=171,255) /
5448 & .7500
d-01, .0000
d+00, .0000
d+00, .5000
d+00, .7500
d-01,
5449 & .1250
d+00, .3000
d+00, .0000
d+00, .0000
d+00, .1800
d+00,
5450 & .3700
d+00, .1300
d+00, .8000
d-01, .4000
d-01, .7000
d-01,
5451 & .1300
d+00, .3700
d+00, .1800
d+00, .4000
d-01, .8000
d-01,
5452 & .1300
d+00, .1300
d+00, .7000
d-01, .7000
d-01, .1300
d+00,
5453 & .2300
d+00, .4700
d+00, .5000
d-01, .2000
d-01, .1000
d-01,
5454 & .2000
d-01, .1300
d+00, .7000
d-01, .4700
d+00, .2300
d+00,
5455 & .5000
d-01, .1000
d-01, .2000
d-01, .2000
d-01, .1000
d+01,
5456 & .6700
d+00, .3300
d+00, .3300
d+00, .6700
d+00, .1000
d+01,
5457 & .2500
d+00, .2700
d+00, .1800
d+00, .3000
d+00, .1700
d+00,
5458 & .8000
d-01, .1800
d+00, .3000
d-01, .2400
d+00, .2000
d+00,
5459 & .1000
d+00, .8000
d-01, .1700
d+00, .2400
d+00, .3000
d-01,
5460 & .1800
d+00, .1000
d+00, .2000
d+00, .2500
d+00, .1800
d+00,
5461 & .2700
d+00, .3000
d+00, .1800
d+00, .3700
d+00, .1300
d+00,
5462 & .8000
d-01, .4000
d-01, .7000
d-01, .1300
d+00, .3700
d+00,
5463 & .1800
d+00, .4000
d-01, .8000
d-01, .1300
d+00, .1300
d+00,
5464 & .7000
d-01, .5000
d+00, .5000
d+00, .1000
d+01, .1000
d+01 /
5465 DATA (wt(k),k=256,340) /
5466 & .1000
d+01, .8000
d+00, .2000
d+00, .6000
d+00, .3000
d+00,
5467 & .1000
d+00, .6000
d+00, .3000
d+00, .1000
d+00, .8000
d+00,
5468 & .2000
d+00, .3300
d+00, .6700
d+00, .6600
d+00, .1700
d+00,
5469 & .1700
d+00, .3200
d+00, .1700
d+00, .3200
d+00, .1900
d+00,
5470 & .3300
d+00, .3300
d+00, .3400
d+00, .3000
d+00, .5000
d-01,
5471 & .6500
d+00, .3800
d+00, .1200
d+00, .3800
d+00, .1200
d+00,
5472 & .3800
d+00, .1200
d+00, .3800
d+00, .1200
d+00, .3000
d+00,
5473 & .5000
d-01, .6500
d+00, .3800
d+00, .2500
d+00, .2500
d+00,
5474 & .2000
d-01, .5000
d-01, .5000
d-01, .2000
d+00, .2000
d+00,
5475 & .1200
d+00, .1000
d+00, .7000
d-01, .7000
d-01, .1400
d+00,
5476 & .5000
d-01, .5000
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 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5479 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5480 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5481 & .4800
d+00, .2400
d+00, .2600
d+00, .2000
d-01, .4700
d+00,
5482 & .3500
d+00, .1500
d+00, .3000
d-01, .1000
d+01, .1000
d+01 /
5483 DATA (wt(k),k=341,425) /
5484 & .5200
d+00, .4800
d+00, .1000
d+01, .1000
d+01, .1000
d+01,
5485 & .1000
d+01, .9000
d+00, .5000
d-01, .5000
d-01, .9000
d+00,
5486 & .5000
d-01, .5000
d-01, .9000
d+00, .5000
d-01, .5000
d-01,
5487 & .3300
d+00, .6700
d+00, .6700
d+00, .3300
d+00, .2500
d+00,
5488 & .2500
d+00, .5000
d+00, .9000
d+00, .5000
d-01, .5000
d-01,
5489 & .9000
d+00, .5000
d-01, .5000
d-01, .9000
d+00, .5000
d-01,
5490 & .5000
d-01, .3300
d+00, .6700
d+00, .6700
d+00, .3300
d+00,
5491 & .2500
d+00, .2500
d+00, .5000
d+00, .1000
d+00, .5000
d+00,
5492 & .1600
d+00, .2400
d+00, .7000
d+00, .3000
d+00, .7000
d+00,
5493 & .3000
d+00, .1000
d+00, .5000
d+00, .1600
d+00, .2400
d+00,
5494 & .3000
d+00, .4000
d+00, .3000
d+00, .3000
d+00, .4000
d+00,
5495 & .3000
d+00, .4900
d+00, .4900
d+00, .2000
d-01, .5500
d+00,
5496 & .4500
d+00, .6800
d+00, .3000
d+00, .2000
d-01, .6800
d+00,
5497 & .3000
d+00, .2000
d-01, .5500
d+00, .4500
d+00, .9000
d+00,
5498 & .1000
d+00, .9000
d+00, .1000
d+00, .6000
d+00, .3000
d+00,
5499 & .1000
d+00, .1000
d+00, .1000
d+00, .8000
d+00, .2800
d+00,
5500 & .2800
d+00, .3500
d+00, .7000
d-01, .2000
d-01, .2800
d+00 /
5501 DATA (wt(k),k=426,510) /
5502 & .2800
d+00, .3500
d+00, .7000
d-01, .2000
d-01, .1000
d+01,
5503 & .1000
d+01, .1000
d+01, .1000
d+01, .2000
d-01, .3000
d-01,
5504 & .7000
d-01, .2000
d-01, .2000
d-01, .4000
d-01, .1300
d+00,
5505 & .7000
d-01, .6000
d-01, .6000
d-01, .2000
d+00, .1400
d+00,
5506 & .4000
d-01, .1000
d+00, .2500
d+00, .3000
d-01, .3000
d+00,
5507 & .4200
d+00, .2200
d+00, .3500
d+00, .1900
d+00, .1600
d+00,
5508 & .8000
d-01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5509 & .1000
d+01, .3700
d+00, .2000
d+00, .3600
d+00, .7000
d-01,
5510 & .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00,
5511 & .5000
d+00, .2000
d-01, .3000
d-01, .7000
d-01, .2000
d-01,
5512 & .2000
d-01, .4000
d-01, .1300
d+00, .7000
d-01, .6000
d-01,
5513 & .6000
d-01, .2000
d+00, .1400
d+00, .4000
d-01, .1000
d+00,
5514 & .2500
d+00, .3000
d-01, .3000
d+00, .4200
d+00, .2200
d+00,
5515 & .3500
d+00, .1900
d+00, .1600
d+00, .8000
d-01, .1000
d+01,
5516 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .3700
d+00,
5517 & .2000
d+00, .3600
d+00, .7000
d-01, .5000
d+00, .5000
d+00,
5518 & .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00, .1000
d+01 /
5519 DATA (wt(k),k=511,540) /
5520 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5521 & .1000
d+01, .1000
d+01, .1000
d+01, .3000
d+00, .3000
d+00,
5522 & .4000
d+00, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5523 & .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01, .1000
d+01,
5524 & .3000
d+00, .3000
d+00, .4000
d+00, .3300
d+00, .3300
d+00,
5525 & .3400
d+00, .5000
d+00, .5000
d+00, .5000
d+00, .5000
d+00 /
5527 DATA (wt(i),i=541,602) / .0
d+00, .3334
d+00, .2083
d+00, 2*.125
d+00,
5528 & .2083
d+00, .0
d+00, .125
d+00, .2083
d+00, .3334
d+00, .2083
d+00,
5529 & .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,
5530 & 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,
5531 & 0.3
d+00, 0.0
d+00, 0.31
d+00, 0.62
d+00, 0.035
d+00, 0.035
d+00,
5532 & 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,
5533 & 0.16
d+00, 2*0.12
d+00, 2*0.05
d+00, 1.
d+00 /
5537 DATA (nzk(k,1),k= 1,170) /
5538 & 1, 2, 3, 4, 5, 6, 7, 1, 2, 4,
5539 & 3, 23, 13, 13, 13, 10, 11, 10, 13, 13,
5540 & 13, 10, 4, 11, 14, 14, 14, 11, 3, 1,
5541 & 8, 1, 1, 2, 9, 2, 2, 13, 23, 8,
5542 & 1, 8, 17, 7, 7, 7, 23, 23, 13, 13,
5543 & 13, 13, 23, 14, 13, 13, 23, 15, 24, 24,
5544 & 15, 16, 25, 25, 16, 15, 24, 24, 15, 16,
5545 & 24, 25, 16, 15, 24, 36, 37, 15, 24, 15,
5546 & 15, 24, 15, 37, 36, 24, 15, 24, 24, 16,
5547 & 24, 38, 39, 16, 25, 16, 16, 25, 16, 39,
5548 & 38, 25, 16, 25, 25, 17, 22, 21, 17, 21,
5549 & 20, 17, 22, 8, 1, 21, 22, 20, 17, 48,
5550 & 50, 49, 8, 1, 17, 17, 17, 21, 20, 22,
5551 & 17, 22, 21, 20, 22, 19, 12, 19, 12, 1,
5552 & 1, 8, 1, 8, 8, 1, 53, 54, 1, 1,
5553 & 8, 53, 54, 55, 1, 8, 1, 8, 54, 55,
5554 & 56, 1, 8, 8, 55, 56, 8, 1, 53, 54 /
5555 DATA (nzk(k,1),k=171,340) /
5556 & 55, 1, 8, 8, 54, 55, 56, 1, 8, 1,
5557 & 8, 53, 54, 55, 1, 8, 1, 8, 54, 55,
5558 & 56, 1, 8, 1, 8, 1, 8, 17, 21, 22,
5559 & 1, 1, 8, 1, 8, 17, 22, 20, 8, 2,
5560 & 2, 9, 2, 9, 9, 2, 67, 68, 2, 2,
5561 & 9, 67, 68, 69, 2, 9, 2, 9, 68, 69,
5562 & 70, 2, 9, 9, 69, 70, 9, 2, 9, 67,
5563 & 68, 69, 2, 9, 2, 9, 68, 69, 70, 2,
5564 & 9, 1, 8, 13, 13, 14, 1, 8, 1, 1,
5565 & 8, 8, 8, 1, 8, 1, 1, 1, 1, 1,
5566 & 8, 2, 2, 9, 9, 2, 2, 9, 15, 15,
5567 & 24, 16, 16, 25, 25, 15, 15, 24, 24, 16,
5568 & 16, 25, 1, 21, 22, 21, 48, 49, 8, 17,
5569 & 20, 17, 22, 20, 20, 22, 20, 0, 0, 0,
5570 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5571 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5572 & 31, 31, 13, 7, 15, 12, 13, 31, 17, 17 /
5573 DATA (nzk(k,1),k=341,510) /
5574 & 2, 9, 18, 9, 18, 18, 17, 21, 22, 17,
5575 & 21, 20, 17, 20, 22, 97, 98, 97, 98, 97,
5576 & 98, 17, 18, 99, 100, 18, 101, 99, 18, 101,
5577 & 100, 102, 103, 102, 103, 102, 103, 18, 16, 16,
5578 & 24, 24, 16, 25, 15, 24, 15, 15, 25, 25,
5579 & 31, 15, 15, 31, 16, 16, 23, 13, 7, 116,
5580 & 116, 116, 117, 117, 119, 118, 118, 119, 119, 120,
5581 & 120, 121, 121, 130, 130, 130, 4, 10, 13, 10,
5582 & 4, 32, 13, 36, 11, 3, 34, 14, 38, 133,
5583 & 134, 135, 136, 21, 21, 1, 97, 104, 54, 53,
5584 & 17, 22, 21, 21, 1, 54, 53, 21, 97, 21,
5585 & 97, 22, 21, 97, 98, 105, 137, 137, 137, 138,
5586 & 139, 97, 97, 109, 109, 140, 138, 137, 139, 138,
5587 & 145, 99, 99, 2, 102, 110, 68, 67, 18, 100,
5588 & 99, 99, 2, 68, 67, 99, 102, 99, 102, 100,
5589 & 99, 102, 103, 111, 149, 149, 149, 150, 151, 113,
5590 & 113, 115, 115, 152, 150, 149, 151, 150, 157, 140 /
5591 DATA (nzk(k,1),k=511,540) /
5592 & 141, 142, 143, 144, 145, 146, 147, 148, 138, 145,
5593 & 140, 152, 153, 154, 155, 156, 157, 158, 159, 160,
5594 & 150, 157, 152, 34, 32, 33, 33, 32, 33, 34 /
5595 DATA (nzk(i,1),i=541,602) / 2, 67, 68, 69, 2, 9, 9, 68, 69,
5596 & 70, 2, 9, 33, 32, 13, 14, 189, 32, 34, 23, 23, 189, 33, 34, 14,
5597 & 14, 189, 23, 13, 15, 24, 36, 38, 37, 39, 194, 195, 196, 197,
5598 & 198, 199, 200, 201, 202, 203, 204, 205, 206, 207, 1, 8, 1, 1, 54,
5599 & 55, 8, 1, 8, 8, 54, 55, 210/
5600 DATA (nzk(k,2),k= 1,170) /
5601 & 0, 0, 0, 0, 0, 0, 0, 3, 4, 6,
5602 & 5, 23, 14, 11, 3, 5, 5, 5, 23, 13,
5603 & 23, 23, 23, 5, 23, 13, 23, 23, 23, 14,
5604 & 23, 3, 11, 13, 23, 4, 10, 14, 23, 14,
5605 & 23, 13, 7, 7, 4, 7, 7, 23, 14, 14,
5606 & 23, 14, 23, 23, 14, 14, 7, 23, 13, 23,
5607 & 14, 23, 14, 23, 13, 23, 13, 23, 14, 23,
5608 & 14, 23, 13, 23, 13, 23, 13, 33, 32, 35,
5609 & 31, 23, 14, 23, 14, 33, 34, 35, 31, 23,
5610 & 14, 23, 14, 33, 34, 35, 31, 23, 13, 23,
5611 & 13, 33, 32, 35, 31, 13, 13, 23, 23, 14,
5612 & 13, 14, 14, 25, 16, 14, 23, 13, 31, 14,
5613 & 13, 23, 25, 16, 23, 35, 33, 34, 32, 33,
5614 & 31, 31, 14, 13, 23, 0, 0, 0, 0, 13,
5615 & 23, 13, 14, 23, 14, 13, 23, 13, 78, 23,
5616 & 13, 14, 23, 13, 79, 78, 14, 23, 14, 23,
5617 & 13, 80, 79, 14, 14, 23, 80, 31, 14, 23 /
5618 DATA (nzk(k,2),k=171,340) /
5619 & 13, 79, 78, 31, 14, 23, 13, 80, 79, 23,
5620 & 13, 14, 23, 13, 79, 78, 14, 23, 14, 23,
5621 & 13, 80, 79, 23, 13, 33, 32, 15, 24, 15,
5622 & 31, 14, 23, 34, 33, 24, 24, 15, 31, 14,
5623 & 23, 14, 13, 23, 13, 14, 23, 14, 80, 23,
5624 & 14, 13, 23, 14, 79, 80, 13, 23, 13, 23,
5625 & 14, 78, 79, 13, 13, 23, 78, 23, 14, 13,
5626 & 23, 14, 79, 80, 13, 23, 13, 23, 14, 78,
5627 & 79, 62, 61, 23, 14, 23, 13, 13, 13, 23,
5628 & 13, 13, 23, 14, 14, 14, 1, 8, 8, 1,
5629 & 8, 1, 8, 8, 1, 8, 1, 8, 1, 8,
5630 & 1, 1, 8, 1, 8, 8, 1, 1, 8, 8,
5631 & 1, 8, 25, 23, 13, 31, 23, 13, 16, 14,
5632 & 35, 34, 34, 33, 31, 14, 23, 0, 0, 0,
5633 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5634 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5635 & 13, 23, 14, 7, 16, 19, 14, 7, 23, 14 /
5636 DATA (nzk(k,2),k=341,510) /
5637 & 23, 14, 7, 13, 23, 13, 13, 23, 13, 23,
5638 & 14, 13, 14, 23, 14, 23, 13, 14, 23, 14,
5639 & 23, 16, 14, 23, 14, 23, 14, 13, 13, 23,
5640 & 13, 23, 14, 13, 23, 13, 23, 15, 13, 13,
5641 & 13, 23, 13, 13, 14, 14, 14, 14, 14, 23,
5642 & 13, 16, 25, 14, 15, 24, 23, 14, 7, 23,
5643 & 7, 13, 23, 7, 14, 23, 7, 23, 7, 7,
5644 & 7, 7, 7, 13, 23, 31, 3, 11, 14, 135,
5645 & 5, 134, 134, 134, 136, 6, 133, 133, 133, 0,
5646 & 0, 0, 0, 31, 95, 25, 15, 31, 95, 16,
5647 & 32, 32, 33, 35, 39, 39, 38, 25, 13, 39,
5648 & 32, 39, 38, 35, 32, 39, 13, 23, 14, 7,
5649 & 7, 25, 37, 32, 13, 25, 13, 25, 13, 25,
5650 & 13, 31, 95, 24, 16, 31, 24, 15, 34, 34,
5651 & 33, 35, 37, 37, 36, 24, 14, 37, 34, 37,
5652 & 36, 35, 34, 37, 14, 23, 13, 7, 7, 24,
5653 & 39, 34, 14, 24, 14, 24, 14, 24, 14, 7 /
5654 DATA (nzk(k,2),k=511,540) /
5655 & 7, 7, 7, 7, 7, 7, 7, 7, 25, 13,
5656 & 25, 7, 7, 7, 7, 7, 7, 7, 7, 7,
5657 & 24, 14, 24, 13, 14, 23, 13, 23, 14, 23 /
5658 DATA (nzk(i,2),i=541,602) / 31, 13, 23, 14, 79, 80, 31, 13, 23,
5659 & 14, 78, 79, 13, 23, 23, 13, 13, 14, 13, 23, 13, 23, 14, 23, 23,
5660 & 14, 14, 23, 14, 16, 25,
5661 & 4*23, 14*0, 23, 13, 23, 13, 23, 13, 23, 14,
5662 & 23, 13, 14, 23, 0 /
5663 DATA (nzk(k,3),k= 1,170) /
5664 & 0, 0, 0, 0, 0, 0, 0, 5, 6, 5,
5665 & 6, 23, 23, 5, 5, 0, 0, 0, 0, 14,
5666 & 23, 5, 5, 0, 0, 14, 23, 5, 5, 0,
5667 & 0, 5, 5, 0, 0, 5, 5, 0, 0, 0,
5668 & 0, 0, 0, 0, 3, 0, 7, 23, 23, 7,
5669 & 0, 0, 0, 0, 23, 0, 0, 0, 0, 0,
5671 DATA (nzk(k,3),k=171,340) /
5673 & 0, 0, 0, 0, 0, 0, 23, 13, 14, 23,
5674 & 23, 14, 23, 23, 23, 14, 23, 13, 23, 14,
5675 & 13, 23, 13, 23, 14, 23, 14, 14, 23, 13,
5676 & 13, 23, 13, 14, 23, 23, 14, 23, 13, 23,
5677 & 14, 14, 0, 0, 0, 0, 0, 0, 0, 0,
5679 & 14, 23, 7, 0, 0, 0, 23, 0, 0, 0 /
5680 DATA (nzk(k,3),k=341,510) /
5682 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 23,
5683 & 14, 0, 13, 0, 14, 0, 0, 23, 13, 0,
5684 & 0, 15, 0, 0, 16, 0, 0, 0, 0, 0,
5685 & 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5686 & 0, 0, 0, 14, 23, 0, 0, 0, 23, 134,
5687 & 134, 0, 0, 0, 133, 133, 0, 0, 0, 0,
5689 DATA (nzk(k,3),k=511,540) /
5690 & 0, 0, 0, 0, 0, 0, 0, 0, 13, 13,
5691 & 25, 0, 0, 0, 0, 0, 0, 0, 0, 0,
5692 & 14, 14, 24, 0, 0, 0, 0, 0, 0, 0 /
5693 DATA (nzk(i,3),i=541,602) / 12*0, 2*0, 23, 13, 0, 2*0, 23, 14, 0,
5694 & 2*0, 23, 13, 0, 4*0, 18*0, 2*0, 23, 14, 2*0, 2*0, 23, 14, 2*0, 0/
5715 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
5717 parameter(lout=6,llook=9)
5719 COMPLEX*16 czero,cone,ctwo
5720 parameter(
zero=0.0d0,
one=1.0d0,two=2.0d0,three=3.0d0,
5721 & onethi=
one/three,tiny25=1.0
d-25)
5722 parameter(twopi = 6.283185307179586454
d+00,
5724 & gev2mb = 0.38938d0,
5725 & gev2fm = 0.1972d0,
5726 & alphem =
one/137.0d0,
5736 COMMON /dpar/ aname(210),aam(210),gam(210),
tau(210),iich(210),
5737 & iibar(210),ka1(210),ka2(210)
5739 parameter(ncompx=1,
neb=50)
5740 COMMON /dshmm/ rash,rbsh(ncompx),bmax(ncompx),bstep(ncompx),
5741 & sigsh,rosh,gsh,bsite(0:
neb,ncompx,ksiteb),
5743 COMMON /glaber/ ecmnn(
neb),ecmnow,
5752 COMMON /vdmpar/ rl2,epspol,intrge(2),idpdf,modega,ishad(3)
5753 COMMON /glapar/ jstatb
5756 COMMON /damp/ ca,ci,ga
5757 COMMON /xsecnu/ecmuu,ecmoo,ngritt,nevtt
5758 COMMON /kglaub/jglaub
5760 parameter(maxncl = 210)
5761 COMPLEX*16 pp11,pp12,pp21,pp22,
5762 & ompp11,ompp12,ompp21,ompp22
5763 dimension coop1(3,maxncl),coot1(3,maxncl),
5764 & coop2(3,maxncl),coot2(3,maxncl),
5765 & bprod(ksiteb),sigshh(
neb),
5769 Write(6,*)
' XSGLAU(NA,NB,IJPROJ,NTARG)',
5771 WRITE(6,*)
'/XSECNU/ECMUU,ECMOO,NGRITT,NEVTT',
5772 &ecmuu,ecmoo,ngritt,nevtt
5776 ctwo = dcmplx(two,
zero)
5780 dellog=(log10(ecmoo)-log10(ecmuu))/(ngritt-1)
5781 deldel=10.d0**dellog
5783 DO 1123 ieee=1,ngritt
5788 WRITE(6,*)
'IE,EC111,S',ie,ec111,
s
5797 rash = rnucle*dble(na)**onethi
5798 rbsh(
ntarg) = rnucle*dble(nb)**onethi
5800 IF(na.EQ.9)rash=2.52d0
5801 IF(na.EQ.10)rash=2.45d0
5802 IF(na.EQ.11)rash=2.37d0
5803 IF(na.EQ.12)rash=2.45d0
5804 IF(na.EQ.13)rash=2.44d0
5805 IF(na.EQ.14)rash=2.55d0
5806 IF(na.EQ.15)rash=2.58d0
5807 IF(na.EQ.16)rash=2.71d0
5808 IF(na.EQ.17)rash=2.66d0
5809 IF(na.EQ.18)rash=2.71d0
5810 IF(nb.EQ.9)rbsh(
ntarg)=2.52d0
5811 IF(nb.EQ.10)rbsh(
ntarg)=2.45d0
5812 IF(nb.EQ.11)rbsh(
ntarg)=2.37d0
5813 IF(nb.EQ.12)rbsh(
ntarg)=2.45d0
5814 IF(nb.EQ.13)rbsh(
ntarg)=2.44d0
5815 IF(nb.EQ.14)rbsh(
ntarg)=2.55d0
5816 IF(nb.EQ.15)rbsh(
ntarg)=2.58d0
5817 IF(nb.EQ.16)rbsh(
ntarg)=2.71d0
5818 IF(nb.EQ.17)rbsh(
ntarg)=2.66d0
5819 IF(nb.EQ.18)rbsh(
ntarg)=2.71d0
5826 IF (ijproj.LE.12)
THEN
5827 bslope = 8.5d0*(1.0d0+0.065d0*
log(
s))
5828 IF (ecmnn(ie).LE.3.0d0)
THEN
5830 ELSEIF ((ecmnn(ie).GT.3.0d0).AND.(ecmnn(ie).LE.50.d0))
THEN
5831 rosh = -0.63d0+0.175d0*
log(ecmnn(ie))
5832 ELSEIF (ecmnn(ie).GT.50.0d0)
THEN
5836 bslope = 6.0d0*(1.0d0+0.065d0*
log(
s))
5841 elab = (
s-aam(ijproj)**2-amp2)/(two*amp)
5842 elabb(ie)=elab/1000.
5843 plab =
sqrt( (elab-aam(ijproj))*(elab+aam(ijproj)) )
5845 sigsh =
dshpto(ijproj,plab)/10.d0
5846 sigshh(ie)=sigsh*10.d0
5847 WRITE(6,*)
' NSTATB,NSITEB,RASH,RBSH(NTARG),BMAX(NTARG),
5848 &BSLOPE,ROSH,SIGSH,ECM ELAB',
5849 & nstatb,nsiteb,rash,rbsh(
ntarg),bmax(
ntarg),
5850 &bslope,rosh,sigsh,ec111,elab
5869 facn =
one/dble(nstatb)
5898 b = dble(ib)*bstep(
ntarg)
5899 facb = 10.0d0*twopi*b*bstep(
ntarg)
5913 gsh = 10.0d0/(two*bslope*gev2mb)
5916 rca = ga*sigsh/twopi
5918 ca = dcmplx(rca,fca)
5926 x11 = b+coot1(1,inb)-coop1(1,ina)
5927 y11 = coot1(2,inb)-coop1(2,ina)
5928 xy11 = ga*(x11*x11+y11*y11)
5929 x12 = b+coot2(1,inb)-coop1(1,ina)
5930 y12 = coot2(2,inb)-coop1(2,ina)
5931 xy12 = ga*(x12*x12+y12*y12)
5932 x21 = b+coot1(1,inb)-coop2(1,ina)
5933 y21 = coot1(2,inb)-coop2(2,ina)
5934 xy21 = ga*(x21*x21+y21*y21)
5935 x22 = b+coot2(1,inb)-coop2(1,ina)
5936 y22 = coot2(2,inb)-coop2(2,ina)
5937 xy22 = ga*(x22*x22+y22*y22)
5938 IF (xy11.LE.15.0d0)
THEN
5939 c = cone-ca*
exp(-xy11)
5942 IF (abs(ar).LT.tiny25) ar =
zero
5943 IF (abs(ai).LT.tiny25) ai =
zero
5944 pp11 = dcmplx(ar,ai)
5948 shi = shi+
log(ar*ar+ai*ai)
5950 IF (xy12.LE.15.0d0)
THEN
5951 c = cone-ca*
exp(-xy12)
5954 IF (abs(ar).LT.tiny25) ar =
zero
5955 IF (abs(ai).LT.tiny25) ai =
zero
5956 pp12 = dcmplx(ar,ai)
5959 IF (xy21.LE.15.0d0)
THEN
5960 c = cone-ca*
exp(-xy21)
5963 IF (abs(ar).LT.tiny25) ar =
zero
5964 IF (abs(ai).LT.tiny25) ai =
zero
5965 pp21 = dcmplx(ar,ai)
5968 IF (xy22.LE.15.0d0)
THEN
5969 c = cone-ca*
exp(-xy22)
5972 IF (abs(ar).LT.tiny25) ar =
zero
5973 IF (abs(ai).LT.tiny25) ai =
zero
5974 pp22 = dcmplx(ar,ai)
5982 ompp11 = ompp11+(cone-pp11)
5983 ompp21 = ompp21+(cone-pp21)
5986 ompp12 = ompp12+(cone-pp12)
5987 ompp22 = ompp22+(cone-pp22)
5989 stotm = dble(ompp11+ompp22)
5990 selam = dble(ompp11*dconjg(ompp22))
5992 sqepm = dble(ompp11*dconjg(ompp21))-selam
5993 sqetm = dble(ompp11*dconjg(ompp12))-selam
5994 sqe2m = dble(ompp11*dconjg(ompp11))-selam-sqepm-sqetm
5996 stotb = stotb+facm*stotm
5997 selab = selab+facm*selam
5998 IF (nb.GT.1) sqepb = sqepb+facm*sqepm
5999 IF (na.GT.1) sqetb = sqetb+facm*sqetm
6000 IF ((na.GT.1).AND.(nb.GT.1)) sqe2b = sqe2b+facm*sqe2m
6001 sprob = sprob+facm*sprom
6005 stotn = stotn+facb*stotb
6006 selan = selan+facb*selab
6007 sqepn = sqepn+facb*sqepb
6008 sqetn = sqetn+facb*sqetb
6009 sqe2n = sqe2n+facb*sqe2b
6010 spron = spron+facb*sprob
6011 bprod(ib+1)= bprod(ib+1)+facn*facb*sprob
6015 stot = stot +facn*stotn
6016 stot2 = stot2+facn*stotn**2
6017 sela = sela +facn*selan
6018 sela2 = sela2+facn*selan**2
6019 sqep = sqep +facn*sqepn
6020 sqep2 = sqep2+facn*sqepn**2
6021 sqet = sqet +facn*sqetn
6022 sqet2 = sqet2+facn*sqetn**2
6023 sqe2 = sqe2 +facn*sqe2n
6024 sqe22 = sqe22+facn*sqe2n**2
6025 spro = spro +facn*spron
6026 spro2 = spro2+facn*spron**2
6043 WRITE(6,*)
' STOT,SELA ,SQEP,SQET,SQE2,SPRO ',
6044 & stot,sela ,sqep,sqet,sqe2,spro
6046 xetot(ie) =
sqrt(abs(stot2-stot**2)/dble(nstatb-1))
6047 xeela(ie) =
sqrt(abs(sela2-sela**2)/dble(nstatb-1))
6048 xeqep(ie) =
sqrt(abs(sqep2-sqep**2)/dble(nstatb-1))
6049 xeqet(ie) =
sqrt(abs(sqet2-sqet**2)/dble(nstatb-1))
6050 xeqe2(ie) =
sqrt(abs(sqe22-sqe2**2)/dble(nstatb-1))
6051 xepro(ie) =
sqrt(abs(spro2-spro**2)/dble(nstatb-1))
6052 WRITE(6,*)
' XETOT,XEELA,XEQEP,XEQET,XEQE2,XEPRO ',
6053 & xetot(ie),xeela(ie),xeqep(ie),
6054 &xeqet(ie),xeqe2(ie),xepro(ie)
6057 bsite(ie,
ntarg,i) = bprod(i)/spro+bsite(ie,
ntarg,i-1)
6059 & bsite(0,
ntarg,i) = bprod(i)/spro+bsite(0,
ntarg,i-1)
6061 WRITE (6,*)
' ECMNN,ELABB,SIGSHH,SIGTO,SIGEL,SIGIN,SIGSD'
6064 sigto(i)=
dshnto(1,1,ecmnn(i))
6065 sigel(i)=
dshnel(1,1,ecmnn(i))
6066 sigin(i)=
siinel(1,1,ecmnn(i))
6067 sigsd(i)=
sippsd(ecmnn(i))
6068 CALL
sihndi(ecmnn(i),1,1,sigdif(i),sigdih)
6069 WRITE (6,
'(2F18.4,6F11.3)')ecmnn(i),elabb(i),sigshh(i),
6070 & sigto(i),sigel(i),sigin(i),sigsd(i),sigdif(i)
6072 WRITE (6,*)
' ECMNN,ELABB,XSQEP,XEQEP,XSQET,XEQET,XSQE2,XEQE2'
6074 WRITE (6,
'(2F18.4,6F11.3)')ecmnn(i),elabb(i),xsqep(i),xeqep(i),
6075 * xsqet(i),xeqet(i),xsqe2(i),xeqe2(i)
6077 WRITE (6,*)
' ECMNN,ELABB,XSTOT,XETOT,XSELA,XEELA,XSPRO,XEPRO'
6079 WRITE (6,
'(2F18.4,6F11.3)')ecmnn(i),elabb(i),xstot(i),xetot(i),
6080 * xsela(i),xeela(i),xspro(i),xepro(i)
6088 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6090 parameter(maxncl = 210)
6091 dimension coop1(3,maxncl)
6092 CALL
conucl(coop1,na,rash)
6100 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6106 COMMON /dinpda/ imps(6,6),imve(6,6),ib08(6,21),ib10(6,21), ia08
6107 +(6,21),ia10(6,21), a1,b1,b2,b3,lt,le,bet,as,b8,ame,diq,isu
6131 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6134 dimension ec(30),sd(30)
6135 DATA ec /0.d0, 5.d0, 20.d0, 50.d0, 100.d0,
6136 * 200.d0, 500.d0, 1000.d0, 1500.d0, 2000.d0,
6137 * 3000.d0, 4000.d0, 6000.d0, 8000.d0, 10000.d0,
6138 * 15000.d0, 20000.d0, 30000.d0, 40000.d0, 60000.d0,
6139 * 80000.d0, 100000.d0, 150000.d0, 200000.d0, 300000.d0,
6140 * 400000.d0, 600000.d0, 800000.d0, 1000000.d0, 2000000.d0/
6141 DATA sd /0.d0, 0.d0, 5.00d0, 6.14d0, 6.93d0,
6142 * 7.64d0, 8.43d0, 8.87d0, 9.07d0, 9.17d0,
6143 * 9.33d0, 9.40d0, 9.49d0, 9.56d0, 9.58d0,
6144 * 9.69d0, 9.72d0, 9.82d0, 9.85d0, 9.97d0,
6145 * 10.02d0, 10.03d0, 10.13d0, 10.16d0, 10.25d0,
6146 * 10.28d0, 10.39d0, 10.42d0, 10.43d0, 10.53d0/
6149 IF((ecm.GE.ec(i)).AND.(ecm.LT.ec(i+1)))
THEN
6151 del=(ecm-ec(i))*(sd(i+1)-sd(i))/(ec(i+1)-ec(i))
6159 DOUBLE PRECISION FUNCTION siinel(KPROJ,KTARG,UMO)
6160 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6173 DOUBLE PRECISION FUNCTION phnsch ( KP, KTARG, PLAB )
6250 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
6251 parameter( kalgnm = 2 )
6252 parameter( anglgb = 5.0
d-16 )
6253 parameter( anglsq = 2.5
d-31 )
6254 parameter( axcssv = 0.2
d+16 )
6255 parameter( andrfl = 1.0
d-38 )
6256 parameter( avrflw = 1.0
d+38 )
6257 parameter( ainfnt = 1.0
d+30 )
6258 parameter( azrzrz = 1.0
d-30 )
6259 parameter( einfnt = +69.07755278982137
d+00 )
6260 parameter( ezrzrz = -69.07755278982137
d+00 )
6261 parameter( onemns = 0.999999999999999
d+00 )
6262 parameter( onepls = 1.000000000000001
d+00 )
6263 parameter( csnnrm = 2.0
d-15 )
6264 parameter( dmxtrn = 1.0
d+08 )
6296 parameter( zerzer = 0.
d+00 )
6297 parameter( oneone = 1.
d+00 )
6298 parameter( twotwo = 2.
d+00 )
6299 parameter( thrthr = 3.
d+00 )
6300 parameter( foufou = 4.
d+00 )
6301 parameter( fivfiv = 5.
d+00 )
6302 parameter( sixsix = 6.
d+00 )
6303 parameter( sevsev = 7.
d+00 )
6304 parameter( eigeig = 8.
d+00 )
6305 parameter( aninen = 9.
d+00 )
6306 parameter( tenten = 10.
d+00 )
6307 parameter( hlfhlf = 0.5
d+00 )
6308 parameter( onethi = oneone / thrthr )
6309 parameter( twothi = twotwo / thrthr )
6310 parameter( pipipi = 3.1415926535897932270
d+00 )
6311 parameter( eneper = 2.7182818284590452354
d+00 )
6312 parameter( sqrent = 1.6487212707001281468
d+00 )
6354 parameter( clight = 2.99792458
d+10 )
6355 parameter( avogad = 6.0221367
d+23 )
6356 parameter( amelgr = 9.1093897
d-28 )
6357 parameter( plckbr = 1.05457266
d-27 )
6358 parameter( elccgs = 4.8032068
d-10 )
6359 parameter( elcmks = 1.60217733
d-19 )
6360 parameter( amugrm = 1.6605402
d-24 )
6361 parameter( ammumu = 0.113428913
d+00 )
6376 parameter( alpfsc = 7.2973530791728595
d-03 )
6377 parameter( fscto2 = 5.3251361962113614
d-05 )
6378 parameter( fscto3 = 3.8859399018437826
d-07 )
6379 parameter( fscto4 = 2.8357075508200407
d-09 )
6380 parameter( plabrc = 0.197327053
d+00 )
6381 parameter( amelct = 0.51099906
d-03 )
6382 parameter( amugev = 0.93149432
d+00 )
6383 parameter( ammuon = 0.105658389
d+00 )
6384 parameter( rclsel = 2.8179409183694872
d-13 )
6385 parameter( gevmev = 1.0
d+03 )
6386 parameter( emvgev = 1.0
d-03 )
6387 parameter( raddeg = 180.
d+00 / pipipi )
6388 parameter( degrad = pipipi / 180.
d+00 )
6418 parameter( mxxrgn = 500 )
6419 parameter( mxxmdf = 56 )
6420 parameter( mxxmde = 50 )
6421 parameter( mfstck = 1000 )
6422 parameter( mestck = 100 )
6423 parameter( nallwp = 39 )
6424 parameter( mpdpdx = 8 )
6425 parameter( icomax = 180 )
6426 parameter( nstbis = 304 )
6427 parameter( idmaxp = 210 )
6456 parameter( lunin = 5 )
6457 parameter( lunout = 6 )
6458 parameter( lunerr = 66 )
6459 parameter( lunber = 14 )
6460 parameter( lunech = 8 )
6461 parameter( lunflu = 86 )
6462 parameter( lungeo = 16 )
6463 parameter( lunpgs = 12 )
6464 parameter( lunran = 2 )
6465 parameter( lunxsc = 81 )
6466 parameter( lunrdb = 1 )
6511 dimension ichrge(39),am(39)
6549 COMMON /dpar/ aname(210),aam(210),ga(210),
tau(210),
6550 + iich(210),iibar(210),k1(210),k2(210)
6551 dimension kptoip(210),iptokp(39)
6567 DATA kptoip / 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
6568 & 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 66*0,
6569 & 34, 36, 31, 32, 33, 35, 37, 5*0, 38, 5*0, 39, 19*0, 27, 28, 74*0/
6573 DATA iptokp / 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15,
6574 & 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 99,
6575 & 100, 101, 97, 102, 98, 103, 109, 115 /
6625 COMMON / qquark / iqechr(-6:6), iqbchr(-6:6), iqichr(-6:6),
6626 & iqschr(-6:6), iqcchr(-6:6), iquchr(-6:6),
6627 & iqtchr(-6:6), mquark(3,39)
6629 dimension sieapp(11), sitapp(16), plaetb(16)
6630 dimension sgtcoe(5,33), plalim(2,33), ihlp(nallwp)
6631 dimension sgtco1(5,10),sgtco2(5,8),sgtco3(5,15)
6632 SAVE plaetb, sieapp, sitapp, sgtcoe, plalim, ihlp
6633 SAVE iqfsc1, iqfsc2, iqbsc1, iqbsc2
6634 equivalence(sgtco1(1,1),sgtcoe(1,1))
6635 equivalence(sgtco2(1,1),sgtcoe(1,11))
6636 equivalence(sgtco3(1,1),sgtcoe(1,19))
6638 DATA ihlp/1,4,5*0,1,4,2*0,3,2*2,2*3,1,4,3,3*1,2,
6639 & 2*3, 2, 4*0, 3*4, 1, 4, 1, 4, 1, 4 /
6643 &0.147
d+00, zerzer , zerzer , 0.0022
d+00, -0.0170
d+00,
6645 &0.300
d+00, zerzer , zerzer , 0.0095
d+00, -0.057
d+00,
6647 &16.4
d+00, 19.3
d+00, -0.42
d+00, 0.19
d+00, zerzer ,
6649 &33.0
d+00, 14.0
d+00, -1.36
d+00, 0.456
d+00, -4.03
d+00,
6651 &56.8
d+00, 42.2
d+00, -1.45
d+00, 0.65
d+00, -5.39
d+00,
6653 &18.1
d+00, zerzer , zerzer , 0.26
d+00, -1.0
d+00,
6655 &18.7
d+00, zerzer , zerzer , 0.21
d+00, -0.89
d+00,
6657 &34.2
d+00, 7.9
d+00, -2.1
d+00, 0.346
d+00, -0.99
d+00,
6659 &32.1
d+00, zerzer , zerzer , 0.66
d+00, -5.6
d+00,
6661 &25.2
d+00, zerzer , zerzer , 0.38
d+00, -2.9
d+00/
6665 &57.6
d+00, zerzer , zerzer , 1.17
d+00, -9.5
d+00,
6667 &48.0
d+00, zerzer , zerzer , 0.522
d+00, -4.51
d+00,
6669 &47.30
d+00, zerzer , zerzer , 0.513
d+00, -4.27
d+00,
6671 &91.3
d+00, zerzer , zerzer , 1.05
d+00, -8.8
d+00,
6673 &38.4
d+00, 77.6
d+00, -0.64
d+00, 0.26
d+00, -1.2
d+00,
6675 &zerzer ,133.6
d+00, -0.70
d+00, -1.22
d+00, 13.7
d+00,
6677 &112.
d+00, 125.
d+00, -1.08
d+00, 1.14
d+00, -12.4
d+00,
6679 &30.4
d+00, zerzer , zerzer , zerzer , 1.6
d+00/
6683 &zerzer , 11.4
d+00, -0.4
d+00, 0.079
d+00, zerzer ,
6685 &1.76
d+00, 11.2
d+00, -0.64
d+00, 0.043
d+00, zerzer ,
6687 &5.0
d+00, 8.1
d+00, -1.8
d+00, 0.16
d+00, -1.3
d+00,
6689 &7.3
d+00, zerzer , zerzer , 0.29
d+00, -2.40
d+00,
6691 &11.9
d+00, 26.9
d+00, -1.21
d+00, 0.169
d+00, -1.85
d+00,
6693 &16.1
d+00, zerzer , zerzer , 0.32
d+00, -3.4
d+00,
6695 &10.2
d+00, 52.7
d+00, -1.16
d+00, 0.125
d+00, -1.28
d+00,
6697 &10.6
d+00, 53.1
d+00, -1.19
d+00, 0.136
d+00, -1.41
d+00,
6699 &36.5
d+00, zerzer , zerzer , zerzer , -11.9
d+00,
6701 &12.3
d+00, zerzer , zerzer , zerzer , -2.4
d+00,
6703 &7.24
d+00, 46.0
d+00, -4.71
d+00, 0.279
d+00, -2.35
d+00,
6705 &zerzer ,0.912
d+00, -1.22
d+00, zerzer , zerzer ,
6707 &zerzer , 3.39
d+00, -1.75
d+00, zerzer , zerzer ,
6709 &zerzer , 7.18
d+00, -2.01
d+00, zerzer , zerzer ,
6711 &zerzer , 18.8
d+00, -2.01
d+00, zerzer , zerzer /
6715 & 3.0
d+00, 183.
d+00, 2.0
d+00, 17.8
d+00, 4.0
d+00, 340.
d+00,
6717 & 2.5
d+00, 370.
d+00, 2.5
d+00, 370.
d+00, 2.0
d+00, 310.
d+00,
6719 & 2.0
d+00, 310.
d+00, 2.0
d+00, 310.
d+00, 3.0
d+00, 310.
d+00,
6721 & 1.8
d+00, 310.
d+00, 3.0
d+00, 310.
d+00, 3.0
d+00, 2100.
d+00,
6723 & 3.0
d+00, 370.
d+00, 3.0
d+00, 370.
d+00, 5.0
d+00, 1.73
d+06,
6725 & 1.1
d+00, 280.
d+00, 2.0
d+00, 280.
d+00, 0.6
d+00, 21.
d+00,
6727 & 2.0
d+00, 200.
d+00, 2.0
d+00, 360.
d+00, 2.0
d+00, 175.
d+00,
6729 & 3.0
d+00, 175.
d+00, 3.0
d+00, 2100.
d+00, 2.0
d+00, 384.
d+00,
6731 & 5.0
d+00, 1.73
d+06, 2.0
d+00, 1.59
d+05, 1.1
d+00, 5.55
d+00,
6733 & 0.6
d+00, 24.
d+00, 2.0
d+00, 175.
d+00, 3.5
d+00, 200.
d+00,
6735 & 2.0
d+00, 40.
d+00, 2.0
d+00, 12.8
d+00, 3.0
d+00, 350.
d+00/
6737 DATA plaetb / 0.1
d+00, 0.2
d+00,
6738 & 0.3
d+00, 0.4
d+00, 0.5
d+00, 0.6
d+00, 0.8
d+00, 1.
d+00,
6739 & 1.2
d+00, 1.5
d+00, 2.
d+00, 2.5
d+00, 3.
d+00, 4.
d+00,
6740 & 4.5
d+00, 5.
d+00 /
6743 DATA sieapp / 142.
d+00, 95.1
d+00,
6744 & 75.0
d+00, 70.0
d+00, 62.0
d+00, 57.0
d+00, 48.0
d+00,
6745 & 44.5
d+00, 43.5
d+00, 38.0
d+00, 33.0
d+00 /
6747 DATA sitapp /1129.
d+00, 424.
d+00,
6748 & 239.
d+00, 195.
d+00, 172.
d+00, 150.
d+00, 124.
d+00,
6749 & 117.
d+00, 109.
d+00, 100.
d+00, 90.2
d+00, 81.5
d+00,
6750 & 78.0
d+00, 72.0
d+00, 67.0
d+00, 64.8
d+00 /
6753 ichrge(ktarg)=iich(ktarg)
6754 am(ktarg)=aam(ktarg)
6756 IF ( kp .NE. 26 )
THEN
6772 IF ( iibar(kp) .GT. 0 )
THEN
6778 ELSE IF ( ip .EQ. 15 )
THEN
6784 ELSE IF ( ip .EQ. 24 )
THEN
6790 ELSE IF ( ip .GE. 38 )
THEN
6806 IF ( plab .GT. 50.
d+00 )
THEN
6809 amtsq = am(ktarg)**2
6810 eproj =
sqrt( plab**2 + ampsq )
6811 umosq = ampsq + amtsq + twotwo * am(ktarg) * eproj
6812 eproj =
sqrt( pla**2 + ampsq )
6813 umo50 = ampsq + amtsq + twotwo * am(ktarg) * eproj
6814 umorat =
sqrt( umosq / umo50 )
6818 ELSE IF ( plab .LT. 3.
d+00 )
THEN
6821 amtsq = am(ktarg)**2
6822 eproj =
sqrt( plab**2 + ampsq )
6823 umosq = ampsq + amtsq + twotwo * am(ktarg) * eproj
6824 eproj =
sqrt( pla**2 + ampsq )
6825 umo50 = ampsq + amtsq + twotwo * am(ktarg) * eproj
6826 umorat =
sqrt( umosq / umo50 )
6839 IF ( ihlp(ip) .EQ. 2 )
THEN
6846 sppptt = acof + bcof * pla**enne + ccof * algpla**2
6854 spppel = acof + bcof * pla**enne + ccof * algpla**2
6857 spppin = sppptt - spppel
6864 spmptt = acof + bcof * pla**enne + ccof * algpla**2
6872 spmpel = acof + bcof * pla**enne + ccof * algpla**2
6875 spmpin = spmptt - spmpel
6876 sigdia = spmpin - spppin
6883 IF ( ichrge(ip) .NE. 0 )
THEN
6885 jreac = 3 + ip - 13 + ichrge(ip) * khelp
6886 acof = sgtcoe(1,jreac)
6887 bcof = sgtcoe(2,jreac)
6888 enne = sgtcoe(3,jreac)
6889 ccof = sgtcoe(4,jreac)
6890 dcof = sgtcoe(5,jreac)
6892 shnctt = acof + bcof * pla**enne + ccof * algpla**2
6894 jreac = 19 + ip - 13 + ichrge(ip) * khelp
6895 acof = sgtcoe(1,jreac)
6896 bcof = sgtcoe(2,jreac)
6897 enne = sgtcoe(3,jreac)
6898 ccof = sgtcoe(4,jreac)
6899 dcof = sgtcoe(5,jreac)
6901 shncel = acof + bcof * pla**enne + ccof * algpla**2
6904 shncin = shnctt - shncel
6906 ndiagr = 1 + ip - 13 + ichrge(ip) * khelp
6908 iqfsc1 = 1 + ip - 13
6911 iqbsc2 = 1 + ip - 13
6919 k2hlp = ( kp - 23 ) / 3
6925 ndiagr = 2 + khelp * ( 2 * k2hlp - 1 ) - k2hlp
6926 shncin = hlfhlf * ( spppin + spmpin )
6938 ELSE IF ( ihlp(ip) .EQ. 3 )
THEN
6945 skpptt = acof + bcof * pla**enne + ccof * algpla**2
6953 skppel = acof + bcof * pla**enne + ccof * algpla**2
6956 skppin = skpptt - skppel
6963 skmptt = acof + bcof * pla**enne + ccof * algpla**2
6971 skmpel = acof + bcof * pla**enne + ccof * algpla**2
6974 skmpin = skmptt - skmpel
6975 sigdia = hlfhlf * ( skmpin - skppin )
6978 IF ( ichrge(ip) .NE. 0 )
THEN
6982 IF ( khelp .EQ. 0 )
THEN
6998 shnctt = acof + bcof * pla**enne + ccof * algpla**2
7003 shncin = shnctt - shncel
7022 IF ( khelp .EQ. 0 )
THEN
7023 shncin = skmpin - sigdia
7040 shnctt = acof + bcof * pla**enne + ccof * algpla**2
7045 shncin = shnctt - shncel + sigdia
7062 ELSE IF ( ihlp(ip) .EQ. 4 .AND. ip .LE. 9 )
THEN
7071 sapptt = acof + bcof * pla**enne + ccof * algpla**2
7073 IF ( pla .LT. fivfiv )
THEN
7078 acof = sgtcoe(1,jreac)
7079 bcof = sgtcoe(2,jreac)
7080 enne = sgtcoe(3,jreac)
7081 ccof = sgtcoe(4,jreac)
7082 dcof = sgtcoe(5,jreac)
7084 sappel = acof + bcof * pla**enne + ccof * algpla**2
7087 sappin = sapptt - sappel
7094 spptot = acof + bcof * pla**enne + ccof * algpla**2
7102 sppela = acof + bcof * pla**enne + ccof * algpla**2
7105 sppine = spptot - sppela
7106 sigdia = ( sappin - sppine ) / fivfiv
7110 IF ( ichrge(ip) .NE. 0 )
THEN
7114 IF ( khelp .EQ. 0 )
THEN
7129 shnctt = acof + bcof * pla**enne + ccof * algpla**2
7134 shncin = shnctt - shncel
7142 rnchen =
rndm(rnchen)
7143 IF ( rnchen .LT. puubar )
THEN
7148 iqbsc1 = -iqfsc1 + khelp
7158 IF ( khelp .EQ. 0 )
THEN
7159 shncin = sappin - sigdia
7171 shncin = shnctt - shncel
7179 rnchen =
rndm(rnchen)
7180 IF ( rnchen .LT. pddbar )
THEN
7185 iqbsc1 = -iqfsc1 + khelp - 1
7202 phnsch = ndiagr * sigdia / shncin
7203 iqechc = iqechr(iqfsc1) + iqechr(iqfsc2) + iqechr(iqbsc1)
7205 iqbchc = iqbchr(iqfsc1) + iqbchr(iqfsc2) + iqbchr(iqbsc1)
7209 iqschc = iqschr(iqfsc1) + iqschr(iqfsc2) + iqschr(iqbsc1)
7211 iqspro = iqschr(mquark(1,ip)) + iqschr(mquark(2,ip))
7212 & + iqschr(mquark(3,ip))
7216 WRITE (lunout,*)
' *** Phnsch,kp,ktarg,pla',
7217 &
phnsch,kp,ktarg,pla,
' ****'
7218 WRITE (lunerr,*)
' *** Phnsch,kp,ktarg,pla',
7219 &
phnsch,kp,ktarg,pla,
' ****'
7227 IF ( iqspro .NE. iqschc .OR. ichrge(ip) + ichrge(ktarg)
7228 & .NE. iqechc .OR. iibar(kp) + iibar(ktarg) .NE. iqbchc)
THEN
7230 &
' *** Phnsch,iqspro,iqschc,ichrge,iqechc,ibar,iqbchc,ktarg',
7231 & iqspro,iqschc,ichrge(ip),iqechc,iibar(kp),iqbchc,ktarg
7233 &
' *** Phnsch,iqspro,iqschc,ichrge,iqechc,ibar,iqbchc,ktarg',
7234 & iqspro,iqschc,ichrge(ip),iqechc,iibar(kp),iqbchc,ktarg
7239 IF ( umorat .GT. onepls )
phnsch = oneone / ( ( oneone /
phnsch
7240 & - oneone ) * umorat + oneone )
7243 entry schqua( jqfsc1, jqfsc2, jqbsc1, jqbsc2 )
7298 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7300 COMMON / qquark / iqechr(-6:6), iqbchr(-6:6), iqichr(-6:6),
7301 & iqschr(-6:6), iqcchr(-6:6), iquchr(-6:6),
7302 & iqtchr(-6:6), mquark(3,39)
7306 DATA iqechr / -2, 1, -2, 1, 1, -2, 0, 2, -1, -1, 2, -1, 2 /
7307 DATA iqbchr / 6*-1, 0, 6*1 /
7308 DATA iqichr / 4*0, 1, -1, 0, 1, -1, 4*0 /
7309 DATA iqschr / 3*0, 1, 5*0, -1, 3*0 /
7310 DATA iqcchr / 2*0, -1, 7*0, 1, 2*0 /
7311 DATA iquchr / 0, 1, 9*0, -1, 0 /
7312 DATA iqtchr / -1, 11*0, 1 /
7313 DATA mquark / 1,1,2, -1,-1,-2,
7314 * 0,0,0, 0,0,0, 0,0,0, 0,0,0, 0,0,0,
7315 * 1,2,2, -1,-2,-2, 0,0,0, 0,0,0, 0,0,0,
7316 * 1,-2,0, 2,-1,0, 1,-3,0, 3,-1,0,
7317 * 1,2,3, -1,-2,-3, 0,0,0,
7318 * 2,2,3, 1,1,3, 1,2,3, 1,-1,0,
7319 * 2,-3,0, 3,-2,0, 2,-2,0, 0,0,0,
7320 * 0,0,0, 0,0,0, 0,0,0,
7321 * -1,-1,-3, -1,-2,-3, -2,-2,-3,
7322 * 1,3,3, -1,-3,-3, 2,3,3, -2,-3,-3,
7329 +ptysa2,plaq2,eaq2, amch1,irej,ikvala,pttq1)
7330 IMPLICIT DOUBLE PRECISION (
a-h,o-
z)
7334 COMMON /dprin/ ipri,ipev,ippa,ipco,
init,iphkk,itopd,ipaupr
7336 LOGICAL intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada, ipadis,
7338 COMMON /droppt/ intpt,fermp,ihadss,ihadsv,ihadvs,ihadvv,ihada,
7339 +ipadis,ishmal,lpauli
7341 COMMON /nucc/ it,itz,ip,ipz,ijproj,ibproj,ijtarg,ibtarg
7357 IF (ikvala.EQ.1)b33=6.0
7362 IF (icount.EQ.10)
THEN
7367 IF (icount.GE.1)
THEN
7369 ptxsq1=qtxsq1+hps*cfe
7370 ptysq1=qtysq1+hps*sfe
7371 ptxsa2=qtxsa2-hps*cfe
7372 ptysa2=qtysa2-hps*sfe
7377 es=-2./(b33**2)*
log(abs(
rndm(v)*
rndm(u))+0.00000001)
7378 hps=
sqrt(es*es+2.*es*0.94)
7381 IF (.NOT.intpt) hps=0.0000001
7385 ptxsq1=qtxsq1+hps*cfe
7386 ptysq1=qtysq1+hps*sfe
7387 ptxsa2=qtxsa2-hps*cfe
7388 ptysa2=qtysa2-hps*sfe
7392 IF (ipev.GE.6)
WRITE(6,1000)ptxsq1,ptysq1,
7394 1000
FORMAT (
' PT S ',8f12.6)
7396 pttq1=ptxsq1**2+ptysq1**2
7397 IF((eq1**2.LE.pttq1)) go to 10
7402 IF (ikvala.EQ.1)b33=6.0
7407 IF (icoun2.EQ.12)
THEN
7414 IF (ipev.GE.6)
WRITE(6,1000)ptxsq1,ptysq1,
7417 ptta2=ptxsa2**2+ptysa2**2
7418 IF((eaq2**2.LE.ptta2)) go to 12
7421 IF(ip.GE.1)go to 1779
7422 plq1=
sqrt(eq1**2-pttq1)
7423 plaq2=-
sqrt(eaq2**2-ptta2)
7428 amch1q=(eq1+eaq2)**2-(ptxsq1+ptxsa2)** 2-(ptysq1+ptysa2)**2-(plq1
7430 IF (amch1q.LE.0.d0)
THEN
7432 301
FORMAT(
' inconsistent Kinematics in SELPT AMCH1Q=',e12.3)
7433 WRITE(6,305) qtxsq1,qtysq1,
7434 +qlq1,qeq1,qtxsa1,qtysa1,qlaq1,qeaq1, qtxsq2,qtysq2,qlq2,qeq2,
7436 +qtysa2,qlaq2,qeaq2, amch1,amch2
7437 305
FORMAT(
'PTXSQ1,PTYSQ1,
7438 +PLQ1,EQ1,PTXSA1,PTYSA1,PLAQ1,EAQ1, PTXSQ2,PTYSQ2,PLQ2,EQ2,PTXSA2,
7439 +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)
DOUBLE PRECISION function rndm(RDUMMY)
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
subroutine dtwopd(UMO, ECM1, ECM2, PCM1, PCM2, COD1, COF1, SIF1, COD2, COF2, SIF2, AM1, AM2)
for(Int_t i=0;i< nentries;i++)
DOUBLE PRECISION function dshnto(KPROJ, KTARG, UMO)
DOUBLE PRECISION function sampex(X1, X2)
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)
void fill(G4double x, G4double weight=1.)
subroutine shmak(ICASE, NN, NNA, NNB, NA, NB, UMO, BIMP)
subroutine shmaki(NA, NCA, NB, NCB, RPROJ, RTARG, PPN)
DOUBLE PRECISION function sigsds(S)
subroutine conucl(X, N, R)
subroutine plot(X, Y, N, M, MM, XO, DX, YO, DY)
subroutine py4ent(IP, KF1, KF2, KF3, KF4, PECM, X1, X2, X4, X12, X14)
DOUBLE PRECISION function phnsch(KP, KTARG, PLAB)
DOUBLE PRECISION function dxlamb(X, Y, Z)
subroutine py2ent(IP, KF1, KF2, PECM)
subroutine profb(BSTEP, NSTAT, NA, RA, NB, RB, BSITE, NSITEB)
DOUBLE PRECISION function dshpto(IT, PO)
subroutine previo(RA, RB, NSTB, BMAX, BSTEP, SIG, RO, G)
DOUBLE PRECISION function dshnel(KPROJ, KTARG, UMO)
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
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)
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
DOUBLE PRECISION function siinel(KPROJ, KTARG, UMO)
subroutine drtran(XO, YO, ZO, CDE, SDE, SFE, CFE, X, Y, Z)
DOUBLE PRECISION function sippsd(ECM)
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)
DOUBLE PRECISION function dbeta(X1, X2, BET)
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