color tracing problem with pythia 8
I am trying to generate some events with monotop final state coming from a minimal extension of the standard model, including heavy colour-triplet scalars interacting with a light Majorana ferminon that can be dark matter (arxiv:1507.02271). The scalars are x1 (particle code 6000001) and x2 (particle code 6000002) and fermion is n (particle code 5000001). When I ask for extra jets in madgraph 2.2.3, some of the events fail in pythia 8210. The error messages are:
9420 Error in ColourTracing:
4529 Error in ColourTracing:
181 Error in HadronLevel::next: failed colour/junction check
4710 Error in Pythia::next: hadronLevel failed; try again
So it seems the problem is in color tracing. The events that fail are all from incoming (anti)quarks with a gluon coming from the incoming quarks or x, the lhe file does not have the complete information, so I cannot match the exact diagram. One suspect diagram is: incoming d~ and s~, d~ emits a gluon, d~ and s~ produce x1, x1 decays into t and n, t decays to W and b. The lhe file looks like this:
<event>
...
-1 -1 0 0 0 501 …
-3 -1 0 0 0 502 …
6000001 2 1 2 504 0 ...
6 2 3 3 504 0 ...
5000001 1 3 3 0 0 ...
24 1 4 4 0 0 ...
5 1 4 4 504 0 ...
21 1 1 2 503 501 ...
</event>
It seems the gluon doesn't carry the right anticolor; changing 501 to 503 solves the issue; the event passes through pythia.
So my question is: is this a problem or am I missing something in the model that is needed by madgraph to work out the color assignments?
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Nikolay Kolev for more information if necessary.