How to genrate 2HDM pseudoscalar A via b- associated production?

Asked by Khawla Jaffel on 2020-08-27

Dear Madgraph experts,

I want to produce a grid pack for 2HDM signal samples( H->ZA->ll +bb).
The process that I want to generate are:
1- quark , antiquark annihilation
2- gluon-gluon fusion
3- b associated production

For know i have somthing like this, but as far i undrestand ( by displayying the diagrams using directly madgraph ) it will give me only the first process( quark-antiquark ) :

set group_subprocesses Auto
set ignore_six_quark_processes False
set loop_optimized_output True
set complex_mass_scheme False
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model 2HDM4MG5-may15
set group_subprocesses Auto
set ignore_six_quark_processes False
set gauge unitary
set complex_mass_scheme False
define l+ = e+ mu+
define l- = e- mu-
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
generate p p > h3 z/h1, h3 > b b~, z > l+ l-
output HToZATo2L2B_1000_200

If I wan to include the gluon fusion, keeping the quark anti-quark as well doing something like this can work?
I mean using
generate p p > h3 z/h1 HG=1, h3 > b b~, z > l+ l-
instead of
generate p p > h3 z/h1, h3 > b b~, z > l+ l-

As for the b-associated production process, I have no idea how to include it, any guide?
is there's a way to include these 3 process or the associated production have to be separate?

Many thanks,
Khawla

Question information

Language:
English Edit question
Status:
Answered
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Last query:
2020-08-28
Last reply:
2020-08-28

If you want to do a five flavor computation (b quark as collinear in the initial state and linked to a PDF).
Then you have to set the masse of the b-quark to zero at the model level see
https://answers.launchpad.net/mg5amcnlo/+faq/2312

for one way of doing that.
with that model the syntax that you used before will include the b as part of the proton automatically

Cheers,

Olivier

> On 27 Aug 2020, at 15:55, Khawla Jaffel <email address hidden> wrote:
>
> New question #692630 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/692630
>
> Dear Madgraph experts,
>
> I want to produce a grid pack for 2HDM signal samples( H->ZA->ll +bb).
> The process that I want to generate are:
> 1- quark , antiquark annihilation
> 2- gluon-gluon fusion
> 3- b associated production
>
> For know i have somthing like this, but as far i undrestand ( by displayying the diagrams using directly madgraph ) it will give me only the first process( quark-antiquark ) :
>
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set loop_optimized_output True
> set complex_mass_scheme False
> import model sm
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s u~ c~ d~ s~
> define l+ = e+ mu+
> define l- = e- mu-
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> import model 2HDM4MG5-may15
> set group_subprocesses Auto
> set ignore_six_quark_processes False
> set gauge unitary
> set complex_mass_scheme False
> define l+ = e+ mu+
> define l- = e- mu-
> define p = g u c d s u~ c~ d~ s~
> define j = g u c d s u~ c~ d~ s~
> define vl = ve vm vt
> define vl~ = ve~ vm~ vt~
> generate p p > h3 z/h1, h3 > b b~, z > l+ l-
> output HToZATo2L2B_1000_200
>
> If I wan to include the gluon fusion, keeping the quark anti-quark as well doing something like this can work?
> I mean using
> generate p p > h3 z/h1 HG=1, h3 > b b~, z > l+ l-
> instead of
> generate p p > h3 z/h1, h3 > b b~, z > l+ l-
>
> As for the b-associated production process, I have no idea how to include it, any guide?
> is there's a way to include these 3 process or the associated production have to be separate?
>
> Many thanks,
> Khawla
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Khawla Jaffel (kjaffel) said : #2

Hello Olivier,

thank you very much for your reply.
I try to follow the instructions that you provide me.
what I did is :
instead of using restrict_XXXX.dat in the UFO model, generated by the python script write_param_card.py in models/sm directory ( since it has to be a valid param_card.dat )
I just use the default one, because when I launch mg5_aMc I see it's already loaded:

[kjaffel@ingrid-ui1 MG5_aMC_v2_5_5]$ ./bin/mg5_aMC
blabla...
Loading default model: sm
INFO: Restrict model sm with file models/sm/restrict_default.dat .

what I change in the default file is the MB value in the Block mass. i set it to zero.
    5 0.000000e+00 # MB

then using the same syntax that I send you before:
I see that I switch to 5F -scheme, which is good for now ...
MG5_aMC>import model sm
INFO: Restrict model sm with file models/sm/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Pass the definition of 'j' and 'p' to 5 flavor scheme.

However, I can't see any difference in the diagrams:
I still get only gg-fusion and the quark-antiquark

here the link of madgraph output, where you can see the diagrams and the card used.
file:///run/user/1002/gvfs/sftp:host=ingrid1-gwceci/home/users/k/j/kjaffel/MG5_aMC_v2_5_5/HToZATo2L2B_200_50_ggfusion_qq_andbassociatedproduction/index.html

can you check with me, please, if what I am doing is correct?

Many thanks in advance,
Khawla

But if you do

import model 2HDM4MG5-may15

this will not applied anymore, you need to use the restrict_default.dat from the
2HDM4MG5-may15
UFO model
(or from 2HDM4MG5 and then use the restrict_may15.dat file)

Cheers,

Olivier
> On 28 Aug 2020, at 14:01, Khawla Jaffel <email address hidden> wrote:
>
> Question #692630 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/692630
>
> Status: Answered => Open
>
> Khawla Jaffel is still having a problem:
> Hello Olivier,
>
> thank you very much for your reply.
> I try to follow the instructions that you provide me.
> what I did is :
> instead of using restrict_XXXX.dat in the UFO model, generated by the python script write_param_card.py in models/sm directory ( since it has to be a valid param_card.dat )
> I just use the default one, because when I launch mg5_aMc I see it's already loaded:
>
> [kjaffel@ingrid-ui1 MG5_aMC_v2_5_5]$ ./bin/mg5_aMC
> blabla...
> Loading default model: sm
> INFO: Restrict model sm with file models/sm/restrict_default.dat .
>
> what I change in the default file is the MB value in the Block mass. i set it to zero.
> 5 0.000000e+00 # MB
>
> then using the same syntax that I send you before:
> I see that I switch to 5F -scheme, which is good for now ...
> MG5_aMC>import model sm
> INFO: Restrict model sm with file models/sm/restrict_default.dat .
> INFO: Run "set stdout_level DEBUG" before import for more information.
> INFO: Change particles name to pass to MG5 convention
> Pass the definition of 'j' and 'p' to 5 flavor scheme.
>
> However, I can't see any difference in the diagrams:
> I still get only gg-fusion and the quark-antiquark
>
> here the link of madgraph output, where you can see the diagrams and the card used.
> file:///run/user/1002/gvfs/sftp:host=ingrid1-gwceci/home/users/k/j/kjaffel/MG5_aMC_v2_5_5/HToZATo2L2B_200_50_ggfusion_qq_andbassociatedproduction/index.html
>
> can you check with me, please, if what I am doing is correct?
>
> Many thanks in advance,
> Khawla
>
> --
> 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 Khawla Jaffel for more information if necessary.

To post a message you must log in.