madgraph5 / madgraph4gpu

GPU development for the Madgraph5_aMC@NLO event generator software package
30 stars 32 forks source link

SUSY sample generation failure in the SYCL plugin (fixed in cudacpp) #694

Open whhopkins opened 1 year ago

whhopkins commented 1 year ago

When attempting to generate a SUSY sample used by the ATLAS SUSY group (stop pair production) an error occurs. Below is the process card:

set group_subprocesses Auto
set ignore_six_quark_processes False
set max_t_for_channel 99
set loop_optimized_output True
set low_mem_multicore_nlo_generation False
set default_unset_couplings 99
set zerowidth_tchannel True
set loop_color_flows False
set gauge unitary
set complex_mass_scheme False
set max_npoint_for_channel 0
import model MSSM_SLHA2
define p = g u c d s u~ c~ d~ s~
define j = g u c d s u~ c~ d~ s~
define pb = g u c d s b u~ c~ d~ s~ b~
define jb = g u c d s b u~ c~ d~ s~ b~
define l+ = e+ mu+
define l- = e- mu-
define vl = ve vm vt
define vl~ = ve~ vm~ vt~
define fu = u c e+ mu+ ta+
define fu~ = u~ c~ e- mu- ta-
define fd = d s ve~ vm~ vt~
define fd~ = d~ s~ ve vm vt
define susystrong = go ul ur dl dr cl cr sl sr t1 t2 b1 b2 ul~ ur~ dl~ dr~ cl~ cr~ sl~ sr~ t1~ t2~ b1~ b2~
define susyweak = el- el+ er- er+ mul- mul+ mur- mur+ ta1- ta1+ ta2- ta2+ n1 n2 n3 n4 x1- x1+ x2- x2+ sve sve~ svm svm~ svt svt~
define susylq = ul ur dl dr cl cr sl sr
define susylq~ = ul~ ur~ dl~ dr~ cl~ cr~ sl~ sr~
define susysq = ul ur dl dr cl cr sl sr t1 t2 b1 b2
define susysq~ = ul~ ur~ dl~ dr~ cl~ cr~ sl~ sr~ t1~ t2~ b1~ b2~
define susysl = el- el+ er- er+ mul- mul+ mur- mur+ ta1- ta1+ ta2- ta2+
define susyv = sve svm svt
define susyv~ = sve~ svm~ svt~
generate p p > t1 t1~ $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @1
add process p p > t1 t1~ j $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @2
add process p p > t1 t1~ j j $ go susylq susylq~ b1 b2 t2 b1~ b2~ t2~ @2
output -f -nojpeg

The error that occurs is:

output madevent CODEGEN_mad_susy_model --hel_recycling=False --vector_size=16384 --me_exporter=standalone_sycl" with error:
ValueError : not enough values to unpack (expected 2, got 1)
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 error even occurs when attempting to generate only events with stop pair production without additional jets, i.e.: import model MSSM_SLHA2; generate g g > t t~

The branch that was used when attempting to generate these SUSY events was gpucpp

nscottnichols commented 1 year ago

@valassi may already be tracking this issue elsewhere (it may even be already fixed in the CUDA branch)

valassi commented 1 year ago

Hi @whhopkins @nscottnichols indeed this was #622 in cudacpp, now fixed - meaning, code generation succeeds. There are still other issues at build time though I am afraid for cudacpp.

Anyway, I reassign this #694 to @nscottnichols as this is sycl specific. Thanks!

valassi commented 1 year ago

For cudacpp it was fixed in MR #624