21 |
|
return sqrt(dEta*dEta+dPhi*dPhi); |
22 |
|
} |
23 |
|
|
24 |
+ |
double DeltaR(const BaseVars & v1, const GenVars & v2) |
25 |
+ |
{ |
26 |
+ |
|
27 |
+ |
double dEta = v1.eta - v2.eta; |
28 |
+ |
double dPhi = v1.phi - v2.phi; |
29 |
+ |
if (dPhi<0) dPhi += 2*TMath::Pi(); |
30 |
+ |
|
31 |
+ |
return sqrt(dEta*dEta+dPhi*dPhi); |
32 |
+ |
} |
33 |
+ |
|
34 |
|
double SameSign(const BaseVars & v1, const BaseVars & v2) |
35 |
|
{ |
36 |
|
|
111 |
|
double py1 = leg1.pT*sin(leg1.phi); |
112 |
|
double py2 = leg2.pT*sin(leg2.phi); |
113 |
|
|
104 |
– |
|
114 |
|
double x1_numerator = px1*py2 - px2*py1; |
115 |
|
double x1_denominator = py2*(px1 + mEx) - px2*(py1 + mEy); |
116 |
|
double x1 = ( x1_denominator != 0. ) ? x1_numerator/x1_denominator : -1.; |