LitAbsorber.h

Go to the documentation of this file.

 
#ifndef LITABSORBER_H_
#define LITABSORBER_H_
 
#include "../LitFieldGrid.h"
 
namespace lit {
  namespace parallel {
 
    template<class T>
    class LitAbsorber {
    public:
      LitAbsorber()
        : fZ(0.)
        , fMaterial()
        , fFieldGridFront()
        , fFieldGridMiddle()
        , fFieldGridBack() {}
 
      virtual ~LitAbsorber() {}
 
      const T& GetZ() const { return fZ; }
 
      void SetZ(const T& z) { fZ = z; }
 
      const LitMaterialInfo<T>& GetMaterial() const { return fMaterial; }
 
      void SetMaterial(const LitMaterialInfo<T>& material) {
        fMaterial = material;
      }
 
      const LitFieldGrid& GetFieldGridFront() const { return fFieldGridFront; }
 
      void SetFieldGridFront(const LitFieldGrid& fieldGrid) {
        fFieldGridFront = fieldGrid;
      }
 
      const LitFieldGrid& GetFieldGridMiddle() const {
        return fFieldGridMiddle;
      }
 
      void SetFieldGridMiddle(const LitFieldGrid& fieldGrid) {
        fFieldGridMiddle = fieldGrid;
      }
 
      const LitFieldGrid& GetFieldGridBack() const { return fFieldGridBack; }
 
      void SetFieldGridBack(const LitFieldGrid& fieldGrid) {
        fFieldGridBack = fieldGrid;
      }
 
      std::string ToString() const {
        return "LitAbsorber: Z=" + lit::parallel::ToString<T>(GetZ())
               + ", material=" + GetMaterial().ToString();
      }
 
      friend std::ostream& operator<<(std::ostream& strm,
                                      const LitAbsorber& absorber) {
        strm << absorber.ToString();
        return strm;
      }
 
    private:
      T fZ;                          // Z position of the absorber in [cm]
      LitMaterialInfo<T> fMaterial;  // Absorber material
      LitFieldGrid
        fFieldGridFront;  // Magnetic field approximation in front of absorber
      LitFieldGrid
        fFieldGridMiddle;  // Magnetic field approximation in the middle of absorber
      LitFieldGrid
        fFieldGridBack;  // Magnetic field approximation in the back of absorber
    } _fvecalignment;
 
    /* Some typedefs for convenience */
    typedef LitAbsorber<fvec> LitAbsorberVec;
    typedef LitAbsorber<fscal> LitAbsorberScal;
 
  }  // namespace parallel
}  // namespace lit
 
#endif /* LITABSORBER_H_ */