* Clear options, settings, and old data
set *
opt *
vec/del *
hist/del *

* Open the hardcopy output file
fort/file 66 klusdat.eps
metafile 66 -113

* Get the data
hfi 33 klusdat.hbk
hist/opera/add 12501 12502 3
hist/opera/add 12503 3 4
hist/opera/add 12505 4 5
hist/opera/add 12509 5 102
hist/opera/add 14501 14502 3
hist/opera/add 14503 3 4
hist/opera/add 14505 4 5
hist/opera/add 14509 5 103
close 33

* Set up the plot
opt nbox
set xsiz 20
set ysiz 10
set xmgl 2.5
set asiz 0.40
set gsiz 0.75
set ygti 0.50
set fwid 5
set *fon -60
set gfon -15

opt utit

igset mtyp 20

* Now add the mass plots
zone 2 1 1
set xlab 1.75
vec/create par(7) R 1200 1.23 0.10 0 0 0 0
igset txci 2 ; igset plci 2
hi/fit 102(1.16:1.32) G+P3 ! 7 par
igset txfp -15 ; igset chhe 0.75 ; igset txal 10 ; itx 1.3 10000 '[D]^++!'
igset txci 1 ; igset plci 1
atitle ' ' 'Events per 10 MeV'
igset txfp -15 ; igset chhe 0.75 ; igset txal 10 ; itx 1.6 5000 'p[p]^+!'
igset chhe 0.5
igset txfp -60 ; igset txal 10 ; itx 1.6 -3000 'Two-Particle Mass (GeV)'

vec/del par
vec/create par(7) R 541 0.77 0.10 14 -14 0 0
igset txci 2 ; igset plci 2
hi/fit 103(0.65:0.90) G+E+P1 ! 7 par
igset txfp -15 ; igset chhe 0.75 ; igset txal 10 ; itx 0.85 800 '[r]^0!'
igset txci 1 ; igset plci 1
igset txfp -15 ; igset chhe 0.75 ; igset txal 10 ; itx 1.0 500 '[p]^+![p]^-!'

* Clean up
close 66