sistrip::LinearFit Class Reference

#include <Utility.h>

Classes
class	Params

Public Member Functions
void	add (const float &value_x, const float &value_y)
void	add (const float &value_x, const float &value_y, const float &error_y)
void	fit (Params &fit_params)
	LinearFit ()
	~LinearFit ()

Private Attributes
std::vector< float >	e_
float	ss_
float	sx_
float	sy_
std::vector< float >	x_
std::vector< float >	y_

Detailed Description

Definition at line 11 of file Utility.h.

Constructor & Destructor Documentation

LinearFit()

LinearFit::LinearFit

(

)

Definition at line 6 of file Utility.cc.

: x_(), y_(), e_(), ss_(0.), sx_(0.), sy_(0.) { ; }

~LinearFit()

sistrip::LinearFit::~LinearFit ( )

inline

Definition at line 15 of file Utility.h.

{ ; }

Member Function Documentation

add() [1/2]

void LinearFit::add	(	const float &	value_x,
		const float &	value_y
	)

Definition at line 10 of file Utility.cc.

                                                         {
  float e = 1.;  // default
  x_.push_back(x);
  y_.push_back(y);
  e_.push_back(e);
  float wt = 1. / sqrt(e);
  ss_ += wt;
  sx_ += x * wt;
  sy_ += y * wt;
}

References MillePedeFileConverter_cfg::e, and mathSSE::sqrt().

Referenced by OptoScanAlgorithm::analyse(), counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

add() [2/2]

void LinearFit::add	(	const float &	value_x,
		const float &	value_y,
		const float &	error_y
	)

Definition at line 23 of file Utility.cc.

                                                                         {
  if (e > 0.) {
    x_.push_back(x);
    y_.push_back(y);
    e_.push_back(e);
    float wt = 1. / sqrt(e);
    ss_ += wt;
    sx_ += x * wt;
    sy_ += y * wt;
  }
}

References MillePedeFileConverter_cfg::e, and mathSSE::sqrt().

Referenced by counter.Counter::register(), SequenceTypes.Task::remove(), and SequenceTypes.Task::replace().

fit()

void LinearFit::fit

(

Params &

fit_params

)

Definition at line 37 of file Utility.cc.

                                         {
  float s2 = 0.;
  float b = 0;
  for (uint16_t i = 0; i < x_.size(); i++) {
    float t = (x_[i] - sx_ / ss_) / e_[i];
    s2 += t * t;
    b += t * y_[i] / e_[i];
  }
 
  // Set parameters
  params.n_ = x_.size();
  params.b_ = b / s2;
  params.a_ = (sy_ - sx_ * params.b_) / ss_;
  params.erra_ = sqrt((1. + (sx_ * sx_) / (ss_ * s2)) / ss_);
  params.errb_ = sqrt(1. / s2);
 
  /*
    params.chi2_ = 0.;
    *q=1.0;
    if (mwt == 0) {
    for (i=1;i<=ndata;i++)
    *chi2 += SQR(y[i]-(*a)-(*b)*x[i]);
    sigdat=sqrt((*chi2)/(ndata-2));
    *sigb *= sigdat;
    */
}

References b, mps_fire::i, CalibrationSummaryClient_cfi::params, indexGen::s2, mathSSE::sqrt(), and OrderedSet::t.

Referenced by OptoScanAlgorithm::analyse(), and trackingPlots.Iteration::modules().

Member Data Documentation

e_

std::vector<float> sistrip::LinearFit::e_

private

Definition at line 44 of file Utility.h.

ss_

float sistrip::LinearFit::ss_

private

Definition at line 45 of file Utility.h.

sx_

float sistrip::LinearFit::sx_

private

Definition at line 46 of file Utility.h.

sy_

float sistrip::LinearFit::sy_

private

Definition at line 47 of file Utility.h.

x_

std::vector<float> sistrip::LinearFit::x_

private

Definition at line 42 of file Utility.h.

y_

std::vector<float> sistrip::LinearFit::y_

private

Definition at line 43 of file Utility.h.