recognizing collinear distribution type

[t7phy]

When producing grids with a FF instead of a PDF, we need to trick PineAPPL faking DIS to make it produce actual result.
Specifically:

1. we have to trick the meta data by making the initial states 211 (pion pid) and 11, instead of the real 11 and -11.
2. we have to trick the subgrids by setting x2 as 1.0 (as in the DIS case) to make pineappl make only single convolution (instead of 2).

this problem will also persist for "higher collinear distribution dimension":

• single spin asymmetries: using one polarized and one unpolarized PDF
• SIDIS: one PDF = one FF
• etc.

sincerely,
@felixhekhorn

PS: the relevant line at which we crashed is

pineappl/pineappl/src/grid.rs

Line 559 in b1036ae

lumi_cache.setup(self, xi).unwrap();

[alecandido]

More than automated recognition, I'd say we need to explicitly specify the number of distributions involved in the PineAPPL grid explicitly.

Do we really need to make PineAPPL aware of the kind of the distribution? It should not matter whether they are two PDFs or one PDF and one FF (or two FFs). For the sake of interpolation, only the number of distributions should be relevant. Correct? @felixhekhorn @cschwan

(of course you could always store the type in metadata, but this would be a convention on PineAPPL usage, no need to make PineAPPL aware and hardcoded)

[cschwan]

For these applications the metadata isn't properly named, so think of initial_state_1 and initial_state_2 not as the initial states, but rather the type of particle that is described by the function with which a convolution is performed. For proton PDFs you obviously have to have it set to 2212 and in general it must match the Particle field of the function's .info file in LHAPDF.

See also #135.

[t7phy]

@cschwan how about also adding a final_state_1 key? This could also make it clear whether the convolution is with a PDF or FF and ultimately the 3 convolutions (2 pdfs and 1 ffs) would also benefit from it.

[alecandido]

I'd rather name them hadron_1, hadron_2, ... otherwise you will have to support identical combinations.
Whether it is used for initial or final state should be up to the user (and you could describe with further metadata).

[felixhekhorn]

1. Should we close this as duplicate of Clarify initial states and PDFs used #135?
2. indeed the way to go would be something along the lines of: initial_state_1: 2212 -> particle_1: {pid: 2212, polarized: False, timelike: False} - if you wish "timelike" can be replaced by "is_outgoing" or similar ... this would have to be repeated eventually up to 3 times

[cschwan]

Closed in favor of #135.