TVacuumPol.h

Go to the documentation of this file.

#ifndef RAD_TVacuumPol
#define RAD_TVacuumPol
 
#include "TRadSpline.h"
#include <complex>
#include <string>
 
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