Nagano2004FluorescenceModel.cc

Go to the documentation of this file.

#include <utl/TabulatedFunction.h>
#include <utl/Reader.h>
#include <utl/ErrorLogger.h>
#include <utl/MathConstants.h>
#include <utl/Vector.h>
#include <utl/PhysicalConstants.h>
#include <utl/Transformation.h>
#include <utl/TabulatedFunctionErrors.h>
#include <det/Detector.h>
#include <atm/Atmosphere.h>
#include <atm/ProfileResult.h>
#include <atm/Nagano2004FluorescenceModel.h>
#include <fwk/CentralConfig.h>
#include <sstream>

using namespace utl;
using namespace atm;
using namespace fwk;
using namespace std;

  

void Nagano2004FluorescenceModel::Init(){
  CentralConfig* cc = CentralConfig::GetInstance(); 
  
  Branch topB = 
    cc->GetTopBranch("Nagano2004FluorescenceModel");
  
  topB.GetChild("wavelength").GetData(fWavelength);
  topB.GetChild("naganoA").GetData(fNaganoA);
  topB.GetChild("naganoAError").GetData(fNaganoAError);
  topB.GetChild("naganoB").GetData(fNaganoB);
  topB.GetChild("naganoBError").GetData(fNaganoBError);
}


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

}


const utl::TabulatedFunction&
Nagano2004FluorescenceModel::EvaluateFluorescenceYield
(const double heightAboveSeaLevel) const{

  double fluorescenceYield = 0;

  const Atmosphere& atmo = det::Detector::GetInstance().GetAtmosphere();
  const ProfileResult& tempProfile    = atmo.EvaluateTemperatureVsHeight();
  const ProfileResult& densityProfile = atmo.EvaluateDensityVsHeight();
  const double temperature = tempProfile.Y(heightAboveSeaLevel);
  const double density     = densityProfile.Y(heightAboveSeaLevel);

  

  fFluorescenceSpectrum.Clear();
  for (unsigned int iwl=0; iwl<fWavelength.size(); iwl++){
    
    double wavelength = fWavelength[iwl];
    
    fluorescenceYield = fNaganoA[iwl]*density/(1.+ density*fNaganoB[iwl] 
                                               *sqrt(temperature));
    
    fFluorescenceSpectrum.PushBack(wavelength, fluorescenceYield);
    
  }
  
  return fFluorescenceSpectrum;

}

double Nagano2004FluorescenceModel::GetdEdX0() const {

    // Energy loss evaluated for an electron of .85 MeV at standard 
    // temperature and pressure.
    // Value taken from the ESTAR program of NIST (2004)
    // http://physics.nist.gov/PhysRefData/Star/Text/ESTAR.html 

    return 1.675*MeV/(g/cm/cm);
}

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k"
// End: