*read in data
nt/cre 1 'data' 10 ! ! nev nevsk diff chi1 chi2 chi mz len dedx enu
nt/read 1 sk3_data_chi.dat
*read in mc
nt/cre 2 'mc' 10 ! ! nev nevsk diff chi1 chi2 chi mz len dedx enu
nt/read 2 sk3_mcrock_chi.dat

*set graphical attributes for 2x2 plots
set ylab .8
set vsiz .3
set tsiz .3
set asiz .3
set xval .4

*book and fill histograms
1d 50 '' 30 0 20
1d 51 '' 30 0 20
nt/proj 50 1.diff
nt/proj 51 2.diff

opt liny
zone 2 2
set dmod 1                *sets line style to solid
norm 51 149               *normalize to the number of data events
h/pl 51
set mscf 1
set mtyp 1
h/pl 50 PSE               *plot data superimposed as points with error bars
graph/hplot/atit 'diff [\143]^2!"\114#10'

*********

1d 60 '' 30 0 800
1d 61 '' 30 0 800
nt/proj 60 1.chi1 $1.and.$2
nt/proj 61 2.chi1 $1.and.$2
norm 60 238
norm 61 238
h/pl 60 PE
h/pl 61 S
graph/hplot/atit '[\143]^2! diff"\114#10, [\143]^2!"\076#50'

1d 52 '' 30 0 20
1d 53 '' 30 0 20
nt/proj 52 1.diff
nt/proj 53 2.diff
norm 52 243
norm 53 243
h/pl 53
h/pl 52 PSE
graph/hplot/atit 'diff'

opt logy
1d 62 '' 30 0 800
1d 63 '' 30 0 800
nt/proj 62 1.chi1
nt/proj 63 2.chi1
norm 62 243
norm 63 243
h/pl 63
h/pl 62 PSE
graph/hplot/atit '[\143]^2!'