ConstantFluorescenceModel.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 <atm/ConstantFluorescenceModel.h>
#include <fwk/CentralConfig.h>
#include <sstream>

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


void
ConstantFluorescenceModel::Init()
{
  Branch topB = 
    CentralConfig::GetInstance()->GetTopBranch("ConstantFluorescenceModel");
  
  topB.GetChild("wavelength").GetData(fWavelength);
  topB.GetChild("fluorescenceYield").GetData(fFluorescenceYield);
  topB.GetChild("totalLight").GetData(fTotalLight);

  fSpectrumNormalization = 0.0;
  for (unsigned int i = 0; i < fFluorescenceYield.size(); ++i)
    fSpectrumNormalization += fFluorescenceYield[i];
}


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


const utl::TabulatedFunction&
ConstantFluorescenceModel::EvaluateFluorescenceYield(const double /*heightAboveSeaLevel*/)
  const
{ 
  // Normalizing to the total light
  vector<double> normalizedFY = fFluorescenceYield;
  double total = 0;
  
  for (int i = 0; i < int(fFluorescenceYield.size()); ++i)
    total += fFluorescenceYield[i];

  fFluorescenceSpectrum.Clear();
  
  for (int i = 0; i < int(normalizedFY.size()); ++i) {
    normalizedFY[i] *= fTotalLight/total;
    fFluorescenceSpectrum.PushBack(fWavelength[i], normalizedFY[i]);
  }

  return fFluorescenceSpectrum;
}


double
ConstantFluorescenceModel::GetdEdX0()
  const
{
  // Energy loss evaluated for an electron of 1.4 MeV at standard 
  // temperature and pressure.
  // Value taken from the ESTAR program of NIST (2004)
  // http://physics.nist.gov/Star 
  return 1.676*MeV/(g/cm2);
}


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