|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
|
#ifndef DataFormats_TauReco_PFTauTransverseImpactParameter_h
#define DataFormats_TauReco_PFTauTransverseImpactParameter_h
/* class PFTauTransverseImpactParameter
*
* Stores information on the impact paramters and flight length of the hadronic decay of a tau lepton
*
* author: Ian M. Nugent
* This work is based on the impact parameter work by Rosamaria Venditti and reconstructing the 3 prong taus.
* The idea of the fully reconstructing the tau using a kinematic fit comes from
* Lars Perchalla and Philip Sauerland Theses under Achim Stahl supervision. This
* work was continued by Ian M. Nugent and Vladimir Cherepanov.
* Thanks goes to Christian Veelken and Evan Klose Friis for their help and suggestions.
*/
#include "DataFormats/Math/interface/Vector3D.h"
#include "DataFormats/Math/interface/Point3D.h"
#include "DataFormats/Math/interface/Error.h"
#include "DataFormats/VertexReco/interface/Vertex.h"
#include "DataFormats/VertexReco/interface/VertexFwd.h"
#include "TVector3.h"
namespace reco {
class PFTauTransverseImpactParameter {
enum { dimension = 3 };
enum { covarianceSize = dimension * (dimension + 1) / 2 };
public:
typedef math::Error<dimension>::type CovMatrix;
typedef math::XYZPoint Point;
typedef math::XYZVector Vector;
PFTauTransverseImpactParameter() {}
/// constructor from values
PFTauTransverseImpactParameter(const Point&, double, double, const Point&, double, double, const VertexRef&);
PFTauTransverseImpactParameter(const Point&,
double,
double,
const Point&,
double,
double,
const VertexRef&,
const Point&,
double,
const VertexRef&);
virtual ~PFTauTransverseImpactParameter() {}
PFTauTransverseImpactParameter* clone() const;
const Point& dxy_PCA() const { return pca_; }
double dxy() const { return dxy_; }
double dxy_error() const { return dxy_error_; }
double dxy_Sig() const { return (dxy_error_ != 0) ? (dxy_ / dxy_error_) : 0.; }
const Point& ip3d_PCA() const { return pca3d_; }
double ip3d() const { return ip3d_; }
double ip3d_error() const { return ip3d_error_; }
double ip3d_Sig() const { return (ip3d_error_ != 0) ? (ip3d_ / ip3d_error_) : 0.; }
const VertexRef& primaryVertex() const { return PV_; }
Point primaryVertexPos() const;
CovMatrix primaryVertexCov() const;
bool hasSecondaryVertex() const { return hasSV_; }
const Vector& flightLength() const;
double flightLengthSig() const;
CovMatrix flightLengthCov() const;
const VertexRef& secondaryVertex() const { return SV_; }
Point secondaryVertexPos() const;
CovMatrix secondaryVertexCov() const;
private:
Point pca_;
double dxy_{0};
double dxy_error_{0};
Point pca3d_;
double ip3d_{0};
double ip3d_error_{0};
VertexRef PV_;
bool hasSV_{false};
Vector FlightLength_;
double FlightLengthSig_{0};
VertexRef SV_;
};
} // namespace reco
#endif
|