mij cuut

Asked by ahh

Hi MG Authors
iam trying to calculate acut for process
p p > w z j j ,z > e- e+ , w > j j
with cuts in run file :

#*********************************************************************
# Minimum and maximum invariant mass for pairs *
#*********************************************************************
 200 = mmjj ! min invariant mass of a jet pair
 0 = mmbb ! min invariant mass of a b pair
 0 = mmaa ! min invariant mass of gamma gamma pair
 200 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
 500 = mmjjmax ! max invariant mass of a jet pair
 -1 = mmbbmax ! max invariant mass of a b pair
 -1 = mmaamax ! max invariant mass of gamma gamma pair
 500 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
#*********************************************************************
# Minimum and maximum invariant mass for all letpons *
#*********************************************************************
 200 = mmnl ! min invariant mass for all letpons (l+- and vl)
 500 = mmnlmax ! max invariant mass for all letpons (l+- and vl)

but it gives me error ???!!!
Cheers
Ahmed

Question information

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

Hi Ahmed,

You have a problem of definition between various cuts in your run_card.dat:

You use decay chain syntax for the w/z (z > e- e+ , w > j j ). Those syntax are valid if the z/w are on shell.
This force by the cut linked to bwcutoff (default 15)
| M_Z^* - M_Z| < bwcutoff * \Gamma_Z
the same for the W mass

So this constraints the MMLL value and MMJJ to a maximum value.

Note that for the Z, you don't have any reason to not include the photon contribution (and the associate interference!)
at such high energy.

Cheers,

Olivier

On Nov 26, 2013, at 3:31 PM, ahmed hammad <email address hidden> wrote:

> New question #239850 on MadGraph5:
> https://answers.launchpad.net/madgraph5/+question/239850
>
> Hi MG Authors
> iam trying to calculate acut for process
> p p > w z j j ,z > e- e+ , w > j j
> with cuts in run file :
>
>
>
> #*********************************************************************
> # Minimum and maximum invariant mass for pairs *
> #*********************************************************************
> 200 = mmjj ! min invariant mass of a jet pair
> 0 = mmbb ! min invariant mass of a b pair
> 0 = mmaa ! min invariant mass of gamma gamma pair
> 200 = mmll ! min invariant mass of l+l- (same flavour) lepton pair
> 500 = mmjjmax ! max invariant mass of a jet pair
> -1 = mmbbmax ! max invariant mass of a b pair
> -1 = mmaamax ! max invariant mass of gamma gamma pair
> 500 = mmllmax ! max invariant mass of l+l- (same flavour) lepton pair
> #*********************************************************************
> # Minimum and maximum invariant mass for all letpons *
> #*********************************************************************
> 200 = mmnl ! min invariant mass for all letpons (l+- and vl)
> 500 = mmnlmax ! max invariant mass for all letpons (l+- and vl)
>
>
>
> but it gives me error ???!!!
> Cheers
> Ahmed
>
> --
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.

Revision history for this message
ahh (ahh) said :
#2

Hi dr Olivier
the way to set the cut , that i want to study the tail of background processes in high energy .
So , by seeting cut 200 ~500 and cut > 500 , i will be sure that MG will generate 100000 events at that region .
So , is there any way to overcome this ????
thanks for your help
Ahmed

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

My suggestion would be to generate:
p p > e+ e- j j j j
with your set of cuts. But for this you need a cluster.

Otherwise, you can do what you did and set bwcutoff to 1000 in the run_card.dat
This should make the generation to go trough but I don't think that the associate results can be consider as physical.

Cheers,

olivier

On Nov 26, 2013, at 4:56 PM, ahmed hammad <email address hidden> wrote:

> Question #239850 on MadGraph5 changed:
> https://answers.launchpad.net/madgraph5/+question/239850
>
> Status: Answered => Open
>
> ahmed hammad is still having a problem:
> Hi dr Olivier
> the way to set the cut , that i want to study the tail of background processes in high energy .
> So , by seeting cut 200 ~500 and cut > 500 , i will be sure that MG will generate 100000 events at that region .
> So , is there any way to overcome this ????
> thanks for your help
> Ahmed
>
> --
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.

Revision history for this message
ahh (ahh) said :
#4

Hi dr Olivier
First of all i want to thank you for your interest
But ,excuse my small information, i just want to ask about things .
i was using calchep for long and i know that the calculation done using NWA
in such aprocess "p p > w z ,z > e- e+ , w > j j"
so CHep calculate it cross section* BR , but it missed the spin correlation
, here W/Z are on-shell ,So , CHep cannot set any cut for it
i tried a lot to know the mechanism of how MG calculate such a process ???
2- whey you set such a cut for the tail of the resonant particles "bwcutoff"
???
3-finally , what did you mean by "non-physical"?
thanks alot for your patient .
Cheers
Ahmed

On Tue, Nov 26, 2013 at 7:16 PM, Olivier Mattelaer <
<email address hidden>> wrote:

> Your question #239850 on MadGraph5 changed:
> https://answers.launchpad.net/madgraph5/+question/239850
>
> Status: Open => Answered
>
> Olivier Mattelaer proposed the following answer:
> My suggestion would be to generate:
> p p > e+ e- j j j j
> with your set of cuts. But for this you need a cluster.
>
> Otherwise, you can do what you did and set bwcutoff to 1000 in the
> run_card.dat
> This should make the generation to go trough but I don't think that the
> associate results can be consider as physical.
>
> Cheers,
>
> olivier
>
>
> On Nov 26, 2013, at 4:56 PM, ahmed hammad <
> <email address hidden>> wrote:
>
> > Question #239850 on MadGraph5 changed:
> > https://answers.launchpad.net/madgraph5/+question/239850
> >
> > Status: Answered => Open
> >
> > ahmed hammad is still having a problem:
> > Hi dr Olivier
> > the way to set the cut , that i want to study the tail of background
> processes in high energy .
> > So , by seeting cut 200 ~500 and cut > 500 , i will be sure that MG
> will generate 100000 events at that region .
> > So , is there any way to overcome this ????
> > thanks for your help
> > Ahmed
> >
> > --
> > You received this question notification because you are a member of
> > MadTeam, which is an answer contact for MadGraph5.
>
> --
> If this answers your question, please go to the following page to let us
> know that it is solved:
>
> https://answers.launchpad.net/madgraph5/+question/239850/+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/madgraph5/+question/239850
>
> You received this question notification because you asked the question.
>

--
Ahmed Hammad
Center for Theoretical Physics (CTP)
*Institute of Basic Sciences (IBS)*
Zewail City of Science and Technology

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

Hi Ahmed,

> , here W/Z are on-shell ,So , CHep cannot set any cut for it
> i tried a lot to know the mechanism of how MG calculate such a process ???

We use the full matrix-element without any particular approximation.
So we have the full width effect and the full spin-correlation. As a correlated effect our
cross-section is not (exactly) give by cross section* BR.

They are however one approximation that you do, is to select a sub-sample of diagram
for a given final state.
In your example
p p > w z, z > e- e+ , w > j j
your final state is
e+ e- j j
they are other diagram which contribute to such final state including
p p > a w > e+ e- j j

Due to your syntax the second process will not be computed.
In general, those other process not only contribute but also interferes and in some case those
interferences restore gauge invariance. So neglect those are in general dangerous.
On the other hand, if the massive particles are on-shell then we know that the on-shell diagram dominates
and therefore the interferences (and the matrix-element square of the others diagrams) are negligible.
This is why we set this bwcutoff since the approximation breaks down if we allow the Z/W to be too-far off-shell.

Please take a look at the following tutorial, where this is explained in detail:
https://cp3.irmp.ucl.ac.be/projects/madgraph/attachment/wiki/TASISchool13/13_06_10_tutomg_tasi.pdf

Cheers,

Olivier

On Nov 27, 2013, at 10:16 AM, ahmed hammad <email address hidden> wrote:

> Question #239850 on MadGraph5 changed:
> https://answers.launchpad.net/madgraph5/+question/239850
>
> Status: Answered => Open
>
> ahmed hammad is still having a problem:
> Hi dr Olivier
> First of all i want to thank you for your interest
> But ,excuse my small information, i just want to ask about things .
> i was using calchep for long and i know that the calculation done using NWA
> in such aprocess "p p > w z ,z > e- e+ , w > j j"
> so CHep calculate it cross section* BR , but it missed the spin correlation
> , here W/Z are on-shell ,So , CHep cannot set any cut for it
> i tried a lot to know the mechanism of how MG calculate such a process ???
> 2- whey you set such a cut for the tail of the resonant particles "bwcutoff"
> ???
> 3-finally , what did you mean by "non-physical"?
> thanks alot for your patient .
> Cheers
> Ahmed
>
>
> On Tue, Nov 26, 2013 at 7:16 PM, Olivier Mattelaer <
> <email address hidden>> wrote:
>
>> Your question #239850 on MadGraph5 changed:
>> https://answers.launchpad.net/madgraph5/+question/239850
>>
>> Status: Open => Answered
>>
>> Olivier Mattelaer proposed the following answer:
>> My suggestion would be to generate:
>> p p > e+ e- j j j j
>> with your set of cuts. But for this you need a cluster.
>>
>> Otherwise, you can do what you did and set bwcutoff to 1000 in the
>> run_card.dat
>> This should make the generation to go trough but I don't think that the
>> associate results can be consider as physical.
>>
>> Cheers,
>>
>> olivier
>>
>>
>> On Nov 26, 2013, at 4:56 PM, ahmed hammad <
>> <email address hidden>> wrote:
>>
>>> Question #239850 on MadGraph5 changed:
>>> https://answers.launchpad.net/madgraph5/+question/239850
>>>
>>> Status: Answered => Open
>>>
>>> ahmed hammad is still having a problem:
>>> Hi dr Olivier
>>> the way to set the cut , that i want to study the tail of background
>> processes in high energy .
>>> So , by seeting cut 200 ~500 and cut > 500 , i will be sure that MG
>> will generate 100000 events at that region .
>>> So , is there any way to overcome this ????
>>> thanks for your help
>>> Ahmed
>>>
>>> --
>>> You received this question notification because you are a member of
>>> MadTeam, which is an answer contact for MadGraph5.
>>
>> --
>> If this answers your question, please go to the following page to let us
>> know that it is solved:
>>
>> https://answers.launchpad.net/madgraph5/+question/239850/+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/madgraph5/+question/239850
>>
>> You received this question notification because you asked the question.
>>
>
>
> --
> Ahmed Hammad
> Center for Theoretical Physics (CTP)
> *Institute of Basic Sciences (IBS)*
> Zewail City of Science and Technology
>
> You received this question notification because you are a member of
> MadTeam, which is an answer contact for MadGraph5.

Revision history for this message
ahh (ahh) said :
#6

Hi Dr Olivier
i really apperciate that , that clearify alot .
Cheers
Ahmed