6 |
|
{ |
7 |
|
|
8 |
|
double dPhi = phi1 - phi2; |
9 |
– |
if (dPhi<0) dPhi += 2*TMath::Pi(); |
10 |
– |
|
11 |
– |
return dPhi; |
12 |
– |
} |
13 |
– |
|
14 |
– |
double DeltaR(const BaseVars & v1, const BaseVars & v2) |
15 |
– |
{ |
16 |
– |
|
17 |
– |
double dEta = v1.eta - v2.eta; |
18 |
– |
double dPhi = v1.phi - v2.phi; |
19 |
– |
if (dPhi<0) dPhi += 2*TMath::Pi(); |
20 |
– |
|
21 |
– |
return sqrt(dEta*dEta+dPhi*dPhi); |
22 |
– |
} |
23 |
– |
|
24 |
– |
double DeltaR(const BaseVars & v1, const GenVars & v2) |
25 |
– |
{ |
9 |
|
|
10 |
< |
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); |
10 |
> |
return dPhi; |
11 |
|
} |
12 |
|
|
13 |
|
double SameSign(const BaseVars & v1, const BaseVars & v2) |
47 |
|
double et = leg1.pT + leg2.pT + TMath::Sqrt(mEx*mEx + mEy*mEy); |
48 |
|
double mt2 = et*et - (px*px + py*py); |
49 |
|
if ( mt2 < 0 ) { |
50 |
< |
std::cout << " --- WARNING : mt2 = " << mt2 << " is negative... Set to 0."; |
50 |
> |
//std::cout << " --- WARNING : mt2 = " << mt2 << " is negative... Set to 0."; |
51 |
|
return 0.; |
52 |
|
} |
53 |
|
return sqrt(mt2); |
62 |
|
double et = leg1.pT + TMath::Sqrt(mEx*mEx + mEy*mEy); |
63 |
|
double mt = et*et - (px*px + py*py); |
64 |
|
if ( mt < 0 ) { |
65 |
< |
std::cout << " --- WARNING : mt = " << mt << " is negative... Set to 0."; |
65 |
> |
//std::cout << " --- WARNING : mt = " << mt << " is negative... Set to 0."; |
66 |
|
return 0.; |
67 |
|
} |
68 |
|
return sqrt(mt); |