PndPlane Class Reference

#include <PndPlane.h>

Public Member Functions
	PndPlane (TVector3 b, TVector3 d1, TVector3 d2)
	PndPlane (PndLine line, Int_t layer)
virtual	~PndPlane ()
PndLine	getIntersection (PndPlane &p)
void	Print ()
TVector3	getBase ()
TVector3	getDir1 ()
TVector3	getDir2 ()
TVector3	getIntersection (PndLine l)

Private Attributes
TVector3	base
TVector3	dir1
TVector3	dir2

Detailed Description

Definition at line 21 of file PndPlane.h.

Constructor & Destructor Documentation

PndPlane::PndPlane ( TVector3  b, TVector3  d1, TVector3  d2  )

inline

Definition at line 23 of file PndPlane.h.

: base(b),dir1(d1),dir2(d2) {}

PndPlane::PndPlane	(	PndLine	line,
		Int_t	layer
	)

Definition at line 11 of file PndPlane.cxx.

References base, CAMath::Cos(), dir1, dir2, PndLine::getP1(), PndLine::getP2(), and CAMath::Sin().

                                           {
        base = line.getP1();
        dir1 = line.getP2()-line.getP1();
        if(layer%4==1){
                dir2 = TVector3(TMath::Sin(-5*TMath::DegToRad()),TMath::Cos(-5*TMath::DegToRad()),0);
        }else if(layer%4==2){
                dir2 = TVector3(TMath::Sin(5*TMath::DegToRad()),TMath::Cos(5*TMath::DegToRad()),0);
        } else {
                dir2 = TVector3(0,1,0);
        }
}

PndPlane::~PndPlane ( )

virtual

Definition at line 23 of file PndPlane.cxx.

                    {
        // TODO Auto-generated destructor stub
}

Member Function Documentation

TVector3 PndPlane::getBase ( )

inline

Definition at line 28 of file PndPlane.h.

{return base;}

TVector3 PndPlane::getDir1 ( )

inline

Definition at line 29 of file PndPlane.h.

{return dir1;}

TVector3 PndPlane::getDir2 ( )

inline

Definition at line 30 of file PndPlane.h.

{return dir2;}

PndLine PndPlane::getIntersection

(

PndPlane &

p

)

Definition at line 28 of file PndPlane.cxx.

References base, dir1, dir2, Double_t, and i.

Referenced by PndTrackCombiner::combine(), PndWayFollower::getBestSkewedApproxTo(), and PndLine::getPerpendicular().

                                            {
        //create the direction vector
        TVector3 norm1 = dir1.Cross(dir2);
        Double_t d1 = norm1.Dot(base);
        TVector3 norm2 = p.dir1.Cross(p.dir2);
        Double_t d2 = norm2.Dot(p.base);
        TVector3 dirLine = norm1.Cross(norm2);
        //create the base vector

        //init array
        Double_t arr[2][3];
        arr[0][0]=norm1[0];
        arr[0][1]=norm1[1];
        arr[0][2]=d1-norm1[2]*base[2];
        arr[1][0]=norm2[0];
        arr[1][1]=norm2[1];
        arr[1][2]=d2-norm2[2]*base[2];

        //swap rows
        if(arr[0][0]==0){
                for(int i=0;i<3;i++){
                        //Double_t temp = arr[0][i]; //[R.K.03/2017] unused variable
                        arr[0][i]=arr[1][i];
                        arr[1][i]=arr[0][i];
                }
        }

        Double_t fac = -arr[1][0]/arr[0][0];
        for(int i=0;i<3;i++){
                arr[1][i]+=fac*arr[0][i];
        }

        Double_t x2 = arr[1][2]/arr[1][1];
        Double_t x1 = (arr[0][2]-x2*arr[0][1])/arr[0][0];

        //create the line
        TVector3 lineBase;
        lineBase.SetXYZ(x1,x2,base[2]);
        PndLine l(lineBase,dirLine);
        return l;
}

TVector3 PndPlane::getIntersection

(

PndLine

l

)

Definition at line 78 of file PndPlane.cxx.

References CAMath::Abs(), base, dir1, dir2, Double_t, PndLine::getDir(), PndLine::getP1(), i, and max().

                                           {
        //create the Matrix
        TVector3 lineDir = l.getDir();
        TVector3 lineBase = l.getP1();
        Double_t matrix[3][4];
        for(int i=0;i<3;i++){
                matrix[i][0]=dir1[i];
                matrix[i][1]=dir2[i];
                matrix[i][2]=-lineDir[i];
                matrix[i][3]=lineBase[i]-base[i];
        }
        //solve the Matrix
        for(int i=0;i<3;i++){
                if(TMath::Abs(matrix[i][i])<0.000000001){//equal 0
                        //search max element
                        Int_t max = i;
                        for(int j=i+1;j<3;j++){
                                if(TMath::Abs(matrix[j][i])>TMath::Abs(matrix[max][i]))max=j;
                        }
                        //no greater element found
                        if(max==i){
                                cout << "\033[1;31mNo Intersection from plane and line\033[0m" << endl;
                                return TVector3(0,0,0);
                        }
                        //swap lines
                        for(int j=i;j<4;j++){
                                Double_t temp = matrix[i][j];
                                matrix[i][j] = matrix[max][j];
                                matrix[max][j] = temp;
                        }
                }
                for(int j=i+1;j<3;j++){
                        Double_t fac = -matrix[j][i]/matrix[i][i];
                        for(int k=i;k<4;k++){
                                matrix[j][k]+=(fac*matrix[i][k]);
                        }
                }
        }
        //create the result
        Double_t lamp = matrix[2][3]/matrix[2][2];
        return lineBase+lamp*lineDir;
}

void PndPlane::Print

(

)

Definition at line 70 of file PndPlane.cxx.

References base, dir1, and dir2.

                    {
        cout << "Plane:";
        cout << "(" << base[0] << "," << base[1] << "," << base[2] << ")";
        cout << "+a*(" << dir1[0] << "," << dir1[1] << "," << dir1[2] << ")";
        cout << "+b*(" << dir2[0] << "," << dir2[1] << "," << dir2[2] << ")";
}

Member Data Documentation

TVector3 PndPlane::base

private

Definition at line 33 of file PndPlane.h.

Referenced by getIntersection(), PndPlane(), and Print().

TVector3 PndPlane::dir1

private

Definition at line 34 of file PndPlane.h.

Referenced by getIntersection(), PndPlane(), and Print().

TVector3 PndPlane::dir2

private

Definition at line 35 of file PndPlane.h.

Referenced by getIntersection(), PndPlane(), and Print().

---

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

• PndPlane.h
• PndPlane.cxx