parameter scan fails after complaining about 3 body decay
Hi Olivier,
I am running a 2 dim parameter scan using a heavy Neutrino model file.
My scan used to run fine in the past but now (2.7.2) it stops after two points.
Both points give me an exception:
Exception : Phase Space generator not yet ready for 3 body decay
Do you have any idea how I can debug this?
I modified the phase space.
But also unmodified scripts fails.
My script starts with:
import model ../mg5amcnlo/
define lep = e+ e- mu+ mu- ta+ ta-
define ferm = lep ve ve~ vm vm~ vt vt~ u u~ d d~ s s~ c c~ b b~
generate p p > ferm n1
output out
launch out
shower = PY8
madspin = none
decay n1 > ferm ferm ferm
set ebeam1 7000
set ebeam2 7000
set nevents 10000
set ptl 0
set etal -1.49
set etalmin -3.44
set drll 0
set mpi False
set use_syst False
set time_of_flight 0.1
set mN2 1000000
set mN3 1000000
set VeN1 0
set VtaN1 0
set VeN2 0
set VmuN2 0
set VtaN2 0
set VeN3 0
set VmuN3 0
set VtaN3 0
# U2 0.0100000000000
# M 1.
set mN1 1.
set VmuN1 0.1
compute_widths --body_decay=3.0025 n1
launch
# U2 0.0016681005372
# M 1.
set mN1 1.
set VmuN1 0.04084238652674522
compute_widths --body_decay=3.0025 n1
launch
# U2 0.0002782559402
# M 1.
set mN1 1.
set VmuN1 0.0166810053720
compute_widths --body_decay=3.0025 n1
launch
# U2 0.0000464158883
# M 1.
set mN1 1.
set VmuN1 0.0068129206905
compute_widths --body_decay=3.0025 n1
launch
# U2 7.7426368268112
# M 1.
set mN1 1.
set VmuN1 0.0027825594022
compute_widths --body_decay=3.0025 n1
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Jan Hajer
- Solved:
- Last query:
- Last reply: