DISPLAY '<<< executing script file = '//STRING(?0)
@extrema.init

infile = 'RANGE_3D.txt'
infile = STRING(?1)

read infile\125 ii xx yy zz

xr = xx/1.e7   ! convert angstroms to mm

generate xbin 0.005 0.010 2.499

bin xr xbin xrb
exrb = sqrt(xrb)
statistics\-mess xrb xrbsum\sum xrbmin\min xrbmax\max
xr1 = 1.
xr2 = 1.125
xr3 = 1.250
xr4 = 1.500
yrl = xrbmin
yrh = xrbmax

statistics\-mess xrb[where(xbin<xr1)] xrgsum[1]\sum
statistics\-mess xrb[where(xbin>=xr1 & xbin<xr2)] xrgsum[2]\sum
statistics\-mess xrb[where(xbin>=xr2 & xbin<xr3)] xrgsum[3]\sum
statistics\-mess xrb[where(xbin>=xr3 & xbin<xr4)] xrgsum[4]\sum
statistics\-mess xrb[where(xbin>=xr4)] xrgsum[5]\sum

scale 0 2.5 5 0 10 2
set autoscale y

set xlabel 'Range (mm)'
set ylabel 'Events'

clear

set histogramtype 1
graph xbin xrb exrb

graph\overlay [xr1;xr1] [yrl;yrh]
graph\overlay [xr2;xr2] [yrl;yrh]
graph\overlay [xr3;xr3] [yrl;yrh]
graph\overlay [xr4;xr4] [yrl;yrh]

set textinteractive 0

set %xtextlocation 60
set %textheight 2
dytext = 4

ytext = 85
set %ytextlocation ytext
text infile

ytext = ytext - dytext
set %ytextlocation ytext
text 'R1: sum='//rchar(xrgsum[1])//', fraction='//rchar(xrgsum[1]/xrbsum)

ytext = ytext - dytext
set %ytextlocation ytext
text 'S1: sum='//rchar(xrgsum[2])//', fraction='//rchar(xrgsum[2]/xrbsum)

ytext = ytext - dytext
set %ytextlocation ytext
text 'S2: sum='//rchar(xrgsum[3])//', fraction='//rchar(xrgsum[3]/xrbsum)

ytext = ytext - dytext
set %ytextlocation ytext
text 'S3: sum='//rchar(xrgsum[4])//', fraction='//rchar(xrgsum[4]/xrbsum)

ytext = ytext - dytext
set %ytextlocation ytext
text 'R2: sum='//rchar(xrgsum[5])//', fraction='//rchar(xrgsum[5]/xrbsum)

display '>>>'