#include <G4KokoulinMuonNuclearXS.hh>

Inheritance diagram for G4KokoulinMuonNuclearXS:

Collaboration diagram for G4KokoulinMuonNuclearXS:

Public Member Functions
	G4KokoulinMuonNuclearXS ()

virtual	~G4KokoulinMuonNuclearXS ()

virtual void	CrossSectionDescription (std::ostream &) const

G4bool	IsElementApplicable (const G4DynamicParticle particle, G4int Z, const G4Material )

G4double	GetElementCrossSection (const G4DynamicParticle particle, G4int Z, const G4Material )

void	BuildPhysicsTable (const G4ParticleDefinition &)

void	BuildCrossSectionTable ()

G4double	ComputeDDMicroscopicCrossSection (G4double incidentKE, G4double Z, G4double A, G4double epsilon)

Public Member Functions inherited from G4VCrossSectionDataSet
	G4VCrossSectionDataSet (const G4String &nam="")

virtual	~G4VCrossSectionDataSet ()

virtual G4bool	IsIsoApplicable (const G4DynamicParticle , G4int Z, G4int A, const G4Element elm=0, const G4Material *mat=0)

G4double	GetCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

G4double	ComputeCrossSection (const G4DynamicParticle , const G4Element , const G4Material *mat=0)

virtual G4double	GetIsoCrossSection (const G4DynamicParticle , G4int Z, G4int A, const G4Isotope iso=0, const G4Element elm=0, const G4Material mat=0)

virtual G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual G4int	GetVerboseLevel () const

virtual void	SetVerboseLevel (G4int value)

G4double	GetMinKinEnergy () const

void	SetMinKinEnergy (G4double value)

G4double	GetMaxKinEnergy () const

void	SetMaxKinEnergy (G4double value)

const G4String &	GetName () const

Static Public Member Functions
static const char *	Default_Name ()

Additional Inherited Members
Protected Member Functions inherited from G4VCrossSectionDataSet
void	SetName (const G4String &)

Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel

Detailed Description

Definition at line 53 of file G4KokoulinMuonNuclearXS.hh.

Constructor & Destructor Documentation

G4KokoulinMuonNuclearXS::G4KokoulinMuonNuclearXS ( )

Definition at line 59 of file G4KokoulinMuonNuclearXS.cc.

   :G4VCrossSectionDataSet(Default_Name()), 
   LowestKineticEnergy(1*GeV), HighestKineticEnergy(1*PeV),
   TotBin(60), CutFixed(0.2*GeV), isInitialized(false), isMaster(false)
 {}

G4KokoulinMuonNuclearXS::~G4KokoulinMuonNuclearXS ( )

virtual

Definition at line 65 of file G4KokoulinMuonNuclearXS.cc.

 {
   if (isMaster) {
     for(G4int i=0; i<MAXZMUN; ++i) {
       delete theCrossSection[i];
       theCrossSection[i] = 0;
     }
   }
 }

Member Function Documentation

void G4KokoulinMuonNuclearXS::BuildCrossSectionTable ( )

Definition at line 108 of file G4KokoulinMuonNuclearXS.cc.

 {
   G4double energy, A, Value;
   G4int Z;
 
   G4int nEl = G4Element::GetNumberOfElements(); 
   const G4ElementTable* theElementTable = G4Element::GetElementTable();
   G4NistManager* nistManager = G4NistManager::Instance();
 
   for (G4int j = 0; j < nEl; j++) {
     Z = G4lrint((*theElementTable)[j]->GetZ());
     A  = nistManager->GetAtomicMassAmu(Z);
     if(Z < MAXZMUN && !theCrossSection[Z]) {
       theCrossSection[Z] = new G4PhysicsLogVector(LowestKineticEnergy,
                           HighestKineticEnergy, 
                           TotBin);
       for (G4int i = 0; i <= TotBin; ++i) {
     energy = theCrossSection[Z]->Energy(i);
     Value = ComputeMicroscopicCrossSection(energy, A);
     theCrossSection[Z]->PutValue(i,Value);
       }
     }
   }
 }

Here is the call graph for this function:

Here is the caller graph for this function:

void G4KokoulinMuonNuclearXS::BuildPhysicsTable ( const G4ParticleDefinition & )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 96 of file G4KokoulinMuonNuclearXS.cc.

 {
   if(!isInitialized) { 
     isInitialized = true; 
     for(G4int i=0; i<MAXZMUN; ++i) {
       if(theCrossSection[i]) { return; }
     }
     isMaster = true; 
   }
   if(isMaster) { BuildCrossSectionTable(); }
 }

Here is the call graph for this function:

G4double G4KokoulinMuonNuclearXS::ComputeDDMicroscopicCrossSection	(	G4double	incidentKE,
		G4double	Z,
		G4double	A,
		G4double	epsilon
	)

Definition at line 179 of file G4KokoulinMuonNuclearXS.cc.

 {
   // Calculate the double-differential microscopic cross section (in muon
   // incident kinetic energy and energy loss) using the cross section formula
   // of R.P. Kokoulin (18/01/98)
 
   static const G4double alam2 = 0.400*GeV*GeV;
   static const G4double alam  = 0.632456*GeV;
   static const G4double coeffn = fine_structure_const/pi;   
 
   G4double ParticleMass = G4MuonMinus::MuonMinus()->GetPDGMass();
   G4double TotalEnergy = KineticEnergy + ParticleMass;
 
   G4double DCrossSection = 0.;
 
   if ((epsilon >= TotalEnergy - 0.5*proton_mass_c2) ||
       (epsilon <= CutFixed) ) { return DCrossSection; }
 
   G4double ep = epsilon/GeV;
   G4double aeff = 0.22*A+0.78*G4Exp(0.89*G4Log(A));       //shadowing 
   G4double sigph = (49.2+11.1*G4Log(ep)+151.8/std::sqrt(ep))*microbarn; 
   
   G4double v = epsilon/TotalEnergy;
   G4double v1 = 1.-v;
   G4double v2 = v*v;
   G4double mass2 = ParticleMass*ParticleMass;
 
   G4double up = TotalEnergy*TotalEnergy*v1/mass2*(1.+mass2*v2/(alam2*v1));
   G4double down = 1.+epsilon/alam*(1.+alam/(2.*proton_mass_c2)+epsilon/alam);
 
   DCrossSection = coeffn*aeff*sigph/epsilon*
                   (-v1+(v1+0.5*v2*(1.+2.*mass2/alam2))*G4Log(up/down));
 
   if (DCrossSection < 0.) { DCrossSection = 0.; }
   return DCrossSection;
 }

Here is the call graph for this function:

void G4KokoulinMuonNuclearXS::CrossSectionDescription ( std::ostream & outFile ) const

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 77 of file G4KokoulinMuonNuclearXS.cc.

 {
     outFile << "G4KokoulinMuonNuclearXS provides the total inelastic\n"
     << "cross section for mu- and mu+ interactions with nuclei.\n"
     << "R. Kokoulin's approximation of the Borog and Petrukhin double\n"
     << "differential cross section at high energy and low Q**2 is integrated\n"
     << "over the muon energy loss to get the total cross section as a\n"
     << "function of muon kinetic energy\n" ;
 }

static const char* G4KokoulinMuonNuclearXS::Default_Name ( )

inlinestatic

Definition at line 60 of file G4KokoulinMuonNuclearXS.hh.

60 {return "KokoulinMuonNuclearXS";}

Here is the caller graph for this function:

G4double G4KokoulinMuonNuclearXS::GetElementCrossSection	(	const G4DynamicParticle *	particle,
		G4int	Z,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 218 of file G4KokoulinMuonNuclearXS.cc.

 {
   return theCrossSection[Z]->Value(aPart->GetKineticEnergy());
 }

Here is the call graph for this function:

G4bool G4KokoulinMuonNuclearXS::IsElementApplicable	(	const G4DynamicParticle *	particle,
		G4int	Z,
		const G4Material *
	)

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 89 of file G4KokoulinMuonNuclearXS.cc.

 {
   return true;
 }

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

geant4.10.03.p01/source/processes/hadronic/cross_sections/include/G4KokoulinMuonNuclearXS.hh
geant4.10.03.p01/source/processes/hadronic/cross_sections/src/G4KokoulinMuonNuclearXS.cc

Public Member Functions

Static Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation