VFluorescenceModel.cc

Go to the documentation of this file.

#include <cmath>
#include <atm/VFluorescenceModel.h>
#include <utl/AugerUnits.h>
using namespace utl;
using namespace atm;

double 
VFluorescenceModel::GetDeExcitationTime(const double height_asl)
  const 
{

  // see GAP-2007-099 Eq. (6.14)

  const double tau_0 = 37.5*ns; // vacuum lifetime
  const double H     = 8005.*m; // (scale height)
  const double a     = 95.;
  
  const double lifetime = (
                           height_asl>=0. ?
                           tau_0 / (a * exp(-height_asl/H) + 1.) :
                           0.
                           );

  return lifetime;

}


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