genf::GFWireHitPolicy Class Reference

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 36 of file GFWireHitPolicy.cxx.

:fMaxdistance(1.E50) {;}

virtual genf::GFWireHitPolicy::~GFWireHitPolicy ( )

inlinevirtual

Definition at line 78 of file GFWireHitPolicy.h.

{;}

Member Function Documentation

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

Check if the detector plane is valid.

Definition at line 67 of file GFWireHitPolicy.cxx.

References e, genf::GFAbsRecoHit::getRawHitCoord(), and genf::GFDetPlane::getV().

Referenced by hitCoord(), and hitCov().

{
  // 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 91 of file GFWireHitPolicy.cxx.

References genf::GFDetPlane::dist(), e, genf::GFAbsTrackRep::extrapolateToLine(), fDetPlane, fMaxdistance, genf::GFAbsRecoHit::getRawHitCoord(), and x.

{

  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.

References fMaxdistance.

{return fMaxdistance;}

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

inline

Definition at line 83 of file GFWireHitPolicy.h.

References fPolicyName.

{return fPolicyName;}

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

Hit coordinates in detector plane.

Definition at line 39 of file GFWireHitPolicy.cxx.

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

{
  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 53 of file GFWireHitPolicy.cxx.

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

{
  checkPlane(hit, plane);

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

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

  return  returnCov;
}

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

inline

Definition at line 81 of file GFWireHitPolicy.h.

References d, and fMaxdistance.

{fMaxdistance=d;}

Member Data Documentation

GFDetPlane genf::GFWireHitPolicy::fDetPlane

private

Definition at line 89 of file GFWireHitPolicy.h.

Referenced by detPlane().

double genf::GFWireHitPolicy::fMaxdistance

private

Definition at line 90 of file GFWireHitPolicy.h.

Referenced by detPlane(), getMaxDistance(), and setMaxDistance().

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

staticprivate

Definition at line 85 of file GFWireHitPolicy.h.

Referenced by getName().

---

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

• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/Genfit/GFWireHitPolicy.h
• /cvmfs/larsoft.opensciencegrid.org/products/larreco/v07_10_02/source/larreco/Genfit/GFWireHitPolicy.cxx