KeyError : 'id (999000006,) is not in mass' when calculating relic density

[Inga S]

I have installed maddm, and the internal installations were successful. I do the following
MadDM>import model MSSM_SLHA2
MadDM>generate relic_density
MadDM>add direct_detection
MadDM>add indirect_detection
MadDM>output RELIC_PROC

Then I copy the SLHA file I want to use to RELIC_PROC/Cards/param_card.dat and do

MadDM>launch

whereupon I get the error message

Command "launch" interrupted with error:
KeyError : 'id (999000006,) is not in mass'
Please report this bug on https://bugs.launchpad.net/mg5amcnlo
More information is found in 'MG5_debug'.
Please attach this file to your report.

Have I misunderstood the workflow?
All I want to do is to calculate the neutralino relic density from an SLHA file where the neutralino is the LSP.
Thank you!

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

Thanks for the report, this is indeed one of the problem of our solution if you use a pre-existing SLHA file.
The problem is that we update your model automatically to be able to do direct detection.
In that process, we unfortunatly add new particles (that atually does not contribute to anything) and you need those to be defined (mass/width) at run time in the param_card.

We will see how we can fix this issue in a future version (of either maddm or mg5amc),
but for the moment a (dirty) solution is to type "update missing" after having defined your SLHA:
(you can see below that it automatically fix the card by adding the additional required --but useless-- line)

Thanks for the valuable information.

Olivier

The following switches determine which programs are run:
/=========== Description ============|=============== values ===============|======== other options =========\
| 1. Compute the Relic Density | relic = ON | OFF |
| 2. Compute direct(ional) detection | direct = directional | direct|OFF |
| 3. Compute indirect detection/flux | indirect = sigmav | flux_source|flux_earth|OFF |
| 4. Run Multinest scan | nestscan = OFF | Please install module |
\============================================================================================================/
 You can also edit the various input card:
 * Enter the name/number to open the editor
 * Enter a path to a file to replace the card
 * Enter set NAME value to change any parameter to the requested value
 /=============================================================================\
 | 5. Edit the model parameters [param] |
 | 6. Edit the MadDM options [maddm] |
 \=============================================================================/

>update missing
INFO: Adding 4 parameter(s) to block mass
INFO: Adding 5 parameter(s) to block decay
INFO: write new param_card with 9 new parameter(s).
The following switches determine which programs are run:
/=========== Description ============|=============== values ===============|======== other options =========\
| 1. Compute the Relic Density | relic = ON | OFF |
| 2. Compute direct(ional) detection | direct = directional | direct|OFF |
| 3. Compute indirect detection/flux | indirect = sigmav | flux_source|flux_earth|OFF |
| 4. Run Multinest scan | nestscan = OFF | Please install module |
\============================================================================================================/
 You can also edit the various input card:
 * Enter the name/number to open the editor
 * Enter a path to a file to replace the card
 * Enter set NAME value to change any parameter to the requested value
 /=============================================================================\
 | 5. Edit the model parameters [param] ...

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

Are you sure that you are using a slha2 format?
From the update missing command it seems that you are using a slha1 format
Can you check that?

In principle we have a command to convert slha1 card to slha2 card but I'm not sure that it will work in this context
(due to the additional mass).

Cheers,

Olivier

> On 27 Oct 2018, at 09:35, Inga S <email address hidden> wrote:
>
> Public bug reported:
>
> I have installed maddm, and the internal installations were successful. I do the following
> MadDM>import model MSSM_SLHA2
> MadDM>generate relic_density
> MadDM>add direct_detection
> MadDM>add indirect_detection
> MadDM>output RELIC_PROC
>
> Then I copy the SLHA file I want to use to
> RELIC_PROC/Cards/param_card.dat and do
>
> MadDM>launch
>
> whereupon I get the error message
>
> Command "launch" interrupted with error:
> KeyError : 'id (999000006,) is not in mass'
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in 'MG5_debug'.
> Please attach this file to your report.
>
> Have I misunderstood the workflow?
> All I want to do is to calculate the neutralino relic density from an SLHA file where the neutralino is the LSP.
> Thank you!
>
> ** Affects: mg5amcnlo
> Importance: Undecided
> Status: New
>
> --
> You received this bug notification because you are subscribed to
> MadGraph5_aMC@NLO.
> https://bugs.launchpad.net/bugs/1800265
>
> Title:
> KeyError : 'id (999000006,) is not in mass' when calculating relic
> density
>
> Status in MadGraph5_aMC@NLO:
> New
>
> Bug description:
> I have installed maddm, and the internal installations were successful. I do the following
> MadDM>import model MSSM_SLHA2
> MadDM>generate relic_density
> MadDM>add direct_detection
> MadDM>add indirect_detection
> MadDM>output RELIC_PROC
>
> Then I copy the SLHA file I want to use to
> RELIC_PROC/Cards/param_card.dat and do
>
> MadDM>launch
>
> whereupon I get the error message
>
> Command "launch" interrupted with error:
> KeyError : 'id (999000006,) is not in mass'
> Please report this bug on https://bugs.launchpad.net/mg5amcnlo
> More information is found in 'MG5_debug'.
> Please attach this file to your report.
>
> Have I misunderstood the workflow?
> All I want to do is to calculate the neutralino relic density from an SLHA file where the neutralino is the LSP.
> Thank you!
>
> To manage notifications about this bug go to:
> https://bugs.launchpad.net/mg5amcnlo/+bug/1800265/+subscriptions

[Olivier Mattelaer (olivier-mattelaer)]

Hi,

This is for sure a SLHA1 format that you are using. Actually your output is INCOMPLETE.
Since a lot of SM information are missing.

Indeed If I use the SLHA1 model, and that I type
update missing, I have the following:

INFO: unknow block SPINFO: block will be ignored
INFO: unknow block MINPAR: block will be ignored
INFO: unknow block EXTPAR: block will be ignored
INFO: Adding 15 parameter(s) to block MASS
INFO: unknow block gauge: block will be ignored
INFO: Adding 16 parameter(s) to block decay

If i do the same from the SLHA2 model:

INFO: unknow block SPINFO: block will be ignored
INFO: unknow block MODSEL: block will be ignored
INFO: unknow block MINPAR: block will be ignored
INFO: unknow block EXTPAR: block will be ignored
INFO: Adding 15 parameter(s) to block MASS
INFO: unknow block alpha: block will be ignored
INFO: unknow block stopmix: block will be ignored
INFO: unknow block sbotmix: block will be ignored
INFO: unknow block staumix: block will be ignored
INFO: unknow block gauge: block will be ignored
INFO: unknow block au: block will be ignored
INFO: unknow block ad: block will be ignored
INFO: unknow block ae: block will be ignored
INFO: Block snumix was missing. Adding the 3 associated parameter(s)
INFO: Block msu2 was missing. Adding the 3 associated parameter(s)
INFO: Block vckm was missing. Adding the 3 associated parameter(s)
INFO: Block fralpha was missing. Adding the 1 associated parameter(s)
INFO: Block upmns was missing. Adding the 3 associated parameter(s)
INFO: Block msd2 was missing. Adding the 3 associated parameter(s)
INFO: Block msl2 was missing. Adding the 3 associated parameter(s)
INFO: Block tu was missing. Adding the 1 associated parameter(s)
INFO: Block selmix was missing. Adding the 8 associated parameter(s)
INFO: Block td was missing. Adding the 1 associated parameter(s)
INFO: Block te was missing. Adding the 1 associated parameter(s)
INFO: Block usqmix was missing. Adding the 8 associated parameter(s)
INFO: Block dsqmix was missing. Adding the 8 associated parameter(s)
INFO: Block mse2 was missing. Adding the 3 associated parameter(s)
INFO: Block msq2 was missing. Adding the 3 associated parameter(s)
INFO: Adding 16 parameter(s) to block decay
INFO: write new param_card with 83 new parameter(s).

But using that command from a slha1 file is certainly not correct since now you have a mix of SLHA1 and SLHA2 card.

So you need
1) to use a converter to slha2 ("update to_slha2" seems to do the job but I'm not 100% confident in that converter)
2) update the missing information ("update missing")
As you can see below

Do you want to edit a card (press enter to bypass editing)?
/------------------------------------------------------------\
| 1. param : param_card.dat |
| 2. run : run_card.dat |
| 3. plot : plot_card.dat |
\------------------------------------------------------------/
 you can also
   - enter the path to a valid card or banner.
   - use the 'set' command to modify a parameter directly.
     The set option works only for param_card and run_card.
     Type 'help set' for more information ...