HT-sliced samples when using FxFx
Dear experts,
I am using MG5_aMC@NLO v2.6.5 and I am studying the possibility of slicing a FxFx Z+jets sample (up to 2 additional patrons in the final state) using the HT variable.
Below you can find the lines I added to the SubProcesses/cuts.f file in order to compute the HT variable and make different HT slices:
right after these lines:
c Apply the jet cuts
if (njet .ne. nQCD .and. njet .ne. nQCD-1) then
return
endif
endif
122 continue
I added:
c
c HT cuts > 1 TeV
c
htj=0d0
do i=1,nexternal
if (istatus(i).eq.1 .and. is_a_j(i)) then ! massless quark/gluon
htj = htj + pt_04(p(0,i))
endif
enddo
if (htj.lt.
return
endif
(this is the example to create the slice 800 < HT < 1500 GeV).
Unfortunately, the cross section I got is identical to the one for the inclusive sample Z+2jets FxFx.
Could you please tell me if what I am trying to do is possible and, if so, could you please advice me on how to do it properly?
Thanks a lot for your time and help.
Best,
Fra
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 Francesco Giuli for more information if necessary.