Trying to generate the process g g > h h , h > b b~ , h > l- l+ l- l+ with little success.

Hi,

> the whole process should be done in NLO

I hope that you confuse LO and NLO here since Higgs pair production is NOT know at NLO precision.
Nobody is able to do such process at NLO precision since nobody is able to compute the two-loop diagram needed in that case.
Now, since this process do not have any leading order diagram the leading order is at one loop. So a lot of people call it a NLO contribution by mistake.

So in the latest version of MG5, you can indeed compute the LO contribution for process starting at one loop (also call loop-induced processes).
If you want approximate NLO computation, you can refer to arXiv:1401.7340 and some of the material can be found here: https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/HiggsPairProduction

> MadGraph5Error : Failed the initialisation of loop-induced matrix element 'PV0_0_1_gg_hh' (trying with a maximum of 15 PS points).

Do you have a correct compiler? What is your compiler?
Do you succeed to compute cross-section for p p > t t~ [QCD] ?

> And if I try to specify the decay modes I get an error immediately.

Decay chain syntax are not supported for loop-induced processes (neither at NLO actually).
The only solution (which is actually perfect for the Higgs) is to use MadSpin (which can be run on the flight) with the following line in the madspin_card:

set spinmode=none
decay h > b b~
decay h > l+ l- l+ l-

Cheers,

Olivier

On 31 Jul 2015, at 14:36, Gerardo <email address hidden> wrote:

> New question #269871 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/269871
>
> Hello Madgraph team,
>
> I am attempting to generate a gluon gluon fusion to higgs pair process for a study on the FCC (Future Circular Collider). I also require that the one higgs decays into b bbar while the other decays into 4 leptons; the whole process should be done in NLO. I am semi-successful in generating the gluon gluon fusion to higgs pair portion of the process with the commands,
>
> generate g g > h h [QCD]
> output FCC_gg_hh
> launch FCC_gg_hh
>
> however, after the launch command I get an error. It looks like:
>
> ************************************************************
> * *
> * W E L C O M E to *
> * M A D G R A P H 5 _ a M C @ N L O *
> * *
> * *
> * * * *
> * * * * * *
> * * * * * 5 * * * * *
> * * * * * *
> * * * *
> * *
> * VERSION 2.3.0 2015-07-01 *
> * *
> * The MadGraph5_aMC@NLO Development Team - Find us at *
> * https://server06.fynu.ucl.ac.be/projects/madgraph *
> * and *
> * http://amcatnlo.web.cern.ch/amcatnlo/ *
> * *
> * Type 'help' for in-line help. *
> * Type 'tutorial' to learn how MG5 works *
> * Type 'tutorial aMCatNLO' to learn how aMC@NLO works *
> * Type 'tutorial MadLoop' to learn how MadLoop works *
> * *
> ************************************************************
> load MG5 configuration from input/mg5_configuration.txt
> set fastjet to /home/jerry/fastjet-3.1.3/fastjet-config
> set lhapdf to lhapdf-config
> Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
> Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
> Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
> Loading default 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
> Defined multiparticle p = g u c d s u~ c~ d~ s~
> Defined multiparticle j = g u c d s u~ c~ d~ s~
> Defined multiparticle l+ = e+ mu+
> Defined multiparticle l- = e- mu-
> Defined multiparticle vl = ve vm vt
> Defined multiparticle vl~ = ve~ vm~ vt~
> 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 w- ta- ta+
> MG5_aMC>generate g g > h h [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 l- / j / vl / l+ / p / 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+
> ------------------------------------------------------------------------
> No Born diagrams found. Now switching to the loop-induced mode.
> Please cite ref. 'arXiv:1507.00020' when using results from this mode.
> ------------------------------------------------------------------------
> 1 processes with 16 diagrams generated in 1.212 s
> Total: 1 processes with 16 diagrams
> MG5_aMC>output FCC_gg_hh
> INFO: initialize a new directory: FCC_gg_hh
> INFO: remove old information in FCC_gg_hh
> INFO: Organizing processes into subprocess groups
> INFO: Generating Helas calls for process: g g > h h WEIGHTED=8 [ noborn = QCD ]
> INFO: Processing color information for loop process: g g > h h [ noborn = QCD ]
> INFO: Creating color matrix loop process: g g > h h WEIGHTED=8 [ noborn = QCD ]
> INFO: Creating files in directory /home/jerry/MG5_aMC_v2_3_0/FCC_gg_hh/SubProcesses/PV0_0_1_gg_hh
> ALOHA: aloha creates FFV1 set of routines with options: L1
> INFO: Computing diagram color coefficients
> INFO: Drawing loop Feynman diagrams for Process: g g > h h WEIGHTED=8 [ noborn = QCD ]
> INFO: Creating files in directory P0_gg_hh
> INFO: Generating Feynman diagrams for Process: g g > h h WEIGHTED=8 [ noborn = QCD ]
> INFO: Finding symmetric diagrams for subprocess group gg_hh
> History written to /home/jerry/MG5_aMC_v2_3_0/FCC_gg_hh/Cards/proc_card_mg5.dat
> Generated helas calls for 1 subprocesses (16 diagrams) in 0.366 s
> Export UFO model to MG4 format
> ALOHA: aloha creates FFS1 set of routines with options: L2
> ALOHA: aloha creates VVS1 routines
> ALOHA: aloha creates SSS1 routines
> ALOHA: aloha creates R2_GGHH routines
> save configuration file to /home/jerry/MG5_aMC_v2_3_0/FCC_gg_hh/Cards/me5_configuration.txt
> INFO: Use Fortran compiler gfortran
> INFO: Generate jpeg diagrams
> INFO: Generate web pages
> Output to directory /home/jerry/MG5_aMC_v2_3_0/FCC_gg_hh done.
> Type "launch" to generate events from this process, or see
> /home/jerry/MG5_aMC_v2_3_0/FCC_gg_hh/README
> Run "open index.html" to see more information about this process.
> MG5_aMC>launch FCC_gg_hh
> ************************************************************
> * *
> * 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 2.3.0 2015-07-01 *
> * *
> * 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/jerry/MG5_aMC_v2_3_0/FCC_gg_hh/Cards/me5_configuration.txt
> INFO: load configuration from /home/jerry/MG5_aMC_v2_3_0/input/mg5_configuration.txt
> INFO: load configuration from /home/jerry/MG5_aMC_v2_3_0/FCC_gg_hh/Cards/me5_configuration.txt
> Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
> No valid SysCalc path found
> generate_events run_01
> The following switches determine which programs are run:
> 1 Run the pythia shower/hadronization: pythia=OFF
> 2 Run PGS as detector simulator: pgs=OFF
> 3 Run Delphes as detector simulator: delphes=OFF
> 4 Decay particles with the MadSpin module: madspin=OFF
> 5 Add weight to events based on coupling parameters: reweight=OFF
> Either type the switch number (1 to 5) to change its default setting,
> or set any switch explicitly (e.g. type 'madspin=ON' at the prompt)
> Type '0', 'auto', 'done' or just press enter when you are done.
> [0, 1, 2, 3, 4, 5, auto, done, pythia=ON, ... ][60s to answer]
>>
> Do you want to edit a card (press enter to bypass editing)?
> 1 / param : param_card.dat
> 2 / run : run_card.dat
> 3 / MadLoopParams : MadLoopParams.dat
> 9 / plot : plot_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, MadLoopParams, 9, enter path, ... ][60s to answer]
>>
> Generating 10000 events with run name run_01
> survey run_01
> INFO: compile directory
> initMadLoop -r -f
>
> (process:4612): GLib-CRITICAL **: g_slice_set_config: assertion 'sys_page_size == 0' failed
> Initializing MadLoop loop-induced matrix elements (this can take some time)...
> root: Error while executing ./check in /home/jerry/MG5_aMC_v2_3_0/FCC_gg_hh/SubProcesses/PV0_0_1_gg_hh
> root: Failed at running the process in /home/jerry/MG5_aMC_v2_3_0/FCC_gg_hh/SubProcesses/PV0_0_1_gg_hh.
> Error detected in "generate_events run_01"
> write debug file /home/jerry/MG5_aMC_v2_3_0/FCC_gg_hh/run_01_tag_1_debug.log
> If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
> MadGraph5Error : Failed the initialization of loop-induced matrix element 'PV0_0_1_gg_hh' (trying with a maximum of 15 PS points).
> quit
> INFO:
>
> INFO:
>
> ***************************************************************************************************************
>
> And if I try to specify the decay modes I get an error immediately. The commands and errors are:
>
> MG5_aMC>define ha = h
> Defined multiparticle ha = h
> MG5_aMC>define hb = h
> Defined multiparticle hb = h
> MG5_aMC>generate g g > ha hb , ha > b b~ , hb > l- l+ l- l+ [QCD]
> Error detected in "generate g g > ha hb , ha > b b~ , hb > l- l+ l- l+ [QCD]"
> write debug file MG5_debug
> If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
> MadGraph5Error : Decay processes cannot be perturbed
> MG5_aMC>
>
> ****************************************************************************************************************
>
> What are my options, if any, in generating these events using Madgraph?
>
> Much appreciated,
> Gerardo
>
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Please discard my last message, which arrive on the wrong thread.

Olivier

Hi,

>Now, since this process do not have any leading order diagram the leading order is at one loop. So a lot of people call it a NLO >contribution by mistake.

I am guilty of the same mistake; I only need the one loop process, so in that case I need the LO contribution.

>Do you have a correct compiler? What is your compiler?
>Do you succeed to compute cross-section for p p > t t~ [QCD] ?

I have the compiler:
gcc version 4.9.2 (Ubuntu 4.9.2-10Ubuntu13)
No, g g > t t~ also fails. I suppose I have the incorrect compiler then. What compiler is needed?

Gerardo

Hi,

 > I suppose I have the incorrect compiler then. What compiler is needed?

I guess yes,
gcc 4.9.2 should be working fine.

What is the error that you get for p p > t t~ [QCD] it might be more clear what the problem is in that case.

Cheers,

Olivier

Hi,

>What is the error that you get for p p > t t~ [QCD] it might be more clear what the problem is in that case.

Well its not exactly and error but, I get a few warnings and then it remains Idle for a very long time. I only ran 100 events but its been Idle for 30+ minutes now.

This is the what is says after launch:

INFO: ...done, continuing with P* directories
INFO: Compiling directories...
INFO: Compiling on 4 cores
INFO: Compiling P0_gg_ttx...
INFO: Compiling P0_uux_ttx...
INFO: Compiling P0_uxu_ttx...
./check_poles: error while loading shared libraries: libgolem.so.0: cannot open shared object file: No such file or directory
./check_poles: error while loading shared libraries: libgolem.so.0: cannot open shared object file: No such file or directory
./check_poles: error while loading shared libraries: libgolem.so.0: cannot open shared object file: No such file or directory
INFO: P0_uxu_ttx done.
INFO: P0_uux_ttx done.
INFO: P0_gg_ttx done.
INFO: Checking test output:
INFO: P0_gg_ttx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_uux_ttx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: P0_uxu_ttx
INFO: Result for test_ME:
INFO: Passed.
INFO: Result for test_MC:
INFO: Passed.
INFO: Result for check_poles:
WARNING: 0 points have been tried
INFO: Starting run
INFO: Using 4 cores
INFO: Cleaning previous results
INFO: Doing NLO matched to parton shower
INFO: Setting up grid
WARNING: program /home/jerry/MG5_aMC_v2_3_0/compiler_test/SubProcesses/P0_uxu_ttx/ajob1 2 F 0 launch ends with non zero status: 127. Stop all computation
INFO: Idle: 1, Running: 0, Completed: 4 [ current time: 08h33 ]

I ran the generation on a non-local computer and the same process ran successfully and in a minute or so. I will try your suggestion about the madspin_card non-localy as it seems to compile correctly and let you know how that goes.

Gerardo

Hi,

I have already seen this golem stuff and Valentin knows how to fix it.
Let's wait that he appears online.

Cheers,

Olivier

Hi,

I ran the generation non-locally and it seemed to work but gave me a warning that Golem was not configured so would not be used. I installed Golem and then re-launched and now it gives me this error:

Generating 100 events with run name run_01
survey run_01
INFO: compile directory

(process:2633): GLib-CRITICAL **: g_slice_set_config: assertion 'sys_page_size == 0' failed
initMadLoop -r -f
Initializing MadLoop loop-induced matrix elements (this can take some time)...
root: Error while executing ./check in /home/jerry/MG5_aMC_v2_3_0/FCC_Test/SubProcesses/PV0_0_1_gg_hh
root: Failed at running the process in /home/jerry/MG5_aMC_v2_3_0/FCC_Test/SubProcesses/PV0_0_1_gg_hh.
Error detected in "generate_events run_01"
write debug file /home/jerry/MG5_aMC_v2_3_0/FCC_Test/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
MadGraph5Error : Failed the initialization of loop-induced matrix element 'PV0_0_1_gg_hh' (trying with a maximum of 15 PS points).

HI,

golem is not mandatory for the run, so just remove the link to golem if this creates trouble.

Cheers,

Olivier
On 03 Aug 2015, at 12:21, Gerardo <email address hidden> wrote:

> Question #269871 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/269871
>
> Status: Answered => Open
>
> Gerardo is still having a problem:
> Hi,
>
> I ran the generation non-locally and it seemed to work but gave me a
> warning that Golem was not configured so would not be used. I installed
> Golem and then re-launched and now it gives me this error:
>
>
> Generating 100 events with run name run_01
> survey run_01
> INFO: compile directory
>
> (process:2633): GLib-CRITICAL **: g_slice_set_config: assertion 'sys_page_size == 0' failed
> initMadLoop -r -f
> Initializing MadLoop loop-induced matrix elements (this can take some time)...
> root: Error while executing ./check in /home/jerry/MG5_aMC_v2_3_0/FCC_Test/SubProcesses/PV0_0_1_gg_hh
> root: Failed at running the process in /home/jerry/MG5_aMC_v2_3_0/FCC_Test/SubProcesses/PV0_0_1_gg_hh.
> Error detected in "generate_events run_01"
> write debug file /home/jerry/MG5_aMC_v2_3_0/FCC_Test/run_01_tag_1_debug.log
> If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
> MadGraph5Error : Failed the initialization of loop-induced matrix element 'PV0_0_1_gg_hh' (trying with a maximum of 15 PS points).
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi,

Okay, I got rid of the Golem path and it fixed the last error. Now it gets to Madspin and decays the particles but it gives me and error and doesn't finish:

WARNING: Order of particle in the event did not agree with parent/child order. This might be problematic for some code.
Command "generate_events run_01" interrupted with error:
IndexError : list index out of range
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in '/home/jerry/MG5_aMC_v2_3_0/FCC_Test/run_01_tag_1_debug.log'.
Please attach this file to your report.

The debug log is:

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.3.0 2015-07-01 *
#* *
#* The MadGraph5_aMC@NLO Development Team - Find us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events run_01
Traceback (most recent call last):
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/interface/extended_cmd.py", line 879, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/interface/madevent_interface.py", line 2019, in do_generate_events
    self.exec_cmd('decay_events -from_cards', postcmd=False)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/interface/extended_cmd.py", line 919, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/interface/common_run_interface.py", line 1828, in do_decay_events
    madspin_cmd.import_command_file(path)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/interface/extended_cmd.py", line 1038, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/interface/extended_cmd.py", line 919, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/interface/extended_cmd.py", line 872, in onecmd_orig
    return func(arg, **opt)
  File "/home/jerry/MG5_aMC_v2_3_0/MadSpin/interface_madspin.py", line 454, in do_launch
    return self.run_bridge(line)
  File "/home/jerry/MG5_aMC_v2_3_0/MadSpin/interface_madspin.py", line 868, in run_bridge
    particle.add_decay(decay)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/various/lhe_parser.py", line 141, in add_decay
    return self.event.add_decay_to_particle(self.event_id, decay_event)
  File "/home/jerry/MG5_aMC_v2_3_0/madgraph/various/lhe_parser.py", line 901, in add_decay_to_particle
    setattr(new_particle, tag, self[nb_part + mother_id -1])
IndexError: list index out of range

Sorry for having so many problems,

Gerardo

Hi,

The patched worked. You have been a great help Olivier, thank you.

cheers,
Gerardo