GFDetPlane.h

Go to the documentation of this file.

//****************************************************************************
//*                   This file is part of PandaRoot.                        *
//*                                                                          *
//*            PandaRoot is distributed under the terms of the               *
//*              GNU General Public License (GPL) version 3,                 *
//*                 copied verbatim in the file "LICENSE".                   *
//*                                                                          *
//*  Copyright (C) 2006 - 2024 FAIR GmbH and copyright holders of PandaRoot  *
//*     The copyright holders are listed in the file "COPYRIGHTHOLDERS".     *
//*               The authors are listed in the file "AUTHORS".              *
//****************************************************************************

/* 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 = nullptr);
  GFDetPlane(const TVector3 &o, const TVector3 &u, const TVector3 &v, GFAbsFinitePlane *finite = nullptr);
  GFDetPlane(const TVector3 &o, const TVector3 &n, GFAbsFinitePlane *finite = nullptr);
  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 = nullptr);

  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 == nullptr)
      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