1 |
khahn |
1.1 |
#ifndef ANGLES
|
2 |
|
|
#define ANGLES
|
3 |
|
|
|
4 |
|
|
#include "TLorentzVector.h"
|
5 |
khahn |
1.2 |
#include "EventData.h"
|
6 |
|
|
|
7 |
|
|
#define VARLIST_ANGLES "costheta1/F:costheta2/F:Phi/F:costhetastar/F:Phi1/F:Phi2/F:phi1/F:phi2/F:theta1/F:theta2/F:bdt/F"
|
8 |
|
|
|
9 |
khahn |
1.1 |
|
10 |
|
|
class Angles {
|
11 |
|
|
public:
|
12 |
khahn |
1.2 |
float costheta1, costheta2, Phi; // decay
|
13 |
|
|
float costhetastar, Phi1; // production
|
14 |
|
|
float Phi2, phi1, phi2; // extra stuff done by JHU
|
15 |
|
|
float theta1, theta2; // for ANL_JHU
|
16 |
|
|
float bdt;
|
17 |
|
|
// float mz1, mz2, m4l;
|
18 |
khahn |
1.1 |
};
|
19 |
|
|
|
20 |
|
|
|
21 |
|
|
|
22 |
|
|
class LabVectors {
|
23 |
|
|
public:
|
24 |
|
|
TLorentzVector vec4l;
|
25 |
|
|
TLorentzVector vecz1;
|
26 |
|
|
TLorentzVector vecz2;
|
27 |
|
|
TLorentzVector vecl1p;
|
28 |
|
|
TLorentzVector vecl1m;
|
29 |
|
|
TLorentzVector vecl2p;
|
30 |
|
|
TLorentzVector vecl2m;
|
31 |
khahn |
1.2 |
|
32 |
|
|
TLorentzVector vec4l_c;
|
33 |
|
|
TLorentzVector vecz1_c;
|
34 |
|
|
TLorentzVector vecz2_c;
|
35 |
|
|
TLorentzVector vecl1p_c;
|
36 |
|
|
TLorentzVector vecl1m_c;
|
37 |
|
|
TLorentzVector vecl2p_c;
|
38 |
|
|
TLorentzVector vecl2m_c;
|
39 |
|
|
|
40 |
|
|
int ncorrections1;
|
41 |
|
|
int ncorrections2;
|
42 |
khahn |
1.1 |
};
|
43 |
|
|
|
44 |
khahn |
1.2 |
|
45 |
khahn |
1.1 |
void getAngles( LabVectors &l, Angles &a );
|
46 |
|
|
void getAngles_JHU( LabVectors &l, Angles &a );
|
47 |
|
|
void getAngles_ANL_NWU( LabVectors &l, Angles &a );
|
48 |
|
|
|
49 |
khahn |
1.2 |
void fillAngles( EventData &, Angles &a );
|
50 |
|
|
|
51 |
khahn |
1.1 |
#endif
|