genf::GFWirepointHitPolicy Class Reference

#include <GFWirepointHitPolicy.h>

Inheritance diagram for genf::GFWirepointHitPolicy:

Public Member Functions
	GFWirepointHitPolicy ()
const GFDetPlane &	detPlane (GFAbsRecoHit *, GFAbsTrackRep *)
	Get detector plane. More...
TMatrixT< Double_t >	hitCoord (GFAbsRecoHit *, const GFDetPlane &)
	Hit coordinates in detector plane. More...
TMatrixT< Double_t >	hitCov (GFAbsRecoHit *, const GFDetPlane &)
	Hit covariances in detector plane. More...
void	checkPlane (GFAbsRecoHit *, const GFDetPlane &)
	Check if the detector plane is valid. More...
virtual	~GFWirepointHitPolicy ()
double	getMaxDistance ()
void	setMaxDistance (double d)
const std::string &	getName ()

Private Attributes
GFDetPlane	fDetPlane
double	fMaxdistance

Static Private Attributes
static const std::string	fPolicyName = "GFWirepointHitPolicy"

Detailed Description

Definition at line 55 of file GFWirepointHitPolicy.h.

Constructor & Destructor Documentation

genf::GFWirepointHitPolicy::GFWirepointHitPolicy

(

)

Definition at line 37 of file GFWirepointHitPolicy.cxx.

: fMaxdistance(1.E50) {;}

virtual genf::GFWirepointHitPolicy::~GFWirepointHitPolicy ( )

inlinevirtual

Definition at line 78 of file GFWirepointHitPolicy.h.

{;}

Member Function Documentation

void genf::GFWirepointHitPolicy::checkPlane	(	GFAbsRecoHit *	hit,
		const GFDetPlane &	plane
	)

Check if the detector plane is valid.

Definition at line 72 of file GFWirepointHitPolicy.cxx.

{
  // raw x1, y1, z1, x2, y2, z2, rdrift, zreco
  TMatrixT<Double_t> rC = hit->getRawHitCoord();

  assert(rC.GetNrows()==8);

  TVector3 wire1(rC[0][0], rC[1][0], rC[2][0]);
  TVector3 wire2(rC[3][0], rC[4][0], rC[5][0]);
  TVector3 wiredirection = wire1 - wire2;

  TVector3 vaxis = plane.getV();
  wiredirection.SetMag(1.);
  vaxis.SetMag(1.);

  if(fabs(TMath::Abs(wiredirection.Dot(vaxis)) - 1) > 1e-3)
    {

      std::cout << "GFWirepointHitPolicy: plane not valid!!" << std::endl;
    }
}

const genf::GFDetPlane & genf::GFWirepointHitPolicy::detPlane	(	GFAbsRecoHit *	hit,
		GFAbsTrackRep *	rep
	)

Get detector plane.

Definition at line 96 of file GFWirepointHitPolicy.cxx.

{

  TMatrixT<Double_t> x=hit->getRawHitCoord();
  assert(x.GetNrows()==8);
  TVector3 wire1(x[0][0],x[1][0],x[2][0]);
  TVector3 wire2(x[3][0],x[4][0],x[5][0]);

  //  distance of one (the first) of the wire extremities from the plane
  Double_t d_from_refplane =  fDetPlane.dist(wire1).Mag();
  if(d_from_refplane < 1e-5) return fDetPlane;


  // point of closest approach
  TVector3 poca, poca_onwire, dirInPoca;

  rep->extrapolateToLine(wire1, wire2, poca, dirInPoca, poca_onwire);


  Double_t distance;
  distance = TMath::Sqrt(fabs(((wire1-poca).Mag2()*(wire2-wire1).Mag2()-pow((wire1-poca).Dot(wire2-wire1),2))/(wire2-wire1).Mag2()));

  // check poca inside tube
  if(distance > fMaxdistance) {
    throw GFException("distance poca-wire > maxdistance", __LINE__,__FILE__)/* .setFatal() */;
  }

  // find plane
  // unitary vector along distance
  // poca (on track), poca_onwire (on wire)
  TVector3 fromwiretoextr = poca - poca_onwire;
  fromwiretoextr.SetMag(1.);
  // unitary vector along the wire
  TVector3 wiredirection = wire2 - wire1;
  wiredirection.SetMag(1.);

  // check orthogonality
  if(fabs(fromwiretoextr * wiredirection) > 1e-3) {
    throw GFException("fromwiretoextr*wiredirection > 1e-3", __LINE__,__FILE__)/* .setFatal() */;
  }

  TVector3 U;
  U = fromwiretoextr;
  TVector3 V;
  V = wiredirection;
  U.SetMag(1.);
  V.SetMag(1.);

  TVector3 O = (wire1 + wire2) * 0.5;


  fDetPlane = GFDetPlane(O, U, V);

  return fDetPlane;
}

double genf::GFWirepointHitPolicy::getMaxDistance ( )

inline

Definition at line 80 of file GFWirepointHitPolicy.h.

{return fMaxdistance;}

const std::string& genf::GFWirepointHitPolicy::getName ( )

inline

Definition at line 83 of file GFWirepointHitPolicy.h.

{return fPolicyName;}

TMatrixT< Double_t > genf::GFWirepointHitPolicy::hitCoord	(	GFAbsRecoHit *	hit,
		const GFDetPlane &	plane
	)

Hit coordinates in detector plane.

Definition at line 40 of file GFWirepointHitPolicy.cxx.

{
  TMatrixT<Double_t> returnMat(2,1);

  checkPlane(hit, plane);

  // raw x1, y1, z1, x2, y2, z2, rdrift, zreco
  TMatrixT<Double_t> rC = hit->getRawHitCoord();

  returnMat[0][0] = rC[6][0];
  returnMat[1][0] = rC[7][0];
  return returnMat;
}

TMatrixT< Double_t > genf::GFWirepointHitPolicy::hitCov	(	GFAbsRecoHit *	hit,
		const GFDetPlane &	plane
	)

Hit covariances in detector plane.

Definition at line 55 of file GFWirepointHitPolicy.cxx.

{
  checkPlane(hit, plane);

  TMatrixT<Double_t> returnCov(2,2);
  TMatrixT<Double_t> rawCov = hit->getRawHitCov();

  returnCov[0][0] = rawCov[6][6];
  returnCov[1][0] = rawCov[7][6];
  returnCov[0][1] = rawCov[6][7];
  returnCov[1][1] = rawCov[7][7];

  return  returnCov;
}

void genf::GFWirepointHitPolicy::setMaxDistance ( double  d )

inline

Definition at line 81 of file GFWirepointHitPolicy.h.

{fMaxdistance=d;}

Member Data Documentation

GFDetPlane genf::GFWirepointHitPolicy::fDetPlane

private

Definition at line 88 of file GFWirepointHitPolicy.h.

double genf::GFWirepointHitPolicy::fMaxdistance

private

Definition at line 89 of file GFWirepointHitPolicy.h.

const std::string genf::GFWirepointHitPolicy::fPolicyName = "GFWirepointHitPolicy"

staticprivate

Definition at line 85 of file GFWirepointHitPolicy.h.

---

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

• larreco/larreco/Genfit/GFWirepointHitPolicy.h
• larreco/larreco/Genfit/GFWirepointHitPolicy.cxx