TRT_PAI_effectiveGas.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2025 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