10 COMMON / / blan(10000)
12 COMMON / inout / inut,iout
18 OPEN(iout,file=
"./tauola.output")
19 OPEN(inut,file=
"./dane.dat")
25 SUBROUTINE dectes(KTORY)
28 DOUBLE PRECISION HH(4)
30 COMMON / jaki / jak1,jak2,jakp,jakm,ktom
33 COMMON / inout / inut,iout
38 COMMON / taurad / xk0dec,itdkrc
61 READ( ninp,3000) testit
62 WRITE(nout,3000) testit
63 READ( ninp,3001) kat1,kat2,kat3,kat4,kat5,kat6
64 READ( ninp,3002) nevt,jak1,jak2,itdkrc
65 READ( ninp,3003) ptau,xk0dec
69 WRITE(nout,
'(6A6/6I6)')
70 $
'KAT1',
'KAT2',
'KAT3',
'KAT4',
'KAT5',
'KAT6',
71 $ kat1 , kat2 , kat3 , kat4 , kat5 , kat6
72 WRITE(nout,
'(4A12/4I12)')
73 $
'NEVT',
'JAK1',
'JAK2',
'ITDKRC',
74 $ nevt, jak1 , jak2 , itdkrc
75 WRITE(nout,
'(2A12/2F12.6)')
92 print *,
'for the sake of these tests KTO has to be 2'
93 print *,
'to change tau- to tau+ change IDFF from -15 to 15'
105 print *,
'NEVTES= ',nevtes
106 WRITE(iout,7011) keya1
109 WRITE(iout,7001) jak,idff,pol(3),ptau
111 WRITE(iout,7004) jak,idff,pol(3),ptau
132 #if defined (history)
142 CALL dekay(kto+10,hh)
154 IF(ipri.EQ.1) print *,
' event no: ',nev,
' NEVTES: ',nevtes
166 7001
FORMAT(//4(/1x,15(5h=====))
167 $ /,
' ', 19x,
' TEST OF RAD. CORR IN ELECTRON DECAY ',9x,1h ,
168 $ /,
' ', 19x,
' TESTS OF TAU DECAY ROUTINES ',9x,1h ,
169 $ /,
' ', 19x,
' INTERFACE OF THE KORAL-Z TYPE ',9x,1h ,
170 $ 2(/,1x,15(5h=====)),
171 $ /,5x ,
'JAK =',i7 ,
' KEY DEFINING DECAY TYPE ',9x,1h ,
172 $ /,5x ,
'IDFF =',i7 ,
' LUND IDENTIFIER FOR FIRST TAU ',9x,1h ,
173 $ /,5x ,
'POL(3)=',f7.2,
' THIRD COMPONENT OF TAU POLARIZ. ',9x,1h ,
174 $ /,5x ,
'PTAU =',f7.2,
' THIRD COMPONENT OF TAU MOM. GEV ',9x,1h ,
175 $ 2(/,1x,15(5h=====))/)
176 7002
FORMAT(///1x,
'===== EVENT NO.',i4,1x,5h=====)
177 7003
FORMAT(5x,
'POLARIMETRIC VECTOR: ',
178 $ 7x,
'HH(1)',7x,
'HH(2)',7x,
'HH(3)',7x,
'HH(4)',
179 $ /, 5x,
' ', 4(1x,f11.8) )
180 7004
FORMAT(//4(/1x,15(5h=====))
181 $ /,
' ', 19x,
' TEST OF RAD. CORR IN ELECTRON DECAY ',9x,1h ,
182 $ /,
' ', 19x,
' TESTS OF TAU DECAY ROUTINES ',9x,1h ,
183 $ /,
' ', 19x,
' INTERFACE OF THE KORAL-B TYPE ',9x,1h ,
184 $ 2(/,1x,15(5h=====)),
185 $ /,5x ,
'JAK =',i7 ,
' KEY DEFINING DECAY TYPE ',9x,1h ,
186 $ /,5x ,
'IDFF =',i7 ,
' LUND IDENTIFIER FOR FIRST TAU ',9x,1h ,
187 $ /,5x ,
'POL(3)=',f7.2,
' THIRD COMPONENT OF TAU POLARIZ. ',9x,1h ,
188 $ /,5x ,
'PTAU =',f7.2,
' THIRD COMPONENT OF TAU MOM. GEV ',9x,1h ,
189 $ 2(/,1x,15(5h=====))/)
190 7011
FORMAT(///1x,
'===== TYPE OF CURRENT',i4,1x,5h=====)
192 SUBROUTINE choice(MNUM,RR,ICHAN,PROB1,PROB2,PROB3,
193 $ AMRX,GAMRX,AMRA,GAMRA,AMRB,GAMRB)
194 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
195 * ,ampiz,ampi,amro,gamro,ama1,gama1
196 * ,amk,amkz,amkst,gamkst
198 real*4 amtau,amnuta,amel,amnue,ammu,amnumu
199 * ,ampiz,ampi,amro,gamro,ama1,gama1
200 * ,amk,amkz,amkst,gamkst
216 ELSEIF(mnum.EQ.1)
THEN
225 ELSEIF(mnum.EQ.2)
THEN
234 ELSEIF(mnum.EQ.3)
THEN
243 ELSEIF(mnum.EQ.4)
THEN
252 ELSEIF(mnum.EQ.5)
THEN
261 ELSEIF(mnum.EQ.6)
THEN
270 ELSEIF(mnum.EQ.7)
THEN
279 ELSEIF(mnum.EQ.8)
THEN
288 ELSEIF(mnum.EQ.101)
THEN
297 ELSEIF(mnum.EQ.102)
THEN
317 IF (rr.LE.prob1)
THEN
319 ELSEIF(rr.LE.(prob1+prob2))
THEN
334 prob3=1.0-prob1-prob2
343 COMMON / decpar / gfermi,gv,ga,ccabib,scabib,gamel
344 real*4 gfermi,gv,ga,ccabib,scabib,gamel
345 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
346 * ,ampiz,ampi,amro,gamro,ama1,gama1
347 * ,amk,amkz,amkst,gamkst
349 real*4 amtau,amnuta,amel,amnue,ammu,amnumu
350 * ,ampiz,ampi,amro,gamro,ama1,gama1
351 * ,amk,amkz,amkst,gamkst
352 COMMON / taubra / gamprt(30),jlist(30),nchan
353 COMMON / taukle / bra1,brk0,brk0b,brks
354 real*4 bra1,brk0,brk0b,brks
356 parameter(nmode=15,nm1=0,nm2=1,nm3=8,nm4=2,nm5=1,nm6=3)
357 COMMON / taudcd /idffin(9,nmode),mulpik(nmode)
359 CHARACTER NAMES(NMODE)*31
361 parameter(nmode=15,nm1=0,nm2=1,nm3=8,nm4=2,nm5=1,nm6=3)
362 COMMON / taudcd /idffin(9,nmode),mulpik(nmode)
364 CHARACTER NAMES(NMODE)*31
366 CHARACTER OLDNAMES(7)*31
369 $ bxinit =
'(1x,1h*,g17.8, 16x, a31,a4,a4, 1x,1h*)'
440 dimension nopik(6,nmode),npik(nmode)
451 DATA nopik / -1,-1, 2, 1, 0, 0, 2, 2, 2,-1, 0, 0,
452 1 -1,-1, 1, 2, 2, 0, -1,-1,-1, 1, 1, 0,
453 2 -1,-1,-1, 1, 1, 2, -1,-1, 1, 2, 2, 2,
454 3 -3,-1, 3, 0, 0, 0, 4,-1,-4, 0, 0, 0,
455 4 -3, 2, 4, 0, 0, 0, 2, 2,-3, 0, 0, 0,
456 5 -3,-1, 1, 0, 0, 0, -1,-4, 2, 0, 0, 0,
457 6 9,-1, 2, 0, 0, 0, -1, 2, 8, 0, 0, 0,
459 DATA nopik / -1,-1, 1, 2, 0, 0, 2, 2, 2,-1, 0, 0,
460 1 -1,-1, 1, 2, 2, 0, -1,-1,-1, 1, 1, 0,
461 2 -1,-1,-1, 1, 1, 2, -1,-1, 1, 2, 2, 2,
462 3 -3,-1, 3, 0, 0, 0, -4,-1, 4, 0, 0, 0,
463 4 -3, 2,-4, 0, 0, 0, 2, 2,-3, 0, 0, 0,
464 5 -3,-1, 1, 0, 0, 0, -1, 4, 2, 0, 0, 0,
465 6 9,-1, 2, 0, 0, 0, -1, 2, 8, 0, 0, 0,
469 7 -3,-4, 0, 0, 0, 0 /
471 7 -3, 4, 0, 0, 0, 0 /
479 IF(i.EQ. 1) gamprt(i) = 1.0000
480 IF(i.EQ. 2) gamprt(i) = 1.0000
481 IF(i.EQ. 3) gamprt(i) = 1.0000
482 IF(i.EQ. 4) gamprt(i) = 1.0000
483 IF(i.EQ. 5) gamprt(i) = 1.0000
484 IF(i.EQ. 6) gamprt(i) = 1.0000
485 IF(i.EQ. 7) gamprt(i) = 1.0000
486 IF(i.EQ. 8) gamprt(i) = 1.0000
487 IF(i.EQ. 9) gamprt(i) = 1.0000
488 IF(i.EQ.10) gamprt(i) = 1.0000
489 IF(i.EQ.11) gamprt(i) = 1.0000
490 IF(i.EQ.12) gamprt(i) = 1.0000
491 IF(i.EQ.13) gamprt(i) = 1.0000
492 IF(i.EQ.14) gamprt(i) = 1.0000
493 IF(i.EQ.15) gamprt(i) = 1.0000
494 IF(i.EQ.16) gamprt(i) = 1.0000
495 IF(i.EQ.17) gamprt(i) = 1.0000
496 IF(i.EQ.18) gamprt(i) = 1.0000
497 IF(i.EQ.19) gamprt(i) = 1.0000
498 IF(i.EQ.20) gamprt(i) = 1.0000
499 IF(i.EQ.21) gamprt(i) = 1.0000
500 IF(i.EQ.22) gamprt(i) = 1.0000
502 IF(i.EQ. 1) gamprt(i) =0.1800
503 IF(i.EQ. 2) gamprt(i) =0.1751
504 IF(i.EQ. 3) gamprt(i) =0.1110
505 IF(i.EQ. 4) gamprt(i) =0.2515
506 IF(i.EQ. 5) gamprt(i) =0.1790
507 IF(i.EQ. 6) gamprt(i) =0.0071
508 IF(i.EQ. 7) gamprt(i) =0.0134
509 IF(i.EQ. 8) gamprt(i) =0.0450
510 IF(i.EQ. 9) gamprt(i) =0.0100
511 IF(i.EQ.10) gamprt(i) =0.0009
512 IF(i.EQ.11) gamprt(i) =0.0004
513 IF(i.EQ.12) gamprt(i) =0.0003
514 IF(i.EQ.13) gamprt(i) =0.0005
515 IF(i.EQ.14) gamprt(i) =0.0015
516 IF(i.EQ.15) gamprt(i) =0.0015
517 IF(i.EQ.16) gamprt(i) =0.0015
518 IF(i.EQ.17) gamprt(i) =0.0005
519 IF(i.EQ.18) gamprt(i) =0.0050
520 IF(i.EQ.19) gamprt(i) =0.0055
521 IF(i.EQ.20) gamprt(i) =0.0017
522 IF(i.EQ.21) gamprt(i) =0.0013
523 IF(i.EQ.22) gamprt(i) =0.0010
524 #elif defined (ALEPH)
525 IF(i.EQ. 1) gamprt(i) = 1.0000
526 IF(i.EQ. 2) gamprt(i) = .9732
527 IF(i.EQ. 3) gamprt(i) = .6217
528 IF(i.EQ. 4) gamprt(i) = 1.4221
529 IF(i.EQ. 5) gamprt(i) = 1.0180
530 IF(i.EQ. 6) gamprt(i) = .0405
531 IF(i.EQ. 7) gamprt(i) = .0781
532 IF(i.EQ. 8) gamprt(i) = .2414
533 IF(i.EQ. 9) gamprt(i) = .0601
534 IF(i.EQ.10) gamprt(i) = .0281
535 IF(i.EQ.11) gamprt(i) = .0045
536 IF(i.EQ.12) gamprt(i) = .0010
537 IF(i.EQ.13) gamprt(i) = .0062
538 IF(i.EQ.14) gamprt(i) = .0096
539 IF(i.EQ.15) gamprt(i) = .0169
540 IF(i.EQ.16) gamprt(i) = .0056
541 IF(i.EQ.17) gamprt(i) = .0045
542 IF(i.EQ.18) gamprt(i) = .0219
543 IF(i.EQ.19) gamprt(i) = .0180
544 IF(i.EQ.20) gamprt(i) = .0096
545 IF(i.EQ.21) gamprt(i) = .0088
546 IF(i.EQ.22) gamprt(i) = .0146
549 IF(i.EQ. 1) oldnames(i)=
' TAU- --> E- '
550 IF(i.EQ. 2) oldnames(i)=
' TAU- --> MU- '
551 IF(i.EQ. 3) oldnames(i)=
' TAU- --> PI- '
552 IF(i.EQ. 4) oldnames(i)=
' TAU- --> PI-, PI0 '
553 IF(i.EQ. 5) oldnames(i)=
' TAU- --> A1- (two subch) '
554 IF(i.EQ. 6) oldnames(i)=
' TAU- --> K- '
555 IF(i.EQ. 7) oldnames(i)=
' TAU- --> K*- (two subch) '
556 IF(i.EQ. 8) names(i-7)=
' TAU- --> 2PI-, PI0, PI+ '
557 IF(i.EQ. 9) names(i-7)=
' TAU- --> 3PI0, PI- '
558 IF(i.EQ.10) names(i-7)=
' TAU- --> 2PI-, PI+, 2PI0 '
559 IF(i.EQ.11) names(i-7)=
' TAU- --> 3PI-, 2PI+, '
560 IF(i.EQ.12) names(i-7)=
' TAU- --> 3PI-, 2PI+, PI0 '
561 IF(i.EQ.13) names(i-7)=
' TAU- --> 2PI-, PI+, 3PI0 '
562 IF(i.EQ.14) names(i-7)=
' TAU- --> K-, PI-, K+ '
563 IF(i.EQ.15) names(i-7)=
' TAU- --> K0, PI-, K0B '
565 IF(i.EQ.16) names(i-7)=
' TAU- --> K- PI0 K0 '
567 IF(i.EQ.16) names(i-7)=
' TAU- --> K-, K0, PI0 '
569 IF(i.EQ.17) names(i-7)=
' TAU- --> PI0 PI0 K- '
570 IF(i.EQ.18) names(i-7)=
' TAU- --> K- PI- PI+ '
571 IF(i.EQ.19) names(i-7)=
' TAU- --> PI- K0B PI0 '
572 IF(i.EQ.20) names(i-7)=
' TAU- --> ETA PI- PI0 '
573 IF(i.EQ.21) names(i-7)=
' TAU- --> PI- PI0 GAM '
574 IF(i.EQ.22) names(i-7)=
' TAU- --> K- K0 '
583 idffin(j,i)=nopik(j,i)
610 scabib = sqrt(1.-ccabib**2)
612 gamel = gfermi**2*amtau**5/(192*pi**3)
618 FUNCTION dcdmas(IDENT)
619 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
620 * ,ampiz,ampi,amro,gamro,ama1,gama1
621 * ,amk,amkz,amkst,gamkst
623 real*4 amtau,amnuta,amel,amnue,ammu,amnumu
624 * ,ampiz,ampi,amro,gamro,ama1,gama1
625 * ,amk,amkz,amkst,gamkst
626 IF (ident.EQ. 1)
THEN
628 ELSEIF (ident.EQ.-1)
THEN
630 ELSEIF (ident.EQ. 2)
THEN
632 ELSEIF (ident.EQ.-2)
THEN
634 ELSEIF (ident.EQ. 3)
THEN
636 ELSEIF (ident.EQ.-3)
THEN
638 ELSEIF (ident.EQ. 4)
THEN
640 ELSEIF (ident.EQ.-4)
THEN
642 ELSEIF (ident.EQ. 8)
THEN
644 ELSEIF (ident.EQ.-8)
THEN
646 ELSEIF (ident.EQ. 9)
THEN
648 ELSEIF (ident.EQ.-9)
THEN
651 print *,
'STOP IN APKMAS, WRONG IDENT=',ident
656 FUNCTION lunpik(ID,ISGN)
657 COMMON / taukle / bra1,brk0,brk0b,brks
658 real*4 bra1,brk0,brk0b,brks
662 IF (ident.EQ. 1)
THEN
664 ELSEIF (ident.EQ.-1)
THEN
666 ELSEIF (ident.EQ. 2)
THEN
668 ELSEIF (ident.EQ.-2)
THEN
670 ELSEIF (ident.EQ. 3)
THEN
672 ELSEIF (ident.EQ.-3)
THEN
675 IF (ident.EQ. 1)
THEN
677 ELSEIF (ident.EQ.-1)
THEN
679 ELSEIF (ident.EQ. 2)
THEN
681 ELSEIF (ident.EQ.-2)
THEN
683 ELSEIF (ident.EQ. 3)
THEN
685 ELSEIF (ident.EQ.-3)
THEN
688 ELSEIF (ident.EQ. 4)
THEN
692 IF (xio(1).GT.brk0)
THEN
697 ELSEIF (ident.EQ.-4)
THEN
701 IF (xio(1).GT.brk0b)
THEN
706 ELSEIF (ident.EQ. 8)
THEN
708 ELSEIF (ident.EQ.-8)
THEN
710 ELSEIF (ident.EQ. 9)
THEN
712 ELSEIF (ident.EQ.-9)
THEN
715 print *,
'STOP IN IPKDEF, WRONG IDENT=',ident
724 SUBROUTINE taurdf(KTO)
728 COMMON / taukle / bra1,brk0,brk0b,brks
729 real*4 bra1,brk0,brk0b,brks
730 COMMON / taubra / gamprt(30),jlist(30),nchan
750 SUBROUTINE taurdf(KTO)
754 COMMON / taukle / bra1,brk0,brk0b,brks
755 real*4 bra1,brk0,brk0b,brks
756 COMMON / taubra / gamprt(30),jlist(30),nchan
764 IF(i.EQ. 1) gamprt(i) = .0000
765 IF(i.EQ. 2) gamprt(i) = .0000
766 IF(i.EQ. 3) gamprt(i) = .0000
767 IF(i.EQ. 4) gamprt(i) = .0000
768 IF(i.EQ. 5) gamprt(i) = .0000
769 IF(i.EQ. 6) gamprt(i) = .0000
770 IF(i.EQ. 7) gamprt(i) = .0000
771 IF(i.EQ. 8) gamprt(i) = 1.0000
772 IF(i.EQ. 9) gamprt(i) = 1.0000
773 IF(i.EQ.10) gamprt(i) = 1.0000
774 IF(i.EQ.11) gamprt(i) = 1.0000
775 IF(i.EQ.12) gamprt(i) = 1.0000
776 IF(i.EQ.13) gamprt(i) = 1.0000
777 IF(i.EQ.14) gamprt(i) = 1.0000
778 IF(i.EQ.15) gamprt(i) = 1.0000
779 IF(i.EQ.16) gamprt(i) = 1.0000
780 IF(i.EQ.17) gamprt(i) = 1.0000
781 IF(i.EQ.18) gamprt(i) = 1.0000
782 IF(i.EQ.19) gamprt(i) = 1.0000
807 IF(i.EQ. 1) gamprt(i) = .0000
808 IF(i.EQ. 2) gamprt(i) = .0000
809 IF(i.EQ. 3) gamprt(i) = .0000
810 IF(i.EQ. 4) gamprt(i) = .0000
811 IF(i.EQ. 5) gamprt(i) = .0000
812 IF(i.EQ. 6) gamprt(i) = .0000
813 IF(i.EQ. 7) gamprt(i) = .0000
814 IF(i.EQ. 8) gamprt(i) = 1.0000
815 IF(i.EQ. 9) gamprt(i) = 1.0000
816 IF(i.EQ.10) gamprt(i) = 1.0000
817 IF(i.EQ.11) gamprt(i) = 1.0000
818 IF(i.EQ.12) gamprt(i) = 1.0000
819 IF(i.EQ.13) gamprt(i) = 1.0000
820 IF(i.EQ.14) gamprt(i) = 1.0000
821 IF(i.EQ.15) gamprt(i) = 1.0000
822 IF(i.EQ.16) gamprt(i) = 1.0000
823 IF(i.EQ.17) gamprt(i) = 1.0000
824 IF(i.EQ.18) gamprt(i) = 1.0000
825 IF(i.EQ.19) gamprt(i) = 1.0000
848 SUBROUTINE iniphy(XK00)
853 COMMON / qedprm /alfinv,alfpi,xk0
854 real*8 alfinv,alfpi,xk0
857 pi8 = 4.d0*datan(1.d0)
859 alfpi = 1d0/(alfinv*pi8)
869 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
870 * ,ampiz,ampi,amro,gamro,ama1,gama1
871 * ,amk,amkz,amkst,gamkst
873 real*4 amtau,amnuta,amel,amnue,ammu,amnumu
874 * ,ampiz,ampi,amro,gamro,ama1,gama1
875 * ,amk,amkz,amkst,gamkst
931 #elif defined (ALEPH)
940 print *,
'INIMAS a1 mass= ',ama1,gama1
954 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
955 * ,ampiz,ampi,amro,gamro,ama1,gama1
956 * ,amk,amkz,amkst,gamkst
958 real*4 amtau,amnuta,amel,amnue,ammu,amnumu
959 * ,ampiz,ampi,amro,gamro,ama1,gama1
960 * ,amk,amkz,amkst,gamkst
964 COMMON /taupos / npa,npb
965 dimension xpb1(4),xpb2(4),aqf1(4),aqf2(4)
977 CALL tralo4(1,aqf1,aqf1,am)
978 CALL tralo4(2,aqf2,aqf2,am)
980 kfb1= 11*idff/iabs(idff)
981 kfb2=-11*idff/iabs(idff)
985 $ xpb1(3)= aqf1(4)*aqf1(3)/abs(aqf1(3))
989 $ xpb2(3)= aqf2(4)*aqf2(3)/abs(aqf2(3))
994 CALL filhep( 1,3, kfb1,0,0,0,0,xpb1, amel,.true.)
995 CALL filhep( 2,3, kfb2,0,0,0,0,xpb2, amel,.true.)
996 CALL filhep(npa,1, idff,1,2,0,0,aqf1,amtau,.true.)
997 CALL filhep(npb,1,-idff,1,2,0,0,aqf2,amtau,.true.)
999 SUBROUTINE tralo4(KTO,P,Q,AM)
1004 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
1005 * ,ampiz,ampi,amro,gamro,ama1,gama1
1006 * ,amk,amkz,amkst,gamkst
1008 real*4 amtau,amnuta,amel,amnue,ammu,amnumu
1009 * ,ampiz,ampi,amro,gamro,ama1,gama1
1010 * ,amk,amkz,amkst,gamkst
1013 etau=sqrt(ptau**2+amtau**2)
1014 exe=(etau+ptau)/amtau
1015 IF(kto.EQ.2) exe=(etau-ptau)/amtau
1016 CALL bostr3(exe,p,q)
1021 SUBROUTINE filhep(N,IST,ID,JMO1,JMO2,JDA1,JDA2,P4,PINV,PHFLAG)
1031 #include "../../include/HEPEVT.h"
1048 ELSE IF (n.GT.0)
THEN
1059 IF ((ihep.LE.0).OR.(ihep.GT.nmxhep))
RETURN
1066 IF(jmo1.LT.0)jmohep(1,ihep)=jmohep(1,ihep)+ihep
1068 IF(jmo2.LT.0)jmohep(2,ihep)=jmohep(2,ihep)+ihep
1083 DO ip=jmohep(1,ihep),jmohep(2,ihep)
1087 IF(isthep(ip).EQ.1)isthep(ip)=2
1090 IF(jdahep(1,ip).EQ.0)
THEN
1094 jdahep(2,ip)=max(ihep,jdahep(2,ip))
1104 IMPLICIT REAL*8(a-h,o-z)
1107 IF(x .LT.-1.0)
GO TO 1
1108 IF(x .LE. 0.5)
GO TO 2
1109 IF(x .EQ. 1.0)
GO TO 3
1110 IF(x .LE. 2.0)
GO TO 4
1114 z=z-0.5* log(abs(x))**2
1120 3 dilogy=1.64493406684822
1124 z=1.64493406684822 - log(x)* log(abs(t))
1125 5 y=2.66666666666666 *t+0.66666666666666
1126 b= 0.00000 00000 00001
1127 a=y*b +0.00000 00000 00004
1128 b=y*a-b+0.00000 00000 00011
1129 a=y*b-a+0.00000 00000 00037
1130 b=y*a-b+0.00000 00000 00121
1131 a=y*b-a+0.00000 00000 00398
1132 b=y*a-b+0.00000 00000 01312
1133 a=y*b-a+0.00000 00000 04342
1134 b=y*a-b+0.00000 00000 14437
1135 a=y*b-a+0.00000 00000 48274
1136 b=y*a-b+0.00000 00001 62421
1137 a=y*b-a+0.00000 00005 50291
1138 b=y*a-b+0.00000 00018 79117
1139 a=y*b-a+0.00000 00064 74338
1140 b=y*a-b+0.00000 00225 36705
1141 a=y*b-a+0.00000 00793 87055
1142 b=y*a-b+0.00000 02835 75385
1143 a=y*b-a+0.00000 10299 04264
1144 b=y*a-b+0.00000 38163 29463
1145 a=y*b-a+0.00001 44963 00557
1146 b=y*a-b+0.00005 68178 22718
1147 a=y*b-a+0.00023 20021 96094
1148 b=y*a-b+0.00100 16274 96164
1149 a=y*b-a+0.00468 63619 59447
1150 b=y*a-b+0.02487 93229 24228
1151 a=y*b-a+0.16607 30329 27855
1152 a=y*a-b+1.93506 43008 6996