VXDGeometry.h

#ifndef VXDGeometry_h
#define VXDGeometry_h 1

//#include "lcio.h"
#include <map>
#include <vector>
//#include <gsl/gsl_rng.h>

#include <gear/GearMgr.h>
#include "CLHEP/Vector/TwoVector.h"

struct VXDLadder{
  double phi{};   // phi of ladder - rotation araound z-axis
  gear::Vector3D trans{}; // translation after rotation
//   CLHEP::Hep2Vector p0 ;  // 'left' end of ladder in r-phi
//   CLHEP::Hep2Vector p1 ;  // 'right' end of ladder in r-phi
//   CLHEP::Hep2Vector u  ;  // unit vector along ladder in r-phi
};
typedef std::vector< std::vector< VXDLadder > > VXDLadders ;

struct VXDLayer{
  double rMin{};
  double rMax{};
  double length{};
  double width{};
  double thickness{};
  double gap{};
  double ladderArea{};
  int nLadders{};
};

typedef std::vector< VXDLayer >  VXDLayers ;



class VXDGeometry  {
  
public:

  VXDGeometry(const VXDGeometry&) = delete;
  VXDGeometry& operator=(const VXDGeometry&) = delete;  

  VXDGeometry(gear::GearMgr* gearMgr) ;

  
  std::pair<int,int> getLadderID( gear::Vector3D labPos, int layerID=-1 ) ;
  
  gear::Vector3D lab2LadderPos( gear::Vector3D labPos, int layerID, int ladderID) ;
  
  gear::Vector3D ladder2LabPos( gear::Vector3D ladderPos, int layerID, int ladderID) ;

  gear::Vector3D lab2LadderDir( gear::Vector3D labPos, int layerID, int ladderID) ;
  
  gear::Vector3D ladder2LabDir( gear::Vector3D ladderPos, int layerID, int ladderID) ;

  void test() ;


protected:

  void init() ;
  VXDGeometry(){}

  gear::GearMgr* _gearMgr{};
  VXDLadders _vxdLadders{};
  VXDLayers  _vxdLayers{};

} ;

#endif