MaterialMap.h

#ifndef GEAR_MATERIALMAP_H
#define GEAR_MATERIALMAP_H 1

#include <string>
#include <vector>

#include "gear/GEAR.h"
#include "gear/GearMgr.h"

namespace gear {
  
  class MaterialMap {
    
  public: 
    MaterialMap(GearMgr *gearMgr,double xmin, double xmax, int nxsteps,
                double ymin, double ymax, int nysteps,
                double zmin, double zmax, int nzsteps, int coord);
    virtual ~MaterialMap() { 
      ;}
    double getInteractionLength(double x, double y, double z) const ;
    double getRadiationLength(double x, double y, double z) const ;
    static const int CARTESIAN = 0 ;
    static const int CYLINDRICAL = 1 ;
    static const int SPHERICAL = 2 ;
 
    
  protected:
    double _xmin{}, _xmax{};
    int _nxsteps{};
    double _ymin{}, _ymax{}; 
    int _nysteps{};
    double _zmin{}, _zmax{};
    int _nzsteps{};
    double _xstep{},_ystep{},_zstep{};
    std::pair<double,double> _values{};
    std::vector< std::vector< std::vector< std::pair<double,double> > > > _myMap{};
    void calculateGridIndex(std::vector<int> & gpos,double x, double y, double z) const ;
    void interpolateOnGrid(std::pair<double,double> &result, double x, double y, double z) const;
  }; // class
} // namespace gear
#endif /* ifndef  GEAR_MATERIALMAPFACTORY_H */