ZPlanarParametersImpl.h

// -*- C++ -*-
#ifndef GEAR_ZPlanarParametersImpl_H
#define GEAR_ZPlanarParametersImpl_H 1

#include "gear/ZPlanarParameters.h"
#include "gearimpl/GearParametersImpl.h"
#include "gearimpl/ZPlanarLayerLayoutImpl.h"


namespace gear {

class ZPlanarLayerLayoutImpl;

class ZPlanarParametersImpl : public GearParametersImpl, public ZPlanarParameters {

public: 
  ZPlanarParametersImpl( int type, double shellInnerRadius, double shellOuterRadius, double shellHalfLength, double shellGap, double shellRadLength ) ;

  // Destructor.
  virtual ~ZPlanarParametersImpl() { /* nop */; }
  
  // static constants enums( for the type of layout )
  
  static const int CCD    = 1 ;
  static const int CMOS   = 2 ;
  static const int HYBRID = 3 ;

  virtual void addLayer(int nLadders, double phi0,
                        double ladderDistance, double ladderOffset, double ladderThickness,
                        double ladderLength, double ladderWidth, double ladderRadLength,
                        double sensitiveDistance, double sensitiveOffset, double sensitiveThickness,
                        double sensitiveLength, double sensitiveWidth, double sensitiveRadLength )
  {
    _layer.addLayer( nLadders, phi0,
                     ladderDistance, ladderOffset, ladderThickness, ladderLength, ladderWidth, ladderRadLength,
                     sensitiveDistance, sensitiveOffset, sensitiveThickness, sensitiveLength, sensitiveWidth, sensitiveRadLength ) ;
    return ;
  }
    
      
  virtual const ZPlanarLayerLayout & getZPlanarLayerLayout() const { return _layer ; }
  
  virtual int getType() const { return _type ; }
  
  virtual double getShellHalfLength() const { return _shellHalfLength ; }

  virtual double getShellGap() const { return _shellGap ; }
  
  virtual double getShellInnerRadius() const { return _shellInnerRadius ; }
  
  virtual double getShellOuterRadius() const { return _shellOuterRadius ; }
  
  virtual double getShellRadLength() const { return _shellRadLength ; }
  
  virtual bool isPointInLadder(Vector3D p) const {
    return isPointInLadder( p, false );
  }
  
  virtual bool isPointInSensitive(Vector3D p, SensorID* sensorID=0 ) const {
    return isPointInLadder( p, true , sensorID ) ;
  }
  
  virtual Vector3D distanceToNearestLadder(Vector3D p) const {
    return distanceToNearestLadder( p, false ) ;
  }

  virtual Vector3D distanceToNearestSensitive(Vector3D p) const {
    return distanceToNearestLadder( p, true ) ;
  }

  virtual Vector3D intersectionLadder( Vector3D p, Vector3D v ) const {
    return intersectionLadder( p, v, false ) ;
  }

  virtual Vector3D intersectionSensitive( Vector3D p, Vector3D v ) const {
    return intersectionLadder( p, v, true ) ;
  }

protected:
  
  ZPlanarLayerLayoutImpl _layer{} ;
  
  int _type{} ;
  
  double _shellInnerRadius{} ;
  
  double _shellOuterRadius{} ;

  double _shellHalfLength{} ;

  double _shellGap{} ;

  double _shellRadLength{} ;

private:

  bool isPointInLadder(Vector3D p , bool sensitive = false, SensorID* sensorID=0 ) const ;

  Vector3D distanceToNearestLadder(Vector3D p, bool sensitive = false) const ;

  Vector3D distanceToPlane(Vector3D p, Vector3D r, Vector3D n, Vector3D u, Vector3D v, float minU, float maxU, float minV, float maxV) const ;

  Vector3D intersectionLadder( Vector3D p, Vector3D v, bool sensitive = false) const ;

  Vector3D planeLineIntersection( Vector3D r, Vector3D n, Vector3D linePoint, Vector3D lineDir) const ;

  double confiningRatio( Vector3D p , Vector3D r, Vector3D n, Vector3D u, Vector3D v, float minU, float maxU, float minV, float maxV ) const ;
  
  Vector3D correctToBorderPoint( Vector3D vPlane , Vector3D u, Vector3D v, float minU, float maxU, float minV, float maxV ) const ;

  bool isEqual( double valueOne , double valueTwo ) const ;

  bool isEqual( Vector3D p1 , Vector3D p2 ) const ;

  bool differsLess( double valueOne , double valueTwo ) const ;

  
  double correctPhiRange( double Phi ) const ;

  double getPhiPoint( Vector3D p ) const ;

}; // class

} // namespace gear

#endif /* ifndef GEAR_ZPlanarPARAMETERS_H */