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

New correlation weights management #39

Closed
eugene274 opened this issue Oct 8, 2020 · 1 comment
Closed

New correlation weights management #39

eugene274 opened this issue Oct 8, 2020 · 1 comment
Labels
help wanted Extra attention is needed

Comments

@eugene274
Copy link
Contributor

Hi!

I've recently found that you refactored management of weights for correlation task. So now I have to pass two functions one for correlation itself and another one which evaluates weights.

I have few questions concerning underlying math.
Previously we had OBSERVABLE type and REFERENCE type of Q-vectors. Is there a representation of those types in new version?
You added weights functions into CorrelationFunctions.hpp. Maybe I'm stupid, but I cannot catch new notation (n2, nd2,etc). Can you clarify this? I see also that there is a common ansatz for weights of the form M(M-1). Can you also say few words on that?

Thanks a lot!

@eugene274 eugene274 added the help wanted Extra attention is needed label Oct 8, 2020
@kreisl
Copy link
Contributor

kreisl commented Oct 8, 2020

The observable and reference type are still there now it is named use_weights_. If use_weights_ is true the resulting correlation will be of type OBSERVABLE.
We now need to pass a second function which takes care of the weights. There are some premade ones for the cumulant methods.
For example:
v_2{2} = sqrt(c_2{2})
c_2{2} = <C_2> <> is event average

Edit:
In the case of my functions c2 = C_2 per event. So n2 is already divided out.
However It is still needed during the event average which is weighted by n2.
<C_2> = (sum_i^N C_2 n2) / (sum_i^N n2)

So n2 is equivalent to the weight / multiplicity m(m-1) the Q-vector
nd2 is supposed to be the same for differential cumulants, which make the distinction between particles of interest p and reference particles r. d2 is the corresponding numerator for differential cumulants.
it is a similar notation to: https://arxiv.org/pdf/1312.3572.pdf
Maybe we should adopt the same notation and link to this paper to make it easier to understand, what the purpose of these functions are.

@kreisl kreisl closed this as completed Oct 29, 2020
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
help wanted Extra attention is needed
Projects
None yet
Development

No branches or pull requests

2 participants