MuCrossSections Class Reference

#include <MuCrossSections.hh>

Public Member Functions
	MuCrossSections ()
	~MuCrossSections ()
G4double	CR_Macroscopic (const G4String &, G4Material *, G4double, G4double)
G4double	CR_PerAtom (const G4String &, G4Element *, G4double, G4double)

Detailed Description

Definition at line 44 of file MuCrossSections.hh.

Constructor & Destructor Documentation

MuCrossSections::MuCrossSections

(

)

Definition at line 44 of file MuCrossSections.cc.

{ }

MuCrossSections::~MuCrossSections

(

)

Definition at line 49 of file MuCrossSections.cc.

{ }

Member Function Documentation

G4double MuCrossSections::CR_Macroscopic	(	const G4String &	process,
		G4Material *	material,
		G4double	tkin,
		G4double	ep
	)

Definition at line 54 of file MuCrossSections.cc.

References G4Material::GetElementVector(), G4Material::GetNumberOfElements(), and G4Material::GetVecNbOfAtomsPerVolume().

{
  const G4ElementVector* theElementVector = material->GetElementVector();
  const G4double* NbOfAtomsPerVolume = material->GetVecNbOfAtomsPerVolume();

  G4double SIGMA = 0 ;

  for ( size_t i=0 ; i < material->GetNumberOfElements() ; i++ )
  {
    G4Element* element = (*theElementVector)[i];
    SIGMA += NbOfAtomsPerVolume[i] * CR_PerAtom(process, element, tkin, ep);
  }
  return SIGMA;
}

G4double MuCrossSections::CR_PerAtom	(	const G4String &	process,
		G4Element *	element,
		G4double	tkin,
		G4double	ep
	)

Definition at line 74 of file MuCrossSections.cc.

References test::a, python.hepunit::Avogadro, python.hepunit::cm2, g(), G4Element::GetA(), G4Element::GetZ(), python.hepunit::GeV, python.hepunit::mole, and z.

{
 G4double z = element->GetZ();
 G4double a = element->GetA();
 
 G4double sigma = 0.;
 if (process == "muBrems")
   sigma = CRB_Mephi(z,a/(g/mole),tkin/GeV,ep/GeV)*(cm2/(g*GeV))*a/Avogadro;
   
 else if (process == "muIoni")
   sigma = CRK_Mephi(z,a/(g/mole),tkin/GeV,ep/GeV)*(cm2/(g*GeV))*a/Avogadro;
   
 //else if (process == "muNucl")
 else if (process == "muonNuclear")
   sigma = CRN_Mephi(z,a/(g/mole),tkin/GeV,ep/GeV)*(cm2/(g*GeV))*a/Avogadro;
   
 else if (process == "muPairProd")
   sigma = CRP_Mephi(z,a/(g/mole),tkin/GeV,ep/GeV)*(cm2/(g*GeV))*a/Avogadro;
   
 return sigma;        
}

---

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

• MuCrossSections.hh
• MuCrossSections.cc