Crash with Phase space cuts inside MadSpin card (MG5_aMC_v2.6.5)
Dear Experts,
I am trying to generate the gridpacks needed to produce NLO MC samples of Higgs dalitz decay to a pair of muons and a photon. The cards used for the purpose can be found here (in particular, the folders *HToMuMuGamma):
https:/
Unfortunately, I am not the original author of the cards and they were created about 3-4 years ago. These cards were written for MG5_aMC_v2.4.2, but now I have to regenerate the gridpacks using necessarily the MG5_aMC_v2.6.5 version.
However, when running the gridpack generation, the pilotrun fails since it is not possible to apply a 0<m_ll<60 GeV phase space cut inside the madspin card. For instance, the problematic lines of the cards are the following:
https:/
set run_card mmll=0
set run_card mmllmax=60
A long discussion has been done with the gen-contacts of the physics group and it can be found here:
https:/
where, the problem is fully described.
Apparently, from the version 2.4.2 to 2.6.5, a lot of changes have been done and it seems that this feature (cuts in the madspin card) are no longer supported, from what I understood. A clean fix could not be found for this problem for the version 2.6.5 and the only way to perform a correct and bug-free gridpack generation was to remove those lines. By this way, a higher number of events would be required in order to take into account the loss in efficiency.
Apparently, a possible way to solve the problem would be to create a second run_card for madspin, but it seems that this solution may lead to bugs that would cause other problems to the gridpack generation.
Thus, I am writing here to ask if there is a clean and simple way to apply this m_ll cut inside the madspin card or somewhere else (like in the run card), avoiding its implementation in the fortran file 'cuts.f'. For a LO generation, this would be simple to implement inside the run_card, since mmll and mmllmax entries are available, but in the NLO card template I could not find a way to set an upper cut for m_ll.
I apologize for the long question and I hope that the description is clear.
Thanks in advance!
Best,
Rocco
Question information
- Language:
- English Edit question
- Status:
- Solved
- Assignee:
- No assignee Edit question
- Solved by:
- Olivier Mattelaer
- Solved:
- Last query:
- Last reply: