S vs B

Asked by SOUAD SEMLALI

Dear MadAnalysis team,

I am using S/B (main.fom.formula=1) to check the signal events vs. the background events after applying some cuts, but the chart shows the sum of weights of signal and background after each cut (myregion.saf) and the ratio of these weights (S/B), instead of the number of events scaled to luminosity. Please correct me if I am missing something.

From myregion.saf:
-----------
Signal:
-------------
initial number of events: 1000000
sum of weights : 482

1st cut:
nentries : 921899
sum of weights: 445.26

------------------
background
------------------
initial number of events: 1000000
sum of weights : 2.95x 10^6

1st cut:
nentries : 133836
sum of weights: 3.94 x10^5

===========================================

------------------------
From the chart:
---------------------
intial no cuts: signal (482); background (2950987) ;; S vs B (0.000164)
first cut: signal (445); background (394948) ;; S vs B (1.13x10^-3)

Many thanks for your help!

Best ,
Souad

Question information

Language:
English Edit question
Status:
Solved
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Solved by:
Jack Y. Araz
Solved:
Last query:
Last reply:
Revision history for this message
Benjamin Fuks (fuks) said :
#1

Hi Souad,

There is nothing wrong here. The weights must be included to account for the different signal and background cross sections.

Cheers,

Benjamin

Revision history for this message
SOUAD SEMLALI (souad-sem) said :
#2

Hi Benjamin,

Thanks a lot for your reply. I am afraid I have another question. The
number of events (nentries) kept after applying the cut of the signal is
higher than the number of events of the background, but the ratio in the
table is very suppressed, which might be explained by the difference in the
cross-section. but, why the chart doesn't show the number of events (or
the normalized events) kept after each cut?
Many thanks for your help to alleviate the confusion.

Best
Souad

Le mer. 5 janv. 2022 à 10:56, Benjamin Fuks <
<email address hidden>> a écrit :

> Your question #700094 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/700094
>
> Status: Open => Answered
>
> Benjamin Fuks proposed the following answer:
> Hi Souad,
>
> There is nothing wrong here. The weights must be included to account for
> the different signal and background cross sections.
>
> Cheers,
>
> Benjamin
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/700094/+confirm?answer_id=0
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/700094
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Revision history for this message
Benjamin Fuks (fuks) said :
#3

Hi Souad,

We need to include cross section information to get a proper information on a number of events (for a given sample). The number of entries is indeed related to how many events you generated, which is not the information we want to see in a cutflow. For this reason, we show the *weighted* number of events passing each cut.

Cheers,

Benjamin

Revision history for this message
SOUAD SEMLALI (souad-sem) said :
#4

Hi Benjamin,

Thank you for your clarification! and apologies for dragging this for so
long. I understand now that the cut-flow chart shows the weighted number
of events after each cut.

I only have one last question. There is an example ( top-antitop, WW, ZZ
=signal with one leptonic and one invisible decay ) in one of the
tutorials
https://indico.cern.ch/event/568875/contributions/2397928/attachments/1459090/2253214/fuks_MC4BSM17.pdf,
slide =20

After applying the cut on MET < 125.0 GeV, the same table as the one in
the slide is produced (signal = 640). but the weighted events were
different (2.31x10^3)
"1_MET < 125.0" # 1st cut
9283 0 # nentries
2.316294e+03 0.000000e+00 # sum of weights
5.779617e+02 0.000000e+00 # sum of weights^2

Could you please clarify why the table doesn't show the number of the
weighted events, which is 2.31x10^3 instead of 640?

Thanks again for your help!

Best
Souad

Le jeu. 6 janv. 2022 à 21:35, Benjamin Fuks <
<email address hidden>> a écrit :

> Your question #700094 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/700094
>
> Status: Open => Answered
>
> Benjamin Fuks proposed the following answer:
> Hi Souad,
>
> We need to include cross section information to get a proper information
> on a number of events (for a given sample). The number of entries is
> indeed related to how many events you generated, which is not the
> information we want to see in a cutflow. For this reason, we show the
> *weighted* number of events passing each cut.
>
> Cheers,
>
> Benjamin
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madanalysis5/+question/700094/+confirm?answer_id=2
>
> If you still need help, you can reply to this email or go to the
> following page to enter your feedback:
> https://answers.launchpad.net/madanalysis5/+question/700094
>
> You received this question notification because you asked the question.
>

--
Souad SEMLALI

Revision history for this message
Best Jack Y. Araz (jackaraz) said :
#5

Hi Suad

The saf files are not final calculation those are the raw info for you to take and calculate what ever you need. For instance you might need to calculate the number of events with different xsec or different lumi. If we give you number of events you wont be able to use saf files again hence you will have to run it multiple times. So refer saf files as raw data where the html/latex report is the processed one where we calculate number of events for you with respect to the xsec and lumi that you have provided.

cheers


Jack Y. Araz

> On Jan 7, 2022, at 10:41 AM, SOUAD SEMLALI <email address hidden> wrote:
>
> [EXTERNAL EMAIL]
>
> Question #700094 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/700094
>
> Status: Answered => Open
>
> SOUAD SEMLALI is still having a problem:
> Hi Benjamin,
>
> Thank you for your clarification! and apologies for dragging this for so
> long. I understand now that the cut-flow chart shows the weighted number
> of events after each cut.
>
> I only have one last question. There is an example ( top-antitop, WW, ZZ
> =signal with one leptonic and one invisible decay ) in one of the
> tutorials
> https://indico.cern.ch/event/568875/contributions/2397928/attachments/1459090/2253214/fuks_MC4BSM17.pdf,
> slide =20
>
> After applying the cut on MET < 125.0 GeV, the same table as the one in
> the slide is produced (signal = 640). but the weighted events were
> different (2.31x10^3)
> "1_MET < 125.0" # 1st cut
> 9283 0 # nentries
> 2.316294e+03 0.000000e+00 # sum of weights
> 5.779617e+02 0.000000e+00 # sum of weights^2
>
> Could you please clarify why the table doesn't show the number of the
> weighted events, which is 2.31x10^3 instead of 640?
>
> Thanks again for your help!
>
> Best
> Souad
>
>
> Le jeu. 6 janv. 2022 à 21:35, Benjamin Fuks <
> <email address hidden>> a écrit :
>
>> Your question #700094 on MadAnalysis 5 changed:
>> https://answers.launchpad.net/madanalysis5/+question/700094
>>
>> Status: Open => Answered
>>
>> Benjamin Fuks proposed the following answer:
>> Hi Souad,
>>
>> We need to include cross section information to get a proper information
>> on a number of events (for a given sample). The number of entries is
>> indeed related to how many events you generated, which is not the
>> information we want to see in a cutflow. For this reason, we show the
>> *weighted* number of events passing each cut.
>>
>> Cheers,
>>
>> Benjamin
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>>
>> https://answers.launchpad.net/madanalysis5/+question/700094/+confirm?answer_id=2
>>
>> If you still need help, you can reply to this email or go to the
>> following page to enter your feedback:
>> https://answers.launchpad.net/madanalysis5/+question/700094
>>
>> You received this question notification because you asked the question.
>>
>
>
> --
> Souad SEMLALI
>
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
SOUAD SEMLALI (souad-sem) said :
#6

Thanks Jack Y. Araz, that solved my question.