G4KokoulinMuonNuclearXS Class Reference

#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.

{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:

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