* KUMAC file to create dm2 fits from dm2.prt output from nuscan1.f

angle 30 -120
opt eah
opt date
opt fit
opt linx
opt liny
opt linz
opt htit
set gsiz 0.4
set asiz 0.2
set tsiz 0.3

tit '[d]m^2! study SK 1489 Days, 100yr os-mc, 0.2-10 GeV'


ve/cr idm(100)      i
ve/cr dm2(100)      r 
ve/cr aldm2(100)    r
ve/cr chisqel(100)  r
ve/cr chisqmu(100)  r
ve/cr chisqmup(100) r
ve/cr rlikel(100)   r
ve/cr rlikmu(100)   r
ve/cr rlikmup(100)  r
ve/cr alrlikel(100) r
ve/cr alrlikmu(100)  r
ve/cr alrlikmup(100) r
ve/cr Rosc(100)      r
ve/cr Roscer(100)    r

ve/read idm,dm2,aldm2,chisqel,chisqmu,chisqmup,rlikel,rlikmu,rlikmup,alrlikel,alrlikmu,alrlikmup,Rosc,Roscer dm2.prt (i10,13e11.2)

1dh 201 'Chisq vs L10(dm2/eV2) electr, mu norm'   100 -5.0 0.0 0.
1dh 202 'Chisq vs L10(dm2/eV2) muons, mu norm'    100 -5.0 0.0 0.
1dh 203 'Chisq vs L10(dm2/eV2) muons, el norm'    100 -5.0 0.0 0.

1dh 301 '-Ln(RLik) vs L10(dm2/eV2) el, mu norm'   100 -5.0 0.0 0.
1dh 302 '-Ln(RLik) vs L10(dm2/eV2) mu, mu norm'   100 -5.0 0.0 0.
1dh 303 '-Ln(RLik) vs L10(dm2/eV2) mu, el norm'   100 -5.0 0.0 0.

1dh 400 'R(osc) vs log_10(dm2/eV2)'   100 -5.0 0.0 0.

hi/put_vect/contents 201 chisqel
hi/put_vect/contents 202 chisqmu
hi/put_vect/contents 203 chisqmup

hi/put_vect/contents 301 alrlikel
hi/put_vect/contents 302 alrlikmu
hi/put_vect/contents 303 alrlikmup

hi/put_vect/contents 400 Rosc
hi/put_vect/errors  400 Roscer

hi/set/max 203 40.0
hi/set/min 203 0.0

zon 2 1
opt liny

hi/pl 203
hi/pl 202 s  
hi/pl 201 s 

opt liny

hi/set/max 303 40.0
hi/set/min 303 20.0

hi/pl 303
hi/pl 302 s
hi/pl 301 s

pict/print 'dm2.ps'
*pict/print 'dm2.gif'
*pict/print 'dm2.eps'

zon 1 1

hi/pl 400 E2
line -5 1 0 1

pict/prin 'rosc.ps'
*pict/prin 'rosc.eps'
*pict/prin 'rosc.gif'