**test against Belle's inclusive BB->pi0 spectrum...

** use the "generic" test - 
** [1] is the histogram id of the inclusive pi0 spectrum (=13 currently)
** [2] is the number of Upsilon decays.
** [3] =1 for David's nice kumac - need stuff in David's envirn though...
hi/fil 1 generic.hst

vec/create mom(30) R 0.05 0.15 0.25 0.35 0.45 0.55 0.65 0.75 0.85 0.95 1.05 1.15 1.25 1.35 1.45 1.55 1.65 1.75 1.85 1.95 2.05 2.15 2.25 2.35 2.45 2.55 2.65 2.75 2.85 2.95
vec/create xbinw(30) R 0.05 0.05 0.05 0.05 0.05  0.05 0.05 0.05 0.05 0.05  0.05 0.05 0.05 0.05 0.05  0.05 0.05 0.05 0.05 0.05  0.05 0.05 0.05 0.05 0.05  0.05 0.05 0.05 0.05 0.05 

vec/create xsect(30) R 4.656 9.936 8.760 7.015 5.640 3.861 2.802 2.087 1.421 1.079 0.766 0.528 0.425 0.249 0.228 0.170 0.147 0.089 0.077 0.061 0.062 0.041 0.048 -0.025 -0.016 -0.001 -0.013 0.010 -0.005 -0.002

vec/create xsecte(30) R 0.272 0.621 0.372 0.247 0.171 0.135 0.101 0.080 0.068 0.054 0.038 0.033 0.029 0.026 0.023 0.020 0.017 0.016 0.014 0.012 0.011 0.012 0.011 0.012 0.10 0.010 0.010 0.009 0.008 0.009

**add also the systematic error
syserr=0.064
syserr=[syserr]*[syserr]

vdiv xsecte xsect xsecte
vmult xsecte xsecte xsecte
vbias xsecte [syserr] xsecte
SIGMA xsecte=sqrt(xsecte)

vmult xsecte xsect xsecte

hiscal=1.0/[2]
hi/copy [1] 1000
hi/oper/add 1000 1000 1000 0 [hiscal]

if ( [3]=1) then
   opt nsta
   opt ndat
   exec beautify.kumac
   set hwid 7
   title ' '
endif
  
max 1000 1.2
set htyp -3
hi/plot 1000
set htyp
vscal xsect 0.1 xsect
vscal xsecte 0.1 xsecte

if ( [3]=1) then
  blue
endif

hplot/err mom xsect xbinw xsecte 30 21

if ( [3]=1) then
  black
endif

atitle '[p]^0! momentum' 'Number of [p]^0! per [Y](4s) per 0.1 GeV/c' ! 220

if ( [3]=1 ) then
   igset chhe 0.6
   igset txal 
   itx 1.0 0.62 'EvtGen'
   blue
   itx 1.0 0.7 'Belle-Conf-0008 (Osaka)'
   black
   pic/print pi0.gif
   exec eps_print pi0.eps
endif