ARICHGeoMasterVolume.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 <string>
#include <Math/Vector3D.h>
#include <Math/Rotation3D.h>
#include <Math/RotationX.h>
#include <Math/RotationY.h>
#include <Math/RotationZ.h>
 
namespace Belle2 {
  class ARICHGeoMasterVolume: public ARICHGeoBase {
 
  public:
 
    ARICHGeoMasterVolume()
    {}
 
    ~ARICHGeoMasterVolume()
    {
      /*if(m_rotation) delete m_rotation;
      if(m_rotationInverse) delete m_rotationInverse;
      if(m_translation) delete m_translation;
      */
    }
 
    bool isConsistent() const override;
 
    void print(const std::string& title = "ARICH Master Volume geometry parameters") const override;
 
    void setPlacement(double x, double y, double z, double rx, double ry, double rz);
 
    void setVolume(double innerR, double outerR, double length, const std::string& material) {m_innerR = innerR; m_outerR = outerR; m_length = length; m_material = material;};
 
    ROOT::Math::XYZVector getPosition() const {return ROOT::Math::XYZVector(m_x / s_unit, m_y / s_unit, m_z / s_unit);}
 
    ROOT::Math::Rotation3D getRotation() const
    {
      ROOT::Math::Rotation3D rot;
      ROOT::Math::RotationX rotX(m_rx);
      ROOT::Math::RotationY rotY(m_ry);
      ROOT::Math::RotationZ rotZ(m_rz);
      rot *= rotZ * rotY * rotX;
      return rot;
    }
 
    double getRotationX() const {return m_rx;}
 
    double getRotationY() const {return m_ry;}
 
    double getRotationZ() const {return m_rz;}
 
    double getInnerRadius() const {return m_innerR / s_unit;};
 
    double getOuterRadius() const {return m_outerR / s_unit;};
 
    double getLength() const {return m_length / s_unit;};
 
    const std::string& getMaterial() const {return m_material;}
 
    ROOT::Math::XYZVector pointToGlobal(const ROOT::Math::XYZVector& point) const;
 
    ROOT::Math::XYZVector momentumToGlobal(const ROOT::Math::XYZVector& momentum) const;
 
    ROOT::Math::XYZVector pointToLocal(const ROOT::Math::XYZVector& point) const;
 
    ROOT::Math::XYZVector momentumToLocal(const ROOT::Math::XYZVector& momentum) const;
 
 
  private:
    // position in global BelleII frame
    double m_x = 0;  
    double m_y = 0;  
    double m_z = 0;  
    // rotations in global BelleII frame
    double m_rx = 0; 
    double m_ry = 0; 
    double m_rz = 0; 
    double m_innerR = 0; 
    double m_outerR = 0; 
    double m_length = 0; 
    std::string m_material; 
    mutable ROOT::Math::Rotation3D* m_rotation = nullptr ; 
    mutable ROOT::Math::Rotation3D* m_rotationInverse = nullptr; 
    mutable ROOT::Math::XYZVector*  m_translation = nullptr; 
    ClassDefOverride(ARICHGeoMasterVolume, 2); 
  };
 
} // end namespace Belle2