unset multiplot
set term postscript enhanced color #size 400,600
set output "hist_out.ps"
#set term aqua size 400, 600
#set term x11 size 400, 600
reset
set key right inside samplen 1
Y_MIN=0
nturns = 100
set multiplot
set size 1.,0.23
set origin 0.,0.
n=100 #number of intervals
max=100. #max value
min=-100. #min value
width=(max-min)/n #interval width
#function used to map a value to the intervals
hist(x,width)=width*floor(x/width)+width/2.0
set xrange [-60:60]
#set yrange [0:]
unset yrange
#to put an empty boundary around the
#data inside an autoscaled graph.
#set offset graph 0.05,0.05,0.05,0.0
set xtics min,(max-min)/5,max
set boxwidth width*0.9
set style fill solid 0.5 #fillstyle
set tics out nomirror
#set xlabel "x[mm]"
set ylabel "Events"
#count and plot
plot filename u (hist(($4)*1.e3,width)):(1.0) smooth freq w boxes lc rgb"green" t "x[mm]"
clear
Y_MAX = (0+floor(GPVAL_Y_MAX/10)) * 10
#set ytics Y_MIN,(Y_MAX-Y_MIN)/4.,GPVAL_Y_MAX 
replot
set origin 0.,0.23
max=100. #max value
min=-100. #min value
width=(max-min)/n #interval width
set xrange [-60:60]
#set yrange [0:]
#to put an empty boundary around the
#data inside an autoscaled graph.
#set offset graph 0.05,0.05,0.05,0.0
set xtics min,(max-min)/5,max
set boxwidth width*0.9
set style fill solid 0.5 #fillstyle
set tics out nomirror
#set xlabel "y[mm]"
set ylabel "Events"
#count and plot
plot filename u (hist(($6)*1.e3,width)):(1.0) smooth freq w boxes lc rgb"green" t "y[mm]"
clear
Y_MAX = (0+floor(GPVAL_Y_MAX/10)) * 10
#set ytics Y_MIN,(Y_MAX-Y_MIN)/4.,GPVAL_Y_MAX 
replot
set origin 0.,0.46
max=3.  #max value
min=-3. #min value
width=(max-min)/n #interval width
set xrange [min:max]
#set yrange [0:]
#to put an empty boundary around the
#data inside an autoscaled graph.
#set offset graph 0.05,0.05,0.05,0.0
set xtics min,(max-min)/5,max
set boxwidth width*0.9
set style fill solid 0.5 #fillstyle
set tics out nomirror
#set xlabel "dp/p[%]"
set ylabel "Events"
#count and plot
plot filename u (hist(($9)*1.e2,width)):(1.0) smooth freq w boxes lc rgb"green" t "dp/p[%]"
clear
Y_MAX = (0+floor(GPVAL_Y_MAX/10)) * 10
#set ytics Y_MIN,(Y_MAX-Y_MIN)/4.,GPVAL_Y_MAX 
replot
set origin 0.,0.69
set size 1.,0.3
max=100. #max value
min=-100. #min value
width=(max-min)/n #interval width
set xrange [min:max]
set autoscale
#set yrange [0:]
#to put an empty boundary around the
#data inside an autoscaled graph.
#set offset graph 0.05,0.05,0.05,0.0
set xtics min,(max-min)/10,max
set xtics 20
set boxwidth width*0.9
set style fill solid 0.5 #fillstyle
set tics out nomirror
#set xlabel "t[ns]"
set ylabel "Events"
#count and plot
set title filename
plot filename u (hist(($10*1.e9),width)):(1.0) smooth freq w boxes lc rgb "green" t "t[ns]"
clear
xdif = (GPVAL_DATA_X_MAX-GPVAL_DATA_X_MIN)/5
set xtics int(GPVAL_DATA_X_MIN/10)*10, int(xdif/10)*10, int(GPVAL_DATA_X_MAX/10)*10 nomirror
Y_MAX = (0+floor(GPVAL_Y_MAX/10)) * 10
#set ytics Y_MIN,(Y_MAX-Y_MIN)/2.,GPVAL_Y_MAX 
incr = (GPVAL_Y_MAX-GPVAL_Y_MIN)/4.
#set ytics autofreq
#set ytics 0,int(incr)
replot
#