Two types of reject Commands

Asked by Kasinath Das

Hi ma5 team,
     What is the difference between the following to commands
                    (1) reject ABSETA(j)>1.5
                    (1) reject (j)ABSETA>1.5

Thanking You,
Kasinath

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
Benjamin Fuks (fuks) said :
#1

Hi Kasinath,

> (1) reject ABSETA(j)>1.5
This rejects all events for which the condition is satisfied

> (1) reject (j)ABSETA>1.5
This removes all jets for which the condition is satisfied from the jet collection. It does not reject any event.

Cheers,

Benjamin

Revision history for this message
Kasinath Das (kasinath-das91) said :
#2

Hi Benjamin,
        I am getting the following results with the use of the above two commands separately.

        (While I was plotting the histograms for ETA(j) ETA(b) and ETA(nb) after applying the cut.)

Case-1 With reject ABSETA(j)>1.5 [rejetcs all the events where at least one jet satisfies the condition.]

                                  For ETA(j) there is no events for ETA(j)>1.5 and ETA(j) < 1.5
                                  For ETA(b) there is no events for ETA(b)>1.5 and ETA(b) < 1.5
                                   For ETA(nb) there is no events for ETA(nb)>1.5 and ETA(nb) < 1.5

 Case -2 With reject (j)ABSETA>1.5 [removes all the jets from events which satisfies the condition.]

                                   For ETA(j) there is no events for ETA(j)>1.5 and ETA(j) < 1.5

               BUT FOR ETA(b) I am getting events for ETA(b)>1.5 and ETA(b) < 1.5
                                   For ETA(nb) there is also events for ETA(nb)>1.5 and ETA(nb) < 1.5

                      If j contains both bjets(b) and non_bjets(nb) then in case-2 we should not get in any event b and nb jets
                      having absolute value of ETA beyond 1.5.

             Here are the links to the histograms for the above two cases.
             https://www.dropbox.com/s/qa8jdp7w1vlln2l/reject_ABSETA%28j%29.pdf?dl=0
            https://www.dropbox.com/s/96o92dvgk4edsvk/reject_%28j%29ABSETA.pdf?dl=0

Thanking You,
Kasinath

Revision history for this message
Benjamin Fuks (fuks) said :
#3

Hi Kasinath,

The command reject (j) ETA > 1.5 does not yield to the rejection of any event. It takes the jet collection, and removes from it all jets for which the condition is not satisfied. This is the expected behavior and this is what you observe.

Cheers,

Benjamin

On 26 Nov 2016, at 08:53 , Kasinath Das <email address hidden> wrote:

> Question #404390 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/404390
>
> Status: Answered => Open
>
> Kasinath Das is still having a problem:
>
> Hi Benjamin,
> I am getting the following results with the use of the above two commands separately.
>
> (While I was plotting the histograms for ETA(j) ETA(b)
> and ETA(nb) after applying the cut.)
>
>
> Case-1 With reject ABSETA(j)>1.5 [rejetcs all the events where at least one jet satisfies the condition.]
>
> For ETA(j) there is no events for ETA(j)>1.5 and ETA(j) < 1.5
> For ETA(b) there is no events for ETA(b)>1.5 and ETA(b) < 1.5
> For ETA(nb) there is no events for ETA(nb)>1.5 and ETA(nb) < 1.5
>
>
>
>
> Case -2 With reject (j)ABSETA>1.5 [removes all the jets from events which satisfies the condition.]
>
> For ETA(j) there is no events for ETA(j)>1.5 and ETA(j) < 1.5
>
> BUT FOR ETA(b) I am getting events for ETA(b)>1.5 and ETA(b) < 1.5
> For ETA(nb) there is also events for ETA(nb)>1.5 and ETA(nb) < 1.5
>
>
> If j contains both bjets(b) and non_bjets(nb) then in case-2 we should not get in any event b and nb jets
> having absolute value of ETA beyond 1.5.
>
>
> Here are the links to the histograms for the above two cases.
> https://www.dropbox.com/s/qa8jdp7w1vlln2l/reject_ABSETA%28j%29.pdf?dl=0
> https://www.dropbox.com/s/96o92dvgk4edsvk/reject_%28j%29ABSETA.pdf?dl=0
>
>
> Thanking You,
> Kasinath
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Kasinath Das (kasinath-das91) said :
#4

Hi Benjamin,

         Yes, I agree that reject (j) ETA > 1.5 removes jets satisfying the condition from the jet collection.
          We know j contains both bjets(b) and non_bjets(nb). Then after applying this cut we should not get
         any jet (wheather it is b or nb) having eta value beyond 1.5.
                                            The problem that I am facing is this. The ETA(j) plot HAS NO EVENTS with jets having eta value greater than 1.5. But the ETA(b) plot HAS EVENTS with bjets having eta value greater than 1.5. ETA(nb) plot also HAS EVENTS with nb jets having eta value more than 1.5. Is not this contradicting the fact that j contains both b and nb ?

Thanking you,
Kasinath

Revision history for this message
Best Benjamin Fuks (fuks) said :
#5

Hi Kasinath,

Those are actually different (anbd independent) collections: the objects are duplicated. You have only applied cuts on the j collection, and not on the b abd nb collections. This is why you still get b and nb jets with eta > 1.5.

Cheers,

Benj

On 28 Nov 2016, at 15:22 , Kasinath Das <email address hidden> wrote:

> Question #404390 on MadAnalysis 5 changed:
> https://answers.launchpad.net/madanalysis5/+question/404390
>
> Status: Answered => Open
>
> Kasinath Das is still having a problem:
> Hi Benjamin,
>
> Yes, I agree that reject (j) ETA > 1.5 removes jets satisfying the condition from the jet collection.
> We know j contains both bjets(b) and non_bjets(nb). Then after applying this cut we should not get
> any jet (wheather it is b or nb) having eta value beyond 1.5.
> The problem that I am facing is this. The ETA(j) plot HAS NO EVENTS with jets having eta value greater than 1.5. But the ETA(b) plot HAS EVENTS with bjets having eta value greater than 1.5. ETA(nb) plot also HAS EVENTS with nb jets having eta value more than 1.5. Is not this contradicting the fact that j contains both b and nb ?
>
>
> Thanking you,
> Kasinath
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.

Revision history for this message
Kasinath Das (kasinath-das91) said :
#6

Thanks Benjamin Fuks, that solved my question.