Problem generating events with third party model

Hi,

I am trying to generate dihiggs events with an external model from:
https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/HiggsPairProduction

I do not see model on that page, only prepared LO process.
Can you give me a direct html link to your model?

However, when I do "ouput pp2hh" I gives me an error :
"""
Command "output pp2hh" interrupted with error:
KeyError : (0,)
"""

Did you ask the model author, if your process is supported?
Since this model contains a lot of hacking, I would not be surprise that it can ONLY do
one process (i.e. t t~ h h) and potentially even MG5aMC version specific.

This issue is actually my first guess.

Could someone please help me out?

I would actually also need the content of the MG5_debug in order to (try to) help you.

Cheers,

Olivier

On 17 Oct 2018, at 11:52, Peter Meiring <<email address hidden><mailto:<email address hidden>>> wrote:

New question #675282 on MadGraph5_aMC@NLO:
https://answers.launchpad.net/mg5amcnlo/+question/675282

Hi MadGraph experts!

I am trying to generate dihiggs events with an external model from:
https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/HiggsPairProduction

I placed the "tthh" model in MG5_aMC_v2.6.3.2/models. Importing the model gives me no problems.
Also the command "generate p p > h h" works fine and I can see the produced diagrams with "display diagrams".

However, when I do "ouput pp2hh" I gives me an error :
"""
Command "output pp2hh" interrupted with error:
KeyError : (0,)
"""

Here is the whole log:
"""
[meiringp@dhcp-135-247 MG5_aMC_v2.6.3.2]$ ./bin/mg5_aMC
************************************************************
* *
* 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.6.3.2 2018-06-22 *
* *
* 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
fastjet-config does not seem to correspond to a valid fastjet-config executable (v3+). We will use fjcore instead.
Please set the 'fastjet'variable to the full (absolute) /PATH/TO/fastjet-config (including fastjet-config).
MG5_aMC> set fastjet /PATH/TO/fastjet-config

set lhapdf to /home/meiringp/LHAPDF6/bin/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
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>import model ./models/tthhUFO/
INFO: Change particles name to pass to MG5 convention
Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchangedDefined 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 p p > h h
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=4: WEIGTHED IS 2*QED+QCD
INFO: Trying coupling order WEIGHTED<=5: WEIGTHED IS 2*QED+QCD
INFO: Trying coupling order WEIGHTED<=6: WEIGTHED IS 2*QED+QCD
INFO: Trying process: g g > h h WEIGHTED<=6 @1
INFO: Process has 2 diagrams
INFO: Trying process: u u~ > h h WEIGHTED<=6 @1
INFO: Trying process: u c~ > h h WEIGHTED<=6 @1
INFO: Trying process: c u~ > h h WEIGHTED<=6 @1
INFO: Trying process: c c~ > h h WEIGHTED<=6 @1
INFO: Trying process: d d~ > h h WEIGHTED<=6 @1
INFO: Trying process: d s~ > h h WEIGHTED<=6 @1
INFO: Trying process: s d~ > h h WEIGHTED<=6 @1
INFO: Trying process: s s~ > h h WEIGHTED<=6 @1
1 processes with 2 diagrams generated in 0.087 s
Total: 1 processes with 2 diagrams

MG5_aMC>display diagrams
Drawing Process: g g > h h WEIGHTED<=6 @1
Wrote file /tmp/diagrams_1_gg_hh.eps
open /tmp/diagrams_1_gg_hh.eps
time to draw 0.00582408905029

MG5_aMC>output pp2hh
INFO: directory /home/meiringp/MG5_aMC_v2.6.3.2/pp2hh already exists.
If you continue this directory will be deleted and replaced.
Do you want to continue? [y, n][60s to answer]
y
INFO: initialize a new directory: pp2hh
INFO: remove old information in pp2hh
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: g g > h h WEIGHTED<=6 @1
INFO: Processing color information for process: g g > h h @1
INFO: Creating files in directory P1_gg_hh
INFO: Generating Feynman diagrams for Process: g g > h h WEIGHTED<=6 @1
INFO: Finding symmetric diagrams for subprocess group gg_hh
Generated helas calls for 1 subprocesses (2 diagrams) in 0.005 s
Wrote files for 7 helas calls in 0.039 s
ALOHA: aloha creates SSS1 routines
ALOHA: aloha creates VVSS5 routines
Command "output pp2hh" interrupted with error:
KeyError : (0,)
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.
"""

I also tried to apply the fixes on https://answers.launchpad.net/mg5amcnlo/+question/670817.
However, for me these didn't work.
Could someone please help me out?

Thanks in advance,

Peter

--
You received this question notification because you are an answer
contact for MadGraph5_aMC@NLO.

Hi Olivier,

Thanks for the quick response!

> Can you give me a direct html link to your model?
The model can be downloaded with:
https://cp3.irmp.ucl.ac.be/projects/madgraph/raw-attachment/wiki/HiggsPairProduction/tthh.tar.gz

> Did you ask the model author, if your process is supported?
Yes, I contacted one of the authors, but have received no response. Indeed it should only contain the LO process, but that is fine for now; I am interfacing MG with Pythia8 inside the MG framework to do the hadronization, but everything is on a very introductory level right now.

> I would actually also need the content of the MG5_debug in order to (try to) help you.
Actually, I found that I had a typo in applying the fixes from https://answers.launchpad.net/mg5amcnlo/+question/670817.
A new error occurs "Error: Symbol ‘mt’ at (1) has no IMPLICIT type". This seems like a clear error statement, but since this is the first time I see Fortran code, I am not sure if I should explicitly set the type (and if so, what type; also should it be set in the model scripts or in the generated process directory?).
Here follows the debug log "run_01_tag_1_debug.log":

"""
#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.3.2 2018-06-22 *
#* *
#* 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/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1501, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/madevent_interface.py", line 2469, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/common_run_interface.py", line 6711, in new_fct
    original_fct(obj, *args, **opts)
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/madevent_interface.py", line 2508, in run_generate_events
    postcmd=False)
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
    return func(arg, **opt)
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/madevent_interface.py", line 3277, in do_survey
    self.configure_directory()
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/madevent_interface.py", line 5603, in configure_directory
    self.compile(arg=[name], cwd=os.path.join(self.me_dir, 'Source'))
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1592, in compile
    return misc.compile(nb_core=self.options['nb_core'], *args, **opts)
  File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/various/misc.py", line 527, in compile
    raise MadGraph5Error, error_text
MadGraph5Error: A compilation Error occurs when trying to compile /home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source.
The compilation fails with the following output message:
    ar cru ../lib/libdsample.a dsample.o ranmar.o DiscreteSampler.o StringCast.o
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o setrun.o setrun.f
    cd DHELAS; make
    cd PDF; make
    cd MODEL; make
    ranlib ../lib/libdsample.a
    make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
    make[1]: Entering directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/DHELAS'
    make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
    make[1]: Entering directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/PDF'
    make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
    make[1]: Entering directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/MODEL'
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o VVSS2_0.o VVSS2_0.f
    ar cru ../../lib/libpdf.a Ctq6Pdf.o pdfwrap.o opendata.o pdf.o PhotonFlux.o pdg2pdf.o NNPDFDriver.o
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o rw_para.o rw_para.f
    ranlib ../../lib/libpdf.a
    make[1]: Leaving directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/PDF'
    VVSS2_0.f:101:12:

           mQ2=MT**2
                1
    Error: Symbol ‘mt’ at (1) has no IMPLICIT type
    make[1]: *** [<builtin>: VVSS2_0.o] Error 1
    make[1]: Leaving directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/DHELAS'
    make: *** [makefile:44: ../lib/libdhelas.a] Error 2
    make: *** Waiting for unfinished jobs....
    ar cru ../../lib/libmodel.a couplings.o lha_read.o printout.o rw_para.o model_functions.o couplings1.o couplings2.o
    ranlib ../../lib/libmodel.a
    make[1]: Leaving directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/MODEL'

Please try to fix this compilations issue and retry.
Help might be found at https://answers.launchpad.net/mg5amcnlo.
If you think that this is a bug, you can report this at https://bugs.launchpad.net/mg5amcnlo
                              Run Options
                              -----------
               stdout_level : 20 (user set)

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
        notification_center : True
          cluster_temp_path : None
             cluster_memory : None (user set)
               cluster_size : 100
              cluster_queue : None (user set)
                    nb_core : 8 (user set)
               cluster_time : 8 (user set)
                   run_mode : 2

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : /home/meiringp/MG5_aMC_v2.6.3.2/HEPTools/pythia8 (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : /home/meiringp/Delphes-3.4.1 (user set)
                thepeg_path : None (user set)
               cluster_type : condor
          madanalysis5_path : None (user set)
           cluster_nb_retry : 1
                 eps_viewer : None
                web_browser : /usr/bin/firefox (user set)
               syscalc_path : None (user set)
           madanalysis_path : None (user set)
                     lhapdf : /home/meiringp/LHAPDF6/bin/lhapdf-config (user set)
              f2py_compiler : None
                 hepmc_path : None (user set)
         cluster_retry_wait : 300
           fortran_compiler : None
                auto_update : 7 (user set)
        exrootanalysis_path : /home/meiringp/MG5_aMC_v2.6.3.2/ExRootAnalysis (user set)
                    timeout : 60
               cpp_compiler : None
"""

Cheers,

Peter

Hi,

>> Did you ask the model author, if your process is supported?
> Yes, I contacted one of the authors, but have received no response. Indeed it should only contain the LO process, but that is fine for now; I am interfacing MG with Pythia8 inside the MG framework to do the hadronization, but everything is on a very introductory level right now.

Note that one of the author left physics. So you should contact Eleni (she is at CERN for the moment, so use her cern email). She has a launchpad account so I will try to assign her to this question.

> A new error occurs "Error: Symbol ‘mt’ at (1) has no IMPLICIT type".

This means that the model is clearly outdated. It will need to be rewritten by changing
MT in all fortran file by mdl_MT (potentially the same for all other parameter of the model).
The best would obviously be that the model author updates their tarball.
Otherwise you should be able to do it yourself.

Cheers,

Olivier

> On 17 Oct 2018, at 13:32, Peter Meiring <email address hidden> wrote:
>
> Question #675282 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/675282
>
> Status: Answered => Open
>
> Peter Meiring is still having a problem:
> Hi Olivier,
>
> Thanks for the quick response!
>
>> Can you give me a direct html link to your model?
> The model can be downloaded with:
> https://cp3.irmp.ucl.ac.be/projects/madgraph/raw-attachment/wiki/HiggsPairProduction/tthh.tar.gz
>
>> Did you ask the model author, if your process is supported?
> Yes, I contacted one of the authors, but have received no response. Indeed it should only contain the LO process, but that is fine for now; I am interfacing MG with Pythia8 inside the MG framework to do the hadronization, but everything is on a very introductory level right now.
>
>> I would actually also need the content of the MG5_debug in order to (try to) help you.
> Actually, I found that I had a typo in applying the fixes from https://answers.launchpad.net/mg5amcnlo/+question/670817.
> A new error occurs "Error: Symbol ‘mt’ at (1) has no IMPLICIT type". This seems like a clear error statement, but since this is the first time I see Fortran code, I am not sure if I should explicitly set the type (and if so, what type; also should it be set in the model scripts or in the generated process directory?).
> Here follows the debug log "run_01_tag_1_debug.log":
>
> """
> #************************************************************
> #* MadGraph5_aMC@NLO/MadEvent *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.3.2 2018-06-22 *
> #* *
> #* 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/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1501, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
> return func(arg, **opt)
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/madevent_interface.py", line 2469, in do_generate_events
> self.run_generate_events(switch_mode, args)
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/common_run_interface.py", line 6711, in new_fct
> original_fct(obj, *args, **opts)
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/madevent_interface.py", line 2508, in run_generate_events
> postcmd=False)
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1528, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1450, in onecmd_orig
> return func(arg, **opt)
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/madevent_interface.py", line 3277, in do_survey
> self.configure_directory()
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/madevent_interface.py", line 5603, in configure_directory
> self.compile(arg=[name], cwd=os.path.join(self.me_dir, 'Source'))
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/interface/extended_cmd.py", line 1592, in compile
> return misc.compile(nb_core=self.options['nb_core'], *args, **opts)
> File "/home/meiringp/MG5_aMC_v2.6.3.2/madgraph/various/misc.py", line 527, in compile
> raise MadGraph5Error, error_text
> MadGraph5Error: A compilation Error occurs when trying to compile /home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source.
> The compilation fails with the following output message:
> ar cru ../lib/libdsample.a dsample.o ranmar.o DiscreteSampler.o StringCast.o
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o setrun.o setrun.f
> cd DHELAS; make
> cd PDF; make
> cd MODEL; make
> ranlib ../lib/libdsample.a
> make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
> make[1]: Entering directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/DHELAS'
> make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
> make[1]: Entering directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/PDF'
> make[1]: warning: jobserver unavailable: using -j1. Add '+' to parent make rule.
> make[1]: Entering directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/MODEL'
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o VVSS2_0.o VVSS2_0.f
> ar cru ../../lib/libpdf.a Ctq6Pdf.o pdfwrap.o opendata.o pdf.o PhotonFlux.o pdg2pdf.o NNPDFDriver.o
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o rw_para.o rw_para.f
> ranlib ../../lib/libpdf.a
> make[1]: Leaving directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/PDF'
> VVSS2_0.f:101:12:
>
> mQ2=MT**2
> 1
> Error: Symbol ‘mt’ at (1) has no IMPLICIT type
> make[1]: *** [<builtin>: VVSS2_0.o] Error 1
> make[1]: Leaving directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/DHELAS'
> make: *** [makefile:44: ../lib/libdhelas.a] Error 2
> make: *** Waiting for unfinished jobs....
> ar cru ../../lib/libmodel.a couplings.o lha_read.o printout.o rw_para.o model_functions.o couplings1.o couplings2.o
> ranlib ../../lib/libmodel.a
> make[1]: Leaving directory '/home/meiringp/MG5_aMC_v2.6.3.2/pp2hh_tthh2/Source/MODEL'
>
> Please try to fix this compilations issue and retry.
> Help might be found at https://answers.launchpad.net/mg5amcnlo.
> If you think that this is a bug, you can report this at https://bugs.launchpad.net/mg5amcnlo
> Run Options
> -----------
> stdout_level : 20 (user set)
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> notification_center : True
> cluster_temp_path : None
> cluster_memory : None (user set)
> cluster_size : 100
> cluster_queue : None (user set)
> nb_core : 8 (user set)
> cluster_time : 8 (user set)
> run_mode : 2
>
> Configuration Options
> ---------------------
> text_editor : None
> cluster_local_path : None
> cluster_status_update : (600, 30)
> pythia8_path : /home/meiringp/MG5_aMC_v2.6.3.2/HEPTools/pythia8 (user set)
> hwpp_path : None (user set)
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : /home/meiringp/Delphes-3.4.1 (user set)
> thepeg_path : None (user set)
> cluster_type : condor
> madanalysis5_path : None (user set)
> cluster_nb_retry : 1
> eps_viewer : None
> web_browser : /usr/bin/firefox (user set)
> syscalc_path : None (user set)
> madanalysis_path : None (user set)
> lhapdf : /home/meiringp/LHAPDF6/bin/lhapdf-config (user set)
> f2py_compiler : None
> hepmc_path : None (user set)
> cluster_retry_wait : 300
> fortran_compiler : None
> auto_update : 7 (user set)
> exrootanalysis_path : /home/meiringp/MG5_aMC_v2.6.3.2/ExRootAnalysis (user set)
> timeout : 60
> cpp_compiler : None
> """
>
> Cheers,
>
> Peter
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Hi,

we do not provide any UFO model on the website:
https://cp3.irmp.ucl.ac.be/projects/madgraph/wiki/HiggsPairProduction#no1

We only give generated processes folders, which would work only with very old versions of madgraph, from the time before loop-induced production was possible in MG5.

If you need to generate Standard Model events for hh, just use the loop_sm model (default model in mg5):

generate p p > h h [QCD]

Cheers,

Eleni