MadGraph5_aMC@NLO

gfortan compliation error with gldgrv UFO

Asked by Joakim Flinckman

Hi,
Im currently trying to look at a specific event ( i.e g g > t1 t1~ ,(t1 > t n1, t > w+ b , n1 > b b~ grv) , (t1~ > t~ n1 , t~ > w- b~ , n1 > b~ b grv) ) with the MSSM goldstino/gravitino model from http://feynrules.irmp.ucl.ac.be/wiki/goldstino. I get the UFO to work and I can replicate the plots described in the article verifying the model (Simulating spin-3/2 particles at colliders, arXiv:1308.1668) and many other types of events works fine.

Both the generating and output runs smoothly, however, when launching the event above I get the following (very long) error message:

INFO: compile directory
compile Source Directory
Error detected in "generate_events run_01"
write debug file /Users/joakimflinckman/MG5_aMC_v2_6_0/slask/run_01_tag_1_debug.log
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
MadGraph5Error : A compilation Error occurs when trying to compile /Users/joakimflinckman/MG5_aMC_v2_6_0/slask/Source.
 The compilation fails with the following output message:
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o StringCast.o StringCast.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o ranmar.o ranmar.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 rw_routines.o rw_routines.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 open_file.o open_file.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
     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 FRS7_3.o FRS7_3.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFS12_3.o FFS12_3.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FRS6_3.o FRS6_3.f
     cd PDF; make
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FRS5_3.o FRS5_3.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FRS1_3.o FRS1_3.f
     cd MODEL; make
     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 couplings.o couplings.f
     cd CERNLIB; make
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFV2_3.o FFV2_3.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o lha_read.o lha_read.f
     cd BIAS/dummy; make
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o abend.o abend.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 FRS2_3.o FRS2_3.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o dummy.o dummy.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o dlsqp2.o dlsqp2.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o Ctq5Pdf.o Ctq5Pdf.f
     ar cru ../../../lib/libbias.a dummy.o
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o lenocc.o lenocc.f
     ranlib ../../../lib/libbias.a
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FFS1C1_2.o FFS1C1_2.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 mtlprt.o mtlprt.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o DiscreteSampler.o DiscreteSampler.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 mtlset.o mtlset.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 model_functions.o model_functions.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o radmul.o radmul.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FRS3_3.o FRS3_3.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o Ctq6Pdf.o Ctq6Pdf.f
     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
     ranlib ../lib/libgeneric.a
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o FRS2C1_2.o FRS2C1_2.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o cteq3.o cteq3.f
     FRS2C1_2.f:63:8:

              R1(I) = R1(I) + RTMP(I)
             1
     Error: Unclassifiable statement at (1)
     FRS2C1_2.f:67:8:

              R1(I) = R1(I) + RTMP(I)
             1
     Error: Unclassifiable statement at (1)
     FRS2C1_2.f:71:8:

              R1(I) = R1(I) + RTMP(I)
             1
     Error: Unclassifiable statement at (1)
     FRS2C1_2.f:60:40:

            CALL FRS2C1_2(R2,S3,COUP1,M1,W1,R1)
                                             1
     Error: Symbol ‘r1’ at (1) has no IMPLICIT type
     make[1]: *** [FRS2C1_2.o] Error 1
     make[1]: *** Waiting for unfinished jobs....
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o couplings1.o couplings1.f
     ar cru ../../lib/libcernlib.a abend.o dlsqp2.o lenocc.o mtlprt.o mtlset.o radmul.o
     ranlib ../../lib/libcernlib.a
     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 couplings2.o couplings2.f
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o mrs98lo.o mrs98lo.f
     make: *** [../lib/libdhelas.a] Error 2
     make: *** Waiting for unfinished jobs....
     gfortran -O -w -fbounds-check -fPIC -ffixed-line-length-132 -c -o couplings3.o couplings3.f
     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 couplings4.o couplings4.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
     ar cru ../../lib/libmodel.a couplings.o lha_read.o printout.o rw_para.o model_functions.o couplings1.o couplings2.o couplings3.o couplings4.o
     ranlib ../../lib/libmodel.a

 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

Any idea of how I could solve this?

Cheers,
Joakim

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Joakim Flinckman
Solved:
Last query:
Last reply:
Revision history for this message
Joakim Flinckman (jofl90) said : 		#1

Not sure if this is relevant but, I have also tried to run each "sub event" and they all run smoothly, however there seams to be a problem in the t1 > t n1 , ( n1 > b b~ grv ) part, since trying to run this results in the same error message as above.

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

Hi,

Thanks for the information,

Here is a patch that we are testing right now and which should solve your issue

=== modified file 'aloha/aloha_writers.py'
--- aloha/aloha_writers.py 2017-09-27 22:16:08 +0000
+++ aloha/aloha_writers.py 2018-04-11 13:07:23 +0000
@@ -862,7 +862,7 @@

         lor_list = (self.routine.name,) + lor_names
         line = " call %(name)s%(addon)s(%(before_coup)s,%(coup)s%(after_coup)s,%(out)s)\n"
- main = '%(spin)s%(id)d' % {'spin': self.particles[self.offshell -1],
+ main = '%(spin)s%(id)d' % {'spin': self.particles[self.outgoing -1],
                            'id': self.outgoing}
         for i, name in enumerate(lor_list):
             data['name'] = name
@@ -885,7 +885,7 @@
                 if not offshell:
                     routine.write( ' vertex = vertex + tmp\n')
                 else:
- size = self.type_to_size[self.particles[offshell -1]] -2
+ size = self.type_to_size[self.particles[self.outgoing -1]] -2
                     routine.write(" do i = %s, %s\n" % (self.momentum_size+1, self.momentum_size+size))
                     routine.write(" %(main)s(i) = %(main)s(i) + %(tmp)s(i)\n" %\
                                {'main': main, 'tmp': data['out']})

> On 10 Apr 2018, at 10:37, Joakim Flinckman <email address hidden> wrote:
>
> Question #667793 on MadGraph5_aMC@NLO changed:
> https://answers.launchpad.net/mg5amcnlo/+question/667793
>
> Joakim Flinckman gave more information on the question:
> Not sure if this is relevant but, I have also tried to run each "sub
> event" and they all run smoothly, however there seams to be a problem in
> the t1 > t n1 , ( n1 > b b~ grv ) part, since trying to run this results
> in the same error message as above.
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

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

Actually,

additionally you need to do this:
=== modified file 'aloha/aloha_writers.py'
--- aloha/aloha_writers.py 2018-04-11 13:12:48 +0000
+++ aloha/aloha_writers.py 2018-04-11 15:25:50 +0000
@@ -874,7 +874,7 @@
                     data['out'] = main
             elif i==1:
                 if self.offshell:
- type = self.particles[self.offshell-1]
+ type = self.particles[self.outgoing-1]
                     self.declaration.add(('list_complex','%stmp' % type))
                 else:
                     type = ''

Cheers and thanks
(thanks to Peter Richardson for his test on this actually)

Revision history for this message
Joakim Flinckman (jofl90) said : 		#4

This seams to work beautifully!
Will this potentially interfere with other processes I run, or was this just a Bug in the code?

Many Cheers
/Joakim

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

Hi,

No this should not interfere with anything else.

Cheers,

Olivier