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

Interface QED and N3LO, n.2 #316

Merged
merged 26 commits into from
Oct 19, 2023
Merged

Interface QED and N3LO, n.2 #316

merged 26 commits into from
Oct 19, 2023

Conversation

niclaurenti
Copy link
Contributor

reopening #311 after rebasing master

@niclaurenti
Copy link
Contributor Author

Now this is more or less ready. @alecandido I wasn't able to reconcile the qed branches as you asked. I think that the uglier branch is the one happening in

if order[1] == 0:

but it is quite complicated to change it since qcd and qed evolution require different ingredients and also the structure is a bit different(in one case we have singlet + non_singlet, in the other also the valence). If you have proposals we can try to restructure it a bit.

Then, whenever we will tag the next eko version, we have to pay attention to the fact that i remover the argument qed from apply_pdf (that information now is contained in the EKO that is passed to the function). The only part that is affected should be
https://github.com/NNPDF/nnpdf/blob/db8b790c1bacf52f71cc3a7fabf0fd9674cbf476/n3fit/src/evolven3fit_new/evolve.py#L178
but maybe some other parts of the NNPDF code are using this function (maybe @andreab1997 knows better)

Copy link
Collaborator

@giacomomagni giacomomagni left a comment

Choose a reason for hiding this comment

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

Looks good to me. to_flavor_basis_tensor still contains a qed=False as default which can be dropped?

src/ekobox/apply.py Show resolved Hide resolved
src/eko/basis_rotation.py Show resolved Hide resolved
@niclaurenti
Copy link
Contributor Author

Looks good to me. to_flavor_basis_tensor still contains a qed=False as default which can be dropped?

Thanks I missed it. I will fix it asap

@felixhekhorn felixhekhorn added physics new physics features refactor Refactor code labels Oct 12, 2023
Copy link
Contributor

@felixhekhorn felixhekhorn left a comment

Choose a reason for hiding this comment

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

Did you check that the apfel benchmark is still fine? else this seems alright!

@niclaurenti
Copy link
Contributor Author

Did you check that the apfel benchmark is still fine? else this seems alright!

@felixhekhorn The apfel bench at order = (3, 2) is still good. However, I cannot test order=(4,2) since apfel doesn't implement N3LO (I think)

@niclaurenti niclaurenti merged commit f0e00e9 into master Oct 19, 2023
5 checks passed
@niclaurenti niclaurenti deleted the n3lo_plus_qed2 branch October 19, 2023 12:06
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
physics new physics features refactor Refactor code
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants