set term pdf fontscale 0.5 size 6in,4in
outfile = '2D_histogram_y.pdf'
set output outfile
#x-px phase space
#xmax=0.086;xmin=0.066
#ymax=0.008;ymin=-0.008

#y-y phase space
xmax=0.030;xmin=-0.030
ymax=0.006;ymin=-0.006

#x_e - Amplitude phase space
#ymax=0.4e-3;ymin=0
#xmin=-0.02;xmax=0.04

ix_bins=100;iy_bins=100
#xmax=0.045;xmin=-0.045;ix_bins=100
#ymax=0.006;ymin=-0.006;iy_bins=100
set key off
set view map
 set pm3d at s hidden3d 100
 set surface
 set dgrid3d 50,50,3
 set view map  scale 1
 set style increment default
# set style line 100 lt 5 lw 0.5
 unset hidden3d
 unset surf
# set pm3d interpolate 0,0
# set autoscale fix
# set contour base
# unset dgrid
set pm3d noborder
set palette defined ( 0 0.05 0.05 0.2, 0.1 0 0 1, 0.25 0.7 0.85 0.9,\
   0.4 0 0.75 0, 0.5 1 1 0, 0.7 1 0 0, 0.9 0.6 0.6 0.6,\
  1 0.95 0.95 0.95 )
    # set palette defined (0 0 0 0, 1 0 0 1, 3 0 1 0, 4 1 0 0, 6 1 1 1)
     set palette color

set palette defined ( 0 "white", 1 "black", 2 "blue" , 3"green")
s=1000
fx=(xmax-xmin)/ix_bins*s
x0=(xmax+xmin)/2*s
#x0=0
fy=(ymax-ymin)/iy_bins*s
y0=(ymax+ymin)/2*s
#y0=0
set xlabel 'y [mm]'
#set xlabel 'x [mm]'
set ylabel "y' [mrad]"
#set ylabel "Amplitude ((x-D_x{/Symbol d})^2/{/Symbol b}+x'^2{/Symbol b}) [mm-mrad]"
set xrange [xmin*s:xmax*s]

#set xtics xmin*s,2
set yrange [ymin*s:ymax*s]
#set zrange [0.00001:]
#set logscale z
#splot '2D_hist.dat' u (($1*fx+x0)):(($2*fy+y0)):($3) w pm3d
splot '2D_hist_y.dat' u (($1*fx+x0)):(($2*fy+y0)):3 w pm3d

print outfile