TVacuumPol.h

Go to the documentation of this file.

#ifndef RAD_TVacuumPol
#define RAD_TVacuumPol
 
#include <string>
#include <complex>
#include "TRadSpline.h"
 
class TVacuumPol{
protected:
  TRadSpline *fsvtr;
  TRadSpline *fsvsr;
  TRadSpline *fsvsi;
  bool fNoVP;
  bool fReadSuccessfullyVP;
  void ReadVacuumPolData(std::string);
public:
  TVacuumPol();
  ~TVacuumPol();
  void Init(std::string, std::string);
  void SetZeroVP(){fNoVP = true;}
  void SetDefaultVP(){fNoVP = !fReadSuccessfullyVP;}
  bool GetVacuumPolarization(){return fNoVP;}
  inline std::complex <double> Evals(const double &s){
    if(fNoVP) return std::complex<double>(1,0);
    return std::complex<double>(fsvsr->Eval(s),fsvsi->Eval(s));
  }
  inline std::complex <double> Evalt(const double &t){
    if(fNoVP) return std::complex<double>(1,0);
    return std::complex<double>(fsvtr->Eval(fabs(t)),0.);
  }
};
 
extern TVacuumPol *gVPol;
 
#endif // #ifndef RAD_TVacuumPol