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 |