21 CHARACTER beam*16,targ*16
22 dimension xsec0(8,0:200),coef0(8,200,20),ini(8),
25 SAVE xsec0, coef0, ini, mint44, mint45
27 common/hijcrdn/
yp(3,300),
yt(3,300)
31 common/histrng/
nfp(300,15),
pp(300,15),
nft(300,15),
pt(300,15)
33 common/hipyint/mint4,mint5,atco(200,20),atxs(0:200)
36 common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
38 common/ludat3/mdcy(500,3),mdme(2000,2),brat(2000),kfdp(2000,5)
40 common/pyhisubs/msel,msub(200),kfin(2,-40:40),ckin(200)
46 common/pyhiint2/iset(200),kfpr(200,2),coef(200,20),icol(40,4,2)
48 common/pyhiint5/
ngen(0:200,3),
xsec(0:200,3)
50 DATA ini/8*0/i_last/-1/
59 ELSE IF(
ihnt2(5).NE.0 .AND.
ihnt2(6).EQ.0)
THEN
61 IF(
nft(jt,4).EQ.2112) i_type=2
62 ELSE IF(
ihnt2(5).EQ.0 .AND.
ihnt2(6).NE.0)
THEN
64 IF(
nfp(jp,4).EQ.2112) i_type=2
66 IF(
nfp(jp,4).EQ.2212 .AND.
nft(jt,4).EQ.2212)
THEN
68 ELSE IF(
nfp(jp,4).EQ.2212 .AND.
nft(jt,4).EQ.2112)
THEN
70 ELSE IF(
nfp(jp,4).EQ.2112 .AND.
nft(jt,4).EQ.2212)
THEN
77 IF(i_trig.NE.0)
go to 160
78 IF(i_trig.EQ.i_last)
go to 150
120 DO 110
j=1,min(8,mdcy(21,3))
121 110 mdme(mdcy(21,2)+
j-1,1)=0
123 IF(
hint1(1).GE.20.0 .and.
ihpr2(18).EQ.1) isel=5
124 mdme(mdcy(21,2)+isel-1,1)=1
130 150
IF(ini(i_type).NE.0)
go to 800
136 IF(i_trig.EQ.i_last)
go to 260
138 ckin(5)=abs(
hipr1(10))-0.25
139 ckin(3)=abs(
hipr1(10))-0.25
140 ckin(4)=abs(
hipr1(10))+0.25
147 IF(
ihpr2(3).EQ.1)
THEN
159 DO 102
j=1,min(8,mdcy(21,3))
160 102 mdme(mdcy(21,2)+
j-1,1)=0
162 IF(
hint1(1).GE.20.0 .and.
ihpr2(18).EQ.1) isel=5
163 mdme(mdcy(21,2)+isel-1,1)=1
165 ELSE IF(
ihpr2(3).EQ.2)
THEN
171 ELSE IF(
ihpr2(3).EQ.3)
THEN
177 DO 105
j=1,min(8,mdcy(21,3))
178 105 mdme(mdcy(21,2)+
j-1,1)=0
180 IF(
hint1(1).GE.20.0 .and.
ihpr2(18).EQ.1) isel=5
181 mdme(mdcy(21,2)+isel-1,1)=1
184 260
IF(ini(i_type).NE.0)
go to 800
189 IF(
nfp(jp,4).EQ.2212)
THEN
191 ELSE IF(
nfp(jp,4).EQ.-2212)
THEN
193 ELSE IF(
nfp(jp,4).EQ.2112)
THEN
195 ELSE IF(
nfp(jp,4).EQ.-2112)
THEN
197 ELSE IF(
nfp(jp,4).EQ.211)
THEN
199 ELSE IF(
nfp(jp,4).EQ.-211)
THEN
201 ELSE IF(
nfp(jp,4).EQ.321)
THEN
203 ELSE IF(
nfp(jp,4).EQ.-321)
THEN
206 WRITE(6,*)
'unavailable beam type',
nfp(jp,4)
208 IF(
nft(jt,4).EQ.2212)
THEN
210 ELSE IF(
nft(jt,4).EQ.-2212)
THEN
212 ELSE IF(
nft(jt,4).EQ.2112)
THEN
214 ELSE IF(
nft(jt,4).EQ.-2112)
THEN
216 ELSE IF(
nft(jt,4).EQ.211)
THEN
218 ELSE IF(
nft(jt,4).EQ.-211)
THEN
220 ELSE IF(
nft(jt,4).EQ.321)
THEN
222 ELSE IF(
nft(jt,4).EQ.-321)
THEN
225 WRITE(6,*)
'unavailable target type',
nft(jt,4)
235 mint44(i_type)=
mint(44)
236 mint45(i_type)=
mint(45)
238 xsec0(i_type,0)=
xsec(0,1)
244 coef0(i_type,
i,
j)=coef(
i,
j)
254 800
mint(44)=mint44(i_type)
255 mint(45)=mint45(i_type)
258 xsec(0,1)=xsec0(i_type,0)
264 coef(
i,
j)=coef0(i_type,
i,
j)