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

#include <GFWirepointHitPolicy.h>

Inheritance diagram for GFWirepointHitPolicy:

Public Member Functions
	GFWirepointHitPolicy ()

const GFDetPlane &	detPlane (GFAbsRecoHit , GFAbsTrackRep )
	Get detector plane. More...

TMatrixT< double >	hitCoord (GFAbsRecoHit *, const GFDetPlane &)
	Hit coordinates in detector plane. More...

TMatrixT< double >	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 ()

	ClassDef (GFWirepointHitPolicy, 1)

Private Attributes
GFDetPlane	fDetPlane

double	fMaxdistance

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

Detailed Description

policy class for hits in wire detectors (STT and DCH) which can 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 8 coordinates: w_x1, w_y1, w_z1, w_x2, w_y2, w_z2, rdrift, zreco 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) and zreco = coordinate along the wire (in the plane reference frame, v coordinate).

Definition at line 52 of file GFWirepointHitPolicy.h.

Constructor & Destructor Documentation

GFWirepointHitPolicy::GFWirepointHitPolicy ( )

Definition at line 36 of file GFWirepointHitPolicy.cxx.

36 : fMaxdistance(1.E50) {;}

GFWirepointHitPolicy::fMaxdistance

double fMaxdistance

Definition: GFWirepointHitPolicy.h:86

virtual GFWirepointHitPolicy::~GFWirepointHitPolicy ( )

inlinevirtual

Definition at line 75 of file GFWirepointHitPolicy.h.

75 {;}

Member Function Documentation

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

Check if the detector plane is valid.

Definition at line 71 of file GFWirepointHitPolicy.cxx.

References CAMath::Abs(), fabs(), GFAbsRecoHit::getRawHitCoord(), and GFDetPlane::getV().

Referenced by hitCoord(), and hitCov().

 {
   // raw x1, y1, z1, x2, y2, z2, rdrift, zreco
   TMatrixT<double> 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;
     }
 }

GFWirepointHitPolicy::ClassDef	(	GFWirepointHitPolicy	,
		1
	)

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

Get detector plane.

Definition at line 95 of file GFWirepointHitPolicy.cxx.

References GFDetPlane::dist(), Double_t, GFAbsTrackRep::extrapolateToLine(), fabs(), fDetPlane, fMaxdistance, GFAbsRecoHit::getRawHitCoord(), CAMath::Sqrt(), and x.

 {
 
   TMatrixT<double> 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) {
     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 = (wire1 + wire2) * 0.5;
 
   
   fDetPlane = GFDetPlane(O, U, V);
   
   return fDetPlane;
 }

double GFWirepointHitPolicy::getMaxDistance ( )

inline

Definition at line 77 of file GFWirepointHitPolicy.h.

References fMaxdistance.

77 {return fMaxdistance;}

GFWirepointHitPolicy::fMaxdistance

double fMaxdistance

Definition: GFWirepointHitPolicy.h:86

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

inline

Definition at line 80 of file GFWirepointHitPolicy.h.

References fPolicyName.

80 {return fPolicyName;}

GFWirepointHitPolicy::fPolicyName

static const std::string fPolicyName

Definition: GFWirepointHitPolicy.h:82

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

Hit coordinates in detector plane.

Definition at line 39 of file GFWirepointHitPolicy.cxx.

References checkPlane(), and GFAbsRecoHit::getRawHitCoord().

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

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

Hit covariances in detector plane.

Definition at line 54 of file GFWirepointHitPolicy.cxx.

References checkPlane(), and GFAbsRecoHit::getRawHitCov().

 {
   checkPlane(hit, plane);
 
   TMatrixT<double> returnCov(2,2);
   TMatrixT<double> 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;
 }