MadAnalysis 5

multi-particle

Asked by audiya fatihah

Hi,

im studying an event which decays to 4 leptons final state,
though i am not sure whether im getting this right,
i defined 4mu as
> define mumu = mu+ mu- mu+ mu-
also 4l as
>define ll = l+ l- l+ l-
>plot M (mumu) 150 0 150
>plot M(ll) 150 0 150
the plot seems odd, i think im not getting this right,

also what is the difference say,
>define mu = mu- mu+
> plot M (mu) 100 0 150
> plot M (mu+ mu-) 100 0 150

can anyone suggest a solution

thanks!

Question information

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

Hi,

>> define mumu = mu+ mu- mu+ mu-
> also 4l as
>> define ll = l+ l- l+ l-
>> plot M (mumu) 150 0 150
This does not do what you aim to do. Here, you ask the code to plot the invariant mass of any single muon present in the event (with one entry for each muon). In order to do what you want, one needs to type
 plot M(mu+[1] mu-[1] mu-[2] mu-[2])

> also what is the difference say,
>> define mu = mu- mu+
>> plot M (mu) 100 0 150
>> plot M (mu+ mu-) 100 0 150
With the second comment, one combines the momentum of the two muons before extracting the invariant mass. With the first comment, one gets one entry for each muon present in the event (with the value being its invariant mass).

I hope this helps.

Regards,

Benjamin

Revision history for this message
audiya fatihah (audiyafatihah) said : 		#2

Hi benjamin thanks for your reply, im so sorry for the basics quest, i already went through the manual & still seem to not get it right

so i applied some cuts,
ma5> plot M( l+[1] l-[1] l+[2] l-[2])
ma5> reject M ( e+[1] e-[1] e+[2] e-[2]) < 5
ma5> reject ETA ( e+[1] e-[1] e+[2] e-[2]) > 2.5
ma5> reject M ( mu+[1] mu-[1] mu+[2] mu-[2]) < 7
ma5> reject ETA ( mu+[1] mu-[1] mu+[2] mu-[2]) > 2.4
ma5> select 120 < M (l+[1] l-[1] l+[2] l-[2]) < 130
ma5>set main.lumi = 30
ma5> set main.stacking_method = normalize2one

the summary table came out fine, but there are no plots, (not a empty plot but no plots at all)

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

Hi,

Do you mind sharing the events and the ma5 script. Note that it is very likely I won't answer before August 26th (am about to leave for vacation).

Regards,

Benjamin

> On 6 Aug 2018, at 14:37 , audiya fatihah <email address hidden> wrote:
>
> Question #671350 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/671350
>
> Status: Answered => Open
>
> audiya fatihah is still having a problem:
> Hi benjamin thanks for your reply, im so sorry for the basics quest, i
> already went through the manual & still seem to not get it right
>
> so i applied some cuts,
> ma5> plot M( l+[1] l-[1] l+[2] l-[2])
> ma5> reject M ( e+[1] e-[1] e+[2] e-[2]) < 5
> ma5> reject ETA ( e+[1] e-[1] e+[2] e-[2]) > 2.5
> ma5> reject M ( mu+[1] mu-[1] mu+[2] mu-[2]) < 7
> ma5> reject ETA ( mu+[1] mu-[1] mu+[2] mu-[2]) > 2.4
> ma5> select 120 < M (l+[1] l-[1] l+[2] l-[2]) < 130
> ma5>set main.lumi = 30
> ma5> set main.stacking_method = normalize2one
>
> the summary table came out fine, but there are no plots, (not a empty
> plot but no plots at all)
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
audiya fatihah (audiyafatihah) said : 		#4

hi, i tried to restart the session, only then the plot appeared
Thanks!