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 throw (Exception, std::exception );
    double getRadiationLength(double x, double y, double z) const throw (Exception, std::exception );
    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 throw (Exception, std::exception );
    void interpolateOnGrid(std::pair<double,double> &result, double x, double y, double z) const;
  }; // class
} // namespace gear
#endif /* ifndef  GEAR_MATERIALMAPFACTORY_H */