#include <G4KokoulinMuonNuclearXS.hh>

Inheritance diagram for G4KokoulinMuonNuclearXS:

Public Member Functions
	G4KokoulinMuonNuclearXS ()

virtual	~G4KokoulinMuonNuclearXS ()

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

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

void	BuildCrossSectionTable ()

G4double	ComputeDDMicroscopicCrossSection (G4double incidentKE, G4double, G4double AtomicWeight, 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	BuildPhysicsTable (const G4ParticleDefinition &)

virtual void	DumpPhysicsTable (const G4ParticleDefinition &)

virtual void	CrossSectionDescription (std::ostream &) const

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

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 55 of file G4KokoulinMuonNuclearXS.hh.

Constructor & Destructor Documentation

G4KokoulinMuonNuclearXS::G4KokoulinMuonNuclearXS ( )

Definition at line 47 of file G4KokoulinMuonNuclearXS.cc.

References BuildCrossSectionTable().

  :G4VCrossSectionDataSet("KokoulinMuonNuclearXS"), theCrossSectionTable(0),
   LowestKineticEnergy(1*GeV), HighestKineticEnergy(1*PeV),
   TotBin(60), CutFixed(0.2*GeV)
 {
   BuildCrossSectionTable();
 }

G4KokoulinMuonNuclearXS::~G4KokoulinMuonNuclearXS ( )

virtual

Definition at line 55 of file G4KokoulinMuonNuclearXS.cc.

References G4PhysicsTable::clearAndDestroy().

 {
   if (theCrossSectionTable) {
     theCrossSectionTable->clearAndDestroy();
     delete theCrossSectionTable;
   }
 }

Member Function Documentation

void G4KokoulinMuonNuclearXS::BuildCrossSectionTable ( )

Definition at line 71 of file G4KokoulinMuonNuclearXS.cc.

References G4PhysicsVector::Energy(), G4Element::GetElementTable(), G4Element::GetNumberOfElements(), G4PhysicsTable::insertAt(), and G4PhysicsVector::PutValue().

Referenced by G4KokoulinMuonNuclearXS().

 {
   G4double lowEdgeEnergy;
   G4PhysicsLogVector* ptrVector;
 
   G4int nEl = G4Element::GetNumberOfElements(); 
   const G4ElementTable* theElementTable = G4Element::GetElementTable();
 
   theCrossSectionTable = new G4PhysicsTable(nEl);
 
   G4double AtomicNumber;
   G4double AtomicWeight;
   G4double Value;
 
   for (G4int j = 0; j < nEl; j++) {
     ptrVector = new G4PhysicsLogVector(LowestKineticEnergy,
                                        HighestKineticEnergy, TotBin);
     AtomicNumber = (*theElementTable)[j]->GetZ();
     AtomicWeight = (*theElementTable)[j]->GetA();
 
     for (G4int i = 0; i <= TotBin; ++i) {
       lowEdgeEnergy = ptrVector->Energy(i);
       Value = ComputeMicroscopicCrossSection(lowEdgeEnergy,
                                              AtomicNumber, AtomicWeight);
       ptrVector->PutValue(i,Value);
     }
 
     theCrossSectionTable->insertAt(j, ptrVector);
   }
 
   // Build (Z,element) look-up table (for use in GetZandACrossSection) 
   for (G4int j = 0; j < nEl; j++) {
     G4int ZZ = G4int((*theElementTable)[j]->GetZ());
     zelMap[ZZ] = j;
   }
 }

G4double G4KokoulinMuonNuclearXS::ComputeDDMicroscopicCrossSection	(	G4double	incidentKE,
		G4double	,
		G4double	AtomicWeight,
		G4double	epsilon
	)

Definition at line 154 of file G4KokoulinMuonNuclearXS.cc.

References test::a, python.hepunit::fine_structure_const, g(), G4ParticleDefinition::GetPDGMass(), python.hepunit::GeV, python.hepunit::microbarn, python.hepunit::mole, G4MuonMinus::MuonMinus(), python.hepunit::pi, python.hepunit::proton_mass_c2, and test::v.

 {
   // 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 a = AtomicWeight/(g/mole);
   G4double aeff = 0.22*a+0.78*std::exp(0.89*std::log(a));       //shadowing 
   G4double sigph = (49.2+11.1*std::log(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))*std::log(up/down));
 
   if (DCrossSection < 0.) DCrossSection = 0.; 
   return DCrossSection;
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 195 of file G4KokoulinMuonNuclearXS.cc.

References G4DynamicParticle::GetKineticEnergy().

 {
   G4int j = zelMap[ZZ];
   return (*theCrossSectionTable)[j]->Value(aPart->GetKineticEnergy());
 }

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 64 of file G4KokoulinMuonNuclearXS.cc.

 {
   return true;
 }

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

Public Member Functions

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation