Error in re-weighting between different groups of operators in EFT model

Asked by Alexandros Marantis on 2021-02-21

Hello Olivier,
I, systematically, get this error when I'm trying to re-weight between two different operators in this EFT model (http://feynrules.irmp.ucl.ac.be/wiki/AnomalousGaugeCoupling#no1), and more specifically, when I'm trying to re-weight between two different "groups" (T0,T1,T2 -> T8,T9 e.g re-weight a sample generated with non-zero T1 to non-zero T8).
There is a KeyError at the end of the debug file. (i'm attaching it)
Can you please help me to understand what "KeyError: ((-1, 2), (-13, -11, -3, 4, 11, 13))" means?
In fact, what I'd like to know is if this error refers to a bug of the program or the inability of the model to re-weight between different groups of operators.

Thank you in advance,
Alexandros

P.S The re-weighting process works successfully and with great precision, between the operators of the same group.

//---------------------------- from ME5_debug ----------------------------//
reweight run_01_0
Traceback (most recent call last):
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/_GENERATION_/pp4ljj/SM_Ltotal_Ind5v2020v2_UFO/pp4ljj_FT2_pureQGC/pp4ljj_pureQGC_FT2_0e-13_0e-13_9e-13_0e-13_0e-13_1_FT2toFT8_woHELICITY/bin/internal/extended_cmd.py", line 1514, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/_GENERATION_/pp4ljj/SM_Ltotal_Ind5v2020v2_UFO/pp4ljj_FT2_pureQGC/pp4ljj_pureQGC_FT2_0e-13_0e-13_9e-13_0e-13_0e-13_1_FT2toFT8_woHELICITY/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/_GENERATION_/pp4ljj/SM_Ltotal_Ind5v2020v2_UFO/pp4ljj_FT2_pureQGC/pp4ljj_pureQGC_FT2_0e-13_0e-13_9e-13_0e-13_0e-13_1_FT2toFT8_woHELICITY/bin/internal/common_run_interface.py", line 2142, in do_reweight
    reweight_cmd.import_command_file(path)
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1660, in import_command_file
    self.exec_cmd(line, precmd=True)
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1543, in exec_cmd
    stop = Cmd.onecmd_orig(current_interface, line, **opt)
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
    return func(arg, **opt)
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/various/misc.py", line 100, in f_with_no_logger
    out = f(self, *args, **opt)
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 546, in do_launch
    weight = self.calculate_weight(event)
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 968, in calculate_weight
    w_new = self.calculate_matrix_element(new_event, 1)
  File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 1157, in calculate_matrix_element
    orig_order, Pdir, hel_dict = self.id_to_path_second[tag]
KeyError: ((-1, 2), (-13, -11, -3, 4, 11, 13))
                              Run Options
                              -----------
               stdout_level : None

                         MadEvent Options
                         ----------------
     automatic_html_opening : False (user set)
        notification_center : True
          cluster_temp_path : /tmp/madgraph (user set)
             cluster_memory : None
               cluster_size : 100
              cluster_queue : None
                    nb_core : 32 (user set)
               cluster_time : None
                   run_mode : 1 (user set)

                      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 : /lustre/home/leisos/marantis/MG5_aMC_v2_6_5/Delphes (user set)
                thepeg_path : None (user set)
               cluster_type : slurm (user set)
          madanalysis5_path : None (user set)
           cluster_nb_retry : 30 (user set)
                 eps_viewer : None
                web_browser : None
               syscalc_path : None (user set)
           madanalysis_path : None (user set)
                     lhapdf : /lustre/home/leisos/marantis/MG5_aMC_v2_6_5/local_6.1.6/bin/lhapdf-config (user set)
              f2py_compiler : None
                 hepmc_path : None (user set)
         cluster_retry_wait : 300 (user set)
           fortran_compiler : None
                auto_update : 7 (user set)
        exrootanalysis_path : /lustre/home/leisos/marantis/MG5_aMC_v2_6_5/ExRootAnalysis (user set)
                    timeout : 60
               cpp_compiler : None

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

The key error means that one set of matrix-element does not have a matrix-element for a given initial/final state
in your case this is for
u d~ > mu+ mu- e+ e- s~ c

Now you use an outdated version of the code (2.6.5) if this is a bug on the reweighting side (which might be) then I would advise to check with a more recent version.

Cheers,

Olivier

> On 21 Feb 2021, at 20:20, Alexandros Marantis <email address hidden> wrote:
>
> New question #695684 on MadGraph5_aMC@NLO:
> https://answers.launchpad.net/mg5amcnlo/+question/695684
>
> Hello Olivier,
> I, systematically, get this error when I'm trying to re-weight between two different operators in this EFT model (http://feynrules.irmp.ucl.ac.be/wiki/AnomalousGaugeCoupling#no1), and more specifically, when I'm trying to re-weight between two different "groups" (T0,T1,T2 -> T8,T9 e.g re-weight a sample generated with non-zero T1 to non-zero T8).
> There is a KeyError at the end of the debug file. (i'm attaching it)
> Can you please help me to understand what "KeyError: ((-1, 2), (-13, -11, -3, 4, 11, 13))" means?
> In fact, what I'd like to know is if this error refers to a bug of the program or the inability of the model to re-weight between different groups of operators.
>
> Thank you in advance,
> Alexandros
>
> P.S The re-weighting process works successfully and with great precision, between the operators of the same group.
>
>
> //---------------------------- from ME5_debug ----------------------------//
> reweight run_01_0
> Traceback (most recent call last):
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/_GENERATION_/pp4ljj/SM_Ltotal_Ind5v2020v2_UFO/pp4ljj_FT2_pureQGC/pp4ljj_pureQGC_FT2_0e-13_0e-13_9e-13_0e-13_0e-13_1_FT2toFT8_woHELICITY/bin/internal/extended_cmd.py", line 1514, in onecmd
> return self.onecmd_orig(line, **opt)
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/_GENERATION_/pp4ljj/SM_Ltotal_Ind5v2020v2_UFO/pp4ljj_FT2_pureQGC/pp4ljj_pureQGC_FT2_0e-13_0e-13_9e-13_0e-13_0e-13_1_FT2toFT8_woHELICITY/bin/internal/extended_cmd.py", line 1463, in onecmd_orig
> return func(arg, **opt)
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/_GENERATION_/pp4ljj/SM_Ltotal_Ind5v2020v2_UFO/pp4ljj_FT2_pureQGC/pp4ljj_pureQGC_FT2_0e-13_0e-13_9e-13_0e-13_0e-13_1_FT2toFT8_woHELICITY/bin/internal/common_run_interface.py", line 2142, in do_reweight
> reweight_cmd.import_command_file(path)
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1660, in import_command_file
> self.exec_cmd(line, precmd=True)
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1543, in exec_cmd
> stop = Cmd.onecmd_orig(current_interface, line, **opt)
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/extended_cmd.py", line 1463, in onecmd_orig
> return func(arg, **opt)
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/various/misc.py", line 100, in f_with_no_logger
> out = f(self, *args, **opt)
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 546, in do_launch
> weight = self.calculate_weight(event)
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 968, in calculate_weight
> w_new = self.calculate_matrix_element(new_event, 1)
> File "/lustre/home/leisos/marantis/MG5_aMC_v2_6_5/madgraph/interface/reweight_interface.py", line 1157, in calculate_matrix_element
> orig_order, Pdir, hel_dict = self.id_to_path_second[tag]
> KeyError: ((-1, 2), (-13, -11, -3, 4, 11, 13))
> Run Options
> -----------
> stdout_level : None
>
> MadEvent Options
> ----------------
> automatic_html_opening : False (user set)
> notification_center : True
> cluster_temp_path : /tmp/madgraph (user set)
> cluster_memory : None
> cluster_size : 100
> cluster_queue : None
> nb_core : 32 (user set)
> cluster_time : None
> run_mode : 1 (user set)
>
> 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 : /lustre/home/leisos/marantis/MG5_aMC_v2_6_5/Delphes (user set)
> thepeg_path : None (user set)
> cluster_type : slurm (user set)
> madanalysis5_path : None (user set)
> cluster_nb_retry : 30 (user set)
> eps_viewer : None
> web_browser : None
> syscalc_path : None (user set)
> madanalysis_path : None (user set)
> lhapdf : /lustre/home/leisos/marantis/MG5_aMC_v2_6_5/local_6.1.6/bin/lhapdf-config (user set)
> f2py_compiler : None
> hepmc_path : None (user set)
> cluster_retry_wait : 300 (user set)
> fortran_compiler : None
> auto_update : 7 (user set)
> exrootanalysis_path : /lustre/home/leisos/marantis/MG5_aMC_v2_6_5/ExRootAnalysis (user set)
> timeout : 60
> cpp_compiler : None
>
> --
> You received this question notification because you are an answer
> contact for MadGraph5_aMC@NLO.

Revision history for this message
Alexandros Marantis (amaranti) said :
#2

Thanks Olivier Mattelaer, that solved my question.