L1MaterialInfo.h

Go to the documentation of this file.

#ifndef L1MaterialInfo_h
#define L1MaterialInfo_h
 
#include "../CbmL1Def.h"
#include <vector>
 
class L1MaterialInfo {
public:
  L1MaterialInfo() : thick(0), RL(0), RadThick(0), logRadThick(0) {};
 
  fvec thick, RL, RadThick, logRadThick;
} _fvecalignment;
 
class L1Material {
public:
  L1Material() : table(0), NBins(0), RMax(0.), iD(0.) {};
 
  std::vector<std::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