TkFittedLasBeam.cc

Go to the documentation of this file.

#include "DataFormats/Alignment/interface/TkFittedLasBeam.h"
#include "FWCore/Utilities/interface/Exception.h"
#include <vector>

TkFittedLasBeam::TkFittedLasBeam()
{
}

TkFittedLasBeam::TkFittedLasBeam(const TkLasBeam &lasBeam) : TkLasBeam(lasBeam)
{
}

void TkFittedLasBeam::setParameters(unsigned int parametrisation,
                    const std::vector<Scalar> &params,
                    const AlgebraicSymMatrix &paramCovariance,
                    const AlgebraicMatrix &derivatives,
                    unsigned int firstFixedParam, float chi2)
{
  parametrisation_ = parametrisation;
  parameters_ = params;
  paramCovariance_ = paramCovariance;
  derivatives_ = derivatives;
  firstFixedParameter_ = firstFixedParam;
  chi2_ = chi2;

  // check integrity
  if (parameters_.size() != static_cast<unsigned int>(derivatives_.num_col()) // # parameters
      || static_cast<unsigned int>(derivatives_.num_row()) != this->getData().size() // # hits
      || firstFixedParameter_ > parameters_.size() // index 'fixed' might be the
      || static_cast<unsigned int>(paramCovariance_.num_row()) != firstFixedParameter_) {
    throw cms::Exception("BadInput")
      << "[TkFittedLasBeam::setParameters] with inconsistent sizes: (parametrisation " 
      << parametrisation << "):\n"
      << parameters_.size() << " parameters,\n"
      << derivatives_.num_row() << "x" << derivatives_.num_col() << " derivatives,\n"
      << "firstFixed = " << firstFixedParameter_ << " (i.e. "
      << static_cast<int>(parameters_.size()) - static_cast<int>(firstFixedParameter_)
      << " global parameters),\n" // cast on line before to allow difference to be < 0, [un]signed!
      << "cov. matrix size " << paramCovariance_.num_row() <<  ".\n";
  }
}