MadGraph5_aMC@NLO

KeyError when running MadSpin in MG5_v2.3.0

Asked by Peyton Rose

Dear experts,

I am trying to generate p p > z at NLO and decay z > b b~ in the events using MadSpin (*note that the eventual goal is to do a more complicated process that contains z bosons in the final state, but I have been using this as a test for the faster iteration time).

The process is generated via the following commands:
-------------------------------------------------------
> import model loop_sm-no_b_mass
> generate p p > z [QCD]
> output ppTOz_QCD
-------------------------------------------------------

Since I have specified the b-quarks to be massless during generation, I use the following commands in my madspin_card.dat (otherwise it complains about the b-quarks being massless):
------------------------------------------------------------------------------------------------
> import model loop_sm < /path/to/param_card_heavyb.dat > --bypass_check
> decay z > b b~
> launch
------------------------------------------------------------------------------------------------

"param_card_heavyb.dat" is a new param card where I have added the b-quark mass under "Block mass" and under "INFORMATION FOR YUKAWA" using the identical syntax as for the other particles in those subsections and a mass of 4.7e+00. I have also removed the b-quark from the list of "dependent parameters"

The event generation proceeds fine, but MadSpin fails on the first event:
-----------------------------------------------------------------------------------------------------------------------
INFO: Running MadSpin
INFO: This functionality allows for the decay of resonances
INFO: in a .lhe file, keeping track of the spin correlation effets.
INFO: BE AWARE OF THE CURRENT LIMITATIONS:
INFO: (1) Only a succession of 2 body decay are currently allowed
************************************************************
* *
* W E L C O M E to M A D S P I N *
* *
************************************************************
INFO: Extracting the banner ...
INFO: process: p p > z
INFO: options:
INFO: detected model: loop_sm-no_b_mass. Loading...
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
import model loop_sm /export/share/data/prose/Run2/Hbb/eventGeneration_forOfficialProduction_mg523/MG5_aMC_v2_3_0/ppTOz_QCD_madspin/Cards/param_card_heavyb.dat --bypass_check
INFO: detected model: loop_sm. Loading...
INFO: Restrict model loop_sm with file ../models/loop_sm/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
Defined multiparticle all = g gh gh~ d u s c d~ u~ s~ c~ a ve vm vt e- m- ve~ vm~ vt~ e+ m+ b t b~ t~ z w+ h w- tt- tt+
decay z > b b~
launch
INFO: Will use seed 157766244
INFO: generating Madspin matrix element
WARNING: set the mass of the b-quark to its value in the param_card.dat: 4.7 GeV
INFO: We need to recalculate the branching fractions for z
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model /export/share/data/prose/Run2/Hbb/eventGeneration_forOfficialProduction_mg523/MG5_aMC_v2_3_0/models/loop_sm with file ../models/loop_sm/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO:
INFO: decay channels for z : ( width = 2.441735 GeV )
INFO: BR d1 d2
INFO: 1.523609e-01 d~ d
INFO: 1.523609e-01 s~ s
INFO: 1.507412e-01 b~ b
INFO: 1.188172e-01 u~ u
INFO: 1.188172e-01 c~ c
INFO: 6.793939e-02 ve~ ve
INFO: 6.793939e-02 vm~ vm
INFO: 6.793939e-02 vt~ vt
INFO: 3.438743e-02 e+ e-
INFO: 3.438743e-02 m+ m-
INFO: 3.430962e-02 tt+ tt-
INFO:
INFO: generating the production square matrix element
INFO: generate p p > z --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > z pert_QCD --no_warning=duplicate;
INFO: Done 11.55
INFO: generating the full square matrix element (with decay)
INFO: generate p p > z , z > b b~ QCD=99 --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > z pert_QCD, z > b b~ QCD=99 --no_warning=duplicate;
INFO: Done 11.53
INFO: generate matrix element for decay only (1 - > N).
INFO: output standalone_msF /export/share/data/prose/Run2/Hbb/eventGeneration_forOfficialProduction_mg523/MG5_aMC_v2_3_0/ppTOz_QCD_madspin/decay_me
INFO: Done 6.908
INFO: Compiling code

Kept definitions of multiparticles l- / j / vl / l+ / p / vl~ unchanged
Defined multiparticle all = g gh gh~ d u s c d~ u~ s~ c~ a ve vm vt e- m- ve~ vm~ vt~ e+ m+ b t b~ t~ z w+ h w- tt- tt+
decay z > b b~
launch
INFO: Will use seed 157766244
INFO: generating Madspin matrix element
WARNING: set the mass of the b-quark to its value in the param_card.dat: 4.7 GeV
INFO: We need to recalculate the branching fractions for z
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model /export/share/data/prose/Run2/Hbb/eventGeneration_forOfficialProduction_mg523/MG5_aMC_v2_3_0/models/loop_sm with file ../models/loop_sm/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
INFO:
INFO: decay channels for z : ( width = 2.441735 GeV )
INFO: BR d1 d2
INFO: 1.523609e-01 d~ d
INFO: 1.523609e-01 s~ s
INFO: 1.507412e-01 b~ b
INFO: 1.188172e-01 u~ u
INFO: 1.188172e-01 c~ c
INFO: 6.793939e-02 ve~ ve
INFO: 6.793939e-02 vm~ vm
INFO: 6.793939e-02 vt~ vt
INFO: 3.438743e-02 e+ e-
INFO: 3.438743e-02 m+ m-
INFO: 3.430962e-02 tt+ tt-
INFO:
INFO: generating the production square matrix element
INFO: generate p p > z --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > z pert_QCD --no_warning=duplicate;
INFO: Done 11.55
INFO: generating the full square matrix element (with decay)
INFO: generate p p > z , z > b b~ QCD=99 --no_warning=duplicate;define pert_QCD = -4 -3 -2 -1 1 2 3 4 21;add process p p > z pert_QCD, z > b b~ QCD=99 --no_warning=duplicate;
INFO: Done 11.53
INFO: generate matrix element for decay only (1 - > N).
INFO: output standalone_msF /export/share/data/prose/Run2/Hbb/eventGeneration_forOfficialProduction_mg523/MG5_aMC_v2_3_0/ppTOz_QCD_madspin/decay_me
INFO: Done 6.908
INFO: Compiling code
INFO: running MadSpin
INFO: detect independant decays
INFO: MadSpin: Estimate the maximum weight
INFO:
INFO: Estimating the maximum weight
INFO: *****************************
INFO: Probing the first 75 events
INFO: with 400 phase space points
INFO:
INFO: Event 1/75 : 0.17s
Command "launch -x -o" interrupted with error:
KeyError : ((-5, 5), (23,))
Please report this bug on https://bugs.launchpad.net/madgraph5
More information is found in '/export/share/data/prose/Run2/Hbb/eventGeneration_forOfficialProduction_mg523/MG5_aMC_v2_3_0/ppTOz_QCD_madspin/run_16_grid_setup_debug.log'.
Please attach this file to your report.

---------------------------------------------------------------------------------------

Do you know what's going wrong here, or have suggestions for how to fix this? Thanks in advance for your help.

Regards,
Peyton

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
Pierre Artoisenet Edit question
Solved by:
Olivier Mattelaer
Solved:
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#1

Hi,

>Since I have specified the b-quarks to be massless during generation, I use the following commands in my madspin_card.dat (otherwise it complains about the b-quarks being massless):

I do not see which kind of complain you refer to. I thought that the code was handling that automatically.
I assign this to Pierre since he plays more than me with the mass of the quarks in MadSpin.

> Command "launch -x -o" interrupted with error:
KeyError : ((-5, 5), (23,))

This error is clearly linked to the fact that in your original model your b were massless and therefore include in the proton definition.
When you change your model to massive, your b are not part of the proton anymore. Therefore when we generate the new matrix element. We do not have any b in the initial state and as soon as you have on event with a b in the initial state it crashes.

The easiest is to do everything with massive b (this is at least fine for your test). (or everything with massless but if Pierre contradict me on this)
If this is not satisfactory, then I have an idea but it would require testing.

Cheers,

Olivier

Revision history for this message
Peyton Rose (pwrose21) said : 		#2

Hi Olivier, Pierre,

Thanks for your quick reply. Just to clarify:

>>Since I have specified the b-quarks to be massless during generation, I use the following commands in my madspin_card.dat (otherwise it complains about the b-quarks being massless):

>I do not see which kind of complain you refer to. I thought that the code was handling that automatically.

The error I am referring to here is in MadSpin. If I use the param card with massless b-quarks, I get the following error:
-----------------------------------------------------------------------------------------------------
INFO: Running MadSpin
INFO: This functionality allows for the decay of resonances
INFO: in a .lhe file, keeping track of the spin correlation effets.
INFO: BE AWARE OF THE CURRENT LIMITATIONS:
INFO: (1) Only a succession of 2 body decay are currently allowed
************************************************************
* *
* W E L C O M E to M A D S P I N *
* *
************************************************************
INFO: Extracting the banner ...
INFO: process: p p > z
INFO: options:
INFO: detected model: loop_sm-no_b_mass. Loading...
set max_weight_ps_point 400 # number of PS to estimate the maximum for each event
decay z > b b~
launch
INFO: Will use seed 619431129
INFO: generating Madspin matrix element
INFO: We need to recalculate the branching fractions for z
INFO: Using MadWidth (arXiv:1402.1178)
INFO: Restrict model /export/share/data/prose/Run2/Hbb/eventGeneration_forOfficialProduction_mg523/MG5_aMC_v2_3_0/models/loop_sm with file ../models/loop_sm/restrict_default.dat .
INFO: Run "set stdout_level DEBUG" before import for more information.
Command "launch -x -o" interrupted with error:
MadGraph5Error : yukawa (5,) not define
------------------------------------------------------------------------------------------------------------

After searching through some older threads where a similar problem was solved, I was led me use a new param card with massive b-quarks. However, the older thread was dealing with decays to massive/massless leptons, so there was no conflict in the proton definition.

Regards,
Peyton

Revision history for this message
Best Olivier Mattelaer (olivier-mattelaer) said : 		#3

Dear Peyton,

Ok I see, this is actually a bug which is fixed in 2.3.2 which is going to be release sometimes today.

Cheers,

Olivier

Revision history for this message
Peyton Rose (pwrose21) said : 		#4

Thanks Olivier Mattelaer, that solved my question.