#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	IsElementApplicable (const G4DynamicParticle , G4int Z, const G4Material mat=0)

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	GetElementCrossSection (const G4DynamicParticle , G4int Z, 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

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::G4KokoulinMuonNuclearXS ( )

Definition at line 47 of file G4KokoulinMuonNuclearXS.cc.

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

◆ ~G4KokoulinMuonNuclearXS()

G4KokoulinMuonNuclearXS::~G4KokoulinMuonNuclearXS ( )

virtual

Definition at line 55 of file G4KokoulinMuonNuclearXS.cc.

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

Member Function Documentation

◆ BuildCrossSectionTable()

void G4KokoulinMuonNuclearXS::BuildCrossSectionTable ( )

Definition at line 71 of file G4KokoulinMuonNuclearXS.cc.

{
  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;
  }
}

Referenced by G4KokoulinMuonNuclearXS().

◆ ComputeDDMicroscopicCrossSection()

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

Definition at line 153 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)
 
  const G4double alam2 = 0.400*GeV*GeV;
  const G4double alam  = 0.632456*GeV;
  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;
}

◆ GetElementCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 194 of file G4KokoulinMuonNuclearXS.cc.

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

◆ IsElementApplicable()

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

◆ G4KokoulinMuonNuclearXS()

◆ ~G4KokoulinMuonNuclearXS()

Member Function Documentation

◆ BuildCrossSectionTable()

◆ ComputeDDMicroscopicCrossSection()

◆ GetElementCrossSection()

◆ IsElementApplicable()