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=3; i<=ZPROJMAX; ++i) {
    fMatData[i] = new std::vector<G4PhysicsLogVector*>;
  }
}

~G4IonICRU73Data()

G4IonICRU73Data::~G4IonICRU73Data

(

)

Definition at line 111 of file G4IonICRU73Data.cc.

{
  delete fVector;
  for(G4int i=3; i<=ZPROJMAX; ++i) {
    auto v = fMatData[i];
    if(nullptr != v) {
      for(auto & dat : *v) {
    delete dat;
      }
      delete v;
    }
    for(G4int j=1; j<=ZTARGMAX; ++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 130 of file G4IonICRU73Data.cc.

{
  // if data does not exist the result is zero
  if (Z > ZPROJMAX) { return 0.0; }
 
  G4PhysicsLogVector* v = nullptr;
  if (1 == mat->GetNumberOfElements()) {
    G4int Z1 = (*(mat->GetElementVector()))[0]->GetZasInt();
    if(Z1 > ZTARGMAX) { return 0.0; }
    v = fElmData[Z][Z1];
  }
  else {
    G4int idx = fMatIndex[mat->GetIndex()];
    if (idx < 0) { return 0.0; }
    v = (*(fMatData[Z]))[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 155 of file G4IonICRU73Data.cc.

{
  // fill directory path
  if(fDataDirectory.empty()) {
    std::ostringstream ost;
    ost << G4EmParameters::Instance()->GetDirLEDATA() << "/ion_stopping_data/";
    fDataDirectory = ost.str();
  }
 
  const std::size_t nmat = G4Material::GetNumberOfMaterials();
 
  // do nothing if for a new run list of materials is not changed
  if(nmat == fMatIndex.size()) { return; }
 
  // look over all materials and recompute all materials
  // do not recompute element data if exist
  if(1 < fVerbose) {
    G4cout << "### G4IonICRU73Data::Initialise() for " << nmat
           << " materials" << G4endl;
  } 
  fMatIndex.resize(nmat, -1);
  for(G4int j=3; j<=ZPROJMAX; ++j) {
    fMatData[j]->resize(nmat, nullptr);
  }
  G4bool useICRU90 = G4EmParameters::Instance()->UseICRU90Data();
  auto mtable = G4Material::GetMaterialTable();
 
  for(G4int i=0; i<(G4int)nmat; ++i) {
    const G4Material* mat = (*mtable)[i];
    const G4String matname = mat->GetName(); 
    G4int idx = (G4int)mat->GetIndex();
    if(1 < fVerbose) {
      G4cout << i << ".  material: " << matname 
             << "  idx=" << idx << "  matIdx=" << fMatIndex[idx] 
             << G4endl;
    }
    if(fMatIndex[idx] == -1) {
      G4bool isOK = false; 
 
      // for material in ICRU90 
      if(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 " << matname << G4endl;
            }
            break;
          }
        }
      }
      // for material in ICRU73
      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 " << matname << G4endl;
            }
            break;
          }
        }
      }
      // dEdx is combined from element stopping power
      // Target element Z <= ZTARGMAX
      if(!isOK) {
    const std::size_t nelm = mat->GetNumberOfElements();
    auto elmv = mat->GetElementVector();
    G4bool isOut = false;
        for (std::size_t j=0; j<nelm; ++j) {
          if ((*elmv)[j]->GetZasInt() > ZTARGMAX) {
        isOut = true;
        break;
      }
    }
    if (!isOut) {
      ReadElementData(mat, useICRU90);
      isOK = true;
      if(1 < fVerbose) {
        G4cout << "Data via elements for " << matname << G4endl;
      }
    }
      }
      if(isOK) { fMatIndex[idx] = i; }
    }
    if(1 < fVerbose) {
      G4cout << "     matData: " << fMatData[i] << 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