Problem with process definition

Asked by Daniel William Busbridge

Dear MadGraph experts,

I am currently trying to generate same sign dimuon events with MadGraph 2.1.2 using the MSSM model generated with SARAH 4.3.0. My process card is

import model MSSM -modelname
# Define multiparticle labels
define p = g u1 u2 d1 d2 u1bar u2bar d1bar d2bar
define j = g u1 u2 d1 d2 u1bar u2bar d1bar d2bar
define l- = e1 e2
define l+ = e1bar e2bar
define vl = nu1 nu2 nu3
define vl~ = nu1bar nu2bar nu3bar
# Specify process(es) to run
generate p p > su2 su2, (su2 > d2 c1bar, c1bar > n1 wmc, wmc > nu2 e2bar)
# Output processes to MadEvent directory
# output
# This will create a directory PROC_$MODELNAME_$X
# If you want to specify the path/name of the directory use
output MSSM_same_sign_dimuon

# To generate events, you can go to the created directory and
# run ./bin/generate_events

In the above notation, su2 is the second lightest up squark, d2 is the strange quark, c1bar is the positively charged lightest chargino, n1 is the lightest neutralino, wmc is the positive W boson, nu2 is the second lightest neutrino and e2bar is the positively charged muon.

The output from running the process card is:

-----------------------------------------------------------------------------

load MG5 configuration from input/mg5_configuration.txt
set lhapdf to /mt/batch/users/busbridge/LHAPDF5/bin/lhapdf-config
set fastjet to /mt/batch/users/busbridge/fastjet/bin/fastjet-config
Using default text editor "vi". Set another one in ./input/mg5_configuration.txt
Using default eps viewer "evince". Set another one in ./input/mg5_configuration.txt
Using default web browser "firefox". Set another one in ./input/mg5_configuration.txt
import /mt/batch/users/busbridge/MG5_aMC_v2_1_2/proc_card_MSSM_same_sign_dimuon.dat
The import format was not given, so we guess it as command
import model MSSM -modelname
impossible to set default multiparticles p because No particle u in model
impossible to set default multiparticles j because No particle u in model
impossible to set default multiparticles l+ because No particle e+ in model
impossible to set default multiparticles l- because No particle e- in model
impossible to set default multiparticles vl because No particle ve in model
impossible to set default multiparticles vl~ because No particle ve~ in model
Defined multiparticle all = g a nu1 nu2 nu3 nu1bar nu2bar nu3bar go sd1 sd2 sd3 sd4 sd5 sd6 su1 su2 su3 su4 su5 su6 sd1c sd2c sd3c sd4c sd5c sd6c su1c su2c su3c su4c su5c su6c d1 d2 d3 u1 u2 u3 d1bar d2bar d3bar u1bar u2bar u3bar sv1 sv2 sv3 se1 se2 se3 se4 se5 se6 h1 h2 ah2 hm2c z wmc sv1c sv2c sv3c se1c se2c se3c se4c se5c se6c hm2 wm n1 n2 n3 n4 c1bar c2bar e1 e2 e3 c1 c2 e1bar e2bar e3bar
define p = g u1 u2 d1 d2 u1bar u2bar d1bar d2bar
Defined multiparticle p = g u1 u2 d1 d2 u1bar u2bar d1bar d2bar
define j = g u1 u2 d1 d2 u1bar u2bar d1bar d2bar
Defined multiparticle j = g u1 u2 d1 d2 u1bar u2bar d1bar d2bar
define l- = e1 e2
Defined multiparticle l- = e1 e2
define l+ = e1bar e2bar
Defined multiparticle l+ = e1bar e2bar
define vl = nu1 nu2 nu3
Defined multiparticle vl = nu1 nu2 nu3
define vl~ = nu1bar nu2bar nu3bar
Defined multiparticle vl~ = nu1bar nu2bar nu3bar
generate p p > su2 su2, (su2 > d2 c1bar, c1bar > n1 wmc, wmc > nu2 e2bar)
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED=2
INFO: Trying process: u1 u1 > su2 su2 WEIGHTED=2
INFO: Process has 2 diagrams
INFO: Trying process: u1 u2 > su2 su2 WEIGHTED=2
INFO: Process has 2 diagrams
INFO: Process u2 u1 > su2 su2 added to mirror process u1 u2 > su2 su2
INFO: Trying process: u2 u2 > su2 su2 WEIGHTED=2
INFO: Process has 2 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: su2 > d2 c1bar WEIGHTED=2
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: c1bar > n1 wmc WEIGHTED=2
INFO: Process has 1 diagrams
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying process: wmc > nu2 e2bar WEIGHTED=2
INFO: Process has 1 diagrams
WARNING: Decay without corresponding particle in core process found.
Decay information for particle(s) wmc is discarded.
Please check your process definition carefully.
This warning usually means that you forgot parentheses in presence of subdecay.
Example of correct syntax: p p > t t~, ( t > w+ b, w+ > l+ vl)
1 processes with 8 diagrams generated in 6.899 s
Total: 1 processes with 8 diagrams
output MSSM_same_sign_dimuon
INFO: initialize a new directory: MSSM_same_sign_dimuon
INFO: remove old information in MSSM_same_sign_dimuon
INFO: Organizing processes into subprocess groups
INFO: Organizing processes into subprocess groups
INFO: Organizing processes into subprocess groups
INFO: Generating Helas calls for process: u1 u1 > su2 su2 WEIGHTED=2
INFO: Generating Helas calls for process: u1 u2 > su2 su2 WEIGHTED=2
INFO: Generating Helas calls for process: u2 u2 > su2 su2 WEIGHTED=2
INFO: Generating Helas calls for process: su2 > d2 c1bar WEIGHTED=2
INFO: Generating Helas calls for process: c1bar > n1 wmc WEIGHTED=2
INFO: Combine su2 > d2 c1bar WEIGHTED=2 with decays c1bar > n1 wmc WEIGHTED=2
INFO: Combine u1 u1 > su2 su2 WEIGHTED=2 with decays su2 > d2 c1bar WEIGHTED=2
  Decay: c1bar > n1 wmc WEIGHTED=2, su2 > d2 c1bar WEIGHTED=2
  Decay: c1bar > n1 wmc WEIGHTED=2
INFO: Combine u1 u2 > su2 su2 WEIGHTED=2 with decays su2 > d2 c1bar WEIGHTED=2
  Decay: c1bar > n1 wmc WEIGHTED=2, su2 > d2 c1bar WEIGHTED=2
  Decay: c1bar > n1 wmc WEIGHTED=2
INFO: Combine u2 u2 > su2 su2 WEIGHTED=2 with decays su2 > d2 c1bar WEIGHTED=2
  Decay: c1bar > n1 wmc WEIGHTED=2, su2 > d2 c1bar WEIGHTED=2
  Decay: c1bar > n1 wmc WEIGHTED=2
INFO: Processing color information for process: u1 u1 > su2 su2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
INFO: Reusing existing color information for process: u1 u2 > su2 su2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
INFO: Reusing existing color information for process: u2 u2 > su2 su2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
INFO: Creating files in directory P0_u1u1_su2su2_su2_d2c1bar_c1bar_n1wmc_su2_d2c1bar_c1bar_n1wmc
INFO: Generating Feynman diagrams for Process: u1 u1 > su2 su2 WEIGHTED=2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
INFO: Finding symmetric diagrams for subprocess group u1u1_su2su2_su2_d2c1bar_c1bar_n1wmc_su2_d2c1bar_c1bar_n1wmc
INFO: Creating files in directory P0_u1u2_su2su2_su2_d2c1bar_c1bar_n1wmc_su2_d2c1bar_c1bar_n1wmc
INFO: Generating Feynman diagrams for Process: u1 u2 > su2 su2 WEIGHTED=2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
INFO: Finding symmetric diagrams for subprocess group u1u2_su2su2_su2_d2c1bar_c1bar_n1wmc_su2_d2c1bar_c1bar_n1wmc
INFO: Creating files in directory P0_u2u2_su2su2_su2_d2c1bar_c1bar_n1wmc_su2_d2c1bar_c1bar_n1wmc
INFO: Generating Feynman diagrams for Process: u2 u2 > su2 su2 WEIGHTED=2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
  Decay: su2 > d2 c1bar WEIGHTED=2
    Decay: c1bar > n1 wmc WEIGHTED=2
INFO: Finding symmetric diagrams for subprocess group u2u2_su2su2_su2_d2c1bar_c1bar_n1wmc_su2_d2c1bar_c1bar_n1wmc
History written to /mt/batch/users/busbridge/MG5_aMC_v2_1_2/MSSM_same_sign_dimuon/Cards/proc_card_mg5.dat
Generated helas calls for 3 subprocesses (6 diagrams) in 0.090 s
Wrote files for 48 helas calls in 0.504 s
Export UFO model to MG4 format
madALOHA: aloha creates FFV2 routines
ALOHA: aloha creates FFV3 routines
ALOHA: aloha creates FFS2 routines
ALOHA: aloha creates FFS1 routines
save configuration file to /mt/batch/users/busbridge/MG5_aMC_v2_1_2/MSSM_same_sign_dimuon/Cards/me5_configuration.txt
INFO: Use Fortran compiler gfortran
INFO: Generate jpeg diagrams
INFO: Generate web pages
Output to directory /mt/batch/users/busbridge/MG5_aMC_v2_1_2/MSSM_same_sign_dimuon done.
Type "launch" to generate events from this process, or see
/mt/batch/users/busbridge/MG5_aMC_v2_1_2/MSSM_same_sign_dimuon/README
Run "open index.html" to see more information about this process.
quit

-----------------------------------------------------------------------------

The diagrams don't include the W decay to a muon and a neutrino. Generating events without the W decay

generate p p > su2 su2, (su2 > d2 c1bar, c1bar > n1 wmc)

and generating the W decay on its own

generate wmc > nu2 e2bar

work fine. Any help you could provide in generating these di-muon events would be appreciated.

Many thanks,

Daniel

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Rikkert Frederix
Solved:
Last query:
Last reply:
Revision history for this message
Best Rikkert Frederix (frederix) said :
#1

Dear Daniel,

The code gives you this warning:

WARNING: Decay without corresponding particle in core process found.
Decay information for particle(s) wmc is discarded.
Please check your process definition carefully.
This warning usually means that you forgot parentheses in presence of subdecay.

Hence, I guess the correct syntax should be:

generate p p > su2 su2, (su2 > d2 c1bar, (c1bar > n1 wmc, wmc > nu2 e2bar))

or

generate p p > su2 su2, ((su2 > d2 c1bar, c1bar > n1 wmc), wmc > nu2 e2bar)

Does that solve the problem?

Best regards,
Rikkert

Revision history for this message
Daniel William Busbridge (d-w-busbridge) said :
#2

Thanks Rikkert Frederix, that solved my question.