2-2, decay 1-2, output to pythia8 program, the program is wrong
Affects | Status | Importance | Assigned to | Milestone | |
---|---|---|---|---|---|
MadGraph5_aMC@NLO |
Fix Released
|
Medium
|
Olivier Mattelaer |
Bug Description
I inputed my model with extra squarks, and wanted to generate pythia8 program with such process:
u~ g > l- squak, squark > l+ u~
and this is a 2->5 process finally. I then typed
output pythia8
After compiled, the generated program showed warning telling me that the colour arrange is wrong.
After I went through the program generated automatically and discovered where the mistake come from:
in each function void ???::setIdColAcol() the code is shown below such as:
if(id1 == 21 && id2 == -2)
{
// Pick one of the flavor combinations (13, -6013)
int flavors[1][3] = {{13, -6013}}; //-6013 is the particle to decay into u~ and lepton
vector<double> probs;
double sum = matrix_element[0];
probs.
int choice = rndmPtr-
id3 = flavors[choice][0];
id4 = flavors[choice][1];
id5 = flavors[choice][2];
}
setId(id1, id2, id3, id4, id5);
// Pick color flow
int ncolor[1] = {1};
if(id1 == 21 && id2 == -2 && id3 == 13 && id4 -6013)
{
vector<double> probs;
double sum = jamp2[0][0];
for(int i = 0; i < ncolor[0]; i++ )
probs.
int ic = rndmPtr-
static int colors[1][10] = {{1, 2, 0, 1, 0, 0, 0, 0, 0, 2}};
setColAcol(
}
-6013 is the new particle I inserted in this model. It is obvious that this program wrongly set id5=0(void), and assigned the wrong color flow since it fully ignored the decay of -6013 as I have inputed in the previous madgraph generate command. So I changed every program by hand into such as:
if(id1 == 21 && id2 == -2)
{
// Pick one of the flavor combinations (13, -6013)
int flavors[1][3] = {{13, -11, -2}}; //-6013
vector<double> probs;
double sum = matrix_element[0];
probs.
int choice = rndmPtr-
id3 = flavors[choice][0];
id4 = flavors[choice][1];
id5 = flavors[choice][2];
}
setId(id1, id2, id3, id4, id5);
// Pick color flow
int ncolor[1] = {1};
if(id1 == 21 && id2 == -2 && id3 == 13 && id4 == -11 && id5 == -2)
{
vector<double> probs;
double sum = jamp2[0][0];
for(int i = 0; i < ncolor[0]; i++ )
probs.
int ic = rndmPtr-
static int colors[1][10] = {{1, 2, 0, 1, 0, 0, 0, 0, 0, 2}};
setColAcol(
}
Finnally, it worked.
Related branches
- Johan Alwall (community): Approve
- Olivier Mattelaer: Approve
-
Diff: 7090 lines (+3695/-1280)50 files modifiedTemplate/Cards/delphes_card_ATLAS.dat (+498/-134)
Template/Cards/delphes_card_CMS.dat (+511/-134)
Template/Cards/delphes_card_default.dat (+69/-53)
Template/Cards/delphes_trigger_ATLAS.dat (+0/-16)
Template/Cards/delphes_trigger_CMS.dat (+0/-20)
Template/Cards/delphes_trigger_default.dat (+0/-20)
Template/Source/PDF/pdf_lhapdf.f (+9/-2)
Template/Source/PDF/pdg2pdf.f (+52/-17)
Template/Source/PDF/pdg2pdf_lhapdf.f (+90/-14)
Template/SubProcesses/cuts.f (+16/-0)
Template/SubProcesses/myamp.f (+3/-1)
Template/bin/internal/run_delphes3 (+38/-0)
UpdateNotes.txt (+46/-0)
aloha/aloha_lib.py (+1/-1)
aloha/aloha_parsers.py (+5/-1)
aloha/aloha_writers.py (+1/-0)
madgraph/VERSION (+1/-1)
madgraph/core/base_objects.py (+30/-10)
madgraph/core/diagram_generation.py (+9/-3)
madgraph/core/helas_objects.py (+74/-45)
madgraph/interface/madevent_interface.py (+178/-27)
madgraph/interface/madgraph_interface.py (+23/-4)
madgraph/iolibs/export_cpp.py (+30/-8)
madgraph/iolibs/export_v4.py (+16/-6)
madgraph/iolibs/template_files/auto_dsig_v4.inc (+4/-1)
madgraph/iolibs/template_files/super_auto_dsig_group_v4.inc (+5/-0)
madgraph/various/cluster.py (+24/-17)
madgraph/various/lhe_parser.py (+246/-0)
models/EWdim6/EWdim6.log (+20/-10)
models/EWdim6/__init__.py (+23/-7)
models/EWdim6/build_restrict.py (+20/-16)
models/EWdim6/coupling_orders.py (+3/-3)
models/EWdim6/couplings.py (+46/-70)
models/EWdim6/function_library.py (+9/-1)
models/EWdim6/lorentz.py (+128/-113)
models/EWdim6/object_library.py (+22/-10)
models/EWdim6/parameters.py (+70/-170)
models/EWdim6/particles.py (+52/-52)
models/EWdim6/vertices.py (+215/-251)
models/EWdim6/write_param_card.py (+2/-2)
models/check_param_card.py (+1/-2)
models/import_ufo.py (+45/-33)
models/write_param_card.py (+2/-2)
tests/acceptance_tests/test_cmd_madevent.py (+60/-1)
tests/unit_tests/core/test_base_objects.py (+3/-1)
tests/unit_tests/core/test_helas_objects.py (+503/-0)
tests/unit_tests/iolibs/test_export_cpp.py (+283/-0)
tests/unit_tests/iolibs/test_export_v4.py (+12/-1)
tests/unit_tests/various/test_aloha.py (+63/-0)
tests/unit_tests/various/test_lhe_parser.py (+134/-0)
Changed in madgraph5: | |
importance: | Undecided → Medium |
Changed in madgraph5: | |
status: | In Progress → Fix Committed |
Changed in madgraph5: | |
status: | Fix Committed → Fix Released |
sorry, I made a mistake in my description, it is a 2 -> 3 process, not 2->5.