ARICHGeoMasterVolume.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Luka Santelj                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <arich/dbobjects/ARICHGeoBase.h>
#include <string>
#include <TVector3.h>
#include <TRotation.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;};
 
    TVector3 getPosition() const {return TVector3(m_x / s_unit, m_y / s_unit, m_z / s_unit);}
 
    TRotation getRotation() const
    {
      TRotation rot;
      rot.RotateX(m_rx).RotateY(m_ry).RotateZ(m_rz);
      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;}
 
    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:
    // 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  TRotation* m_rotation = 0 ;
    mutable TRotation* m_rotationInverse = 0;
    mutable TVector3*  m_translation = 0;
 
    ClassDefOverride(ARICHGeoMasterVolume, 1); 
  };
 
} // end namespace Belle2