genie::mueloss::PetrukhinShestakovModel Class Reference

Bethe-Heitler, Petrukhin-Shestakov model for the energy loss of muons due to bremsstrahlung. Concrete implementation of the MuELossI interface. More...

#include <PetrukhinShestakovModel.h>

Inheritance diagram for genie::mueloss::PetrukhinShestakovModel:

Collaboration diagram for genie::mueloss::PetrukhinShestakovModel:

Public Member Functions
	PetrukhinShestakovModel ()
	PetrukhinShestakovModel (string config)
virtual	~PetrukhinShestakovModel ()
double	dE_dx (double E, MuELMaterial_t material) const
	Implement the MuELossI interface. More...
MuELProcess_t	Process (void) const
Public Member Functions inherited from genie::mueloss::MuELossI
virtual	~MuELossI ()
Public Member Functions inherited from genie::Algorithm
virtual	~Algorithm ()
virtual void	Configure (const Registry &config)
virtual void	Configure (string config)
virtual void	FindConfig (void)
virtual const Registry &	GetConfig (void) const
Registry *	GetOwnedConfig (void)
virtual const AlgId &	Id (void) const
	Get algorithm ID. More...
virtual AlgStatus_t	GetStatus (void) const
	Get algorithm status. More...
virtual bool	AllowReconfig (void) const
virtual AlgCmp_t	Compare (const Algorithm *alg) const
	Compare with input algorithm. More...
virtual void	SetId (const AlgId &id)
	Set algorithm ID. More...
virtual void	SetId (string name, string config)
const Algorithm *	SubAlg (const RgKey &registry_key) const
void	AdoptConfig (void)
void	AdoptSubstructure (void)
virtual void	Print (ostream &stream) const
	Print algorithm info. More...

Additional Inherited Members
Static Public Member Functions inherited from genie::Algorithm
static string	BuildParamVectKey (const std::string &comm_name, unsigned int i)
static string	BuildParamVectSizeKey (const std::string &comm_name)
static string	BuildParamMatKey (const std::string &comm_name, unsigned int i, unsigned int j)
static string	BuildParamMatRowSizeKey (const std::string &comm_name)
static string	BuildParamMatColSizeKey (const std::string &comm_name)
Protected Member Functions inherited from genie::mueloss::MuELossI
	MuELossI ()
	MuELossI (string name)
	MuELossI (string name, string config)
Protected Member Functions inherited from genie::Algorithm
	Algorithm ()
	Algorithm (string name)
	Algorithm (string name, string config)
void	Initialize (void)
void	DeleteConfig (void)
void	DeleteSubstructure (void)
Registry *	ExtractLocalConfig (const Registry &in) const
Registry *	ExtractLowerConfig (const Registry &in, const string &alg_key) const
	Split an incoming configuration Registry into a block valid for the sub-algo identified by alg_key. More...
template<class T >
bool	GetParam (const RgKey &name, T &p, bool is_top_call=true) const
template<class T >
bool	GetParamDef (const RgKey &name, T &p, const T &def) const
template<class T >
int	GetParamVect (const std::string &comm_name, std::vector< T > &v, bool is_top_call=true) const
	Handle to load vectors of parameters. More...
int	GetParamVectKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
template<class T >
int	GetParamMat (const std::string &comm_name, TMatrixT< T > &mat, bool is_top_call=true) const
	Handle to load matrix of parameters. More...
template<class T >
int	GetParamMatSym (const std::string &comm_name, TMatrixTSym< T > &mat, bool is_top_call=true) const
int	GetParamMatKeys (const std::string &comm_name, std::vector< RgKey > &k, bool is_top_call=true) const
int	AddTopRegistry (Registry *rp, bool owns=true)
	add registry with top priority, also update ownership More...
int	AddLowRegistry (Registry *rp, bool owns=true)
	add registry with lowest priority, also update ownership More...
int	MergeTopRegistry (const Registry &r)
int	AddTopRegisties (const vector< Registry * > &rs, bool owns=false)
	Add registries with top priority, also udated Ownerships. More...
Protected Attributes inherited from genie::Algorithm
bool	fAllowReconfig
bool	fOwnsSubstruc
	true if it owns its substructure (sub-algs,...) More...
AlgId	fID
	algorithm name and configuration set More...
vector< Registry * >	fConfVect
vector< bool >	fOwnerships
	ownership for every registry in fConfVect More...
AlgStatus_t	fStatus
	algorithm execution status More...
AlgMap *	fOwnedSubAlgMp
	local pool for owned sub-algs (taken out of the factory pool) More...

Detailed Description

Bethe-Heitler, Petrukhin-Shestakov model for the energy loss of muons due to bremsstrahlung. Concrete implementation of the MuELossI interface.

References:

W.Lohmann, R.Kopp and R.Voss, Energy Loss of Muons in the Energy Range 1-10000 GeV, CERN 85-03

Author

Costas Andreopoulos <c.andreopoulos cern.ch> University of Liverpool

Created:

December 10, 2003

License:

Copyright (c) 2003-2024, The GENIE Collaboration For the full text of the license visit http://copyright.genie-mc.org

Definition at line 39 of file PetrukhinShestakovModel.h.

Constructor & Destructor Documentation

PetrukhinShestakovModel::PetrukhinShestakovModel

(

)

Definition at line 28 of file PetrukhinShestakovModel.cxx.

                                                 :
MuELossI("genie::mueloss::PetrukhinShestakovModel")
{

}

PetrukhinShestakovModel::PetrukhinShestakovModel

(

string

config

)

Definition at line 34 of file PetrukhinShestakovModel.cxx.

                                                              :
MuELossI("genie::mueloss::PetrukhinShestakovModel", config)
{

}

PetrukhinShestakovModel::~PetrukhinShestakovModel ( )

virtual

Definition at line 40 of file PetrukhinShestakovModel.cxx.

{

}

Member Function Documentation

double PetrukhinShestakovModel::dE_dx ( double  E, MuELMaterial_t  material  ) const

virtual

Implement the MuELossI interface.

Implements genie::mueloss::MuELossI.

Definition at line 45 of file PetrukhinShestakovModel.cxx.

References genie::units::A, genie::mueloss::MuELMaterial::A(), genie::mueloss::eMuUndefined, genie::units::g, genie::utils::gsl::Integration1DimTypeFromString(), genie::mueloss::kMaxMuE, genie::constants::kMuonMass, genie::constants::kNA, genie::constants::kSqrtNapierConst, Process(), genie::mueloss::MuELProcess::Threshold(), and genie::mueloss::MuELMaterial::Z().

{
// Calculate the muon -dE/dx due to muon bremsstrahlung (in GeV^-2).
// To convert the result to more handly units, eg MeV/(gr/cm^2), just write:
// dE_dx /= (units::MeV/(units::g/units::cm2));

  if(material == eMuUndefined) return 0;
  if(E<=MuELProcess::Threshold(this->Process()) || E>=kMaxMuE) return 0;

  // material Z,E
  double Z = MuELMaterial::Z(material);
  double A = MuELMaterial::A(material);

  // calculate (the min,max) fraction of energy, v,  carried to the photon
  double Vmin = 0.;
  double Vmax = 1. - 0.75*kSqrtNapierConst* (kMuonMass/E) * TMath::Power(Z,1/3.);

  // integrate the Bethe-Heitler muon bremsstrahlung 
  // cross section v*ds/dv over v

  ROOT::Math::IBaseFunctionOneDim * integrand =
          new gsl::PetrukhinShestakovIntegrand(E,Z);
  ROOT::Math::IntegrationOneDim::Type ig_type =
          utils::gsl::Integration1DimTypeFromString("adaptive");

  double abstol   = 1;    // We mostly care about relative tolerance   
  double reltol   = 1E-4;
  int    nmaxeval = 100000;
  ROOT::Math::Integrator ig(*integrand,ig_type,abstol,reltol,nmaxeval);

  // calculate the b factor (bE = -dE/dx) in GeV^-3
  A *= units::g;
  double bbrem = (kNA/A) * ig.Integral(Vmin, Vmax);

  delete integrand;

  // calculate the dE/dx due to muon bremsstrahlung in GeV^-2
  double de_dx = bbrem*E;
  return de_dx;
}

MuELProcess_t genie::mueloss::PetrukhinShestakovModel::Process ( void  ) const

inlinevirtual

Implements genie::mueloss::MuELossI.

Definition at line 48 of file PetrukhinShestakovModel.h.

References genie::mueloss::eMupBremsstrahlung.

Referenced by dE_dx().

{ return eMupBremsstrahlung; }

---

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

• PetrukhinShestakovModel.h
• PetrukhinShestakovModel.cxx