#include <IORTInteractionParameters.hh>

Inheritance diagram for IORTInteractionParameters:

Collaboration diagram for IORTInteractionParameters:

Public Member Functions
	IORTInteractionParameters (G4bool)

	~IORTInteractionParameters ()

G4bool	GetStoppingTable (const G4String &vararg)

G4double	GetStopping (G4double energy, const G4ParticleDefinition , const G4Material , G4double density=0.)

void	ListOfNistMaterials (const G4String &vararg)

void	BeamOn ()

bool	ParseArg (const G4String &vararg)

Public Member Functions inherited from G4EmCalculator
	G4EmCalculator ()

	~G4EmCalculator ()

G4double	GetDEDX (G4double kinEnergy, const G4ParticleDefinition , const G4Material , const G4Region *r=0)

G4double	GetDEDX (G4double kinEnergy, const G4String &part, const G4String &mat, const G4String &s="world")

G4double	GetRangeFromRestricteDEDX (G4double kinEnergy, const G4ParticleDefinition , const G4Material , const G4Region *r=0)

G4double	GetRangeFromRestricteDEDX (G4double kinEnergy, const G4String &part, const G4String &mat, const G4String &s="world")

G4double	GetCSDARange (G4double kinEnergy, const G4ParticleDefinition , const G4Material , const G4Region *r=0)

G4double	GetCSDARange (G4double kinEnergy, const G4String &part, const G4String &mat, const G4String &s="world")

G4double	GetRange (G4double kinEnergy, const G4ParticleDefinition , const G4Material , const G4Region *r=0)

G4double	GetRange (G4double kinEnergy, const G4String &part, const G4String &mat, const G4String &s="world")

G4double	GetKinEnergy (G4double range, const G4ParticleDefinition , const G4Material , const G4Region *r=0)

G4double	GetKinEnergy (G4double range, const G4String &part, const G4String &mat, const G4String &s="world")

G4double	GetCrossSectionPerVolume (G4double kinEnergy, const G4ParticleDefinition , const G4String &processName, const G4Material , const G4Region *r=0)

G4double	GetCrossSectionPerVolume (G4double kinEnergy, const G4String &part, const G4String &proc, const G4String &mat, const G4String &s="world")

G4double	GetShellIonisationCrossSectionPerAtom (const G4String &part, G4int Z, G4AtomicShellEnumerator shell, G4double kinEnergy)

G4double	GetMeanFreePath (G4double kinEnergy, const G4ParticleDefinition , const G4String &processName, const G4Material , const G4Region *r=0)

G4double	GetMeanFreePath (G4double kinEnergy, const G4String &part, const G4String &proc, const G4String &mat, const G4String &s="world")

void	PrintDEDXTable (const G4ParticleDefinition *)

void	PrintRangeTable (const G4ParticleDefinition *)

void	PrintInverseRangeTable (const G4ParticleDefinition *)

G4double	ComputeDEDX (G4double kinEnergy, const G4ParticleDefinition , const G4String &processName, const G4Material , G4double cut=DBL_MAX)

G4double	ComputeDEDX (G4double kinEnergy, const G4String &part, const G4String &proc, const G4String &mat, G4double cut=DBL_MAX)

G4double	ComputeElectronicDEDX (G4double kinEnergy, const G4ParticleDefinition , const G4Material mat, G4double cut=DBL_MAX)

G4double	ComputeElectronicDEDX (G4double kinEnergy, const G4String &part, const G4String &mat, G4double cut=DBL_MAX)

G4double	ComputeDEDXForCutInRange (G4double kinEnergy, const G4ParticleDefinition , const G4Material mat, G4double rangecut=DBL_MAX)

G4double	ComputeDEDXForCutInRange (G4double kinEnergy, const G4String &part, const G4String &mat, G4double rangecut=DBL_MAX)

G4double	ComputeNuclearDEDX (G4double kinEnergy, const G4ParticleDefinition , const G4Material )

G4double	ComputeNuclearDEDX (G4double kinEnergy, const G4String &part, const G4String &mat)

G4double	ComputeTotalDEDX (G4double kinEnergy, const G4ParticleDefinition , const G4Material , G4double cut=DBL_MAX)

G4double	ComputeTotalDEDX (G4double kinEnergy, const G4String &part, const G4String &mat, G4double cut=DBL_MAX)

G4double	ComputeCrossSectionPerVolume (G4double kinEnergy, const G4ParticleDefinition , const G4String &processName, const G4Material , G4double cut=0.0)

G4double	ComputeCrossSectionPerVolume (G4double kinEnergy, const G4String &part, const G4String &proc, const G4String &mat, G4double cut=0.0)

G4double	ComputeCrossSectionPerAtom (G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, G4double Z, G4double A, G4double cut=0.0)

G4double	ComputeCrossSectionPerAtom (G4double kinEnergy, const G4String &part, const G4String &processName, const G4Element *, G4double cut=0.0)

G4double	ComputeCrossSectionPerShell (G4double kinEnergy, const G4ParticleDefinition *, const G4String &processName, G4int Z, G4int shellIdx, G4double cut=0.0)

G4double	ComputeCrossSectionPerShell (G4double kinEnergy, const G4String &part, const G4String &processName, const G4Element *, G4int shellIdx, G4double cut=0.0)

G4double	ComputeGammaAttenuationLength (G4double kinEnergy, const G4Material *)

G4double	ComputeShellIonisationCrossSectionPerAtom (const G4String &part, G4int Z, G4AtomicShellEnumerator shell, G4double kinEnergy, const G4Material *mat=0)

G4double	ComputeMeanFreePath (G4double kinEnergy, const G4ParticleDefinition , const G4String &processName, const G4Material , G4double cut=0.0)

G4double	ComputeMeanFreePath (G4double kinEnergy, const G4String &, const G4String &, const G4String &processName, G4double cut=0.0)

G4double	ComputeEnergyCutFromRangeCut (G4double range, const G4ParticleDefinition , const G4Material )

G4double	ComputeEnergyCutFromRangeCut (G4double range, const G4String &, const G4String &)

const G4ParticleDefinition *	FindParticle (const G4String &)

const G4ParticleDefinition *	FindIon (G4int Z, G4int A)

const G4Material *	FindMaterial (const G4String &)

const G4Region *	FindRegion (const G4String &)

const G4MaterialCutsCouple *	FindCouple (const G4Material , const G4Region r=0)

G4VProcess *	FindProcess (const G4ParticleDefinition *part, const G4String &processName)

void	SetupMaterial (const G4Material *)

void	SetupMaterial (const G4String &)

void	SetVerbose (G4int val)

Private Member Functions
G4Material *	GetNistMaterial (G4String material)

Private Attributes
G4NistElementBuilder *	nistEle

G4NistMaterialBuilder *	nistMat

std::ofstream	outfile

std::ostream	data

G4Material *	Pmaterial

IORTParameterMessenger *	pMessenger

bool	beamFlag

G4double	kinEmin

G4double	kinEmax

G4double	npoints

G4String	particle

G4String	material

G4String	filename

G4double	dedxtot

G4double	density

std::vector< G4double >	energy

std::vector< G4double >	massDedx

Detailed Description

Definition at line 64 of file IORTInteractionParameters.hh.

Constructor & Destructor Documentation

◆ IORTInteractionParameters()

IORTInteractionParameters::IORTInteractionParameters ( G4bool wantMessenger )

Definition at line 65 of file IORTInteractionParameters.cc.

                                                                         : 
     nistEle(new G4NistElementBuilder(0)),                                         
     nistMat(new G4NistMaterialBuilder(nistEle, 0)),                                   
     data(G4cout.rdbuf()), 
     pMessenger(0),
     beamFlag(false)
 #ifdef G4ANALYSIS_USE_ROOT 
     ,theRootCanvas(0),
     theRootGraph(0)
 #endif
 {
     if (wantMessenger) pMessenger = new IORTParameterMessenger(this); 
 }

◆ ~IORTInteractionParameters()

IORTInteractionParameters::~IORTInteractionParameters ( )

Definition at line 79 of file IORTInteractionParameters.cc.

 {
     if (pMessenger) delete pMessenger; 
     delete nistMat; 
     delete nistEle; 
 }

Member Function Documentation

◆ BeamOn()

void IORTInteractionParameters::BeamOn ( )

Definition at line 256 of file IORTInteractionParameters.cc.

 {
     // first check if RunManager is above G4State_Idle 
     G4StateManager* mState = G4StateManager::GetStateManager();
     G4ApplicationState  aState = mState -> GetCurrentState(); 
     if ( aState <= G4State_Idle && beamFlag == false)
      {
         G4cout << "Issuing a G4RunManager::beamOn()... "; 
         G4cout << "Current Run State is " << mState -> GetStateString( aState ) << G4endl; 
         G4RunManager::GetRunManager() -> BeamOn(0);
         beamFlag = true;
      }
 
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetNistMaterial()

G4Material * IORTInteractionParameters::GetNistMaterial ( G4String material )

private

Definition at line 195 of file IORTInteractionParameters.cc.

 {
     Pmaterial = G4NistManager::Instance()->FindOrBuildMaterial(mat);
     if (Pmaterial) density = Pmaterial -> GetDensity(); 
     return Pmaterial;
 }

Here is the call graph for this function:

Here is the caller graph for this function:

◆ GetStopping()

G4double IORTInteractionParameters::GetStopping	(	G4double	energy,
		const G4ParticleDefinition *	pDef,
		const G4Material *	pMat,
		G4double	density = `0.`
	)

Definition at line 86 of file IORTInteractionParameters.cc.

 {
     if (dens) return ComputeTotalDEDX(ene, pDef, pMat)/dens;
     return ComputeTotalDEDX(ene, pDef, pMat);
 }

Here is the call graph for this function:

◆ GetStoppingTable()

bool IORTInteractionParameters::GetStoppingTable ( const G4String & vararg )

Definition at line 94 of file IORTInteractionParameters.cc.

 {
   // Check arguments
   if ( !ParseArg(vararg)) return false;
   // Clear previous energy & mass sp vectors
   energy.clear(); 
   massDedx.clear();
   // log scale 
   if (kinEmin != kinEmax && npoints >1)
     {
       G4double logmin = std::log10(kinEmin);
       G4double logmax = std::log10(kinEmax); 
       G4double en;
       // uniform log space
       for (G4double c = 0.; c < npoints; c++)
     {
       en = std::pow(10., logmin + ( c*(logmax-logmin)  / (npoints - 1.)) );  
       energy.push_back(en/MeV);
       dedxtot =  ComputeTotalDEDX (en, particle, material);
       massDedx.push_back ( (dedxtot / density)/(MeV*cm2/g) );
     }
     }
   else // one point only
     {
       energy.push_back(kinEmin/MeV);
       dedxtot =  ComputeTotalDEDX (kinEmin, particle, material);
       massDedx.push_back ( (dedxtot / density)/(MeV*cm2/g) );
     }
   
   G4cout.precision(6);  
   data <<  "MeV             " << "MeV*cm2/g      " << particle << " (into " << 
     material << ", density = " << G4BestUnit(density,"Volumic Mass") << ")" << G4endl;
   data << G4endl;
   data << std::left << std::setfill(' ');
   for (size_t i=0; i<energy.size(); i++){
     data << std::setw(16) << energy[i] << massDedx[i] << G4endl;
   }
   outfile.close();
   // This will plot 
 #ifdef G4ANALYSIS_USE_ROOT 
   PlotStopping("pdf");
 #endif
 
   // Info to user
   G4String ofName = (filename == "") ? "User terminal": filename;
   G4cout << "User choice:\n";
   G4cout << "Kinetic energy lower limit= "<< G4BestUnit(kinEmin,"Energy") << 
     ", Kinetic energy upper limit= " << G4BestUnit(kinEmax,"Energy") << 
     ", npoints= "<< npoints << ", particle= \"" << particle << 
     "\", material= \"" << material << "\", filename= \""<< 
     ofName << "\"" << G4endl;
   return true;
 }

Here is the call graph for this function:

◆ ListOfNistMaterials()

void IORTInteractionParameters::ListOfNistMaterials ( const G4String & vararg )

Definition at line 271 of file IORTInteractionParameters.cc.

 {
 /*
  $G4INSTALL/source/materials/src/G4NistElementBuilder.cc
  You can also construct a new material by the ConstructNewMaterial method:
  see $G4INSTALL/source/materials/src/G4NistMaterialBuilder.cc
 */
   // Get simplest full list
   if (vararg =="list")
     {
       const std::vector<G4String>& vec =  nistMat -> GetMaterialNames(); 
       for (size_t i=0; i<vec.size(); i++)
     {
       G4cout << std::setw(12) << std::left << i+1 << vec[i] << G4endl;
     }
       G4cout << G4endl;
     }
   else if (vararg =="all" || vararg =="simple" || 
        vararg =="compound" || vararg =="hep" )
     {
       nistMat -> ListMaterials(vararg);
     }
 }

◆ ParseArg()

bool IORTInteractionParameters::ParseArg ( const G4String & vararg )

Definition at line 202 of file IORTInteractionParameters.cc.

 {
   kinEmin = kinEmax = npoints = 0.;
   particle = material = filename = "";
   // set internal variables
   std::istringstream strParam(vararg);
   // TODO here check for number and parameters consistency 
   strParam >> std::skipws >> material >> kinEmin >> kinEmax >> npoints >> particle >> filename;
   // npoints must be an integer!
   npoints = std::floor(npoints); 
 
 // Check that  kinEmax >= kinEmin > 0 &&  npoints >= 1 
 // TODO NIST points and linear scale
    if (kinEmax == 0. && kinEmin > 0. ) kinEmax = kinEmin;
    if (kinEmax == 0. && kinEmin == 0. ) kinEmax = kinEmin = 1.*MeV;
    if (kinEmax < kinEmin) 
    {
      G4cout << "WARNING: kinEmin must not exceed kinEmax!" << G4endl;
      G4cout << "Usage: /parameter/command  material kinetic Emin kinetic Emax nPoints [particle] [output filename]" << G4endl;  
      return false;
    }
   if (npoints < 1) npoints = 1;
 
     // check if element/material is into database
     if (!GetNistMaterial(material) )
        {
          G4cout << "WARNING: material \"" << material << "\" doesn't exist in NIST elements/materials"
                " table [$G4INSTALL/source/materials/src/G4NistMaterialBuilder.cc]" << G4endl; 
          G4cout << "Use command \"/parameter/nist\" to see full materials list" << G4endl; 
          return false;
        }
     // Check for particle
     if (particle == "") particle = "proton"; // default to "proton"
     else if ( !FindParticle(particle) )
        {
         G4cout << "WARNING: Particle \"" << particle << "\" isn't supported." << G4endl;
         G4cout << "Try the command \"/particle/list\" to get full supported particles list." << G4endl;
         G4cout << "If you are interested in an ion that isn't in this list you must give it to the particle gun."
                   "\nTry the commands:\n/gun/particle ion"
               "\n/gun/ion <atomic number> <mass number> <[charge]>" << G4endl << G4endl;
         return false;
        }
     // start physics by forcing a G4RunManager::BeamOn(): 
     BeamOn();
     // Set output file
     if( filename != "" ) 
        {
           outfile.open(filename,std::ios_base::trunc); // overwrite existing file
           data.rdbuf(outfile.rdbuf());
        }
     else data.rdbuf(G4cout.rdbuf());    // output is G4cout!                
     return true;
 }