Public Member Functions | Private Member Functions | Private Attributes | List of all members
bpfit::Lutz Class Reference

Break-point track fitter. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-15/BreakPointFitter/func/Lutz.h"

Public Member Functions

 Lutz (unsigned int n, unsigned int N)
 Construct a fitter for n measurements and N scattering planes. More...
 
void SetMeasurement (unsigned int i, double z, double x, double sigx)
 Add measurements. More...
 
void SetScatteringPlane (unsigned int j, double Z, double sigSJ)
 Set Z location of a scattering plane j. More...
 
void Fit (double *a, double *b, double *alpha, double *chi2)
 Do the fit. More...
 
double Chi2XIi (unsigned int i) const
 After fit, contribution to chi^2 of the ith measurement. More...
 
double Chi2BetaJ (unsigned int j) const
 After fit, contribution to chi^2 of the jth scattering plane. More...
 
double Chi2XI () const
 
double Chi2Beta () const
 
double X (double z) const
 After fit, best-fit x location of track at location z. More...
 

Private Member Functions

void CalcZeta ()
 
void CalcD ()
 
bool CalcAinv ()
 
void Calcq ()
 
void Calcx ()
 
void NullResult ()
 

Private Attributes

unsigned int fn
 Number of measurements. More...
 
unsigned int fN
 Number of scattering planes. More...
 
std::vector< double > fz
 Measurement plane longitudinal locations. More...
 
TVectorT< double > fXI
 Measurements in transverse location. More...
 
std::vector< double > fwx
 Measurement weight (1/sigx^2) More...
 
std::vector< double > fZ
 Location of scattering planes. More...
 
std::vector< double > fwSJ
 Scattering angle weight (1/sigSJ^2) More...
 
std::vector< std::vector< double > > fZeta
 measurement-scattering plane distances More...
 
TMatrixT< double > fD
 Eqn. 31. More...
 
TMatrixTSym< double > fAinv
 System of equations, Eqn. 29. More...
 
double fa
 Fitted intercept. More...
 
double fb
 Fitted slope. More...
 
std::vector< double > fAlpha
 Fitted scatters. More...
 
std::vector< double > fx
 Fitted track locations. More...
 

Detailed Description

Break-point track fitter.

Implementation of ................................................................... OPTIMUM TRACK FITTING IN THE PRESENCE OF MULTIPLE SCATTERING G. LUTZ Nuclear Instruments and Methods in Physics Research A273 (1998) pp. 349-361 ...................................................................

Definition at line 20 of file Lutz.h.

Constructor & Destructor Documentation

bpfit::Lutz::Lutz ( unsigned int  n,
unsigned int  N 
)

Construct a fitter for n measurements and N scattering planes.

Parameters
n- number of measurements to be added
N- number of scattering surfaces to be added

Definition at line 13 of file Lutz.cxx.

References fZeta, MECModelEnuComparisons::i, bpfit::kNn, getGoodRuns4SAM::n, and stan::math::resize().

13  :
14  fn(n),
15  fN(N),
16  fz(n),
17  fXI(n),
18  fwx(n),
19  fZ(N),
20  fwSJ(N),
21  fZeta(n),
22  fD(N+2,n),
23  fAinv(N+2),
24  fa(0),
25  fb(0),
26  fAlpha(fN),
27  fx(n)
28  {
29  if(N >= n)
30  throw BPException(__FILE__, __LINE__, kNn);
31 
32  for(unsigned int i = 0; i < n; i++)
33  fZeta[i].resize(N);
34  }
TMatrixTSym< double > fAinv
System of equations, Eqn. 29.
Definition: Lutz.h:99
std::vector< double > fz
Measurement plane longitudinal locations.
Definition: Lutz.h:92
std::vector< double > fZ
Location of scattering planes.
Definition: Lutz.h:95
unsigned int fn
Number of measurements.
Definition: Lutz.h:90
double fa
Fitted intercept.
Definition: Lutz.h:100
const unsigned int kNn
number of scat surf >= number of measurements
Definition: BPException.h:20
TMatrixT< double > fD
Eqn. 31.
Definition: Lutz.h:98
TVectorT< double > fXI
Measurements in transverse location.
Definition: Lutz.h:93
void resize(T &x, std::vector< int > dims)
Definition: resize.hpp:41
double fb
Fitted slope.
Definition: Lutz.h:101
std::vector< double > fwx
Measurement weight (1/sigx^2)
Definition: Lutz.h:94
std::vector< double > fwSJ
Scattering angle weight (1/sigSJ^2)
Definition: Lutz.h:96
std::vector< double > fAlpha
Fitted scatters.
Definition: Lutz.h:102
std::vector< std::vector< double > > fZeta
measurement-scattering plane distances
Definition: Lutz.h:97
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
std::vector< double > fx
Fitted track locations.
Definition: Lutz.h:103

Member Function Documentation

bool bpfit::Lutz::CalcAinv ( )
private

Inverting the matrix Sometimes this operation fails, producing error: matrix is singular Thus we need to put it into try/catch block to avoid the hard exit.

Definition at line 100 of file Lutz.cxx.

References om::cerr, fillBadChanDBTables::det, allTimeWatchdog::endl, fAinv, fn, fN, fwSJ, fwx, fz, fZeta, MECModelEnuComparisons::i, genie::utils::mec::J(), K, and sum.

Referenced by Fit().

101  {
102  unsigned int i, J, K;
103 
104  double sum00 = 0;
105  double sum01 = 0;
106  double sum11 = 0;
107  double fwxfz;
108  for (i=0; i<fn; ++i) {
109  fwxfz = fwx[i]*fz[i];
110  sum00 += fwx[i];
111  sum01 += fwxfz;
112  sum11 += fwxfz*fz[i];
113  }
114  fAinv[0][0] = sum00;
115  fAinv[0][1] = sum01;
116  fAinv[1][0] = sum01;
117  fAinv[1][1] = sum11;
118 
119  double sum0J2;
120  double sum1J2;
121  double fwxZeta;
122  for (J=0; J<fN; ++J) {
123  sum0J2 = 0;
124  sum1J2 = 0;
125  for (i=0; i<fn; ++i) {
126  fwxZeta = fwx[i]*fZeta[i][J];
127  sum0J2 += fwxZeta;
128  sum1J2 += fz[i]*fwxZeta;
129  }
130  fAinv[0][J+2] = sum0J2;
131  fAinv[J+2][0] = sum0J2;
132  fAinv[1][J+2] = sum1J2;
133  fAinv[J+2][1] = sum1J2;
134  }
135 
136  double sum;
137  for (J=0; J<fN; ++J) {
138  for (K=J; K<fN; ++K) {
139  sum = 0;
140  for (i=0; i<fn; ++i) sum += fwx[i]*fZeta[i][J]*fZeta[i][K];
141  if (J==K) sum += fwSJ[J];
142  fAinv[J+2][K+2] = sum;
143  fAinv[K+2][J+2] = sum;
144  }
145  }
146 
147  double det = 0.0;
148 
149  /// Inverting the matrix
150  /// Sometimes this operation fails, producing error: matrix is singular
151  /// Thus we need to put it into try/catch block to avoid the hard exit.
152  try {
153  fAinv.InvertFast(&det);
154  }
155  catch(...) {
156  std::cerr << __FILE__ << ":" << __LINE__
157  << " Matrix inversion failed" << std::endl;
158  return false;
159  } // end of inverting the matrix
160 
161  if (det<=0.0) {
162  std::cerr << __FILE__ << ":" << __LINE__
163  << " Matrix inversion failed" << std::endl;
164  return false;
165  }
166  return true;
167  }
TMatrixTSym< double > fAinv
System of equations, Eqn. 29.
Definition: Lutz.h:99
std::vector< double > fz
Measurement plane longitudinal locations.
Definition: Lutz.h:92
double J(double q0, double q3, double Enu, double ml)
Definition: MECUtils.cxx:141
OStream cerr
Definition: OStream.cxx:7
unsigned int fn
Number of measurements.
Definition: Lutz.h:90
Double_t K
std::vector< double > fwx
Measurement weight (1/sigx^2)
Definition: Lutz.h:94
std::vector< double > fwSJ
Scattering angle weight (1/sigSJ^2)
Definition: Lutz.h:96
std::vector< std::vector< double > > fZeta
measurement-scattering plane distances
Definition: Lutz.h:97
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
Double_t sum
Definition: plot.C:31
void bpfit::Lutz::CalcD ( )
private

Definition at line 84 of file Lutz.cxx.

References fD, fn, fN, fwx, fz, fZeta, MECModelEnuComparisons::i, and calib::j.

Referenced by Fit().

85  {
86  unsigned int i, j;
87  for (i=0; i<fn; ++i) {
88  fD[0][i] = fwx[i];
89  fD[1][i] = fwx[i]*fz[i];
90  for (j=0; j<fN; ++j) {
91  fD[j+2][i] = fwx[i]*fZeta[i][j];
92  }
93  }
94  }
std::vector< double > fz
Measurement plane longitudinal locations.
Definition: Lutz.h:92
unsigned int fn
Number of measurements.
Definition: Lutz.h:90
TMatrixT< double > fD
Eqn. 31.
Definition: Lutz.h:98
std::vector< double > fwx
Measurement weight (1/sigx^2)
Definition: Lutz.h:94
const double j
Definition: BetheBloch.cxx:29
std::vector< std::vector< double > > fZeta
measurement-scattering plane distances
Definition: Lutz.h:97
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
void bpfit::Lutz::Calcq ( )
private

Definition at line 173 of file Lutz.cxx.

References fa, fAinv, fAlpha, fb, fD, fN, fXI, and calib::j.

Referenced by Fit().

174  {
175  unsigned int j;
176  //
177  // Solve according to equation 35. Don't bother with the second term
178  // which is by construction zero.
179  //
181  q *= fD;
182  q *= fAinv;
183 
184  //
185  // Unpack the resutls
186  //
187  fa = q[0];
188  fb = q[1];
189  for (j=0; j<fN; ++j) fAlpha[j] = q[2+j];
190  }
TMatrixTSym< double > fAinv
System of equations, Eqn. 29.
Definition: Lutz.h:99
double fa
Fitted intercept.
Definition: Lutz.h:100
TMatrixT< double > fD
Eqn. 31.
Definition: Lutz.h:98
TVectorT< double > fXI
Measurements in transverse location.
Definition: Lutz.h:93
double fb
Fitted slope.
Definition: Lutz.h:101
const double j
Definition: BetheBloch.cxx:29
std::vector< double > fAlpha
Fitted scatters.
Definition: Lutz.h:102
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
void bpfit::Lutz::Calcx ( )
private

Definition at line 194 of file Lutz.cxx.

References fa, fAlpha, fb, fn, fN, fx, fz, fZeta, MECModelEnuComparisons::i, and calib::j.

Referenced by Fit().

195  {
196  unsigned int i, j;
197  for (i=0; i<fn; ++i) {
198  fx[i] = fa + fb*fz[i];
199  for (j=0; j<fN; ++j) {
200  if (fZeta[i][j]>0.0) fx[i] += fAlpha[j]*fZeta[i][j];
201  }
202  }
203  }
std::vector< double > fz
Measurement plane longitudinal locations.
Definition: Lutz.h:92
unsigned int fn
Number of measurements.
Definition: Lutz.h:90
double fa
Fitted intercept.
Definition: Lutz.h:100
double fb
Fitted slope.
Definition: Lutz.h:101
const double j
Definition: BetheBloch.cxx:29
std::vector< double > fAlpha
Fitted scatters.
Definition: Lutz.h:102
std::vector< std::vector< double > > fZeta
measurement-scattering plane distances
Definition: Lutz.h:97
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
std::vector< double > fx
Fitted track locations.
Definition: Lutz.h:103
void bpfit::Lutz::CalcZeta ( )
private

Definition at line 69 of file Lutz.cxx.

References fn, fN, fz, fZ, fZeta, MECModelEnuComparisons::i, and calib::j.

Referenced by Fit().

70  {
71  unsigned int i, j;
72  for (i=0; i<fn; ++i) {
73  for (j=0; j<fN; ++j) {
74  fZeta[i][j] = fz[i]-fZ[j];
75  if (fZeta[i][j]<0.0) fZeta[i][j] = 0.0;
76  }
77  }
78  }
std::vector< double > fz
Measurement plane longitudinal locations.
Definition: Lutz.h:92
std::vector< double > fZ
Location of scattering planes.
Definition: Lutz.h:95
unsigned int fn
Number of measurements.
Definition: Lutz.h:90
const double j
Definition: BetheBloch.cxx:29
std::vector< std::vector< double > > fZeta
measurement-scattering plane distances
Definition: Lutz.h:97
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
double bpfit::Lutz::Chi2Beta ( ) const

Definition at line 287 of file Lutz.cxx.

References Chi2BetaJ(), fN, calib::j, and sum.

Referenced by bpfit::BreakPoint::FillTrackNt(), and Fit().

288  {
289  unsigned int j;
290  double sum = 0.0;
291  for (j=0; j<fN; ++j) sum += this->Chi2BetaJ(j);
292  return sum;
293  }
const double j
Definition: BetheBloch.cxx:29
double Chi2BetaJ(unsigned int j) const
After fit, contribution to chi^2 of the jth scattering plane.
Definition: Lutz.cxx:264
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
Double_t sum
Definition: plot.C:31
double bpfit::Lutz::Chi2BetaJ ( unsigned int  j) const

After fit, contribution to chi^2 of the jth scattering plane.

Definition at line 264 of file Lutz.cxx.

References fAlpha, fwSJ, and calib::j.

Referenced by Chi2Beta().

265  {
266  double chi;
267 
268  chi = 0.0-fAlpha[j];
269  chi *= chi;
270  chi *= fwSJ[j];
271 
272  return chi;
273  }
std::vector< double > fwSJ
Scattering angle weight (1/sigSJ^2)
Definition: Lutz.h:96
const double j
Definition: BetheBloch.cxx:29
std::vector< double > fAlpha
Fitted scatters.
Definition: Lutz.h:102
double bpfit::Lutz::Chi2XI ( ) const

Definition at line 277 of file Lutz.cxx.

References Chi2XIi(), fn, MECModelEnuComparisons::i, and sum.

Referenced by bpfit::BreakPoint::FillTrackNt(), and Fit().

278  {
279  unsigned int i;
280  double sum = 0.0;
281  for (i=0; i<fn; ++i) sum += this->Chi2XIi(i);
282  return sum;
283  }
double Chi2XIi(unsigned int i) const
After fit, contribution to chi^2 of the ith measurement.
Definition: Lutz.cxx:251
unsigned int fn
Number of measurements.
Definition: Lutz.h:90
Double_t sum
Definition: plot.C:31
double bpfit::Lutz::Chi2XIi ( unsigned int  i) const

After fit, contribution to chi^2 of the ith measurement.

Definition at line 251 of file Lutz.cxx.

References fwx, fx, fXI, and MECModelEnuComparisons::i.

Referenced by Chi2XI().

252  {
253  double chi;
254 
255  chi = fXI[i]-fx[i];
256  chi *= chi;
257  chi *= fwx[i];
258 
259  return chi;
260  }
TVectorT< double > fXI
Measurements in transverse location.
Definition: Lutz.h:93
std::vector< double > fwx
Measurement weight (1/sigx^2)
Definition: Lutz.h:94
std::vector< double > fx
Fitted track locations.
Definition: Lutz.h:103
void bpfit::Lutz::Fit ( double *  a,
double *  b,
double *  alpha,
double *  chi2 
)

Do the fit.

Parameters
a- on return, fitted intercept. OK to pass 0.
b- on return, fitted initial slope. OK to pass 0.
alpha- on return, fitted scattering slopes. OK to pass 0.
chi2- on return, the chi^2 result of the fit. OK to pass 0.

Definition at line 221 of file Lutz.cxx.

References CalcAinv(), CalcD(), Calcq(), Calcx(), CalcZeta(), Chi2Beta(), Chi2XI(), fa, fAlpha, fb, fN, calib::j, and NullResult().

222  {
223  bool ok;
224  unsigned int j;
225  //
226  // Make the fit
227  //
228  this->CalcZeta();
229  this->CalcD();
230  ok = this->CalcAinv();
231  if (ok) {
232  this->Calcq();
233  this->Calcx();
234  }
235  else this->NullResult();
236  //
237  // Unpack the results
238  //
239  if (a!=0) *a = fa;
240  if (b!=0) *b = fb;
241  if (alpha!=0) {
242  for (j=0; j<fN; ++j) alpha[j] = fAlpha[j];
243  }
244  if (chi2!=0) {
245  *chi2 = this->Chi2XI() + this->Chi2Beta();
246  }
247  }
void Calcx()
Definition: Lutz.cxx:194
double fa
Fitted intercept.
Definition: Lutz.h:100
double chi2()
const double a
double fb
Fitted slope.
Definition: Lutz.h:101
double Chi2XI() const
Definition: Lutz.cxx:277
const double j
Definition: BetheBloch.cxx:29
std::vector< double > fAlpha
Fitted scatters.
Definition: Lutz.h:102
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
void Calcq()
Definition: Lutz.cxx:173
const hit & b
Definition: hits.cxx:21
void CalcZeta()
Definition: Lutz.cxx:69
bool CalcAinv()
Definition: Lutz.cxx:100
double Chi2Beta() const
Definition: Lutz.cxx:287
void CalcD()
Definition: Lutz.cxx:84
void NullResult()
Definition: Lutz.cxx:297
void bpfit::Lutz::NullResult ( )
private

Definition at line 297 of file Lutz.cxx.

References fa, fAlpha, fb, fn, fN, fx, MECModelEnuComparisons::i, and calib::j.

Referenced by Fit().

298  {
299  unsigned int i, j;
300  fa = 0.0;
301  fb = 0.0;
302  for (i=0; i<fn; ++i) fx[i] = 0.0;
303  for (j=0; j<fN; ++j) fAlpha[j] = 0.0;
304  }
unsigned int fn
Number of measurements.
Definition: Lutz.h:90
double fa
Fitted intercept.
Definition: Lutz.h:100
double fb
Fitted slope.
Definition: Lutz.h:101
const double j
Definition: BetheBloch.cxx:29
std::vector< double > fAlpha
Fitted scatters.
Definition: Lutz.h:102
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
std::vector< double > fx
Fitted track locations.
Definition: Lutz.h:103
void bpfit::Lutz::SetMeasurement ( unsigned int  i,
double  z,
double  x,
double  sigx 
)

Add measurements.

Parameters
i- measurement index [0...n)
z- z location of measurement (longitudinal location)
x- x location of measurement (transverse location)
sigx- measurement error on x position

Definition at line 38 of file Lutz.cxx.

References fn, fwx, fXI, fz, MECModelEnuComparisons::i, bpfit::kINDEX, bpfit::kNAN, bpfit::kSIGX, submit_syst::x, and test::z.

Referenced by bpfit::BreakPoint::FitTracks(), and bpfit::DimuonFitter::FitView().

39  {
40  if(i >= fn)
41  throw BPException(__FILE__, __LINE__, kINDEX);
42  if(z != z || x != x || sigx != sigx)
43  throw BPException(__FILE__, __LINE__, kNAN);
44  if(sigx <= 0.0)
45  throw BPException(__FILE__, __LINE__, kSIGX, sigx);
46 
47  fz[i] = z;
48  fXI[i] = x;
49  fwx[i] = 1.0/(sigx*sigx);
50  }
std::vector< double > fz
Measurement plane longitudinal locations.
Definition: Lutz.h:92
const unsigned int kSIGX
non-positive value of sigx
Definition: BPException.h:23
const unsigned int kNAN
nan encountered
Definition: BPException.h:22
unsigned int fn
Number of measurements.
Definition: Lutz.h:90
TVectorT< double > fXI
Measurements in transverse location.
Definition: Lutz.h:93
std::vector< double > fwx
Measurement weight (1/sigx^2)
Definition: Lutz.h:94
z
Definition: test.py:28
const unsigned int kINDEX
index problem
Definition: BPException.h:21
void bpfit::Lutz::SetScatteringPlane ( unsigned int  j,
double  Z,
double  sigSJ 
)

Set Z location of a scattering plane j.

Parameters
j- index of scattering plane (start from 0)
Z- longitudinal location of scattering plane
sigSJ- RMS scattering of track slope at this plane

Definition at line 54 of file Lutz.cxx.

References fN, fwSJ, fZ, calib::j, bpfit::kINDEX, bpfit::kNAN, bpfit::kSIGSJ, and Z.

55  {
56  if(j >= fN)
57  throw BPException(__FILE__, __LINE__, kINDEX);
58  if(Z != Z || sigSJ != sigSJ)
59  throw BPException(__FILE__, __LINE__, kNAN);
60  if(sigSJ <= 0.0)
61  throw BPException(__FILE__, __LINE__, kSIGSJ, sigSJ);
62 
63  fZ[j] = Z;
64  fwSJ[j] = 1.0/(sigSJ*sigSJ);
65  }
const unsigned int kNAN
nan encountered
Definition: BPException.h:22
std::vector< double > fZ
Location of scattering planes.
Definition: Lutz.h:95
Float_t Z
Definition: plot.C:38
std::vector< double > fwSJ
Scattering angle weight (1/sigSJ^2)
Definition: Lutz.h:96
const double j
Definition: BetheBloch.cxx:29
const unsigned int kINDEX
index problem
Definition: BPException.h:21
const unsigned int kSIGSJ
non-positive value of sigSJ
Definition: BPException.h:24
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91
double bpfit::Lutz::X ( double  z) const

After fit, best-fit x location of track at location z.

Definition at line 207 of file Lutz.cxx.

References fa, fAlpha, fb, fN, fZ, calib::j, submit_syst::x, and test::z.

Referenced by bpfit::BreakPoint::LogResiduals().

208  {
209  unsigned int j;
210  double zeta;
211  double x = fa + fb*z;
212  for (j=0; j<fN; ++j) {
213  zeta = z-fZ[j];
214  if (zeta>0.0) x += fAlpha[j]*zeta;
215  }
216  return x;
217  }
std::vector< double > fZ
Location of scattering planes.
Definition: Lutz.h:95
double fa
Fitted intercept.
Definition: Lutz.h:100
double fb
Fitted slope.
Definition: Lutz.h:101
const double j
Definition: BetheBloch.cxx:29
z
Definition: test.py:28
std::vector< double > fAlpha
Fitted scatters.
Definition: Lutz.h:102
unsigned int fN
Number of scattering planes.
Definition: Lutz.h:91

Member Data Documentation

double bpfit::Lutz::fa
private

Fitted intercept.

Definition at line 100 of file Lutz.h.

Referenced by Calcq(), Calcx(), Fit(), NullResult(), and X().

TMatrixTSym<double> bpfit::Lutz::fAinv
private

System of equations, Eqn. 29.

Definition at line 99 of file Lutz.h.

Referenced by CalcAinv(), and Calcq().

std::vector<double> bpfit::Lutz::fAlpha
private

Fitted scatters.

Definition at line 102 of file Lutz.h.

Referenced by Calcq(), Calcx(), Chi2BetaJ(), Fit(), NullResult(), and X().

double bpfit::Lutz::fb
private

Fitted slope.

Definition at line 101 of file Lutz.h.

Referenced by Calcq(), Calcx(), Fit(), NullResult(), and X().

TMatrixT<double> bpfit::Lutz::fD
private

Eqn. 31.

Definition at line 98 of file Lutz.h.

Referenced by CalcD(), and Calcq().

unsigned int bpfit::Lutz::fn
private

Number of measurements.

Definition at line 90 of file Lutz.h.

Referenced by CalcAinv(), CalcD(), Calcx(), CalcZeta(), Chi2XI(), NullResult(), and SetMeasurement().

unsigned int bpfit::Lutz::fN
private

Number of scattering planes.

Definition at line 91 of file Lutz.h.

Referenced by CalcAinv(), CalcD(), Calcq(), Calcx(), CalcZeta(), Chi2Beta(), Fit(), NullResult(), SetScatteringPlane(), and X().

std::vector<double> bpfit::Lutz::fwSJ
private

Scattering angle weight (1/sigSJ^2)

Definition at line 96 of file Lutz.h.

Referenced by CalcAinv(), Chi2BetaJ(), and SetScatteringPlane().

std::vector<double> bpfit::Lutz::fwx
private

Measurement weight (1/sigx^2)

Definition at line 94 of file Lutz.h.

Referenced by CalcAinv(), CalcD(), Chi2XIi(), and SetMeasurement().

std::vector<double> bpfit::Lutz::fx
private

Fitted track locations.

Definition at line 103 of file Lutz.h.

Referenced by Calcx(), Chi2XIi(), and NullResult().

TVectorT<double> bpfit::Lutz::fXI
private

Measurements in transverse location.

Definition at line 93 of file Lutz.h.

Referenced by Calcq(), Chi2XIi(), and SetMeasurement().

std::vector<double> bpfit::Lutz::fz
private

Measurement plane longitudinal locations.

Definition at line 92 of file Lutz.h.

Referenced by CalcAinv(), CalcD(), Calcx(), CalcZeta(), and SetMeasurement().

std::vector<double> bpfit::Lutz::fZ
private

Location of scattering planes.

Definition at line 95 of file Lutz.h.

Referenced by CalcZeta(), SetScatteringPlane(), and X().

std::vector< std::vector<double> > bpfit::Lutz::fZeta
private

measurement-scattering plane distances

Definition at line 97 of file Lutz.h.

Referenced by CalcAinv(), CalcD(), Calcx(), CalcZeta(), and Lutz().


The documentation for this class was generated from the following files: