MeasuredStateOnPlane.h

Go to the documentation of this file.

/* Copyright 2008-2010, Technische Universitaet Muenchen,
   Authors: Christian Hoeppner & Sebastian Neubert & Johannes Rauch

   This file is part of GENFIT.

   GENFIT is free software: you can redistribute it and/or modify
   it under the terms of the GNU Lesser General Public License as published
   by the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.

   GENFIT is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public License
   along with GENFIT.  If not, see <http://www.gnu.org/licenses/>.
*/

#ifndef genfit_MeasuredStateOnPlane_h
#define genfit_MeasuredStateOnPlane_h

#include "StateOnPlane.h"
#include "AbsTrackRep.h"

#include <TMatrixDSym.h>

namespace genfit {

class MeasuredStateOnPlane : public StateOnPlane {

 public:
  MeasuredStateOnPlane(const AbsTrackRep *rep = nullptr);
  MeasuredStateOnPlane(const TVectorD &state, const TMatrixDSym &cov, const genfit::SharedPlanePtr &plane, const AbsTrackRep *rep);
  MeasuredStateOnPlane(const TVectorD &state, const TMatrixDSym &cov, const genfit::SharedPlanePtr &plane, const AbsTrackRep *rep, const TVectorD &auxInfo);
  MeasuredStateOnPlane(const MeasuredStateOnPlane &o);
  MeasuredStateOnPlane(const StateOnPlane &state, const TMatrixDSym &cov);

  MeasuredStateOnPlane &operator=(MeasuredStateOnPlane other);
  void swap(MeasuredStateOnPlane &other); // nothrow

  virtual ~MeasuredStateOnPlane() {}

  const TMatrixDSym &getCov() const { return cov_; }
  TMatrixDSym &getCov() { return cov_; }

  void blowUpCov(double blowUpFac, bool resetOffDiagonals = true, double maxVal = -1.);

  void setStateCov(const TVectorD &state, const TMatrixDSym &cov)
  {
    setState(state);
    setCov(cov);
  }
  void setStateCovPlane(const TVectorD &state, const TMatrixDSym &cov, const SharedPlanePtr &plane)
  {
    setStatePlane(state, plane);
    setCov(cov);
  }
  void setCov(const TMatrixDSym &cov)
  {
    if (cov_.GetNrows() == 0)
      cov_.ResizeTo(cov);
    cov_ = cov;
  }

  // Shortcuts to TrackRep functions
  TMatrixDSym get6DCov() const { return getRep()->get6DCov(*this); };
  void getPosMomCov(TVector3 &pos, TVector3 &mom, TMatrixDSym &cov) const { getRep()->getPosMomCov(*this, pos, mom, cov); }
  void get6DStateCov(TVectorD &stateVec, TMatrixDSym &cov) const { getRep()->get6DStateCov(*this, stateVec, cov); }
  double getMomVar() const { return getRep()->getMomVar(*this); }

  void setPosMomErr(const TVector3 &pos, const TVector3 &mom, const TVector3 &posErr, const TVector3 &momErr) { getRep()->setPosMomErr(*this, pos, mom, posErr, momErr); }
  void setPosMomCov(const TVector3 &pos, const TVector3 &mom, const TMatrixDSym &cov6x6) { getRep()->setPosMomCov(*this, pos, mom, cov6x6); }
  void setPosMomCov(const TVectorD &state6, const TMatrixDSym &cov6x6) { getRep()->setPosMomCov(*this, state6, cov6x6); }

  virtual void Print(Option_t *option = "") const;

 protected:
  TMatrixDSym cov_;

 public:
  ClassDef(MeasuredStateOnPlane, 1)
};

MeasuredStateOnPlane calcAverageState(const MeasuredStateOnPlane &forwardState, const MeasuredStateOnPlane &backwardState);

inline void MeasuredStateOnPlane::swap(MeasuredStateOnPlane &other)
{
  StateOnPlane::swap(other);
  this->cov_.ResizeTo(other.cov_);
  std::swap(this->cov_, other.cov_);
}

inline MeasuredStateOnPlane::MeasuredStateOnPlane(const AbsTrackRep *rep) : StateOnPlane(rep), cov_(0, 0)
{
  if (rep != nullptr) {
    cov_.ResizeTo(rep->getDim(), rep->getDim());
  }
}

inline MeasuredStateOnPlane::MeasuredStateOnPlane(const TVectorD &state, const TMatrixDSym &cov, const SharedPlanePtr &plane, const AbsTrackRep *rep)
  : StateOnPlane(state, plane, rep), cov_(cov)
{
  assert(rep != nullptr);
  // assert(cov_.GetNcols() == (signed)rep->getDim());
}

inline MeasuredStateOnPlane::MeasuredStateOnPlane(const TVectorD &state, const TMatrixDSym &cov, const SharedPlanePtr &plane, const AbsTrackRep *rep, const TVectorD &auxInfo)
  : StateOnPlane(state, plane, rep, auxInfo), cov_(cov)
{
  assert(rep != nullptr);
  // assert(cov_.GetNcols() == (signed)rep->getDim());
}

inline MeasuredStateOnPlane::MeasuredStateOnPlane(const MeasuredStateOnPlane &o) : StateOnPlane(o), cov_(o.cov_) {}

inline MeasuredStateOnPlane::MeasuredStateOnPlane(const StateOnPlane &state, const TMatrixDSym &cov) : StateOnPlane(state), cov_(cov)
{
  // assert(cov_.GetNcols() == (signed)getRep()->getDim());
}

inline MeasuredStateOnPlane &MeasuredStateOnPlane::operator=(MeasuredStateOnPlane other)
{
  swap(other);
  return *this;
}

} /* End of namespace genfit */
#endif // genfit_MeasuredStateOnPlane_h