Exclude syntax could work better (i.e. generate relic / X)
Hi,
I am trying to calculate the velocity averaged cross-section of DM annihilation. I would like to split this calculation by final state, to see what the individual contributions are. As an example let's say I want annihilation into everything except bbar.
Naively, one would expect this to be possible by calling "generate relic / b b~". Unfortunately, this does not work. The reason it does not work is the way madgraph is called in the generate_relic function in maddm_interface.py:
if excluded_particles:
proc = "dm_particles dm_particles > fs_particles fs_particles / %s %s @DM2SM" \
Since the b quarks are part of fs_particles as well as excluded_particles, the slash notation to exclude b quarks is simply ignored by madgraph. E.g. "generate dm dm > b b~ / b b~" does give diagrams (Why MG decides to handle it like that is another question...)
A simple solution to this would be to simply kick the excluded particles out of the fs_particles group, as demonstrated here: https:/
I think this would be a useful addition/fix of functionality to include in MadDM.
Best
Andreas
Question information
- Language:
- English Edit question
- Status:
- Answered
- For:
- MadDM Edit question
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Andreas for more information if necessary.