problem in getting distribution of electron for four leptonic final state

Asked by Deepanjali Goswami

Hi Benjamin,

I am using Minimal Left Right Symmeric model, where I have generated some number of events for purely four leptonic final state in Madgraph. Now I am using this generated events in MA5, also want to plot like M(l1, l2) , M(l2, l3).....etc.
I want to see distributions for only M(e-, e+) or M(mu+, mu-). For this purpose I have defined like

define mu = mu+ mu-
define e = e+ e-
plot M(mu[1] mu[2])
plot M(e[1] e[2])

now, I am getting distribution for M(mu[1] mu[2]) but not for M(e[1] e[2]).

When I checked with
plot N(mu), i can see a distribution with 5 bins and N(e) also gives one bin for electron/positron. I am also not getting any other distribution for electron/positron. I have kept Pyhtia on and Standard cut on in Madgraph while generating no of events. Where the minimum PT for leptons is 10 GeV in standard cut.

I have done it in two steps: first generated root files in MA5 and then used root files for distributions.

I am not understanding what may be the problem. Is the standard cut a problem for not getting a distibution for electron?

And one more thing, I am not getting plot for M(e, mu) also.

Thanking you.
Regards,
Deepanjali

Question information

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

Hi Deepanjali,

From your message, it is clear that you have from 0 to 4 muons/antimuons in each event, and at most one electron/positron in each event. Therefore, it is normal that
  - M(e[1] e[2]) returns an empty plot
  - If M(e[1] mu[1]) is empty, that means that you have no event featuring both electrons and muons.
To me, the issue lies at the level of event generation (including the simulation of the detector) where you may loose many of your electrons/positrons. Please investigate the delphes card.

Cheers,

Benjamin

On 15 Mar 2016, at 12:02 , Deepanjali Goswami <email address hidden> wrote:

> Question #288648 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/288648
>
> Description changed to:
> Hi Benjamin,
>
> I am using Minimal Left Right Symmeric model, where I have generated some number of events for purely four leptonic final state in Madgraph. Now I am using this generated events in MA5, also want to plot like M(l1, l2) , M(l2, l3).....etc.
> I want to see distributions for only M(e-, e+) or M(mu+, mu-). For this purpose I have defined like
>
> define mu = mu+ mu-
> define e = e+ e-
> plot M(mu[1] mu[2])
> plot M(e[1] e[2])
>
> now, I am getting distribution for M(mu[1] mu[2]) but not for M(e[1]
> e[2]).
>
> When I checked with
> plot N(mu), i can see a distribution with 5 bins and N(e) also gives one bin for electron/positron. I am also not getting any other distribution for electron/positron. I have kept Pyhtia on and Standard cut on in Madgraph while generating no of events. Where the minimum PT for leptons is 10 GeV in standard cut.
>
> I have done it in two steps: first generated root files in MA5 and then
> used root files for distributions.
>
> I am not understanding what may be the problem. Is the standard cut a
> problem for not getting a distibution for electron?
>
> And one more thing, I am not getting plot for M(e, mu) also.
>
>
> Thanking you.
> Regards,
> Deepanjali
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Deepanjali Goswami (g-deepanjali-j) said :
#2

Hi Benjamin,

Thanks for your reply.

If one electron or positron is present in one event, then atleast I should get distribution for any kin observable for e+/e- . Am I not correct?
But I am not getting any kin. distribution for e+/e-.

Thanking you.
Regards,
Deepanjali

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

Hi,

No this is not correct. If you want to plot M(e[1] e[2]), you need at least two of these guys.

Cheers,

Benjamin

On 15 Mar 2016, at 12:52 , Deepanjali Goswami <email address hidden> wrote:

> Question #288648 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/288648
>
> Deepanjali Goswami posted a new comment:
> Hi Benjamin,
>
> Thanks for your reply.
>
> If one electron or positron is present in one event, then atleast I should get distribution for any kin observable for e+/e- . Am I not correct?
> But I am not getting any kin. distribution for e+/e-.
>
>
> Thanking you.
> Regards,
> Deepanjali
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Deepanjali Goswami (g-deepanjali-j) said :
#4

Hi Benjamin,

I mean observable like Pt[e+], PT[e-] or E[e+/e-] i should get.

Is it wrong?

Thanks.
Regards,
Deepanjali

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

Hi,

Indeed, you should get something for those. Could you please share the sample? I would like to check what we have in there, in order to understand whether what you get could be expected.

Cheers,

Benjamin

On 16 Mar 2016, at 05:31 , Deepanjali Goswami <email address hidden> wrote:

> Question #288648 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/288648
>
> Deepanjali Goswami posted a new comment:
> Hi Benjamin,
>
> I mean observable like Pt[e+], PT[e-] or E[e+/e-] i should get.
>
> Is it wrong?
>
> Thanks.
> Regards,
> Deepanjali
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Deepanjali Goswami (g-deepanjali-j) said :
#6

Hi Benjamin,

Thanks for your reply.
I am sending you my sample files to your email.

Thanking you.
Regards,
Deepanjali

Revision history for this message
Deepanjali Goswami (g-deepanjali-j) said :
#7

Thanks Benjamin Fuks, that solved my question.