TPCGas.h

#ifndef TPCGAS_H
#define TPCGAS_H 1

#include <string>

namespace marlintpc
{
  
  class TPCGas;
  
  class TPCGas
  {
  public:
    TPCGas(std::string gasname);
    ~TPCGas();
    void SetEField(double EField){_EField=EField;}
    double GetEField(){return _EField;}
    void SetBField(double BField){_BField=BField;}
    double GetBField(){return _BField;}
    double GetDriftVelocity(){return _vdrift;}
    double GetTransverseDiffusion(){return _diffTrans;}
    double GetLongitudinalDiffusion(){return _diffLong;}
    double GetAttachment(double EField, double dist);
    void CalculateGasParameters();
    std::string GetGasName(){return _gasName;}
 
  protected:
    void CalculateTransverseDiffusion();
    double EvaluateEFieldDependence(double a, double b ,double c, double d,
                                    double e, double f, double g, double h,
                                    double i, double j, double k);
    double EvaluateBFieldDependence(double a, double b, double c, double d,double e, double f);
    
    void initTDRGas();
    void initP5Gas();
    void initP10Gas();
    void initT2KGas();

    double _BField;
    double _EField;
    double _vdrift;
    double _diffLong;
    double _diffTrans;
    double _attachment;

    double _dL_params[11];
    double _vd_params[11];
    double _dT_params[10][6];
    double _att[7];

    std::string _gasName;

   };
  
} //namespace marlintpc 

#endif // TPCGAS_H