AirFluorescenceModel.h

Go to the documentation of this file.

#ifndef _atm_AirFluorescenceModel_h_
#define _atm_AirFluorescenceModel_h_

#include <atm/VFluorescenceModel.h>
#include <utl/TabulatedFunction.h>


namespace utl {
  class TabulatedFunction;
  class Reader;
}

namespace atm {

  class AirFluorescenceModel : public atm::VFluorescenceModel {

  public:
    AirFluorescenceModel() = default;

    virtual ~AirFluorescenceModel() = default;

    void Init();

    const utl::TabulatedFunction&
    EvaluateFluorescenceYield(const double heightAboveSeaLevel) const;

    const std::vector<double>& GetWavelengths() const { return fWavelength; }

    double GetdEdX0() const;

    bool HasData() const { return true; }

  private:
    enum FluorescenceDataSet {
      eAuger,
      eKakimoto,
      eNagano,
      eKeilhauer,
      eUserDefined
    };

    enum TemperatureParameterization {
      eAIRtemp,
      eNoTempParam
    };

    enum TemperatureParameterizationWater {
      eAIRWatertemp,
      eNoWaterTempParam
    };

    enum HumidityParameterization {
      eAIRhum,
      eNoHumParam
    };

    TemperatureParameterization fTempParam;
    TemperatureParameterizationWater fTempWaterParam;
    HumidityParameterization fHumParam;

    mutable utl::TabulatedFunction fFluorescenceSpectrum;

    std::vector<double> fWavelength;
    std::vector<double> fRelativeIntensity;
    std::vector<double> fPPrime;
    std::vector<double> fRelativeIntensityError;
    std::vector<double> fPPrimeErrorUncorrelated;
    std::vector<double> fPPrimeErrorCorrelated;
    std::vector<double> fAlpha;
    std::vector<double> fAlphaWater;
    std::vector<double> fPPrimeHum;

    std::string NaganoAirflyFileName;
    std::string NaganoFileName;

    double fTemperature0 = 0;
    double fPressure0 = 0;
    double fYield337 = 0;

    REGISTER_MODEL(atm::VFluorescenceModel, "AirFluorescence", AirFluorescenceModel);

  };

}


#endif