IMPLICIT REAL*8(A-H,O-Z) PARAMETER(NX=5,NY=6,NPO=34,NRO=26*NX,NFIT=7*NY) CHARACTER*1 OEXT,OVAL,ORAD,OCHAN CHARACTER*2 OF2 CHARACTER*3 OF3 CHARACTER*4 OF4 CHARACTER*5 OF5 CHARACTER*6 OF6 CHARACTER*7 OF7 * COMMON/TK/IFK COMMON/TFITO/APO(9) COMMON/TFIT/FIT(NFIT) * DIMENSION RS(NX) DIMENSION OTPPO(NPO),OTPRO(NRO) DIMENSION SET(NY),AET(NY),SMUT(NY),AMUT(NY),STAUT(NY), # ATAUT(NY),SHADT(NY),ESET(NY),EAET(NY), # ESMUT(NY),EAMUT(NY),ESTAUT(NY),EATAUT(NY), # ESHADT(NY),ACT(NY),ABT(NY),EACT(NY),EABT(NY), # SBT(NY),ESBT(NY),SCT(NY),ESCT(NY),BRAT(NY), # EBRAT(NY),CRAT(NY),ECRAT(NY) * IFK= 0 NRS= 5 NRST= 6 ZM= 91.1867D0 TQM= 173.8D0 ALS= 0.119D0 HM= 100.D0 RS(1)= ZM-3.D0 RS(2)= ZM-1.8D0 RS(3)= ZM RS(4)= ZM+1.8D0 RS(5)= ZM+3.D0 * OEXT= 'E' CALL TINIT(NRS,OEXT) * SE= 1.D-1 CALL TOPAZ0(SE,NRS,RS,ZM,TQM,HM,ALS,OTPPO,OTPRO,OTPDRO) CALL TCINFO(0) CALL TCINFO(1) PRINT 1,(OTPPO(I),I=1,NPO) * IF(OEXT.NE.'P') THEN DO I=1,NRS DO J=1,2 DO L=1,13 K= 13*(2*(I-1)+J-1)+L IF(L.EQ.1) THEN IF(J.EQ.1) THEN SET(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN ESET(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.2) THEN IF(J.EQ.1) THEN SMUT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN ESMUT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.3) THEN IF(J.EQ.1) THEN STAUT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN ESTAUT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.4) THEN IF(J.EQ.1) THEN SHADT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN ESHADT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.5) THEN IF(J.EQ.1) THEN SBT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN ESBT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.6) THEN IF(J.EQ.1) THEN BRAT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN EBRAT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.7) THEN IF(J.EQ.1) THEN SCT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN ESCT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.8) THEN IF(J.EQ.1) THEN CRAT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN ECRAT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.9) THEN IF(J.EQ.1) THEN AET(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN EAET(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.10) THEN IF(J.EQ.1) THEN AMUT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN EAMUT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.11) THEN IF(J.EQ.1) THEN ATAUT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN EATAUT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.12) THEN IF(J.EQ.1) THEN ACT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN EACT(I)= OTPRO(K) ENDIF ELSE IF(L.EQ.13) THEN IF(J.EQ.1) THEN ABT(I)= OTPRO(K) ELSE IF(J.EQ.2) THEN EABT(I)= OTPRO(K) ENDIF ENDIF ENDDO ENDDO PRINT 3,RS(I),SET(I),SMUT(I),STAUT(I),SHADT(I), # SBT(I),SCT(I),AET(I),AMUT(I),ATAUT(I), # ACT(I),ABT(I) ENDDO ENDIF * 1 FORMAT( ' W MASS (GEV) = ',E20.6/ # ' NU = ',E20.6/ # ' ELECTRON = ',E20.6/ # ' MUON = ',E20.6/ # ' TAU = ',E20.6/ # ' UP = ',E20.6/ # ' DOWN(STRANGE) = ',E20.6/ # ' CHARM = ',E20.6/ # ' BOTTOM = ',E20.6/ # ' SIN^2(E) = ',E20.6/ # ' SIN^2(B) = ',E20.6/ # ' A_FB(L) EFF. = ',E20.6/ # ' A_LR EFF. = ',E20.6/ # ' TOTAL WIDTH (GEV) = ',E20.6/ # ' G_H/G_E = ',E20.6/ # ' SIGMA0_H (NB) = ',E20.6/ # ' G(B)/G(HAD) = ',E20.6/ # ' A_FB(B) = ',E20.6/ # ' HADRONIC WIDTH (GEV) = ',E20.6/ # ' INVISIBLE = ',E20.6/ # ' A_FB(C) = ',E20.6/ # ' G(C)/G(HAD) = ',E20.6/ # ' A_LR(B) EFF. = ',E20.6/ # ' A_LR(C) EFF. = ',E20.6/ # ' SIN^2(C) = ',E20.6/ # ' RHO(E) = ',E20.6/ # ' RHO(C) = ',E20.6/ # ' RHO(B) = ',E20.6/ # ' 1/ALPHA^5 = ',E20.6/ # ' 1/ALPHA = ',E20.6/ # ' SIGMA0_L (NB) = ',E20.6/ # ' ALPHA_S(M_T) = ',E20.6/ # ' G_H/G_M = ',E20.6/ # ' G_H/G_T = ',E20.6///) * 3 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/ # ' SIGMA(B) (NB) = ',E20.7/ # ' SIGMA(C) (NB) = ',E20.7/ # ' A_FB(E) = ',E20.7/ # ' A_FB(MU) = ',E20.7/ # ' A_FB(TAU) = ',E20.7/ # ' A_FB(C) = ',E20.7/ # ' A_FB(B) = ',E20.7///) * STOP END