21 SUBROUTINE inietc(ITAUXPAR,xpar)
29 COMMON / taurad / xk0dec,itdkrc
30 DOUBLE PRECISION XK0DEC
31 COMMON / jaki / jak1,jak2,jakp,jakm,ktom
44 idff = xpar(itauxpar+3)
46 xk0dec = xpar(itauxpar+5)
48 itdkrc = xpar(itauxpar+4)
50 jak1 = xpar(itauxpar+1)
51 jak2 = xpar(itauxpar+2)
55 keya1 = xpar(itauxpar+6)
58 WRITE(iout,bxtxt)
' Parameters passed from KK to Tauola: ' 59 WRITE(iout,bxl1i) jak1,
'dec. type 1-st tau ',
'Jak1 ',
't01' 60 WRITE(iout,bxl1i) jak2,
'dec. type 2-nd tau ',
'Jak2 ',
't02' 61 WRITE(iout,bxl1i) keya1,
'current type a1 dec.',
'KeyA1 ',
't03' 62 WRITE(iout,bxl1i) idff,
'PDG id 1-st tau ',
'idff ',
't04' 63 WRITE(iout,bxl1i) itdkrc,
'R.c. switch lept dec',
'itdkRC',
't05' 64 WRITE(iout,bxl1g) xk0dec,
'IR-cut for lept r.c.',
'xk0dec',
't06' 69 SUBROUTINE initdk(ITAUXPAR,xpar)
82 COMMON / decpar / gfermi,gv,ga,ccabib,scabib,gamel
83 REAL*4 GFERMI,GV,GA,CCABIB,SCABIB,GAMEL
84 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
85 * ,ampiz,ampi,amro,gamro,ama1,gama1
86 * ,amk,amkz,amkst,gamkst
88 REAL*4 AMTAU,AMNUTA,AMEL,AMNUE,AMMU,AMNUMU
89 * ,AMPIZ,AMPI,AMRO,GAMRO,AMA1,GAMA1
90 * ,AMK,AMKZ,AMKST,GAMKST
91 COMMON / taubra / gamprt(500),jlist(500),nchan
92 COMMON / taukle / bra1,brk0,brk0b,brks
93 REAL*4 BRA1,BRK0,BRK0B,BRKS
95 parameter(nmode=86,nm1=0,nm2=11,nm3=19,nm4=22,nm5=21,nm6=13)
96 COMMON / taudcd /idffin(9,nmode),mulpik(nmode)
98 CHARACTER NAMES(NMODE)*31
100 CHARACTER OLDNAMES(7)*31
103 $ bxinit =
'(1x,1h*,g17.8, 16x, a31,a4,a4, 1x,1h*)' 127 dimension nopik(9,nmode),npik(nmode)
175 DATA nopik / -1,-1, 1, 2, 0, 0,3*0, 2, 2, 2,-1, 0, 0,3*0,
176 a 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
177 b 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
178 c 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
179 d 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
180 e 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
181 a 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
182 b 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
183 c 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
184 d 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
185 e 4, 2, 2,-1, 0, 0,3*0, 4, 2, 2,-1, 0, 0,3*0,
186 1 -1,-1, 1, 2, 2, 0,3*0,
187 a -1,-1, 1, 2, 4, 0,3*0, -1,-1, 1, 2, 4, 0,3*0,
188 a 1,-1,-1, 2, 2, 0,3*0, -1, 2, 2, 2, 2, 0,3*0,
189 a -1, 1, 1,-1,-1, 0,3*0, -1,-1, 1, 2, 4, 0,3*0,
190 a -1,-1, 1, 2, 4, 0,3*0, -1,-1, 1, 2, 4, 0,3*0,
191 a -1,-1, 1, 2, 4, 0,3*0, -1,-1, 1, 2, 4, 0,3*0,
192 a -1,-1, 1, 2, 4, 0,3*0, -1,-1, 1, 2, 4, 0,3*0,
193 a -1,-1, 1, 2, 4, 0,3*0, -1,-1, 1, 2, 4, 0,3*0,
194 a -1,-1, 1, 2, 4, 0,3*0, -1,-1, 1, 2, 4, 0,3*0,
195 a -1,-1, 1, 2, 4, 0,3*0, -1,-1, 1, 2, 4, 0,3*0,
196 a -1,-1, 1, 2, 4, 0,3*0, -1,-1, 1, 2, 4, 0,3*0,
197 x -1,-1,-1, 1, 1, 0,3*0,
198 2 -1,-1,-1, 1, 1, 2,3*0, -1,-1, 1, 2, 2, 2,3*0,
199 a -1,-1,-1, 1, 1, 1,3*0, -1,-1, 1, 2, 2, 1,3*0,
200 b -1,-1,-1, 1, 1, 1,3*0, -1,-1, 1, 2, 2, 1,3*0,
201 c -1,-1,-1, 1, 1, 1,3*0, -1,-1, 1, 2, 2, 1,3*0,
202 d -1,-1,-1, 1, 1, 1,3*0, -1,-1, 1, 2, 2, 1,3*0,
203 e -1,-1,-1, 1, 1, 1,3*0, -1,-1, 1, 2, 2, 1,3*0,
204 3 -3,-1, 3, 0, 0, 0,3*0, -4,-1, 4, 0, 0, 0,3*0,
205 4 -3, 2,-4, 0, 0, 0,3*0, 2, 2,-3, 0, 0, 0,3*0,
206 5 -3,-1, 1, 0, 0, 0,3*0, -1, 4, 2, 0, 0, 0,3*0,
207 6 9,-1, 2, 0, 0, 0,3*0, -1, 2, 8, 0, 0, 0,3*0,
211 7 2, 2,-1, 0, 0, 0,3*0,
212 7 2, 2, 2, 0, 0, 0,3*0, 2, 2, 2, 0, 0, 0,3*0,
213 7 2, 2, 2, 0, 0, 0,3*0, 2, 2, 2, 0, 0, 0,3*0,
214 7 2, 2, 2, 0, 0, 0,3*0, 2, 2, 2, 0, 0, 0,3*0,
215 7 2, 2, 2, 0, 0, 0,3*0, 2, 2, 2, 0, 0, 0,3*0,
216 7 2, 2, 2, 0, 0, 0,3*0, 2, 2, 2, 0, 0, 0,3*0,
218 8 -3,-4, 0, 0, 0, 0,3*0,
219 8 -3,-3, 0, 0, 0, 0,3*0, -3,-3, 0, 0, 0, 0,3*0,
220 8 -3,-3, 0, 0, 0, 0,3*0, -3,-3, 0, 0, 0, 0,3*0,
221 8 -3,-3, 0, 0, 0, 0,3*0, -3,-3, 0, 0, 0, 0,3*0,
222 8 -3,-3, 0, 0, 0, 0,3*0, -3,-3, 0, 0, 0, 0,3*0,
223 8 -3,-3, 0, 0, 0, 0,3*0, -3,-3, 0, 0, 0, 0,3*0 /
232 IF(i.EQ. 1) gamprt(i) = 1.0000
233 IF(i.EQ. 2) gamprt(i) = 1.0000
234 IF(i.EQ. 3) gamprt(i) = 1.0000
235 IF(i.EQ. 4) gamprt(i) = 1.0000
236 IF(i.EQ. 5) gamprt(i) = 1.0000
237 IF(i.EQ. 6) gamprt(i) = 1.0000
238 IF(i.EQ. 7) gamprt(i) = 1.0000
239 IF(i.EQ. 8) gamprt(i) = 1.0000
240 IF(i.EQ. 9) gamprt(i) = 1.0000
241 IF(i.EQ.10) gamprt(i) = 1.0000
242 IF(i.EQ.11) gamprt(i) = 1.0000
243 IF(i.EQ.12) gamprt(i) = 1.0000
244 IF(i.EQ.13) gamprt(i) = 1.0000
245 IF(i.EQ.14) gamprt(i) = 1.0000
246 IF(i.EQ.15) gamprt(i) = 1.0000
247 IF(i.EQ.16) gamprt(i) = 1.0000
248 IF(i.EQ.17) gamprt(i) = 1.0000
249 IF(i.EQ.18) gamprt(i) = 1.0000
250 IF(i.EQ.19) gamprt(i) = 1.0000
251 IF(i.EQ.20) gamprt(i) = 1.0000
252 IF(i.EQ.21) gamprt(i) = 1.0000
253 IF(i.EQ.22) gamprt(i) = 1.0000
254 IF(i.GT.22.AND.i.LE.93) gamprt(i) = 1.0000
256 IF(i.GT.0.AND.i.LE.93) gamprt(i) = 0.0000
257 IF(i.EQ. 1) gamprt(i) =0.1800
258 IF(i.EQ. 2) gamprt(i) =0.1751
259 IF(i.EQ. 3) gamprt(i) =0.1110
260 IF(i.EQ. 4) gamprt(i) =0.2515
261 IF(i.EQ. 5) gamprt(i) =0.1790 /2
262 IF(i.EQ. 6) gamprt(i) =0.0071
263 IF(i.EQ. 7) gamprt(i) =0.0134
264 IF(i.EQ. 8) gamprt(i) =0.0450
265 IF(i.EQ. 9) gamprt(i) =0.0100
267 IF(i.EQ.30) gamprt(i) =0.0009
268 IF(i.EQ.33) gamprt(i) =0.004
269 IF(i.EQ.34) gamprt(i) =0.002
270 IF(i.EQ.35) gamprt(i) =0.001
272 IF(i.EQ.51) gamprt(i) =0.0004
273 IF(i.EQ.52) gamprt(i) =0.0003
274 IF(i.EQ.53) gamprt(i) =0.0005
276 IF(i.EQ.64) gamprt(i) =0.0015
277 IF(i.EQ.65) gamprt(i) =0.0015
278 IF(i.EQ.66) gamprt(i) =0.0015
279 IF(i.EQ.67) gamprt(i) =0.0005
280 IF(i.EQ.68) gamprt(i) =0.0050
281 IF(i.EQ.69) gamprt(i) =0.0055
282 IF(i.EQ.70) gamprt(i) =0.0017
283 IF(i.EQ.71) gamprt(i) =0.0013
284 IF(i.EQ.72) gamprt(i) =0.1790 /2
286 IF(i.EQ.83) gamprt(i) =0.0010
288 IF(i.EQ. 1) oldnames(i)=
' TAU- --> E- ' 289 IF(i.EQ. 2) oldnames(i)=
' TAU- --> MU- ' 290 IF(i.EQ. 3) oldnames(i)=
' TAU- --> PI- ' 291 IF(i.EQ. 4) oldnames(i)=
' TAU- --> PI-, PI0 ' 292 IF(i.EQ. 5) oldnames(i)=
' TAU- --> PI-, PI-, PI+ ' 293 IF(i.EQ. 6) oldnames(i)=
' TAU- --> K- ' 294 IF(i.EQ. 7) oldnames(i)=
' TAU- --> K*- (two subch) ' 295 IF(i.EQ. 8) names(i-7)=
' TAU- --> 2PI-, PI0, PI+ ' 296 IF(i.EQ. 9) names(i-7)=
' TAU- --> 3PI0, PI- ' 298 IF(i.EQ.10) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 299 IF(i.EQ.11) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 300 IF(i.EQ.12) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 301 IF(i.EQ.13) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 302 IF(i.EQ.14) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 303 IF(i.EQ.15) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 304 IF(i.EQ.16) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 305 IF(i.EQ.17) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 306 IF(i.EQ.18) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 307 IF(i.EQ.19) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 308 IF(i.EQ.20) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 309 IF(i.EQ.21) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 310 IF(i.EQ.22) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 311 IF(i.EQ.23) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 312 IF(i.EQ.24) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 313 IF(i.EQ.25) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 314 IF(i.EQ.26) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 315 IF(i.EQ.27) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 316 IF(i.EQ.28) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 317 IF(i.EQ.29) names(i-7)=
' TAU- --> xxxxxxx4xxxxxxxx ' 320 IF(i.EQ.30) names(i-7)=
' TAU- --> 2PI-, PI+, 2PI0 ' 322 IF(i.EQ.31) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 323 IF(i.EQ.32) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 324 IF(i.EQ.33) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 325 IF(i.EQ.34) names(i-7)=
' TAU- --> PI- 4PI0 ' 326 IF(i.EQ.35) names(i-7)=
' TAU- --> 3PI- 2PI+ ' 327 IF(i.EQ.36) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 328 IF(i.EQ.37) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 329 IF(i.EQ.38) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 330 IF(i.EQ.39) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 331 IF(i.EQ.40) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 333 IF(i.EQ.41) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 334 IF(i.EQ.42) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 335 IF(i.EQ.43) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 336 IF(i.EQ.44) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 337 IF(i.EQ.45) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 338 IF(i.EQ.46) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 339 IF(i.EQ.47) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 340 IF(i.EQ.48) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 341 IF(i.EQ.49) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 342 IF(i.EQ.50) names(i-7)=
' TAU- --> xxxxxxxxx5xxxxxx ' 344 IF(i.EQ.51) names(i-7)=
' TAU- --> 3PI-, 2PI+, ' 345 IF(i.EQ.52) names(i-7)=
' TAU- --> 3PI-, 2PI+, PI0 ' 346 IF(i.EQ.53) names(i-7)=
' TAU- --> 2PI-, PI+, 3PI0 ' 347 IF(i.EQ.54) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 348 IF(i.EQ.55) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 349 IF(i.EQ.56) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 350 IF(i.EQ.57) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 351 IF(i.EQ.58) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 352 IF(i.EQ.59) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 353 IF(i.EQ.60) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 354 IF(i.EQ.61) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 355 IF(i.EQ.62) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 356 IF(i.EQ.63) names(i-7)=
' TAU- --> xxxxxxxxxnxxxxxx ' 358 IF(i.EQ.64) names(i-7)=
' TAU- --> K-, PI-, K+ ' 359 IF(i.EQ.65) names(i-7)=
' TAU- --> K0, PI-, K0B ' 361 IF(i.EQ.66) names(i-7)=
' TAU- --> K-, K0, PI0 ' 363 IF(i.EQ.67) names(i-7)=
' TAU- --> PI0 PI0 K- ' 364 IF(i.EQ.68) names(i-7)=
' TAU- --> K- PI- PI+ ' 365 IF(i.EQ.69) names(i-7)=
' TAU- --> PI- K0B PI0 ' 366 IF(i.EQ.70) names(i-7)=
' TAU- --> ETA PI- PI0 ' 367 IF(i.EQ.71) names(i-7)=
' TAU- --> PI- PI0 GAM ' 368 IF(i.EQ.72) names(i-7)=
' TAU- --> PI- PI0 PI0 ' 369 IF(i.EQ.73) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 370 IF(i.EQ.74) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 371 IF(i.EQ.75) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 372 IF(i.EQ.76) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 373 IF(i.EQ.77) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 374 IF(i.EQ.78) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 375 IF(i.EQ.79) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 376 IF(i.EQ.80) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 377 IF(i.EQ.81) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 378 IF(i.EQ.82) names(i-7)=
' TAU- --> xxxxxxxxx3xxxxxx ' 381 IF(i.EQ.83) names(i-7)=
' TAU- --> K- K0 ' 382 IF(i.EQ.84) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 383 IF(i.EQ.85) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 384 IF(i.EQ.86) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 385 IF(i.EQ.87) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 386 IF(i.EQ.88) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 387 IF(i.EQ.89) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 388 IF(i.EQ.90) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 389 IF(i.EQ.91) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 390 IF(i.EQ.92) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 391 IF(i.EQ.93) names(i-7)=
' TAU- --> xxxxxxxxx2xxxxxx ' 401 idffin(j,i)=nopik(j,i)
428 IF (xpar(itauxpar+100+1).GT.-1d0)
THEN 430 ccabib = xpar(itauxpar+7)
431 gv = xpar(itauxpar+8)
432 ga = xpar(itauxpar+9)
435 brks = xpar(itauxpar+11)
436 brk0 = xpar(itauxpar+12)
437 brk0b = xpar(itauxpar+13)
441 IF(k.EQ. 1) gamprt(k) =xpar(itauxpar+100+k)
442 IF(k.EQ. 2) gamprt(k) =xpar(itauxpar+100+k)
443 IF(k.EQ. 3) gamprt(k) =xpar(itauxpar+100+k)
444 IF(k.EQ. 4) gamprt(k) =xpar(itauxpar+100+k)
445 IF(k.EQ. 5) gamprt(k) =xpar(itauxpar+100+k) /2
446 IF(k.EQ. 6) gamprt(k) =xpar(itauxpar+100+k)
447 IF(k.EQ. 7) gamprt(k) =xpar(itauxpar+100+k)
448 IF(k.EQ. 8) gamprt(k) =xpar(itauxpar+100+k)
449 IF(k.EQ. 9) gamprt(k) =xpar(itauxpar+100+k)
451 IF(k.EQ.30) gamprt(k) =xpar(itauxpar+100+10)
453 IF(k.EQ.51) gamprt(k) =xpar(itauxpar+100+11)
454 IF(k.EQ.52) gamprt(k) =xpar(itauxpar+100+12)
455 IF(k.EQ.53) gamprt(k) =xpar(itauxpar+100+13)
457 IF(k.EQ.64) gamprt(k) =xpar(itauxpar+100+14)
458 IF(k.EQ.65) gamprt(k) =xpar(itauxpar+100+15)
459 IF(k.EQ.66) gamprt(k) =xpar(itauxpar+100+16)
460 IF(k.EQ.67) gamprt(k) =xpar(itauxpar+100+17)
461 IF(k.EQ.68) gamprt(k) =xpar(itauxpar+100+18)
462 IF(k.EQ.69) gamprt(k) =xpar(itauxpar+100+19)
463 IF(k.EQ.70) gamprt(k) =xpar(itauxpar+100+20)
464 IF(k.EQ.71) gamprt(k) =xpar(itauxpar+100+21)
465 IF(k.EQ.72) gamprt(k) =xpar(itauxpar+100+5) /2
467 IF(k.EQ.83) gamprt(k) =xpar(itauxpar+100+22)
469 IF(k.EQ.33) gamprt(i) =xpar(itauxpar+100+23)
470 IF(k.EQ.34) gamprt(i) =xpar(itauxpar+100+24)
471 IF(k.EQ.35) gamprt(i) =xpar(itauxpar+100+25)
476 scabib = sqrt(1.-ccabib**2)
478 gamel = gfermi**2*amtau**5/(192*pi**3)
491 WRITE(iout,bxtxt)
' TAUOLA Initialization SUBROUTINE INITDK: ' 492 WRITE(iout,bxtxt)
' Adopted to read from KK ' 493 WRITE(iout,bxtxt)
' ' 494 WRITE(iout,bxtxt)
' Choice Probability -- Decay Channel' 496 WRITE(iout,bxinit) gamprt(k)/sum, oldnames(k),
'****',
'***' 499 WRITE(iout,bxinit) gamprt(k)/sum, names(k-7),
'****',
'***' 501 WRITE(iout,bxtxt)
' In addition:' 502 WRITE(iout,bxinit) gv,
'Vector W-tau-nu coupl. ',
'****',
'***' 503 WRITE(iout,bxinit) ga,
'Axial W-tau-nu coupl. ',
'****',
'***' 504 WRITE(iout,bxinit) gfermi,
'Fermi Coupling ',
'****',
'***' 505 WRITE(iout,bxinit) ccabib,
'cabibo angle ',
'****',
'***' 506 WRITE(iout,bxinit) bra1,
'a1 br ratio (massless) ',
'****',
'***' 507 WRITE(iout,bxinit) brks,
'K* br ratio (massless) ',
'****',
'***' 513 SUBROUTINE iniphy(XK00)
518 COMMON / qedprm /alfinv,alfpi,xk0
519 REAL*8 ALFINV,ALFPI,XK0
522 pi8 = 4.d0*datan(1.d0)
524 alfpi = 1d0/(alfinv*pi8)
528 SUBROUTINE inimas(ITAUXPAR,xpar)
540 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
541 * ,ampiz,ampi,amro,gamro,ama1,gama1
542 * ,amk,amkz,amkst,gamkst
544 REAL*4 AMTAU,AMNUTA,AMEL,AMNUE,AMMU,AMNUMU
545 * ,AMPIZ,AMPI,AMRO,GAMRO,AMA1,GAMA1
546 * ,AMK,AMKZ,AMKST,GAMKST
549 $ bxinit =
'(1x,1h*,g17.8, 16x, a31,a4,a4, 1x,1h*)' 576 WRITE(iout,bxtxt)
' TAUOLA Initialization SUBROUTINE INIMAS: ' 577 WRITE(iout,bxtxt)
' Adopted to read from KK ' 578 WRITE(iout,bxinit) amtau,
'AMTAU tau-mass ',
'****',
'***' 579 WRITE(iout,bxinit) amel ,
'AMEL electron-mass ',
'****',
'***' 580 WRITE(iout,bxinit) ammu ,
'AMMU muon-mass ',
'****',
'***' 584 SUBROUTINE choice(MNUM,RR,ICHAN,PROB1,PROB2,PROB3,
585 $ AMRX,GAMRX,AMRA,GAMRA,AMRB,GAMRB)
586 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
587 * ,ampiz,ampi,amro,gamro,ama1,gama1
588 * ,amk,amkz,amkst,gamkst
590 REAL*4 AMTAU,AMNUTA,AMEL,AMNUE,AMMU,AMNUMU
591 * ,ampiz,ampi,amro,gamro,ama1,gama1
592 * ,amk,amkz,amkst,gamkst
608 ELSEIF(mnum.EQ.1)
THEN 617 ELSEIF(mnum.EQ.2)
THEN 626 ELSEIF(mnum.EQ.3)
THEN 635 ELSEIF(mnum.EQ.4)
THEN 644 ELSEIF(mnum.EQ.5)
THEN 653 ELSEIF(mnum.EQ.6)
THEN 662 ELSEIF(mnum.EQ.7)
THEN 671 ELSEIF(mnum.EQ.8)
THEN 680 ELSEIF(mnum.EQ.9)
THEN 689 ELSEIF(mnum.EQ.101)
THEN 698 ELSEIF(mnum.EQ.102)
THEN 707 ELSEIF(mnum.GE.103.AND.mnum.LE.112)
THEN 729 IF (rr.LE.prob1)
THEN 731 ELSEIF(rr.LE.(prob1+prob2))
THEN 746 prob3=1.0-prob1-prob2
748 FUNCTION dcdmas(IDENT)
749 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
750 * ,ampiz,ampi,amro,gamro,ama1,gama1
751 * ,amk,amkz,amkst,gamkst
753 REAL*4 AMTAU,AMNUTA,AMEL,AMNUE,AMMU,AMNUMU
754 * ,ampiz,ampi,amro,gamro,ama1,gama1
755 * ,amk,amkz,amkst,gamkst
756 IF (ident.EQ. 1)
THEN 758 ELSEIF (ident.EQ.-1)
THEN 760 ELSEIF (ident.EQ. 2)
THEN 762 ELSEIF (ident.EQ.-2)
THEN 764 ELSEIF (ident.EQ. 3)
THEN 766 ELSEIF (ident.EQ.-3)
THEN 768 ELSEIF (ident.EQ. 4)
THEN 770 ELSEIF (ident.EQ.-4)
THEN 772 ELSEIF (ident.EQ. 8)
THEN 774 ELSEIF (ident.EQ.-8)
THEN 776 ELSEIF (ident.EQ. 9)
THEN 778 ELSEIF (ident.EQ.-9)
THEN 781 print *,
'STOP IN APKMAS, WRONG IDENT=',ident
786 FUNCTION lunpik(ID,ISGN)
787 COMMON / taukle / bra1,brk0,brk0b,brks
788 REAL*4 BRA1,BRK0,BRK0B,BRKS
792 IF (ident.EQ. 1)
THEN 794 ELSEIF (ident.EQ.-1)
THEN 796 ELSEIF (ident.EQ. 2)
THEN 798 ELSEIF (ident.EQ.-2)
THEN 800 ELSEIF (ident.EQ. 3)
THEN 802 ELSEIF (ident.EQ.-3)
THEN 805 ELSEIF (ident.EQ. 4)
THEN 809 IF (xio(1).GT.brk0)
THEN 814 ELSEIF (ident.EQ.-4)
THEN 818 IF (xio(1).GT.brk0b)
THEN 823 ELSEIF (ident.EQ. 8)
THEN 825 ELSEIF (ident.EQ.-8)
THEN 827 ELSEIF (ident.EQ. 9)
THEN 829 ELSEIF (ident.EQ.-9)
THEN 832 print *,
'STOP IN IPKDEF, WRONG IDENT=',ident
841 SUBROUTINE taurdf(KTO)
845 COMMON / taukle / bra1,brk0,brk0b,brks
846 REAL*4 BRA1,BRK0,BRK0B,BRKS
847 COMMON / taubra / gamprt(500),jlist(500),nchan
856 IF(i.EQ. 1) gamprt(i) = .0000
857 IF(i.EQ. 2) gamprt(i) = .0000
858 IF(i.EQ. 3) gamprt(i) = .0000
859 IF(i.EQ. 4) gamprt(i) = .0000
860 IF(i.EQ. 5) gamprt(i) = .0000
861 IF(i.EQ. 6) gamprt(i) = .0000
862 IF(i.EQ. 7) gamprt(i) = .0000
863 IF(i.EQ. 8) gamprt(i) = 1.0000
864 IF(i.EQ. 9) gamprt(i) = 1.0000
865 IF(i.EQ.10) gamprt(i) = 1.0000
866 IF(i.EQ.11) gamprt(i) = 1.0000
867 IF(i.EQ.12) gamprt(i) = 1.0000
868 IF(i.EQ.13) gamprt(i) = 1.0000
869 IF(i.EQ.14) gamprt(i) = 1.0000
870 IF(i.EQ.15) gamprt(i) = 1.0000
871 IF(i.EQ.16) gamprt(i) = 1.0000
872 IF(i.EQ.17) gamprt(i) = 1.0000
873 IF(i.EQ.18) gamprt(i) = 1.0000
874 IF(i.EQ.19) gamprt(i) = 1.0000
899 IF(i.EQ. 1) gamprt(i) = .0000
900 IF(i.EQ. 2) gamprt(i) = .0000
901 IF(i.EQ. 3) gamprt(i) = .0000
902 IF(i.EQ. 4) gamprt(i) = .0000
903 IF(i.EQ. 5) gamprt(i) = .0000
904 IF(i.EQ. 6) gamprt(i) = .0000
905 IF(i.EQ. 7) gamprt(i) = .0000
906 IF(i.EQ. 8) gamprt(i) = 1.0000
907 IF(i.EQ. 9) gamprt(i) = 1.0000
908 IF(i.EQ.10) gamprt(i) = 1.0000
909 IF(i.EQ.11) gamprt(i) = 1.0000
910 IF(i.EQ.12) gamprt(i) = 1.0000
911 IF(i.EQ.13) gamprt(i) = 1.0000
912 IF(i.EQ.14) gamprt(i) = 1.0000
913 IF(i.EQ.15) gamprt(i) = 1.0000
914 IF(i.EQ.16) gamprt(i) = 1.0000
915 IF(i.EQ.17) gamprt(i) = 1.0000
916 IF(i.EQ.18) gamprt(i) = 1.0000
917 IF(i.EQ.19) gamprt(i) = 1.0000