IMPLICIT REAL*8(A-H,O-Z) PARAMETER(NX=5,NPO=32,NRO=26*NX) CHARACTER*1 OEXT,OVAL,ORAD,OCHAN CHARACTER*2 OV2 CHARACTER*3 OF3 CHARACTER*4 OF4 CHARACTER*5 OF5 CHARACTER*7 OF7 * DIMENSION RS(NX) DIMENSION OTPPO(NPO),OTPRO(NRO),OTPROD(9,NX),OTPRODD(9,NX), # OTPROZZ(9,NX),OTPROWI(9,NX) DIMENSION SET(NX),AET(NX),SMUT(NX),AMUT(NX),STAUT(NX), # ATAUT(NX),SHADT(NX),ESET(NX),EAET(NX), # ESMUT(NX),EAMUT(NX),ESTAUT(NX),EATAUT(NX), # ESHADT(NX),ACT(NX),ABT(NX),EACT(NX),EABT(NX), # SBT(NX),ESBT(NX),SCT(NX),ESCT(NX),BRAT(NX), # EBRAT(NX),CRAT(NX),ECRAT(NX),DECXSH(NX),DDECXSH(NX), # ZZECXSH(NX),WIECXSH(NX) * NRS= 5 ZM= 91.1867D0 TQM= 173.8D0 ALS= 0.119D0 HM= 76.D0 RS(1)= ZM-3.D0 RS(2)= ZM-1D0 RS(3)= ZM RS(4)= ZM+1.D0 RS(5)= ZM+3.D0 * OEXT= 'P' OEXT= 'E' CALL TINIT(NRS,OEXT) * OVAL= 'Y' * OF4= 'ORAD' * CALL TCFLAG(OF4,OVAL) * * OV2= 'FC' * OF5= 'OCUTF' * CALL TCFLAG(OF5,OV2) * * DMY= -1.D0 * IND= 1 * XCUTF= 0.01D0 * CALL TCUTSET(IND,DMY,DMY,DMY,DMY,DMY,DMY,XCUTF,DMY) * *---------------------------------------------------------------------- * SE= 1.D0 CALL TOPAZ0(SE,NRS,RS,ZM,TQM,HM,ALS,OTPPO,OTPRO,OTPROD, # OTPRODD,OTPROZZ,OTPROWI) * 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 2,RS(I),SET(I),ESET(I),SMUT(I),ESMUT(I), # STAUT(I),ESTAUT(I),SHADT(I),ESHADT(I), # SBT(I),ESBT(I),BRAT(I),EBRAT(I), # SCT(I),ESCT(I),CRAT(I),ECRAT(I),AET(I), # EAET(I),AMUT(I),EAMUT(I),ATAUT(I), # EATAUT(I),ACT(I),EACT(I),ABT(I),EABT(I) ENDDO ENDIF * IF(OEXT.EQ.'E') THEN DO I=1,NRS DECXSH(I)= 0.D0 DO J= 2,6 DECXSH(I)= DECXSH(I)+OTPROD(J,I) ENDDO PRINT 3,RS(I),OTPROD(1,I),OTPROD(2,I),OTPROD(3,I), # OTPROD(4,I),OTPROD(5,I),OTPROD(6,I), # DECXSH(I), # OTPROD(7,I),OTPROD(8,I),OTPROD(9,I) ENDDO DO I=1,NRS DDECXSH(I)= 0.D0 DO J= 2,6 DDECXSH(I)= DDECXSH(I)+OTPRODD(J,I) ENDDO PRINT 4,RS(I),OTPRODD(1,I),OTPRODD(2,I),OTPRODD(3,I), # OTPRODD(4,I),OTPRODD(5,I),OTPRODD(6,I), # DDECXSH(I), # OTPRODD(7,I),OTPRODD(8,I),OTPRODD(9,I) ENDDO DO I=1,NRS ZZECXSH(I)= 0.D0 DO J= 2,6 ZZECXSH(I)= ZZECXSH(I)+OTPROZZ(J,I) ENDDO PRINT 5,RS(I),OTPROZZ(1,I),OTPROZZ(2,I),OTPROZZ(3,I), # OTPROZZ(4,I),OTPROZZ(5,I),OTPROZZ(6,I), # ZZECXSH(I), # OTPROZZ(7,I),OTPROZZ(8,I),OTPROZZ(9,I) ENDDO DO I=1,NRS WIECXSH(I)= 0.D0 DO J= 2,6 WIECXSH(I)= WIECXSH(I)+OTPROWI(J,I) ENDDO PRINT 6,RS(I),OTPROWI(1,I),OTPROWI(2,I),OTPROWI(3,I), # OTPROWI(4,I),OTPROWI(5,I),OTPROWI(6,I), # WIECXSH(I), # OTPROWI(7,I),OTPROWI(8,I),OTPROWI(9,I) ENDDO ENDIF IF(OEXT.EQ.'E') THEN DO I=1,NRS PRINT 7,RS(I),1.D2*(SMUT(I)/OTPROD(1,I)-1.D0), # 1.D2*(SHADT(I)/DECXSH(I)-1.D0), # AMUT(I)-OTPROD(7,I) ENDDO ENDIF IF(OEXT.EQ.'E') THEN DO I=1,NRS PRINT 8,RS(I),1.D2*(OTPRODD(1,I)/OTPROWI(1,I)-1.D0), # 1.D2*(DDECXSH(I)/WIECXSH(I)-1.D0), # OTPRODD(7,I)-OTPROWI(7,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///) * 2 FORMAT(' E_CM (GEV) = ',E20.5// # ' SIGMA(E) (NB) = ',E20.7,' +/- ',E20.7/ # ' SIGMA(MU) (NB) = ',E20.7,' +/- ',E20.7/ # ' SIGMA(TAU) (NB) = ',E20.7,' +/- ',E20.7/ # ' SIGMA(HAD) (NB) = ',E20.7,' +/- ',E20.7/ # ' SIGMA(B) (NB) = ',E20.7,' +/- ',E20.7/ # ' R_B = ',E20.7,' +/- ',E20.7/ # ' SIGMA(C) (NB) = ',E20.7,' +/- ',E20.7/ # ' R_C = ',E20.7,' +/- ',E20.7/ # ' A_FB(E) = ',E20.7,' +/- ',E20.7/ # ' A_FB(MU) = ',E20.7,' +/- ',E20.7/ # ' A_FB(TAU) = ',E20.7,' +/- ',E20.7/ # ' A_FB(C) = ',E20.7,' +/- ',E20.7/ # ' A_FB(B) = ',E20.7,' +/- ',E20.7///) * 3 FORMAT(' E_CM (GEV) = ',E20.5// # ' DECONVOLUTED QUANTITIES '/ # ' SIGMA(MU) (NB) = ',E20.7/ # ' SIGMA(U) (NB) = ',E20.7/ # ' SIGMA(D) (NB) = ',E20.7/ # ' SIGMA(C) (NB) = ',E20.7/ # ' SIGMA(S) (NB) = ',E20.7/ # ' SIGMA(B) (NB) = ',E20.7/ # ' SIGMA(HAD) (NB) = ',E20.7/ # ' A_FB(MU) = ',E20.7/ # ' A_FB(C) = ',E20.7/ # ' A_FB(B) = ',E20.7///) * 4 FORMAT(' E_CM (GEV) = ',E20.5// # ' DDECONVOLUTED QUANTITIES '/ # ' SIGMA(MU) (NB) = ',E20.7/ # ' SIGMA(U) (NB) = ',E20.7/ # ' SIGMA(D) (NB) = ',E20.7/ # ' SIGMA(C) (NB) = ',E20.7/ # ' SIGMA(S) (NB) = ',E20.7/ # ' SIGMA(B) (NB) = ',E20.7/ # ' SIGMA(HAD) (NB) = ',E20.7/ # ' A_FB(MU) = ',E20.7/ # ' A_FB(C) = ',E20.7/ # ' A_FB(B) = ',E20.7///) 5 FORMAT(' E_CM (GEV) = ',E20.5// # ' ZZ QUANTITIES '/ # ' SIGMA(MU) (NB) = ',E20.7/ # ' SIGMA(U) (NB) = ',E20.7/ # ' SIGMA(D) (NB) = ',E20.7/ # ' SIGMA(C) (NB) = ',E20.7/ # ' SIGMA(S) (NB) = ',E20.7/ # ' SIGMA(B) (NB) = ',E20.7/ # ' SIGMA(HAD) (NB) = ',E20.7/ # ' A_FB(MU) = ',E20.7/ # ' A_FB(C) = ',E20.7/ # ' A_FB(B) = ',E20.7///) * 6 FORMAT(' E_CM (GEV) = ',E20.5// # ' NO ZG QUANTITIES '/ # ' SIGMA(MU) (NB) = ',E20.7/ # ' SIGMA(U) (NB) = ',E20.7/ # ' SIGMA(D) (NB) = ',E20.7/ # ' SIGMA(C) (NB) = ',E20.7/ # ' SIGMA(S) (NB) = ',E20.7/ # ' SIGMA(B) (NB) = ',E20.7/ # ' SIGMA(HAD) (NB) = ',E20.7/ # ' A_FB(MU) = ',E20.7/ # ' A_FB(C) = ',E20.7/ # ' A_FB(B) = ',E20.7///) * 7 FORMAT(' E_CM (GEV) = ',E20.5// # ' DSIGMA(MU) (PERCENT) = ',E20.7/ # ' DSIGMA(HAD) (PERCENT) = ',E20.7/ # ' DA_FB(MU) (ABS) = ',E20.7///) 8 FORMAT(' E_CM (GEV) = ',E20.5// # ' WISIGMA(MU) (PERCENT) = ',E20.7/ # ' WISIGMA(HAD) (PERCENT) = ',E20.7/ # ' WIA_FB(MU) (ABS) = ',E20.7///) STOP END