MadAnalysis 5

lepton PT cut seems didn,t work

Asked by lukeskywalker

ma5>import /home/luke/MG5_aMC_v2_7_31/3t1b/Events/run_01/tag_1_delphes_events.root
ma5>set main.stacking_method = normalize2one
ma5>plot N(b)
ma5>plot N(l)
ma5>select N(b)>=3
ma5>select N(l)>=2
ma5>plot N(b)
ma5>plot N(l)
ma5>plot PT(b)
ma5>plot PT(l)
ma5>plot PT(b) 100 0 100
ma5>plot PT(l) 100 0 100
ma5>select PT(b)>20
ma5>select PT(l)>20
ma5>plot PT(b)
ma5>plot PT(l)
ma5>plot PT(b) 100 0 100
ma5>plot PT(l) 100 0 100
ma5>submit c4

This is my flow of cuts and plots,but after select PT(l)>20 ,the PT plot for leptons generated by"plot PT(l) 100 0 100"seems doesn't change. It still have parts of PT<20. But after I use reject PT(l)<20 , it works.
So I am confused is there something difference between select PT(l)>20 and reject PT(l)<20 .
Best Regards.
LUKE.

Question information

Language:
English Edit question
Status:
Solved
For:
MadAnalysis 5 Edit question
Assignee:
No assignee Edit question
Solved by:
Benjamin Fuks
Solved:
Last query:
Last reply:
Revision history for this message
Best Benjamin Fuks (fuks) said : 		#1

Hi Luke,

The cut ```select PT(l)>20``` selects those events for which at least one lepton has a pT larger than 20 GeV, whereas the cut ```reject PT(l)<20``` rejects those events for which at last one lepton with a pT smaller than 20 GeV is found. What you observe is related to the fact that there are several leptons in each event.

Maybe the cut you actually want to implement is ```select (l) PT > 20```, that corresponds to keeping in each events only the leptons with a pT larger than 20 GeV (and ignoring any other lepton).

May the force be with you,

Benjamin

Revision history for this message
lukeskywalker (luke19941001) said : 		#2

Thanks Benjamin Fuks, that solved my question.

Revision history for this message
lukeskywalker (luke19941001) said : 		#3

Thanks for your reply,that's helps a lot.
Best Regards.
luke