Question #818653 “ggF Higgs LO and NLO in SM” : Questions : MadGraph5

Revision history for this message

Olivier Mattelaer (olivier-mattelaer) said on 2024-09-04:

#1

Concerning:
> InvalidCmd : No particle h, in model
you are missing a space between the "h" and the ","

Now the "," (onshell decay) are not supported for loop induced processes.
For this particular case (higgs decay), the solution is to use madspin with the line "set spinmode=none"
WHile this mode is bad for Z/top decay (and so many particle) it is perfect for Higgs (spin 0 and small width).

Concerning the NLO contribution, we are not able to do two loop conputation within our framework and therefore we do not have NLO.

What we can do is "an approximate NLO" like in this paper: https://inspirehep.net/literature/1279493
So here you got to NLO in the heft framework and then you do a special reweighting to introduce back the top mass effect.
Obviously this is only approximate since we do miss the two loop computation and have not the correct reweighting factor for term related to the loop.

Cheers,

Olivier

Revision history for this message

Kalash (r-b-99) said on 2024-09-05:

#2

Hello,
Thank you so much for responding.

> you are missing a space between the "h" and the ","

I fixed it, and still gave me an error (even after set spinmode=none):

MG5_aMC>generate g g > h , h > a a QED<=1 [noborn = QCD]
The current model sm does not allow to generate loop corrections of type ['QCD'].
MG5_aMC now loads 'loop_sm'.
import model loop_sm
INFO: Restrict model loop_sm with file models/loop_sm/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
Defined multiparticle all = g gh gh~ d u s c d~ u~ s~ c~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ b t b~ t~ z w+ h w- ta- ta+
Command "generate g g > h , h > a a QED<=1 [noborn = QCD]" interrupted with error:
InvalidCmd : No particle , in model

Tried this other notation:

MG5_aMC>generate g g > h > a a QED<=1 [noborn = QCD]
Command "generate g g > h > a a QED<=1 [noborn = QCD]" interrupted with error:
NoDiagramException : No amplitudes generated from process Process: g g > h> a a QED=1 [ noborn = QCD ]. Please enter a valid process
MG5_aMC>

Regarding the NLO calculation, thank you so much. Just to check:
I need to get to the NLO in heft, then do the reweighting to get the NLO in sm.
In this tutorial (https://www.niu.edu/spmartin/madgraph/madhiggs.html) it is LO heft, correct?
and when trying to add [QCD] to the heft process, I get this:

MG5_aMC>generate p p > h , h > a a [QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate p p > h , h > a a [QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes
MG5_aMC>generate p p > h [QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate p p > h [QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes
MG5_aMC>generate p p > h [noborn=QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate p p > h [noborn=QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes
MG5_aMC>generate p p > h [QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate p p > h [QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes
MG5_aMC>generate g g > h [QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate g g > h [QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes

Can you advise on how to generate NLO processes in HEFT?

Also, when I perform the same heft generating process in the tutorial but for WW decay channel instead for the diphoton, I get an error after the "launch":

MG5_aMC>import model heft_raghad
INFO: Restrict model heft_raghad with file models/heft_raghad/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h h1 w- ta- ta+
MG5_aMC>generate p p > h, h > w+ w-
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: g g > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Process has 1 diagrams
INFO: Trying process: u u~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Trying process: u c~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Trying process: c u~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Trying process: c c~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Trying process: d d~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Trying process: d s~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Trying process: s d~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Trying process: s s~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2
INFO: Process has 1 diagrams
1 processes with 2 diagrams generated in 0.035 s
Total: 1 processes with 2 diagrams
MG5_aMC>output pp_h_ww_heft
INFO: initialize a new directory: pp_h_ww_heft
INFO: remove old information in pp_h_ww_heft
INFO: Organizing processes into subprocess groups
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
INFO: Generating Helas calls for process: h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2
INFO: Combine g g > h HIG<=1 HIW<=1 WEIGHTED<=2 @1 with decays h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2
INFO: Processing color information for process: g g > h HIG<=1 HIW<=1 @1
  Decay: h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2
INFO: Creating files in directory P1_gg_h_h_wpwm
INFO: Generating Feynman diagrams for Process: g g > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
  Decay: h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2
INFO: Finding symmetric diagrams for subprocess group gg_h_h_wpwm
Generated helas calls for 1 subprocesses (1 diagrams) in 0.005 s
Wrote files for 6 helas calls in 0.031 s
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 2 routines in 0.353 s
save configuration file to /home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft/Cards/me5_configuration.txt
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
INFO: Generating MadAnalysis5 default cards tailored to this process
^[[AOutput to directory /home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft done.
Type "launch" to generate events from this process, or see
/home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft/README
Run "open index.html" to see more information about this process.
MG5_aMC>launch pp_h_ww_heft
************************************************************
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 3.5.5 2024-06-17 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft/Cards/me5_configuration.txt
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/input/mg5_configuration.txt
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft/Cards/me5_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
generate_events run_01
The following switches determine which programs are run:
/=================== Description ===================|================== values ===================|===== other options =====\
| 1. Choose the shower/hadronization program | shower = OFF | Pythia8 |
| 2. Choose the detector simulation program | detector = OFF | Delphes |
| 3. Choose an analysis package (plot/convert) | analysis = MadAnalysis5 | Rivet|OFF |
| 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
| 5. Add weights to events for new hypp. | reweight = OFF | ON |
\===========================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=Pythia8' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>0
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. madanalysis5_parton : madanalysis5_parton_card.dat |
\------------------------------------------------------------/
you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
[0, done, 1, param, 2, run, 3, madanalysis5_parton, enter path][90s to answer]
>0
INFO: Update the dependent parameter of the param_card.dat
WARNING: update the strong coupling value (alpha_s) to the value from the pdf selected: 0.13
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
compile Source Directory
Using random number seed offset = 21
INFO: Running Survey
Creating Jobs
Working on SubProcesses
INFO: Compiling for process 1/1.
INFO: P1_gg_h_h_wpwm
INFO: P1_gg_h_h_wpwm
Zero result detected: No Phase Space. Please check particle masses.

INFO:
quit
INFO:
Generation failed (no results.dat file found)
MG5_aMC>

Sorry for the long reply, and thank you in advance,

Raghad.

Hello,
Thank you so much for responding.

> you are missing a space between the "h" and the ","

I fixed it, and still gave me an error (even after set spinmode=none):

MG5_aMC>generate g g > h , h > a a QED<=1 [noborn = QCD]
The current model sm does not allow to generate loop corrections of type ['QCD'].
MG5_aMC now loads 'loop_sm'.
 import model loop_sm
INFO: Restrict model loop_sm with file models/loop_sm/restrict_default.dat . 
INFO: Run "set stdout_level DEBUG" before import for more information. 
INFO: Change particles name to pass to MG5 convention 
Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
Defined multiparticle all = g gh gh~ d u s c d~ u~ s~ c~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ b t b~ t~ z w+ h w- ta- ta+
Command "generate g g > h , h > a a QED<=1 [noborn = QCD]" interrupted with error:
InvalidCmd : No particle , in model

Tried this other notation:

MG5_aMC>generate g g > h > a a QED<=1 [noborn = QCD]
Command "generate g g > h > a a QED<=1 [noborn = QCD]" interrupted with error:
NoDiagramException : No amplitudes generated from process Process: g g > h> a a QED=1 [ noborn = QCD ]. Please enter a valid process
MG5_aMC>

Regarding the NLO calculation, thank you so much. Just to check:
I need to get to the NLO in heft, then do the reweighting to get the NLO in sm. 
In this tutorial (https://www.niu.edu/spmartin/madgraph/madhiggs.html) it is LO heft, correct?
and when trying to add [QCD] to the heft process, I get this:

MG5_aMC>generate p p > h , h > a a [QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate p p > h , h > a a [QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes
MG5_aMC>generate p p > h [QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate p p > h [QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes
MG5_aMC>generate p p > h [noborn=QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate p p > h [noborn=QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes
MG5_aMC>generate p p > h [QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate p p > h [QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes
MG5_aMC>generate g g > h [QCD]
The current model heft_raghad does not allow to generate loop corrections of type ['QCD'].
Command "generate g g > h [QCD]" interrupted with error:
InvalidCmd : The model heft_raghad cannot handle loop processes

Can you advise on how to generate NLO processes in HEFT?

Also, when I perform the same heft generating process in the tutorial but for WW decay channel instead for the diphoton, I get an error after the "launch":

MG5_aMC>import model heft_raghad
INFO: Restrict model heft_raghad with file models/heft_raghad/restrict_default.dat . 
INFO: Run "set stdout_level DEBUG" before import for more information. 
INFO: Change particles name to pass to MG5 convention 
Kept definitions of multiparticles p / j / l+ / l- / vl / vl~ unchanged
Defined multiparticle all = g u c d s u~ c~ d~ s~ a ve vm vt e- mu- ve~ vm~ vt~ e+ mu+ t b t~ b~ z w+ h h1 w- ta- ta+
MG5_aMC>generate p p > h, h > w+ w-
INFO: Checking for minimal orders which gives processes. 
INFO: Please specify coupling orders to bypass this step. 
INFO: Trying process: g g > h HIG<=1 HIW<=1 WEIGHTED<=2 @1  
INFO: Process has 1 diagrams 
INFO: Trying process: u u~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1  
INFO: Trying process: u c~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1  
INFO: Trying process: c u~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1  
INFO: Trying process: c c~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1  
INFO: Trying process: d d~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1  
INFO: Trying process: d s~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1  
INFO: Trying process: s d~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1  
INFO: Trying process: s s~ > h HIG<=1 HIW<=1 WEIGHTED<=2 @1  
INFO: Checking for minimal orders which gives processes. 
INFO: Please specify coupling orders to bypass this step. 
INFO: Trying process: h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2  
INFO: Process has 1 diagrams 
1 processes with 2 diagrams generated in 0.035 s
Total: 1 processes with 2 diagrams
MG5_aMC>output pp_h_ww_heft
INFO: initialize a new directory: pp_h_ww_heft 
INFO: remove old information in pp_h_ww_heft 
INFO: Organizing processes into subprocess groups 
INFO: Organizing processes into subprocess groups 
INFO: Generating Helas calls for process: g g > h HIG<=1 HIW<=1 WEIGHTED<=2 @1 
INFO: Generating Helas calls for process: h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2 
INFO: Combine g g > h HIG<=1 HIW<=1 WEIGHTED<=2 @1 with decays h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2 
INFO: Processing color information for process: g g > h HIG<=1 HIW<=1 @1
  Decay: h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2 
INFO: Creating files in directory P1_gg_h_h_wpwm 
INFO: Generating Feynman diagrams for Process: g g > h HIG<=1 HIW<=1 WEIGHTED<=2 @1
  Decay: h > w+ w- HIG<=1 HIW<=1 WEIGHTED<=2 
INFO: Finding symmetric diagrams for subprocess group gg_h_h_wpwm 
Generated helas calls for 1 subprocesses (1 diagrams) in 0.005 s
Wrote files for 6 helas calls in 0.031 s
ALOHA: aloha starts to compute helicity amplitudes
ALOHA: aloha creates 2 routines in  0.353 s
save configuration file to /home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft/Cards/me5_configuration.txt
INFO: Use Fortran compiler gfortran 
INFO: Use c++ compiler g++ 
INFO: Generate jpeg diagrams 
INFO: Generate web pages 
INFO: Generating MadAnalysis5 default cards tailored to this process 
^[[AOutput to directory /home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft done.
Type "launch" to generate events from this process, or see
/home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft/README
Run "open index.html" to see more information about this process.
MG5_aMC>launch pp_h_ww_heft
************************************************************
*                                                          *
*                      W E L C O M E to                    *
*             M A D G R A P H 5 _ a M C @ N L O            *
*                      M A D E V E N T                     *
*                                                          *
*                 *                       *                *
*                   *        * *        *                  *
*                     * * * * 5 * * * *                    *
*                   *        * *        *                  *
*                 *                       *                *
*                                                          *
*         VERSION 3.5.5                 2024-06-17         *
*                                                          *
*    The MadGraph5_aMC@NLO Development Team - Find us at   *
*    https://server06.fynu.ucl.ac.be/projects/madgraph     *
*                                                          *
*               Type 'help' for in-line help.              *
*                                                          *
************************************************************
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft/Cards/me5_configuration.txt  
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/input/mg5_configuration.txt  
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/pp_h_ww_heft/Cards/me5_configuration.txt  
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
generate_events run_01
The following switches determine which programs are run:
/=================== Description ===================|================== values ===================|===== other options =====\
| 1. Choose the shower/hadronization program        |        shower = OFF                         |     Pythia8             |
| 2. Choose the detector simulation program         |      detector = OFF                         |     Delphes             |
| 3. Choose an analysis package (plot/convert)      |      analysis = MadAnalysis5                |     Rivet|OFF           |
| 4. Decay onshell particles                        |       madspin = OFF                         |     ON|onshell|full     |
| 5. Add weights to events for new hypp.            |      reweight = OFF                         |     ON                  |
\===========================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=Pythia8' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer] 
>0
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
|  1. param               : param_card.dat                   |
|  2. run                 : run_card.dat                     |
|  3. madanalysis5_parton : madanalysis5_parton_card.dat     |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, madanalysis5_parton, enter path][90s to answer] 
>0
INFO: Update the dependent parameter of the param_card.dat 
WARNING: update the strong coupling value (alpha_s) to the value from the pdf selected: 0.13 
Generating 10000 events with run name run_01
survey  run_01 
INFO: compile directory 
compile Source Directory
Using random number seed offset = 21
INFO: Running Survey 
Creating Jobs
Working on SubProcesses
INFO: Compiling for process 1/1. 
INFO:     P1_gg_h_h_wpwm  
INFO:     P1_gg_h_h_wpwm  
Zero result detected:  No Phase Space. Please check particle masses.

INFO:  
quit
INFO:  
Generation failed (no results.dat file found)
MG5_aMC>

Sorry for the long reply, and thank you in advance,

Raghad.

Revision history for this message

Sihyun Jeon (shjeon) said on 2024-09-05:

#3

The first question where you get error

InvalidCmd : No particle , in model
generate g g > h , h > a a QED<=1 [noborn = QCD]

this needs to change into
generate g g > h, h > a a QED<=1 [noborn = QCD]

"h , " -> "h, " without space between higgs and the comma

Revision history for this message

Sihyun Jeon (shjeon) said on 2024-09-05:

#4

With the second error,

MG5_aMC>generate g g > h > a a QED<=1 [noborn = QCD]
Command "generate g g > h > a a QED<=1 [noborn = QCD]" interrupted with error:
NoDiagramException : No amplitudes generated from process Process: g g > h> a a QED=1 [ noborn = QCD ]. Please enter a valid process

This just implies there is no interaction modeled for "h -> aa". [noborn=QCD] is used for calculating the loop for the "production part of the higgs", it's not possible to model the loop for the "decay of higgs", is what I think is happening.

Olivier can correct me if I am wrong.

Revision history for this message

Sihyun Jeon (shjeon) said on 2024-09-05:

#5

Third question,

Regarding the NLO calculation, thank you so much. Just to check:
I need to get to the NLO in heft, then do the reweighting to get the NLO in sm.
In this tutorial (https://www.niu.edu/spmartin/madgraph/madhiggs.html) it is LO heft, correct?
and when trying to add [QCD] to the heft process, I get this:

I am not sure what model this is, same model in http://madgraph.phys.ucl.ac.be/Downloads/models/heft.tgz ?
But in any case, not all UFO models are NLO available, you have to find the model that allows NLO calculations at NLO and the model you are plugging in seems to be not providing it.

If you look around google, I see
https://feynrules.irmp.ucl.ac.be/wiki/HiggsCharacterisation

Revision history for this message

Sihyun Jeon (shjeon) said on 2024-09-05:

#6

Last question ( hope I didn't miss anything in the middle)

Also, when I perform the same heft generating process in the tutorial but for WW decay channel instead for the diphoton, I get an error after the "launch":

Answer is same as this.
https://answers.launchpad.net/mg5amcnlo/+question/818646

Your madspin should be written
spinmode=none
decay h > w+ w-

You should not give the full matrix element " p p > h, h > w+ w- "
You should first give "generate p p > h" production part
and then decay higgs using madspin as above

Revision history for this message

Sihyun Jeon (shjeon) said on 2024-09-05:

#7

-

Revision history for this message

Olivier Mattelaer (olivier-mattelaer) said on 2024-09-06:

#8

Concerning
MG5_aMC>generate g g > h , h > a a QED<=1 [noborn = QCD]

The issue here is that such computation is a two loop computation, which can not be evaluated by our code.
The solution is to generate
generate g g > h [noborn=QCD]
and let madspin do the decay (potentially in the heft model or again in the sm)
decay h > a a [noborn=QCD]
or
decay h > a a [noborn=QED]
or
decay h > a a [noborn=QED QCD]

Concerning
h> w+ w-
In top of what Sihyun mentioned, note that you need to decay the W.
This is important since, like this, you ask an onshell Higgs to decay to two onshell W which is impossible.
so you need to use the syntax:
decay h > l+l- vl vl~

Cheers,

Olivier

Revision history for this message

Sihyun Jeon (shjeon) said on 2024-09-06:

#9

ah yes thanks for the correction, my bad

h > w+ w- should not work in madspin for the same reason = mass restrictions

Revision history for this message

Kalash (r-b-99) said on 2024-09-07:

#10

Hello,

Thank you so much for responding.

When trying the decay command, I get this:

MG5_aMC>generate g g > h [noborn=QCD]
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=2: WEIGTHED IS QCD+2*QED
INFO: Trying coupling order WEIGHTED<=3: WEIGTHED IS QCD+2*QED
INFO: Contributing diagrams generated: 0 Born, 4 loops, 2 R2, 0 UV
1 processes with 4 diagrams generated in 0.071 s
Total: 1 processes with 4 diagrams
MG5_aMC> decay h > a a
Command "decay" not recognized, please try again
MG5_aMC>

I hope I am using it as requested.
Can you help?

Many thanks,

Raghad.

Revision history for this message

Sihyun Jeon (shjeon) said on 2024-09-07:

#11

No "decay h > a a" should be written in the madspin stage

https://workshops.ift.uam-csic.es/files/187/MG5_Tutorial.pdf

look at slide 15 and 35

you should turn on madspin module after you do "launch" and write the decay when it asks to edit the cards

Revision history for this message

Kalash (r-b-99) said on 2024-09-08:

#12

Hello,

Thank you so much for answering and clarifying.

I still get the error at "launch" stage the after editing the card:

MG5_aMC>generate g g > h [noborn=QCD]
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=2: WEIGTHED IS QCD+2*QED
INFO: Trying coupling order WEIGHTED<=3: WEIGTHED IS QCD+2*QED
INFO: Contributing diagrams generated: 0 Born, 4 loops, 2 R2, 0 UV
1 processes with 4 diagrams generated in 0.045 s
Total: 1 processes with 4 diagrams
MG5_aMC>output gg_h_aa
INFO: initialize a new directory: gg_h_aa
INFO: remove old information in gg_h_aa
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Processing color information for loop process: g g > h [ noborn = QCD ]
INFO: Creating color matrix loop process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Creating files in directory /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/SubProcesses/PV0_0_1_gg_h
INFO: Computing diagram color coefficients
INFO: Drawing loop Feynman diagrams for Process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Creating files in directory P0_gg_h
INFO: Generating Feynman diagrams for Process: g g > h WEIGHTED<=4 [ noborn = QCD ]
INFO: Finding symmetric diagrams for subprocess group gg_h
Generated helas calls for 1 subprocesses (4 diagrams) in 0.014 s
save configuration file to /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/Cards/me5_configuration.txt
INFO: Use Fortran compiler gfortran
INFO: Use c++ compiler g++
INFO: Generate jpeg diagrams
INFO: Generate web pages
Output to directory /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa done.
Type "launch" to generate events from this process, or see
/home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/README
Run "open index.html" to see more information about this process.
MG5_aMC>launch gg_h_aa
************************************************************
* *
* W E L C O M E to *
* M A D G R A P H 5 _ a M C @ N L O *
* M A D E V E N T *
* *
* * * *
* * * * * *
* * * * * 5 * * * * *
* * * * * *
* * * *
* *
* VERSION 3.5.5 2024-06-17 *
* *
* The MadGraph5_aMC@NLO Development Team - Find us at *
* https://server06.fynu.ucl.ac.be/projects/madgraph *
* *
* Type 'help' for in-line help. *
* *
************************************************************
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/Cards/me5_configuration.txt
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/input/mg5_configuration.txt
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/Cards/me5_configuration.txt
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
generate_events run_01
The following switches determine which programs are run:
/=================== Description ===================|================== values ===================|====== other options =======\
| 1. Choose the shower/hadronization program | shower = OFF | Pythia8 |
| 2. Choose the detector simulation program | detector = OFF | Delphes |
| 3. Choose an analysis package (plot/convert) | analysis = OFF | MadAnalysis5|Rivet |
| 4. Decay onshell particles | madspin = OFF | ON|onshell|full |
| 5. Add weights to events for new hypp. | reweight = OFF | ON |
\==============================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=Pythia8' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer]
>4
The following switches determine which programs are run:
/=================== Description ===================|================== values ===================|====== other options =======\
| 1. Choose the shower/hadronization program | shower = OFF | Pythia8 |
| 2. Choose the detector simulation program | detector = OFF | Delphes |
| 3. Choose an analysis package (plot/convert) | analysis = OFF | MadAnalysis5|Rivet |
| 4. Decay onshell particles | madspin = ON | onshell|full|OFF |
| 5. Add weights to events for new hypp. | reweight = OFF | ON |
\==============================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=Pythia8' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>0
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. madspin : madspin_card.dat |
| 4. MadLoopParams : MadLoopParams.dat |
\------------------------------------------------------------/
you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
[0, done, 1, param, 2, run, 3, madspin, 4, enter path, ... ][90s to answer]
>3
open /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/Cards/madspin_card.dat
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. madspin : madspin_card.dat |
| 4. MadLoopParams : MadLoopParams.dat |
\------------------------------------------------------------/
you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
[0, done, 1, param, 2, run, 3, madspin, 4, enter path, ... ]
>0
INFO: Update the dependent parameter of the param_card.dat
WARNING: update the strong coupling value (alpha_s) to the value from the pdf selected: 0.13
WARNING: Failed to update dependent parameter. This might create trouble for external program (like MadSpin/shower/...)
Generating 10000 events with run name run_01
survey run_01
INFO: compile directory
initMadLoop -r -f
Initializing MadLoop loop-induced matrix elements (this can take some time)...
MadLoop initialization finished.
compile Source Directory
Using random number seed offset = 21
INFO: Running Survey
Creating Jobs
Working on SubProcesses
INFO: P0_gg_h
INFO: Idle: 1, Running: 1, Completed: 0 [ current time: 10h05 ]
INFO: Idle: 0, Running: 1, Completed: 1 [ 0.3s ]
INFO: P0_gg_h/G1 is at 16 +- 0.267 pb. Now submitting iteration #2.
INFO: P0_gg_h/G1 is at 15.85 +- 0.05 pb. Now submitting iteration #3.
INFO: Survey finished for P0_gg_h/G1 at 15.72 +- 0.0483 pb
INFO: Idle: 0, Running: 0, Completed: 9 [ 2.4s ]
INFO: End survey
refine 10000
Creating Jobs
INFO: Refine results to 10000
INFO: Generating 10000.0 unweighted events.
sum of cpu time of last step: 0 second
INFO: Effective Luminosity 731.5056294128872 pb^-1
INFO: need to improve 1 channels
- Current estimate of cross-section: 15.721 +- 0.048287
INFO: Idle: 1, Running: 7, Completed: 0 [ current time: 10h05 ]
INFO: P0_gg_h/G1 is at 3803/11413 (33%) event. Resubmit 6 job at iteration 1.
INFO: Idle: 0, Running: 5, Completed: 10 [ 1.1s ]
INFO: P0_gg_h/G1 is at 6679/11475 (58%) event. Resubmit 10 job at iteration 2.
INFO: found enough event for P0_gg_h/G1
INFO: Idle: 0, Running: 0, Completed: 27 [ 2.3s ]
INFO: Combining runs
sum of cpu time of last step: 0 second
INFO: finish refine
refine 10000 --treshold=0.9
No need for second refine due to stability of cross-section
INFO: Combining Events
combination of events done in 1.4699656963348389 s

=== Results Summary for run: run_01 tag: tag_1 ===

     Cumulative sequential time for this run: 17 seconds
     Cross-section : 15.69 +- 0.05053 pb
     Nb of events : 10000

fail
Failed to access python version of LHAPDF: If the python interface to LHAPDF is available on your system, try adding its location to the PYTHONPATH environment variable and theLHAPDF library location to LD_LIBRARY_PATH (linux) or DYLD_LIBRARY_PATH (mac os x).The required LD_LIBRARY_PATH is /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lhapdf6_py3//lib
INFO: can not run systematics since can not link python to lhapdf
store_events
INFO: Storing parton level results
INFO: End Parton
reweight -from_cards
decay_events -from_cards
INFO: Running MadSpin
INFO: This functionality allows for the decay of resonances
INFO: in a .lhe file, keeping track of the spin correlation effets.
INFO: BE AWARE OF THE CURRENT LIMITATIONS:
INFO: (1) Only a succession of 2 body decay are currently allowed
************************************************************
* *
* W E L C O M E to M A D S P I N *
* *
************************************************************
INFO: Extracting the banner ...
INFO: process: g g > h
INFO: options:
INFO: detected model: loop_sm. Loading...
No MadAnalysis5 version number could be read from the path supplied '/home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/madanalysis5/madanalysis5'.
The specified version of MadAnalysis5 will not be active in your session.
set fastjet to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/fastjet/bin/fastjet-config
set ninja to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lib
set collier to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lib
set lhapdf to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lhapdf6_py3/bin/lhapdf-config
set lhapdf to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lhapdf6_py3/bin/lhapdf-config
Set group_subprocesses to Auto
Note that you need to regenerate all processes
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
decay h > a a [noborn=QCD]
WARNING: Note that coupling order restriction are not associated to specific Branching Ratio. The total cross-section might therefore use the wrong branching ratio.
launch
INFO: Will use seed 788771334
INFO: We need to recalculate the branching fractions for t~,w+,z,w-,t,h
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model /home/raghad/Downloads/MG5_aMC_v3_5_4/models/loop_sm with file models/loop_sm/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
WARNING: The LO estimate for the width of particle 25
WARNING: differs from the one in the banner by 12 percent
INFO:
INFO: decay channels for t : ( width = 1.4915 GeV )
INFO: BR d1 d2
INFO: 1.000000e+00 b w+
INFO:
INFO:
INFO: decay channels for z : ( width = 2.441756 GeV )
INFO: BR d1 d2
INFO: 1.523658e-01 d~ d
INFO: 1.523658e-01 s~ s
INFO: 1.507440e-01 b~ b
INFO: 1.188141e-01 u~ u
INFO: 1.188141e-01 c~ c
INFO: 6.793744e-02 ve~ ve
INFO: 6.793744e-02 vm~ vm
INFO: 6.793744e-02 vt~ vt
INFO: 3.438713e-02 e+ e-
INFO: 3.438713e-02 mu+ mu-
INFO: 3.430972e-02 ta+ ta-
INFO:
INFO:
INFO: decay channels for w+ : ( width = 2.04793 GeV )
INFO: BR d1 d2
INFO: 3.333610e-01 d~ u
INFO: 3.333610e-01 s~ c
INFO: 1.111195e-01 e+ ve
INFO: 1.111195e-01 mu+ vm
INFO: 1.110390e-01 ta+ vt
INFO:
INFO:
INFO: decay channels for h : ( width = 0.00564955 GeV )
INFO: BR d1 d2
INFO: 9.541999e-01 b~ b
INFO: 4.580011e-02 ta+ ta-
INFO:
INFO:
INFO: decay channels for t~ : ( width = 1.4915 GeV )
INFO: BR d1 d2
INFO: 1.000000e+00 b~ w-
INFO:
INFO:
INFO: decay channels for w- : ( width = 2.04793 GeV )
INFO: BR d1 d2
INFO: 3.333610e-01 d u~
INFO: 3.333610e-01 s c~
INFO: 1.111195e-01 e- ve~
INFO: 1.111195e-01 mu- vm~
INFO: 1.110390e-01 ta- vt~
INFO:
INFO: generating the production square matrix element
INFO: generate g g > h [sqrvirt=QCD] ;
Command "generate_events run_01" interrupted with error:
InvalidCmd : not supported format standalone_msP
quit
INFO: storing files of previous run
gzipping output file: unweighted_events.lhe
INFO: Done
INFO:
more information in /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/index.html
MG5_aMC>

Hello,

Thank you so much for answering and clarifying.

I still get the error at "launch" stage the after editing the card:

MG5_aMC>generate g g > h [noborn=QCD]
INFO: Checking for minimal orders which gives processes. 
INFO: Please specify coupling orders to bypass this step. 
INFO: Trying coupling order WEIGHTED<=2: WEIGTHED IS QCD+2*QED 
INFO: Trying coupling order WEIGHTED<=3: WEIGTHED IS QCD+2*QED 
INFO: Contributing diagrams generated: 0 Born, 4 loops, 2 R2, 0 UV 
1 processes with 4 diagrams generated in 0.045 s
Total: 1 processes with 4 diagrams
MG5_aMC>output gg_h_aa
INFO: initialize a new directory: gg_h_aa 
INFO: remove old information in gg_h_aa 
INFO: Organizing processes into subprocess groups 
INFO: Generating Helas calls for process: g g > h WEIGHTED<=4 [ noborn = QCD ] 
INFO: Processing color information for loop process: g g > h [ noborn = QCD ] 
INFO: Creating color matrix  loop process: g g > h WEIGHTED<=4 [ noborn = QCD ] 
INFO: Creating files in directory /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/SubProcesses/PV0_0_1_gg_h 
INFO: Computing diagram color coefficients 
INFO: Drawing loop Feynman diagrams for Process: g g > h WEIGHTED<=4 [ noborn = QCD ] 
INFO: Creating files in directory P0_gg_h 
INFO: Generating Feynman diagrams for Process: g g > h WEIGHTED<=4 [ noborn = QCD ] 
INFO: Finding symmetric diagrams for subprocess group gg_h 
Generated helas calls for 1 subprocesses (4 diagrams) in 0.014 s
save configuration file to /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/Cards/me5_configuration.txt
INFO: Use Fortran compiler gfortran 
INFO: Use c++ compiler g++ 
INFO: Generate jpeg diagrams 
INFO: Generate web pages 
Output to directory /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa done.
Type "launch" to generate events from this process, or see
/home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/README
Run "open index.html" to see more information about this process.
MG5_aMC>launch gg_h_aa
************************************************************
*                                                          *
*                      W E L C O M E to                    *
*             M A D G R A P H 5 _ a M C @ N L O            *
*                      M A D E V E N T                     *
*                                                          *
*                 *                       *                *
*                   *        * *        *                  *
*                     * * * * 5 * * * *                    *
*                   *        * *        *                  *
*                 *                       *                *
*                                                          *
*         VERSION 3.5.5                 2024-06-17         *
*                                                          *
*    The MadGraph5_aMC@NLO Development Team - Find us at   *
*    https://server06.fynu.ucl.ac.be/projects/madgraph     *
*                                                          *
*               Type 'help' for in-line help.              *
*                                                          *
************************************************************
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/Cards/me5_configuration.txt  
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/input/mg5_configuration.txt  
INFO: load configuration from /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/Cards/me5_configuration.txt  
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
generate_events run_01
The following switches determine which programs are run:
/=================== Description ===================|================== values ===================|====== other options =======\
| 1. Choose the shower/hadronization program        |        shower = OFF                         |     Pythia8                |
| 2. Choose the detector simulation program         |      detector = OFF                         |     Delphes                |
| 3. Choose an analysis package (plot/convert)      |      analysis = OFF                         |     MadAnalysis5|Rivet     |
| 4. Decay onshell particles                        |       madspin = OFF                         |     ON|onshell|full        |
| 5. Add weights to events for new hypp.            |      reweight = OFF                         |     ON                     |
\==============================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=Pythia8' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.[60s to answer] 
>4
The following switches determine which programs are run:
/=================== Description ===================|================== values ===================|====== other options =======\
| 1. Choose the shower/hadronization program        |        shower = OFF                         |     Pythia8                |
| 2. Choose the detector simulation program         |      detector = OFF                         |     Delphes                |
| 3. Choose an analysis package (plot/convert)      |      analysis = OFF                         |     MadAnalysis5|Rivet     |
| 4. Decay onshell particles                        |       madspin = ON                          |     onshell|full|OFF       |
| 5. Add weights to events for new hypp.            |      reweight = OFF                         |     ON                     |
\==============================================================================================================================/
Either type the switch number (1 to 5) to change its setting,
Set any switch explicitly (e.g. type 'shower=Pythia8' at the prompt)
Type 'help' for the list of all valid option
Type '0', 'auto', 'done' or just press enter when you are done.
>0
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
|  1. param         : param_card.dat                         |
|  2. run           : run_card.dat                           |
|  3. madspin       : madspin_card.dat                       |
|  4. MadLoopParams : MadLoopParams.dat                      |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, madspin, 4, enter path, ... ][90s to answer] 
>3
open /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/Cards/madspin_card.dat
Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
|  1. param         : param_card.dat                         |
|  2. run           : run_card.dat                           |
|  3. madspin       : madspin_card.dat                       |
|  4. MadLoopParams : MadLoopParams.dat                      |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information on this command.
   - call an external program (ASperGE/MadWidth/...).
     Type 'help' for the list of available command
 [0, done, 1, param, 2, run, 3, madspin, 4, enter path, ... ] 
>0
INFO: Update the dependent parameter of the param_card.dat 
WARNING: update the strong coupling value (alpha_s) to the value from the pdf selected: 0.13 
WARNING: Failed to update dependent parameter. This might create trouble for external program (like MadSpin/shower/...) 
Generating 10000 events with run name run_01
survey  run_01 
INFO: compile directory 
initMadLoop -r -f
Initializing MadLoop loop-induced matrix elements (this can take some time)...
MadLoop initialization finished.
compile Source Directory
Using random number seed offset = 21
INFO: Running Survey 
Creating Jobs
Working on SubProcesses
INFO:     P0_gg_h  
INFO:  Idle: 1,  Running: 1,  Completed: 0 [ current time: 10h05 ] 
INFO:  Idle: 0,  Running: 1,  Completed: 1 [  0.3s  ] 
INFO: P0_gg_h/G1 is at 16 +- 0.267 pb. Now submitting iteration #2. 
INFO: P0_gg_h/G1 is at 15.85 +- 0.05 pb. Now submitting iteration #3. 
INFO: Survey finished for P0_gg_h/G1 at 15.72 +- 0.0483 pb 
INFO:  Idle: 0,  Running: 0,  Completed: 9 [  2.4s  ] 
INFO: End survey 
refine 10000
Creating Jobs
INFO: Refine results to 10000 
INFO: Generating 10000.0 unweighted events. 
sum of cpu time of last step: 0 second
INFO: Effective Luminosity 731.5056294128872 pb^-1 
INFO: need to improve 1 channels 
- Current estimate of cross-section: 15.721 +- 0.048287
INFO:  Idle: 1,  Running: 7,  Completed: 0 [ current time: 10h05 ] 
INFO: P0_gg_h/G1 is at 3803/11413 (33%) event. Resubmit 6 job at iteration 1. 
INFO:  Idle: 0,  Running: 5,  Completed: 10 [  1.1s  ] 
INFO: P0_gg_h/G1 is at 6679/11475 (58%) event. Resubmit 10 job at iteration 2. 
INFO: found enough event for P0_gg_h/G1 
INFO:  Idle: 0,  Running: 0,  Completed: 27 [  2.3s  ] 
INFO: Combining runs 
sum of cpu time of last step: 0 second
INFO: finish refine 
refine 10000 --treshold=0.9
No need for second refine due to stability of cross-section
INFO: Combining Events 
combination of events done in 1.4699656963348389 s 
 
 
 
  === Results Summary for run: run_01 tag: tag_1 ===

Cumulative sequential time for this run: 17 seconds
     Cross-section :   15.69 +- 0.05053 pb
     Nb of events :  10000
 
fail
Failed to access python version of LHAPDF: If the python interface to LHAPDF is available on your system, try adding its location to the PYTHONPATH environment variable and theLHAPDF library location to LD_LIBRARY_PATH (linux) or DYLD_LIBRARY_PATH (mac os x).The required LD_LIBRARY_PATH is /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lhapdf6_py3//lib
INFO: can not run systematics since can not link python to lhapdf 
store_events
INFO: Storing parton level results 
INFO: End Parton 
reweight -from_cards
decay_events -from_cards
INFO: Running MadSpin 
INFO: This functionality allows for the decay of resonances 
INFO: in a .lhe file, keeping track of the spin correlation effets. 
INFO: BE AWARE OF THE CURRENT LIMITATIONS: 
INFO:   (1) Only a succession of 2 body decay are currently allowed 
************************************************************
*                                                          *
*           W E L C O M E  to  M A D S P I N               *
*                                                          *
************************************************************
INFO: Extracting the banner ... 
INFO: process: g g > h  
INFO: options:  
INFO: detected model: loop_sm. Loading... 
No MadAnalysis5 version number could be read from the path supplied '/home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/madanalysis5/madanalysis5'.
The specified version of MadAnalysis5 will not be active in your session.
set fastjet to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/fastjet/bin/fastjet-config
set ninja to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lib
set collier to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lib
set lhapdf to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lhapdf6_py3/bin/lhapdf-config
set lhapdf to /home/raghad/Downloads/MG5_aMC_v3_5_4/HEPTools/lhapdf6_py3/bin/lhapdf-config
Set group_subprocesses to Auto
Note that you need to regenerate all processes
set max_weight_ps_point 400  # number of PS to estimate the maximum for each event
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
decay h > a a [noborn=QCD]
WARNING: Note that coupling order restriction are not associated to specific Branching Ratio. The total cross-section might therefore use the wrong branching ratio. 
launch
INFO: Will use seed 788771334 
INFO: We need to recalculate the branching fractions for t~,w+,z,w-,t,h 
INFO: Using MadWidth (arXiv:1402.1178) 
INFO: Restrict model /home/raghad/Downloads/MG5_aMC_v3_5_4/models/loop_sm with file models/loop_sm/restrict_default.dat . 
INFO: Run "set stdout_level DEBUG" before import for more information. 
WARNING: The LO estimate for the width of particle 25  
WARNING: differs from the one in the banner by 12 percent  
INFO:    
INFO: decay channels for t : ( width = 1.4915 GeV ) 
INFO:        BR                 d1  d2 
INFO:    1.000000e+00            b  w+  
INFO:    
INFO:    
INFO: decay channels for z : ( width = 2.441756 GeV ) 
INFO:        BR                 d1  d2 
INFO:    1.523658e-01            d~  d  
INFO:    1.523658e-01            s~  s  
INFO:    1.507440e-01            b~  b  
INFO:    1.188141e-01            u~  u  
INFO:    1.188141e-01            c~  c  
INFO:    6.793744e-02            ve~  ve  
INFO:    6.793744e-02            vm~  vm  
INFO:    6.793744e-02            vt~  vt  
INFO:    3.438713e-02            e+  e-  
INFO:    3.438713e-02            mu+  mu-  
INFO:    3.430972e-02            ta+  ta-  
INFO:    
INFO:    
INFO: decay channels for w+ : ( width = 2.04793 GeV ) 
INFO:        BR                 d1  d2 
INFO:    3.333610e-01            d~  u  
INFO:    3.333610e-01            s~  c  
INFO:    1.111195e-01            e+  ve  
INFO:    1.111195e-01            mu+  vm  
INFO:    1.110390e-01            ta+  vt  
INFO:    
INFO:    
INFO: decay channels for h : ( width = 0.00564955 GeV ) 
INFO:        BR                 d1  d2 
INFO:    9.541999e-01            b~  b  
INFO:    4.580011e-02            ta+  ta-  
INFO:    
INFO:    
INFO: decay channels for t~ : ( width = 1.4915 GeV ) 
INFO:        BR                 d1  d2 
INFO:    1.000000e+00            b~  w-  
INFO:    
INFO:    
INFO: decay channels for w- : ( width = 2.04793 GeV ) 
INFO:        BR                 d1  d2 
INFO:    3.333610e-01            d  u~  
INFO:    3.333610e-01            s  c~  
INFO:    1.111195e-01            e-  ve~  
INFO:    1.111195e-01            mu-  vm~  
INFO:    1.110390e-01            ta-  vt~  
INFO:    
INFO: generating the production square matrix element 
INFO: generate g g > h  [sqrvirt=QCD] ; 
Command "generate_events run_01" interrupted with error:
InvalidCmd : not supported format standalone_msP
quit
INFO: storing files of previous run 
gzipping output file: unweighted_events.lhe
INFO: Done 
INFO:  
more information in /home/raghad/Downloads/MG5_aMC_v3_5_4/gg_h_aa/index.html
MG5_aMC>

Revision history for this message

Kalash (r-b-99) said on 2024-09-08:

#13

My madspin cards is as follows:

#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#* Manual: *
#* cp3.irmp.ucl.ac.be/projects/madgraph/wiki/MadSpin *
#* *
#************************************************************
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weight 75 # number of events for the estimate of the max. weight
# set BW_cut 15 # cut on how far the particle can be off-shell
# set spinmode onshell # Use one of the madspin special mode
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event

# specify the decay for the final state particles
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
decay h > a a [noborn=QCD]

# running the actual code
launch

Many thank for your help,

Raghad.

Revision history for this message

Olivier Mattelaer (olivier-mattelaer) said on 2024-09-08:

#14

You forgot to set madspin in the
spinmode=none
mode.

Madspin does not work for loop induced decay/production otherwise.

Cheers,

Olivier

Revision history for this message

Kalash (r-b-99) said on 2024-09-09:

#15

Thank you so much.

Just a note: this is the command that worked:

set spinmode none

Many thanks for your help,

Raghad.

Revision history for this message

Kalash (r-b-99) said on 2024-09-09:

#16

Thanks Olivier Mattelaer, that solved my question.

Revision history for this message

Kalash (r-b-99) said on 2024-11-09 (last edit on 2024-11-09):

#17

Hello,

Sorry to bother you again, but the branching ratios for the diphoton and lvlv (and evmuv) branching ratios are inconsistent with the SM predictions.

The when generating ggf with lvlv decay I get the cross sections (ggf= 15.7 ± 0.033 pb, to lvlv=0.1048 ± 0.00022 pb) which results in BR of ~0.667% which is aprox half the SM prediction (1.06%).

Same thing for the ggf to diphoton, I get the cross sections (ggf=15.72 ± 0.0085 pb, to aa=0.02396 ± 1.3e-05 pb) which results in a BR of ~0.152%, which is inconsistent with the SM prediction of 0.228%.

Both madspin cards are listed here, do you know what might be causing this?

The lvlv:

#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#* Manual: *
#* cp3.irmp.ucl.ac.be/projects/madgraph/wiki/MadSpin *
#* *
#************************************************************
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weight 75 # number of events for the estimate of the max. weight
# set BW_cut 15 # cut on how far the particle can be off-shell
# set spinmode onshell # Use one of the madspin special mode
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
set spinmode none

# specify the decay for the final state particles

decay h > l+ l- vl vl~
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
# running the actual code
launch

The diphoton:

#************************************************************
#* MadSpin *
#* *
#* P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer *
#* *
#* Part of the MadGraph5_aMC@NLO Framework: *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#* Manual: *
#* cp3.irmp.ucl.ac.be/projects/madgraph/wiki/MadSpin *
#* *
#************************************************************
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weight 75 # number of events for the estimate of the max. weight
# set BW_cut 15 # cut on how far the particle can be off-shell
# set spinmode onshell # Use one of the madspin special mode
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
set spinmode none

# specify the decay for the final state particles
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
decay h > a a [noborn= QCD QED]

# running the actual code
launch

Thank you so much for your work and help.

Hello,

Sorry to bother you again, but the branching ratios for the diphoton and lvlv (and evmuv) branching ratios are inconsistent with the SM predictions.

The when generating ggf with lvlv decay I get the cross sections (ggf= 15.7 ± 0.033 pb, to lvlv=0.1048 ± 0.00022 pb) which results in BR of ~0.667% which is aprox half the SM prediction  (1.06%).

Same thing for the ggf to diphoton, I get the cross sections (ggf=15.72 ± 0.0085 pb, to aa=0.02396 ± 1.3e-05 pb) which results in a BR of ~0.152%, which is inconsistent with the SM prediction of 0.228%.

Both madspin cards are listed here, do you know what might be causing this?

The lvlv:

#************************************************************
#*                        MadSpin                           *
#*                                                          *
#*    P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer * 
#*                                                          *
#*    Part of the MadGraph5_aMC@NLO Framework:              *
#*    The MadGraph5_aMC@NLO Development Team - Find us at   *
#*    https://server06.fynu.ucl.ac.be/projects/madgraph     *
#*                                                          *
#*    Manual:                                               *
#*    cp3.irmp.ucl.ac.be/projects/madgraph/wiki/MadSpin     *
#*                                                          *
#************************************************************
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weight 75 # number of events for the estimate of the max. weight
# set BW_cut 15                 # cut on how far the particle can be off-shell
# set spinmode onshell          # Use one of the madspin special mode
 set max_weight_ps_point 400  # number of PS to estimate the maximum for each event
 set spinmode none

# specify the decay for the final state particles

decay h > l+ l- vl vl~
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
# running the actual code
launch

The diphoton:

#************************************************************
#*                        MadSpin                           *
#*                                                          *
#*    P. Artoisenet, R. Frederix, R. Rietkerk, O. Mattelaer * 
#*                                                          *
#*    Part of the MadGraph5_aMC@NLO Framework:              *
#*    The MadGraph5_aMC@NLO Development Team - Find us at   *
#*    https://server06.fynu.ucl.ac.be/projects/madgraph     *
#*                                                          *
#*    Manual:                                               *
#*    cp3.irmp.ucl.ac.be/projects/madgraph/wiki/MadSpin     *
#*                                                          *
#************************************************************
#Some options (uncomment to apply)
#
# set seed 1
# set Nevents_for_max_weight 75 # number of events for the estimate of the max. weight
# set BW_cut 15                 # cut on how far the particle can be off-shell
# set spinmode onshell          # Use one of the madspin special mode
 set max_weight_ps_point 400  # number of PS to estimate the maximum for each event
 set spinmode none

# specify the decay for the final state particles
decay t > w+ b, w+ > all all
decay t~ > w- b~, w- > all all
decay w+ > all all
decay w- > all all
decay z > all all
decay h > a a [noborn= QCD QED]

# running the actual code
launch

Thank you so much for your work and help.

Revision history for this message

Kalash (r-b-99) said on 2024-11-09:

#18

Revision history for this message

Sihyun Jeon (shjeon) said on 2024-11-09:

#19

Your BR is at LO (so there is no contributions such as h>gg). It is meant to be different from SM predictions

Revision history for this message

Kalash (r-b-99) said on 2024-11-09:

#20

Thank you so much for responding.

Sorry to bother you with this, but I need to compare this value (from MG) to an SM prediction value, so I cannot do that?

Revision history for this message

Sihyun Jeon (shjeon) said on 2024-11-09 (last edit on 2024-11-09):

#21

You of course you "can" in a technical sense. But the decay is computed at LO which makes it far from reality especially for higgs case where loops get involved.

Revision history for this message

Kalash (r-b-99) said on 2024-11-09:

#22

Thanks Sihyun Jeon, that solved my question.

MadGraph5_aMC@NLO

ggF Higgs LO and NLO in SM

Question information

This question was reopened

Subscribers

MadGraph5_aMC@NLO

ggF Higgs LO and NLO in SM

Question information

Related bugs

Related FAQ:

This question was reopened

Subscribers