Comparing GEANT to Data
q2 = 5, E = 644
post script file on monsoon  = /data/aniol/geant_cp_e89044/data_geant_5_644.ps
geant
data
peak location
epeak
elo,ehi:  in this example (ehi - elo) = 4 MeV


608.3
606.5  610.5
full width half maximum
FWHM


1.67 MeV
total counts in peak
pg_fitted= peak geant counts in fitted spectrum
xnorm = normalization factor for geant
peak shift = peak_data - peak_geant (in MeV)
pg = pg_fitted/xnorm
pdt = total peak data counts
fit region = 606.5 to 610.5 MeV
chi2pdf = 2.61

0.2377e6
3.02
0.725
78709






0.2346e6
background counts
background region ebglo to ebghi = 610.5 to 612. :  in this example (ebghi - ebglo)=1.5 MeV
data_bgd/MeV = data_bgd/(ebghi - ebglo)
background in data =bd= data_bgd/MeV*(ehi - elo)
geant_bgd/MeV/xnorm = geant_bgd/(ebghi - ebglo)/xnorm
background in geant = bg/xnorm = geant_bgd/MeV/xnorm*(ehi - elo)




222
888


757
3028
delta background = uncertainty in background
in data uncertainty=delbd : this is a judgment call
in geant = 0


0

2000
net peak = peak - background
in data, pd = pdt - bd
in geant pg = pg


77821

231572
error in data net peak = delpd
delpd =sqrt( pdt + delbd**2)


2060
geant total
epkt

151429

ratio = epkt/pg
1.946

uncertainty in ratio
delratio/ratio = sqrt(1/epkt + 1/pg)

0.004

data total = datatot = ratio*pd

4.488e5
(deldatatot/datatot)**2 = 1/epkt + 1/pg + (delpd/pd)**2 + (delratio/ratio)**2

9.3e-5
deldatatot/datatot = fractional uncertainty in the total data counts

1e-2