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

defination of entries per event #176

Open
Karabo-Mosala opened this issue Feb 8, 2023 · 5 comments
Open

defination of entries per event #176

Karabo-Mosala opened this issue Feb 8, 2023 · 5 comments
Labels
❓question Further information is requested

Comments

@Karabo-Mosala
Copy link

Question

Hi,
I was reading lhco file in MadAnalysis5. Ma5 report gives me 2 numbers "events kept" (given in table of Cut section) and "Integral" (given in table of Histogram section). "Event kept" refers to remaining events after cuts which are generated by MadGraph. These two numbers are related by "entries per event". How conversion between these quantities is done, can you please help?

Thanks

@Karabo-Mosala Karabo-Mosala added the ❓question Further information is requested label Feb 8, 2023
@jackaraz
Copy link
Member

jackaraz commented Feb 8, 2023

Hi @Karabo-Mosala, I'm not sure if I understand your question. Ma5 counts MC events that are passing the cuts, but MC events have no physical meaning, but each has a weight, $w_i$, attached to it. So after cuts, we calculate the probability of events passing the cuts $p = \sum w^{pass}_i / \sum w_i$. So physical events that are passing each cut are $N = p\sigma\mathcal{L}$ where $\sigma$ is the cross-section, and $\mathcal{L}$ is the luminosity. The integral of the histogram is simply the integral of the curve that you are getting.

@Karabo-Mosala
Copy link
Author

Hi,

My question is, looking at the picture attached below - for the histogram, we have entries per event on distribution. How is that value calculated.

Screenshot 2023-02-08 at 11 44 20

@jackaraz
Copy link
Member

jackaraz commented Feb 9, 2023

Hi @Karabo-Mosala entries per event means "sum of event-weights over entries" divided by "sum of event-weights over events".

@Karabo-Mosala
Copy link
Author

Hi Jack,

Thanks for the help.

In some processes this ratio "entries per events" is 1 and in some it is some number.
I have some confusion like;

  1. In my example case, entries per event = 234527/171330 = 1.37
    Here sum of events-weights over enteries = 234527 is physical number, right?
  2. We get this physical number by formula as you explained in first reply. MA5 doesn't show at what luminosity we get this number.
    Thanks for your patience.

@jackaraz
Copy link
Member

Hi @Karabo-Mosala

  1. No, this number is just for you to estimate how sufficient your MC simulation is. In case your particular event needs more events or not. This is not a physical number. The only physical number is the ratio between the initial and final sum of event weights which gives you the efficiency of the cut or histogram bin.
  2. It's a ratio, luminosity cancels out. Columns with events kept and rejected are computed with luminosity.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
❓question Further information is requested
Projects
None yet
Development

No branches or pull requests

2 participants