GENIEGenerator
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
Public Member Functions | List of all members
genie::mueloss::KokoulinPetrukhinModel Class Reference

Kokoulin-Petrukhin model for the energy loss of muons due to direct e+e- pair production. Concrete implementation of the MuELossI interface. More...

#include <KokoulinPetrukhinModel.h>

Inheritance diagram for genie::mueloss::KokoulinPetrukhinModel:
Inheritance graph
[legend]
Collaboration diagram for genie::mueloss::KokoulinPetrukhinModel:
Collaboration graph
[legend]

Public Member Functions

 KokoulinPetrukhinModel ()
 
 KokoulinPetrukhinModel (string config)
 
virtual ~KokoulinPetrukhinModel ()
 
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 RegistryGetConfig (void) const
 
RegistryGetOwnedConfig (void)
 
virtual const AlgIdId (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 AlgorithmSubAlg (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)
 
RegistryExtractLocalConfig (const Registry &in) const
 
RegistryExtractLowerConfig (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...
 
AlgMapfOwnedSubAlgMp
 local pool for owned sub-algs (taken out of the factory pool) More...
 

Detailed Description

Kokoulin-Petrukhin model for the energy loss of muons due to direct e+e- pair production. 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 37 of file KokoulinPetrukhinModel.h.

Constructor & Destructor Documentation

KokoulinPetrukhinModel::KokoulinPetrukhinModel ( )

Definition at line 22 of file KokoulinPetrukhinModel.cxx.

22  :
23 MuELossI("genie::mueloss::KokoulinPetrukhinModel")
24 {
25 
26 }
KokoulinPetrukhinModel::KokoulinPetrukhinModel ( string  config)

Definition at line 28 of file KokoulinPetrukhinModel.cxx.

28  :
29 MuELossI("genie::mueloss::KokoulinPetrukhinModel", config)
30 {
31 
32 }
KokoulinPetrukhinModel::~KokoulinPetrukhinModel ( )
virtual

Definition at line 34 of file KokoulinPetrukhinModel.cxx.

35 {
36 
37 }

Member Function Documentation

double KokoulinPetrukhinModel::dE_dx ( double  E,
MuELMaterial_t  material 
) const
virtual

Implement the MuELossI interface.

Implements genie::mueloss::MuELossI.

Definition at line 39 of file KokoulinPetrukhinModel.cxx.

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

40 {
41 // Calculate the muon -dE/dx due to e+e- pair production (in GeV^-2).
42 // To convert the result to more handly units, eg MeV/(gr/cm^2), just write:
43 // dE_dx /= (units::MeV/(units::g/units::cm2));
44 
45  if(material == eMuUndefined) return 0;
46  if(E<=MuELProcess::Threshold(this->Process()) || E>=kMaxMuE) return 0;
47 
48  // material Z,E
49  double Z = MuELMaterial::Z(material);
50  double A = MuELMaterial::A(material);
51 
52  // calculate (the min,max) fraction of energy, v, carried to the photon
53  double Vmin = 4.*kElectronMass/E;
54  double Vmax = 1. - 0.75*kSqrtNapierConst* (kMuonMass/E) * TMath::Power(Z,1/3.);
55 
56  // claculate the limits of the asymmetry parameter of the e+e- pair
57  // p = (E(+) - E(-)) / (E(+) + E(-))
58  double Pmin = 0.;
59  double Pmax = 1.;
60 
61  // integrate the Kokulin-Petrukhin differential cross section v*ds/dv for
62  // muon e+e- pair production over v and p
63 
64  ROOT::Math::IBaseFunctionMultiDim * integrand =
66  ROOT::Math::IntegrationMultiDim::Type ig_type =
68 
69  double abstol = 1; // We mostly care about relative tolerance.
70  double reltol = 1E-4;
71  int nmaxeval = 500000;
72 
73  ROOT::Math::IntegratorMultiDim ig(*integrand, ig_type, abstol, reltol, nmaxeval);
74 
75  double min[2] = { Vmin, Pmin };
76  double max[2] = { Vmax, Pmax };
77 
78  // calculate the b factor (bE = -dE/dx) in GeV^-3
79  A *= units::g;
80  double bpair = (2*kNA/A) * ig.Integral(min, max);
81 
82  delete integrand;
83 
84  // calculate the dE/dx due to muon bremsstrahlung in GeV^-2
85  double de_dx = bpair*E;
86  return de_dx;
87 }
static constexpr double g
Definition: Units.h:144
ROOT::Math::IntegrationMultiDim::Type IntegrationNDimTypeFromString(string type)
Definition: GSLUtils.cxx:59
static const double kSqrtNapierConst
static double Z(MuELMaterial_t material)
Definition: MuELMaterial.h:236
static const double kElectronMass
static double Threshold(MuELProcess_t p)
Definition: MuELProcess.h:58
static constexpr double A
Definition: Units.h:74
const double kMaxMuE
Definition: MuELossI.h:29
static double A(MuELMaterial_t material)
Definition: MuELMaterial.h:304
MuELProcess_t genie::mueloss::KokoulinPetrukhinModel::Process ( void  ) const
inlinevirtual

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