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

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.

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

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

Thanks a lot,

Ilaria

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.

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

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