10 COMMON / / blan(10000)
12 COMMON / inout / inut,iout
18 OPEN(iout,file=
"./tauola.output")
19 OPEN( 16,file=
"./tauola.lund")
20 OPEN(inut,file=
"./dane.dat")
26 SUBROUTINE dectes(KTORY)
29 DOUBLE PRECISION hh(4)
31 COMMON / jaki / jak1,jak2,jakp,jakm,ktom
34 COMMON / inout / inut,iout
39 COMMON / taurad / xk0dec,itdkrc
62 READ( ninp,3000) testit
63 WRITE(nout,3000) testit
64 READ( ninp,3001) kat1,kat2,kat3,kat4,kat5,kat6
65 READ( ninp,3002) nevt,jak1,jak2,itdkrc
66 READ( ninp,3003) ptau,xk0dec
70 WRITE(nout,
'(6A6/6I6)')
71 $
'KAT1',
'KAT2',
'KAT3',
'KAT4',
'KAT5',
'KAT6',
72 $ kat1 , kat2 , kat3 , kat4 , kat5 , kat6
73 WRITE(nout,
'(4A12/4I12)')
74 $
'NEVT',
'JAK1',
'JAK2',
'ITDKRC',
75 $ nevt, jak1 , jak2 , itdkrc
76 WRITE(nout,
'(2A12/2F12.6)')
93 print *,
'for the sake of these tests KTO has to be 2'
94 print *,
'to change tau- to tau+ change IDFF from -15 to 15'
106 print *,
'NEVTES= ',nevtes
107 WRITE(iout,7011) keya1
110 WRITE(iout,7001) jak,idff,pol(3),ptau
112 WRITE(iout,7004) jak,idff,pol(3),ptau
133 #if defined (history)
143 CALL dekay(kto+10,hh)
155 IF(ipri.EQ.1) print *,
' event no: ',nev,
' NEVTES: ',nevtes
167 7001
FORMAT(//4(/1x,15(5h=====))
168 $ /,
' ', 19x,
' TEST OF RAD. CORR IN ELECTRON DECAY ',9x,1h ,
169 $ /,
' ', 19x,
' TESTS OF TAU DECAY ROUTINES ',9x,1h ,
170 $ /,
' ', 19x,
' INTERFACE OF THE KORAL-Z TYPE ',9x,1h ,
171 $ 2(/,1x,15(5h=====)),
172 $ /,5x ,
'JAK =',i7 ,
' KEY DEFINING DECAY TYPE ',9x,1h ,
173 $ /,5x ,
'IDFF =',i7 ,
' LUND IDENTIFIER FOR FIRST TAU ',9x,1h ,
174 $ /,5x ,
'POL(3)=',f7.2,
' THIRD COMPONENT OF TAU POLARIZ. ',9x,1h ,
175 $ /,5x ,
'PTAU =',f7.2,
' THIRD COMPONENT OF TAU MOM. GEV ',9x,1h ,
176 $ 2(/,1x,15(5h=====))/)
177 7002
FORMAT(///1x,
'===== EVENT NO.',i4,1x,5h=====)
178 7003
FORMAT(5x,
'POLARIMETRIC VECTOR: ',
179 $ 7x,
'HH(1)',7x,
'HH(2)',7x,
'HH(3)',7x,
'HH(4)',
180 $ /, 5x,
' ', 4(1x,f11.8) )
181 7004
FORMAT(//4(/1x,15(5h=====))
182 $ /,
' ', 19x,
' TEST OF RAD. CORR IN ELECTRON DECAY ',9x,1h ,
183 $ /,
' ', 19x,
' TESTS OF TAU DECAY ROUTINES ',9x,1h ,
184 $ /,
' ', 19x,
' INTERFACE OF THE KORAL-B TYPE ',9x,1h ,
185 $ 2(/,1x,15(5h=====)),
186 $ /,5x ,
'JAK =',i7 ,
' KEY DEFINING DECAY TYPE ',9x,1h ,
187 $ /,5x ,
'IDFF =',i7 ,
' LUND IDENTIFIER FOR FIRST TAU ',9x,1h ,
188 $ /,5x ,
'POL(3)=',f7.2,
' THIRD COMPONENT OF TAU POLARIZ. ',9x,1h ,
189 $ /,5x ,
'PTAU =',f7.2,
' THIRD COMPONENT OF TAU MOM. GEV ',9x,1h ,
190 $ 2(/,1x,15(5h=====))/)
191 7011
FORMAT(///1x,
'===== TYPE OF CURRENT',i4,1x,5h=====)
193 SUBROUTINE choice(MNUM,RR,ICHAN,PROB1,PROB2,PROB3,
194 $ amrx,gamrx,amra,gamra,amrb,gamrb)
195 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
196 * ,ampiz,ampi,amro,gamro,ama1,gama1
197 * ,amk,amkz,amkst,gamkst
199 REAL*4 amtau,amnuta,amel,amnue,ammu,amnumu
200 * ,ampiz,ampi,amro,gamro,ama1,gama1
201 * ,amk,amkz,amkst,gamkst
217 ELSEIF(mnum.EQ.1)
THEN
226 ELSEIF(mnum.EQ.2)
THEN
235 ELSEIF(mnum.EQ.3)
THEN
244 ELSEIF(mnum.EQ.4)
THEN
253 ELSEIF(mnum.EQ.5)
THEN
262 ELSEIF(mnum.EQ.6)
THEN
271 ELSEIF(mnum.EQ.7)
THEN
280 ELSEIF(mnum.EQ.8)
THEN
289 ELSEIF(mnum.EQ.101)
THEN
298 ELSEIF(mnum.EQ.102)
THEN
318 IF (rr.LE.prob1)
THEN
320 ELSEIF(rr.LE.(prob1+prob2))
THEN
335 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 / decomp /idffin(9,nmode),mulpik(nmode)
364 CHARACTER names(nmode)*31
434 dimension nopik(6,nmode),npik(nmode)
445 DATA nopik / -1,-1, 2, 1, 0, 0, 2, 2, 2,-1, 0, 0,
446 1 -1,-1, 1, 2, 2, 0, -1,-1,-1, 1, 1, 0,
447 2 -1,-1,-1, 1, 1, 2, -1,-1, 1, 2, 2, 2,
448 3 -3,-1, 3, 0, 0, 0, 4,-1,-4, 0, 0, 0,
449 4 -3, 2, 4, 0, 0, 0, 2, 2,-3, 0, 0, 0,
450 5 -3,-1, 1, 0, 0, 0, -1,-4, 2, 0, 0, 0,
451 6 9,-1, 2, 0, 0, 0, -1, 2, 8, 0, 0, 0,
453 DATA nopik / -1,-1, 1, 2, 0, 0, 2, 2, 2,-1, 0, 0,
454 1 -1,-1, 1, 2, 2, 0, -1,-1,-1, 1, 1, 0,
455 2 -1,-1,-1, 1, 1, 2, -1,-1, 1, 2, 2, 2,
456 3 -3,-1, 3, 0, 0, 0, -4,-1, 4, 0, 0, 0,
457 4 -3, 2,-4, 0, 0, 0, 2, 2,-3, 0, 0, 0,
458 5 -3,-1, 1, 0, 0, 0, -1, 4, 2, 0, 0, 0,
459 6 9,-1, 2, 0, 0, 0, -1, 2, 8, 0, 0, 0,
463 7 -3,-4, 0, 0, 0, 0 /
465 7 -3, 4, 0, 0, 0, 0 /
472 IF(i.EQ. 1) gamprt(i) = 1.0000
473 IF(i.EQ. 2) gamprt(i) = 1.0000
474 IF(i.EQ. 3) gamprt(i) = 1.0000
475 IF(i.EQ. 4) gamprt(i) = 1.0000
476 IF(i.EQ. 5) gamprt(i) = 1.0000
477 IF(i.EQ. 6) gamprt(i) = 1.0000
478 IF(i.EQ. 7) gamprt(i) = 1.0000
479 IF(i.EQ. 8) gamprt(i) = 1.0000
480 IF(i.EQ. 9) gamprt(i) = 1.0000
481 IF(i.EQ.10) gamprt(i) = 1.0000
482 IF(i.EQ.11) gamprt(i) = 1.0000
483 IF(i.EQ.12) gamprt(i) = 1.0000
484 IF(i.EQ.13) gamprt(i) = 1.0000
485 IF(i.EQ.14) gamprt(i) = 1.0000
486 IF(i.EQ.15) gamprt(i) = 1.0000
487 IF(i.EQ.16) gamprt(i) = 1.0000
488 IF(i.EQ.17) gamprt(i) = 1.0000
489 IF(i.EQ.18) gamprt(i) = 1.0000
490 IF(i.EQ.19) gamprt(i) = 1.0000
491 IF(i.EQ.20) gamprt(i) = 1.0000
492 IF(i.EQ.21) gamprt(i) = 1.0000
493 IF(i.EQ.22) gamprt(i) = 1.0000
495 IF(i.EQ. 1) gamprt(i) = 1.0000
496 IF(i.EQ. 2) gamprt(i) = 1.0000
497 IF(i.EQ. 3) gamprt(i) = 1.0000
498 IF(i.EQ. 4) gamprt(i) = 1.0000
499 IF(i.EQ. 5) gamprt(i) = 1.0000
500 IF(i.EQ. 6) gamprt(i) = 1.0000
501 IF(i.EQ. 7) gamprt(i) = 1.0000
502 IF(i.EQ. 8) gamprt(i) = 1.0000
503 IF(i.EQ. 9) gamprt(i) = 1.0000
504 IF(i.EQ.10) gamprt(i) = 1.0000
505 IF(i.EQ.11) gamprt(i) = 1.0000
506 IF(i.EQ.12) gamprt(i) = 1.0000
507 IF(i.EQ.13) gamprt(i) = 1.0000
508 IF(i.EQ.14) gamprt(i) = 1.0000
509 IF(i.EQ.15) gamprt(i) = 1.0000
510 IF(i.EQ.16) gamprt(i) = 1.0000
511 IF(i.EQ.17) gamprt(i) = 1.0000
512 IF(i.EQ.18) gamprt(i) = 1.0000
513 IF(i.EQ.19) gamprt(i) = 1.0000
514 IF(i.EQ.20) gamprt(i) = 1.0000
515 IF(i.EQ.21) gamprt(i) = 1.0000
516 IF(i.EQ.22) gamprt(i) = 1.0000
518 IF(i.EQ. 1) gamprt(i) =0.1800
519 IF(i.EQ. 2) gamprt(i) =0.1751
520 IF(i.EQ. 3) gamprt(i) =0.1110
521 IF(i.EQ. 4) gamprt(i) =0.2515
522 IF(i.EQ. 5) gamprt(i) =0.1790
523 IF(i.EQ. 6) gamprt(i) =0.0071
524 IF(i.EQ. 7) gamprt(i) =0.0134
525 IF(i.EQ. 8) gamprt(i) =0.0450
526 IF(i.EQ. 9) gamprt(i) =0.0100
527 IF(i.EQ.10) gamprt(i) =0.0009
528 IF(i.EQ.11) gamprt(i) =0.0004
529 IF(i.EQ.12) gamprt(i) =0.0003
530 IF(i.EQ.13) gamprt(i) =0.0005
531 IF(i.EQ.14) gamprt(i) =0.0015
532 IF(i.EQ.15) gamprt(i) =0.0015
533 IF(i.EQ.16) gamprt(i) =0.0015
534 IF(i.EQ.17) gamprt(i) =0.0005
535 IF(i.EQ.18) gamprt(i) =0.0050
536 IF(i.EQ.19) gamprt(i) =0.0055
537 IF(i.EQ.20) gamprt(i) =0.0017
538 IF(i.EQ.21) gamprt(i) =0.0013
539 IF(i.EQ.22) gamprt(i) =0.0010
540 #elif defined (ALEPH)
541 IF(i.EQ. 1) gamprt(i) = 1.0000
542 IF(i.EQ. 2) gamprt(i) = .9732
543 IF(i.EQ. 3) gamprt(i) = .6217
544 IF(i.EQ. 4) gamprt(i) = 1.4221
545 IF(i.EQ. 5) gamprt(i) = 1.0180
546 IF(i.EQ. 6) gamprt(i) = .0405
547 IF(i.EQ. 7) gamprt(i) = .0781
548 IF(i.EQ. 8) gamprt(i) = .2414
549 IF(i.EQ. 9) gamprt(i) = .0601
550 IF(i.EQ.10) gamprt(i) = .0281
551 IF(i.EQ.11) gamprt(i) = .0045
552 IF(i.EQ.12) gamprt(i) = .0010
553 IF(i.EQ.13) gamprt(i) = .0062
554 IF(i.EQ.14) gamprt(i) = .0096
555 IF(i.EQ.15) gamprt(i) = .0169
556 IF(i.EQ.16) gamprt(i) = .0056
557 IF(i.EQ.17) gamprt(i) = .0045
558 IF(i.EQ.18) gamprt(i) = .0219
559 IF(i.EQ.19) gamprt(i) = .0180
560 IF(i.EQ.20) gamprt(i) = .0096
561 IF(i.EQ.21) gamprt(i) = .0088
562 IF(i.EQ.22) gamprt(i) = .0146
565 IF(i.EQ. 8) names(i-7)=
' TAU- --> 2PI-, PI0, PI+ '
566 IF(i.EQ. 9) names(i-7)=
' TAU- --> 3PI0, PI- '
567 IF(i.EQ.10) names(i-7)=
' TAU- --> 2PI-, PI+, 2PI0 '
568 IF(i.EQ.11) names(i-7)=
' TAU- --> 3PI-, 2PI+, '
569 IF(i.EQ.12) names(i-7)=
' TAU- --> 3PI-, 2PI+, PI0 '
570 IF(i.EQ.13) names(i-7)=
' TAU- --> 2PI-, PI+, 3PI0 '
571 IF(i.EQ.14) names(i-7)=
' TAU- --> K-, PI-, K+ '
572 IF(i.EQ.15) names(i-7)=
' TAU- --> K0, PI-, K0B '
574 IF(i.EQ.16) names(i-7)=
' TAU- --> K- PI0 K0 '
576 IF(i.EQ.16) names(i-7)=
' TAU- --> K-, K0, PI0 '
578 IF(i.EQ.17) names(i-7)=
' TAU- --> PI0, PI0, K- '
579 IF(i.EQ.18) names(i-7)=
' TAU- --> K-, PI-, PI+ '
580 IF(i.EQ.19) names(i-7)=
' TAU- --> PI-, K0B, PI0 '
581 IF(i.EQ.20) names(i-7)=
' TAU- --> ETA, PI-, PI0 '
582 IF(i.EQ.21) names(i-7)=
' TAU- --> PI-, PI0, GAM '
583 IF(i.EQ.22) names(i-7)=
' TAU- --> K-, K0 '
592 idffin(j,i)=nopik(j,i)
617 scabib = sqrt(1.-ccabib**2)
618 gamel = gfermi**2*amtau**5/(192*pi**3)
624 FUNCTION dcdmas(IDENT)
625 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
626 * ,ampiz,ampi,amro,gamro,ama1,gama1
627 * ,amk,amkz,amkst,gamkst
629 REAL*4 amtau,amnuta,amel,amnue,ammu,amnumu
630 * ,ampiz,ampi,amro,gamro,ama1,gama1
631 * ,amk,amkz,amkst,gamkst
632 IF (ident.EQ. 1)
THEN
634 ELSEIF (ident.EQ.-1)
THEN
636 ELSEIF (ident.EQ. 2)
THEN
638 ELSEIF (ident.EQ.-2)
THEN
640 ELSEIF (ident.EQ. 3)
THEN
642 ELSEIF (ident.EQ.-3)
THEN
644 ELSEIF (ident.EQ. 4)
THEN
646 ELSEIF (ident.EQ.-4)
THEN
648 ELSEIF (ident.EQ. 8)
THEN
650 ELSEIF (ident.EQ.-8)
THEN
652 ELSEIF (ident.EQ. 9)
THEN
654 ELSEIF (ident.EQ.-9)
THEN
657 print *,
'STOP IN APKMAS, WRONG IDENT=',ident
662 FUNCTION lunpik(ID,ISGN)
663 COMMON / taukle / bra1,brk0,brk0b,brks
664 REAL*4 bra1,brk0,brk0b,brks
669 IF (ident.EQ. 1)
THEN
671 ELSEIF (ident.EQ.-1)
THEN
673 ELSEIF (ident.EQ. 2)
THEN
675 ELSEIF (ident.EQ.-2)
THEN
677 ELSEIF (ident.EQ. 3)
THEN
679 ELSEIF (ident.EQ.-3)
THEN
682 IF (ident.EQ. 1)
THEN
684 ELSEIF (ident.EQ.-1)
THEN
686 ELSEIF (ident.EQ. 2)
THEN
688 ELSEIF (ident.EQ.-2)
THEN
690 ELSEIF (ident.EQ. 3)
THEN
692 ELSEIF (ident.EQ.-3)
THEN
695 ELSEIF (ident.EQ. 4)
THEN
699 IF (xio(1).GT.brk0)
THEN
704 ELSEIF (ident.EQ.-4)
THEN
708 IF (xio(1).GT.brk0b)
THEN
713 ELSEIF (ident.EQ. 8)
THEN
715 ELSEIF (ident.EQ.-8)
THEN
717 ELSEIF (ident.EQ. 9)
THEN
719 ELSEIF (ident.EQ.-9)
THEN
722 print *,
'STOP IN IPKDEF, WRONG IDENT=',ident
729 SUBROUTINE taurdf(KTO)
733 COMMON / taukle / bra1,brk0,brk0b,brks
734 REAL*4 bra1,brk0,brk0b,brks
735 COMMON / taubra / gamprt(30),jlist(30),nchan
755 SUBROUTINE taurdf(KTO)
759 COMMON / taukle / bra1,brk0,brk0b,brks
760 REAL*4 bra1,brk0,brk0b,brks
761 COMMON / taubra / gamprt(30),jlist(30),nchan
769 IF(i.EQ. 1) gamprt(i) = .0000
770 IF(i.EQ. 2) gamprt(i) = .0000
771 IF(i.EQ. 3) gamprt(i) = .0000
772 IF(i.EQ. 4) gamprt(i) = .0000
773 IF(i.EQ. 5) gamprt(i) = .0000
774 IF(i.EQ. 6) gamprt(i) = .0000
775 IF(i.EQ. 7) gamprt(i) = .0000
776 IF(i.EQ. 8) gamprt(i) = 1.0000
777 IF(i.EQ. 9) gamprt(i) = 1.0000
778 IF(i.EQ.10) gamprt(i) = 1.0000
779 IF(i.EQ.11) gamprt(i) = 1.0000
780 IF(i.EQ.12) gamprt(i) = 1.0000
781 IF(i.EQ.13) gamprt(i) = 1.0000
782 IF(i.EQ.14) gamprt(i) = 1.0000
783 IF(i.EQ.15) gamprt(i) = 1.0000
784 IF(i.EQ.16) gamprt(i) = 1.0000
785 IF(i.EQ.17) gamprt(i) = 1.0000
786 IF(i.EQ.18) gamprt(i) = 1.0000
787 IF(i.EQ.19) gamprt(i) = 1.0000
812 IF(i.EQ. 1) gamprt(i) = .0000
813 IF(i.EQ. 2) gamprt(i) = .0000
814 IF(i.EQ. 3) gamprt(i) = .0000
815 IF(i.EQ. 4) gamprt(i) = .0000
816 IF(i.EQ. 5) gamprt(i) = .0000
817 IF(i.EQ. 6) gamprt(i) = .0000
818 IF(i.EQ. 7) gamprt(i) = .0000
819 IF(i.EQ. 8) gamprt(i) = 1.0000
820 IF(i.EQ. 9) gamprt(i) = 1.0000
821 IF(i.EQ.10) gamprt(i) = 1.0000
822 IF(i.EQ.11) gamprt(i) = 1.0000
823 IF(i.EQ.12) gamprt(i) = 1.0000
824 IF(i.EQ.13) gamprt(i) = 1.0000
825 IF(i.EQ.14) gamprt(i) = 1.0000
826 IF(i.EQ.15) gamprt(i) = 1.0000
827 IF(i.EQ.16) gamprt(i) = 1.0000
828 IF(i.EQ.17) gamprt(i) = 1.0000
829 IF(i.EQ.18) gamprt(i) = 1.0000
830 IF(i.EQ.19) gamprt(i) = 1.0000
852 SUBROUTINE iniphy(XK00)
857 COMMON / qedprm /alfinv,alfpi,xk0
858 REAL*8 alfinv,alfpi,xk0
861 pi8 = 4.d0*datan(1.d0)
863 alfpi = 1d0/(alfinv*pi8)
872 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
873 * ,ampiz,ampi,amro,gamro,ama1,gama1
874 * ,amk,amkz,amkst,gamkst
876 REAL*4 amtau,amnuta,amel,amnue,ammu,amnumu
877 * ,ampiz,ampi,amro,gamro,ama1,gama1
878 * ,amk,amkz,amkst,gamkst
944 #elif defined (ALEPH)
953 print *,
'INIMAS a1 mass= ',ama1,gama1
967 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
968 * ,ampiz,ampi,amro,gamro,ama1,gama1
969 * ,amk,amkz,amkst,gamkst
971 REAL*4 amtau,amnuta,amel,amnue,ammu,amnumu
972 * ,ampiz,ampi,amro,gamro,ama1,gama1
973 * ,amk,amkz,amkst,gamkst
977 COMMON /taupos / npa,npb
978 dimension xpb1(4),xpb2(4),aqf1(4),aqf2(4)
990 CALL tralo4(1,aqf1,aqf1,am)
991 CALL tralo4(2,aqf2,aqf2,am)
993 kfb1= 11*idff/iabs(idff)
994 kfb2=-11*idff/iabs(idff)
998 $ xpb1(3)= aqf1(4)*aqf1(3)/abs(aqf1(3))
1002 $ xpb2(3)= aqf2(4)*aqf2(3)/abs(aqf2(3))
1007 CALL filhep( 1,3, kfb1,0,0,0,0,xpb1, amel,.true.)
1008 CALL filhep( 2,3, kfb2,0,0,0,0,xpb2, amel,.true.)
1009 CALL filhep(npa,1, idff,1,2,0,0,aqf1,amtau,.true.)
1010 CALL filhep(npb,1,-idff,1,2,0,0,aqf2,amtau,.true.)
1012 SUBROUTINE tralo4(KTO,P,Q,AM)
1017 COMMON / parmas / amtau,amnuta,amel,amnue,ammu,amnumu
1018 * ,ampiz,ampi,amro,gamro,ama1,gama1
1019 * ,amk,amkz,amkst,gamkst
1021 REAL*4 amtau,amnuta,amel,amnue,ammu,amnumu
1022 * ,ampiz,ampi,amro,gamro,ama1,gama1
1023 * ,amk,amkz,amkst,gamkst
1026 etau=sqrt(ptau**2+amtau**2)
1027 exe=(etau+ptau)/amtau
1028 IF(kto.EQ.2) exe=(etau-ptau)/amtau
1029 CALL bostr3(exe,p,q)
1034 SUBROUTINE filhep(N,IST,ID,JMO1,JMO2,JDA1,JDA2,P4,PINV,PHFLAG)
1044 #include "../../include/HEPEVT.h"
1061 ELSE IF (n.GT.0)
THEN
1072 IF ((ihep.LE.0).OR.(ihep.GT.nmxhep))
RETURN
1079 IF(jmo1.LT.0)jmohep(1,ihep)=jmohep(1,ihep)+ihep
1081 IF(jmo2.LT.0)jmohep(2,ihep)=jmohep(2,ihep)+ihep
1096 DO ip=jmohep(1,ihep),jmohep(2,ihep)
1100 IF(isthep(ip).EQ.1)isthep(ip)=2
1103 IF(jdahep(1,ip).EQ.0)
THEN
1107 jdahep(2,ip)=max(ihep,jdahep(2,ip))
1117 IMPLICIT REAL*8(a-h,o-z)
1120 IF(x .LT.-1.0) go to 1
1121 IF(x .LE. 0.5) go to 2
1122 IF(x .EQ. 1.0) go to 3
1123 IF(x .LE. 2.0) go to 4
1127 z=z-0.5* log(abs(x))**2
1133 3 dilogy=1.64493406684822
1137 z=1.64493406684822 - log(x)* log(abs(t))
1138 5 y=2.66666666666666 *t+0.66666666666666
1139 b= 0.00000 00000 00001
1140 a=y*b +0.00000 00000 00004
1141 b=y*a-b+0.00000 00000 00011
1142 a=y*b-a+0.00000 00000 00037
1143 b=y*a-b+0.00000 00000 00121
1144 a=y*b-a+0.00000 00000 00398
1145 b=y*a-b+0.00000 00000 01312
1146 a=y*b-a+0.00000 00000 04342
1147 b=y*a-b+0.00000 00000 14437
1148 a=y*b-a+0.00000 00000 48274
1149 b=y*a-b+0.00000 00001 62421
1150 a=y*b-a+0.00000 00005 50291
1151 b=y*a-b+0.00000 00018 79117
1152 a=y*b-a+0.00000 00064 74338
1153 b=y*a-b+0.00000 00225 36705
1154 a=y*b-a+0.00000 00793 87055
1155 b=y*a-b+0.00000 02835 75385
1156 a=y*b-a+0.00000 10299 04264
1157 b=y*a-b+0.00000 38163 29463
1158 a=y*b-a+0.00001 44963 00557
1159 b=y*a-b+0.00005 68178 22718
1160 a=y*b-a+0.00023 20021 96094
1161 b=y*a-b+0.00100 16274 96164
1162 a=y*b-a+0.00468 63619 59447
1163 b=y*a-b+0.02487 93229 24228
1164 a=y*b-a+0.16607 30329 27855
1165 a=y*a-b+1.93506 43008 6996