genf::PointHit Class Reference

#include "PointHit.h"

Inheritance diagram for genf::PointHit:

Public Member Functions
	PointHit ()
	PointHit (TVector3 point, double res)
	PointHit (TVector3 point, std::vector< double > &ref)
	PointHit (TVector3 point, TVector3 res)
virtual	~PointHit ()
virtual GFAbsRecoHit *	clone ()
	Get clone of this object. More...
virtual TMatrixT< Double_t >	getHMatrix (const GFAbsTrackRep *stateVector)
	Get transformation matrix. Transformation between hit coordinates and track representation coordinates. More...
virtual TMatrixT< Double_t >	getHMatrix (const GFAbsTrackRep *stateVector, const Double_t &, const Double_t &)
virtual const GFDetPlane &	getDetPlane (GFAbsTrackRep *rep)
	Returns the detector plane object for this hit and a given track representation. More...
virtual TMatrixT< Double_t >	getHitCoord (const GFDetPlane &plane, const GFDetPlane &planePrev)
	Get hit coordinates in a specific detector plane. More...
virtual TMatrixT< Double_t >	getHitCoord (const GFDetPlane &plane)
virtual TMatrixT< Double_t >	getHitCov (const GFDetPlane &plane)
	Get hit covariances in a specific detector plane. More...
virtual TMatrixT< Double_t >	getHitCov (const GFDetPlane &plane, const GFDetPlane &planePrev, const TMatrixT< Double_t > &state, const Double_t &mass)
const std::string &	getPolicyName ()
virtual TMatrixT< Double_t >	residualVector (const GFAbsTrackRep *stateVector, const TMatrixT< Double_t > &state, const GFDetPlane &d)
	Calculate residual with respect to a track representation. More...
virtual TMatrixT< Double_t >	residualVector (const GFAbsTrackRep *stateVector, const TMatrixT< Double_t > &state, const GFDetPlane &d, const GFDetPlane &dPrev, const double &mass)
TMatrixT< Double_t >	getRawHitCov () const
	Get raw hit covariances. More...
TMatrixT< Double_t >	getRawHitCoord () const
	Get raw hit coordinates. More...
virtual void	Print (std::ostream &out=std::cout) const
	Print raw hit coordinates. More...
int	getNparHit ()

Protected Attributes
HitPolicy	fPolicy
TMatrixT< Double_t >	fHitCoord
	Vector of raw coordinates of hit. More...
TMatrixT< Double_t >	fHitCov
	Covariance of raw hit coordinates. More...

Static Private Attributes
static const int	NparHitRep = 3

Detailed Description

Definition at line 17 of file PointHit.h.

Constructor & Destructor Documentation

genf::PointHit::PointHit

(

)

Definition at line 23 of file PointHit.cxx.

Referenced by clone().

  : SpacepointRecoHit(NparHitRep)
{}

genf::PointHit::PointHit	(	TVector3	point,
		double	res
	)

Definition at line 27 of file PointHit.cxx.

References d, genf::GFAbsRecoHit::fHitCoord, and genf::GFAbsRecoHit::fHitCov.

  : SpacepointRecoHit(NparHitRep){

  fHitCov[0][0] = res*res;
  fHitCov[1][1] = res*res;
  fHitCov[2][2] = res*res;
  GFDetPlane d;

  fHitCoord[0][0] = point.X();
  fHitCoord[1][0] = point.Y();
  fHitCoord[2][0] = point.Z();

}

genf::PointHit::PointHit	(	TVector3	point,
		std::vector< double > &	ref
	)

Definition at line 42 of file PointHit.cxx.

References d, genf::GFAbsRecoHit::fHitCoord, and genf::GFAbsRecoHit::fHitCov.

  : SpacepointRecoHit(NparHitRep){

  //assert (res.size()==4);

  fHitCov[0][0] = res.at(0);
  fHitCov[1][1] = res.at(1);
  fHitCov[2][2] = res.at(3);
  fHitCov[1][2] = res.at(2); // yz cov element.
  fHitCov[2][1] = res.at(2); // yz cov element.

  GFDetPlane d;

  fHitCoord[0][0] = point.X();
  fHitCoord[1][0] = point.Y();
  fHitCoord[2][0] = point.Z();
}

genf::PointHit::PointHit	(	TVector3	point,
		TVector3	res
	)

genf::PointHit::~PointHit ( )

virtual

Definition at line 20 of file PointHit.cxx.

{}

Member Function Documentation

genf::GFAbsRecoHit * genf::PointHit::clone ( )

virtual

Get clone of this object.

Virtual abstract method. Has to be implemented by inherting classes. Creates a deep copy of this object. Ownership is trandsferred to the caller!

Implements genf::GFAbsRecoHit.

Definition at line 60 of file PointHit.cxx.

References PointHit().

                                     {
  return new PointHit(*this);
}

template<class HitPolicy >

virtual const GFDetPlane& genf::GFRecoHitIfc< HitPolicy >::getDetPlane ( GFAbsTrackRep *  rep )

inlinevirtualinherited

Returns the detector plane object for this hit and a given track representation.

The actutal code for this method depends on the hit geometry and is implemented in the HitPolicy

See also

PlanarHitPolicy

SpacepointHitPolicy

WirepointHitPolicy

Implements genf::GFAbsRecoHit.

Definition at line 82 of file GFRecoHitIfc.h.

{return fPolicy.detPlane(this,rep);}

template<class HitPolicy >

virtual TMatrixT<Double_t> genf::GFRecoHitIfc< HitPolicy >::getHitCoord ( const GFDetPlane &  plane, const GFDetPlane &  planePrev  )

inlinevirtualinherited

Get hit coordinates in a specific detector plane.

Implementation in the HitPolicy

Implements genf::GFAbsRecoHit.

Definition at line 88 of file GFRecoHitIfc.h.

  {return fPolicy.hitCoord(this,plane,planePrev);}

template<class HitPolicy >

virtual TMatrixT<Double_t> genf::GFRecoHitIfc< HitPolicy >::getHitCoord ( const GFDetPlane &  plane )

inlinevirtualinherited

Implements genf::GFAbsRecoHit.

Definition at line 90 of file GFRecoHitIfc.h.

  {return fPolicy.hitCoord(this,plane);}

template<class HitPolicy >

virtual TMatrixT<Double_t> genf::GFRecoHitIfc< HitPolicy >::getHitCov ( const GFDetPlane &  plane )

inlinevirtualinherited

Get hit covariances in a specific detector plane.

Implementation in the HitPolicy

Implements genf::GFAbsRecoHit.

Definition at line 97 of file GFRecoHitIfc.h.

  {return fPolicy.hitCov(this,plane);}

template<class HitPolicy >

virtual TMatrixT<Double_t> genf::GFRecoHitIfc< HitPolicy >::getHitCov ( const GFDetPlane &  plane, const GFDetPlane &  planePrev, const TMatrixT< Double_t > &  state, const Double_t &  mass  )

inlinevirtualinherited

Implements genf::GFAbsRecoHit.

Definition at line 99 of file GFRecoHitIfc.h.

  {return fPolicy.hitCov(this,plane,planePrev, state, mass);}

TMatrixT< Double_t > genf::PointHit::getHMatrix ( const GFAbsTrackRep *  stateVector )

virtual

Get transformation matrix. Transformation between hit coordinates and track representation coordinates.

This is a virtual abstract method which has to be implemented in the child classes.

In general there is a linear transformation between the coordinate system of the hit (which is defined by the detector plane) and the coordinates of the track representation in that plane. In the most simple case the track representation has 5 parameters (space + momentum) while a hit usually has less (one to three space coordinates).

The transformation matrix is then simply projecting out the space-components of the track representation.

Its dimensions are NxM. Where N is the number of dimensions of the hit in the detector plane (usually 2 or 1) and M is the dimension of the track representation.

In this method a hit has to define with which track representations it can work together. It should be the only point where this explicit coordination is necessary.

For example code see implementing classes below:

Implements genf::GFAbsRecoHit.

Definition at line 141 of file PointHit.cxx.

{
  if ((dynamic_cast<const genf::GeaneTrackRep2*>(stateVector) != NULL) ||
      (dynamic_cast<const genf::RKTrackRep*>(stateVector) != NULL)){
    //I know, since this is the same everytime, it could be done in the
    //the constructor, but I do it here anyway, to make clear that in the
    //case of several track-reps per hit, it would have to be done here
   //    fHMatrix.ResizeTo(NparHitRep,5);


    //   TMatrixT<Double_t> HMatrix(2,6);
    // WTF? Making this 2,5 not 2,6. EC, 3-Jan-2011.
    TMatrixT<Double_t> HMatrix(2,5);

    HMatrix[0][0] = 0.;
    HMatrix[0][1] = 0.;
    HMatrix[0][2] = 0.;
    HMatrix[0][3] = 1.;
    HMatrix[0][4] = 0.;
    //HMatrix[0][5] = 0.;

    HMatrix[1][0] = 0.;
    HMatrix[1][1] = 0.;
    HMatrix[1][2] = 0.;
    HMatrix[1][3] = 0.;
    HMatrix[1][4] = 1.;
    //HMatrix[1][5] = 0.;

    return HMatrix;
  }
  else if(dynamic_cast<const genf::SlTrackRep*>(stateVector)){
    TMatrixT<Double_t> HMatrix(2,4);
    
    HMatrix[0][0] = 1.;
    HMatrix[0][1] = 0.;
    HMatrix[0][2] = 0.;
    HMatrix[0][3] = 0.;

    
    HMatrix[1][0] = 0.;
    HMatrix[1][1] = 1.;
    HMatrix[1][2] = 0.;
    HMatrix[1][3] = 0.;
    return HMatrix;
  }

 else {
   std::cerr << "PointHit can only handle state"
             << " vectors of type GeaneTrackRep2 -> abort" << std::endl;
   throw;
 }
 
}

TMatrixT< Double_t > genf::PointHit::getHMatrix ( const GFAbsTrackRep *  stateVector, const Double_t &  betac, const Double_t &  dist  )

virtual

Implements genf::GFAbsRecoHit.

Definition at line 66 of file PointHit.cxx.

{
  if ((dynamic_cast<const genf::GeaneTrackRep2*>(stateVector) != NULL) ||
      (dynamic_cast<const genf::RKTrackRep*>(stateVector) != NULL)){
    //I know, since this is the same everytime, it could be done in the
    //the constructor, but I do it here anyway, to make clear that in the
    //case of several track-reps per hit, it would have to be done here
   //    fHMatrix.ResizeTo(NparHitRep,5);


    //   TMatrixT<Double_t> HMatrix(2,6);
    // WTF? Making this 2,5 not 2,6. EC, 3-Jan-2011.
    //    TMatrixT<Double_t> HMatrix(2,5);
    TMatrixT<Double_t> HMatrix(5,5);

    Double_t C = 0.0136/betac*sqrt(dist/14.0)*(1+0.038*log(dist/14.0)); // EC, 2-Jan-2012.

    HMatrix[0][0] = C; 
    HMatrix[0][1] = 0.;
    HMatrix[0][2] = 0.;
    HMatrix[0][3] = 0.;
    HMatrix[0][4] = 0.;

    HMatrix[1][0] = 0.; 
    HMatrix[1][1] = 1.;
    HMatrix[1][2] = 0.;
    HMatrix[1][3] = 0.;
    HMatrix[1][4] = 0.;

    HMatrix[2][0] = 0.; 
    HMatrix[2][1] = 0.;
    HMatrix[2][2] = 1.;
    HMatrix[2][3] = 0.;
    HMatrix[2][4] = 0.;

    HMatrix[3][0] = 0.; 
    HMatrix[3][1] = 0.;
    HMatrix[3][2] = 0.;
    HMatrix[3][3] = 1.;
    HMatrix[3][4] = 0.;

    HMatrix[4][0] = 0.; 
    HMatrix[4][1] = 0.;
    HMatrix[4][2] = 0.;
    HMatrix[4][3] = 0.;
    HMatrix[4][4] = 1.;


    return HMatrix;
  }
  else if(dynamic_cast<const genf::SlTrackRep*>(stateVector)){
    TMatrixT<Double_t> HMatrix(2,4);
    
    HMatrix[0][0] = 1.;
    HMatrix[0][1] = 0.;
    HMatrix[0][2] = 0.;
    HMatrix[0][3] = 0.;

    
    HMatrix[1][0] = 0.;
    HMatrix[1][1] = 1.;
    HMatrix[1][2] = 0.;
    HMatrix[1][3] = 0.;
    return HMatrix;
  }

 else {
   std::cerr << "PointHit can only handle state"
             << " vectors of type GeaneTrackRep2 -> abort" << std::endl;
   throw;
 }
 
}

int genf::GFAbsRecoHit::getNparHit ( )

inlineinherited

Definition at line 238 of file GFAbsRecoHit.h.

References genf::GFAbsRecoHit::fNparHit.

{return fNparHit;}

template<class HitPolicy >

const std::string& genf::GFRecoHitIfc< HitPolicy >::getPolicyName ( )

inlinevirtualinherited

Reimplemented from genf::GFAbsRecoHit.

Definition at line 102 of file GFRecoHitIfc.h.

{return fPolicy.getName();}

TMatrixT<Double_t> genf::GFAbsRecoHit::getRawHitCoord ( ) const

inlineinherited

Get raw hit coordinates.

Definition at line 182 of file GFAbsRecoHit.h.

References genf::GFAbsRecoHit::clone(), genf::GFAbsRecoHit::fHitCoord, genf::GFAbsRecoHit::getDetPlane(), genf::GFAbsRecoHit::getHitCoord(), and genf::GFAbsRecoHit::getHitCov().

Referenced by genf::GFWireHitPolicy::checkPlane(), genf::GFWirepointHitPolicy::checkPlane(), genf::GFWirepointHitPolicy::detPlane(), genf::GFWireHitPolicy::detPlane(), genf::GFSpacepointHitPolicy::detPlane(), genf::GFWirepointHitPolicy::hitCoord(), genf::GFWireHitPolicy::hitCoord(), genf::GFSpacepointHitPolicy::hitCoord(), genf::GFPlanarHitPolicy::hitCoord(), and genf::GFKalman::processHit().

{return fHitCoord;}

TMatrixT<Double_t> genf::GFAbsRecoHit::getRawHitCov ( ) const

inlineinherited

Get raw hit covariances.

Definition at line 177 of file GFAbsRecoHit.h.

References genf::GFAbsRecoHit::fHitCov.

Referenced by genf::GFWireHitPolicy::hitCov(), genf::GFWirepointHitPolicy::hitCov(), genf::GFSpacepointHitPolicy::hitCov(), and genf::GFPlanarHitPolicy::hitCov().

{return fHitCov;}

virtual void genf::GFAbsRecoHit::Print ( std::ostream &  out = std::cout ) const

inlinevirtualinherited

Print raw hit coordinates.

Definition at line 233 of file GFAbsRecoHit.h.

References genf::GFAbsRecoHit::getPolicyName(), and genf::PrintROOTobject().

    { PrintROOTobject(out, fHitCoord); }

virtual TMatrixT<Double_t> genf::GFAbsRecoHit::residualVector ( const GFAbsTrackRep *  stateVector, const TMatrixT< Double_t > &  state, const GFDetPlane &  d  )

inlinevirtualinherited

Calculate residual with respect to a track representation.

Returns the N-dimensional residual of this vector to a given track representation.

This method is not doing any extrapolation. But it creates the necessary detector plane object. See GFAbsRecoHit::getGFDetPlane

Parameters

stateVector	pointer to track representation - used to synchronize with the track repesentation
state	parameter vector of the track representation

See also

setHMatrix

getGFDetPlane

Definition at line 147 of file GFAbsRecoHit.h.

References genf::GFAbsRecoHit::getHitCoord(), and genf::GFAbsRecoHit::getHMatrix().

Referenced by genf::GFKalman::getChi2Hit(), genf::GFTrack::getResiduals(), and genf::GFKalman::processHit().

                                                                 {
    std::cout << "GFAbsRecoHit::residualVector(3args) Not correctly Using theta -- multiple scattering -- information !!! Fix this if you really want to use getChi2Hit" << std::endl;
    TMatrixT<Double_t> H = getHMatrix(stateVector);
    return ( getHitCoord(d) - (H*state ));
  }

virtual TMatrixT<Double_t> genf::GFAbsRecoHit::residualVector ( const GFAbsTrackRep *  stateVector, const TMatrixT< Double_t > &  state, const GFDetPlane &  d, const GFDetPlane &  dPrev, const double &  mass  )

inlinevirtualinherited

Definition at line 155 of file GFAbsRecoHit.h.

References beta, genf::GFAbsRecoHit::getHitCoord(), genf::GFAbsRecoHit::getHMatrix(), and genf::GFDetPlane::getO().

  {
    Double_t dist = (d.getO()-dPrev.getO()).Mag();
    //    Double_t mass = 0.104; // close enough for muons, pions, I think.
    Double_t mom = fabs(1.0/state[0][0]);
    Double_t beta = mom/sqrt(mass*mass+mom*mom);
    if (std::isnan(dist) || dist<0.2) dist=0.2; // don't allow 0s here.
    if (std::isnan(beta) || beta<0.04) beta=0.04;

    TMatrixT<Double_t> H = getHMatrix(stateVector,beta,dist);
    return ( getHitCoord(d,dPrev) - (H*state ));
  }

Member Data Documentation

TMatrixT<Double_t> genf::GFAbsRecoHit::fHitCoord

protectedinherited

Vector of raw coordinates of hit.

Definition at line 80 of file GFAbsRecoHit.h.

Referenced by genf::GFAbsRecoHit::getRawHitCoord(), and PointHit().

TMatrixT<Double_t> genf::GFAbsRecoHit::fHitCov

protectedinherited

Covariance of raw hit coordinates.

Definition at line 83 of file GFAbsRecoHit.h.

Referenced by genf::GFAbsRecoHit::getRawHitCov(), and PointHit().

template<class HitPolicy >

HitPolicy genf::GFRecoHitIfc< HitPolicy >::fPolicy

protectedinherited

Definition at line 64 of file GFRecoHitIfc.h.

const int genf::PointHit::NparHitRep = 3

staticprivate

Definition at line 37 of file PointHit.h.

---

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

• larreco/v06_64_02/source/larreco/Genfit/PointHit.h
• larreco/v06_64_02/source/larreco/Genfit/PointHit.cxx