* * -- see "MACRO SET_VALUES" for Mw, Mtop values Used. * * =========== MACRO UPDATE * =========== * EXEC OLD_NEW_MT ; WAIT EXEC WORLD_NUTEV ; WAIT EXEC WORLD_NUTEV ; WAIT EXEC WORLD_NUTEV_CCFR ; WAIT EXEC 40_5_PLOT ; WAIT EXEC WORLD ; WAIT EXEC CDF_WORLD ; WAIT EXEC TEV_WORLD ; WAIT EXEC CDF_D0_LEP2 ; WAIT EXEC CDF_D0_LEP2_LEP1 ; WAIT EXEC WORLD_POINT ; WAIT EXEC CDF_RUN2 ; WAIT EXEC WORLD_PRE_LHC ; WAIT EXEC TEVATRON ; WAIT EXEC TEVATRON_LEP2 ; WAIT EXEC TEVATRON_LEP2_LEP1 ; WAIT EXEC CDF_D0 ; WAIT EXEC CDF ; WAIT EXEC WORLD_MSSM ; WAIT RETURN *================ MACRO SET_VALUES *================ VEC/DEL MW VEC/CREATE MW(20) R VEC/INPUT MW(1:2) 80.470 0.089 | CDF 1B VEC/INPUT MW(3:4) 80.433 0.079 | CDF ALL VEC/INPUT MW(5:6) 80.482 0.091 | D0 ALL VEC/INPUT MW(7:8) 80.427 0.046 | LEP2 ALL VEC/INPUT MW(9:10) 80.454 0.063 | TEVATRON VEC/INPUT MW(11:12) 80.452 0.062 | TEVATRON+UA2 VEC/INPUT MW(13:14) 80.250 0.110 | NUTEV VEC/INPUT MW(15:16) 80.436 0.037 | DIRECT:ALL VEC/INPUT MW(17:18) 80.417 0.035 | INDIRECT+DIRECT VEC/INPUT MW(19:20) 80.430 0.030 | CDF-RUN2 VEC/PRI MW VEC/DEL MT VEC/CREATE MT(8) R VEC/INPUT MT(1:2) 176.1 6.6 | CDF VEC/INPUT MT(3:4) 172.1 7.1 | D0 VEC/INPUT MT(5:6) 174.3 5.1 | TEVATRON VEC/INPUT MT(7:8) 176.0 3.0 | CDF-RUN2 RETURN * ================ MACRO WORLD_NUTEV * ================ * EXEC OPS mwmt_world_direct_nutev.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC SET_VALUES EXEC DRAW_NUTEV MW(13:13) MW(14:14) 100. 1000. 130. 200. | nuTeV band EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 6 MW(15:15) MW(16:16) MT(5:5) MT(6:6) IGSET TXAL 31 | RIGHT ALIGN ITX 160. 80.460 'DIRECT M?W!' ARROW 170. 160. 80.4 80.4 ; ITX 159. 80.415 'pp^&&"\276#! + LEP2' * EXEC DRAW_POINT 21 2 MW(15:15) MW(16:16) MT(5:5) MT(6:6) | world-av:direct EXEC DRAW_LEP_SLC EXEC FONTS 2 SET FAIS 3 SET FACI 5 SET BORD 1. EXEC COL 1 IGSET TXAL 10 BOX 135. 140. 80.59 80.62 ITX 142. 80.595 'NuTeV : 100 [\074] M?H! [\074] 1000 GeV/c^2!' * EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ==================== MACRO WORLD_NUTEV_CCFR * ===================== * EXEC OPS mwmt_world_direct_nutev_ccfr.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.70 1 EXEC SET_VALUES EXEC DRAW_CCFR EXEC DRAW_NUTEV MW(13:13) MW(14:14) 100. 1000. 130. 200. | nuTeV band EXEC DRAW_HIGGS EXEC LABEL_HIGGS 27.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 6 MW(15:15) MW(16:16) MT(5:5) MT(6:6) IGSET TXAL 31 | RIGHT ALIGN ITX 160. 80.460 'DIRECT M?W!' ARROW 170. 160. 80.4 80.4 ; ITX 159. 80.415 'pp^&&"\276#! + LEP2' * EXEC DRAW_POINT 21 2 MW(15:15) MW(16:16) MT(5:5) MT(6:6) | world-av:direct EXEC DRAW_LEP_SLC EXEC FONTS 2 SET FAIS 3 SET FACI 5 SET BORD 1. EXEC COL 1 IGSET TXAL 10 BOX 135. 140. 80.59 80.62 ITX 142. 80.595 'NuTeV : 100 [\074] M?H! [\074] 1000 GeV/c^2!' * EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.70 2 CLOSE 66 RETURN * ================ MACRO CDF_D0_LEP2 * ================ * EXEC OPS mwmt_cdf_d0_lep2.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC SHADE_LEP $SIGMA(MW(7)) $SIGMA(MW(8)) | LEP-2 EXEC DRAW_ELLIPSE 3 MW(5:5) MW(6:6) MT(3:3) MT(4:4) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 170. 150. 80.58 80.58 ; ITX 149. 80.595 'D0' EXEC DRAW_ELLIPSE 4 MW(3:3) MW(4:4) MT(1:1) MT(2:2) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 175. 150. 80.5 80.5 ; ITX 149. 80.515 'CDF' EXEC 68 EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.655 2 CLOSE 66 RETURN * ===================== MACRO CDF_D0_LEP2_LEP1 * ===================== * EXEC OPS mwmt_cdf_d0_lep2_lep1.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC SHADE_LEP $SIGMA(MW(7)) $SIGMA(MW(8)) | LEP-2 EXEC DRAW_ELLIPSE 3 MW(5:5) MW(6:6) MT(3:3) MT(4:4) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 170. 150. 80.58 80.58 ; ITX 149. 80.595 'D0' EXEC DRAW_ELLIPSE 4 MW(3:3) MW(4:4) MT(1:1) MT(2:2) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 175. 150. 80.5 80.5 ; ITX 149. 80.515 'CDF' EXEC DRAW_LEP_SLC EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.655 2 CLOSE 66 RETURN * =============== MACRO 40_5_PLOT * ============== * EXEC OPS mwmt_40_5.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 6 80.41 0.040 MT(5:5) MT(6:6) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 159. 175. 80.41 80.41 ; ITX 158. 80.45 '[D]Mw = 40 MeV' ITX 158. 80.40 '[D]Mt = 5.1 GeV' * EXEC DRAW_LEP_SLC EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ========== MACRO WORLD * ========== * EXEC OPS mwmt_world_direct.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 6 MW(15:15) MW(16:16) MT(5:5) MT(6:6) IGSET TXAL 31 | RIGHT ALIGN ARROW 170. 160. 80.4 80.4 ; ITX 159. 80.415 'pp^&&"\276#! + LEP2' EXEC DRAW_LEP_SLC EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ============== MACRO CDF_WORLD * ============== * EXEC OPS mwmt_cdf_world_direct.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 4 MW(3:3) MW(4:4) MT(1:1) MT(2:2) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 170. 160. 80.52 80.52 ; ITX 159. 80.535 'CDF' EXEC DRAW_ELLIPSE 6 MW(15:15) MW(16:16) MT(5:5) MT(6:6) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 170. 160. 80.4 80.4 ; ITX 159. 80.415 'WORLD: Direct' EXEC DRAW_LEP_SLC EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ============== MACRO TEV_WORLD * ============== * EXEC OPS mwmt_tev_world_direct.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 4 MW(9:9) MW(10:10) MT(5:5) MT(6:6) | cdf+do combined IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 175. 160. 80.51 80.51 ; ITX 159. 80.525 'TEVATRON' EXEC DRAW_ELLIPSE 6 MW(15:15) MW(16:16) MT(5:5) MT(6:6) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 170. 160. 80.4 80.4 ; ITX 159. 80.415 'WORLD: Direct' EXEC DRAW_LEP_SLC EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ================ MACRO WORLD_POINT * ================ * EXEC OPS mwmt_world_direct_point.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_POINT 20 1 MW(15:15) MW(16:16) MT(5:5) MT(6:6) EXEC DRAW_LEP_SLC EXEC COL 1 ; IGSET TXAL 20 ; IGSET CHHE 0.7 ; ITX 174.3 80.50 'World Av.' EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ============= MACRO CDF_RUN2 * ============= * EXEC OPS mwmt_cdfrun2.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 4 MW(19:19) MW(20:20) MT(7:7) MT(8:8) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 176. 160. 80.43 80.43 ; ITX 159. 80.445 'CDF RUN-II' EXEC 68 EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ================== MACRO WORLD_PRE_LHC * ================== * EXEC OPS mwmt_world_prelhc.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 4 80.400 0.020 175.0 2.5 IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 175. 165. 80.4 80.5 ITX 187. 80.6 'Indirect Higgs Constraint' ITX 189. 80.55 'After LEP2 + Tevatron Run-II' EXEC 68 EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ============= MACRO TEVATRON * ============= * EXEC OPS mwmt_tevatron.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS 10.0 EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 4 MW(9:9) MW(10:10) MT(5:5) MT(6:6) | cdf+do combined IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 175. 160. 80.45 80.45 ; ITX 159. 80.465 'TEVATRON' EXEC 68 EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ================== MACRO TEVATRON_LEP2 * ================== * EXEC OPS mwmt_tevatron_lep2.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS 10.0 EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC SET_VALUES EXEC DRAW_ELLIPSE 4 MW(9:9) MW(10:10) MT(5:5) MT(6:6) | cdf+do combined IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 175. 160. 80.51 80.51 ; ITX 159. 80.525 'TEVATRON' EXEC 68 EXEC SHADE_LEP $SIGMA(MW(7)) $SIGMA(MW(8)) | LEP-2 EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ======================= MACRO TEVATRON_LEP2_LEP1 * ======================= * EXEC OPS mwmt_tevatron_lep2_lep1.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS 10.0 EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC SET_VALUES EXEC DRAW_ELLIPSE 4 MW(9:9) MW(10:10) MT(5:5) MT(6:6) | cdf+do combined IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 175. 160. 80.51 80.51 ; ITX 159. 80.525 'TEVATRON' EXEC 68 EXEC SHADE_LEP $SIGMA(MW(7)) $SIGMA(MW(8)) | LEP-2 EXEC DRAW_LEP_SLC EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * =========== MACRO CDF_D0 * =========== * EXEC OPS mwmt_cdf_d0.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS 10.0 EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 3 MW(5:5) MW(6:6) MT(3:3) MT(4:4) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 170. 150. 80.58 80.58 ; ITX 149. 80.595 'D0' EXEC DRAW_ELLIPSE 4 MW(3:3) MW(4:4) MT(1:1) MT(2:2) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 175. 150. 80.5 80.5 ; ITX 149. 80.515 'CDF' EXEC DRAW_LEP_SLC EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ======== MACRO CDF * ======== * EXEC OPS mwmt_cdf.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 4 MW(3:3) MW(4:4) MT(1:1) MT(2:2) IGSET TXAL 31 | RIGHT ALIGN SET LWID 1. ; ARROW 175. 150. 80.5 80.5 ; ITX 149. 80.515 'CDF' EXEC DRAW_LEP_SLC EXEC FONTS 1 EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * =============== MACRO WORLD_MSSM * =============== * EXEC OPS mwmt_world_mssm.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.8 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 24.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 6 MW(15:15) MW(16:16) MT(5:5) MT(6:6) IGSET TXAL 31 | RIGHT ALIGN ARROW 170. 150. 80.4 80.4 ; ITX 149. 80.415 'pp^&&"\276#! + LEP2' EXEC DRAW_LEP_SLC EXEC DRAW_MSSM IGSET CHHE 0.6 ; SET CSIZ 0.6 EXEC SETUP_AXIS 130. 200.0 80.001 80.8 2 CLOSE 66 RETURN * ========== MACRO FONTS * ========== IF [1] = 1 THEN SET GFON -20. FONT=-130 SET VFON [FONT] SET TFON [FONT] SET CFON [FONT] SET LFON [FONT] IGSET TXFP [FONT] IGSET CHHE 0.6 ENDIF IF [1] = 2 THEN SET GFON -20. FONT=-20 SET VFON [FONT] SET TFON [FONT] SET CFON [FONT] SET LFON [FONT] IGSET TXFP [FONT] IGSET CHHE 0.6 ENDIF RETURN * ================ MACRO SETUP_AXIS * =============== * EXEC FONTS 2 * IGSET TXAL 10 SET YSIZ 20.0 SET XMGL 5.0 SET YMGL 3.0 EXEC COL 1 SET XLAB 2.4 IF [5] = 1 THEN SET VSIZ 0.6 SET LWID 1. NULL [1] [2] [3] [4] SET XLAB 2.7 ; SET YLAB 1.5 SET ASIZ 0.6 ATITLE 'M?top! (GeV/c^2!)' 'M?W! (GeV/c^2!)' ELSE SET LWID 1. EXEC REDRAW-AXIS * LINE [2] [3] [2] [4] ENDIF RETURN * =============== MACRO REDRAW-AXIS * =============== X1 = $GRAFINFO('WNXMIN') X2 = $GRAFINFO('WNXMAX') Y1 = $GRAFINFO('WNYMIN') Y2 = $GRAFINFO('WNYMAX') SET FAIS 0 ; SET FACI 1 SET LWID 1.0 * BOX [X1] [X2] [Y1] [Y2] TICK XL [X1] ! TICK YB ! [Y1] RETURN * =============== MACRO DRAW_HIGGS * =============== v/create t(1) v/create w(1) v/create ew(1) v/create et(1) v/create mw(20) v/create mt(20) v/create mwu(20) v/create mwd(20) v/create dum(20) v/create mwall(40) v/create mtall(40) * set ftyp -3 * ----------- COLOR ------- COL = 7 * ------------------------ set faci [COL] set fais 3 set fasi 158 set bord 1. v/read mt,dum,dum,mw,mwu,mwd higgs_100_new_bk_mz911884_fj.data vect/copy mt(1:20) mtall(1:20) vect/copy mwu(1:20) mwall(1:20) call mwmtop_fill.for(mt,mwd) vect/copy xout(1:20) mtall(21:40) vect/copy yout(1:20) mwall(21:40) farea 40 mtall mwall v/read mt,dum,dum,mw,mwu,mwd higgs_250_new_bk_mz911884_fj.data vect/copy mt(1:20) mtall(1:20) vect/copy mwu(1:20) mwall(1:20) call mwmtop_fill.for(mt,mwd) vect/copy xout(1:20) mtall(21:40) vect/copy yout(1:20) mwall(21:40) farea 40 mtall mwall v/read mt,dum,dum,mw,mwu,mwd higgs_500_new_bk_mz911884_fj.data vect/copy mt(1:20) mtall(1:20) vect/copy mwu(1:20) mwall(1:20) call mwmtop_fill.for(mt,mwd) vect/copy xout(1:20) mtall(21:40) vect/copy yout(1:20) mwall(21:40) farea 40 mtall mwall v/read mt,dum,dum,mw,mwu,mwd higgs_1000_new_bk_mz911884_fj.data vect/copy mt(1:20) mtall(1:20) vect/copy mwu(1:20) mwall(1:20) call mwmtop_fill.for(mt,mwd) vect/copy xout(1:20) mtall(21:40) vect/copy yout(1:20) mwall(21:40) farea 40 mtall mwall set fais 1 set dmod 0 RETURN * ================ MACRO LABEL_HIGGS * ================ * * exec col 1 IGSET CHHE 0.6 IGSET TANG [1] IGSET CHHE 0.6 ITX 190 80.485 '100' ITX 190.0 80.42 '250' ITX 190.0 80.360 '500' ITX 190.0 80.305 '1000' IGSET CHHE 0.6 ; ITX 147.6 80.05 'Higgs Mass (GeV/c^2!)' IGSET TANG 0.0 * RETURN * =============== MACRO DRAW_POINT * =============== * Inputs : MTYP, COLOR, MW, MW_E, MT, MT_E SET MTYP [1] EXEC COL [2] INPUT W [3] INPUT EW [4] INPUT T [5] INPUT ET [6] MESSAGE MW = $SIGMA(W(1)) ERROR = $SIGMA(EW(1)) MESSAGE MT = $SIGMA(T(1)) ERROR = $SIGMA(ET(1)) SET HWID 5. SET LWID 5. * LINE $SIGMA([5]-[6]) [3] $SIGMA([5]+[6]) [3] * LINE [5] $SIGMA([3]-[4]) [5] $SIGMA([3]+[4]) GRAPH/HPLOT/ERRORS T W ET EW 4 [1] 0.50 '1' RETURN * ================= MACRO DRAW_ELLIPSE * ================= * * Inputs : COLOR, MW, MW_E, MT, MT_E * INPUT W [2] INPUT EW [3] INPUT T [4] INPUT ET [5] MESSAGE ----- ELLIPSE ------- MESSAGE MW = $SIGMA(W(1)) ERROR = $SIGMA(EW(1)) MESSAGE MT = $SIGMA(T(1)) ERROR = $SIGMA(ET(1)) * * DRAW ELLIPSE FOR 68% CL I.E SCALE ERRORS BY 1.515 * SET HWID 5. SET LWID 5. APPLIC SIGMA EET = ET*1.515 EEW = EW*1.515 X1 = ARRAY(100,T-EET#T+EET) DD = (X1-T)**2 DD = DD/EET**2 YY = 1-DD Y1 = W+EEW*SQRT(YY) Y2 = W-EEW*SQRT(YY) EXIT SET LWID 5. SET FTYP -3 EXEC COL [1] * ----------- COLOR ------- COL = [1] * ------------------------ SET FACI [COL] SET FAIS 3 SET FASI 244 SET BORD 1. FAREA 100 X1 Y1 FAREA 100 X1 Y2 PLINE 100 X1 Y1 PLINE 100 X1 Y2 RETURN * ================ MACRO LEP_CONTOUR * ================ EXEC OPS lep_contours.ps NULL 140.0 190.0 80.26 80.43 EXEC DRAW_LEP_SLC EXEC DRAW_LEP_SLC mw_mtop_sld_lep1_nutev_moriond99.dat CLOSE 66 RETURN * =============== MACRO OLD_NEW_MT * =============== EXEC OPS mwmt_world_direct_old_new_contour.ps EXEC SETUP_AXIS 130. 200.0 80.001 80.65 1 EXEC DRAW_HIGGS EXEC LABEL_HIGGS 29.0 EXEC SET_VALUES EXEC DRAW_ELLIPSE 6 MW(15:15) MW(16:16) MT(5:5) MT(6:6) IGSET TXAL 31 | RIGHT ALIGN ARROW 170. 160. 80.4 80.4 ; ITX 159. 80.415 'pp^&&"\276#! + LEP2' EXEC DRAW_LEP_SLC EXEC DRAW_LEP_SLC mw_mtop_sld_lep1_nutev_moriond99.dat EXEC SETUP_AXIS 130. 200.0 80.001 80.65 2 CLOSE 66 RETURN * ================ MACRO LEP_CONTOUR * ================ EXEC OPS lep_contours.ps NULL 140.0 190.0 80.26 80.43 EXEC DRAW_LEP_SLC EXEC DRAW_LEP_SLC mw_mtop_sld_lep1_nutev_moriond99.dat CLOSE 66 RETURN * ================= MACRO DRAW_LEP_SLC * ================= VEC/DEL WMASS VEC/DEL TMASS VEC/DEL INDEX IGSET CHHE 0.5 IGSET TXAL 10 SET LWID 5. SET FTYP -3 SET FAIS 3 SET FASI 244 SET BORD 1. SET BORD 1 IF [1] .NE. ' ' THEN COL = 4 EXEC COL [COL] SET FACI [COL] VEC/READ INDEX,TMASS,WMASS [1] VEC/CREATE XB(5) R 168. 173. 173. 168. 168. VEC/CREATE YB(5) R 80.06 80.06 80.08 80.08 80.06 PLINE 5 XB YB FAREA 5 XB YB ITX 174. 80.061 'LEP1,SLD,[n]N data' ELSE COL = 2 EXEC COL [COL] SET FACI [COL] VEC/READ INDEX,TMASS,WMASS mw_mtop_sld_lep1_nutev_lp99.dat VEC/CREATE XB(5) R 168. 173. 173. 168. 168. VEC/CREATE YB(5) R 80.10 80.10 80.12 80.12 80.10 PLINE 5 XB YB FAREA 5 XB YB ITX 174. 80.101 'LEP1,SLD,[n]N data' ENDIF GRAPH/PRIMI/GRAPH $VDIM(TMASS,1) TMASS WMASS 'CF' EXEC 68 RETURN *======== MACRO 68 * ======== IGSET CHHE 0.4 IGSET TXAL 31 EXEC COL 1 ; ITX 198. 80.037 'M?W!-M?top! contours : 68% CL' RETURN * ============== MACRO DRAW_MSSM * ============== exec col 4 vec/input wmass(1:7) 80.14 80.19 80.3 80.42 80.55 80.693 80.84 vec/input tmass(1:7) 130 140 160 180 200 220 240 graph/primi/graph 7 tmass(1:7) wmass(1:7) 'c' vec/input wmass(1:7) 80.24 80.3 80.45 80.63 80.84 81.05 81.28 vec/input tmass(1:7) 130 140 160 180 200 220 240 graph/primi/graph 7 tmass(1:7) wmass(1:7) 'c' set lwid 5. igset chhe 0.5 igset txal 10 * itx 135 81.1 'P. Chankowski, et al. Nucl. Phys. B417(1994) 101' * itx 135 81.05 'A. Dabelstein, et al. KA-THEP-5-1995 ' itx 132 80.76 '90 [\074] M?Higgs! [\074] 1000 GeV' itx 132 80.715 ' no SUSY particles at LEP II' igset chhe 0.7 itx 185 80.60 'MSSM' set lwid 1. arrow 171. 183.3 80.733 80.628 RETURN * ========= MACRO OPS * ========= SET * IGSET * OPT NGRI OPT NBOX * SET XSIZ 25. SET YSIZ 28. SET GSIZ 0.8 SET YGTI 0.5 SET YLAB 1.3 SET YMGU 2. SET YMGL 2. SET XMGL 3.0 SET XMGR 3.0 SET XVAL 0.5 SET GFON -20. FONT=-130 SET VFON [FONT] SET TFON [FONT] SET CFON [FONT] SET LFON [FONT] IGSET TXFP [FONT] IGSET CHHE 0.6 * SET CHHE 0.5 SET CSIZ 0.5 SET VSIZ 0.5 SET ASIZ 0.5 SET YHTI 100. * SET YLAB 1.0 SET XLAB 1.8 SET XWIN 3.0 * SET YHTI 1. SET TSIZ 0.5 close 66 ; fort/file 66 [1] ; meta 66 -111 RETURN * ======== MACRO COL * ======== IGSET TXCI [1] SET HCOL [1] SET PMCI [1] SET PLCI [1] RETURN * ============== MACRO SHADE_LEP * ============== COL = 6 EXEC COL [COL] ; LINE 130. [1] 200. [1] SET FTYP -3 SET FACI [COL] SET FAIS 3 SET FASI 533 SET BORD 1. SET BORD 1 VEC/CREATE X(5) R 130. 200. 200. 130. 130. VEC/CREATE Y(5) R $SIGMA([1]-[2]) $SIGMA([1]-[2]) $SIGMA([1]+[2]) $SIGMA([1]+[2]) $SIGMA([1]-[2]) FAREA 5 X Y IGSET CHHE 0.5 IGSET TXAL 10 EXEC FONTS 2 ; EXEC COL 1 ; ITX 135. $SIGMA([1]+0.01) 'LEP2' RETURN * =============== MACRO DRAW_NUTEV * =============== * * NUTEV MW DEPENDS ON MT AND MHIGGS : DRAW CONTOUR. * * INPUTS ARE : [1] = MW, [2] = MW-ERROR , [3],[4]:MH RANGE, [5],[6]:MT RANGE * VEC/CREATE XP(5) R [5] [6] [6] [5] [5] VEC/CREATE YP(5) R * APPLIC COMIS QUIT SUBROUTINE GRANGE(MW,MWE,MT1,MT2) VECTOR YP REAL MW,MWE,MT1,MT2,MT,MW1,MW2 MT = 174.3 MH1 = 100.0 MH2 = 1000.0 WRITE(*,90) MW,MWE,MT1,MT2,IFLAG 90 FORMAT(4(1X,F8.4),1X,I2) MW1 = MW-MWE MW2 = MW+MWE YP(1) = MW1 + 0.073*(MT1**2-MT**2)/(100**2) - 0.025*LOG(MH2/150.0) YP(2) = MW1 + 0.073*(MT2**2-MT**2)/(100**2) - 0.025*LOG(MH2/150.0) YP(3) = MW2 + 0.073*(MT2**2-MT**2)/(100**2) - 0.025*LOG(MH1/150.0) YP(4) = MW2 + 0.073*(MT1**2-MT**2)/(100**2) - 0.025*LOG(MH1/150.0) YP(5) = YP(1) END QUIT CALL GRANGE([1],[2],[5],[6]) MESSAGE 'NUTEV POINTS --' VEC/PRI XP VEC/PRI YP SET FAIS 3 SET FACI 5 SET FASI 244 EXEC COL 5 FAREA 5 XP YP RETURN * =============== MACRO DRAW_CCFR * =============== * * -- Just draw as shaded area -- region taken from Druckers talk at SLAC-98. * VEC/DEL XP; VEC/DEL YP MT_MIN = 130. MT_MAX = 200. VEC/CREATE XP(5) R [MT_MIN] [MT_MAX] [MT_MAX] [MT_MIN] [MT_MIN] VEC/CREATE YP(5) R * APPLIC COMIS QUIT SUBROUTINE CCFR VECTOR YP,XP WRITE(*,*) 'MT1,MT2 = ',MT1,MT2 YP(1) = (-0.150/110.0)*XP(1) + 80.377 YP(2) = (-0.150/110.0)*XP(2) + 80.377 YP(3) = (-0.081/110.0)*XP(3) + 80.670 YP(4) = (-0.081/110.0)*XP(4) + 80.670 YP(5) = YP(1) END QUIT CALL CCFR MESSAGE 'CCFR POINTS --' VEC/PRI XP VEC/PRI YP SET FAIS 3 SET FACI 3 SET FASI 244 FAREA 5 XP YP SET BORD 1. EXEC COL 1 IGSET TXAL 10 SET LWID 5. BOX 135. 140. 80.64 80.67 ITX 142. 80.645 'CCFR' RETURN