TPCParametersImpl.h

#ifndef GEAR_TPCParametersImpl_H
#define GEAR_TPCParametersImpl_H 1

#include "gear/TPCParameters.h"
#include <map>
#include "gearimpl/GearParametersImpl.h"
#include "gearimpl/GlobalPadIndex.h"



namespace gear {

class TPCModule;

class TPCParametersImpl :  public TPCParameters,    public GearParametersImpl {
    
  public: 
    
    TPCParametersImpl( double maxDriftLength=-1., int coordinateType=-1 ) ;
    TPCParametersImpl( const TPCParametersImpl & ) ;

    TPCParametersImpl & operator = ( const TPCParametersImpl & ) ;

//     /** Copy constructor for base parameters */
//     TPCParametersImpl( GearParameters& gearParameters ) ;

    virtual ~TPCParametersImpl() ;
    
    virtual const TPCModule& getModule(int moduleID) const; 

    virtual int getNModules() const;

    virtual const TPCModule & getNearestModule(double c0, double c1) const;

    virtual const TPCModule & getNearestModule(double c0, double c1, double z) const;

    virtual double getMaxDriftLength() const;

    virtual bool isInsideModule(double c0, double c1) const;

    virtual bool isInsideModule(double c0, double c1, double z) const;

    virtual bool isInsidePad(double c0, double c1) const;

    virtual bool isInsidePad(double c0, double c1, double z) const;

    virtual GlobalPadIndex getNearestPad(double c0, double c1) const;

    virtual GlobalPadIndex getNearestPad(double c0, double c1, double z) const;

    virtual const std::vector<double>  & getPlaneExtent() const;

    virtual int getCoordinateType() const;

    virtual const std::vector<TPCModule *> & getModules() const;
   
    virtual void  addModule( TPCModule * TPCModule ) ;

    virtual const PadRowLayout2D& getPadLayout() const ;
    
    virtual double getDriftVelocity() const ;
    
    virtual double getReadoutFrequency() const ;

    virtual void setPadLayout( PadRowLayout2D * padLayout ) ;
    
    virtual void  setMaxDriftLength( double maxDriftLength ) ;
    
    virtual void  setDriftVelocity( double driftVelocity ) ;
    
    virtual void  setReadoutFrequency( double readoutFrequency ) ;

protected:

  // one vector with all the modules, and one for each half tpc
    std::vector<TPCModule *>  _TPCModules{}, _modulesPositiveHalfTPC{}, _modulesNegativeHalfTPC{};

  const TPCModule & getNearestModule(double c0, double c1, 
                                     std::vector<TPCModule *> const & modulesVector) const;
  
  bool isInsideModule(double c0, double c1, std::vector<TPCModule *> const & modulesVector) const;

  bool isInsidePad(double c0, double c1, std::vector<TPCModule *> const & modulesVector ) const;

  GlobalPadIndex getNearestPad(double c0, double c1,
                               std::vector<TPCModule *> const & modulesVector ) const;
 

  double _maxDriftLength {};

  int _coordinateType{};

  double _cathodePosition{};
    
  std::map<int,int> _moduleIDMap {};

  std::vector<double> _planeExtent {};

  double _driftVelocity {};

  void copy_and_assign(const  TPCParametersImpl & );

  void cleanup();

  // Internal function to avoid code duplication
  // Recalculation of the z postion is necessary when a module is added and when the 
  // maximal drift distance is changed.
  void setCathodePosition();

}; // class

} // namespace gear


#endif /* ifndef GEAR_TPCParametersImpl_H */