#include <G4IonICRU73Data.hh>

Public Member Functions
	G4IonICRU73Data ()

	G4IonICRU73Data (const G4IonICRU73Data &)=delete

G4double	GetDEDX (const G4Material *, const G4int Z, const G4double e, const G4double loge) const

void	Initialise ()

G4IonICRU73Data &	operator= (const G4IonICRU73Data &right)=delete

	~G4IonICRU73Data ()

Private Member Functions
G4PhysicsLogVector *	FindOrBuildElementData (const G4int Z, const G4int Z1, G4bool useICRU90)

void	ReadElementData (const G4Material *mat, G4bool type)

void	ReadMaterialData (const G4String &name, G4bool type)

G4PhysicsLogVector *	RetrieveVector (std::ostringstream &in, G4bool warn)

Private Attributes
G4String	fDataDirectory = ""

G4PhysicsLogVector *	fElmData [81][81] = {{nullptr}}

G4double	fEmax

G4double	fEmin

std::vector< G4PhysicsLogVector * > *	fMatData [81] = {nullptr}

std::vector< G4int >	fMatIndex

G4int	fNbins = 0

G4int	fNbinsPerDecade = 10

G4int	fNmat = 0

G4bool	fSpline = false

G4PhysicsFreeVector *	fVector = nullptr

G4int	fVerbose = 0

Detailed Description

Definition at line 55 of file G4IonICRU73Data.hh.

Constructor & Destructor Documentation

◆ G4IonICRU73Data() [1/2]

G4IonICRU73Data::G4IonICRU73Data ( )

explicit

Definition at line 90 of file G4IonICRU73Data.cc.

{
  fEmin = 0.025*CLHEP::MeV;
  fEmax = 2.5*CLHEP::MeV;
  fNbins = fNbinsPerDecade*G4lrint(std::log10(fEmax/fEmin) + 0.000001);
  fVector = new G4PhysicsFreeVector(fSpline);
  for(G4int i=0; i<81; ++i) {
    fMatData[i] = new std::vector<G4PhysicsLogVector*>;
  }
}

References fEmax, fEmin, fMatData, fNbins, fNbinsPerDecade, fSpline, fVector, G4lrint(), and CLHEP::MeV.

◆ ~G4IonICRU73Data()

G4IonICRU73Data::~G4IonICRU73Data ( )

Definition at line 103 of file G4IonICRU73Data.cc.

{
  delete fVector;
  for(G4int i=0; i<81; ++i) {
    auto v = fMatData[i];
    for(G4int j=0; j<fNmat; ++j) {
      delete (*v)[j];
    }
    delete v;
    for(G4int j=0; j<81; ++j) { delete fElmData[i][j]; }
  }
}

References fElmData, fMatData, fNmat, and fVector.

◆ G4IonICRU73Data() [2/2]

G4IonICRU73Data::G4IonICRU73Data ( const G4IonICRU73Data & )

delete

Member Function Documentation

◆ FindOrBuildElementData()

G4PhysicsLogVector * G4IonICRU73Data::FindOrBuildElementData	(	const G4int	Z,
		const G4int	Z1,
		G4bool	useICRU90
	)

private

Definition at line 262 of file G4IonICRU73Data.cc.

{
  G4PhysicsLogVector* v = nullptr;
  if(Z <= 80 && Z1 <= 80) {
    v = fElmData[Z][Z1];
    if(nullptr == v) {
      G4bool isICRU90 = (useICRU90 && Z <= 18 && 
             (Z1 == 1 || Z1 == 6 || Z1 == 7 || Z1 == 8)); 
      std::ostringstream ost;
      ost << fDataDirectory << "icru";
      if(isICRU90) { ost << "90"; }
      else { ost << "73"; }
      ost << "/z" << Z << "_" << Z1 << ".dat"; 
      v = RetrieveVector(ost, false);
      fElmData[Z][Z1] = v;
    }
  }
  return v;
}

References fDataDirectory, fElmData, RetrieveVector(), Z, and anonymous_namespace{paraMaker.cc}::Z1.

Referenced by ReadElementData().

◆ GetDEDX()

G4double G4IonICRU73Data::GetDEDX	(	const G4Material *	mat,
		const G4int	Z,
		const G4double	e,
		const G4double	loge
	)		const

Definition at line 118 of file G4IonICRU73Data.cc.

{
  G4PhysicsLogVector* v = nullptr;
  if(1 == mat->GetNumberOfElements()) {
    G4int Z1 = (*(mat->GetElementVector()))[0]->GetZasInt();
    if(Z <= 80 && Z1 <= 80) {
      v = fElmData[Z][Z1];
    }
  } else {
    G4int idx = fMatIndex[mat->GetIndex()];
    if(idx < fNmat && Z <= 80) {
      v = (*(fMatData[Z]))[idx];
    }
  }
  G4double res = (nullptr != v) ? v->LogVectorValue(e, loge) : 0.0;
  return res;
}

References fElmData, fMatData, fMatIndex, fNmat, G4Material::GetElementVector(), G4Material::GetIndex(), G4Material::GetNumberOfElements(), G4PhysicsVector::LogVectorValue(), Z, and anonymous_namespace{paraMaker.cc}::Z1.

Referenced by G4LindhardSorensenIonModel::CorrectionsAlongStep().

◆ Initialise()

void G4IonICRU73Data::Initialise ( )

Definition at line 139 of file G4IonICRU73Data.cc.

{
  // fill directory path
  if(fDataDirectory.empty()) {
    char* path = std::getenv("G4LEDATA");
    if (nullptr != path) {
      std::ostringstream ost;
      ost << path << "/ion_stopping_data/";
      fDataDirectory = ost.str();
    } else {
      G4Exception("G4IonICRU73Data::Initialise(..)","em013",
                  FatalException,
                  "Environment variable G4LEDATA is not defined");
    }
  }
 
  G4bool useICRU90 = G4EmParameters::Instance()->UseICRU90Data();
  size_t nmat = G4Material::GetNumberOfMaterials();
  fMatIndex.resize(nmat, -1);
 
  const G4ProductionCutsTable* theCoupleTable=
        G4ProductionCutsTable::GetProductionCutsTable();
  size_t numOfCouples = theCoupleTable->GetTableSize();
 
  for(size_t i=0; i<numOfCouples; ++i) {
    const G4MaterialCutsCouple* couple = 
      theCoupleTable->GetMaterialCutsCouple(i);
    const G4Material* mat = couple->GetMaterial();
    G4int idx = mat->GetIndex();
    if(fMatIndex[idx] == -1) {
      G4String matname = mat->GetName();
      G4bool isOK = false; 
      if(1 == mat->GetNumberOfElements()) {
    ReadElementData(mat, useICRU90);
        isOK = true;
      }
      if(!isOK && useICRU90) {
    for(G4int j=0; j<3; ++j) {
      if(matname == namesICRU90[j]) {
        ReadMaterialData(matname, true);
            isOK = true;
        fMatIndex[idx] = fNmat; 
        ++fNmat;
        break;
      }
    }
      }
      if(!isOK) {
    for(G4int j=0; j<31; ++j) {
      if(matname == namesICRU73[j]) {
        ReadMaterialData(matname, false);
            isOK = true;
        fMatIndex[idx] = fNmat; 
        ++fNmat;
        break;
      }
    }
      }
      if(!isOK) {
    ReadElementData(mat, useICRU90);
    fMatIndex[idx] = fNmat; 
    ++fNmat;
      }
    }
  }
}

References FatalException, fDataDirectory, fMatIndex, fNmat, G4Exception(), G4Material::GetIndex(), G4MaterialCutsCouple::GetMaterial(), G4ProductionCutsTable::GetMaterialCutsCouple(), G4Material::GetName(), G4Material::GetNumberOfElements(), G4Material::GetNumberOfMaterials(), G4ProductionCutsTable::GetProductionCutsTable(), G4ProductionCutsTable::GetTableSize(), G4EmParameters::Instance(), namesICRU73, namesICRU90, ReadElementData(), ReadMaterialData(), and G4EmParameters::UseICRU90Data().

Referenced by G4LindhardSorensenIonModel::Initialise().

◆ operator=()

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

delete

◆ ReadElementData()

void G4IonICRU73Data::ReadElementData	(	const G4Material *	mat,
		G4bool	type
	)

private

Definition at line 223 of file G4IonICRU73Data.cc.

{
  const G4ElementVector* elmv = mat->GetElementVector();
  const G4double* dens = mat->GetFractionVector();
  const G4int nelm = mat->GetNumberOfElements();
  for(G4int Z=3; Z<81; ++Z) {
    if(1 == nelm) {
      FindOrBuildElementData(Z, (*elmv)[0]->GetZasInt(), useICRU90);
    } else {
      G4PhysicsLogVector* v2 = nullptr; 
      for(G4int j=0; j<nelm; ++j) {
    v2 = FindOrBuildElementData(Z, (*elmv)[j]->GetZasInt(), useICRU90);
        if(nullptr == v2) { break; }
      }
      // for one or more elements there is no data
      if(nullptr == v2) {
    fMatData[Z]->push_back(v2);
        continue;
      }
      G4PhysicsLogVector* v =
    new G4PhysicsLogVector(fEmin, fEmax, fNbins, fSpline);
      for(G4int i=0; i<=fNbins; ++i) {
    G4double dedx = 0;
        for(G4int j=0; j<nelm; ++j) {
      G4PhysicsLogVector* v1 = 
        FindOrBuildElementData(Z, (*elmv)[j]->GetZasInt(), useICRU90);
          dedx += (*v1)[i]*dens[j];
    }
        v->PutValue(i, dedx);
      }
      if(fSpline) { v->FillSecondDerivatives(); }
      fMatData[Z]->push_back(v);
    }
  }
}

References fEmax, fEmin, G4PhysicsVector::FillSecondDerivatives(), FindOrBuildElementData(), fMatData, fNbins, fSpline, G4Material::GetElementVector(), G4Material::GetFractionVector(), G4Material::GetNumberOfElements(), G4PhysicsVector::PutValue(), and Z.

Referenced by Initialise().

◆ ReadMaterialData()

void G4IonICRU73Data::ReadMaterialData	(	const G4String &	name,
		G4bool	type
	)

private

Definition at line 208 of file G4IonICRU73Data.cc.

{
  for(G4int Z=3; Z<81; ++Z) {
    std::ostringstream ost;
    ost << fDataDirectory << "icru";
    if(isICRU90 && Z <= 18) { ost << "90"; }
    else { ost << "73"; }
    ost << "/z" << Z << "_" << name << ".dat"; 
    G4PhysicsLogVector* v = RetrieveVector(ost, false);
    (fMatData[Z])->push_back(v);
  }
}

References fDataDirectory, fMatData, G4InuclParticleNames::name(), RetrieveVector(), and Z.

Referenced by Initialise().

◆ RetrieveVector()

G4PhysicsLogVector * G4IonICRU73Data::RetrieveVector	(	std::ostringstream &	in,
		G4bool	warn
	)

private

Definition at line 286 of file G4IonICRU73Data.cc.

{
  G4PhysicsLogVector* v = nullptr;
  std::ifstream filein(ost.str().c_str());
  if (!filein.is_open()) {
    if(warn) {
      G4ExceptionDescription ed;
      ed << "Data file <" << ost.str().c_str()
         << "> is not opened";
      G4Exception("G4IonICRU73Data::RetrieveVector(..)","em013",
                  FatalException, ed, "Check G4LEDATA");
    }
  } else {
    if(fVerbose > 0) {
      G4cout << "File " << ost.str() 
             << " is opened by G4IonICRU73Data" << G4endl;
    }
    // retrieve data from DB
    if(!fVector->Retrieve(filein, true)) {
      G4ExceptionDescription ed;
      ed << "Data file <" << ost.str().c_str()
         << "> is not retrieved!";
      G4Exception("G4IonICRU73Data::RetrieveVector(..)","had015",
                  FatalException, ed, "Check G4LEDATA");
    } else {
      if(fSpline) { fVector->FillSecondDerivatives(); }
      v = new G4PhysicsLogVector(fEmin, fEmax, fNbins, fSpline);
      for(G4int i=0; i<=fNbins; ++i) {
    G4double e = v->Energy(i);
    G4double dedx = fVector->Value(e);
    v->PutValue(i, dedx);
      }
      const G4double fact = CLHEP::MeV * CLHEP::cm2 /( 0.001 * CLHEP::g); 
      v->ScaleVector(CLHEP::MeV, fact);
      if(fSpline) { v->FillSecondDerivatives(); }
      if(fVerbose > 1) { G4cout << *v << G4endl; }
    }
  }
  return v;
}

References CLHEP::cm2, G4PhysicsVector::Energy(), FatalException, fEmax, fEmin, G4PhysicsVector::FillSecondDerivatives(), fNbins, fSpline, fVector, fVerbose, CLHEP::g, G4cout, G4endl, G4Exception(), CLHEP::MeV, G4PhysicsVector::PutValue(), G4PhysicsVector::Retrieve(), G4PhysicsVector::ScaleVector(), and G4PhysicsVector::Value().

Referenced by FindOrBuildElementData(), and ReadMaterialData().

Field Documentation

◆ fDataDirectory

G4String G4IonICRU73Data::fDataDirectory = ""

private

Definition at line 101 of file G4IonICRU73Data.hh.

Referenced by FindOrBuildElementData(), Initialise(), and ReadMaterialData().

◆ fElmData

G4PhysicsLogVector* G4IonICRU73Data::fElmData[81][81] = {{nullptr}}

private

Definition at line 93 of file G4IonICRU73Data.hh.

Referenced by FindOrBuildElementData(), GetDEDX(), and ~G4IonICRU73Data().

◆ fEmax

G4double G4IonICRU73Data::fEmax

private

Definition at line 89 of file G4IonICRU73Data.hh.

Referenced by G4IonICRU73Data(), ReadElementData(), and RetrieveVector().

◆ fEmin

G4double G4IonICRU73Data::fEmin

private

Definition at line 88 of file G4IonICRU73Data.hh.

Referenced by G4IonICRU73Data(), ReadElementData(), and RetrieveVector().

◆ fMatData

std::vector<G4PhysicsLogVector*>* G4IonICRU73Data::fMatData[81] = {nullptr}

private

Definition at line 92 of file G4IonICRU73Data.hh.

Referenced by G4IonICRU73Data(), GetDEDX(), ReadElementData(), ReadMaterialData(), and ~G4IonICRU73Data().

◆ fMatIndex

std::vector<G4int> G4IonICRU73Data::fMatIndex

private

Definition at line 91 of file G4IonICRU73Data.hh.

Referenced by GetDEDX(), and Initialise().

◆ fNbins

G4int G4IonICRU73Data::fNbins = 0

private

Definition at line 97 of file G4IonICRU73Data.hh.

Referenced by G4IonICRU73Data(), ReadElementData(), and RetrieveVector().

◆ fNbinsPerDecade

G4int G4IonICRU73Data::fNbinsPerDecade = 10

private

Definition at line 98 of file G4IonICRU73Data.hh.

Referenced by G4IonICRU73Data().

◆ fNmat

G4int G4IonICRU73Data::fNmat = 0

private

Definition at line 96 of file G4IonICRU73Data.hh.

Referenced by GetDEDX(), Initialise(), and ~G4IonICRU73Data().

◆ fSpline

G4bool G4IonICRU73Data::fSpline = false

private

Definition at line 100 of file G4IonICRU73Data.hh.

Referenced by G4IonICRU73Data(), ReadElementData(), and RetrieveVector().

◆ fVector

G4PhysicsFreeVector* G4IonICRU73Data::fVector = nullptr

private

Definition at line 94 of file G4IonICRU73Data.hh.

Referenced by G4IonICRU73Data(), RetrieveVector(), and ~G4IonICRU73Data().

◆ fVerbose

G4int G4IonICRU73Data::fVerbose = 0

private

Definition at line 99 of file G4IonICRU73Data.hh.

Referenced by RetrieveVector().

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

source/processes/electromagnetic/standard/include/G4IonICRU73Data.hh
source/processes/electromagnetic/standard/src/G4IonICRU73Data.cc

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4IonICRU73Data() [1/2]

◆ ~G4IonICRU73Data()

◆ G4IonICRU73Data() [2/2]

Member Function Documentation

◆ FindOrBuildElementData()

◆ GetDEDX()

◆ Initialise()

◆ operator=()

◆ ReadElementData()

◆ ReadMaterialData()

◆ RetrieveVector()

Field Documentation

◆ fDataDirectory

◆ fElmData

◆ fEmax

◆ fEmin

◆ fMatData

◆ fMatIndex

◆ fNbins

◆ fNbinsPerDecade

◆ fNmat

◆ fSpline

◆ fVector

◆ fVerbose