GFDetPlane.h

Go to the documentation of this file.

/* Copyright 2008-2010, Technische Universitaet Muenchen,
   Authors: Christian Hoeppner & Sebastian Neubert

   This file is part of GENFIT.

   GENFIT is free software: you can redistribute it and/or modify
   it under the terms of the GNU Lesser General Public License as published
   by the Free Software Foundation, either version 3 of the License, or
   (at your option) any later version.

   GENFIT is distributed in the hope that it will be useful,
   but WITHOUT ANY WARRANTY; without even the implied warranty of
   MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
   GNU Lesser General Public License for more details.

   You should have received a copy of the GNU Lesser General Public License
   along with GENFIT.  If not, see <http://www.gnu.org/licenses/>.
*/
// Description:
//      Detector plane - a geometric object
#ifndef GFDETPLANE_H
#define GFDETPLANE_H

#include"GFAbsFinitePlane.h"

#include "TObject.h"

#include "TVector3.h"
#include "TPolyMarker3D.h"
#include "TPolyLine3D.h"

class GFDetPlane : public TObject {
public:

  // Constructors/Destructors ---------
  GFDetPlane(GFAbsFinitePlane* finite=NULL);
  GFDetPlane(const TVector3& o,
             const TVector3& u,
             const TVector3& v,
             GFAbsFinitePlane* finite=NULL);
  GFDetPlane(const TVector3& o,
             const TVector3& n,
             GFAbsFinitePlane* finite=NULL);
  virtual ~GFDetPlane();
  GFDetPlane(const GFDetPlane&);
  GFDetPlane& operator=(const GFDetPlane&);
  // Accessors -----------------------
  TVector3 getO() const {return fO;}
  TVector3 getU() const {return fU;}
  TVector3 getV() const {return fV;}

  // Modifiers -----------------------
  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(GFAbsFinitePlane* finite){fFinitePlane=finite;}

  // Operations ----------------------
  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;
  TVector2 LabToPlane(const TVector3& x) const;
  TVector3 toLab(const TVector2& x) const;
  // get vector from point to plane (normal)
  TVector3 dist(const TVector3& point) const;
  
  TVector2 straightLineToPlane(const TVector3& point,const TVector3& dir) const;


  void Print() const;

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

  friend bool operator== (const GFDetPlane& lhs, const GFDetPlane& rhs);
  friend bool operator!= (const GFDetPlane& lhs, const GFDetPlane& rhs);

  double distance(TVector3&) const;
  double distance(double,double,double) const;


  bool inActive(const TVector3& point, const TVector3& dir) const{
    return this->inActive( this->straightLineToPlane(point,dir));
  }

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

  bool inActive(const TVector2& v) const{
    return inActive(v.X(),v.Y());
  }
  
private:

  // Private Data Members ------------
  // origin
  TVector3 fO;
  // Vectors spanning the plane
  TVector3 fU;
  TVector3 fV;

  GFAbsFinitePlane* fFinitePlane;

  // Private Methods -----------------

  void sane(); // ensures orthnormal coordinates
public:
  ClassDef(GFDetPlane,2)

};

#endif