ARICHGlobalAlignment.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <arich/dbobjects/ARICHGeoBase.h>
#include <arich/dbobjects/ARICHPositionElement.h>
#include <string>
#include <TVector3.h>
#include <TRotation.h>
 
namespace Belle2 {
  class ARICHGlobalAlignment: public ARICHGeoBase {
 
  public:
 
    ARICHGlobalAlignment()
    {}
 
    ARICHGlobalAlignment(const ARICHGlobalAlignment& align): ARICHGeoBase()
    {
      *this = align;
      m_rotation = 0;
      m_rotationInverse = 0;
      m_translation = 0;
    }
 
 
    ARICHGlobalAlignment& operator=(const ARICHGlobalAlignment& align)
    {
      if (this != &align) {
        ARICHGeoBase::operator=(align);
        m_alignPars = align.getAlignmentElement();
        if (m_rotation) delete m_rotation;
        if (m_rotationInverse) delete m_rotationInverse;
        if (m_translation) delete m_translation;
        m_rotation = 0;
        m_rotationInverse = 0;
        m_translation = 0;
      }
      return *this;
    }
 
    ~ARICHGlobalAlignment()
    {
      if (m_rotation) delete m_rotation;
      if (m_rotationInverse) delete m_rotationInverse;
      if (m_translation) delete m_translation;
    }
 
    void setAlignmentElement(const ARICHPositionElement& align)
    {
      m_alignPars = align;
    }
 
    const ARICHPositionElement& getAlignmentElement() const
    {
      return m_alignPars;
    }
 
 
    void print(const std::string& title = "ARICH global alignment parameters") const;
 
 
    const TVector3& getTranslation() const {if (!m_translation) setTransformation(); return *m_translation;}
 
    const TRotation& getRotation() const
    {
      if (!m_rotation) setTransformation();
      return *m_rotation;
    }
 
    TVector3 pointToGlobal(const TVector3& point) const;
    TVector3 momentumToGlobal(const TVector3& momentum) const;
    TVector3 pointToLocal(const TVector3& point) const;
    TVector3 momentumToLocal(const TVector3& momentum) const;
 
  private:
 
    void setTransformation() const;
 
    ARICHPositionElement m_alignPars;
 
    mutable TRotation* m_rotation = 0 ;
    mutable TRotation* m_rotationInverse = 0;
    mutable TVector3*  m_translation = 0;
 
    ClassDef(ARICHGlobalAlignment, 1); 
  };
 
} // end namespace Belle2