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.

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

systematic uncertainties from fix parameters. #111

Merged
merged 4 commits into from
Sep 13, 2023
Merged

Conversation

jiangyi15
Copy link
Owner

@jiangyi15 jiangyi15 commented Sep 12, 2023

Propagate uncertainties from fix params to fit params using AD.
The minimal condition

$$- \frac{\partial \ln L(a,b)}{\partial a} = 0$$

is a implicit function $a(b)$. The jacobian is

$$\frac{\partial a}{\partial b} = - (\frac{\partial^2 \ln L}{\partial a \partial a})^{-1} \frac{\partial^2 \ln L}{\partial a \partial b}.$$

The uncertainties from fix paramters $b$, can be propagete to $a$ as

$$V_a = \frac{\partial a}{\partial b} V_b \frac{\partial a}{\partial b}.$$

V = config.attach_fix_params_error({"A_mass": 0.01})

Then V
will be added to config.inv_he.
The config.params_trans() will include the uncertainties of fix parameters.
To solve the systematic uncertainties only, set config.inv_he=None first.

@codecov
Copy link

codecov bot commented Sep 12, 2023

Codecov Report

Patch coverage: 88.09% and project coverage change: +0.02% 🎉

Comparison is base (c2ea2d5) 76.18% compared to head (11cde8a) 76.21%.

Additional details and impacted files
@@            Coverage Diff             @@
##              dev     #111      +/-   ##
==========================================
+ Coverage   76.18%   76.21%   +0.02%     
==========================================
  Files         108      108              
  Lines       15855    15894      +39     
  Branches     2971     2981      +10     
==========================================
+ Hits        12079    12113      +34     
- Misses       3081     3084       +3     
- Partials      695      697       +2     
Flag Coverage Δ
unittests 76.21% <88.09%> (+0.02%) ⬆️

Flags with carried forward coverage won't be shown. Click here to find out more.

Files Changed Coverage Δ
tf_pwa/amp/core.py 76.41% <80.00%> (+<0.01%) ⬆️
tf_pwa/config_loader/config_loader.py 68.37% <88.88%> (+0.66%) ⬆️
tf_pwa/tests/test_amp.py 100.00% <100.00%> (ø)
tf_pwa/tests/test_full.py 100.00% <100.00%> (ø)

☔ View full report in Codecov by Sentry.
📢 Have feedback on the report? Share it here.

@jiangyi15 jiangyi15 merged commit 20ac860 into dev Sep 13, 2023
12 checks passed
@jiangyi15 jiangyi15 deleted the syst_fix_params branch September 13, 2023 14:48
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

1 participant