MACRO FINAL_SPRING97_XSECPLOT EPS=-113
* IF EPS=-111, PS FILE, WILL BE COLORED.
*IF EPS=-113, EPS FILE, WILL BE BW.
IF [EPS]=-113 THEN
f/file 56 TOP_XSEC_1997.eps
ELSEIF [EPS]=-111 THEN
f/file 56 TOP_XSEC_1997.ps
ENDIF
Meta 56 [eps]
set ygti 1.0
set txci 001
set xsiz 8.57
set ysiz 8.57
set xmgr 0.30
set xmgl 1.00
set ymgu 0.15
set ymgl 1.2
set ytic 0.2
set xval 0.15
set xtic 0.2
set yval 0.15
zone 1 1
set vsiz 0.22
set asiz 0.25
set ylab 1.2
igset cshi 0.02
igset chhe 0.3
igset txfp -112
igset lwid 1.0
opt liny
opt ndat
opt nbox
opt nsta
vector/delete *
dim=7
vector/create x([dim])
vector/create y([dim])
dim1=1
vector/create x1([dim1])
vector/create y1([dim1])
dim3=7
vector/create x3([dim3])
vector/create y3([dim3])
dim4=4
vector/create x4([dim4])
vector/create y4([dim4])
dim5=3
vector/create x5([dim5])
vector/create y5([dim5])
*
set ndvx 206
set ndvy 210
set bwid 3
set vfon -130
null 140 200 0 20
*title 'D0&/ Prelim Run 1 Top Cross Section'
*atitle 'Top Mass (GeV/c^2!)' 'Cross Section (pb)'
set txal 20
set chhe 0.25
set txfp -130
set tang 0
itx 170 -2.6 'Top Quark Mass (GeV/c^2!)'
set tang 90
itx 135 10 'Cross Section (pb)'
set tang 0
***************************
* D0 cross section values *
***************************
vector/input x3 140. 150. 160. 170. 180. 190. 200.
vector/input y3 10.52 8.25 6.80 5.74 5.01 4.67 4.36
vector/create dx3([dim3])
vector/create dy3([dim3])
vector/input dx3 0.01 0.01 0.01 0.01 0.01 0.01 0.01
vector/input dy3 4.17 3.05 2.34 1.88 1.57 1.43 1.31
*
*itx 166 17.0 '[s] = 7.08 "a# 2.23 at 160 GeV/c^2!'
*itx 166 16.0 '[s] = 5.98 "a# 1.82 at 170 GeV/c^2!'
*itx 166 15.0 '[s] = 5.22 "a# 1.55 at 180 GeV/c^2!'
*itx 166 15.0 '[s] = 5.76 "a# 1.75 at 172.7 "a# 8.5 GeV/c^2!'
*
vector/create xarea(14)
vector/create yarea(14)
vector/input xarea 140.05 150.  160. 170. 180. 190. 199.95 199.95 190. 180. 170. 160. 150. 140.05
vector/input yarea  14.69 11.30 9.14 7.62 6.58 6.10 5.66   3.05   3.24 3.44 3.87 4.46 5.19 6.35
* GTW -- Change fill style to color
set fais 1
IF [EPS]=-113 THEN
* GTW -- Change fill area to something grey
igset faci 7
ELSEIF [EPS]=-111 THEN
* GTW
igset faci 7
ENDIF
igset fasi -120
igset plci 001
igset txci 001
*farea 14 xarea yarea   |<--SK
* GTW -- Redraw the axis to get the tick marks back
igset plci 001
axis 140.0 140.0 0.0 20.0 0 20 210 U
* GTW -- Redraw the boundry box incase the farea has clobbered any of it
*        The box primitive didn't seem to work here, had to draw lines.
set lwid 3
line 140 0 200 0
line 200 0 200 20
line 200 20 140 20
line 140 20 140 0
*
igset txci 001
igset plci 001
igset faci 1
igset fais 0
igset fasi 1
igset plci 004
igset lwid 3.0
set dmod 1
*graph [dim3] x3 y3 c   |<--SK
igset lwid 1.0
set dmod 1
igset plci 001
*
IF [EPS]=-113 THEN
set txci 001
ELSEIF [EPS]=-111 THEN
set txci 001
ENDIF
*
set faci 1
set plci 1
set txci 001
***************************
* background subtr at 95% *
***************************
igset pmci 001
set txci 001
set lwid 3
* GTW -- Changed arrow to line; arrow
*        wasn't showing up for some *odd*
*        reason.  I hate paw.
*arrow 164.9 181.7  5.53  5.53 0.001
*arrow 172.8 173.8  7.20  7.20 0.001
*arrow 172.8 173.8  3.86  3.86 0.001
*arrow 173.3 173.3  3.86  7.20 0.001
*arrow 164.9 164.9  5.38  5.68 0.001
*arrow 181.7 181.7  5.38  5.68 0.001
************************************************
* D0 cross section aproved on Apr 4,97
*
*       sigt =   5.48 +/- 1.76 pb
* at      mt = 173.3  +/- 8.4  GeV/c^2
************************************************
igset plci 002
line 164.9 5.48  181.7  5.48
*line 172.8 7.24  173.8  7.24    |<--SK
*line 172.8 3.72  173.8  3.72    |<--SK
line 173.3 3.72  173.3  7.24
*line 164.9 5.33  164.9  5.63    |<--SK
*line 181.7 5.33  181.7  5.63    |<--SK
igset pmci 002
set txci 002
vector/input y1  5.48
vector/input x1 173.3
hplot/symbols x1 y1 1 20 0.2
*
vector/input x1 174
vector/input y1 17.1
hplot/symbols x1 y1 1 20 0.2
itx 182.0 16.8 'DO&/'
*
*hplot/symbols x1 y1 1 20 0.28
*key 175 17.0 20 'DO&/'
************************************************
* CDF cross section from FERMILAB-CONF-97/004-E
*   "Recent Results on Top Quark Physics at CDF"
*    by Michele Gallinaro
*
*       sigt =   7.5  +1.9 -1.6 pb
* at      mt = 176.8  +/- 6.5  GeV/c^2
************************************************
igset plci 007
line 170.3  7.5  183.3  7.5
line 176.8  5.9  176.8  9.4
igset pmci 007
set txci 007
vector/input y1  7.5
vector/input x1 176.8
hplot/symbols x1 y1 1 21 0.2
*
vector/input x1 174
vector/input y1 19
hplot/symbols x1 y1 1 21 0.2
itx 182.5 18.7 'CDF'
*
*hplot/symbols x1 y1 1 21 0.28
*key 175 19.0 21 'CDF'
************************************************
igset pmci 001
set txci 001
SET FAIS 1
* GTW -- Change to a good grey-scale value
SET FACI 7
SET PLCI 4
vector/create bt(4)
vector/create bs(4)
vector/input bt(1:4) 170   170   178   178
vector/input bs(1:4) 17.5  18.3  18.3  17.5
set lwid 3
* GTW -- Set to color fill area.
set fais 1
set fasi -120
*graph/primitive/farea 4 bt bs  |<--SK
* GTW -- moved line after the farea so that it shows up
*line 170 17.9 178 17.9     |<--SK
set faci 1
set plci 1
set fais 1
*goto done
igset chhe 0.36
igset txal 0
igset tang 0
igset chhe 0.25
*itx 180 17.5 'DO&/'    |<--SK
*
*itx 180 15.0 Laenen et al.
igset txci 003
itx 180 13.5 Laenen et al.
set dmod 1
*line 170 15.4 178 15.4
igset plci 003
line 170 13.9 178 13.9
*
*itx 180 13.5 Berger et al.
igset txci 004
itx 180 15.0 Berger et al.
set dmod 2
*line 170 13.9 178 13.9
igset plci 004
line 170 15.4 178 15.4
*
igset txci 006
itx 180 12.0 Catani et al.
set dmod 3
igset plci 006
line 170 12.4 178 12.4
*
*goto done
************************
* theoretical xsections
************************
IF [EPS]=-113 THEN
set txci 001
ELSEIF [EPS]=-111 THEN
set txci 001
ENDIF
*text 132 34 'Theory' 0.3
set txci 001
*
* laenen etal central cross section
*
* GTW- Lost the clipping; had to cut off 210 and 220
*vector/input x 140. 150. 160. 170. 180. 190. 200. 210. 220.
*vector/input y 16.9 11.7 8.16 5.83 4.21 3.06 2.26 1.68 1.25
vector/input x 140. 150. 160. 170. 180. 190. 200.
vector/input y 16.9 11.7 8.16 5.83 4.21 3.06 2.26
*igset plci 001
igset plci 003
igset lwid 3.0
set dmod 1
graph [dim] x y c
set dmod 1
igset lwid 1
igset plci 001
*
* berger etal central cross section
*
* GTW- Lost the clipping; had to cut off 210 and 220
vector/input x 140. 150. 160. 170. 180. 190. 200.
vector/input y 17.9 12.5 8.87 6.40 4.66 3.43 2.55
*vector/input x 140. 150. 160. 170. 180. 190. 200. 210. 220.
*vector/input y 17.9 12.5 8.87 6.40 4.66 3.43 2.55 1.91 1.43
*igset plci 001
igset plci 004
igset lwid 3.0
set dmod 2
graph [dim] x y c
set dmod 1
igset lwid 1
igset plci 001
*
*
* catani etal central cross section
*
* GTW- Lost the clipping; had to cut off 210 and 220
*vector/input x 140. 150. 160. 170. 180. 190. 200. 210. 220.
*vector/input y 15.8 11.1 7.78 5.62 4.05 2.98 2.20 1.64 1.23
vector/input x 140. 150. 160. 170. 180. 190. 200.
vector/input y 15.8 11.1 7.78 5.62 4.05 2.98 2.20
*igset plci 001
igset plci 006
igset lwid 3.0
set dmod 3
graph [dim] x y c
set dmod 1
igset lwid 1
igset plci 001
*
*****************************
OPT LINY
opt date
*wait
done:
close 56