Decay width 1 -> 3 calculation crashes

[Adam Alloul]

Hi guys,

I'm trying to calculate the decay widths of some susy particles. Some of these particles only decay into 3 bodies.

After the launch command and after some time I get the following error log

calculate_decay_widths run_01 -f
Traceback (most recent call last):
  File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/extended_cmd.py", line 819, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/extended_cmd.py", line 812, in onecmd_orig
    return func(arg, **opt)
  File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/madevent_interface.py", line 2315, in do_calculate_decay_widths
    postcmd=False)
  File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/extended_cmd.py", line 859, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/extended_cmd.py", line 812, in onecmd_orig
    return func(arg, **opt)
  File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/madevent_interface.py", line 2584, in do_survey
    misc.compile(['gensym'], cwd=Pdir)
  File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/various/misc.py", line 223, in compile
    raise MadGraph5Error, error_text
MadGraph5Error: A compilation Error occurs when trying to compile /opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/decay_widths_III/SubProcesses/P0_x2p_x
2ptaumtaup.
The compilation fails with the following output message:
    gfortran -O -w -ffixed-line-length-132 -w -c -o symmetry.o symmetry.f
    touch qmass.inc
    gfortran -O -w -ffixed-line-length-132 -w -c -o cuts.o cuts.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o cluster.o cluster.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o myamp.o myamp.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o genps.o genps.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o initcluster.o initcluster.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o setscales.o setscales.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o reweight.o reweight.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o get_color.o get_color.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o idenparts.o idenparts.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o matrix.o matrix.f
    gfortran -O -w -ffixed-line-length-132 -w -c -o setcuts.o setcuts.f
    gfortran -O -w -ffixed-line-length-132 -w -o gensym symmetry.o setcuts.o cuts.o cluster.o myamp.o genps.o initcluster.o setscales.o reweight.o get_color.o
 idenparts.o matrix.o -L../../lib/ -ldhelas -ldsample -lmodel -lgeneric -lpdf -lcernlib
    matrix.o: In function `matrix_':
    matrix.f:(.text+0x1bc): undefined reference to `ffv1c1p0_3_'
    collect2: ld returned 1 exit status
    make: *** [gensym] Error 1

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 : 20 (user set)

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
          cluster_temp_path : None
              cluster_queue : None (user set)
                    nb_core : 16 (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 : None (user set)
                    td_path : None (user set)
               delphes_path : None (user set)
                auto_update : 7 (user set)
               cluster_type : condor
           fortran_compiler : None (user set)
        exrootanalysis_path : None (user set)
                 eps_viewer : None
                    timeout : 60

Do you have any idea from where this may come from?

My model has been generated by FeynRules through the UFO interface and I have set all the unknown decay widths to 0 and I'm using the version 1.5.8 of MadGraph5 on Scientific Linux SL release 5.5 (Boron)

Cheers,
Adam.

[Pierre Artoisenet (partois)]

Dear Adam,

Thanks for your bug report.
Olivier Mattelaer is away for a week and he aksed me to inform
the mg5 users that he will answer bug reports on next Sunday.

Thanks for your patience,

Pierre

[Olivier Mattelaer (olivier-mattelaer)]

Hi Adam,

Could you send me by email or attach to this bug report (whichever you prefer)
1) your model
2) your processes.

I think that this is easy to fix, but I need those to be able to reproduce your error and be 100% sure that this is correctly fixed.

Cheers,

Olivier

On Apr 18, 2013, at 12:39 PM, Adam Alloul <email address hidden> wrote:

> Public bug reported:
>
> Hi guys,
>
> I'm trying to calculate the decay widths of some susy particles. Some of
> these particles only decay into 3 bodies.
>
> After the launch command and after some time I get the following error
> log
>
> calculate_decay_widths run_01 -f
> Traceback (most recent call last):
> File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/extended_cmd.py", line 819, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/extended_cmd.py", line 812, in onecmd_orig
> return func(arg, **opt)
> File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/madevent_interface.py", line 2315, in do_calculate_decay_widths
> postcmd=False)
> File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/extended_cmd.py", line 859, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/extended_cmd.py", line 812, in onecmd_orig
> return func(arg, **opt)
> File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/interface/madevent_interface.py", line 2584, in do_survey
> misc.compile(['gensym'], cwd=Pdir)
> File "/opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/madgraph/various/misc.py", line 223, in compile
> raise MadGraph5Error, error_text
> MadGraph5Error: A compilation Error occurs when trying to compile /opt/sbg/grand-est/ui7_data1/aalloul/MadGraph5_v1_5_8/decay_widths_III/SubProcesses/P0_x2p_x
> 2ptaumtaup.
> The compilation fails with the following output message:
> gfortran -O -w -ffixed-line-length-132 -w -c -o symmetry.o symmetry.f
> touch qmass.inc
> gfortran -O -w -ffixed-line-length-132 -w -c -o cuts.o cuts.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o cluster.o cluster.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o myamp.o myamp.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o genps.o genps.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o initcluster.o initcluster.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o setscales.o setscales.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o reweight.o reweight.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o get_color.o get_color.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o idenparts.o idenparts.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o matrix.o matrix.f
> gfortran -O -w -ffixed-line-length-132 -w -c -o setcuts.o setcuts.f
> gfortran -O -w -ffixed-line-length-132 -w -o gensym symmetry.o setcuts.o cuts.o cluster.o myamp.o genps.o initcluster.o setscales.o reweight.o get_color.o
> idenparts.o matrix.o -L../../lib/ -ldhelas -ldsample -lmodel -lgeneric -lpdf -lcer...

[Olivier Mattelaer (olivier-mattelaer)]

Hi Adam,

Thanks so much to have found this problem.
I've fixed it in 1.5.10 (the next version of MG5).
In the mean time, you can use the patch in attachment in order to have a working code.

Cheers,

Olivier

[Olivier Mattelaer (olivier-mattelaer)]

In order to use this patch,
you have to do
1) copy it in the MG5 directory
2) patch -p1 c1p0.patch

Cheers,

Olivier

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

I've update the patch since I was not 100% sure that the previous one was correct in all cases.

Cheers,

Olivier

[Adam Alloul (adam-alloul)]

Hi Olivier,

thank you for the update. I'll run the tests right now and keep you updated.

Cheers,

Adam.

[Adam Alloul (adam-alloul)]

Hi Olivier,

unfortunately I have a bad news. The patch you gave me doesn't seem to solve the whole problem. I still get exactly the same problem at the same place.

On the other hand, all the other decay processes seem to work so I'm wondering if the problem is from my model? (I can send it to you if you ever want and have time to look at it?)

Cheers,

Adam.

[Olivier Mattelaer (olivier-mattelaer)]

Hi Adam,

Are you sure that the patch works correctly?

If you look at th file UpdatesNotes do you have the following line:
1.5.10 (xx/05/13) OM: Fix Bug #1170417: fix crash for conjugate routine in presence of
                      massless propagator. (introduce in 1.5.9)
                  OM: Fix question #226810: checking that patch program exists before
                      trying to update MG5 code.
                  OM: Fix Bug #1171049: an error in the order of wavefuntions
                      making the code to crash (introduce in 1.5.7)

If not this means that the patch fails to be done.
Looks like that my instruction for the patch was wrong: this is the correct syntax.
patch -p1 < c1p0.patch

Cheers,

Olivier

[Adam Alloul (adam-alloul)]

Thank you for the prompt answer!

Indeed, the syntax was incorrect but I managed to find the correct one. I do have patched correctly the files. I checked one by one even the UpdateNotes.txt file and they agree with the patch file you sent me. Here are the first lines of UpdateNotes.txt that I have

1.5.10 (xx/05/13) OM: Fix Bug #1170417: fix crash for conjugate routine in presence of
                      massless propagator. (introduce in 1.5.9)
                  OM: Fix question #226810: checking that patch program exists before
                      trying to update MG5 code.
                  OM: Fix Bug #1171049: an error in the order of wavefuntions
                      making the code to crash (introduce in 1.5.7)

Cheers,

Adam.

[Adam Alloul (adam-alloul)]

Hi Olivier.

I just wanted to let you know that I somehow solved the problem. I noticed that the crash was due to a process of the form

p1 > p1 tau- tau+

where p1 is some massive particle. All I had to do was to remove p1 from the definition of "all" and everything went smooth.

Thank you very much for your help!

Cheers,
Adam.

[Adam Alloul (adam-alloul)]

Of course I removed it because such process is supposed to give a zero xsection