#ifndef L1MaterialInfo_h
#define L1MaterialInfo_h

#include "CbmL1Def.h"
#include <vector>
using std::vector;

class L1MaterialInfo{
 public:
 L1MaterialInfo(): thick(0), RL(0),  RadThick(0), logRadThick(0) {};
   
  fvec thick, RL,  RadThick, logRadThick;
} _fvecalignment;


class L1Material{
 public:

 L1Material() : table(), NBins(0), RMax(0.), iD(0.) {}

  vector<vector<float> > table;
  
  // static const int NBins = 100; // TODO file?
  // static const float RMax = 60.f;
  // static const float iD = 0.5*NBins/60.f;//RMax!;

  void SetBins( int n, float r ) {
    NBins = n;
    RMax = r;
    iD = 0.5*NBins/RMax;
  }

  float GetRadThick( float x, float y ) {
    x = ( x < RMax && x >= -RMax ) ? x : 0;
    y = ( y < RMax && y >= -RMax ) ? y : 0;
    int i = static_cast<int>( (x + RMax)*iD );
    int j = static_cast<int>(  (y + RMax)*iD );
    i = ( i < NBins && i >= 0 ) ? i : NBins/2;
    j = ( j < NBins && j >= 0 ) ? j : NBins/2;
    return table[ i ][ j ];
  }


  fvec GetRadThick( fvec x, fvec y ) {
    fvec r;
    for ( int i = 0; i < fvecLen; i++ )
      r[i] = GetRadThick(x[i], y[i]);
    return r;
  }

  int NBins;
  float RMax;
  float iD;

} _fvecalignment;

#endif