Problem generating VLQ tp+j

Dear experts

I have a question regarding VLQ NLO tp+j production. In particular, we need to generate some single tp +j events and we tried to follow the instructions on the Appendix of https://arxiv.org/pdf/2105.08745.pdf which uses something like

import model VLQ_v4_4FNS_UFO-3rdL_noCKM
set complex_mass_scheme

generate p p > tpp bb j / bp x y [QCD] QED=1 QCD=1

but that fails with the error [1]

anyway, we also noticed that as explained in the text, the NLO diagrams need special modifications of the MG5 core files and it refers the 1610.04622 where it is stated in pg11:

"As for electroweak contributions to vector-like quark pair production, mixed QCD and electroweak loop diagrams appear at the NLO level and must be treated consistently for getting ultravioletfinite results. This step being not automated, we provide scripts that steer the event generation process on the UFO model webpage "

I guess “model webpage” refers to the feyrules/NLOModels, and there it seems like that there are like VLQ-V5_full, 4FNS_UFO, 5FNS_UFO, generation_scripts tarball.

So, I took the event_generation_scipts (scripts.tar.gz) and I saw that there is a generation_pp_tpj but there, the proc_card calls the VLQ_NLO_UFO model and not the VLQ_v4_4FNS_UFO one. In any case, I tried to use the generate.sh script from the generate_pp_tpj folder which I guess handles the mixed QCD/EWK loop diagrams first and then tries to generate events. However, I then get this error I(we have to use MG2.9.9) [2].

So, in short, we followed the syntax/model of a more recent paper (2105.08745) than 1610.04622 that discusses the VLQ production but then when tried to treat properly the loop diagrams by using the scripts that are provided in the feynrules pages we got stuck.

That said, could you please help us understand what model (ie VLQ_v4_4FNS_UFO, VLQ_v5_UFO or the VLQ_NLO_UFO for example?) and what modifications are needed and finally which syntax is the right one to get the single tp+j production @NLO (like p p > tpp bb j / bp x y [QCD] QED=1 QCD=1 or generate p p > ttp j [QCD]?)

I hope my email details all that you need to know to provide a solution.

Regards and thanks

Alexis

[1]

generate p p > tpp bb j / bp x y [QCD] QED=1 QCD=1
Interpreting 'QCD=1' as 'QCD<=1'
Interpreting 'QED=1' as 'QED<=1'
INFO: Generating FKS-subtracted matrix elements for born process: g u > tp b~ d QCD<=1 QED<=1 [ all = QCD ] / bp x y (1 / 24)
Command "import /afs/cern.ch/work/a/alkaloge/Gridpacks/VLQ/genproductions/bin/MadGraph5_aMCatNLO/bW_t_NLO/bW_t_NLO_gridpack/work/bW_t_NLO_proc_card.dat" interrupted in sub-command:
"generate p p > tpp bb j / bp x y [QCD] QED=1 QCD=1" with error:
FKSProcessError : Cannot map born/real configurations between g u > tp b~ d QCD<=1 QED<=1 [ all = QCD ] / bp x y and g u > tp b~ d g QCD<=2 QED<=1 WEIGHTED<=4 [ all = QCD ] / bp x y (i,j=6,1): not same number of configurations: 4 0
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

[2]
 File "/afs/cern.ch/work/a/alkaloge/Gridpacks/VLQ/genproductions/bin/MadGraph5_aMCatNLO/bW_t_NLO/bW_t_NLO_gridpack/work/./MG5_aMC_v2_9_9/bin/mg5_aMC", line 148, in <module>
    import madgraph.interface.master_interface as interface
  File "/afs/cern.ch/work/a/alkaloge/Gridpacks/VLQ/genproductions/bin/MadGraph5_aMCatNLO/bW_t_NLO/bW_t_NLO_gridpack/work/MG5_aMC_v2_9_9/madgraph/interface/master_interface.py", line 40, in <module>
    import madgraph.core.diagram_generation as diagram_generation
  File "/afs/cern.ch/work/a/alkaloge/Gridpacks/VLQ/genproductions/bin/MadGraph5_aMCatNLO/bW_t_NLO/bW_t_NLO_gridpack/work/MG5_aMC_v2_9_9/madgraph/core/diagram_generation.py", line 32, in <module>
    import madgraph.core.base_objects as base_objects
  File "/afs/cern.ch/work/a/alkaloge/Gridpacks/VLQ/genproductions/bin/MadGraph5_aMCatNLO/bW_t_NLO/bW_t_NLO_gridpack/work/MG5_aMC_v2_9_9/madgraph/core/base_objects.py", line 84
    raise self.PhysicsObjectError, \
                                 ^
SyntaxError: invalid syntax