something about h>zz>jjll

Asked by duc

Hi,

    I am trying to generate p p > H > h h, h > b b~, h > z z > q q l l process in 2HDM(contain heft), and there are something strange in the distribution of z mass in parton level.

    Here is my script:
"define l = e+ mu+ ta+ e- mu- ta-
define l+ = e+ mu+ ta+
define l- = e- mu- ta-
define q = u d c s u~ d~ c~ s~
generate p p > h2, ( h2 > h h, h > b b~, ( h > q q z /q, ( z > l+ l-)))"

where I use "h > q q z /q, ( z > l+ l-)" to make sure on shell z decays to ll, off shell z decays to qq. Now the question is, when I reconstruct the invariant mass of l+l-, in addition to a pick around 90GeV as expect, there is another higher pick near 55GeV, I don't understand what had happened...

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

What is your value of bwcutoff?
With the default at 15, you should not be able to go down up to 55 GeV. (But if you change the Z width).

Cheers,

Olivier

On Nov 3, 2014, at 9:46 AM, duc <email address hidden> wrote:

> New question #256695 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Hi,
>
> I am trying to generate p p > H > h h, h > b b~, h > z z > q q l l process in 2HDM(contain heft), and there are something strange in the distribution of z mass in parton level.
>
> Here is my script:
> "define l = e+ mu+ ta+ e- mu- ta-
> define l+ = e+ mu+ ta+
> define l- = e- mu- ta-
> define q = u d c s u~ d~ c~ s~
> generate p p > h2, ( h2 > h h, h > b b~, ( h > q q z /q, ( z > l+ l-)))"
>
> where I use "h > q q z /q, ( z > l+ l-)" to make sure on shell z decays to ll, off shell z decays to qq. Now the question is, when I reconstruct the invariant mass of l+l-, in addition to a pick around 90GeV as expect, there is another higher pick near 55GeV, I don't understand what had happened...
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
duc (chun-thazen-du) said :
#2

hi,

I use the default value of the Z width and bwcutoff(15). The shape of
the invariant
mass of l+l- is like this in fact: there is no value below about
55GeV, suddenly
there is a stage at 55GeV, then it goes down until 90GeV where there is a
peak as expect...

2014-11-03 20:16 GMT+08:00 Olivier Mattelaer <
<email address hidden>>:

> Your question #256695 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
> What is your value of bwcutoff?
> With the default at 15, you should not be able to go down up to 55 GeV.
> (But if you change the Z width).
>
> Cheers,
>
> Olivier
>
> On Nov 3, 2014, at 9:46 AM, duc <email address hidden>
> wrote:
>
> > New question #256695 on MadGraph5_aMC@NLO:
> > https://answers.launchpad.net/mg5amcnlo/+question/256695
> >
> > Hi,
> >
> > I am trying to generate p p > H > h h, h > b b~, h > z z > q q l l
> process in 2HDM(contain heft), and there are something strange in the
> distribution of z mass in parton level.
> >
> > Here is my script:
> > "define l = e+ mu+ ta+ e- mu- ta-
> > define l+ = e+ mu+ ta+
> > define l- = e- mu- ta-
> > define q = u d c s u~ d~ c~ s~
> > generate p p > h2, ( h2 > h h, h > b b~, ( h > q q z /q, ( z > l+ l-)))"
> >
> > where I use "h > q q z /q, ( z > l+ l-)" to make sure on shell z decays
> to ll, off shell z decays to qq. Now the question is, when I reconstruct
> the invariant mass of l+l-, in addition to a pick around 90GeV as expect,
> there is another higher pick near 55GeV, I don't understand what had
> happened...
> >
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadGraph5_aMC@NLO.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/mg5amcnlo/+question/256695/+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/mg5amcnlo/+question/256695
>
> You received this question notification because you asked the question.
>

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#3

Hi,

Then what you observe is the cut in the “onshell” definition
indeed 55 = 91(MZ) - 15 (bwcutoff) * 2.44 (Width)

So in your case you probably want to have a smaller value for bwcutoff.

Cheers,

Olivier

On Nov 4, 2014, at 2:07 AM, duc <email address hidden> wrote:

> Question #256695 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Status: Answered => Open
>
> duc is still having a problem:
> hi,
>
> I use the default value of the Z width and bwcutoff(15). The shape of
> the invariant
> mass of l+l- is like this in fact: there is no value below about
> 55GeV, suddenly
> there is a stage at 55GeV, then it goes down until 90GeV where there is a
> peak as expect...
>
>
> 2014-11-03 20:16 GMT+08:00 Olivier Mattelaer <
> <email address hidden>>:
>
>> Your question #256695 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/256695
>>
>> Status: Open => Answered
>>
>> Olivier Mattelaer proposed the following answer:
>> Hi,
>>
>> What is your value of bwcutoff?
>> With the default at 15, you should not be able to go down up to 55 GeV.
>> (But if you change the Z width).
>>
>> Cheers,
>>
>> Olivier
>>
>> On Nov 3, 2014, at 9:46 AM, duc <email address hidden>
>> wrote:
>>
>>> New question #256695 on MadGraph5_aMC@NLO:
>>> https://answers.launchpad.net/mg5amcnlo/+question/256695
>>>
>>> Hi,
>>>
>>> I am trying to generate p p > H > h h, h > b b~, h > z z > q q l l
>> process in 2HDM(contain heft), and there are something strange in the
>> distribution of z mass in parton level.
>>>
>>> Here is my script:
>>> "define l = e+ mu+ ta+ e- mu- ta-
>>> define l+ = e+ mu+ ta+
>>> define l- = e- mu- ta-
>>> define q = u d c s u~ d~ c~ s~
>>> generate p p > h2, ( h2 > h h, h > b b~, ( h > q q z /q, ( z > l+ l-)))"
>>>
>>> where I use "h > q q z /q, ( z > l+ l-)" to make sure on shell z decays
>> to ll, off shell z decays to qq. Now the question is, when I reconstruct
>> the invariant mass of l+l-, in addition to a pick around 90GeV as expect,
>> there is another higher pick near 55GeV, I don't understand what had
>> happened...
>>>
>>>
>>> --
>>> You received this question notification because you are an answer
>>> contact for MadGraph5_aMC@NLO.
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>>
>> https://answers.launchpad.net/mg5amcnlo/+question/256695/+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/mg5amcnlo/+question/256695
>>
>> You received this question notification because you asked the question.
>>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
duc (chun-thazen-du) said :
#4

Hi,

Thank you for last replay, I regenerated this process these days. As you said, when I used bwcutoff=5GeV, there is just one pick around 91GeV for the invariant mass of l+l-. But at the same time, there is also one pick around 91GeV for the invariant mass of qq(remember that my script is "h2 > h h, h > b b~, ( h > q q z /q, ( z > l+ l-))", z to qq should be off shell, so there may not be a pick around 91GeV... ). What do you think?

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#5

What is the width of h and h2 for your model?

Cheers,

Olivier
On Nov 9, 2014, at 12:46 PM, duc <email address hidden> wrote:

> Question #256695 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Status: Answered => Open
>
> duc is still having a problem:
> Hi,
>
> Thank you for last replay, I regenerated this process these days. As you
> said, when I used bwcutoff=5GeV, there is just one pick around 91GeV for
> the invariant mass of l+l-. But at the same time, there is also one pick
> around 91GeV for the invariant mass of qq(remember that my script is
> "h2 > h h, h > b b~, ( h > q q z /q, ( z > l+ l-))", z to qq should be
> off shell, so there may not be a pick around 91GeV... ). What do you
> think?
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
duc (chun-thazen-du) said :
#6

The mass of h2 is 300GeV, the width is 5GeV; the mass of h is 125GeV, the width is 5GeV.

Thanks!

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#7

Hi,

Looks like that you have a very large width for your Higgs.
Then your second peak should be linked to a combination of the Higgs off-shel (above it’s mass) and the second Z off shell (below it’s mass).

Cheers,

Olivier

On Nov 11, 2014, at 3:12 PM, duc <email address hidden> wrote:

> Question #256695 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Status: Answered => Open
>
> duc is still having a problem:
> The mass of h2 is 300GeV, the width is 5GeV; the mass of h is 125GeV,
> the width is 5GeV.
>
> Thanks!
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
duc (chun-thazen-du) said :
#8

sorry, the width of 125GeV Higgs is 0.0058GeV in my model, I made a typo... The width of 300GeV Higgs is 5GeV.
so you still suggest me to reduce that width?...

Thanks a lot!

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#9

Hi,

Then I do not understand. I have made the following:
import heft
generate g g > h, (h > j j z , z > l+ l- )

and the invariant mass of the two jet is strongly peaked at around 45GeV as expected.
Did you make your plot at parton level?

Cheers,

Olivier

On Nov 12, 2014, at 1:51 AM, duc <email address hidden> wrote:

> Question #256695 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Status: Answered => Open
>
> duc is still having a problem:
> sorry, the width of 125GeV Higgs is 0.0058GeV in my model, I made a typo... The width of 300GeV Higgs is 5GeV.
> so you still suggest me to reduce that width?...
>
> Thanks a lot!
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
duc (chun-thazen-du) said :
#10

ok, I will try this and plot the parton level first, and what's the
bwcutoff you used?

2014-11-12 17:51 GMT+08:00 Olivier Mattelaer <
<email address hidden>>:

> Your question #256695 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
> Then I do not understand. I have made the following:
> import heft
> generate g g > h, (h > j j z , z > l+ l- )
>
> and the invariant mass of the two jet is strongly peaked at around 45GeV
> as expected.
> Did you make your plot at parton level?
>
> Cheers,
>
> Olivier
>
>
> On Nov 12, 2014, at 1:51 AM, duc <email address hidden>
> wrote:
>
> > Question #256695 on MadGraph5_aMC@NLO changed:
> > https://answers.launchpad.net/mg5amcnlo/+question/256695
> >
> > Status: Answered => Open
> >
> > duc is still having a problem:
> > sorry, the width of 125GeV Higgs is 0.0058GeV in my model, I made a
> typo... The width of 300GeV Higgs is 5GeV.
> > so you still suggest me to reduce that width?...
> >
> > Thanks a lot!
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadGraph5_aMC@NLO.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/mg5amcnlo/+question/256695/+confirm?answer_id=8
>
> 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/mg5amcnlo/+question/256695
>
> You received this question notification because you asked the question.
>

Revision history for this message
duc (chun-thazen-du) said :
#11

hi, I have generate that process as you did, I use the default bwcutoff.
But the invariant mass of l+l- doesn't have the on shell pick around 91GeV
as my plot at parton level...

2014-11-12 17:51 GMT+08:00 Olivier Mattelaer <
<email address hidden>>:

> Your question #256695 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi,
>
> Then I do not understand. I have made the following:
> import heft
> generate g g > h, (h > j j z , z > l+ l- )
>
> and the invariant mass of the two jet is strongly peaked at around 45GeV
> as expected.
> Did you make your plot at parton level?
>
> Cheers,
>
> Olivier
>
>
> On Nov 12, 2014, at 1:51 AM, duc <email address hidden>
> wrote:
>
> > Question #256695 on MadGraph5_aMC@NLO changed:
> > https://answers.launchpad.net/mg5amcnlo/+question/256695
> >
> > Status: Answered => Open
> >
> > duc is still having a problem:
> > sorry, the width of 125GeV Higgs is 0.0058GeV in my model, I made a
> typo... The width of 300GeV Higgs is 5GeV.
> > so you still suggest me to reduce that width?...
> >
> > Thanks a lot!
> >
> > --
> > You received this question notification because you are an answer
> > contact for MadGraph5_aMC@NLO.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/mg5amcnlo/+question/256695/+confirm?answer_id=8
>
> 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/mg5amcnlo/+question/256695
>
> You received this question notification because you asked the question.
>

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#12

Hi Duc,

Did you mean that the plot are fine at parton level?
If yes this means that MG is correct and that you need to understand if your problem comes from the shower and/or of your detector simulation program and then get in contact with the associate author.

Cheers,

Olivier

Revision history for this message
duc (chun-thazen-du) said :
#13

Hi Olivier, the invariant mass of the two jet is fine, but the invariant
mass of l+ l- is just like the two jet case, only one peak around 45 GeV,
they should be on shell in our case...

2014-11-13 0:31 GMT+08:00 Olivier Mattelaer <
<email address hidden>>:

> Your question #256695 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> Hi Duc,
>
> Did you mean that the plot are fine at parton level?
> If yes this means that MG is correct and that you need to understand if
> your problem comes from the shower and/or of your detector simulation
> program and then get in contact with the associate author.
>
> Cheers,
>
> Olivier
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/mg5amcnlo/+question/256695/+confirm?answer_id=11
>
> 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/mg5amcnlo/+question/256695
>
> You received this question notification because you asked the question.
>

Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said :
#14

Sorry but I’m completely confuse by your question.
What is the problem?

Did you understand the effect of the bwcutoff parameter?
If no please take a look at one of the MG5 tutorial (to the answer of one of them)
for example:
https://cp3.irmp.ucl.ac.be/projects/madgraph/attachment/wiki/Lund2014/14_09_23_tutomg_Lund.pdf

Cheers,

Olivier
On Nov 14, 2014, at 4:36 AM, duc <email address hidden> wrote:

> Question #256695 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/256695
>
> Status: Answered => Open
>
> duc is still having a problem:
> Hi Olivier, the invariant mass of the two jet is fine, but the invariant
> mass of l+ l- is just like the two jet case, only one peak around 45 GeV,
> they should be on shell in our case...
>
> 2014-11-13 0:31 GMT+08:00 Olivier Mattelaer <
> <email address hidden>>:
>
>> Your question #256695 on MadGraph5_aMC@NLO changed:
>> https://answers.launchpad.net/mg5amcnlo/+question/256695
>>
>> Status: Open => Answered
>>
>> Olivier Mattelaer proposed the following answer:
>> Hi Duc,
>>
>> Did you mean that the plot are fine at parton level?
>> If yes this means that MG is correct and that you need to understand if
>> your problem comes from the shower and/or of your detector simulation
>> program and then get in contact with the associate author.
>>
>> Cheers,
>>
>> Olivier
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>>
>> https://answers.launchpad.net/mg5amcnlo/+question/256695/+confirm?answer_id=11
>>
>> 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/mg5amcnlo/+question/256695
>>
>> You received this question notification because you asked the question.
>>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Can you help with this problem?

Provide an answer of your own, or ask duc for more information if necessary.

To post a message you must log in.