MadGraph5_aMC@NLO

How to apply "HT cut" for heavy particles?

Asked by Li

Hi,
    I try to generate different pT range for process like pp->ww, however, the pt_min/max_pdg cut applies for both W so not all the pT range is include ( for example one W with pT<100GeV and the other W with pT>100GeV ), I wonder if you have and idea.

    What I think is to modify "htj" from counting jet pt to final particles pt, but I don't know how to modify in the cuts.f,
    I guess the following is the main code,

      htj=ptjet(1)

      do i=2,njets
         htj=htj+ptjet(i)

      if(htj.lt.htjmin.or.htjmax.ge.0d0.and.htj.gt.htjmax)then
         if(debug) write (*,*) i, ' htj -> fails'
         passcuts=.false.
         return

    What can I do to modify the htj to the "ht" for other final state particles like W?

    Thank you!

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Li
Solved:
Last query:
Last reply:
Revision history for this message
Olivier Mattelaer (olivier-mattelaer) said : 		#1

Hi,

Changing existing cut is dangerous since you not only need to change cuts.f but also all the phase-space optimization which are done based on the cut.

So my advise would be to add an additional cut within the function which is designed for that in dummy_fcts.f.

Cheers,

Olivier

Revision history for this message
Li (huangli-itp) said : 		#2

Hi Oliviver,
    Thank you! I tried what you said these days and it works perfect!
    Thank you very much!

Best,
Li