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<stdexcept> // std::logic_error
#include"GFAbsFinitePlane.h"
#include "TObject.h"

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

namespace genf {

  class GFDetPlane : public TObject {
  public:
    
    // Constructors/Destructors ---------
    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&);
    // 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(genf::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(std::ostream& out = std::cout) const;

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

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

    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;

    genf::GFAbsFinitePlane* fFinitePlane;

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

    void sane(); // ensures orthnormal coordinates
    
  private:
    virtual void Print(Option_t*) const
      { throw std::logic_error(std::string(__func__) + "::Print(Option_t*) not available"); }

    //public:
    //ClassDef(GFDetPlane,2)

  };

  bool operator==(const genf::GFDetPlane&, const genf::GFDetPlane&);
  bool operator!=(const genf::GFDetPlane&, const genf::GFDetPlane&);
  
} // namespace genf

#endif