JRadiationSource.hh

Go to the documentation of this file.

#ifndef __JPHYSICS__JRADIATIONSOURCE__
#define __JPHYSICS__JRADIATIONSOURCE__
 
#include <memory>
 
#include "TRandom3.h"
 
#include "JLang/JObjectID.hh"
 
#include "JPhysics/JConstants.hh"
#include "JPhysics/JRadiation.hh"
#include "JPhysics/JDIS.hh"
#include "JPhysics/JDeltaRays.hh"
 
 
/**
 * \author mdejong
 */
 
namespace JPHYSICS {}
namespace JPP { using namespace JPHYSICS; }
 
namespace JPHYSICS {
 
  using JLANG::JObjectID;
 
  /**
   * Interface for calculation of inverse interaction length and shower energy.
   */
  class JRadiationInterface :
    public JObjectID
  {
  public:
    /**
     * Constructor.
     *
     * \param  id     radiation identifier
     */
    JRadiationInterface(const int id) :
      JObjectID(id)
    {}
 
 
    /**
     * Virtual destructor.
     */
    virtual ~JRadiationInterface()
    {}
 
 
    /**
     * Get inverse interaction length.
     *
     * \param  E      muon energy [GeV]
     * \return        inverse interaction length [m^-1]
     */
    virtual double getInverseInteractionLength(const double E) const = 0;
 
 
    /**
     * Get energy of shower.
     *
     * \param  E      muon   energy [GeV]
     * \return        shower energy [GeV]
     */
    virtual double getEnergyOfShower(const double E) const = 0;
 
 
    /**
     * Get RMS of scattering angle.
     *
     * \param  E      muon   energy [GeV]
     * \param  Es     shower energy [GeV]
     * \return        RMS scattering angle [rad]
     */
    virtual double getThetaRMS(const double E, const double Es) const = 0;
  };
 
 
  /**
   * Implementation for calculation of inverse interaction length and shower energy.
   *
   * This class implements the JRadiationInterface interface.
   */
  class JRadiationSource :
    public JRadiationInterface,
    public std::shared_ptr<JRadiation>,
    public JRadiation::radiation_type
  {
  public:
    /**
     * Constructor.
     *
     * \param  id              radiation identifier
     * \param  radiation       radiation
     * \param  density         mass density of radiation material [gr/cm³]
     * \param  source          radiation type
     */
    JRadiationSource(const int                          id,
                     const std::shared_ptr<JRadiation>& radiation,
                     const double                       density,
                     const JRadiation::radiation_type   source) :
      JRadiationInterface(id),
      std::shared_ptr<JRadiation>(radiation),
      radiation_type(source),
      rho(density)
    {}
 
 
    /**
     * Get inverse interaction length.
     *
     * \param  E      muon energy [GeV]
     * \return        inverse interaction length [m^-1]
     */
    virtual double getInverseInteractionLength(const double E) const override 
    { 
      return (get()->*sigma)(E) * rho * 1.0e6;
    }
 
 
    /**
     * Get energy of shower.
     *
     * \param  E      muon   energy [GeV]
     * \return        shower energy [GeV]
     */
    virtual double getEnergyOfShower(const double E) const override 
    { 
      return (get()->*eloss)(E);
    }
 
 
    /**
     * Get RMS of scattering angle.
     *
     * \param  E      muon   energy [GeV]
     * \param  Es     shower energy [GeV]
     * \return        RMS scattering angle [rad]
     */
    virtual double getThetaRMS(const double E, const double Es) const override
    {
      return (get()->*theta)(E, Es/E);
    }
 
  protected:
    const double rho;
  };
 
 
  /**
   * Implementation for calculation of inverse interaction length and shower energy due to deep-inelastic muon-nucleon scattering.
   *
   * This class implements the JRadiationInterface interface.
   */
  class JDISSource :
    public JRadiationInterface,
    public JDIS
  {
  public:
    /**
     * Constructor.
     *
     * \param  id              radiation identifier
     * \param  density         mass density of radiation material [gr/cm³]
     */
    JDISSource(const int    id,
               const double density) :
      JRadiationInterface(id),
      rho(density)
    {}
    
    
    /**
     * Get inverse interaction length.
     *
     * \param  E      muon energy [GeV]
     * \return        inverse interaction length [m^-1]
     */
    virtual double getInverseInteractionLength(const double E) const override 
    {
      return this->getCrossSection(E) * (rho / NUCLEON_MOLAR_MASS) * AVOGADRO * 1.0e2;
    }
 
 
    /**
     * Get energy of shower.
     *
     * \param  E      muon   energy [GeV]
     * \return        shower energy [GeV]
     */
    virtual double getEnergyOfShower(const double E) const override 
    { 
      return this->getE(E);
    }
 
 
    /**
     * Get RMS of scattering angle.
     *
     * \param  E      muon   energy [GeV]
     * \param  Es     shower energy [GeV]
     * \return        RMS scattering angle [rad]
     */
    virtual double getThetaRMS(const double E, const double Es) const override
    {
      return static_cast<const JDIS*>(this)->getThetaRMS(E, Es/E);
    }
    
  protected:
    const double rho;
  };
 
 
  /**
   * Implementation for calculation of inverse interaction length and shower energy due to delta-rays from muon.
   *
   * This class implements the JRadiationInterface interface.
   */
  class JDeltaRaysSource :
    public JRadiationInterface
  {
  public:
    /**
     * Constructor.
     *
     * \param  id              radiation identifier
     * \param  density         mass density of radiation material [gr/cm³]
     * \param  A               atomic mass
     * \param  Z               atomic number
     * \param  Tmin            minimal kunetic energy [GeV]
     */
    JDeltaRaysSource(const int    id,
                     const double density,
                     const double Tmin) :
      JRadiationInterface(id),
      rho(density),
      Tmin(Tmin)
    {}
 
 
    /**
     * Get inverse interaction length.
     *
     * \param  E      muon energy [GeV]
     * \return        inverse interaction length [m^-1]
     */
    virtual double getInverseInteractionLength(const double E) const override 
    {
      const double Tmax = JDeltaRays::getTmax(E, MASS_MUON);
 
      return JDeltaRays::getCount(E, MASS_MUON, Tmin, Tmax) * rho * 1.0e2;
    }
 
 
    /**
     * Get energy of shower.
     *
     * \param  E      muon   energy [GeV]
     * \return        shower energy [GeV]
     */
    virtual double getEnergyOfShower(const double E) const override 
    { 
      const double Tmax = JDeltaRays::getTmax(E, MASS_MUON);
 
      const double xmin = 1.0 / Tmax;
      const double xmax = 1.0 / Tmin;
 
      return 1.0 / gRandom->Uniform(xmin, xmax);
    }
 
 
    /**
     * Get RMS of scattering angle.
     *
     * \param  E      muon   energy [GeV]
     * \param  Es     shower energy [GeV]
     * \return        RMS scattering angle [rad]
     */
    virtual double getThetaRMS(const double E, const double Es) const override
    {
      return 0.0;
    }
 
  protected:
    const double  rho;
    const double  Tmin;    //!< minimal kinetic energy [GeV]
  };
}    
 
#endif