Forced semi-leptonic decay of b-hadrons
Dear MG5 Team members,
I'm facing a problem with the MG5-PY8 framework while generating heavy flavor pair production process, especially p p > b b~. I want to study di-leptonic events from this process only and want to switch ON the leptonic modes of the b-hadrons to generate events efficiently. To do so, I write the following lines at the end of the PY8 card in MG5 Cards directory,
!======
511:onMode = off
511:onIfAny = 11 13 15 -11 -13 -15
521:onMode = off
521:onIfAny = 11 13 15 -11 -13 -15
531:onMode = off
531:onIfAny = 11 13 15 -11 -13 -15
541:onMode = off
541:onIfAny = 11 13 15 -11 -13 -15
!======
Now, after PY8 run is over, the log and cmnd files print the following
!======
541:onMode=off
541:onIfAny=11 13 15 -11 -13 -15
521:onIfAny=11 13 15 -11 -13 -15
521:onMode=off
531:onIfAny=11 13 15 -11 -13 -15
531:onMode=off
511:onMode=off
511:onIfAny=11 13 15 -11 -13 -15
!======
i.e, the order of decay statements is changed while the order is very important for such forced decay. This happens almost for any given list of b-hadrons; some are correctly printed while the order is changed for others. Of course, it changes the expected final state lepton multiplicity significantly. For example, the decay of particles with ID 521 and 531 are OFF completely (I have checked the Pythia log file too, it is indeed a stable particle) in stead of leptonic decays.
It would be nice if you kindly provide some help to solve this.
Thank you in advance.
Best regards,
Amit Chakraborty
Question information
- Language:
- English Edit question
- Status:
- Open
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Amit Chakraborty for more information if necessary.