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(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