MadGraph5_aMC@NLO

Error detected in "import model SDMV2" write debug file MG5_debug If you need help with this issue pleasstr : b has not coherent incoming/outcoming status between interactions

Asked by Debajyoti

Hi!

I tried to load a model into Madgraph after generating its UFO file. However on typing the import command, I get the following message:

MG5_aMC>import model SDMV2
INFO: load particles
Error detected in "import model SDMV2"
write debug file MG5_debug
If you need help with this issue please contact us on https://answers.launchpad.net/mg5amcnlo
str : b has not coherent incoming/outcoming status between interactions

The MG5_debug file reads:

set group_subprocesses Auto
set ignore_six_quark_processes False
set low_mem_multicore_nlo_generation False
set complex_mass_scheme False
set include_lepton_initiated_processes False
set gauge unitary
set loop_optimized_output True
set loop_color_flows False
set max_npoint_for_channel 0
set default_unset_couplings 99
set max_t_for_channel 99
set zerowidth_tchannel True
set nlo_mixed_expansion True
import model sm
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
import model SDMV2
Traceback (most recent call last):
  File "/home/djb/Downloads/env3.8/Software/MG5_aMC_v3.5.3/MG5_aMC_v3_5_3/madgraph/interface/extended_cmd.py", line 1546, in onecmd
    return self.onecmd_orig(line, **opt)
  File "/home/djb/Downloads/env3.8/Software/MG5_aMC_v3.5.3/MG5_aMC_v3_5_3/madgraph/interface/extended_cmd.py", line 1495, in onecmd_orig
    return func(arg, **opt)
  File "/home/djb/Downloads/env3.8/Software/MG5_aMC_v3.5.3/MG5_aMC_v3_5_3/madgraph/interface/master_interface.py", line 281, in do_import
    self.cmd.do_import(self, *args, **opts)
  File "/home/djb/Downloads/env3.8/Software/MG5_aMC_v3.5.3/MG5_aMC_v3_5_3/madgraph/interface/madgraph_interface.py", line 5619, in do_import
    self._curr_model = import_ufo.import_model(args[1], prefix=prefix,
  File "/home/djb/Downloads/env3.8/Software/MG5_aMC_v3.5.3/MG5_aMC_v3_5_3/models/import_ufo.py", line 250, in import_model
    model = import_full_model(model_path, decay, prefix)
  File "/home/djb/Downloads/env3.8/Software/MG5_aMC_v3.5.3/MG5_aMC_v3_5_3/models/import_ufo.py", line 407, in import_full_model
    model = ufo2mg5_converter.load_model()
  File "/home/djb/Downloads/env3.8/Software/MG5_aMC_v3.5.3/MG5_aMC_v3_5_3/models/import_ufo.py", line 562, in load_model
    self.detect_incoming_fermion()
  File "/home/djb/Downloads/env3.8/Software/MG5_aMC_v3.5.3/MG5_aMC_v3_5_3/models/import_ufo.py", line 1336, in detect_incoming_fermion
    raise InvalidModel('%s has not coherent incoming/outcoming status between interactions' %\
models.import_ufo.InvalidModel: b has not coherent incoming/outcoming status between interactions
                          MadGraph5_aMC@NLO Options
                          ----------------
                auto_update : 7
        complex_mass_scheme : False
    default_unset_couplings : 99
                      gauge : unitary
         group_subprocesses : Auto
  ignore_six_quark_processes : False
  include_lepton_initiated_processes : False
           loop_color_flows : False
      loop_optimized_output : True
  low_mem_multicore_nlo_generation : False
     max_npoint_for_channel : 0
          max_t_for_channel : 99
        nlo_mixed_expansion : True
               stdout_level : 20 (user set)
         zerowidth_tchannel : True

                         MadEvent Options
                          ----------------
     automatic_html_opening : True
                    nb_core : None
        notification_center : True
                   run_mode : 2

I hope I have posted the relevant part of the debug file.

I hope this information suffices in pointing out the problem I am facing. Any help in figuring this out will be highly appreciated. Thanks.

Sincerely,
djb

Question information

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

Hi,

This means that your UFO model is not valid and that some of the lorentz structure for interactions with the b are not consistent with some other lorentz interactions with the "b".

Cheers,

Olivier

Revision history for this message
Debajyoti (djb24) said : 		#2

Since I have also treated other quarks the same, I do not see any particular reason for b-interactions to be different. Anyway, I would be happy to be corrected.

However, on trying to import the same model another time, I receive a slightly different message. Maybe this could help.

MG5_aMC>import model SDMV2
fail to load model but auto_convert_model is on True. Trying to convert the model
convert model /home/djb/Downloads/env3.8/Software/MG5_aMC_v3.5.3/MG5_aMC_v3_5_3/models/SDMV2
retry the load of the model
import model SDMV2
INFO: load particles
Command "import model SDMV2" interrupted with error:
UFOError : invalid syntax (object_library.py, line 268)
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

The object_library.py file reads:

def pole(self, x):
        """ the self.value attribute can be a dictionary directly specifying the Laurent serie using normal
        parameter or just a string which can possibly contain CTparameter defining the Laurent serie."""

        if isinstance(self.value,dict):
            if -x in self.value.keys():
                return self.value[-x]
            else:
                return 'ZERO'

        CTparam=None
        for param in all_CTparameters:
           pattern=re.compile(r"(?P<first>\A|\*|\+|\-|\()(?P<name>"+param.name+r")(?P<second>\Z|\*|\+|\-|\))")
           numberOfMatches=len(pattern.findall(self.value))
           if numberOfMatches==1:
               if not CTparam:
                   CTparam=param
               else:
                   raise UFOError("UFO does not support yet more than one occurence of CTParameters in the couplings values.")
           elif numberOfMatches>1:
               raise UFOError("UFO does not support yet more than one occurence of CTParameters in the couplings values.")

        if not CTparam:
            if x==0:
                return self.value
            else:
                return 'ZERO'
        else:
            if CTparam.pole(x)=='ZERO':
                return 'ZERO'
            else:
                def substitution(matchedObj):
                    return matchedObj.group('first')+"("+CTparam.pole(x)+")"+matchedObj.group('second')
                pattern=re.compile(r"(?P<first>\A|\*|\+|\-|\()(?P<name>"+CTparam.name+r")(?P<second>\Z|\*|\+|\-|\))")
                return pattern.sub(substitution,self.value)

NOTE: Line 268: raise UFOError("UFO does not support yet more than one occurence of CTParameters in the couplings values.") *after else:*

I hope this offers a deeper insight into this issue. Thanks a lot.

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

Well maybe your model is wrong for all quark and the model crash on "b" because this is the first one that MG5aMC is testing.

When developping UFO model I would advise to turn auto_convert_model to False to
1) be sure that your model is python3 compatible
2) avoid that the converter runs multiple times on your model and corrupt it (likely what happened here)

Cheers,

Olivier

Revision history for this message
Debajyoti (djb24) said : 		#4

Thanks for the response.

I realized the problem was in not specifying the interaction order of new physics. Also, there were syntax errors in the proposed Lagrangian.

Anyway, thanks !