FairRoot/PandaRoot
MeasuredStateOnPlane.h
Go to the documentation of this file.
1 /* Copyright 2008-2010, Technische Universitaet Muenchen,
2  Authors: Christian Hoeppner & Sebastian Neubert & Johannes Rauch
3 
4  This file is part of GENFIT.
5 
6  GENFIT is free software: you can redistribute it and/or modify
7  it under the terms of the GNU Lesser General Public License as published
8  by the Free Software Foundation, either version 3 of the License, or
9  (at your option) any later version.
10 
11  GENFIT is distributed in the hope that it will be useful,
12  but WITHOUT ANY WARRANTY; without even the implied warranty of
13  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
14  GNU Lesser General Public License for more details.
15 
16  You should have received a copy of the GNU Lesser General Public License
17  along with GENFIT. If not, see <http://www.gnu.org/licenses/>.
18 */
19 
24 #ifndef genfit_MeasuredStateOnPlane_h
25 #define genfit_MeasuredStateOnPlane_h
26 
27 #include "StateOnPlane.h"
28 #include "AbsTrackRep.h"
29 
30 #include <TMatrixDSym.h>
31 
32 
33 namespace genfit {
34 
39 
40  public:
41 
42  MeasuredStateOnPlane(const AbsTrackRep* rep = NULL);
43  MeasuredStateOnPlane(const TVectorD& state, const TMatrixDSym& cov, const genfit::SharedPlanePtr& plane, const AbsTrackRep* rep);
44  MeasuredStateOnPlane(const TVectorD& state, const TMatrixDSym& cov, const genfit::SharedPlanePtr& plane, const AbsTrackRep* rep, const TVectorD& auxInfo);
46  MeasuredStateOnPlane(const StateOnPlane& state, const TMatrixDSym& cov);
47 
49  void swap(MeasuredStateOnPlane& other); // nothrow
50 
51  virtual ~MeasuredStateOnPlane() {}
52 
53  const TMatrixDSym& getCov() const {return cov_;}
54  TMatrixDSym& getCov() {return cov_;}
55 
57  void blowUpCov(double blowUpFac, bool resetOffDiagonals = true, double maxVal = -1.);
58 
59  void setStateCov(const TVectorD& state, const TMatrixDSym& cov) {setState(state); setCov(cov);}
60  void setStateCovPlane(const TVectorD& state, const TMatrixDSym& cov, const SharedPlanePtr& plane) {setStatePlane(state, plane); setCov(cov);}
61  void setCov(const TMatrixDSym& cov) {if(cov_.GetNrows() == 0) cov_.ResizeTo(cov); cov_ = cov;}
62 
63  // Shortcuts to TrackRep functions
64  TMatrixDSym get6DCov() const {return getRep()->get6DCov(*this);};
65  void getPosMomCov(TVector3& pos, TVector3& mom, TMatrixDSym& cov) const {getRep()->getPosMomCov(*this, pos, mom, cov);}
66  void get6DStateCov(TVectorD& stateVec, TMatrixDSym& cov) const {getRep()->get6DStateCov(*this, stateVec, cov);}
67  double getMomVar() const {return getRep()->getMomVar(*this);}
68 
69  void setPosMomErr(const TVector3& pos, const TVector3& mom, const TVector3& posErr, const TVector3& momErr) {getRep()->setPosMomErr(*this, pos, mom, posErr, momErr);}
70  void setPosMomCov(const TVector3& pos, const TVector3& mom, const TMatrixDSym& cov6x6) {getRep()->setPosMomCov(*this, pos, mom, cov6x6);}
71  void setPosMomCov(const TVectorD& state6, const TMatrixDSym& cov6x6) {getRep()->setPosMomCov(*this, state6, cov6x6);}
72 
73 
74  virtual void Print(Option_t* option = "") const;
75 
76  protected:
77 
78  TMatrixDSym cov_;
79 
80  public:
81  ClassDef(MeasuredStateOnPlane,1)
82 
83 };
84 
85 
89 MeasuredStateOnPlane calcAverageState(const MeasuredStateOnPlane& forwardState, const MeasuredStateOnPlane& backwardState);
90 
91 
93  StateOnPlane::swap(other);
94  this->cov_.ResizeTo(other.cov_);
95  std::swap(this->cov_, other.cov_);
96 }
97 
99  StateOnPlane(rep), cov_(0,0)
100 {
101  if (rep != NULL) {
102  cov_.ResizeTo(rep->getDim(), rep->getDim());
103  }
104 }
105 
106 inline MeasuredStateOnPlane::MeasuredStateOnPlane(const TVectorD& state, const TMatrixDSym& cov, const SharedPlanePtr& plane, const AbsTrackRep* rep) :
107  StateOnPlane(state, plane, rep), cov_(cov)
108 {
109  assert(rep != NULL);
110  //assert(cov_.GetNcols() == (signed)rep->getDim());
111 }
112 
113 inline MeasuredStateOnPlane::MeasuredStateOnPlane(const TVectorD& state, const TMatrixDSym& cov, const SharedPlanePtr& plane, const AbsTrackRep* rep, const TVectorD& auxInfo) :
114  StateOnPlane(state, plane, rep, auxInfo), cov_(cov)
115 {
116  assert(rep != NULL);
117  //assert(cov_.GetNcols() == (signed)rep->getDim());
118 }
119 
121  StateOnPlane(o), cov_(o.cov_)
122 {
123 }
124 
125 inline MeasuredStateOnPlane::MeasuredStateOnPlane(const StateOnPlane& state, const TMatrixDSym& cov) :
126  StateOnPlane(state), cov_(cov)
127 {
128  //assert(cov_.GetNcols() == (signed)getRep()->getDim());
129 }
130 
132  swap(other);
133  return *this;
134 }
135 
136 
137 } /* End of namespace genfit */
140 #endif // genfit_MeasuredStateOnPlane_h
void setPosMomCov(const TVectorD &state6, const TMatrixDSym &cov6x6)
TVector3 pos
void swap(MeasuredStateOnPlane &other)
void get6DStateCov(TVectorD &stateVec, TMatrixDSym &cov) const
boost::shared_ptr< genfit::DetPlane > SharedPlanePtr
Shared Pointer to a DetPlane.
void setPosMomErr(const TVector3 &pos, const TVector3 &mom, const TVector3 &posErr, const TVector3 &momErr)
void setStateCovPlane(const TVectorD &state, const TMatrixDSym &cov, const SharedPlanePtr &plane)
virtual unsigned int getDim() const =0
Get the dimension of the state vector used by the track representation.
Abstract base class for a track representation.
Definition: AbsTrackRep.h:66
const TMatrixDSym & getCov() const
virtual void get6DStateCov(const MeasuredStateOnPlane &state, TVectorD &stateVec, TMatrixDSym &cov) const
Translates MeasuredStateOnPlane into 6D state vector (x, y, z, p_x, p_y, p_z) and 6x6 covariance...
StateOnPlane with additional covariance matrix.
void setPosMomCov(const TVector3 &pos, const TVector3 &mom, const TMatrixDSym &cov6x6)
void getPosMomCov(TVector3 &pos, TVector3 &mom, TMatrixDSym &cov) const
void setState(const TVectorD &state)
Definition: StateOnPlane.h:67
void setStatePlane(const TVectorD &state, const SharedPlanePtr &plane)
Definition: StateOnPlane.h:69
Double_t mom
Definition: plot_dirc.C:14
virtual void getPosMomCov(const MeasuredStateOnPlane &state, TVector3 &pos, TVector3 &mom, TMatrixDSym &cov) const =0
Translates MeasuredStateOnPlane into 3D position, momentum and 6x6 covariance.
void setStateCov(const TVectorD &state, const TMatrixDSym &cov)
A state with arbitrary dimension defined in a DetPlane.
Definition: StateOnPlane.h:45
virtual void setPosMomErr(MeasuredStateOnPlane &state, const TVector3 &pos, const TVector3 &mom, const TVector3 &posErr, const TVector3 &momErr) const =0
Set position and momentum and error of state.
MeasuredStateOnPlane calcAverageState(const MeasuredStateOnPlane &forwardState, const MeasuredStateOnPlane &backwardState)
Calculate weighted average between two MeasuredStateOnPlanes.
void blowUpCov(double blowUpFac, bool resetOffDiagonals=true, double maxVal=-1.)
Blow up covariance matrix with blowUpFac. Per default, off diagonals are reset to 0 and the maximum v...
void swap(StateOnPlane &other)
Definition: StateOnPlane.h:174
void setCov(const TMatrixDSym &cov)
virtual void setPosMomCov(MeasuredStateOnPlane &state, const TVector3 &pos, const TVector3 &mom, const TMatrixDSym &cov6x6) const =0
Set position, momentum and covariance of state.
const AbsTrackRep * getRep() const
Definition: StateOnPlane.h:65
MeasuredStateOnPlane & operator=(MeasuredStateOnPlane other)
virtual void Print(Option_t *option="") const
MeasuredStateOnPlane(const AbsTrackRep *rep=NULL)
virtual double getMomVar(const MeasuredStateOnPlane &state) const =0
get the variance of the absolute value of the momentum .
virtual TMatrixDSym get6DCov(const MeasuredStateOnPlane &state) const =0
Get the 6D covariance.
Matrix inversion tools.
Definition: AbsBField.h:29