genf::GFDetPlane Class Reference

#include <GFDetPlane.h>

Inheritance diagram for genf::GFDetPlane:

Public Member Functions
	GFDetPlane (genf::GFAbsFinitePlane *finite=NULL)
	GFDetPlane (const TVector3 &o, const TVector3 &u, const TVector3 &v, genf::GFAbsFinitePlane *finite=NULL)
	GFDetPlane (const TVector3 &o, const TVector3 &n, genf::GFAbsFinitePlane *finite=NULL)
virtual	~GFDetPlane ()
	GFDetPlane (const genf::GFDetPlane &)
GFDetPlane &	operator= (const genf::GFDetPlane &)
TVector3	getO () const
TVector3	getU () const
TVector3	getV () const
void	set (const TVector3 &o, const TVector3 &u, const TVector3 &v)
void	setO (const TVector3 &o)
void	setO (double, double, double)
void	setU (const TVector3 &u)
void	setU (double, double, double)
void	setV (const TVector3 &v)
void	setV (double, double, double)
void	setUV (const TVector3 &u, const TVector3 &v)
void	setON (const TVector3 &o, const TVector3 &n)
void	setFinitePlane (genf::GFAbsFinitePlane *finite)
TVector3	getNormal () const
void	setNormal (TVector3 n)
void	setNormal (double, double, double)
void	setNormal (const double &theta, const double &phi)
TVector2	project (const TVector3 &x) const
	projecting a direction onto the plane: More...
TVector2	LabToPlane (const TVector3 &x) const
	transform from Lab system into plane More...
TVector3	toLab (const TVector2 &x) const
	transform from plane coordinates to lab system More...
TVector3	dist (const TVector3 &point) const
TVector2	straightLineToPlane (const TVector3 &point, const TVector3 &dir) const
	gives u,v coordinates of the intersection point of a straight line with plane More...
void	Print (std::ostream &out=std::cout) const
void	getGraphics (double mesh, double length, TPolyMarker3D **pl, TPolyLine3D **plLine, TPolyLine3D **u, TPolyLine3D **v, TPolyLine3D **n=NULL)
	for poor attempts of making an event display. There is a lot of room for improvements. More...
double	distance (TVector3 &) const
double	distance (double, double, double) const
bool	inActive (const TVector3 &point, const TVector3 &dir) const
	intersect in the active area? C.f. GFAbsFinitePlane More...
bool	inActive (double u, double v) const
	inActive methods refer to finite plane. C.f. GFAbsFinitePlane More...
bool	inActive (const TVector2 &v) const
	inActive methods refer to finite plane. C.f. GFAbsFinitePlane More...
void	sane ()

Public Attributes
TVector3	fO
TVector3	fU
TVector3	fV
genf::GFAbsFinitePlane *	fFinitePlane

Private Member Functions
virtual void	Print (Option_t *) const

Friends
bool	operator== (const GFDetPlane &, const GFDetPlane &)
bool	operator!= (const GFDetPlane &, const GFDetPlane &)
	returns NOT == More...

Detailed Description

Definition at line 65 of file GFDetPlane.h.

Constructor & Destructor Documentation

genf::GFDetPlane::GFDetPlane

(

genf::GFAbsFinitePlane *

finite = NULL

)

Definition at line 44 of file GFDetPlane.cxx.

  :fFinitePlane(finite)
{
  static TRandom3 r(0);
  fO.SetXYZ(0.,0.,0.);
  fU.SetXYZ(r.Uniform(),r.Uniform(),0.);
  fV.SetXYZ(r.Uniform(),r.Uniform(),0.);
  sane();
}

genf::GFDetPlane::GFDetPlane	(	const TVector3 &	o,
		const TVector3 &	u,
		const TVector3 &	v,
		genf::GFAbsFinitePlane *	finite = NULL
	)

Definition at line 36 of file GFDetPlane.cxx.

:fO(o),fU(u),fV(v),fFinitePlane(finite)
{
  sane();
}

genf::GFDetPlane::GFDetPlane	(	const TVector3 &	o,
		const TVector3 &	n,
		genf::GFAbsFinitePlane *	finite = NULL
	)

Definition at line 54 of file GFDetPlane.cxx.

  :fO(o),fFinitePlane(finite){
  setNormal(n);
}

genf::GFDetPlane::~GFDetPlane ( )

virtual

Definition at line 61 of file GFDetPlane.cxx.

                           {
  if(fFinitePlane!=NULL) delete fFinitePlane;
}

genf::GFDetPlane::GFDetPlane

(

const genf::GFDetPlane &

rhs

)

Definition at line 65 of file GFDetPlane.cxx.

                                                   : TObject(rhs) {
  if(rhs.fFinitePlane != NULL) fFinitePlane = rhs.fFinitePlane->clone();
  else fFinitePlane = NULL;
  fO = rhs.fO;
  fU = rhs.fU;
  fV = rhs.fV;
}

Member Function Documentation

TVector3 genf::GFDetPlane::dist

(

const TVector3 &

point

)

const

Definition at line 214 of file GFDetPlane.cxx.

{
  TVector2 p=LabToPlane(x);
  TVector3 xplane=toLab(p);
  return xplane-x;
}

double genf::GFDetPlane::distance

(

TVector3 &

v

)

const

Definition at line 354 of file GFDetPlane.cxx.

                                                {
  double s = (v - fO)*fU;
  double t = (v - fO)*fV;
  TVector3 distanceVector = v - fO - (s*fU) - (t*fV);
  return distanceVector.Mag();
}

double genf::GFDetPlane::distance	(	double	x,
		double	y,
		double	z
	)		const

Definition at line 360 of file GFDetPlane.cxx.

                                                               {
  TVector3 v(x,y,z);
  double s = (v - fO)*fU;
  double t = (v - fO)*fV;
  TVector3 distanceVector = v - fO - (s*fU) - (t*fV);
  return distanceVector.Mag();
}

void genf::GFDetPlane::getGraphics	(	double	mesh,
		double	length,
		TPolyMarker3D **	pl,
		TPolyLine3D **	plLine,
		TPolyLine3D **	u,
		TPolyLine3D **	v,
		TPolyLine3D **	n = NULL
	)

for poor attempts of making an event display. There is a lot of room for improvements.

Definition at line 294 of file GFDetPlane.cxx.

                                                                                                                                                     {
  *pl = new TPolyMarker3D(21*21,24);
  (*pl)->SetMarkerSize(0.1);
  (*pl)->SetMarkerColor(kBlue);
  int nI=10;
  int nJ=10;
  *plLine = new TPolyLine3D(5);

  {
        TVector3 linevec;
        int i,j;
        i=-1*nI;j=-1*nJ;
        linevec=(fO+(mesh*i)*fU+(mesh*j)*fV);
        (*plLine)->SetPoint(0,linevec.X(),linevec.Y(),linevec.Z());
        i=-1*nI;j=1*nJ;
        linevec=(fO+(mesh*i)*fU+(mesh*j)*fV);
        (*plLine)->SetPoint(0,linevec.X(),linevec.Y(),linevec.Z());
        i=1*nI;j=-1*nJ;
        linevec=(fO+(mesh*i)*fU+(mesh*j)*fV);
        (*plLine)->SetPoint(2,linevec.X(),linevec.Y(),linevec.Z());
        i=1*nI;j=1*nJ;
        linevec=(fO+(mesh*i)*fU+(mesh*j)*fV);
        (*plLine)->SetPoint(1,linevec.X(),linevec.Y(),linevec.Z());
        i=-1*nI;j=-1*nJ;
        linevec=(fO+(mesh*i)*fU+(mesh*j)*fV);
        (*plLine)->SetPoint(4,linevec.X(),linevec.Y(),linevec.Z());

  }
  for (int i=-1*nI;i<=nI;++i){
        for (int j=-1*nJ;j<=nJ;++j){
          TVector3 vec(fO+(mesh*i)*fU+(mesh*j)*fV);
          int id=(i+10)*21+j+10;
          (*pl)->SetPoint(id,vec.X(),vec.Y(),vec.Z());
        }
  }


  *u = new TPolyLine3D(2);
  (*u)->SetPoint(0,fO.X(),fO.Y(),fO.Z());
  (*u)->SetPoint(1,fO.X()+length*fU.X(),fO.Y()+length*fU.Y(),fO.Z()+length*fU.Z());
  (*u)->SetLineWidth(2);
  (*u)->SetLineColor(kGreen);


  *v = new TPolyLine3D(2);
  (*v)->SetPoint(0,fO.X(),fO.Y(),fO.Z());
  (*v)->SetPoint(1,fO.X()+length*fV.X(),fO.Y()+length*fV.Y(),fO.Z()+length*fV.Z());
  (*v)->SetLineWidth(2);
  (*v)->SetLineColor(kRed);

  if(n!=NULL){
    *n = new TPolyLine3D(2);
    TVector3 _n=getNormal();
    (*n)->SetPoint(0,fO.X(),fO.Y(),fO.Z());
    (*n)->SetPoint(1,fO.X()+length*_n.X(),fO.Y()+length*_n.Y(),fO.Z()+length*_n.Z());
    (*n)->SetLineWidth(2);
    (*n)->SetLineColor(kBlue);
  }
}

TVector3 genf::GFDetPlane::getNormal

(

)

const

Definition at line 145 of file GFDetPlane.cxx.

{
  TVector3 result=fU.Cross(fV);
  result.SetMag(1.);
  return result;
}

TVector3 genf::GFDetPlane::getO ( ) const

inline

Definition at line 81 of file GFDetPlane.h.

{return fO;}

TVector3 genf::GFDetPlane::getU ( ) const

inline

Definition at line 82 of file GFDetPlane.h.

{return fU;}

TVector3 genf::GFDetPlane::getV ( ) const

inline

Definition at line 83 of file GFDetPlane.h.

{return fV;}

bool genf::GFDetPlane::inActive ( const TVector3 &  point, const TVector3 &  dir  ) const

inline

intersect in the active area? C.f. GFAbsFinitePlane

Definition at line 138 of file GFDetPlane.h.

                                                                   {
      return this->inActive( this->straightLineToPlane(point,dir));
    }

bool genf::GFDetPlane::inActive ( double  u, double  v  ) const

inline

inActive methods refer to finite plane. C.f. GFAbsFinitePlane

Definition at line 143 of file GFDetPlane.h.

                                           {
      if(fFinitePlane==NULL) return true;
      return fFinitePlane->inActive(u,v);
    }

bool genf::GFDetPlane::inActive ( const TVector2 &  v ) const

inline

inActive methods refer to finite plane. C.f. GFAbsFinitePlane

Definition at line 149 of file GFDetPlane.h.

                                          {
      return inActive(v.X(),v.Y());
    }

TVector2 genf::GFDetPlane::LabToPlane

(

const TVector3 &

x

)

const

transform from Lab system into plane

Definition at line 196 of file GFDetPlane.cxx.

{
  TVector3 d=x-fO;
  return project(d);
}

genf::GFDetPlane & genf::GFDetPlane::operator=

(

const genf::GFDetPlane &

rhs

)

Definition at line 72 of file GFDetPlane.cxx.

                                                                  {
  if (this == &rhs) return *this;
  if(fFinitePlane!=NULL) {
    delete fFinitePlane;
  }
  if(rhs.fFinitePlane != NULL){
    fFinitePlane = rhs.fFinitePlane->clone();
  }
  else{
    fFinitePlane = NULL;
  }
  fO = rhs.fO;
  fU = rhs.fU;
  fV = rhs.fV;
  return *this;
}

void genf::GFDetPlane::Print

(

std::ostream &

out = std::cout

)

const

Definition at line 247 of file GFDetPlane.cxx.

{
  out<<"GFDetPlane: "
           <<"O("<<fO.X()<<","<<fO.Y()<<","<<fO.Z()<<") "
           <<"u("<<fU.X()<<","<<fU.Y()<<","<<fU.Z()<<") "
           <<"v("<<fV.X()<<","<<fV.Y()<<","<<fV.Z()<<") "
           <<"n("<<getNormal().X()<<","<<getNormal().Y()<<","<<getNormal().Z()<<") "
                   <<std::endl;
  out << fFinitePlane << std::endl;
  if(fFinitePlane!=NULL) fFinitePlane->Print(out);

}

virtual void genf::GFDetPlane::Print ( Option_t *  ) const

inlineprivatevirtual

Definition at line 169 of file GFDetPlane.h.

      { throw std::logic_error(std::string(__func__) + "::Print(Option_t*) not available"); }

TVector2 genf::GFDetPlane::project

(

const TVector3 &

x

)

const

projecting a direction onto the plane:

Definition at line 188 of file GFDetPlane.cxx.

{
  Double_t xfU=fU*x;
  Double_t xfV=fV*x;
  return TVector2(xfU,xfV);
}

void genf::GFDetPlane::sane

(

)

Definition at line 224 of file GFDetPlane.cxx.

                    {
  if (fU==fV)
    throw GFException("genf::GFDetPlane::sane() sanity check failed", __LINE__, __FILE__).setFatal();
  // ensure unit vectors
  fU.SetMag(1.);
  fV.SetMag(1.);
  // ensure orthogonal system
  // fV should be reached by
  // rotating fU around _n in a counterclockwise direction

  TVector3 n=getNormal();
  fV=n.Cross(fU);

  TVector3 v=fU;
  v.Rotate(TMath::Pi()*0.5,n);
  TVector3 null=v-fV;
  if (null.Mag() >= 1E-6)
    throw GFException("genf::GFDetPlane::sane(): non orthogonal!", __LINE__, __FILE__).setFatal();

}

void genf::GFDetPlane::set	(	const TVector3 &	o,
		const TVector3 &	u,
		const TVector3 &	v
	)

Definition at line 90 of file GFDetPlane.cxx.

{
  fO=o;fU=u;fV=v;
  sane();
}

void genf::GFDetPlane::setFinitePlane ( genf::GFAbsFinitePlane *  finite )

inline

Optionally, set the finite plane definition. This is most important for avoiding fake intersection points in fitting of loopers. This should be implemented for silicon detectors most importantly.

Definition at line 102 of file GFDetPlane.h.

{fFinitePlane=finite;}

void genf::GFDetPlane::setNormal

(

TVector3

n

)

Definition at line 163 of file GFDetPlane.cxx.

                                   {
  n.SetMag(1.);
  if( fabs(n.X()) > 0.1 ){
        fU.SetXYZ(1./n.X()*(-1.*n.Y()-1.*n.Z()),1.,1.);
        fU.SetMag(1.);
  }
  else {
        if(fabs(n.Y()) > 0.1){
          fU.SetXYZ(1.,1./n.Y()*(-1.*n.X()-1.*n.Z()),1.);
          fU.SetMag(1.);
        }
        else{
          fU.SetXYZ(1.,1.,1./n.Z()*(-1.*n.X()-1.*n.Y()));
          fU.SetMag(1.);
        }
  }
  fV=n.Cross(fU);
}

void genf::GFDetPlane::setNormal	(	double	X,
		double	Y,
		double	Z
	)

Definition at line 158 of file GFDetPlane.cxx.

                                                   {
  TVector3 N(X,Y,Z);
  setNormal(N);
}

void genf::GFDetPlane::setNormal	(	const double &	theta,
		const double &	phi
	)

Definition at line 182 of file GFDetPlane.cxx.

                                                                    {
  TVector3 n(TMath::Sin(theta)*TMath::Cos(phi),TMath::Sin(theta)*TMath::Sin(phi),TMath::Cos(theta));
  setNormal(n);
}

void genf::GFDetPlane::setO

(

const TVector3 &

o

)

Definition at line 99 of file GFDetPlane.cxx.

{
  fO=o;
  sane();
}

void genf::GFDetPlane::setO	(	double	X,
		double	Y,
		double	Z
	)

Definition at line 105 of file GFDetPlane.cxx.

{
  fO.SetXYZ(X,Y,Z);
  sane();
}

void genf::GFDetPlane::setON	(	const TVector3 &	o,
		const TVector3 &	n
	)

Definition at line 152 of file GFDetPlane.cxx.

                                                             {
  fO = o;
  setNormal(n);
}

void genf::GFDetPlane::setU

(

const TVector3 &

u

)

Definition at line 112 of file GFDetPlane.cxx.

{
  fU=u;
  sane();
}

void genf::GFDetPlane::setU	(	double	X,
		double	Y,
		double	Z
	)

Definition at line 118 of file GFDetPlane.cxx.

{
  fU.SetXYZ(X,Y,Z);
  sane();
}

void genf::GFDetPlane::setUV	(	const TVector3 &	u,
		const TVector3 &	v
	)

Definition at line 137 of file GFDetPlane.cxx.

{
  fU=u;
  fV=v;
  sane();
}

void genf::GFDetPlane::setV

(

const TVector3 &

v

)

Definition at line 125 of file GFDetPlane.cxx.

{
  fV=v;
  sane();
}

void genf::GFDetPlane::setV	(	double	X,
		double	Y,
		double	Z
	)

Definition at line 131 of file GFDetPlane.cxx.

{
  fV.SetXYZ(X,Y,Z);
  sane();
}

TVector2 genf::GFDetPlane::straightLineToPlane	(	const TVector3 &	point,
		const TVector3 &	dir
	)		const

gives u,v coordinates of the intersection point of a straight line with plane

Definition at line 368 of file GFDetPlane.cxx.

                                                                                            {
  TVector3 dirNorm(dir);
  dirNorm.SetMag(1.);
  TVector3 normal = getNormal();
  double dirTimesN = dirNorm*normal;
  if(fabs(dirTimesN)<1.E-6){//straight line is parallel to plane, so return infinity
    //doesnt parallel mean that they cross in infinity ;-)
    return TVector2(1.E100,1.E100);
  }
  double t = 1/dirTimesN * ((fO-point)*normal);
  return project(point - fO + t * dirNorm);
}

TVector3 genf::GFDetPlane::toLab

(

const TVector2 &

x

)

const

transform from plane coordinates to lab system

Definition at line 203 of file GFDetPlane.cxx.

{
  TVector3 d(fO);
  d+=x.X()*fU;
  d+=x.Y()*fV;
  return d;
}

Friends And Related Function Documentation

bool operator!= ( const GFDetPlane &  , const GFDetPlane &    )

friend

returns NOT ==

bool operator== ( const GFDetPlane &  , const GFDetPlane &    )

friend

this operator is called very often in Kalman filtering. It checks equality of planes by comparing the 9 double values that define them.

Member Data Documentation

genf::GFAbsFinitePlane* genf::GFDetPlane::fFinitePlane

Definition at line 162 of file GFDetPlane.h.

TVector3 genf::GFDetPlane::fO

Definition at line 157 of file GFDetPlane.h.

TVector3 genf::GFDetPlane::fU

Definition at line 159 of file GFDetPlane.h.

TVector3 genf::GFDetPlane::fV

Definition at line 160 of file GFDetPlane.h.

---

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

• larreco/larreco/Genfit/GFDetPlane.h
• larreco/larreco/Genfit/GFDetPlane.cxx