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

Private Member Functions
G4double	ComputeMicroscopicCrossSection (G4double incidentKE, G4double A)
G4KokoulinMuonNuclearXS &	operator= (const G4KokoulinMuonNuclearXS &right)
	G4KokoulinMuonNuclearXS (const G4KokoulinMuonNuclearXS &)

Private Attributes
G4double	LowestKineticEnergy
G4double	HighestKineticEnergy
G4int	TotBin
G4double	CutFixed
G4bool	isInitialized
G4bool	isMaster

Static Private Attributes
static G4PhysicsVector *	theCrossSection [MAXZMUN] = {0}

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() [1/2]

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

Here is the caller graph for this function:

~G4KokoulinMuonNuclearXS()

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

G4KokoulinMuonNuclearXS() [2/2]

G4KokoulinMuonNuclearXS::G4KokoulinMuonNuclearXS ( const G4KokoulinMuonNuclearXS &  )

private

Member Function Documentation

BuildCrossSectionTable()

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:

BuildPhysicsTable()

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:

Here is the caller graph for this function:

ComputeDDMicroscopicCrossSection()

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:

Here is the caller graph for this function:

ComputeMicroscopicCrossSection()

G4double G4KokoulinMuonNuclearXS::ComputeMicroscopicCrossSection ( G4double  incidentKE, G4double  A  )

private

Definition at line 134 of file G4KokoulinMuonNuclearXS.cc.

{
  // Calculate cross section (differential in muon incident kinetic energy) by 
  // integrating the double differential cross section over the energy loss

  static const G4double xgi[] = 
    {0.0199,0.1017,0.2372,0.4083,0.5917,0.7628,0.8983,0.9801};
  static const G4double wgi[] = 
    {0.0506,0.1112,0.1569,0.1813,0.1813,0.1569,0.1112,0.0506};
  static const G4double ak1 = 6.9;
  static const G4double ak2 = 1.0;

  G4double Mass = G4MuonMinus::MuonMinus()->GetPDGMass();

  G4double CrossSection = 0.0;
  if (KineticEnergy <= CutFixed) return CrossSection; 

  G4double epmin = CutFixed;
  G4double epmax = KineticEnergy + Mass - 0.5*proton_mass_c2;
  if (epmax <= epmin) return CrossSection; // NaN bug correction

  G4double aaa = G4Log(epmin);
  G4double bbb = G4Log(epmax);
  G4int kkk = std::max(1,G4int((bbb-aaa)/ak1 +ak2));
  G4double hhh = (bbb-aaa)/kkk ;
  G4double epln;
  G4double ep;
  G4double x;

  for (G4int l = 0; l < kkk; ++l) {
    x = aaa + hhh*l;
    for (G4int ll = 0; ll < 8; ++ll) {
      epln=x+xgi[ll]*hhh;
      ep = G4Exp(epln);
      CrossSection += 
    ep*wgi[ll]*ComputeDDMicroscopicCrossSection(KineticEnergy, 0, A, ep);
    }
  }

  CrossSection *= hhh ;
  if (CrossSection < 0.) { CrossSection = 0.; }
  return CrossSection;
}

Here is the call graph for this function:

Here is the caller graph for this function:

CrossSectionDescription()

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

Here is the caller graph for this function:

Default_Name()

static const char* G4KokoulinMuonNuclearXS::Default_Name ( )

inlinestatic

Definition at line 60 of file G4KokoulinMuonNuclearXS.hh.

{return "KokoulinMuonNuclearXS";}

Here is the call graph for this function:

Here is the caller graph for this function:

GetElementCrossSection()

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:

Here is the caller graph for this function:

IsElementApplicable()

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

virtual

Reimplemented from G4VCrossSectionDataSet.

Definition at line 89 of file G4KokoulinMuonNuclearXS.cc.

{
  return true;
}

Here is the caller graph for this function:

operator=()

G4KokoulinMuonNuclearXS& G4KokoulinMuonNuclearXS::operator= ( const G4KokoulinMuonNuclearXS &  right )

private

Here is the caller graph for this function:

Member Data Documentation

CutFixed

G4double G4KokoulinMuonNuclearXS::CutFixed

private

Definition at line 91 of file G4KokoulinMuonNuclearXS.hh.

HighestKineticEnergy

G4double G4KokoulinMuonNuclearXS::HighestKineticEnergy

private

Definition at line 89 of file G4KokoulinMuonNuclearXS.hh.

isInitialized

G4bool G4KokoulinMuonNuclearXS::isInitialized

private

Definition at line 92 of file G4KokoulinMuonNuclearXS.hh.

isMaster

G4bool G4KokoulinMuonNuclearXS::isMaster

private

Definition at line 93 of file G4KokoulinMuonNuclearXS.hh.

LowestKineticEnergy

G4double G4KokoulinMuonNuclearXS::LowestKineticEnergy

private

Definition at line 88 of file G4KokoulinMuonNuclearXS.hh.

theCrossSection

G4PhysicsVector * G4KokoulinMuonNuclearXS::theCrossSection = {0}

staticprivate

Definition at line 86 of file G4KokoulinMuonNuclearXS.hh.

TotBin

G4int G4KokoulinMuonNuclearXS::TotBin

private

Definition at line 90 of file G4KokoulinMuonNuclearXS.hh.

---

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

• G4KokoulinMuonNuclearXS.hh
• G4KokoulinMuonNuclearXS.cc