* ===========
 MACRO UPDATE
* ===========
 EXEC 1
 WAIT
 EXEC 2
 WAIT
 EXEC 3
 WAIT
 EXEC 4
 WAIT
RETURN   
 
*================ 
 MACRO SET_VALUES
*================ 
   
   VEC/DEL MW
   VEC/CREATE MW(40) R
   VEC/INPUT MW(1:2)    80.360 0.370 | UA2
   
   VEC/INPUT MW(3:4)    79.910 0.390 | CDF RUN 0
   VEC/INPUT MW(5:6)    80.410 0.180 | CDF RUN 1A
   VEC/INPUT MW(7:8)    80.470 0.089 | CDF RUN 1B
   VEC/INPUT MW(9:10)   80.433 0.079 | CDF COMBINED
   
   VEC/INPUT MW(11:12)  80.350 0.270 | D0 1A
   VEC/INPUT MW(13:14)  80.498 0.095 | D0 1B
   VEC/INPUT MW(15:16)  80.482 0.091 | D0 COMBINED
   
   VEC/INPUT MW(17:18)  80.454 0.063 | TEVATRON
   
   VEC/INPUT MW(19:20)  80.452 0.062 | TEVATRON+UA2
   
   VEC/INPUT MW(21:22)  80.440 0.064 | ALEPH
   VEC/INPUT MW(23:24)  80.380 0.071 | DELPHI
   VEC/INPUT MW(25:26)  80.375 0.077 | L3
   VEC/INPUT MW(27:28)  80.485 0.065 | OPAL
   
   VEC/INPUT MW(29:30)  80.427 0.046 | LEP COMBINED
   
   VEC/INPUT MW(31:32)  80.250 0.110 | NUTEV
   
   VEC/INPUT MW(33:34)  80.436 0.037 | DIRECT:ALL
   VEC/INPUT MW(35:36)  80.417 0.035 | INDIRECT+DIRECT
   
   VEC/INPUT MW(37:38)  80.430 0.030 | CDF-RUN2
   VEC/INPUT MW(39:40)  80.376 0.034 | LEP-1 INDIRECT
   
   VEC/PRI MW
   
   RETURN
   
* ========
 MACRO 1
* ========
 EXEC SET_VALUES
 set xsiz 20.
 set ysiz 16.
 set lwid 2.
 exec oeps mw_worldav_4exp.ps
 max = 7			| ---->> this is # of mw points + 1
 igset chhe 0.5
 igset txfp -20  
 exec col 2  
 null 80.0 80.9 0. [max] 'AB' 
 AXIS 80.0 80.9 1.0 1.0 80.0 80.9 510 'U'
 igset txal 20
 yv = 0.5
 exec col 6
 itx 80.0 [yv] '80.0'
 itx 80.2 [yv] '80.2'
 itx 80.4 [yv] '80.4'
 itx 80.6 [yv] '80.6'
 itx 80.8 [yv] '80.8'
* 
 exec col 4
 igset txal 20 
 itx 80.4 0.1 'Mw (GeV)'  
 set dmod 1 ; set ltyp 1
 mk = 20
 
 exec band MW(33:33) MW(34:34) [max]  | draw the world average band in yellow
* 
*            MASS  ERR  MARKER COLOR POSITION  TEXT
*
 exec draw MW(9:9) MW(10:10)   [mk]       2      6.5 'CDF'
 exec draw MW(15:15) MW(16:16) [mk]       2      5.5 'D0'
 exec draw MW(19:19) MW(20:20) [mk]       4      4.5 'Hadron Colliders'
 
 exec draw  MW(29:29) MW(30:30) [mk]       4      4.0 'LEP II'
 exec draw  MW(33:33) MW(34:34) [mk]       2      3.0 'World Average'
 
 exec draw  MW(35:35) MW(36:36) [mk]       1      2.0 'LEP1/SLD/[n]N/Mt'
* 
 
 VEC/CREATE MWX(6) 
 VEC/COPY MW(9:9)   MWX(1:1) 
 VEC/COPY MW(15:15) MWX(2:2) 
 VEC/COPY MW(21:21) MWX(3:3) 
 VEC/COPY MW(23:23) MWX(4:4) 
 VEC/COPY MW(25:25) MWX(5:5) 
 VEC/COPY MW(27:27) MWX(6:6) 
* 
 VEC/CREATE EMW(6) 
 VEC/COPY MW(10:10)   EMW(1:1) 
 VEC/COPY MW(16:16) EMW(2:2) 
 VEC/COPY MW(22:22) EMW(3:3) 
 VEC/COPY MW(24:24) EMW(4:4) 
 VEC/COPY MW(26:26) EMW(5:5) 
 VEC/COPY MW(28:28) EMW(6:6) 
 
 applic sigma
   mws  = vsum(mwx)
   mwa  = mws/6
   yd   = (mwx-mwa)**2/emw**2
   chi2 = vsum(yd)
   ndf  = 6
 exit
 VEC/WRITE CHI2,ndf ! '' 'chi2/df = ',f4.1,'/',i2 ''
* 
* 
* write on some numbers ..
* 
 igset txal 13
 EXEC COL 4 ; ITX 80.0 4.5  $FORMAT($SIGMA(MW(19)),F5.3) '+/-' $FORMAT($SIGMA(MW(20)),F5.3)  | HADRON AVERAGE
 EXEC COL 4; ITX 80.0 4.0   $FORMAT($SIGMA(MW(29)),F5.3) '+/-' $FORMAT($SIGMA(MW(30)),F5.3)  | LEP-2
 EXEC COL 2; ITX 80.0 3.0   $FORMAT($SIGMA(MW(33)),F5.3) '+/-' $FORMAT($SIGMA(MW(34)),F5.3)  | WORLD AV DIRECT
 close 66
 RETURN
   
 
 
* ========
 MACRO 2
* ========
 EXEC SET_VALUES
 set xsiz 20.
 set ysiz 16.
 set lwid 2.
 exec ops mw_worldav_all_lep.ps
 max = 17			| ---->> this is # of mw points + 1
 igset chhe 0.4
 igset txfp -20  
 exec col 2  
 null 80.0 80.9 0. 17 'AB' 
 AXIS 80.0 80.9 1.0 1.0 80.0 80.9 510 'U'
 igset txal 20
 yv = 0.1
 exec col 6
 itx 80.0 [yv] '80.0'
 itx 80.2 [yv] '80.2'
 itx 80.4 [yv] '80.4'
 itx 80.6 [yv] '80.6'
 itx 80.8 [yv] '80.8'
* 
 exec col 4
 igset txal 20 
 itx 80.4 -0.8 'Mw (GeV)'  
 set dmod 1 ; set ltyp 1
 mk = 20
 
 exec band MW(33:33) MW(34:34) [max]  | draw the world average band in yellow

 
*            MASS  ERR  MARKER COLOR POSITION  TEXT
* exec draw  80.360 0.370 [mk]       2      16  'UA2 (W [\256] e[n])'
*
 exec draw MW(9:9) MW(10:10)   [mk]       2      15 'CDF'
 exec draw MW(15:15) MW(16:16) [mk]       2      14 'D0'
 exec draw MW(19:19) MW(20:20) [mk]       4      13 'Hadron Colliders'
 
 exec draw  MW(27:27) MW(28:28) [mk]       2      11  'OPAL'
 exec draw  MW(21:21) MW(22:22) [mk]       2      10  'ALEPH'
 exec draw  MW(25:25) MW(26:26) [mk]       2      9   'L3'
 exec draw  MW(23:23) MW(24:24) [mk]       2      8   'DELPHI'
 exec draw  MW(29:29) MW(30:30) [mk]       4      7   'LEP II'
 
 exec draw  MW(33:33) MW(34:34) [mk]       2      5   'World Average'
 
 exec col 1
 set ltyp 2 ; set dmod 2; line 80.0 4. 80.2 4.; line 80.6 4. 80.8 4.
 itx 80.27 4. 'INDIRECT Mw'
 set dmod 1; set ltyp 1
 igset chhe 0.4
 
 exec draw MW(31:31) MW(32:32) [mk]       1      3 'NUTEV'
 exec draw MW(39:39) MW(40:40) [mk]       1      2 'LEP1/SLD/[n]N/Mt'
* 
 VEC/CREATE MWX(6) 
 VEC/COPY MW(9:9)   MWX(1:1) 
 VEC/COPY MW(15:15) MWX(2:2) 
 VEC/COPY MW(21:21) MWX(3:3) 
 VEC/COPY MW(23:23) MWX(4:4) 
 VEC/COPY MW(25:25) MWX(5:5) 
 VEC/COPY MW(27:27) MWX(6:6) 
* 
 VEC/CREATE EMW(6) 
 VEC/COPY MW(10:10)   EMW(1:1) 
 VEC/COPY MW(16:16) EMW(2:2) 
 VEC/COPY MW(22:22) EMW(3:3) 
 VEC/COPY MW(24:24) EMW(4:4) 
 VEC/COPY MW(26:26) EMW(5:5) 
 VEC/COPY MW(28:28) EMW(6:6) 
 
 applic sigma
   mws  = vsum(mwx)
   mwa  = mws/6
   yd   = (mwx-mwa)**2/emw**2
   chi2 = vsum(yd)
   ndf  = 6
 exit
 VEC/WRITE CHI2,ndf ! '' 'chi2/df = ',f4.1,'/',i2 ''
* 
* 
* write on some numbers ..
* 
 igset txal 13
 exec col 4 ; itx 80.0 13.  $FORMAT($SIGMA(MW(19)),F5.3) '+/-' $FORMAT($SIGMA(MW(20)),F5.3)  | HADRON AVERAGE
 EXEC COL 4; ITX 80.0 7.    $FORMAT($SIGMA(MW(29)),F5.3) '+/-' $FORMAT($SIGMA(MW(30)),F5.3)  | LEP-2
 EXEC COL 2; ITX 80.0 5.    $FORMAT($SIGMA(MW(33)),F5.3) '+/-' $FORMAT($SIGMA(MW(34)),F5.3)  | world av direct
 exec col 1 ; itx 80.5 16.6 [h]^2!/Nexp = $UNQUOTE($FORMAT($sigma(chi2(1)),F3.1))/6  
 close 66
 RETURN
   
 
* ========
 MACRO 3
* ========
 EXEC SET_VALUES
 exec ops mw_worldav_all.ps
 max = 14 | ---->> this is # of mw points + 1
 exec setup [max] 0.1
 set dmod 1 ; set ltyp 1
 mk = 20
 
 exec band MW(35:35) MW(36:36) [max]  | WORLD AVERAGE BAND (INDIRECT+DIRECT)

 
*            MASS  ERR  MARKER COLOR POSITION     TEXT
 exec draw  MW(1:1) MW(2:2) [mk]      1      13      'UA2 (W [\256] e[n])' 0.2
*
 exec draw MW(3:3) MW(4:4) [mk] 2 12 'CDF(Run 0,  W [\256] e[n],[m][n])'
 exec draw MW(5:5) MW(6:6) [mk] 2 11 'CDF(Run 1A, W [\256] e[n],[m][n])'
 exec draw MW(7:7) MW(8:8) [mk] 2 10 'CDF(Run 1B, W [\256] e[n],[m][n])'
 exec draw MW(9:9) MW(10:10) [mk] 2 9  'CDF combined'
* 
 exec draw MW(11:11) MW(12:12) [mk] 1 8 'D0(Run 1A,  W [\256] e[n])'  0.05
 exec draw MW(13:13) MW(14:14) [mk] 1 7 'D0(Run 1B, W [\256] e[n])'   0.05
 exec draw MW(15:15) MW(16:16) [mk] 1 6 'D0 combined'                 0.05
 
 exec draw MW(19:19) MW(20:20) [mk] 4 5 'Hadron Collider Average'
* itx 80.85 4.6 '(25 MeV Common Error)'
 
 exec draw MW(31:31) MW(32:32) [mk] 1 4 'NUTEV  ([n]Fe)'
 exec draw MW(29:29) MW(30:30) [mk] 1 3 'LEP II (ee [\256] WW)'
 exec draw MW(35:35) MW(36:36) [mk] 1 2 'World Average'
* 
 VEC/CREATE MWX(5) 
 VEC/COPY MW(9:9)   MWX(1:1) 
 VEC/COPY MW(15:15) MWX(2:2) 
 VEC/COPY MW(1:1)   MWX(3:3) 
 VEC/COPY MW(29:29) MWX(4:4) 
 VEC/COPY MW(31:31) MWX(5:5) 
* 
 VEC/CREATE EMW(5) 
 VEC/COPY MW(10:10)   EMW(1:1) 
 VEC/COPY MW(16:16)   EMW(2:2) 
 VEC/COPY MW(2:2)     EMW(3:3) 
 VEC/COPY MW(30:30)   EMW(4:4) 
 VEC/COPY MW(32:32)   EMW(5:5) 
 
 applic sigma
   mws  = vsum(mwx)
   mwa  = mws/5
   yd   = (mwx-mwa)**2/emw**2
   chi2 = vsum(yd)
   ndf  = 5
 exit
 VEC/WRITE CHI2,ndf ! '' 'chi2/df = ',f4.1,'/',i2 ''

* 
* 
* write on some numbers ..
* 
 exec col 1
 itx 79.5 13. $FORMAT($SIGMA(MW(1)),F5.3) '+/-' $FORMAT($SIGMA(MW(2)),F5.3)    | UA2
 exec col 2
 itx 79.5 9.  $FORMAT($SIGMA(MW(9)),F5.3)  '+/-' $FORMAT($SIGMA(MW(10)),F5.3)  | CDF
 exec col 1
 itx 79.5 6.  $FORMAT($SIGMA(MW(15)),F5.3) '+/-' $FORMAT($SIGMA(MW(16)),F5.3)  | D0 average
 exec col 4
 itx 79.5 5.  $FORMAT($SIGMA(MW(19)),F5.3) '+/-' $FORMAT($SIGMA(MW(20)),F5.3)  | hadron average
 itx 79.5 4.  $FORMAT($SIGMA(MW(31)),F5.3) '+/-' $FORMAT($SIGMA(MW(32)),F5.3)  | Nutev
 itx 79.5 3.  $FORMAT($SIGMA(MW(29)),F5.3) '+/-' $FORMAT($SIGMA(MW(30)),F5.3)  | lep-2
 itx 79.5 2.  $FORMAT($SIGMA(MW(35)),F5.3) '+/-' $FORMAT($SIGMA(MW(36)),F5.3)  | world av
 exec col 2
 exec col 1; itx 80.8 14.0 [h]^2!/Nexp = $UNQUOTE($FORMAT($sigma(chi2(1)),F3.1))/5  
 close 66
 RETURN
 
* ========
 MACRO 4
* ========
 EXEC SET_VALUES
 exec ops mw_worldav_direct.ps
 max = 16 | ---->> this is # of mw points + 1
 exec setup [max] 0.1
 set dmod 1 ; set ltyp 1
 mk = 20
 
 exec band MW(33:33) MW(34:34) [max]  | DIRECT WORLD AVERAGE BAND
 
*            MASS  ERR  MARKER COLOR POSITION     TEXT
 exec draw  MW(1:1) MW(2:2) [mk]      1      15.5      'UA2 (W [\256] e[n])' 0.2
*
 exec draw MW(5:5) MW(6:6)   [mk] 2 14 'CDF(Run 1A, W [\256] e[n],[m][n])'
 exec draw MW(7:7) MW(8:8)   [mk] 2 13 'CDF(Run 1B, W [\256] e[n],[m][n])'
 exec draw MW(9:9) MW(10:10) [mk] 2 12 'CDF combined'
* 
 exec draw MW(11:11) MW(12:12) [mk] 4 10 'D0(Run 1A,  W [\256] e[n])'  0.05
 exec draw MW(13:13) MW(14:14) [mk] 4 9 'D0(Run 1B, W [\256] e[n])'   0.05
 exec draw MW(15:15) MW(16:16) [mk] 4 8 'D0 combined'                 0.05
 
 exec draw MW(19:19) MW(20:20) [mk] 1 6 'Hadron Collider Average'
* itx 80.85 5.6 '(25 MeV Common Error)'
 
 exec draw MW(29:29) MW(30:30) [mk] 2 4 'LEP II (ee [\256] WW)'
 exec draw MW(33:33) MW(34:34) [mk] 4 2 'World Average'
 
 VEC/CREATE MWX(4) 
 VEC/COPY MW(9:9)   MWX(1:1) 
 VEC/COPY MW(15:15) MWX(2:2) 
 VEC/COPY MW(1:1)   MWX(3:3) 
 VEC/COPY MW(29:29) MWX(4:4) 
* 
 VEC/CREATE EMW(4) 
 VEC/COPY MW(10:10)   EMW(1:1) 
 VEC/COPY MW(16:16)   EMW(2:2) 
 VEC/COPY MW(2:2)     EMW(3:3) 
 VEC/COPY MW(30:30)   EMW(4:4) 
 
 applic sigma
   mws  = vsum(mwx)
   mwa  = mws/4
   yd   = (mwx-mwa)**2/emw**2
   chi2 = vsum(yd)
   ndf  = 4
 exit
 VEC/WRITE CHI2,ndf ! '' 'chi2/df = ',f4.1,'/',i2 ''

* 
* write on some numbers ..
* 
 exec col 1
 itx 79.5 15.5 $FORMAT($SIGMA(MW(1)),F5.3) '+/-' $FORMAT($SIGMA(MW(2)),F5.3) | UA2
 exec col 2
 itx 79.5 14.  $FORMAT($SIGMA(MW(5)),F5.3) '+/-' $FORMAT($SIGMA(MW(6)),F5.3)  | cdf run 1a
 itx 79.5 13.  $FORMAT($SIGMA(MW(7)),F5.3) '+/-' $FORMAT($SIGMA(MW(8)),F5.3)  | cdf run 1b
 itx 79.5 12.  $FORMAT($SIGMA(MW(9)),F5.3) '+/-' $FORMAT($SIGMA(MW(10)),F5.3)  | cdf average
 exec col 4
 itx 79.5 10. $FORMAT($SIGMA(MW(11)),F5.3) '+/-' $FORMAT($SIGMA(MW(12)),F5.3) | D0 1a
 itx 79.5 9.  $FORMAT($SIGMA(MW(13)),F5.3) '+/-' $FORMAT($SIGMA(MW(14)),F5.3)   | D0 1b
 itx 79.5 8.  $FORMAT($SIGMA(MW(15)),F5.3) '+/-' $FORMAT($SIGMA(MW(16)),F5.3)   | D0 average
 exec col 1
 itx 79.5 6.  $FORMAT($SIGMA(MW(19)),F5.3) '+/-' $FORMAT($SIGMA(MW(20)),F5.3)   | hadron average
 exec col 2
 itx 79.5 4.  $FORMAT($SIGMA(MW(29)),F5.3) '+/-' $FORMAT($SIGMA(MW(30)),F5.3)   | lep-2
 exec col 4
 itx 79.5 2.  $FORMAT($SIGMA(MW(33)),F5.3) '+/-' $FORMAT($SIGMA(MW(34)),F5.3)   | world av
 exec col 2
 exec col 1 ; itx 80.8 16.5 [h]^2!/Nexp = $UNQUOTE($FORMAT($sigma(chi2(1)),F3.1))/4
 close 66
 RETURN 

* =========
 MACRO BAND 
* =========
   
* ---- world averge line ...
   
 exec col 1 ; set dmod 1; 
 VEC/CREATE XAV(1) R [1]  
 VEC/CREATE XER(1) R [2]
 MESSAGE MW = $SIGMA(XAV(1)) ERROR = $SIGMA(XER(1))
 
 VEC/CREATE Y(5) R 1. 1. 0. 0. 1.
 VEC/INPUT  Y(3:3) [3]
 VEC/INPUT  Y(4:4) [3]
 VEC/PRI Y
 
 APPLIC SIGMA
   XLOW=XAV-XER
   XHI=XAV+XER
 EXIT
 VEC/PRI XLOW  
 VEC/PRI XHI
 VEC/CREATE X(5)  
 VEC/COPY  XLOW(1:1)  X(1:1)
 VEC/COPY  XHI(1:1)   X(2:2)
 VEC/COPY  XHI(1:1)   X(3:3)
 VEC/COPY  XLOW(1:1)  X(4:4)
 VEC/COPY  XLOW(1:1)  X(5:5)
 
 SET FTYP 1
 SET FACI 5
 SET FAIS 1
 SET FASI 533
 SET BORD 1
 VEC/PRI X
 FAREA 5 X Y
 SET DMOD 2 ; SET LTYP 2; LINE $SIGMA(XAV(1)) 1. $SIGMA(XAV(1)) [3] ; SET DMOD 1 ; SET LTYP 1
 RETURN

* =========
 MACRO DRAW
* =========
   
   exec vcre x 1 r
   exec vcre y 1 r
   exec vcre ex 1 r
   exec vcre ey 1 r
   
   exec col [4]
   vec/input x [1]
   vec/input ex [2]
   vec/input ey 0.000001
   vec/input y [5]
   gra/hplo/err x y ex ey 1 [3] 0.4 '1'
   
   igset txal 13
* itx 80.80 [5] [6]
  itx $sigma(80.60+[7]) [5] [6]
 RETURN  
* ===========
 MACRO SETUP
* ===========
   
 igset txfp -20  
 exec col 2  
 null 79.5 81.5 0. [1] 'AB' 
 AXIS 79.5 81.5 1.0 1.0 79.4 81.4 510 'U'
 igset txal 20
 yv = [2]
 itx 79.5 [yv] '79.5'
 itx 79.7 [yv] '79.7'
 itx 79.9 [yv] '79.9'
 itx 80.1 [yv] '80.1'
 itx 80.3 [yv] '80.3'
 itx 80.5 [yv] '80.5'
 itx 80.7 [yv] '80.7'
 itx 80.9 [yv] '80.9'
 itx 81.1 [yv] '81.1'
 itx 81.3 [yv] '81.3'
 itx 81.5 [yv] '81.5'
 
 igset txal 20 
 itx 80.5 $sigma([yv]-0.7) 'Mw (GeV)'  
* igset txfp -130
return

*========
 MACRO COL
*========
   IGSET TXCI [1]
   SET HCOL [1]
   SET PMCI [1]
   SET PLCI [1]
 RETURN  

*=========
MACRO VCRE
*=========
  if $vexist([1])<>0 then
    vec/del [1]
  endif
  vec/create [1]([2]) [3]
return
*========
macro ops
*========
close 66 ; fort/file 66 [1] ; meta 66 -111
return
*========
macro oeps
*========
close 66 ; fort/file 66 [1] ; meta 66 -113
return