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