16 |
|
{ |
17 |
|
double deta = eta1 - eta2; |
18 |
|
double dphi = deltaPhi(phi1, phi2); |
19 |
< |
return std::sqrt(deta*deta + dphi*dphi); |
19 |
> |
return TMath::Sqrt(deta*deta + dphi*dphi); |
20 |
|
} |
21 |
|
|
22 |
|
|
43 |
|
H1.SetMag( ( - n2.Py() * V.Px() + n2.Px() * V.Py() ) / (sin(n1.Theta()) *det ) ); |
44 |
|
H2.SetMag( ( + n1.Py() * V.Px() - n1.Px() * V.Py() ) / (sin(n2.Theta()) *det ) ); |
45 |
|
|
46 |
< |
float mass=std::sqrt( TMath::Power( (H1.Mag()+H2.Mag()),2 ) - TMath::Power(( ( H1+H2 ).Mag()),2) ); |
46 |
> |
float mass=TMath::Sqrt( TMath::Power( (H1.Mag()+H2.Mag()),2 ) - TMath::Power(( ( H1+H2 ).Mag()),2) ); |
47 |
|
|
48 |
|
return mass; |
49 |
|
|