Skip to content

ccbar asymmetry#1629

Merged
tgiani merged 14 commits into
masterfrom
ccbar
Dec 1, 2022
Merged

ccbar asymmetry#1629
tgiani merged 14 commits into
masterfrom
ccbar

Conversation

@tgiani
Copy link
Copy Markdown
Contributor

@tgiani tgiani commented Nov 11, 2022

addressing #1627

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 11, 2022

when adding the new basis in vp2 I ve reformatted the file using black, not sure if that was a good idea. The only actual modification I did to pdfbasis is at the end, introducing the basis ccbar_asymm

@scarlehoff scarlehoff self-requested a review November 14, 2022 09:08
Comment thread validphys2/src/validphys/pdfbases.py Outdated
"""
if fitbasis == 'NN31IC':
if fitbasis == "NN31IC":
return np.identity(8)
Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

If the output is to be 9 elements now it has to be 9 elements everywhere.

I'm guessing this will require changes anywhere the 8 is set by hand (which is potentially more places than just here)

@scarlehoff
Copy link
Copy Markdown
Member

scarlehoff commented Nov 14, 2022

While I'm usually in favour of passing all files through black, the matrix-like shape of the transformation at the end of the file would be helpful to keep... black creates a bit of a mess there (and just because of the size of the line).

I'm happy with black but make sure it doesn't modify those lines (might be enough with -l 110 or whatever)

edit: ups, clicked close instead of comment :_

@scarlehoff scarlehoff closed this Nov 14, 2022
@scarlehoff scarlehoff reopened this Nov 14, 2022
@Zaharid
Copy link
Copy Markdown
Contributor

Zaharid commented Nov 14, 2022

@tgiani indeed please don't reformat lines that haven't otherwise being touched. It makes reviewing this a pain and messes up the history. Not to mention that the formatting gets worse in a number of places.

@scarlehoff
Copy link
Copy Markdown
Member

Have you tried doing a first fit with this basis for Theory 200 ?

@scarlehoff scarlehoff added the run-fit-bot Starts fit bot from a PR. label Nov 16, 2022
@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 16, 2022

@scarlehoff I m doing it now

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 16, 2022

should be ready soon, if not today tomorrow

@scarlehoff
Copy link
Copy Markdown
Member

Great! (I was going to offer doing it myself otherwise :P)

Other than that the PR looks good to me. We can probably merge once we have both the fitbot and your fit.

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 16, 2022

ok, maybe I could add also the corresponding flavour basis where we fit c and cbar?

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 16, 2022

also, the report for validphys should be modified to add v15 I guess, should i do that in this PR?

@scarlehoff
Copy link
Copy Markdown
Member

also, the report for validphys should be modified to add v15 I guess, should i do that in this PR?

Since (for now) the 8-flavour basis is still the default I would leave the comparefits like it is (at least for this PR). This is what we do already with perturvative charm (the plot is in the 8 basis).

We might want to add the plot in whatever the fit basis was whenever the fit is done with a non-default basis though, but I would separate it from this.

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 16, 2022

ok, but here we will have to look at 9 pdfs at the end, so we do need to have an option for plotting V15 I would say? But yes, this can be done in another PR

@scarlehoff
Copy link
Copy Markdown
Member

ok, but here we will have to look at 9 pdfs at the end, so we do need to have an option for plotting V15 I would say? But yes, this can be done in another PR

Personally I would do evolution / flavour / fitbasis (only when it is different) rather than adding V15 everywhere.
For this fit in particular I guess we are only interested in the plot_pdf action with basis the new basis.

Copy link
Copy Markdown
Member

@Radonirinaunimi Radonirinaunimi left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Some minor comments.

Comment thread validphys2/src/validphys/pdfbases.py Outdated
Comment thread validphys2/src/validphys/pdfbases.py Outdated
Comment thread validphys2/src/validphys/pdfbases.py Outdated
@github-actions
Copy link
Copy Markdown

Greetings from your nice fit 🤖 !
I have good news for you, I just finished my tasks:

Check the report carefully, and please buy me a ☕ , or better, a GPU 😉!

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 17, 2022

here s the report for the fit https://vp.nnpdf.science/UMM4lusgR-S8ktV-ie-5iA==/

@scarlehoff
Copy link
Copy Markdown
Member

I think the sum rules need to be actually updated https://vp.nnpdf.science/UMM4lusgR-S8ktV-ie-5iA==/#sum-rules

(but the fit itself is great, the change in the sum rules is small and the chi2 is equivalent!)

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 17, 2022

uhm shouldn t the sumrules be the same tho? Even if we are fitting v15 we should still have \int V = \int V15 = 3 no?

@felixhekhorn
Copy link
Copy Markdown
Contributor

uhm shouldn t the sumrules be the same tho? Even if we are fitting v15 we should still have \int V = \int V15 = 3 no?

yes! there are no charm valence quarks in the proton (same as there are no strange)

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 22, 2022

and here there are some additional plots for c, cbar, c+ and c-
https://vp.nnpdf.science/5FPOitRhQGeXd6RFCuoSnQ==/
https://vp.nnpdf.science/8oaLv7p3QcGZ5ZTW6xx48A==/

i ll fix the tests, need to update the sumrules test in vp2

@scarlehoff
Copy link
Copy Markdown
Member

Oh, wow. It looks great!

We might want to redo it with Th 400 (but maybe we need to undo some of the optimisations done to the fktables!)

Anyway, thanks!

@RoyStegeman
Copy link
Copy Markdown
Member

Do you understand why the agreement with pseudodata (E_val and E_tr) has improved significantly while agreement to the central data remained largely unchanged?

@scarlehoff
Copy link
Copy Markdown
Member

The improvements are within the error though

@RoyStegeman
Copy link
Copy Markdown
Member

RoyStegeman commented Nov 22, 2022

Oh actually nvm, it's still being compared to the fit from 2107.

Compared to a more recent baseline fit, none of the chi2 values have really changed (as you would expect I guess): https://vp.nnpdf.science/y_AjuiH2RHO3QGI9Z97tBQ==/

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 25, 2022

how do I fix the vp2 regression test for the sumrules? I mean, I have to add the etries for c and cbar in the regression test files like test_sum_rules_MC.csv, but where do I take the values for cvalence and cm momentum fraction? From the fit above? Or we need first some baseline fit?

@scarlehoff
Copy link
Copy Markdown
Member

The sumrules are computed so it should be fine just regenerating the .csv file.
Remove the csv files and run pytest and they will be regenerated with the new values. Then you can commit them.

If you want to add (also) a regression test with the ccbar asymmetry (I'd be in favour) you need to create a new low-precision fit for that with theory 162.

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 25, 2022

ok thank you I ll do that. A part from this in this PR there s no other things to do right? How do we proceede to create theory 400?

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Nov 25, 2022

also I was thinking that it might be interesting to implement an additional sumrule for x(c-cbar) using the lattice result from https://www.sciencedirect.com/science/article/pii/S0370269320304366?via%3Dihub and see how this affect the result. But this should probably go in a separate PR I guess

@scarlehoff
Copy link
Copy Markdown
Member

A part from this in this PR there s no other things to do right?

I would say so.

How do we proceede to create theory 400?

@cschwan @andreab1997 the optimizations for the non-independent flavours (when we collapsed V15 into V) is done only to the fktables right? The ekos and the grids contain all the information, so nothing needs to be recomputed, just reapply the ekos with a different set of optimization rules. Right?

@cschwan
Copy link
Copy Markdown
Contributor

cschwan commented Nov 25, 2022

@scarlehoff I believe this is true; we discussed at some point to implement these optimizations at the level of the EKOs, but this was never done (probably better).

@andreab1997
Copy link
Copy Markdown
Contributor

A part from this in this PR there s no other things to do right?

I would say so.

How do we proceede to create theory 400?

@cschwan @andreab1997 the optimizations for the non-independent flavours (when we collapsed V15 into V) is done only to the fktables right? The ekos and the grids contain all the information, so nothing needs to be recomputed, just reapply the ekos with a different set of optimization rules. Right?

Yes I believe this to be the case as well.

@scarlehoff scarlehoff added the run-fit-bot Starts fit bot from a PR. label Nov 25, 2022
Copy link
Copy Markdown
Member

@scarlehoff scarlehoff left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

lgtm, don't know whether @Radonirinaunimi wants to have a last look

from my side once the bot is done I'd be happy with merging it

@Radonirinaunimi
Copy link
Copy Markdown
Member

I've had a second look and everything LGTM.

@github-actions
Copy link
Copy Markdown

Greetings from your nice fit 🤖 !
I have good news for you, I just finished my tasks:

Check the report carefully, and please buy me a ☕ , or better, a GPU 😉!

@scarlehoff
Copy link
Copy Markdown
Member

@tgiani feel free to merge it if you are happy with it

@Zaharid
Copy link
Copy Markdown
Contributor

Zaharid commented Nov 25, 2022

Is this going to change the fits by default? If so, can it be enabled conditionally?

@scarlehoff
Copy link
Copy Markdown
Member

scarlehoff commented Nov 25, 2022

The default fit is going to be (8,) and EVOLUTION. This just adds another option.

Also, now the sum rule for V15 is computed separated from the rest but most of the time the entry for V15 is the entry of V.

(edit, note that the only regression test that has been updated is the one for the sum rules because now there is an extra entry)

@scarlehoff scarlehoff removed the run-fit-bot Starts fit bot from a PR. label Dec 1, 2022
@scarlehoff
Copy link
Copy Markdown
Member

Can we merge this? There's nothing missing right?

@tgiani
Copy link
Copy Markdown
Contributor Author

tgiani commented Dec 1, 2022

ah yes sure, forgot to merge. Doing now

@tgiani tgiani merged commit d7a6dbc into master Dec 1, 2022
@tgiani tgiani deleted the ccbar branch December 1, 2022 08:19
@scarlehoff scarlehoff changed the title [WIP] ccbar asymmetry ccbar asymmetry Dec 1, 2022
@scarlehoff scarlehoff mentioned this pull request Dec 1, 2022
5 tasks
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.

8 participants