G4IonICRU73Data Class Reference

#include <G4IonICRU73Data.hh>

Public Member Functions
	G4IonICRU73Data ()
	~G4IonICRU73Data ()
void	Initialise ()
G4double	GetDEDX (const G4Material *, const G4int Z, const G4double e, const G4double loge) const
G4IonICRU73Data &	operator= (const G4IonICRU73Data &right)=delete
	G4IonICRU73Data (const G4IonICRU73Data &)=delete

Detailed Description

Definition at line 55 of file G4IonICRU73Data.hh.

Constructor & Destructor Documentation

G4IonICRU73Data() [1/2]

G4IonICRU73Data::G4IonICRU73Data ( )

explicit

Definition at line 98 of file G4IonICRU73Data.cc.

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

~G4IonICRU73Data()

G4IonICRU73Data::~G4IonICRU73Data

(

)

Definition at line 111 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<93; ++j) { delete fElmData[i][j]; }
  }
}

G4IonICRU73Data() [2/2]

G4IonICRU73Data::G4IonICRU73Data ( const G4IonICRU73Data &  )

delete

Member Function Documentation

GetDEDX()

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

Definition at line 126 of file G4IonICRU73Data.cc.

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

Referenced by G4LindhardSorensenIonModel::CorrectionsAlongStep().

Initialise()

void G4IonICRU73Data::Initialise

(

)

Definition at line 148 of file G4IonICRU73Data.cc.

{
  // fill directory path
  if(fDataDirectory.empty()) {
    const char* path = G4FindDataDir("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");
    }
  }
 
  std::size_t nmat = G4Material::GetNumberOfMaterials();
  if(nmat == fMatIndex.size()) { return; }
 
  if(0 < fVerbose) {
    G4cout << "### G4IonICRU73Data::Initialise() for " << nmat
           << " materials" << G4endl;
  } 
  fMatIndex.resize(nmat, -1);
  for(G4int j=0; j<81; ++j) {
    fMatData[j]->resize(nmat, nullptr);
  }
  G4bool useICRU90 = G4EmParameters::Instance()->UseICRU90Data();
  auto mtable = G4Material::GetMaterialTable();
 
  for(G4int i=0; i<(G4int)nmat; ++i) {
    fNmat = i;
    const G4Material* mat = (*mtable)[i];
    G4int idx = (G4int)mat->GetIndex();
    if(1 < fVerbose) {
      G4cout << i << ".  material:" << mat->GetName() 
             << "  idx=" << idx << G4endl;
    }
    if(fMatIndex[idx] == -1) {
      fMatIndex[idx] = i;
      G4String matname = mat->GetName();
      G4bool isOK = false; 
      if(1 == mat->GetNumberOfElements()) {
        ReadElementData(mat, useICRU90);
        isOK = true;
        if(1 < fVerbose) {
          G4cout << "Material from single element" << G4endl;
        }
      }
      if(!isOK && useICRU90) {
        for(G4int j=0; j<3; ++j) {
          if(matname == namesICRU90[j]) {
            ReadMaterialData(mat, densityCoef[j], true);
            isOK = true;
            if(1 < fVerbose) {
              G4cout << "ICRU90 material" << G4endl;
            }
            break;
          }
        }
      }
      if(!isOK) {
        for(G4int j=0; j<31; ++j) {
          if(matname == namesICRU73[j]) {
            ReadMaterialData(mat, 1.0, false);
            isOK = true;
            if(1 < fVerbose) {
              G4cout << "ICRU73 material" << G4endl;
            }
            break;
          }
        }
      }
      if(!isOK) {
        ReadElementData(mat, useICRU90);
        if(1 < fVerbose) {
          G4cout << "Read element data" << G4endl;
        }
      }
    }
  }
}

Referenced by G4LindhardSorensenIonModel::Initialise().

operator=()

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

delete

---

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

• G4IonICRU73Data.hh
• G4IonICRU73Data.cc