Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Generic simple_calibration #108

Draft
wants to merge 16 commits into
base: main
Choose a base branch
from
+231 −95
Draft

Generic simple_calibration #108

Changes from 1 commit
Commits
Show all changes
16 commits
Select commit Hold shift + click to select a range
a0b628f
add more generic verison of simple calibration, e.g. for co60 data
LisaSchlueter Dec 12, 2024
fa21052
take most prominent peak not not largest energy
LisaSchlueter Dec 12, 2024
1639f38
add minimal peakshape without tail
LisaSchlueter Dec 14, 2024
888996b
update docstring and add fall-back solution for fwhm calculation
LisaSchlueter Dec 16, 2024
421e2c3
add e_unit as kwarg
LisaSchlueter Dec 16, 2024
692c1de
merge commong function for th228 and generic gamma
LisaSchlueter Jan 10, 2025
b8dda3f
clean up th228 remainders and docstrings
LisaSchlueter Jan 10, 2025
8e924d9
merge simple calibration plto recipes
LisaSchlueter Jan 10, 2025
13a85f9
add simple calibration test
LisaSchlueter Jan 11, 2025
4851690
bug fix tests
LisaSchlueter Jan 11, 2025
da82248
account for case when dof = 0. uncertainties can still be propagated
LisaSchlueter Mar 7, 2025
a84f399
add fhagemann suggestions
LisaSchlueter Mar 7, 2025
8ff7002
Merge branch 'main' into co60
LisaSchlueter Mar 7, 2025
db57f33
bug fix
LisaSchlueter Mar 11, 2025
35f9b55
Merge branch 'co60' of github.com:legend-exp/LegendSpecFits.jl into co60
LisaSchlueter Mar 11, 2025
cace095
remove broadcasting for now; causes problems
LisaSchlueter Mar 11, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
8 changes: 7 additions & 1 deletion src/chi2fit.jl
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -94,7 +94,13 @@ function chi2fit(f_fit::Function, x::AbstractVector{<:Union{Real,Measurement{<:R
# gof
chi2min = res.objective
dof = length(x) - length(v_chi2)
pvalue = ccdf(Chisq(dof), chi2min)
if dof == 0
pvalue = NaN
@warn "The number of fit parameters is equivalent to number of data points --> dof = 0 ; p-value not meaningful, set to NaN"
else
pvalue = ccdf(Chisq(dof), chi2min)
end

function get_y_pred_err(f_fit, x_val, x_err, pars) # get final uncertainties for normalized residuals
dual_x = ForwardDiff.Dual{UncertTag}(x_val, x_err)
dual_y = f_fit(dual_x, pars...)
Expand Down