Decay chain with diagram filter returns NoDiagramException
Hello MG5 experts,
I'm trying to
generate p p > t t~ sig8 NP=2, sig8 > t t~ --diagram_filter
where sig8 is a top-philic particle; the diagram filter is discarding all the diagrams with two sig8 in the s-channel, as I want to ban all the p p > sig8 sig8 > ... contributions from my sample (see the filter below).
It seems the filter is removing the decay diagram, as no diagrams are found in the end:
Interpreting 'NP=2' as 'NP<=2'
INFO: Trying process: g g > t t~ sig8 NP<=2 @1
WARNING: Diagram filter is ON and removed 4 diagrams for this subprocess.
INFO: Process has 23 diagrams
INFO: Trying process: u u~ > t t~ sig8 NP<=2 @1
WARNING: Diagram filter is ON and removed 1 diagrams for this subprocess.
INFO: Process has 7 diagrams
INFO: Trying process: u c~ > t t~ sig8 NP<=2 @1
INFO: Trying process: c u~ > t t~ sig8 NP<=2 @1
INFO: Trying process: c c~ > t t~ sig8 NP<=2 @1
WARNING: Diagram filter is ON and removed 1 diagrams for this subprocess.
INFO: Process has 7 diagrams
INFO: Trying process: d d~ > t t~ sig8 NP<=2 @1
WARNING: Diagram filter is ON and removed 1 diagrams for this subprocess.
INFO: Process has 7 diagrams
INFO: Trying process: d s~ > t t~ sig8 NP<=2 @1
INFO: Trying process: s d~ > t t~ sig8 NP<=2 @1
INFO: Trying process: s s~ > t t~ sig8 NP<=2 @1
WARNING: Diagram filter is ON and removed 1 diagrams for this subprocess.
INFO: Process has 7 diagrams
INFO: Process u~ u > t t~ sig8 added to mirror process u u~ > t t~ sig8
INFO: Process c~ c > t t~ sig8 added to mirror process c c~ > t t~ sig8
INFO: Process d~ d > t t~ sig8 added to mirror process d d~ > t t~ sig8
INFO: Process s~ s > t t~ sig8 added to mirror process s s~ > t t~ sig8
INFO: Checking for minimal orders which gives processes.
INFO: Please specify coupling orders to bypass this step.
INFO: Trying coupling order WEIGHTED<=1: WEIGTHED IS NP+QCD+2*QED
WARNING: Diagram filter is ON and removed 1 diagrams for this subprocess.
INFO: Trying coupling order WEIGHTED<=2: WEIGTHED IS NP+QCD+2*QED
WARNING: Diagram filter is ON and removed 1 diagrams for this subprocess.
INFO: Trying coupling order WEIGHTED<=3: WEIGTHED IS NP+QCD+2*QED
WARNING: Diagram filter is ON and removed 1 diagrams for this subprocess.
INFO: Trying process: sig8 > t t~ WEIGHTED<=4
WARNING: Diagram filter is ON and removed 1 diagrams for this subprocess.
Command "generate p p > t t~ sig8 NP=2, sig8 > t t~ --diagram_filter" interrupted with error:
NoDiagramException : No amplitudes generated from process Process: sig8 > t t~ WEIGHTED=4. Please enter a valid process
How is this possible? There is only one sig8 in the decay process, what am I missing?
The same error appears when trying to decay the p p > t t~ sig8 sample through MadSpin, as the diagram filter is still active at that level.
Is there a way to solve this?
Thank you,
Matteo
P.S. Here's my diagram filter:
import madgraph.
def remove_diag (diag,model):
'''
Remove all diagrams with two sig8 in s-channel from
generate p p > t t~ sig8 NP=2 QED=0
'''
draw = drawing.
draw.
draw.
sig8_sch = 0
for p in draw.lineList:
if p.id==9000001 and p.begin.
if sig8_sch > 1:
return True
else:
return False
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- Last query:
- Last reply: