G4KokoulinMuonNuclearXS Class Reference

#include <G4KokoulinMuonNuclearXS.hh>

Inheritance 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=nullptr, const G4Material *mat=nullptr)
G4double	GetCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
G4double	ComputeCrossSection (const G4DynamicParticle *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	ComputeCrossSectionPerElement (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, const G4Element *, const G4Material *mat=nullptr)
virtual G4double	GetIsoCrossSection (const G4DynamicParticle *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
virtual G4double	ComputeIsoCrossSection (G4double kinEnergy, G4double loge, const G4ParticleDefinition *, G4int Z, G4int A, const G4Isotope *iso=nullptr, const G4Element *elm=nullptr, const G4Material *mat=nullptr)
virtual const G4Isotope *	SelectIsotope (const G4Element *, G4double kinEnergy, G4double logE)
virtual void	DumpPhysicsTable (const G4ParticleDefinition &)
virtual void	SetVerboseLevel (G4int value)
G4double	GetMinKinEnergy () const
void	SetMinKinEnergy (G4double value)
G4double	GetMaxKinEnergy () const
void	SetMaxKinEnergy (G4double value)
bool	ForAllAtomsAndEnergies () const
void	SetForAllAtomsAndEnergies (G4bool val)
const G4String &	GetName () const
void	SetName (const G4String &nam)
G4VCrossSectionDataSet &	operator= (const G4VCrossSectionDataSet &right)=delete
	G4VCrossSectionDataSet (const G4VCrossSectionDataSet &)=delete

Static Public Member Functions
static const char *	Default_Name ()

Additional Inherited Members
Protected Attributes inherited from G4VCrossSectionDataSet
G4int	verboseLevel
G4String	name

Detailed Description

Definition at line 52 of file G4KokoulinMuonNuclearXS.hh.

Constructor & Destructor Documentation

G4KokoulinMuonNuclearXS()

G4KokoulinMuonNuclearXS::G4KokoulinMuonNuclearXS

(

)

Definition at line 58 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::~G4KokoulinMuonNuclearXS ( )

virtual

Definition at line 64 of file G4KokoulinMuonNuclearXS.cc.

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

Member Function Documentation

BuildCrossSectionTable()

void G4KokoulinMuonNuclearXS::BuildCrossSectionTable

(

)

Definition at line 107 of file G4KokoulinMuonNuclearXS.cc.

{
  G4double energy, A, Value;
  G4int Z;
 
  std::size_t nEl = G4Element::GetNumberOfElements(); 
  const G4ElementTable* theElementTable = G4Element::GetElementTable();
  G4NistManager* nistManager = G4NistManager::Instance();
 
  for (std::size_t j = 0; j < nEl; ++j) {
    Z = G4lrint((*theElementTable)[j]->GetZ());
 
    //AR-24Apr2018 Switch to treat transuranic elements as uranium  
    const G4bool isHeavyElementAllowed = true; if ( isHeavyElementAllowed && Z>92 ) Z=92;
 
    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);
      }
    }
  }
}

Referenced by BuildPhysicsTable().

BuildPhysicsTable()

void G4KokoulinMuonNuclearXS::BuildPhysicsTable ( const G4ParticleDefinition & )

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 95 of file G4KokoulinMuonNuclearXS.cc.

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

ComputeDDMicroscopicCrossSection()

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

Definition at line 181 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;
}

CrossSectionDescription()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 76 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" ;
}

Default_Name()

static const char * G4KokoulinMuonNuclearXS::Default_Name ( )

inlinestatic

Definition at line 59 of file G4KokoulinMuonNuclearXS.hh.

{return "KokoulinMuonNuclearXS";}

Referenced by G4MuonVDNuclearModel::G4MuonVDNuclearModel().

GetElementCrossSection()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 220 of file G4KokoulinMuonNuclearXS.cc.

{
  //AR-24Apr2018 Switch to treat transuranic elements as uranium  
  G4int Z = std::min(ZZ, 92);
  return theCrossSection[Z]->LogVectorValue(aPart->GetKineticEnergy(),
                        aPart->GetLogKineticEnergy());
}

IsElementApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 88 of file G4KokoulinMuonNuclearXS.cc.

{
  return true;
}

---

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

• G4KokoulinMuonNuclearXS.hh
• G4KokoulinMuonNuclearXS.cc