Cut Decays on HepMC
Dear MadGraph, I have a question regarding cut_decays. I am trying to generate events that will automatically pass a 27 GeV lepton trigger. For this, I run:
generate p p > Z H, (Z > l+ l-)
On the run card, I set:
27.0 = xptl ! minimum pt for at least one charged lepton
Of course, I set cut_decays to True. If I inspect the generated lhe file, everything is ok. However, after running Pythia, the leptons on the HepMC output do not necessarily pass the cuts. For instance, one of the events on the lhe has an electron that passes the cut:
11 1 3 3 0 0 -2.8638742914e+01 -1.6432462349e+00 +1.2197184318e+02 1.2529967431e+02 0.0000000000e+00 0.0000e+00 -1.0000e+00
However, on the HepMC this electron will end up with lower pT (it appears in two "vertices", one of which has only one particle):
P 901 11 -2.500547280160
P 904 11 -2.500547055315
As you can see, the pT is reduced from 28.6 GeV to 25.2 GeV. This makes me lose about 2% of the events, which I guess I can live with. A similar thing happens if I run t t~ H production. Still, first request, could you shed light on why this happens?
The second question is much more worrying. I now run:
generate p p > W+ H, (W+ > l+ nu)
add process p p > W- H, (W- > l- nu)
Here, MadGraph does not give me the option for setting xptl, so I set:
27.0 = ptl ! minimum pt for the charged leptons
In this case, I end up losing about 25% of the generated events! Here I have again traced an example. On the lhe, everything is good, look at this electron:
11 1 3 3 0 0 -2.3672448034e+01 -4.7028479844e+01 -8.5434356876e+00 5.3339038290e+01 0.0000000000e+00 0.0000e+00 -1.0000e+00
On the HepMC, one finds that a lot of energy is lost by photon emission. Let me copy the two relevant vertices:
V -306 0 -4.411290175414
P 399 11 -2.314560986254
P 400 12 -1.836148988103
V -307 0 -4.411290175414
P 401 11 -1.175117303160
P 402 22 -1.190031601169
Here, one sees that the photon takes 33.4 GeV as pT, leaving the original electron with 17.9 GeV (so it does not pass the cut). I understand that these radiative processes are to be expected, but what I do not understand is why they do not happen (or do not happen so frequently) for the other processes I generate?
Many thanks!
Question information
- Language:
- English Edit question
- Status:
- Answered
- Assignee:
- No assignee Edit question
- Last query:
- Last reply:
Can you help with this problem?
Provide an answer of your own, or ask Joel Jones for more information if necessary.