Problem with process definition
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_
# To generate events, you can go to the created directory and
# run ./bin/generate_
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_
set lhapdf to /mt/batch/
set fastjet to /mt/batch/
Using default text editor "vi". Set another one in ./input/
Using default eps viewer "evince". Set another one in ./input/
Using default web browser "firefox". Set another one in ./input/
import /mt/batch/
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_
INFO: initialize a new directory: MSSM_same_
INFO: remove old information in MSSM_same_
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_
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_
INFO: Creating files in directory P0_u1u2_
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_
INFO: Creating files in directory P0_u2u2_
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_
History written to /mt/batch/
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/
INFO: Use Fortran compiler gfortran
INFO: Generate jpeg diagrams
INFO: Generate web pages
Output to directory /mt/batch/
Type "launch" to generate events from this process, or see
/mt/batch/
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
- Assignee:
- No assignee Edit question
- Solved by:
- Rikkert Frederix
- Solved:
- Last query:
- Last reply: