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 <string>

#include "GFAbsFinitePlane.h"
#include "TObject.h"

#include "RtypesCore.h"
#include "TVector2.h"
#include "TVector3.h"

class TPolyMarker3D;
class TPolyLine3D;

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