8 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
10 common/pyhipars/mstp(200),parp(200),msti(200),pari(200)
12 common/pyhiint1/mint(400),vint(400)
14 common/pyhiint5/ngen(0:200,3),xsec(0:200,3)
16 dimension bcs(5,8),bcb(2,5),bcc(3)
21 DATA ((bcs(i,j),j=1,8),i=1,5)/
22 1 41.74, 0.66, 0.0000, 337., 0.0, 0.0, -39.3, 0.48,
23 2 41.66, 0.60, 0.0000, 306., 0.0, 0.0, -34.6, 0.51,
24 3 41.36, 0.63, 0.0000, 299., 7.3, 0.5, -40.4, 0.47,
25 4 41.68, 0.63, 0.0083, 330., 0.0, 0.0, -39.0, 0.48,
26 5 41.13, 0.59, 0.0074, 278., 10.5, 0.5, -41.2, 0.46/
27 DATA ((bcb(i,j),j=1,5),i=1,2)/
28 1 10.79, -0.049, 0.040, 21.5, 1.23,
29 2 9.92, -0.027, 0.013, 18.9, 1.07/
30 DATA bcc/2.0164346,-0.5590311,0.0376279/
33 nfit=
min(5,
max(1,mstp(31)))
34 sigp=bcs(nfit,1)+bcs(nfit,2)*(-0.25*paru(1)**2*
35 &(1.-0.25*bcs(nfit,3)*paru(1)**2)+(1.+0.5*bcs(nfit,3)*paru(1)**2)*
36 &(
log(vint(2)/bcs(nfit,4)))**2+bcs(nfit,3)*
37 &(
log(vint(2)/bcs(nfit,4)))**4)/
38 &((1.-0.25*bcs(nfit,3)*paru(1)**2)**2+2.*bcs(nfit,3)*
39 &(1.+0.25*bcs(nfit,3)*paru(1)**2)*(
log(vint(2)/bcs(nfit,4)))**2+
40 &bcs(nfit,3)**2*(
log(vint(2)/bcs(nfit,4)))**4)+bcs(nfit,5)*
41 &vint(2)**(bcs(nfit,6)-1.)*
sin(0.5*paru(1)*bcs(nfit,6))
42 sigm=-bcs(nfit,7)*vint(2)**(bcs(nfit,8)-1.)*
43 &
cos(0.5*paru(1)*bcs(nfit,8))
44 refp=bcs(nfit,2)*paru(1)*
log(vint(2)/bcs(nfit,4))/
45 &((1.-0.25*bcs(nfit,3)*paru(1)**2)**2+2.*bcs(nfit,3)*
46 &(1.+0.25*bcs(nfit,3)*paru(1)**2)+(
log(vint(2)/bcs(nfit,4)))**2+
47 &bcs(nfit,3)**2*(
log(vint(2)/bcs(nfit,4)))**4)-bcs(nfit,5)*
48 &vint(2)**(bcs(nfit,6)-1.)*
cos(0.5*paru(1)*bcs(nfit,6))
49 refm=-bcs(nfit,7)*vint(2)**(bcs(nfit,8)-1.)*
50 &
sin(0.5*paru(1)*bcs(nfit,8))
51 sigma=sigp-isign(1,mint(11)*mint(12))*sigm
52 rho=(refp-isign(1,mint(11)*mint(12))*refm)/sigma
56 IF(mstp(31).GE.4) nfit=2
57 bp=bcb(nfit,1)+bcb(nfit,2)*
log(vint(2))+
58 &bcb(nfit,3)*(
log(vint(2)))**2
59 bm=bcb(nfit,4)+bcb(nfit,5)*
log(vint(2))
60 b=
bp-isign(1,mint(11)*mint(12))*sigm/sigp*(bm-
bp)
62 c=-0.5*bcc(2)/bcc(3)*(1.-
sqrt(
max(0.,1.+4.*bcc(3)/bcc(2)**2*
63 &(1.
e-03*vint(1)-bcc(1)))))
67 sigel=sigma**2*(1.+
rho**2)/(16.*paru(1)*paru(5)*
b)
70 sigsd=2.*0.68*(1.+36./vint(2))*
log(0.6+0.1*vint(2))
74 sigdd=sigsd**2/(3.*sigel)
77 signd=sigma-sigdd-sigsd-sigel
80 IF(iabs(mint(11)).EQ.211.AND.iabs(mint(12)).EQ.211)
THEN
86 ELSEIF(iabs(mint(11)).EQ.211.OR.iabs(mint(12)).EQ.211)
THEN