Definitions of cuts

Asked by Sarah

Hello
I am new user of Madgraph and I have questions
can you answer my questions please,
1- I don't know what are these cuts in Madgraph pt/E/eta/dr/mij ?
2- I want to generate these processes in 2HDM type II q q~ > A > (Z* off shell ) h, Z*> e- e+ and from gluons g g~ > A > (Z* off shell ) h, Z*> e- e+.
with applying cuts on charged leptons |ηl| < 2.5 and plT > 10 GeV. my questions are
Do I have to write the full process or it is ok to write q q~or (g g~)> A > h e- e+ and Madgraph will understand there is z* ?
Also, because I don't know more about cuts, how can I apply these cuts ?

I am really Sorry if I ask questions they were answered before.

Best regards
Sarah

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

 1- I don't know what are these cuts in Madgraph pt/E/eta/dr/mij ?

pt: transverse momenta of the particle
E: Energy of the particle in the PARTONIC center of mass frame of the collision (not the lab frame for p p)
eta: rapidity
mij : invariant mass cut between two particles

2- I want to generate these processes in 2HDM type II q q~ > A > (Z* off shell ) h, Z*> e- e+ and from gluons g g~ > A > (Z* off shell ) h, Z*> e- e+.

Here it is difficult to know for me the exact syntax that you want since I do not know the model and the hyppothesis that you can use for your process.
I would do
generate p p > A, A > h e+ e-
or
generate p p > A, A > Z > h e+ e-
if you want (and can) to force the presence of a S-channel Z in the A onshell decay. (Please check the gauge invariance and lorentz invariance when you use such dangerous syntax).

To apply cuts, you have to defined those in the file run_card.dat.
Be carefull about the cut_decays parameter and be sure to have it set on True otherwise since your electron came from the onshell decay of the "A", your cut will not be applied.

Cheers,

Olivier

Revision history for this message
Sarah (sarahwarad) said :
#2

thank you Olivier
I noticed there is no A particle in the 2HDM particles. is it h1?

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

I guess it does depend of the 2hdm model that you use. I can only suggest you to check the model reference paper.

Cheers,

Olivier

Revision history for this message
Sarah (sarahwarad) said :
#4

(Here it is difficult to know for me the exact syntax that you want since I do not know the model and the hyppothesis that you can use for your process.)

I meant 2HDM model

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

I guess that you mean this one then:
http://feynrules.irmp.ucl.ac.be/wiki/2HDM

Since if i do
import model 2HDM I see:
MG5_aMC>import model 2HDM
INFO: download model from http://feynrules.irmp.ucl.ac.be/raw-attachment/wiki/2HDM/2HDM_UFO.tar.gz

Cheers,

Olivier

Revision history for this message
Scott Dressler (gressler000) said :
#6

Often, one wants to combine more than one process in the same event sample. For example, you can do combined top and bottom production:

    generate p p > t t~
    add process p p > b b~

Note that if you do another "generate", it will just overwrite the first generate. That's why you have to do "add process" instead.

Now you can output the resulting combined processes:

    output MY_SECOND_RUN

and launch:

    launch MY_SECOND_RUN

The resulting events generated will be correctly weighted according to their individual cross-sections. This time we will only be interested in the cross-sections, and not want to look at plots of kinematic variables. So we won't run Pythia, or either of the detector simulators. So when Madgraph asks what switches you want to turn on, you say

    0

which means none. Madgraph will now give you the option of editing the parameter card, the run card, or the plot card. We'll just be looking at the total cross-section, so we just need a very quick run. So let's edit the run card to only ask for 1000 events. Type

    2

and then change "10000" to "1000" on line 32. Save the file, and then type

    0

to start the run.

When Madgraph is finished, you should see that the cross-section is something like 4.41e+08 pb, which is huge! To see the individual contributions, you can type:https://www.atketo.com

    open index.html

and click on "Results and Event Database". If you click on the cross-section number link, it will bring up a page that lists all of the individual partonic subprocesses, and how much they contributed. You can see that almost all of the cross-section is due to gluon+gluon > bottom+antibottom.

Can you help with this problem?

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

To post a message you must log in.