Error when using "Auto" for decay width in param_card.dat

Asked by Pawin Ittisamai

Hello,

When I set the decay width to "Auto"; for example,
DECAY 6 Auto # WT
I got an error message that includes
"""
                Included at param_read.inc:5:
                Included at rw_para.f:26:

                  WT = AUTO
                           1
            Error: Symbol 'auto' at (1) has no IMPLICIT type
            make[1]: *** [rw_para.o] Error 1
"""
[The complete log, from version 1.5.7, will be posted below this message].

This kind of errors occurs in the default settings/processes in the Template directory, in both 1.3.30 and 1.5.7. MadGraph runs fine if I don't use that "Auto" width. Do you know what happened and how to fix this?

Thank you very much and have a good day!
Pawin

[Log file]

#************************************************************
#* MadGraph/MadEvent 5 *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.1.5.7 *
#* *
#* The MadGraph Development Team - Please visit us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events -f --multicore default-mg5-processes-autowidth --nb_core=16
Traceback (most recent call last):
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/bin/internal/extended_cmd.py", line 819, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/bin/internal/extended_cmd.py", line 812, in onecmd_orig
    return func(arg, **opt)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/bin/internal/madevent_interface.py", line 2114, in do_generate_events
    postcmd=False)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/bin/internal/extended_cmd.py", line 859, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/bin/internal/extended_cmd.py", line 812, in onecmd_orig
    return func(arg, **opt)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/bin/internal/madevent_interface.py", line 2431, in do_survey
    self.configure_directory()
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/bin/internal/madevent_interface.py", line 3708, in configure_directory
    misc.compile(arg=[name], cwd=os.path.join(self.me_dir, 'Source'))
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/bin/internal/misc.py", line 223, in compile
    raise MadGraph5Error, error_text
MadGraph5Error: A compilation Error occurs when trying to compile /mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/Source.
The compilation fails with the following output message:
    ar cru ../lib/libdsample.a dsample.o ranmar.o
    /usr/bin/gfortran -I/opt/software/ClusterStudio/2011.0/composerxe-2011.5.220/mkl/include/intel64/lp64 -I/opt/software/ClusterStudio/2011.0/composerxe-2011.5.220/mkl/include/fftw -I/opt/software/ClusterStudio/2011.0/composerxe-2011.5.220/mkl/include -ffixed-line-length-132 -c -o setrun.o setrun.f
    ranlib ../lib/libdsample.a
    cd MODEL; make
    ar cru ../lib/libgeneric.a alfas_functions.o transpole.o invarients.o hfill.o pawgraphs.o ran1.o rw_events.o rw_routines.o kin_functions.o open_file.o basecode.o setrun.o run_printout.o dgauss.o readgrid.o getissud.o
    make[1]: Entering directory `/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/Source/MODEL'
    make[1]: warning: jobserver unavailable: using -j1. Add `+' to parent make rule.
    /usr/bin/gfortran -I/opt/software/ClusterStudio/2011.0/composerxe-2011.5.220/mkl/include/intel64/lp64 -I/opt/software/ClusterStudio/2011.0/composerxe-2011.5.220/mkl/include/fftw -I/opt/software/ClusterStudio/2011.0/composerxe-2011.5.220/mkl/include -ffixed-line-length-132 -ffixed-line-length-132 -c -o rw_para.o rw_para.f
    ../param_card.inc:14.15:
        Included at param_read.inc:5:
        Included at rw_para.f:26:

          WT = AUTO
                   1
    Error: Symbol 'auto' at (1) has no IMPLICIT type
    make[1]: *** [rw_para.o] Error 1
    make[1]: Leaving directory `/mnt/home/ittisama/Work/MadGraph5_v1_5_7/TestMG157/Source/MODEL'
    make: *** [../lib/libmodel.a] Error 2
    make: *** Waiting for unfinished jobs....
    ranlib ../lib/libgeneric.a

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

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
          cluster_temp_path : None
              cluster_queue : madgraph
                    nb_core : 32 (user set)
                   run_mode : 2

                      Configuration Options
                      ---------------------
                web_browser : None
                text_editor : None
           madanalysis_path : None (user set)
               pythia8_path : None (user set)
            pythia-pgs_path : /mnt/home/ittisama/Work/MadGraph5_v1_5_7/pythia-pgs (user set)
                    td_path : None (user set)
               delphes_path : None (user set)
                auto_update : 7 (user set)
               cluster_type : pbs (user set)
           fortran_compiler : gfortran (user set)
        exrootanalysis_path : None (user set)
                 eps_viewer : None
                    timeout : 60

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#1

Thanks so much for reporting this.
I've created a bug report linked to this behavior.

A temporary fix is to avoid of the '-f' options.
Which bypass the step to compute the width on the fly.

If you don't use this options directly, this means that you way you submit your job is equivalent to ask for the -f.

one way to submit which is equivalent but going to work:
./bin/generate_events default-mg5-processes-autowidth
then you will have two question: just answer by pressing enter.

If you need to script this:
write a file with the following line:
generate_events default-mg5-processes-autowidth
and run it like:
./bin/madevent FILE

Note that by default the script is using multi-core (and all available).

Finally, concerning 1.3.X, this version indeed didn't support the Auto syntax (which was introduced in 1.5.x)

Cheers,

Olivier

Revision history for this message
Pawin Ittisamai (ittisama) said :
#2

Thanks Olivier Mattelaer, that solved my question.

Revision history for this message
Pawin Ittisamai (ittisama) said :
#3

Olivier, I just wanted to thank you again (in addition to the automatic message generated above).

A quick question that's probably related to this: the option "Auto" does not work with models generated in "v4" format, right?

When I tried to do so, I got the following messages:

"""
>
Computing the width set on auto in the param_card.dat
Be carefull automatic computation of the width is
ONLY valid if all three (or more) body decay are negligeable. In doubt use a
calculator.
In a future version of MG5 those mode will also be taken into account
Command "generate_events TEST-axi-pp-C-jj-lhc14-gl1gr1-autowidth-cteq5m-defaultcuts-cutT" interrupted with error:
UFOImportError : Path /mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj is not a valid pathname
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in '/mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj-cteq5m/TEST-axi-pp-C-jj-lhc14-gl1gr1-cteq5m-defaultcuts-cutT_tag_1_debug.log'.
Please attach this file to your report.
quit
"""

And here is the log file
"""
#************************************************************
#* MadGraph/MadEvent 5 *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 5.1.5.7 *
#* *
#* The MadGraph Development Team - Please visit us at *
#* https://server06.fynu.ucl.ac.be/projects/madgraph *
#* *
#************************************************************
#* *
#* Command File for MadEvent *
#* *
#* run as ./bin/madevent.py filename *
#* *
#************************************************************
generate_events TEST-axi-pp-C-jj-lhc14-gl1gr1-autowidth-cteq5m-defaultcuts-cutT
Traceback (most recent call last):
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj-cteq5m/bin/internal/extended_cmd.py", line 819, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj-cteq5m/bin/internal/extended_cmd.py", line 812, in onecmd_orig
    return func(arg, **opt)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj-cteq5m/bin/internal/madevent_interface.py", line 2085, in do_generate_events
    self.ask_run_configuration(mode)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj-cteq5m/bin/internal/madevent_interface.py", line 4217, in ask_run_configuration
    self.check_param_card(pjoin(self.me_dir,'Cards','param_card.dat' ))
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj-cteq5m/bin/internal/madevent_interface.py", line 4495, in check_param_card
    self.do_compute_widths('%s %s' % (' '.join(pdg), path))
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj-cteq5m/bin/internal/madevent_interface.py", line 2680, in do_compute_widths
    args = self.check_compute_widths(args)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj-cteq5m/bin/internal/madevent_interface.py", line 906, in check_compute_widths
    decay=True)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/models/import_ufo.py", line 80, in import_model
    model_path = find_ufo_path(model_name)
  File "/mnt/home/ittisama/Work/MadGraph5_v1_5_7/models/import_ufo.py", line 65, in find_ufo_path
    raise UFOImportError("Path %s is not a valid pathname" % model_name)
UFOImportError: Path /mnt/home/ittisama/Work/MadGraph5_v1_5_7/work-axigluon-pp-C-jj is not a valid pathname
                              Run Options
                              -----------
               stdout_level : None

                         MadEvent Options
                         ----------------
     automatic_html_opening : True
          cluster_temp_path : None
              cluster_queue : madgraph
                    nb_core : 32 (user set)
                   run_mode : 2

                      Configuration Options
                      ---------------------
                web_browser : None
                text_editor : None
           madanalysis_path : None (user set)
               pythia8_path : None (user set)
            pythia-pgs_path : /mnt/home/ittisama/Work/MadGraph5_v1_5_7/pythia-pgs (user set)
                    td_path : None (user set)
               delphes_path : None (user set)
                auto_update : 7 (user set)
               cluster_type : pbs (user set)
           fortran_compiler : gfortran (user set)
        exrootanalysis_path : None (user set)
                 eps_viewer : None
                    timeout : 60
"""

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

>Olivier, I just wanted to thank you again (in addition to the automatic message generated above).

You are welcome.

No you right,

> A quick question that's probably related to this: the option "Auto" does not work with models generated in "v4" format, right?

This is only working for some specific UFO model (not even all of them).
In order to generate a UFO model supporting this syntax, you need to use the current development version of FeynRules.
This version is not yet public but you can have it if you ask any of the FeynRules author to have it.

In fact FeynRules provides us the analytical formula for 1>2 body decay. (so this method is not safe for the Higgs width where 1>3 are not negligeable). We have in development a version of MG5 where FR will compute 1>2 and MG5 will automatically estimate if the 1>3 are negligeable or not and if not will use madevent in order to have the precise value of the 1>3 decay. (This will probably part of MG5 1.6.0)

Cheers,

Olivier

Revision history for this message
Pawin Ittisamai (ittisama) said :
#5

Thanks again : )

Pawin