compilation error topBSM model

Dear Mohammad, (cc Claude)

The solution is to edit the file:
functions_library.py and change the following line from:

cot = Function(name = 'cot',
               arguments = ('z',),
               expression = '1./cmath.tan(z)')

to
cot = Function(name = 'cot',
               arguments = ('z',),
               expression = '1./cmath.tan(z.real)')

I thought that such change was include in the latest version of FeynRules already.
Claude can you check? (the error reported is actually compiler dependent)

Cheers,

Olivier

On 15 Apr 2015, at 10:11, Mohammad <email address hidden> wrote:

> New question #265236 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/265236
>
> Hello,
> I want to use MG5_aMC for calculations of some processes using the topBSM model ( http://feynrules.irmp.ucl.ac.be/wiki/topBSM ) . I used the latest version of FeynRules to create UFO model files which I imported in MG5_aMC.
> ---------------------------------------------------------------------------------------
> $FeynRulesPath = "/theo/altakach/Mohammad/mathematica/feynrules"; \
> SetDirectory[$FeynRulesPath];
> << FeynRules`;
>
> SetDirectory[$FeynRulesPath <> "/Models/topBSM"];
> LoadModel["SM.fr", "topBSM.fr"];
> FeynmanRules[LSM, LS0, LO0, LS1, LO1, ScreenOutput -> False];
> WriteUFO[LSM + LS0 + LO0 + LS1 + LO1];
> ------------------------------------------------------------------------
> Generating the process works fine and I can also view the Feynman diagram nicely.
> But when I want to launch the calculation, I get a compilation error.
> --------------------------------------------------------------
> MG5_aMC>import model topBSM_UFO
> MG5_aMC>generate p p > t t~
> MG5_aMC>output realtest
> MG5_aMC>launch realtest
> -----------------------------------------------------------
> Here is what I get:
> ************************************************************
> * *
> * 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.2.3 2015-02-10 *
> * *
> * 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 /theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/Cards/me5_configuration.txt
> INFO: load configuration from /theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/input/mg5_configuration.txt
> INFO: load configuration from /theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/Cards/me5_configuration.txt
> Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
> No valid eps viewer found. Please set in ./input/mg5_configuration.txt
> 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=NOT INSTALLED
> 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, 4, 5, auto, done, pythia=ON, pythia=OFF, ... ][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
> 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, enter path][60s to answer]
>>
> Generating 10000 events with run name run_01
> survey run_01
> INFO: compile directory
> Traceback (most recent call last):
> File "write_param_card.py", line 205, in <module>
> ParamCardWriter('./param_card.dat', generic=True)
> File "write_param_card.py", line 30, in __init__
> self.write_card(list_of_parameters)
> File "write_param_card.py", line 90, in write_card
> self.write_dep_param_block(lhablock)
> File "write_param_card.py", line 123, in write_dep_param_block
> exec("%s = %s" % (parameter.name, parameter.value))
> File "<string>", line 1, in <module>
> ZeroDivisionError: complex division
>
>
> Error detected in "generate_events run_01"
> write debug file /theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/run_01_tag_1_debug.log
> If you need help with this issue please contact us on https://answers.launchpad.net/madgraph5
> MadGraph5Error : A compilation Error occurs when trying to compile /theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/Source.
> The compilation fails with the following output message:
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o dsample.o dsample.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o alfas_functions.o alfas_functions.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o transpole.o transpole.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o invarients.o invarients.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o hfill.o hfill.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o pawgraphs.o pawgraphs.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o ran1.o ran1.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o rw_events.o rw_events.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o kin_functions.o kin_functions.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o basecode.o basecode.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o setrun.o setrun.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o run_printout.o run_printout.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o dgauss.o dgauss.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o readgrid.o readgrid.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o getissud.o getissud.f
> cd DHELAS; make
> cd PDF; make
> make[1]: Entering directory `/theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/Source/DHELAS'
> make[1]: warning: jobserver unavailable: using -j1. Add `+' to parent make rule.
> make[1]: Entering directory `/theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/Source/PDF'
> make[1]: warning: jobserver unavailable: using -j1. Add `+' to parent make rule.
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o Ctq4Fn.o Ctq4Fn.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o aloha_functions.o aloha_functions.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o Ctq5Par.o Ctq5Par.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o Ctq5Pdf.o Ctq5Pdf.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o Partonx5.o Partonx5.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o Ctq6Pdf.o Ctq6Pdf.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o cteq3.o cteq3.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o mrs98.o mrs98.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV1_0.o FFV1_0.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV1_2.o FFV1_2.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV1_1.o FFV1_1.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o VVV1P0_1.o VVV1P0_1.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV1P0_3.o FFV1P0_3.f
> ar cru ../../lib/libdhelas.a aloha_functions.o FFV1_0.o FFV1_2.o FFV1_1.o VVV1P0_1.o FFV1P0_3.o
> ranlib ../../lib/libdhelas.a
> make[1]: Leaving directory `/theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/Source/DHELAS'
> cd MODEL; make
> make[1]: Entering directory `/theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/Source/MODEL'
> make[1]: warning: jobserver unavailable: using -j1. Add `+' to parent make rule.
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o couplings.o couplings.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o lha_read.o lha_read.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o mrs98lo.o mrs98lo.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o printout.o printout.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o rw_para.o rw_para.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o model_functions.o model_functions.f
> model_functions.f:51.29:
>
> COT = 1.000000D+00/TAN(Z)
> 1
> Error: 'x' argument of 'tan' intrinsic at (1) must be REAL
> make[1]: *** [model_functions.o] Error 1
> make[1]: Leaving directory `/theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/Source/MODEL'
> make: *** [../lib/libmodel.a] Error 2
> make: *** Waiting for unfinished jobs....
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o mrs98ht.o mrs98ht.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o mrs99.o mrs99.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o mrst2001.o mrst2001.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o mrst2002.o mrst2002.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o jeppe02.o jeppe02.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o pdfwrap.o pdfwrap.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o opendata.o opendata.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o pdf.o pdf.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o PhotonFlux.o PhotonFlux.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o pdg2pdf.o pdg2pdf.f
> gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o NNPDFDriver.o NNPDFDriver.f
> ar cru ../../lib/libpdf.a Ctq4Fn.o Ctq5Par.o Ctq5Pdf.o Partonx5.o Ctq6Pdf.o cteq3.o mrs98.o mrs98lo.o mrs98ht.o mrs99.o mrst2001.o mrst2002.o jeppe02.o pdfwrap.o opendata.o pdf.o PhotonFlux.o pdg2pdf.o NNPDFDriver.o
> ranlib ../../lib/libpdf.a
> make[1]: Leaving directory `/theo/altakach/Documents/madgraph5/MG5_aMC_v2_2_3/bin/realtest/Source/PDF'
>
> 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
>
>
>
>
>
>
> Thank you for your help.
>
> Cheers.
>
>
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.