Error with four-fermion coupling at tree level --- Symbol ‘gc_10’ at (1) has no IMPLICIT type

Asked by Ilaria Brivio on 2017-08-29

Hi, I am testing a UFO model with four fermion operators and I get a weird error for processes in which four identical particles enter the vertex.

I am using version 2.6 and generating e+ e- > e+ e- / z h a with an operator ( l \gamma_mu l )^2. The diagrams are generated correctly but an error appears after launch or after calling madevent survey. I copy the debug.log file below.
I noticed that the coupling gc_10 does not appear in the the coup.inc file, as if it hadn't been read.

What can that be?
Thanks a lot in advance!

----------------

#************************************************************
#* MadGraph5_aMC@NLO/MadEvent *
#* *
#* * * *
#* * * * * *
#* * * * * 5 * * * * *
#* * * * * *
#* * * *
#* *
#* *
#* VERSION 2.6.0 20xx-xx-xx *
#* *
#* 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 *
#* *
#************************************************************
survey
Traceback (most recent call last):
  File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/extended_cmd.py", line 1438, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
    return func(arg, **opt)
  File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/madevent_interface.py", line 2948, in do_survey
    jobs, P_zero_result = ajobcreator.launch()
  File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/gen_ximprove.py", line 184, in launch
    self.cmd.compile(['madevent'], cwd=Pdir)
  File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/extended_cmd.py", line 1529, in compile
    return misc.compile(nb_core=self.options['nb_core'], *args, **opts)
  File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/misc.py", line 479, in compile
    raise MadGraph5Error, error_text
MadGraph5Error: A compilation Error occurs when trying to compile /home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/SubProcesses/P1_epem_epem.
The compilation fails with the following output message:
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c driver.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c myamp.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c genps.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c unwgt.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c setcuts.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c get_color.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c cuts.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c cluster.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c reweight.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c initcluster.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c addmothers.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c setscales.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c dummy_fct.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c auto_dsig.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c auto_dsig1.f -I../../Source/
    gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c matrix1.f -I../../Source/
    matrix1.f:315:55:

           CALL FFFF4_5_0(W(1,2),W(1,1),W(1,3),W(1,4),-GC_10,GC_606,AMP(1))
                                                           1
    Error: Symbol ‘gc_10’ at (1) has no IMPLICIT type
    makefile:60: set di istruzioni per l'obiettivo "matrix1.o" non riuscito
    make: *** [matrix1.o] Errore 1
    make: *** Attesa per i processi non terminati....

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 : None

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

                      Configuration Options
                      ---------------------
                text_editor : None
         cluster_local_path : None
      cluster_status_update : (600, 30)
               pythia8_path : None (user set)
                  hwpp_path : None (user set)
            pythia-pgs_path : None (user set)
                    td_path : None (user set)
               delphes_path : None (user set)
                thepeg_path : None (user set)
               cluster_type : condor
          madanalysis5_path : None (user set)
           cluster_nb_retry : 1
                 eps_viewer : None
                web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : None (user set)
                     lhapdf : lhapdf-config
              f2py_compiler : None
                 hepmc_path : None (user set)
         cluster_retry_wait : 300
           fortran_compiler : None
                auto_update : 7 (user set)
        exrootanalysis_path : None (user set)
                    timeout : 60
               cpp_compiler : None

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Ilaria Brivio
Solved:
2017-09-08
Last query:
2017-09-08
Last reply:
2017-09-07

This question was reopened

Hi,

Is the model available somewhere online?
Without it, it will be difficult to understand what the problem is.

Cheers,

Olivier

> On 29 Aug 2017, at 12:18, Ilaria Brivio <email address hidden> wrote:
>
> New question #657094 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/657094
>
> Hi, I am testing a UFO model with four fermion operators and I get a weird error for processes in which four identical particles enter the vertex.
>
> I am using version 2.6 and generating e+ e- > e+ e- / z h a with an operator ( l \gamma_mu l )^2. The diagrams are generated correctly but an error appears after launch or after calling madevent survey. I copy the debug.log file below.
> I noticed that the coupling gc_10 does not appear in the the coup.inc file, as if it hadn't been read.
>
>
> What can that be?
> Thanks a lot in advance!
>
> ----------------
>
> #************************************************************
> #* MadGraph5_aMC@NLO/MadEvent *
> #* *
> #* * * *
> #* * * * * *
> #* * * * * 5 * * * * *
> #* * * * * *
> #* * * *
> #* *
> #* *
> #* VERSION 2.6.0 20xx-xx-xx *
> #* *
> #* 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 *
> #* *
> #************************************************************
> survey
> Traceback (most recent call last):
> File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/extended_cmd.py", line 1438, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/extended_cmd.py", line 1392, in onecmd_orig
> return func(arg, **opt)
> File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/madevent_interface.py", line 2948, in do_survey
> jobs, P_zero_result = ajobcreator.launch()
> File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/gen_ximprove.py", line 184, in launch
> self.cmd.compile(['madevent'], cwd=Pdir)
> File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/extended_cmd.py", line 1529, in compile
> return misc.compile(nb_core=self.options['nb_core'], *args, **opts)
> File "/home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/bin/internal/misc.py", line 479, in compile
> raise MadGraph5Error, error_text
> MadGraph5Error: A compilation Error occurs when trying to compile /home/ilaria/Work/SMEFTsim_tests/MFV_alphaScheme/eeee_ll/SubProcesses/P1_epem_epem.
> The compilation fails with the following output message:
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c driver.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c myamp.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c genps.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c unwgt.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c setcuts.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c get_color.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c cuts.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c cluster.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c reweight.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c initcluster.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c addmothers.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c setscales.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c dummy_fct.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c auto_dsig.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c auto_dsig1.f -I../../Source/
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -w -c matrix1.f -I../../Source/
> matrix1.f:315:55:
>
> CALL FFFF4_5_0(W(1,2),W(1,1),W(1,3),W(1,4),-GC_10,GC_606,AMP(1))
> 1
> Error: Symbol ‘gc_10’ at (1) has no IMPLICIT type
> makefile:60: set di istruzioni per l'obiettivo "matrix1.o" non riuscito
> make: *** [matrix1.o] Errore 1
> make: *** Attesa per i processi non terminati....
>
> 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 : None
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> notification_center : True
> cluster_temp_path : None
> cluster_memory : None
> cluster_size : 100
> cluster_queue : None
> nb_core : 4 (user set)
> cluster_time : None
> run_mode : 2
>
> Configuration Options
> ---------------------
> text_editor : None
> cluster_local_path : None
> cluster_status_update : (600, 30)
> pythia8_path : None (user set)
> hwpp_path : None (user set)
> pythia-pgs_path : None (user set)
> td_path : None (user set)
> delphes_path : None (user set)
> thepeg_path : None (user set)
> cluster_type : condor
> madanalysis5_path : None (user set)
> cluster_nb_retry : 1
> eps_viewer : None
> web_browser : None
> syscalc_path : None (user set)
> madanalysis_path : None (user set)
> lhapdf : lhapdf-config
> f2py_compiler : None
> hepmc_path : None (user set)
> cluster_retry_wait : 300
> fortran_compiler : None
> auto_update : 7 (user set)
> exrootanalysis_path : None (user set)
> timeout : 60
> cpp_compiler : None
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Ilaria Brivio (ilariab) said : #2

Hi,

I have the UFO here:

https://www.dropbox.com/sh/e3jkx5t15gr8d7l/AADkY6-qsaZcDX525R6ObFdRa?dl=0

there are a lot of parameters. the one giving me trouble here is cll.

Thanks!
Ilaria

Hi,

Thanks this was a problem of the model restriction which can happen only for 4 fermion interaction model.
You can find the associated patch here:
http://bazaar.launchpad.net/~mg5core2/mg5amcnlo/2.6.1/diff/277?context=3

Thanks a lot for your help to make MG5aMC a bug free program,

Cheers,

Olivier

Ilaria Brivio (ilariab) said : #4

Hi Olivier, that's great, it completely solved my problem!

Thanks a lot,

Ilaria

Ilaria Brivio (ilariab) said : #5

Hi Olivier,

actually I realized that that doesn't solve *all* the problems.. the trouble is that not only the "main" coupling in a list of identical ones can appear with a minus sign in the vertices, but also the other couplings in this list.

because they appear only as -GC_XXX in the vertices, they are not appended to self.coupling_pos and they exit the loop in merge_iden_couplings() without being replaced with +- main.

Is there a simple fix?

Thanks again!
Ilaria

I would need a pratical case to understand the problem and resolve it.
Do you have one of your restriction that have such problem?

Cheers,

Olivier
> On 29 Aug 2017, at 18:14, Ilaria Brivio <email address hidden> wrote:
>
> Question #657094 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/657094
>
> Status: Solved => Open
>
> Ilaria Brivio is still having a problem:
> Hi Olivier,
>
> actually I realized that that doesn't solve *all* the problems.. the
> trouble is that not only the "main" coupling in a list of identical ones
> can appear with a minus sign in the vertices, but also the other
> couplings in this list.
>
> because they appear only as -GC_XXX in the vertices, they are not
> appended to self.coupling_pos and they exit the loop in
> merge_iden_couplings() without being replaced with +- main.
>
> Is there a simple fix?
>
> Thanks again!
> Ilaria
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Ilaria Brivio (ilariab) said : #7

yes,

 keeping the parameter cquqd10 for example gives me trouble with the coupling gc_1940, which is supposed to be the same as 1939 but doesn't get replaced because it appears with a minus sign.

Thanks again!
Ilaria

Hi Ilaria,

Sorry for the late reply, I only saw your reply today.
Not sure what happen in practise, since your reply arrive correctly to my inbox. But I missed it (not sure what happens).
I remember discussing this thread with Ken, and telling him, that I was waiting for your reply, while looks like you did it.

So i have just pushed a second patch to handle that problem:
http://bazaar.launchpad.net/~mg5core2/mg5amcnlo/2.6.1/revision/281

Cheers and sorry to have miss your email,

Olivier

Ilaria Brivio (ilariab) said : #9

Hi Olivier,

that worked like a charm, thanks a lot! No worries for the delay. ;)
I actually mentioned the thing to Fabio as well yesterday before receiving you answer.. but that is completely solved now.

Say hi to Ken and thank you very much again!

Cheers,
Ilaria