← Back to branch summary

~christopher-hunt08/maus/maus_integrated_kalman

« back to all changes in this revision

Viewing changes to src/common_cpp/DataStructure/SciFiHelicalPRTrack.cc

Merged Updates

Show diffs side-by-side

added added

removed removed

Lines of Context:
src/common_cpp/DataStructure/SciFiHelicalPRTrack.cc
36
36
    _circle_x0(-1.0),
37
37
    _circle_y0(-1.0),
38
38
    _circle_chisq(-1.0),
39
 
    _chisq(-1.0),
40
 
    _chisq_dof(-1.0),
41
39
    _point_spread(-1.0),
42
40
    _pos0(-1.0, -1.0, -1.0),
43
41
    _phi(0) {
45
43
 
46
44
SciFiHelicalPRTrack::SciFiHelicalPRTrack(int tracker, int num_points, int charge, ThreeVector pos0,
47
45
                                         double phi0, SimpleCircle circle, SimpleLine line_sz,
48
 
                                         double chisq, double chisq_dof, double point_spread,
49
 
                                         DoubleArray phi, SciFiSpacePointPArray spoints,
50
 
                                         const DoubleArray& covariance) :
51
 
    SciFiBasePRTrack(covariance, spoints),
52
 
    _tracker(tracker),
53
 
    _num_points(num_points),
54
 
    _charge(charge),
55
 
    _R(circle.get_R()),
56
 
    _phi0(phi0),
57
 
    _dsdz(line_sz.get_m()),
58
 
    _line_sz_c(line_sz.get_c()),
59
 
    _line_sz_chisq(line_sz.get_chisq()),
60
 
    _circle_x0(circle.get_x0()),
61
 
    _circle_y0(circle.get_y0()),
62
 
    _circle_chisq(circle.get_chisq()),
63
 
    _chisq(chisq),
64
 
    _chisq_dof(chisq_dof),
65
 
    _point_spread(point_spread),
66
 
    _pos0(pos0),
67
 
    _phi(phi) {
68
 
}
69
 
 
70
 
SciFiHelicalPRTrack::SciFiHelicalPRTrack(int tracker, int num_points, int charge, ThreeVector pos0,
71
 
                                         double phi0, SimpleCircle circle, SimpleLine line_sz,
72
46
                                         double point_spread, DoubleArray phi,
73
47
                                         SciFiSpacePointPArray spoints,
74
48
                                         const DoubleArray& covariance) :
84
58
    _circle_x0(circle.get_x0()),
85
59
    _circle_y0(circle.get_y0()),
86
60
    _circle_chisq(circle.get_chisq()),
87
 
    _chisq(0.0),
88
 
    _chisq_dof(0.0),
89
61
    _point_spread(point_spread),
90
62
    _pos0(pos0),
91
63
    _phi(phi) {
92
 
 
93
 
  _chisq = std::sqrt(_circle_chisq*_circle_chisq + _line_sz_chisq*_line_sz_chisq);
94
 
  _chisq_dof = _chisq / (_num_points - 2);
 
64
  this->set_chi_squared(_circle_chisq*_circle_chisq + _line_sz_chisq*_line_sz_chisq);
 
65
  this->set_ndf((2*_num_points) - 5);
95
66
}
96
67
 
97
68
// Destructor
114
85
  _circle_chisq = htrk.get_circle_chisq();
115
86
  _circle_x0 = htrk.get_circle_x0();
116
87
  _circle_y0 = htrk.get_circle_y0();
117
 
  _chisq = htrk.get_chisq();
118
 
  _chisq_dof = htrk.get_chisq_dof();
119
88
  _num_points = htrk.get_num_points();
120
89
  _charge = htrk.get_charge();
121
90
  _tracker = htrk.get_tracker();
138
107
    _circle_x0(htrk.get_circle_x0()),
139
108
    _circle_y0(htrk.get_circle_y0()),
140
109
    _circle_chisq(htrk.get_circle_chisq()),
141
 
    _chisq(htrk.get_chisq()),
142
 
    _chisq_dof(htrk.get_chisq_dof()),
143
110
    _point_spread(htrk.get_point_spread()),
144
111
    _pos0(htrk.get_pos0()),
145
112
    _phi(htrk.get_phi()) {