G4LossTableBuilder Class Reference

#include <G4LossTableBuilder.hh>

Public Member Functions
	G4LossTableBuilder (G4bool master)
	~G4LossTableBuilder ()
void	BuildDEDXTable (G4PhysicsTable *dedxTable, const std::vector< G4PhysicsTable * > &)
void	BuildRangeTable (const G4PhysicsTable *dedxTable, G4PhysicsTable *rangeTable)
void	BuildInverseRangeTable (const G4PhysicsTable *rangeTable, G4PhysicsTable *invRangeTable)
G4PhysicsTable *	BuildTableForModel (G4PhysicsTable *table, G4VEmModel *model, const G4ParticleDefinition *, G4double emin, G4double emax, G4bool spline)
void	InitialiseBaseMaterials (const G4PhysicsTable *table=nullptr)
void	SetSplineFlag (G4bool flag)
void	SetBaseMaterialActive (G4bool flag)
G4LossTableBuilder &	operator= (const G4LossTableBuilder &right)=delete
	G4LossTableBuilder (const G4LossTableBuilder &)=delete

Static Public Member Functions
static const std::vector< G4int > *	GetCoupleIndexes ()
static const std::vector< G4double > *	GetDensityFactors ()
static G4bool	GetFlag (std::size_t idx)
static G4bool	GetBaseMaterialFlag ()

Detailed Description

Definition at line 59 of file G4LossTableBuilder.hh.

Constructor & Destructor Documentation

G4LossTableBuilder() [1/2]

G4LossTableBuilder::G4LossTableBuilder ( G4bool master )

explicit

Definition at line 76 of file G4LossTableBuilder.cc.

  : isInitializer(master)
{
  theParameters = G4EmParameters::Instance();
  if (nullptr == theFlag) {
    theDensityFactor = new std::vector<G4double>;
    theDensityIdx = new std::vector<G4int>;
    theFlag = new std::vector<G4bool>;
  }
}

Referenced by G4LossTableBuilder(), and operator=().

~G4LossTableBuilder()

G4LossTableBuilder::~G4LossTableBuilder

(

)

Definition at line 89 of file G4LossTableBuilder.cc.

{
  if (isInitializer) {
    delete theDensityFactor;
    delete theDensityIdx;
    delete theFlag;
    theDensityFactor = nullptr;
    theDensityIdx = nullptr;
    theFlag = nullptr;
  }
}

G4LossTableBuilder() [2/2]

G4LossTableBuilder::G4LossTableBuilder ( const G4LossTableBuilder & )

delete

Member Function Documentation

BuildDEDXTable()

void G4LossTableBuilder::BuildDEDXTable	(	G4PhysicsTable *	dedxTable,
		const std::vector< G4PhysicsTable * > &	list )

Definition at line 132 of file G4LossTableBuilder.cc.

{
  InitialiseBaseMaterials(dedxTable);
  std::size_t n_processes = list.size();
  if(1 >= n_processes) { return; }
 
  std::size_t nCouples = dedxTable->size();
  //G4cout << "Nproc= " << n_processes << " nCouples=" << nCouples << " Nv= " 
  //     << dedxTable->size() << G4endl;
  if(0 >= nCouples) { return; }
 
  for (std::size_t i=0; i<nCouples; ++i) {
    auto pv0 = static_cast<G4PhysicsLogVector*>((*(list[0]))[i]);
    //if (0 == i) G4cout << i << ". pv0=" << pv0 << "  t:" << list[0] << G4endl;
    if(pv0 == nullptr) { continue; } 
    std::size_t npoints = pv0->GetVectorLength();
    auto pv = new G4PhysicsLogVector(*pv0);
    for (std::size_t j=0; j<npoints; ++j) {
      G4double dedx = 0.0;
      for (std::size_t k=0; k<n_processes; ++k) {
    const G4PhysicsVector* pv1 = (*(list[k]))[i];
    //if (0 == i) G4cout << "     " << k << ". pv1=" << pv1 << "  t:" << list[k] << G4endl;
    dedx += (*pv1)[j];
      }
      pv->PutValue(j, dedx);
    }
    if(splineFlag) { pv->FillSecondDerivatives(); }
    G4PhysicsTableHelper::SetPhysicsVector(dedxTable, i, pv);
  }
  //G4cout << "### G4LossTableBuilder::BuildDEDXTable " << G4endl; 
  //G4cout << *dedxTable << G4endl;
}

BuildInverseRangeTable()

void G4LossTableBuilder::BuildInverseRangeTable	(	const G4PhysicsTable *	rangeTable,
		G4PhysicsTable *	invRangeTable )

Definition at line 250 of file G4LossTableBuilder.cc.

{
  std::size_t nCouples = rangeTable->size();
  if(0 >= nCouples) { return; }
 
  for (std::size_t i=0; i<nCouples; ++i) {
    G4PhysicsVector* pv = (*rangeTable)[i];
    if((pv == nullptr) || (isBaseMatActive && !(*theFlag)[i])) { continue; } 
    std::size_t npoints = pv->GetVectorLength();
      
    delete (*invRangeTable)[i];
    auto v = new G4PhysicsFreeVector(npoints, splineFlag);
 
    for (std::size_t j=0; j<npoints; ++j) {
      G4double e  = pv->Energy(j);
      G4double r  = (*pv)[j];
      v->PutValues(j,r,e);
    }
    if (splineFlag) { v->FillSecondDerivatives(); }
    v->EnableLogBinSearch(theParameters->NumberForFreeVector());
 
    G4PhysicsTableHelper::SetPhysicsVector(invRangeTable, i, v);
  }
  //G4cout << "### Inverse range table" << G4endl; 
  //G4cout << *invRangeTable << G4endl;
}

BuildRangeTable()

void G4LossTableBuilder::BuildRangeTable	(	const G4PhysicsTable *	dedxTable,
		G4PhysicsTable *	rangeTable )

Definition at line 168 of file G4LossTableBuilder.cc.

{
  //G4cout << "### G4LossTableBuilder::BuildRangeTable: DEDX table" << G4endl; 
  //G4cout << *const_cast<G4PhysicsTable*>(dedxTable) << G4endl;
  const std::size_t nCouples = dedxTable->size();
  if(0 >= nCouples) { return; }
 
  const std::size_t n = 100;
  const G4double del = 1.0/(G4double)n;
 
  for (std::size_t i=0; i<nCouples; ++i) {
    auto pv = static_cast<G4PhysicsLogVector*>((*dedxTable)[i]);
    if((pv == nullptr) || (isBaseMatActive && !(*theFlag)[i])) { continue; } 
    std::size_t npoints = pv->GetVectorLength();
    std::size_t bin0    = 0;
    G4double elow  = pv->Energy(0);
    G4double ehigh = pv->Energy(npoints-1);
    G4double dedx1 = (*pv)[0];
 
    // protection for specific cases dedx=0
    if(dedx1 == 0.0) {
      for (std::size_t k=1; k<npoints; ++k) {
    ++bin0;
    elow  = pv->Energy(k);
    dedx1 = (*pv)[k];
    if(dedx1 > 0.0) { break; }
      }
      npoints -= bin0;
    }
 
    // initialisation of a new vector
    if(npoints < 3) { npoints = 3; }
 
    delete (*rangeTable)[i];
    G4PhysicsLogVector* v;
    if(0 == bin0) { v = new G4PhysicsLogVector(*pv); }
    else { v = new G4PhysicsLogVector(elow, ehigh, npoints-1, splineFlag); }
 
    // assumed dedx proportional to beta
    G4double energy1 = v->Energy(0);
    G4double range   = 2.*energy1/dedx1;
    /*
    G4cout << "New Range vector Npoints=" << v->GetVectorLength()
       << " coupleIdx=" << i << " spline=" << v->GetSpline() 
       << " Elow=" << v->GetMinEnergy() <<" Ehigh=" << v->GetMinEnergy()
       << " DEDX(Elow)=" << dedx1 << " R(Elow)=" << range << G4endl;
    */
    v->PutValue(0,range);
 
    for (std::size_t j=1; j<npoints; ++j) {
 
      G4double energy2 = v->Energy(j);
      G4double de      = (energy2 - energy1) * del;
      G4double energy  = energy2 + de*0.5;
      G4double sum = 0.0;
      std::size_t idx = j - 1;
      for (std::size_t k=0; k<n; ++k) {
    energy -= de;
    dedx1 = pv->Value(energy, idx);
    if(dedx1 > 0.0) { sum += de/dedx1; }
      }
      range += sum;
      /*
      if(energy < 10.) 
    G4cout << "j= " << j << " e1= " << energy1 << " e2= " << energy2 
           << " n= " << n << " range=" << range<< G4endl;
      */
      v->PutValue(j,range);
      energy1 = energy2;
    }
    if(splineFlag) { v->FillSecondDerivatives(); }
    G4PhysicsTableHelper::SetPhysicsVector(rangeTable, i, v);
  }
  //G4cout << "### Range table" << G4endl; 
  //G4cout << *rangeTable << G4endl;
}

BuildTableForModel()

G4PhysicsTable * G4LossTableBuilder::BuildTableForModel	(	G4PhysicsTable *	table,
		G4VEmModel *	model,
		const G4ParticleDefinition *	part,
		G4double	emin,
		G4double	emax,
		G4bool	spline )

Definition at line 357 of file G4LossTableBuilder.cc.

{
  // check input
  G4PhysicsTable* table = G4PhysicsTableHelper::PreparePhysicsTable(aTable);
  if (nullptr == table) { return table; }
  if (aTable != nullptr && aTable != table) {
    aTable->clearAndDestroy();
    delete aTable;
  }
 
  InitialiseBaseMaterials(table);
  G4int nbins = theParameters->NumberOfBinsPerDecade();
 
  // Access to materials
  const G4ProductionCutsTable* theCoupleTable=
        G4ProductionCutsTable::GetProductionCutsTable();
  std::size_t numOfCouples = theCoupleTable->GetTableSize();
  /*
  G4cout << "   G4LossTableBuilder::BuildTableForModel Ncouple=" << numOfCouples
     << " isMaster=" << isInitializer << " model:" << model->GetName()
     << "  " << part->GetParticleName() << G4endl; 
  */
  G4PhysicsLogVector* aVector = nullptr;
 
  for(G4int i=0; i<(G4int)numOfCouples; ++i) {
    //G4cout << i << ".  " << (*theFlag)[i] << "  "  << table->GetFlag(i) << G4endl;
    if (table->GetFlag(i)) {
 
      // create physics vector and fill it
      auto couple = theCoupleTable->GetMaterialCutsCouple(i);
      delete (*table)[i];
 
      // if start from zero then change the scale
      const G4Material* mat = couple->GetMaterial();
 
      G4double tmin = std::max(emin, model->MinPrimaryEnergy(mat,part));
      if(0.0 >= tmin) { tmin = CLHEP::eV; }
      G4int n = nbins;
 
      if(tmin >= emax) {
        aVector = nullptr;
      } else {
        n *= G4lrint(std::log10(emax/tmin));
        n = std::max(n, 3);
        aVector = new G4PhysicsLogVector(tmin, emax, n, spline);
      }
 
      if(nullptr != aVector) {
        //G4cout << part->GetParticleName() << " in " << mat->GetName() 
    //     << " emin= " << tmin << " emax=" << emax << " n=" << n << G4endl;
        for(G4int j=0; j<=n; ++j) {
      G4double e = aVector->Energy(j);
      G4double y = model->Value(couple, part, e);
      //G4cout << "      " << j << ") E=" << e << "  y=" << y << G4endl;
          aVector->PutValue(j, y);
        }
        if(spline) { aVector->FillSecondDerivatives(); }
      }
      G4PhysicsTableHelper::SetPhysicsVector(table, i, aVector);
    }
  }
  /*
  G4cout << "G4LossTableBuilder::BuildTableForModel done for "
         << part->GetParticleName() << " and "<< model->GetName()
         << "  " << table << G4endl;
  */
  //G4cout << *table << G4endl;
  return table; 
}

Referenced by G4EmTableUtil::BuildMscProcess().

GetBaseMaterialFlag()

G4bool G4LossTableBuilder::GetBaseMaterialFlag ( )

static

Definition at line 124 of file G4LossTableBuilder.cc.

{
  return baseMatFlag;
}

Referenced by G4EmTableUtil::BuildMscProcess(), G4EmUtility::FillPeaksStructure(), G4GammaGeneralProcess::PreparePhysicsTable(), G4TransportationWithMsc::PreparePhysicsTable(), G4VEmProcess::PreparePhysicsTable(), and G4VEnergyLossProcess::PreparePhysicsTable().

GetCoupleIndexes()

const std::vector< G4int > * G4LossTableBuilder::GetCoupleIndexes ( )

static

Definition at line 103 of file G4LossTableBuilder.cc.

{
  return theDensityIdx;
}

Referenced by G4EmUtility::FillPeaksStructure(), G4VEmModel::G4VEmModel(), G4VEmProcess::G4VEmProcess(), and G4VEnergyLossProcess::G4VEnergyLossProcess().

GetDensityFactors()

const std::vector< G4double > * G4LossTableBuilder::GetDensityFactors ( )

static

Definition at line 110 of file G4LossTableBuilder.cc.

{
  return theDensityFactor;
}

Referenced by G4VEmModel::G4VEmModel(), G4VEmProcess::G4VEmProcess(), and G4VEnergyLossProcess::G4VEnergyLossProcess().

GetFlag()

G4bool G4LossTableBuilder::GetFlag ( std::size_t idx )

static

Definition at line 117 of file G4LossTableBuilder.cc.

{
  return (idx < theFlag->size()) ? (*theFlag)[idx] : false;
}

Referenced by G4EmTableUtil::BuildDEDXTable(), G4EmTableUtil::BuildLambdaTable(), G4EmTableUtil::BuildLambdaTable(), G4GammaGeneralProcess::BuildPhysicsTable(), and G4TransportationWithMsc::BuildPhysicsTable().

InitialiseBaseMaterials()

void G4LossTableBuilder::InitialiseBaseMaterials

(

const G4PhysicsTable *

table = nullptr

)

Definition at line 281 of file G4LossTableBuilder.cc.

{
  if(!isInitializer) { return; }
  const G4ProductionCutsTable* theCoupleTable=
    G4ProductionCutsTable::GetProductionCutsTable();
  std::size_t nCouples = theCoupleTable->GetTableSize();
  std::size_t nFlags = theFlag->size();
  /*
  G4cout << "### InitialiseBaseMaterials: nCouples=" << nCouples
     << " nFlags=" << nFlags << " isInit:" << isInitialized
     << " baseMat:" << baseMatFlag << G4endl; 
  */
  // define base material flag
  if(isBaseMatActive && !baseMatFlag) {
    for(G4int i=0; i<(G4int)nCouples; ++i) {
      if(nullptr != theCoupleTable->GetMaterialCutsCouple(i)->GetMaterial()->GetBaseMaterial()) {
    baseMatFlag = true;
        isInitialized = false;
    break;
      }
    }
  }
 
  if(nFlags != nCouples) { isInitialized = false; }
  if(isInitialized) { return; }
 
  // reserve memory
  theFlag->resize(nCouples, true);
  theDensityFactor->resize(nCouples,1.0);
  theDensityIdx->resize(nCouples, 0);
 
  // define default flag and index of used material cut couple
  for (G4int i=0; i<(G4int)nCouples; ++i) {
    (*theFlag)[i] = (nullptr == table) ? true : table->GetFlag(i); 
    (*theDensityIdx)[i] = i;
  }
  isInitialized = true;
  if (!baseMatFlag) { return; }
 
  // use base materials
  for (G4int i=0; i<(G4int)nCouples; ++i) { 
    // base material is needed only for a couple which is not
    // initialised and for which tables will be computed
    auto couple = theCoupleTable->GetMaterialCutsCouple(i);
    auto pcuts = couple->GetProductionCuts();
    auto mat  = couple->GetMaterial();
    auto bmat = mat->GetBaseMaterial();
 
    // base material exists - find it and check if it can be reused
    if(nullptr != bmat) {
      for(G4int j=0; j<(G4int)nCouples; ++j) {
    if(j == i) { continue; }
    auto bcouple = theCoupleTable->GetMaterialCutsCouple(j);
 
    if(bcouple->GetMaterial() == bmat && 
       bcouple->GetProductionCuts() == pcuts) {
 
      // based couple exist in the same region
      (*theDensityFactor)[i] = mat->GetDensity()/bmat->GetDensity();
      (*theDensityIdx)[i] = j;
      (*theFlag)[i] = false;
 
      // ensure that there will no double initialisation
      (*theDensityFactor)[j] = 1.0;
      (*theDensityIdx)[j] = j;
      (*theFlag)[j] = true;
      break; 
    }
      }
    }
  }
}

Referenced by BuildDEDXTable(), BuildTableForModel(), G4GammaGeneralProcess::PreparePhysicsTable(), G4TransportationWithMsc::PreparePhysicsTable(), G4VEmProcess::PreparePhysicsTable(), and G4VEnergyLossProcess::PreparePhysicsTable().

operator=()

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

delete

SetBaseMaterialActive()

void G4LossTableBuilder::SetBaseMaterialActive ( G4bool flag )

inline

Definition at line 131 of file G4LossTableBuilder.hh.

{
  isBaseMatActive = flag;
}

SetSplineFlag()

void G4LossTableBuilder::SetSplineFlag ( G4bool flag )

inline

Definition at line 121 of file G4LossTableBuilder.hh.

{
  splineFlag = flag;
}

---

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

• G4LossTableBuilder.hh
• G4LossTableBuilder.cc