DELTAR(b[1],b[2]) returns wrong results in the normal mode

Hi Lily,

I do not find any bug. My calculation of the deltaR also gives the 5.36 number for the first event. Note that the deltaphi has to be taken between -pi and pi.

Cheers,

Benjamin

On 08 Nov 2015, at 12:26 , lily <email address hidden> wrote:

> New question #273881 on MadAnalysis 5:
> https://answers.launchpad.net/madanalysis5/+question/273881
>
> Hi,
> Using the user-friendly fashion of Madananlysis5 (the latest version), I find that the command DELTAR(b[1],b[2]) and plot DELTAR(b[1],b[2]) returns the wrong results.
>
>> plot DELTAR(b[1],b[2]) 50 5 6
>> select DELTAR(b[1],b[2]) > 5.5
>
> The lhco. file contain 4 events.
> According the definition of "DELTAR = ((b1["eta"] - b2["eta"])**2 + (b1["phi"] - b2["phi"])**2)**0.5, all the values of "DELTAR(b[1],b[2])" for the 4 events are larger than 5.5. However, "plot DELTAR(b[1],b[2]) 50 5 6" shows a event with DELTAR(b[1],b[2]) = 5.36. "select DELTAR(b[1],b[2]) > 5.5" reject 4 events, and no event is kept.
>
> The reasons maybe from:
>
> 1. The definition of "DELTAR" in ma5 is different from the definition given here?
> 2. There is something wrong with the b-jets ordering of transverse momentum?
>
> The lhco. file is:
> # Event number: 0
> # number trig word
> 0 602 0
> # typ eta phi pt jmass ntrk btag had/em dummy dummy
> 1 1 -0.884 -0.126 120.83 0.0 1.0 0.0 0.0 0.0 0.0
> 2 1 -1.075 -3.018 17.03 0.0 -1.0 0.0 0.0 0.0 0.0
> 3 4 -2.829 2.25 81.0 10.56 3.0 1.0 4.68 0.0 0.0
> 4 4 2.033 -1.742 77.27 11.88 17.0 1.0 1.26 0.0 0.0
> 5 4 -1.846 2.097 73.05 18.62 13.0 1.0 2.51 0.0 0.0
> 6 4 1.993 2.915 44.2 6.33 8.0 1.0 1.65 0.0 0.0
> 7 4 0.611 -1.12 37.15 8.57 8.0 0.0 3.47 0.0 0.0
> 8 4 1.184 -1.202 32.25 4.96 5.0 0.0 2.23 0.0 0.0
> 9 6 0.0 1.676 23.67 0.0 0.0 0.0 0.0 0.0 0.0
> # Event number: 1
> # number trig word
> 0 1610 0
> # typ eta phi pt jmass ntrk btag had/em dummy dummy
> 1 1 1.303 -0.857 288.77 0.0 1.0 0.0 0.0 0.0 0.0
> 2 1 1.553 -0.216 59.8 0.0 -1.0 0.0 0.0 0.0 0.0
> 3 4 -0.114 2.631 188.5 22.26 13.0 1.0 0.12 0.0 0.0
> 4 4 0.766 -2.924 45.79 6.31 4.0 1.0 1.08 0.0 0.0
> 5 4 0.497 1.851 27.38 7.56 6.0 0.0 1.35 0.0 0.0
> 6 4 3.209 -0.58 26.9 5.84 0.0 0.0 6.27 0.0 0.0
> 7 4 1.017 0.243 23.97 4.42 3.0 1.0 0.08 0.0 0.0
> 8 6 0.0 2.149 182.16 0.0 0.0 0.0 0.0 0.0 0.0
> # Event number: 2
> # number trig word
> 0 4013 0
> # typ eta phi pt jmass ntrk btag had/em dummy dummy
> 1 1 0.072 -0.329 56.9 0.0 1.0 0.0 0.0 0.0 0.0
> 2 1 0.684 -2.791 39.67 0.0 -1.0 0.0 0.0 0.0 0.0
> 3 4 2.239 2.451 119.67 11.9 7.0 1.0 0.39 0.0 0.0
> 4 4 -1.466 1.808 40.24 8.24 3.0 0.0 4.18 0.0 0.0
> 5 4 -0.531 -2.559 31.23 10.17 6.0 1.0 0.65 0.0 0.0
> 6 4 1.961 -2.351 27.01 5.24 6.0 0.0 0.57 0.0 0.0
> 7 4 -3.97 -2.58 25.44 2.31 0.0 1.0 0.98 0.0 0.0
> 8 4 -0.198 0.748 22.86 4.14 4.0 1.0 1.29 0.0 0.0
> 9 6 0.0 -0.438 137.99 0.0 0.0 0.0 0.0 0.0 0.0
> # Event number: 3
> # number trig word
> 0 9653 0
> # typ eta phi pt jmass ntrk btag had/em dummy dummy
> 1 1 0.701 -2.551 64.26 0.0 -1.0 0.0 0.0 0.0 0.0
> 2 1 -0.096 0.145 28.09 0.0 1.0 0.0 0.0 0.0 0.0
> 3 4 -0.732 1.798 46.6 8.66 15.0 0.0 2.52 0.0 0.0
> 4 4 -1.042 -1.449 45.09 4.44 6.0 0.0 81.12 0.0 0.0
> 5 4 -0.502 1.133 41.37 8.18 16.0 0.0 0.03 0.0 0.0
> 6 4 -0.249 2.526 34.3 6.58 7.0 1.0 0.21 0.0 0.0
> 7 4 -2.505 -3.023 26.86 8.02 6.0 1.0 2.73 0.0 0.0
> 8 4 -2.217 2.526 25.11 6.63 3.0 1.0 1.13 0.0 0.0
> 9 4 -1.041 0.58 23.69 5.31 6.0 1.0 0.59 0.0 0.0
> 10 4 -2.105 1.361 21.38 3.57 7.0 0.0 0.65 0.0 0.0
> 11 6 0.0 -0.926 96.7 0.0 0.0 0.0 0.0 0.0 0.0
>
> --
> You received this question notification because you are an answer
> contact for MadAnalysis 5.