CsiGeometryPar.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.                  *
 **************************************************************************/
#ifndef CSIGEOMETRYPAR_H
#define CSIGEOMETRYPAR_H
 
#include <vector>
#include <G4String.hh>
#include <G4Material.hh>
 
#include <CLHEP/Geometry/Transform3D.h>
#include <CLHEP/Vector/ThreeVector.h>
#include <TVector3.h>
 
typedef int CsiIdentifier ;
 
typedef CLHEP::HepRotation RotationMatrix;
typedef CLHEP::Hep3Vector ThreeVector;
 
typedef HepGeom::Transform3D Transform3D;
//typedef HepGeom::Rotate3D  Rotate3D;
//typedef HepGeom::RotateX3D RotateX3D;
typedef HepGeom::RotateY3D RotateY3D;
typedef HepGeom::RotateZ3D RotateZ3D;
 
typedef HepGeom::Translate3D Translate3D;
//typedef HepGeom::TranslateX3D TranslateX3D;
//typedef HepGeom::TranslateY3D TranslateY3D;
//typedef HepGeom::TranslateZ3D TranslateZ3D;
 
 
namespace Belle2 {
  namespace csi {
 
 
    class CsiGeometryPar {
 
    public:
 
      CsiGeometryPar();
 
      virtual ~CsiGeometryPar();
 
 
      static CsiGeometryPar* Instance();
 
      void clear();
 
      void Print(const int cid, int debuglevel = 80);
 
      void PrintAll(int debuglevel = 80);
 
 
      void read();
 
      int CsiVolNameToCellID(const G4String VolumeName);//Mapping from VolumeName to Crystal CellID
 
      G4Material* GetMaterial(int cid);
 
      ThreeVector GetPosition(int cid) { return m_Position.at(cid); };
 
      ThreeVector GetOrientation(int cid) { return m_Orientation.at(cid); };
 
      TVector3 GetPositionTV3(int cid) { return m_PositionTV3.at(cid); };
 
      TVector3 GetOrientationTV3(int cid) { return m_OrientationTV3.at(cid); };
 
      TVector3 ConvertToTVector3(ThreeVector _hepTV)
      {
        TVector3 pos(_hepTV.x(), _hepTV.y(), _hepTV.z());
        return pos;
      };
 
      int GetEnclosureID(int cid) { return m_BoxID.at(cid); };
 
      int GetSlotID(int cid) { return m_SlotID.at(cid); };
 
      double GetMaterialProperty(int cid, const char* propertyname);
 
      double GetTauFast(int cid) { return GetMaterialProperty(cid, "FASTTIMECONSTANT"); };
 
      double GetTauSlow(int cid) { return GetMaterialProperty(cid, "SLOWTIMECONSTANT"); };
 
 
    private:
 
 
      int m_cellID;
 
      std::vector<int> m_thetaID;
 
      std::vector<int> m_phiID;
 
      std::vector<int> m_BoxID;
 
      std::vector<int> m_SlotID;
 
      std::vector<ThreeVector> m_Position;
      std::vector<TVector3> m_PositionTV3;
 
      std::vector<ThreeVector> m_Orientation;
      std::vector<TVector3> m_OrientationTV3;
 
      static CsiGeometryPar* m_B4CsiGeometryParDB;
    };
 
  } // end of namespace csi
} // end of namespace Belle2
 
 
#endif