policy class for hits in wire detectors (STT and DCH) which do not measure the coordinate along the wire More...

#include <GFWireHitPolicy.h>

Inheritance diagram for genf::GFWireHitPolicy:

Public Member Functions
	GFWireHitPolicy ()

const GFDetPlane &	detPlane (GFAbsRecoHit , GFAbsTrackRep )
	Get detector plane Calls GFAbsTrackRep::extrapolateToLine for POCA. 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	~GFWireHitPolicy ()

double	getMaxDistance ()

void	setMaxDistance (double d)

const std::string &	getName ()

Private Attributes
GFDetPlane	fDetPlane

double	fMaxdistance

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

Detailed Description

policy class for hits in wire detectors (STT and DCH) which do not measure the coordinate along the wire

Author: Christian Höppner (Technische Universität München, original author); Lia Lavezzi (INFN Pavia, original author); Sebastian Neubert (Technische Universität München, original author)

This policy is not valid for any kind of plane orientation choice: to use it you MUST choose a plane described by u and v axes with v coincident with the wire (and u orthogonal to it, obviously). The hit will be described by 7 coordinates: w_x1, w_y1, w_z1, w_x2, w_y2, w_z2, rdrift where w_ji (with j = x, y, z and i = 1, 2) are the wire extremities coordinates; rdrift = distance from the wire (u coordinate in the plane)

Definition at line 54 of file GFWireHitPolicy.h.

Constructor & Destructor Documentation

genf::GFWireHitPolicy::GFWireHitPolicy ( )

Definition at line 38 of file GFWireHitPolicy.cxx.

38 :fMaxdistance(1.E50) {;}

genf::GFWireHitPolicy::fMaxdistance

double fMaxdistance

Definition: GFWireHitPolicy.h:90

virtual genf::GFWireHitPolicy::~GFWireHitPolicy ( )

inlinevirtual

Definition at line 78 of file GFWireHitPolicy.h.

78 {;}

Member Function Documentation

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

Check if the detector plane is valid.

Definition at line 69 of file GFWireHitPolicy.cxx.

 {
   // raw x1, y1, z1, x2, y2, z2, rdrift, zreco
   TMatrixT<Double_t> rC = hit->getRawHitCoord();
 
   assert(rC.GetNrows()==7);
 
   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 << "GFWireHitPolicy: plane not valid!!" << std::endl;
     }
 }

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

Get detector plane Calls GFAbsTrackRep::extrapolateToLine for POCA.

Definition at line 93 of file GFWireHitPolicy.cxx.

 {
 
   TMatrixT<Double_t> x=hit->getRawHitCoord();
   assert(x.GetNrows()==7);
   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) {
     GFException exc("distance poca-wire > maxdistance", __LINE__,__FILE__);
     throw exc;
   }
 
   // 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) {
     GFException exc("fromwiretoextr*wiredirection > 1e-3", __LINE__,__FILE__);
     throw exc;
   }
 
   TVector3 U;
   U = fromwiretoextr;
   TVector3 V;
   V = wiredirection;
   U.SetMag(1.);
   V.SetMag(1.);
 
   TVector3 O = poca_onwire;
 
   fDetPlane = GFDetPlane(O, U, V);
 
   return fDetPlane;
 }

double genf::GFWireHitPolicy::getMaxDistance ( )

inline

Definition at line 80 of file GFWireHitPolicy.h.

80 {return fMaxdistance;}

genf::GFWireHitPolicy::fMaxdistance

double fMaxdistance

Definition: GFWireHitPolicy.h:90

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

inline

Definition at line 83 of file GFWireHitPolicy.h.

83 {return fPolicyName;}

genf::GFWireHitPolicy::fPolicyName

static const std::string fPolicyName

Definition: GFWireHitPolicy.h:85

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

Hit coordinates in detector plane.

Definition at line 41 of file GFWireHitPolicy.cxx.

 {
   TMatrixT<Double_t> returnMat(1,1);
 
   checkPlane(hit, plane);
 
   // raw x1, y1, z1, x2, y2, z2, rdrift
   TMatrixT<Double_t> rC = hit->getRawHitCoord();
 
   returnMat[0][0] = rC[6][0];
   return returnMat;
 }

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

Hit covariances in detector plane.

Definition at line 55 of file GFWireHitPolicy.cxx.

 {
   checkPlane(hit, plane);
 
   TMatrixT<Double_t> returnCov(1,1);
   TMatrixT<Double_t> rawCov = hit->getRawHitCov();
 
   returnCov[0][0] = rawCov[6][6];
 
   return  returnCov;
 }