JACoeffSource.hh

Go to the documentation of this file.

#ifndef __JPHYSICS__JACOEFFSOURCE__
#define __JPHYSICS__JACOEFFSOURCE__
 
#include <memory>
 
#include "JPhysics/JRadiation.hh"
 
 
/**
 * \author mdejong
 */
 
namespace JPHYSICS {}
namespace JPP { using namespace JPHYSICS; }
 
namespace JPHYSICS {
  
 
  /**
   * Interface for calculation of ionization constant.
   */
  class JACoeffInterface {
  public:
    /**
     * Virtual destructor.
     */
    virtual ~JACoeffInterface()
    {}
 
    /**
     * Get ionization constant a.
     *
     * \param  E      muon energy [GeV]
     * \return        ionization constant [GeV/m]
     */
    virtual double getA(const double E) const = 0;
  };
  
  
  /**
   * Implementation for calculation of ionization constant.
 
   * This class implements the JACoeffInterface interface.\n
   * N.B: This class owns the object pointed to using JLANG::JSharedPointer.
   */
  class JACoeffSource : 
    public JACoeffInterface,
    public std::shared_ptr<JRadiation>
  {
  public:
    /**
     * Constructor.
     *
     * \param  radiation       pointer to valid JRadition object
     * \param  density         mass density of radiation material [gr/cm³]
     */
    JACoeffSource(const std::shared_ptr<JRadiation>& radiation,
                  const double                       density) :
      std::shared_ptr<JRadiation>(radiation),
      rho(density)
    {}
 
    
    /**
     * Get ionization constant a.
     *
     * \param  E      muon energy [GeV]
     * \return        ionization constant [GeV/m]
     */
    virtual double getA(const double E) const override 
    {
      return get()->CalculateACoeff(E) * rho * 1.0e6; 
    }
 
  protected:
    const double rho;
  };
}    
 
#endif