JPHYSICS::JGeaneWater Class Reference

Function object for energy dependent energy loss of the muon. More...

#include <JGeane.hh>

Inheritance diagram for JPHYSICS::JGeaneWater:

Public Member Functions
	JGeaneWater ()
	Default constructor.
virtual double	getA () const override
	Get energy loss constant.
virtual double	getB () const override
	Get energy loss constant.
virtual double	getE (const double E, const double dx) const override
	Get energy of muon after specified distance.
double	getEa (const double E, const double dx) const
	Get energy loss due to ionisation.
double	getEb (const double E, const double dx) const
	Get energy loss due to pair production and bremsstrahlung.
virtual double	getX (const double E0, const double E1) const override
	Get distance traveled by muon.
double	operator() (const double E, const double dx) const
	Energy of muon after specified distance.
double	operator() (const double E) const
	Range of muon.
double	operator() () const
	Equivalent unit track length per unit shower energy and per unit track length.

Detailed Description

Function object for energy dependent energy loss of the muon.

Approximate values of energy loss parameters taken from reference: Proceedings of ICRC 2001, "Precise parametrizations of muon energy losses in water", S. Klimushin, E. Bugaev and I. Sokalski.

Definition at line 206 of file JGeane.hh.

Constructor & Destructor Documentation

JGeaneWater()

JPHYSICS::JGeaneWater::JGeaneWater ( )

inline

Default constructor.

Definition at line 214 of file JGeane.hh.

    {
      using namespace std;
 
      this->insert(make_pair( 0.0e0,  JGeane_t( 2.30e-1 * DENSITY_SEA_WATER, 15.50e-4 * DENSITY_SEA_WATER)));
      this->insert(make_pair(30.0e0,  JGeane_t( 2.67e-1 * DENSITY_SEA_WATER,  3.40e-4 * DENSITY_SEA_WATER)));
      this->insert(make_pair(35.3e3,  JGeane_t(-6.50e-1 * DENSITY_SEA_WATER,  3.66e-4 * DENSITY_SEA_WATER)));
    }

Member Function Documentation

getA()

virtual double JPHYSICS::JGeaneWater::getA ( ) const

inlineoverridevirtual

Get energy loss constant.

N.B. The return value corresponds to the low-energy regime.

Returns

Energy loss due to ionisation [GeV/m]

Implements JPHYSICS::JGeane.

Definition at line 231 of file JGeane.hh.

    {
      return 2.30e-1 * DENSITY_SEA_WATER; //This is the value for low energy (<30 GeV), the value used for step(ds) and getRange())
    }

getB()

virtual double JPHYSICS::JGeaneWater::getB ( ) const

inlineoverridevirtual

Get energy loss constant.

N.B. The return value corresponds to the medium-energy regime.

Returns

Energy loss due to pair production and bremsstrahlung [m^-1]

Implements JPHYSICS::JGeane.

Definition at line 244 of file JGeane.hh.

    {
      return 3.40e-4 * DENSITY_SEA_WATER;
    }

getE()

virtual double JPHYSICS::JGeaneWater::getE ( const double E, const double dx ) const

inlineoverridevirtual

Get energy of muon after specified distance.

Parameters

E	Energy of muon [GeV]
dx	distance traveled [m]

Returns

Energy of muon [GeV]

Implements JPHYSICS::JGeane.

Definition at line 257 of file JGeane.hh.

    {
      double E1 = E;
      double x1 = dx;
 
      if (E1 > MASS_MUON / getSinThetaC()) {
 
        const_iterator p = this->lower_bound(E1);
 
        do {
 
          --p;
 
          const double x2 = p->second.getX(E1, p->first);
 
          if (x2 > x1) {
            return p->second.getE(E1, x1);
          }
 
          E1  = p->first;
          x1 -= x2;
 
        } while (p != this->begin());
      }
 
      return E1;
    }

getEa()

double JPHYSICS::JGeaneWater::getEa ( const double E, const double dx ) const

inline

Get energy loss due to ionisation.

Parameters

E	initial energy [GeV]
dx	distance traveled [m]

Returns

energy loss due to ionisation [GeV]

Definition at line 295 of file JGeane.hh.

    {
      using namespace std;
      using namespace JPP;
      
      double Ea = 0.0;
      
      double E1 = E;
      double x1 = dx;
 
      if (E1 > MASS_MUON / getSinThetaC()) {
 
        map<double, JGeane_t>::const_iterator p = this->lower_bound(E1);
 
        do {
 
          --p;
 
          const double x2 = p->second.getX(E1, p->first);
 
          Ea += (x2 > x1 ? x1 : x2) * p->second.getA();
          E1  = p->first;
          
          x1 -= x2;
          
        } while (p != this->cbegin() && x1 > 0.0);
      }
 
      return Ea;
    }

getEb()

double JPHYSICS::JGeaneWater::getEb ( const double E, const double dx ) const

inline

Get energy loss due to pair production and bremsstrahlung.

Parameters

E	initial energy [GeV]
dx	distance traveled [m]

Returns

energy loss due to pair production and bremsstrahlung [GeV]

Definition at line 334 of file JGeane.hh.

    {
      const double dE = E - getE(E, dx);
      
      return dE - getEa(E, dx);
    }

getX()

virtual double JPHYSICS::JGeaneWater::getX ( const double E0, const double E1 ) const

inlineoverridevirtual

Get distance traveled by muon.

Parameters

E0	Energy of muon at start [GeV]
E1	Energy of muon at end [GeV]

Returns

distance traveled [m]

Implements JPHYSICS::JGeane.

Definition at line 349 of file JGeane.hh.

    {
      double E  = E0;
      double dx = 0.0;
 
      if (E > MASS_MUON / getSinThetaC()) {
 
        const_iterator p = this->lower_bound(E);
 
        do {
 
          --p;
 
          if (E1 > p->first) {
            return dx + p->second.getX(E, E1);
          }
 
          dx += p->second.getX(E, p->first);
          E   = p->first;
 
        } while (p != this->begin());
      }
 
      return dx;
    }

operator()() [1/3]

double JPHYSICS::JGeane::operator() ( const double E, const double dx ) const

inlineinherited

Energy of muon after specified distance.

Parameters

E	Energy of muon [GeV]
dx	distance traveled [m]

Returns

Energy of muon [GeV]

Definition at line 85 of file JGeane.hh.

    {
      return this->getE(E, dx);
    }

operator()() [2/3]

double JPHYSICS::JGeane::operator() ( const double E ) const

inlineinherited

Range of muon.

Parameters

E	Energy of muon [GeV]

Returns

range [m]

Definition at line 97 of file JGeane.hh.

    {
      return this->getX(E, 0.0);
    }

operator()() [3/3]

double JPHYSICS::JGeane::operator() ( ) const

inlineinherited

Equivalent unit track length per unit shower energy and per unit track length.

Returns

equivalent unit track length [Gev^-1]

Definition at line 108 of file JGeane.hh.

    {
      return this->getB() * geanc();
    }

---

The documentation for this class was generated from the following file:

• JGeane.hh