Skip to content

Write the right nf in the LHAPDF .info file.#281

Merged
felixhekhorn merged 5 commits into
masterfrom
fix-for-NumFlavors-in-info-file
Jun 14, 2023
Merged

Write the right nf in the LHAPDF .info file.#281
felixhekhorn merged 5 commits into
masterfrom
fix-for-NumFlavors-in-info-file

Conversation

@scarlehoff
Copy link
Copy Markdown
Member

@scarlehoff scarlehoff commented Jun 6, 2023

Currently the NumFlavors written by eko in the info file correspond to the number of possible partons, while this field is actually the maximum nf used for the evolution of alpha_s.

Ref: https://lhapdf.hepforge.org/config.html

@scarlehoff scarlehoff requested a review from andreab1997 June 6, 2023 14:46
Comment thread src/ekobox/info_file.py Outdated
@felixhekhorn felixhekhorn added the bug Something isn't working label Jun 6, 2023
@andreab1997
Copy link
Copy Markdown
Collaborator

Ok for me but what about the Flavors list? It should be consistent with NumFlavors right?

@felixhekhorn
Copy link
Copy Markdown
Collaborator

Currently the NumFlavors written by eko in the info file correspond to the number of possible partons, while this field is actually the maximum nf used for the evolution of alpha_s.

Ref: https://lhapdf.hepforge.org/config.html

not quite - quoting said HTML page

Maximum/fixed number of flavors used in the fit.

so this refers to the PDFs in the first place - so it should be theory_card.heavy.num_flavs_max_pdf (of course 14 was just wrong)

alpha_s evolution will use this value if AlphaS_NumFlavors is not explicitly set.

this is only a fallback

PS: @alecandido I know you want to give up on that specific variable (did we do it already?), so we would need to provide an alternative

@felixhekhorn
Copy link
Copy Markdown
Collaborator

Ok for me but what about the Flavors list? It should be consistent with NumFlavors right?

No - again quoting from the officals:

List of supported (user-visible) parton flavors, using the PDG Monte Carlo numbering scheme (quarks from 1 .. 6, antiquarks from -6 .. -1, and gluon = 21, as well as any "optional extras" such as photon = 22).

so you can stuff what ever there (also BSM) - NumFlavors refers to DGLAP

@scarlehoff
Copy link
Copy Markdown
Member Author

so this refers to the PDFs in the first place

Oh, I missed the . I understood "on the fit evolution" as one object.

In any case, it is not 14 :P

Ok for me but what about the Flavors list? It should be consistent with NumFlavors right?

Not really since we are just filling with 0s the numbers we don't care about in Flavors.

Comment thread src/ekobox/info_file.py Outdated
Comment thread src/ekobox/info_file.py Outdated
@felixhekhorn felixhekhorn merged commit e81f315 into master Jun 14, 2023
@felixhekhorn felixhekhorn deleted the fix-for-NumFlavors-in-info-file branch June 14, 2023 16:05
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

bug Something isn't working

Projects

None yet

Development

Successfully merging this pull request may close these issues.

4 participants