| 1 |
– |
#include "Angles.h" |
| 1 |
|
#include <iostream> |
| 2 |
+ |
#include "Angles.h" |
| 3 |
+ |
#include "EventData.h" |
| 4 |
+ |
|
| 5 |
+ |
void fillAngles( EventData &data, Angles &a ) { |
| 6 |
+ |
|
| 7 |
+ |
LabVectors l; |
| 8 |
+ |
|
| 9 |
+ |
if( data.Z1leptons[0].charge > 0 ) { |
| 10 |
+ |
l.vecl1p = data.Z1leptons[0].vec; |
| 11 |
+ |
l.vecl1m = data.Z1leptons[1].vec; |
| 12 |
+ |
} else { |
| 13 |
+ |
l.vecl1p = data.Z1leptons[1].vec; |
| 14 |
+ |
l.vecl1m = data.Z1leptons[0].vec; |
| 15 |
+ |
} |
| 16 |
+ |
|
| 17 |
+ |
if( data.Z2leptons[0].charge > 0 ) { |
| 18 |
+ |
l.vecl2p = data.Z2leptons[0].vec; |
| 19 |
+ |
l.vecl2m = data.Z2leptons[1].vec; |
| 20 |
+ |
} else { |
| 21 |
+ |
l.vecl2p = data.Z2leptons[1].vec; |
| 22 |
+ |
l.vecl2m = data.Z2leptons[0].vec; |
| 23 |
+ |
} |
| 24 |
+ |
|
| 25 |
+ |
l.vecz1 = (l.vecl1p+l.vecl1m); |
| 26 |
+ |
l.vecz2 = (l.vecl2p+l.vecl2m); |
| 27 |
+ |
l.vec4l = (l.vecz1+l.vecz2); |
| 28 |
+ |
|
| 29 |
+ |
getAngles( l, a ); |
| 30 |
+ |
} |
| 31 |
|
|
| 32 |
|
void getAngles( LabVectors &l, |
| 33 |
|
Angles &a ) { |
| 177 |
|
TVector3 vecpp_in_Z1coords_Xframe = vecpp.Transform(rotation); |
| 178 |
|
a.Phi1 = vecpp_in_Z1coords_Xframe.Phi(); |
| 179 |
|
|
| 180 |
< |
a.mz1 = (l.vecl1m+l.vecl1p).M(); |
| 181 |
< |
a.mz2 = (l.vecl2m+l.vecl2p).M(); |
| 182 |
< |
a.m4l = (l.vec4l).M(); |
| 180 |
> |
// a.mz1 = (l.vecl1m+l.vecl1p).M(); |
| 181 |
> |
// a.mz2 = (l.vecl2m+l.vecl2p).M(); |
| 182 |
> |
// a.m4l = (l.vec4l).M(); |
| 183 |
|
|
| 184 |
|
}; |
| 185 |
|
|
