Adding off-shell contribution Scalar Higgs model
Affects | Status | Importance | Assigned to | Milestone | |
---|---|---|---|---|---|
MadDM |
Fix Committed
|
High
|
Daniele Massaro |
Bug Description
I am trying to calculate the contribution to the cross section from the 2 --> 3 annihilation processes via an off-shell W or Z-boson, DMDM --> WW* and ZZ*.
I am doing this for the ScalarHiggsPortal model.
Following the suggestion from this paper of MadDM https:/
###############
python2.7 bin/maddm.py
MadDM > import model ScalarHiggsPort
MadDM > define darkmatter n1
MadDM > generate relic_density
MadDM > define ferm = ve vm vt ve~ vm~ vt~ u c t d s b u~ c~ t~ d~ s~ b~ e- mu- ta- e+ mu+ ta+
MadDM > define wz = w+ w- z
MadDM > add indirect_detection wz ferm ferm
MadDM > output test
MadDM > launch test
###############
The diagrams indeed contain the off-shell production of Z* and W*.
However, during the calculation of the sigmav, I obtain this error:
###############
INFO: generating the square matrix element for reweighting
INFO: generate n1 n1 > w+ b t~ ..... ;add process n1 n1 > w- c s~ SIEFFV=0 SIEFFS=0 SIEFFF=0 SDEFFV=0 SDEFFF=0 ;
INFO: Done 5.209
WARNING: REWEIGHTING: original card and new card are identical.
INFO: starts to compute weight for events with the following modification to the param_card:
INFO:
INFO: Event nb 0 current time: 13h42
.........
INFO: Event nb 9000 8.7s
INFO: All event done (nb_event: 10000) 9.7s
INFO: Eventfiles is/are now created with new central weight
INFO: new cross-section is : 3.76501e+07 pb (indicative error: 124739 pb)
INFO: Original cross-section: 37649450.0 +- 120937.7 pb (cross-section from sum of weights: 37649453.0)
INFO: Computed cross-section:
INFO: rwgt_1 : 37650051.4121 +- 124738.969479 pb
quit
Command "launch test" interrupted with error:
KeyError : 'n1n1_wpvltam'
###############
I send in attachment the zip of the Card folder with also the Debug file.
A similar error occurred for the model import model DMsimp_s_spin1. See this bug:
https:/
Thanks for the consideration!
Related branches
Changed in maddm: | |
importance: | Undecided → High |
status: | New → Triaged |
assignee: | nobody → Daniele Massaro (dmassaro) |
Changed in maddm: | |
status: | Triaged → Fix Committed |
Hi Mattia,
I committed a fix in the lp:~maddm/maddm/process_keys which you can download.
The main issue is that MadGraph automatically groups different subprocesses when it is faster to compute them together. So they appear like dmdm_wmvll or dmdm_zvlvl. So far it is not possible to translate back the key name to one of a generated process.
Cheers.
--
Daniele