ptl1min cut not working

Asked by rojalin padhan on 2018-11-18

hiii..

       I'm using the cuts ptl1min . I set ptl1min to 100. so with this cut we expect that PT of leading lepton to be greater than 100.
But it seems to be not working, because the histogram for leading lepton (PT(l+[1]) ) is not starting from 100. It is only taking the default PT cut.
 I want to generate event with PT of leading lepton > 100. So could you please help me.

Here is my script:
import model SM_plus_Type_II_Seesaw_UFO
define pb = p b b~
generate pb pb > z z , z > l+ l-, z > l+ l-
output bkg
launch bkg
set vevT 1e-8
set lam1 1
set lam3 1
set lam4 1
set lam 0.52
set muT 8e-8
set width 9900041 Auto
set cut_decays=on
set ptl1min 200.0

thanks and regards,
rojalin

Question information

Language:
English Edit question
Status:
Solved
For:
MadGraph5_aMC@NLO Edit question
Assignee:
No assignee Edit question
Solved by:
Olivier Mattelaer
Solved:
2018-12-06
Last query:
2018-12-06
Last reply:
2018-12-06
Launchpad Janitor (janitor) said : #1

This question was expired because it remained in the 'Open' state without activity for the last 15 days.

Hi,

Sorry for the long time for this answer. I have missed flag it. I guess that you found the reason for this issue right now but just in case.

So in short everything is fine and behaves as it should.

The apparent issue is that you are plotting the pt of the positively charged lepton (l+)
and that madgraph applies that cut on all lepton (positively and negatively charged).

Therefore when you have events where a l- has a pt larger than 100, they are no cut anymore on the l+ and in that case both can be lower than 100GeV.

Cheers,

Olivier

rojalin padhan (rossy.mg) said : #3

hii..

Thanks for your answer.

 I have one question here,

Does PT(l+[1]) mean PT of leading plus sign lepton in an event ?

regards,

rojalin

Hi,

I'm not a MadAnalysis5 author. So you should consider my answer accordingly.

>Does PT(l+[1]) mean PT of leading plus sign lepton in an event ?

I would say yes (otherwise what would be the difference compare to PT(l-[1]).
As far as i know, it should also be possible to define your own class of particles within MadAnalysis5.

Cheers,

Olivier

rojalin padhan (rossy.mg) said : #5

Thanks Olivier Mattelaer, that solved my question.