G4ICRU90StoppingData Class Reference

#include <G4ICRU90StoppingData.hh>

Public Member Functions
	G4ICRU90StoppingData ()
	~G4ICRU90StoppingData ()
void	Initialise ()
G4double	GetElectronicDEDXforProton (const G4Material *, G4double kinEnergy) const
G4double	GetElectronicDEDXforAlpha (const G4Material *, G4double scaledKinEnergy) const
G4int	GetIndex (const G4Material *) const
G4int	GetIndex (const G4String &) const
G4double	GetElectronicDEDXforProton (G4int idx, G4double kinEnergy) const
G4double	GetElectronicDEDXforAlpha (G4int idx, G4double scaledKinEnergy) const
G4bool	IsApplicable (const G4Material *) const

Detailed Description

Definition at line 56 of file G4ICRU90StoppingData.hh.

Constructor & Destructor Documentation

G4ICRU90StoppingData()

G4ICRU90StoppingData::G4ICRU90StoppingData ( )

explicit

Definition at line 53 of file G4ICRU90StoppingData.cc.

                                           : isInitialized(false)
{
  // 1st initialisation 
  for(size_t i=0; i<nvectors; ++i) { 
    materials[i]    = nullptr; 
    sdata_proton[i] = nullptr; 
    sdata_alpha[i]  = nullptr;
  }
  FillData();
 
  Initialise();
}

~G4ICRU90StoppingData()

G4ICRU90StoppingData::~G4ICRU90StoppingData

(

)

Definition at line 68 of file G4ICRU90StoppingData.cc.

{
  for(size_t i=0; i<nvectors; ++i) { 
    delete sdata_proton[i]; 
    delete sdata_alpha[i];
  }
}

Member Function Documentation

GetElectronicDEDXforAlpha() [1/2]

G4double G4ICRU90StoppingData::GetElectronicDEDXforAlpha	(	const G4Material *	mat,
		G4double	scaledKinEnergy
	)		const

Definition at line 125 of file G4ICRU90StoppingData.cc.

{
  G4int idx = GetIndex(mat);
  return (idx < 0) ? 0.0 : GetDEDX(sdata_alpha[idx], scaledKinEnergy);
}

Referenced by G4BetheBlochModel::ComputeDEDXPerVolume().

GetElectronicDEDXforAlpha() [2/2]

G4double G4ICRU90StoppingData::GetElectronicDEDXforAlpha ( G4int  idx, G4double  scaledKinEnergy  ) const

inline

Definition at line 151 of file G4ICRU90StoppingData.hh.

{
  return (idx < 0 || idx >= nvectors) ? 0.0 
    : GetDEDX(sdata_alpha[idx], scaledKinEnergy);
}

GetElectronicDEDXforProton() [1/2]

G4double G4ICRU90StoppingData::GetElectronicDEDXforProton	(	const G4Material *	mat,
		G4double	kinEnergy
	)		const

Definition at line 116 of file G4ICRU90StoppingData.cc.

{
  G4int idx = GetIndex(mat);
  return (idx < 0) ? 0.0 : GetDEDX(sdata_proton[idx], kinEnergy);
}

Referenced by G4BetheBlochModel::ComputeDEDXPerVolume().

GetElectronicDEDXforProton() [2/2]

G4double G4ICRU90StoppingData::GetElectronicDEDXforProton ( G4int  idx, G4double  kinEnergy  ) const

inline

Definition at line 142 of file G4ICRU90StoppingData.hh.

{
  return (idx < 0 || idx >= nvectors) ? 0.0 
    : GetDEDX(sdata_proton[idx], kinEnergy);
}

GetIndex() [1/2]

G4int G4ICRU90StoppingData::GetIndex ( const G4Material *  mat ) const

inline

Definition at line 111 of file G4ICRU90StoppingData.hh.

{  
  G4int idx = -1;
  if (mat == materials[0]) { idx = 0; }
  else if(mat == materials[1]) { idx = 1; }
  else if(mat == materials[2]) { idx = 2; }
  return idx;
}

Referenced by G4BetheBlochModel::ComputeDEDXPerVolume(), G4BetheBlochModel::CorrectionsAlongStep(), GetElectronicDEDXforAlpha(), and GetElectronicDEDXforProton().

GetIndex() [2/2]

G4int G4ICRU90StoppingData::GetIndex ( const G4String &  nam ) const

inline

Definition at line 122 of file G4ICRU90StoppingData.hh.

{
  G4int idx = -1;
  if (nam == materials[0]->GetName()) { idx = 0; }
  else if(nam == materials[1]->GetName()) { idx = 1; }
  else if(nam == materials[2]->GetName()) { idx = 2; }
  return idx;
}

Initialise()

void G4ICRU90StoppingData::Initialise

(

)

Definition at line 78 of file G4ICRU90StoppingData.cc.

{
  if(isInitialized) { return; }
  // this method may be called several times during initialisation
  G4int nmat = G4Material::GetNumberOfMaterials();
  if(nmat == (G4int)nvectors) { return; }
 
  static const G4String nameNIST_ICRU90[3] = 
    {"G4_AIR","G4_WATER","G4_GRAPHITE"};
 
  // loop via material list to add extra data
  for(G4int i=0; i<nmat; ++i) {
    const G4Material* mat = (*(G4Material::GetMaterialTable()))[i];
 
    G4bool isThere = false;  
    for(G4int j=0; j<nvectors; ++j) {
      if(mat == materials[j]) {
    isThere = true;
    break;
      }
    }
    if(!isThere) {
      // check list of NIST materials
      G4String mname = mat->GetName();
      for(G4int j=0; j<nvectors; ++j) {
        if(mname == nameNIST_ICRU90[j]) {
          materials[j] = mat;
          break;
    }
      }
    }
    isInitialized = (materials[0] && materials[1] && materials[2]);
    if(isInitialized) { return; }
  }
}

Referenced by G4ICRU90StoppingData(), G4BetheBlochModel::Initialise(), G4BraggIonModel::Initialise(), and G4BraggModel::Initialise().

IsApplicable()

G4bool G4ICRU90StoppingData::IsApplicable ( const G4Material *  mat ) const

inline

Definition at line 106 of file G4ICRU90StoppingData.hh.

{
  return (mat == materials[0] || mat == materials[1] || mat == materials[2]);
}

---

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

• G4ICRU90StoppingData.hh
• G4ICRU90StoppingData.cc