* Use samosatest.kumac to load in uniform phase space distribution into ntuple 1 and real s5 data into ntuple2 

h/cr/2d 101 uniformPaseSpace 50 -4 6 50 0 150
h/cr/2d 102 realS5Data       50 -4 6 50 0 150
nt/proj 101 1.sqrt(l)%log(1.+foldedsamosadist.f(maxnf,minnf,l,79.e3,9.9e3,29.5e3,336.3e3))
nt/proj 102 2.sqrt(l)%log(1.+foldedsamosadist.f(maxnf,minnf,l,79.e3,9.9e3,29.5e3,336.3e3))

h/op/div 102 101 103 ! ! e
h/plot 103 surf3