Skip to content

Commit

Permalink
adjust to uncertainties in literature values
Browse files Browse the repository at this point in the history
  • Loading branch information
LisaSchlueter committed Dec 11, 2024
1 parent 8133691 commit 96d918e
Showing 1 changed file with 1 addition and 2 deletions.
3 changes: 1 addition & 2 deletions src/simple_calibration.jl
Original file line number Diff line number Diff line change
Expand Up @@ -100,8 +100,7 @@ function get_peakhists_th228(e::Vector{<:Unitful.Energy{<:Real}}, th228_lines::V
# get histograms around calibration lines and peakstats
peakenergies = [ustrip.(e_unit, filter(in(peak - first(window)..peak + last(window)), e)) for (peak, window) in zip(th228_lines, window_sizes)]
peakstats = StructArray(estimate_single_peak_stats.(peakenergies, ustrip.(e_unit, bin_widths)))
peakhists = [fit(Histogram, ustrip.(e_unit, e), ustrip(e_unit, peak - first(window)):ps.bin_width:ustrip(e_unit, peak + last(window))) for (peak, window, ps) in zip(th228_lines, window_sizes, peakstats)]

peakhists = [fit(Histogram, ustrip.(e_unit, e), ustrip(e_unit, peak - first(window)):ps.bin_width:ustrip(e_unit, peak + last(window))) for (peak, window, ps) in zip(Measurements.value.(th228_lines), window_sizes, peakstats)]
peakhists, peakstats, h, peakstats.bin_width .* e_unit
end
export get_peakhists_th228
Expand Down

0 comments on commit 96d918e

Please sign in to comment.