Problem with spin-2 particles

Asked by Eric

Hi,

I am trying to generate events for the production of spin-2 particles with MG5_aMC 2.9.9. I am using the DMspin2 UFO model from https://feynrules.irmp.ucl.ac.be/wiki/Spin2#no1, but I run into the same problem with the RS model that is automatically downloaded when trying to import it in MG. After letting MG convert the UFO to Python3, I run
generate p p > Y2
output test
launch
I then get the following error:

Command "generate_events run_01" interrupted with error:
ValueError : (1, -1) is not in list

And in the log file:
Traceback (most recent call last):
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/extended_cmd.py", line 1541, in onecmd
    return self.onecmd_orig(line, **opt)
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/madevent_interface.py", line 2474, in do_generate_events
    self.run_generate_events(switch_mode, args)
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/common_run_interface.py", line 7111, in new_fct
    original_fct(obj, *args, **opts)
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/madevent_interface.py", line 2513, in run_generate_events
    self.exec_cmd('survey %s %s' % (self.run_name,' '.join(args)),
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
    return func(arg, **opt)
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/madevent_interface.py", line 3371, in do_survey
    jobs, P_zero_result = ajobcreator.get_helicity()
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/gen_ximprove.py", line 288, in get_helicity
    recycler.generate_output_file()
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 738, in generate_output_file
    self.read_orig()
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 694, in read_orig
    self.template_dict['helas_calls'] += self.unfold_helicities(
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 588, in unfold_helicities
    new_objs = Amplitude.generate_amps(line, self.dag, self.all_hel, self.bad_amps_perhel)
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 359, in generate_amps
    new_amps = [cls.get_obj(line, wavs, graph, diag_num)
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 361, in <listcomp>
    if cls.good_helicity(wavs, graph, diag_num, all_hel,all_bad_hel)]
  File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 150, in good_helicity
    hel_number = 1 + all_hel.index(tuple(this_hel))
ValueError: (1, -1) is not in list

Is this a bug in MadGraph, a problem in the model, or something that went wrong in the translation from py2 to py3?

Thanks for your help.

Best,
Eric

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said :
#1

Hi,

You can put hel_recycling to False in the run_card.dat and this should fix the problem.
I will investigate why this is not done automatically (it should)

Cheers,

Olivier

> On 20 Apr 2022, at 19:05, Eric <email address hidden> wrote:
>
> New question #701439 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/701439
>
> Hi,
>
> I am trying to generate events for the production of spin-2 particles with MG5_aMC 2.9.9. I am using the DMspin2 UFO model from https://feynrules.irmp.ucl.ac.be/wiki/Spin2#no1, but I run into the same problem with the RS model that is automatically downloaded when trying to import it in MG. After letting MG convert the UFO to Python3, I run
> generate p p > Y2
> output test
> launch
> I then get the following error:
>
> Command "generate_events run_01" interrupted with error:
> ValueError : (1, -1) is not in list
>
> And in the log file:
> Traceback (most recent call last):
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/extended_cmd.py", line 1541, in onecmd
> return self.onecmd_orig(line, **opt)
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
> return func(arg, **opt)
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/madevent_interface.py", line 2474, in do_generate_events
> self.run_generate_events(switch_mode, args)
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/common_run_interface.py", line 7111, in new_fct
> original_fct(obj, *args, **opts)
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/madevent_interface.py", line 2513, in run_generate_events
> self.exec_cmd('survey %s %s' % (self.run_name,' '.join(args)),
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/extended_cmd.py", line 1570, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/extended_cmd.py", line 1490, in onecmd_orig
> return func(arg, **opt)
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/interface/madevent_interface.py", line 3371, in do_survey
> jobs, P_zero_result = ajobcreator.get_helicity()
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/gen_ximprove.py", line 288, in get_helicity
> recycler.generate_output_file()
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 738, in generate_output_file
> self.read_orig()
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 694, in read_orig
> self.template_dict['helas_calls'] += self.unfold_helicities(
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 588, in unfold_helicities
> new_objs = Amplitude.generate_amps(line, self.dag, self.all_hel, self.bad_amps_perhel)
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 359, in generate_amps
> new_amps = [cls.get_obj(line, wavs, graph, diag_num)
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 361, in <listcomp>
> if cls.good_helicity(wavs, graph, diag_num, all_hel,all_bad_hel)]
> File "~/MadGraph/MG5_aMC_v2_9_9/madgraph/madevent/hel_recycle.py", line 150, in good_helicity
> hel_number = 1 + all_hel.index(tuple(this_hel))
> ValueError: (1, -1) is not in list
>
> Is this a bug in MadGraph, a problem in the model, or something that went wrong in the translation from py2 to py3?
>
> Thanks for your help.
>
> Best,
> Eric
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Eric (e-r-i-c) said :
#2

Thanks Olivier Mattelaer, that solved my question.

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

just for information:
this patch ensures to automatically remove the problematic optimization for spin 3/2 and spin2:

https://bazaar.launchpad.net/~maddevelopers/mg5amcnlo/LTS_dev/revision/404

Thanks,

Olivier