IMPLICIT REAL*8(A-H,O-Z) DIMENSION RS(27),PO(9),RO(4,27) * NRS= 27 IC= 1 ICN= 4 ZM= 91.1867D0 TQM= 173.8D0 ALS= 0.119D0 HM= 100.D0 DO I=1,27 RS(I)= ZM-1.3D0+(I-1)*0.1D0 ENDDO PO(1)= ZM PO(2)= 0.249549E+01 PO(3)= 0.414761E+02 PO(4)= 0.207435E+02 PO(5)= 0.207437E+02 PO(6)= 0.207906E+02 PO(7)= 0.160835E-01 PO(8)= 0.160835E-01 PO(9)= 0.160835E-01 * CALL TTHEORY(IC,ICN,NRS,RS,ZM,TQM,ALS,HM,PO,RO) * DO I=1,NRS PRINT 1,RS(I),RO(1,I),RO(2,I),RO(3,I),RO(4,I) ENDDO 1 FORMAT(' E_CM (GEV) = ',E20.5// # ' SIGMA(E) (NB) = ',E20.7/ # ' SIGMA(MU) (NB) = ',E20.7/ # ' SIGMA(TAU) (NB) = ',E20.7/ # ' SIGMA(HAD) (NB) = ',E20.7///) * STOP END * *--------------------------------------------------------------------- * SUBROUTINE TTHEORY(IC,ICN,NRS,RS,ZMT,TQMT,ALST,HMT,PO,RO) IMPLICIT REAL*8(A-H,O-Z) * PARAMETER(NX=27,NY=28,NPO=34,NRO=26*NX,NFIT=7*NY) * CHARACTER*1 OEXT CHARACTER*4 OF4 CHARACTER*6 OF6 * COMMON/TFITO/APO(9) COMMON/TFIT/FIT(NFIT) * DIMENSION RS(NRS),PO(9),RO(ICN,NRS) DIMENSION OTPPO(NPO),OTPRO(NRO),OTPDRO(5,9,NX) DIMENSION SET(NY),AET(NY),SMUT(NY),AMUT(NY),STAUT(NY), # ATAUT(NY),SHADT(NY) * NRST= NRS+1 ZM= ZMT TQM= TQMT ALS= ALST HM= HMT * DO I=1,9 APO(I)= PO(I) ENDDO * OEXT= 'E' CALL TINIT(NRS,OEXT) OF4= 'FITE' OF6= 'OMODES' CALL TCFLAG(OF6,OF4) * SE= 1.D2 CALL TOPAZ0(SE,NRS,RS,ZM,TQM,HM,ALS,OTPPO,OTPRO,OTPDRO) * DO IT=1,7 DO IF=1,NRST K= NRST*(IT-1)+IF IF(IT.EQ.1) THEN SET(IF)= FIT(K) ELSE IF(IT.EQ.2) THEN SMUT(IF)= FIT(K) ELSE IF(IT.EQ.3) THEN STAUT(IF)= FIT(K) ELSE IF(IT.EQ.4) THEN SHADT(IF)= FIT(K) ELSE IF(IT.EQ.5) THEN AET(IF)= FIT(K) ELSE IF(IT.EQ.6) THEN AMUT(IF)= FIT(K) ELSE IF(IT.EQ.7) THEN ATAUT(IF)= FIT(K) ENDIF ENDDO ENDDO * IF(IC.EQ.1) THEN DO I=1,NRS RO(1,I)= SET(I) RO(2,I)= SMUT(I) RO(3,I)= STAUT(I) RO(4,I)= SHADT(I) ENDDO ELSE IF(IC.EQ.2) THEN DO I=1,NRS RO(1,I)= AET(I) RO(2,I)= AMUT(I) RO(3,I)= ATAUT(I) ENDDO ENDIF * RETURN END