|
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
|
#include "DataFormats/Candidate/interface/VertexCompositePtrCandidate.h"
using namespace reco;
VertexCompositePtrCandidate::VertexCompositePtrCandidate(Charge q,
const LorentzVector& p4,
const Point& vtx,
const CovarianceMatrix& err,
double chi2,
double ndof,
int pdgId,
int status,
bool integerCharge)
: CompositePtrCandidate(q, p4, vtx, pdgId, status, integerCharge), chi2_(chi2), ndof_(ndof), time_(0.) {
setCovariance(err);
}
VertexCompositePtrCandidate::VertexCompositePtrCandidate(Charge q,
const LorentzVector& p4,
const Point& vtx,
double time,
const CovarianceMatrix4D& err,
double chi2,
double ndof,
int pdgId,
int status,
bool integerCharge)
: CompositePtrCandidate(q, p4, vtx, pdgId, status, integerCharge), chi2_(chi2), ndof_(ndof), time_(time) {
setCovariance(err);
}
VertexCompositePtrCandidate::~VertexCompositePtrCandidate() {}
VertexCompositePtrCandidate* VertexCompositePtrCandidate::clone() const {
return new VertexCompositePtrCandidate(*this);
}
void VertexCompositePtrCandidate::fillVertexCovariance(CovarianceMatrix& err) const {
Error4D temp;
fillVertexCovariance(temp);
err = temp.Sub<Error>(0, 0);
}
void VertexCompositePtrCandidate::fillVertexCovariance(CovarianceMatrix4D& err) const {
index idx = 0;
for (index i = 0; i < dimension4D; ++i)
for (index j = 0; j <= i; ++j)
err(i, j) = covariance_[idx++];
}
void VertexCompositePtrCandidate::setCovariance(const CovarianceMatrix& err) {
index idx = 0;
for (index i = 0; i < dimension4D; ++i) {
for (index j = 0; j <= i; ++j) {
if (i == dimension || j == dimension) {
covariance_[idx++] = 0.0;
} else {
covariance_[idx++] = err(i, j);
}
}
}
}
void VertexCompositePtrCandidate::setCovariance(const CovarianceMatrix4D& err) {
index idx = 0;
for (index i = 0; i < dimension4D; ++i)
for (index j = 0; j <= i; ++j)
covariance_[idx++] = err(i, j);
}
|