TRT_PAI_effectiveGas.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRT_PAI_effectiveGas_h
#define TRT_PAI_effectiveGas_h
 
#include "AthenaBaseComps/AthMessaging.h"
 
#include <vector>
 
class TRT_PAI_gasMixture;

class TRT_PAI_effectiveGas : public AthMessaging {
public:
  TRT_PAI_effectiveGas(TRT_PAI_gasMixture * gm,
               double Emin,
               double Emax,
               double tempK,
               double eps);
  void GasTab(const std::vector<float> & gamvec,
          std::vector<float>& EArray,
          std::vector< std::vector<float> >& fnArray,
          std::vector<float>& dndx);
  
private:
  double XSigma(double lnE, double dummy);
  double XFReal(double lnD, double lnE);
  double XGInt(double (TRT_PAI_effectiveGas::*pt2Func)(double,double),
           double  lnLo,
           double  lnHi,
           const double  eps,
           double  extraParameter);
  double dndedx(double lgE, double gamma);
  std::vector<float> m_lnELvls;
  std::vector<float> m_lnFosc;
  std::vector<float> m_lnIntegratedSigmas;
  std::vector<float> m_lnEpsR;
  std::vector<float> m_lnEpsI;
  const double m_lnEmin;
  const double m_lnEmax;
  const double m_eps;        // epsilon for numerical integration.
  double m_S1 ;              // x-section to F.osc
  double m_Wp2;              // plasma freq**2 {ev}
  double m_S2 ;              // dN/dx scale
  double m_ne ;              // Electron density
};
 
#endif