macro plot_pi0_photons
* 
* Plot histograms obtained using pi_plot from the bggen.nt file
* 02/09/11 ES
*
*
* set options
*
* hi/create/title_global 'BGGEN Photons from [p]^0! /10 MeV/s E"g#0.1,0.06GeV'
hi/create/title_global 'BGGEN Photons from [p]^0! /10 MeV/s E"g#0.2GeV'
* hi/create/title_global 'BGGEN Photons from [c] /10 MeV/s E"g#0.5GeV'
* hi/create/title_global 'BGGEN Photons from [c] /10 MeV/s E"g#0.1,0.06GeV'
option ndate
option nbox
* set stat 1111111
set stat 1111
option stat
set fit 111
option fit
option ngrid
*
* plotting options
*
* set *
set xmgl 4.
set ymgl 4.
set asiz 0.4
* set xlab 2.
set xlab 1.25
set ylab 1.
set xsiz 20.
set xmgl 3.
set ymgl 3.
set ysiz 20.
set gsiz 0.4
*
* set font definitions to bold roman 
*
set CFON -21
set GFON -21
set LFON -21
set TFON -21
set VFON -21
set txfp -21
set SMGU 0.02
set SMGR 0.02
set CSIZ 0.33
set VSIZ 0.3
set TSIZ 0.35
set YHTI 0.9
set HWID 3.0
set BWID 3.0
*
zone
*
* get ntuple
*
h/file 2 bggen.nt 0 -x
*
* define histograms per 10 MeV
*
* momentum
*
* 1dh 100 'fcal only photons' 400 0 4
* 1dh 101 'bcal only photons' 400 0 4
* 1dh 102 'fcal and bcal photons' 400 0 4
* 1dh 103 'failed photons' 400 0 4
*
* angle
*
1dh 100 'fcal only photons' 180 0 180
1dh 101 'bcal only photons' 180 0 180
1dh 102 'fcal and bcal photons' 180 0 180
1dh 103 'failed photons' 180 0 180
*
npts1 = 1001
*
* message PDE ntubes dt DR sqrtDR pedr atten angle90 angle45 angle15 
* message [PDE] [ntubes] [dt] [DR] [sqrtDR] [pedr] [atten] [angle90] [angle45] [angle15] 
* wait
*
emin=0
emax=4
ymin=0
ymax=10000
*  KPHYTH flags gamma=22, pi0=111, eta=221
nevents=999999999
*
* nt/plot //lun2/9.pi_plot.f(0,0,111,100) pin(3,1).gt.2 [nevents]
* nt/plot //lun2/9.pi_plot.f(1,0,221,100)  !  [nevents]
nt/plot //lun2/9.pi_plot.f(1,0,111,100)  !  [nevents]
*
zone 2 2
*
* open metafile
*
* for/file 66 plot_pi0_photons.ps
* meta 66 -111
for/file 66 plot_pi0_photons_angle_cut2.eps
meta 66 -113
*
*
option stat
option logy
csize = 0.05
igset chhe [csize]
* hplot/null  [emin] [emax] [ymin] [ymax]
* hplot/atitle   'Fcal only p (GeV)'  'Photons/10 MeV/s' 
* hi/plot 100 's'
hi/plot 100
* igset chhe 0.3
* exe window#push
* itx 0 1.1 'Angle "G#  10 deg'
* exe window#pop
csize = 0.05
igset chhe [csize]
* hplot/null  [emin] [emax] [ymin] [ymax]
* hplot/atitle   'Bcal only p (GeV)'  'Photons/10 MeV/s' 
* hi/plot 101 's'
hi/plot 101 
csize = 0.05
igset chhe [csize]
* hplot/null  [emin] [emax] [ymin] [ymax]
* hplot/atitle  'Fcal and Bcal p (GeV)'  'Photons/10 MeV/s' 
* hi/plot 102 's'
hi/plot 102 
csize = 0.05
igset chhe [csize]
* hplot/null  [emin] [emax] [ymin] [ymax]
* hplot/atitle 'Failed p (GeV)'  'Photons/10 MeV/s' 
* hi/plot 103 's'
hi/plot 103 

*
*
*
 close 66
*
exitm 
return