G4PhysicsFreeVector Class Reference

#include <G4PhysicsFreeVector.hh>

Inheritance diagram for G4PhysicsFreeVector:

Public Member Functions
	G4PhysicsFreeVector (G4bool spline=false)
	G4PhysicsFreeVector (G4int length)
	G4PhysicsFreeVector (std::size_t length, G4bool spline=false)
	G4PhysicsFreeVector (std::size_t length, G4double emin, G4double emax, G4bool spline=false)
	G4PhysicsFreeVector (const std::vector< G4double > &energies, const std::vector< G4double > &values, G4bool spline=false)
	G4PhysicsFreeVector (const G4double *energies, const G4double *values, std::size_t length, G4bool spline=false)
	~G4PhysicsFreeVector () override=default
void	PutValues (const std::size_t index, const G4double energy, const G4double value)
void	InsertValues (const G4double energy, const G4double value)
void	EnableLogBinSearch (const G4int n=1)
void	PutValue (const std::size_t index, const G4double e, const G4double value)
Public Member Functions inherited from G4PhysicsVector
	G4PhysicsVector (G4bool spline=false)
	G4PhysicsVector (const G4PhysicsVector &)=default
G4PhysicsVector &	operator= (const G4PhysicsVector &)=default
	G4PhysicsVector (const G4PhysicsVector &&)=delete
G4PhysicsVector &	operator= (const G4PhysicsVector &&)=delete
G4bool	operator== (const G4PhysicsVector &right) const =delete
G4bool	operator!= (const G4PhysicsVector &right) const =delete
virtual	~G4PhysicsVector ()=default
G4double	Value (const G4double energy, std::size_t &lastidx) const
G4double	Value (const G4double energy) const
G4double	GetValue (const G4double energy, G4bool &isOutRange) const
G4double	LogVectorValue (const G4double energy, const G4double theLogEnergy) const
G4double	LogFreeVectorValue (const G4double energy, const G4double theLogEnergy) const
G4double	operator[] (const std::size_t index) const
G4double	operator() (const std::size_t index) const
void	PutValue (const std::size_t index, const G4double value)
G4double	Energy (const std::size_t index) const
G4double	GetLowEdgeEnergy (const std::size_t index) const
G4double	GetMinEnergy () const
G4double	GetMaxEnergy () const
G4double	GetMinValue () const
G4double	GetMaxValue () const
std::size_t	GetVectorLength () const
std::size_t	ComputeLogVectorBin (const G4double logenergy) const
G4PhysicsVectorType	GetType () const
G4bool	GetSpline () const
void	SetVerboseLevel (G4int value)
G4double	FindLinearEnergy (const G4double rand) const
std::size_t	FindBin (const G4double energy, std::size_t idx) const
void	ScaleVector (const G4double factorE, const G4double factorV)
void	FillSecondDerivatives (const G4SplineType=G4SplineType::Base, const G4double dir1=0.0, const G4double dir2=0.0)
G4double	GetEnergy (const G4double value) const
G4bool	Store (std::ofstream &fOut, G4bool ascii=false) const
G4bool	Retrieve (std::ifstream &fIn, G4bool ascii=false)
void	DumpValues (G4double unitE=1.0, G4double unitV=1.0) const

Additional Inherited Members
Protected Member Functions inherited from G4PhysicsVector
virtual void	Initialise ()
void	PrintPutValueError (std::size_t index, G4double value, const G4String &text)
Protected Attributes inherited from G4PhysicsVector
G4double	edgeMin = 0.0
G4double	edgeMax = 0.0
G4double	invdBin = 0.0
G4double	logemin = 0.0
G4double	iBin1 = 0.0
G4double	lmin1 = 0.0
G4int	verboseLevel = 0
std::size_t	idxmax = 0
std::size_t	imax1 = 0
std::size_t	numberOfNodes = 0
std::size_t	nLogNodes = 0
G4PhysicsVectorType	type = T_G4PhysicsFreeVector
std::vector< G4double >	binVector
std::vector< G4double >	dataVector
std::vector< G4double >	secDerivative
std::vector< std::size_t >	scale

Detailed Description

Definition at line 53 of file G4PhysicsFreeVector.hh.

Constructor & Destructor Documentation

G4PhysicsFreeVector() [1/6]

G4PhysicsFreeVector::G4PhysicsFreeVector ( G4bool spline = false )

explicit

Definition at line 40 of file G4PhysicsFreeVector.cc.

  : G4PhysicsVector(spline)
{}

G4PhysicsFreeVector() [2/6]

G4PhysicsFreeVector::G4PhysicsFreeVector ( G4int length )

explicit

Definition at line 45 of file G4PhysicsFreeVector.cc.

  : G4PhysicsFreeVector(static_cast<std::size_t>(length), false)
{}

G4PhysicsFreeVector() [3/6]

G4PhysicsFreeVector::G4PhysicsFreeVector ( std::size_t length, G4bool spline = false )

explicit

Definition at line 50 of file G4PhysicsFreeVector.cc.

  : G4PhysicsVector(spline)
{
  numberOfNodes = length;
 
  if (0 < length) {
    binVector.resize(numberOfNodes, 0.0);
    dataVector.resize(numberOfNodes, 0.0);
  }
  Initialise();
}

G4PhysicsFreeVector() [4/6]

G4PhysicsFreeVector::G4PhysicsFreeVector ( std::size_t length, G4double emin, G4double emax, G4bool spline = false )

explicit

Definition at line 63 of file G4PhysicsFreeVector.cc.

  : G4PhysicsFreeVector(length, spline)
{}

G4PhysicsFreeVector() [5/6]

G4PhysicsFreeVector::G4PhysicsFreeVector ( const std::vector< G4double > & energies, const std::vector< G4double > & values, G4bool spline = false )

explicit

Definition at line 69 of file G4PhysicsFreeVector.cc.

  : G4PhysicsVector(spline)
{
  numberOfNodes = energies.size();
 
  if (numberOfNodes != values.size())
  {
    G4ExceptionDescription ed;
    ed << "The size of energy vector " << numberOfNodes << " != " << values.size();
    G4Exception("G4PhysicsFreeVector constructor: ","glob04", FatalException, ed);
  }
 
  binVector = energies;
  dataVector = values;
  Initialise();
}

G4PhysicsFreeVector() [6/6]

G4PhysicsFreeVector::G4PhysicsFreeVector ( const G4double * energies, const G4double * values, std::size_t length, G4bool spline = false )

explicit

Definition at line 89 of file G4PhysicsFreeVector.cc.

  : G4PhysicsVector(spline)
{
  numberOfNodes = length;
 
  if (0 < numberOfNodes) 
  {
    binVector.resize(numberOfNodes);
    dataVector.resize(numberOfNodes);
 
    for(std::size_t i = 0; i < numberOfNodes; ++i)
    {
      binVector[i] = energies[i];
      dataVector[i] = values[i];
    }
  }
  Initialise();
}

~G4PhysicsFreeVector()

G4PhysicsFreeVector::~G4PhysicsFreeVector ( )

overridedefault

Member Function Documentation

EnableLogBinSearch()

void G4PhysicsFreeVector::EnableLogBinSearch

(

const G4int

n = 1

)

Definition at line 149 of file G4PhysicsFreeVector.cc.

{
  // check if log search is applicable
  if (0 >= n || edgeMin <= 0.0 || edgeMin == edgeMax || numberOfNodes < 3)
  {
    return;
  }
  nLogNodes = static_cast<std::size_t>(static_cast<G4int>(numberOfNodes)/n);
  if (nLogNodes < 3) { nLogNodes = 3; }
  scale.resize(nLogNodes, 0);
  imax1 = nLogNodes - 2;
  iBin1 = (imax1 + 1) / G4Log(edgeMax/edgeMin);
  lmin1 = G4Log(edgeMin);
  scale[0] = 0;
  scale[imax1 + 1] = idxmax;
  std::size_t j = 0;
  for (std::size_t i = 1; i <= imax1; ++i)
  {
    G4double e = edgeMin*G4Exp(i/iBin1);
    for (; j <= idxmax; ++j)
    {
      if (binVector[j] <= e && e < binVector[j+1])
      {
        scale[i] = j;
        break;
      }
    }
  }
}

InsertValues()

void G4PhysicsFreeVector::InsertValues	(	const G4double	energy,
		const G4double	value )

Definition at line 134 of file G4PhysicsFreeVector.cc.

{
  auto binLoc = std::lower_bound(binVector.cbegin(), binVector.cend(), energy);
  auto dataLoc = dataVector.cbegin();
  dataLoc += binLoc - binVector.cbegin(); 
 
  binVector.insert(binLoc, energy);
  dataVector.insert(dataLoc, value);
 
  ++numberOfNodes;
  Initialise();
}

Referenced by G4MaterialPropertiesTable::AddEntry(), G4SPSEneDistribution::ArbEnergyHisto(), G4SPSEneDistribution::ArbEnergyHistoFile(), G4Cerenkov::BuildPhysicsTable(), G4AdjointPrimaryGenerator::ComputeAccumulatedDepthVectorAlongBackRay(), G4SPSEneDistribution::EpnEnergyHisto(), G4SPSRandomGenerator::GenRandEnergy(), G4SPSRandomGenerator::GenRandPhi(), G4SPSRandomGenerator::GenRandPosPhi(), G4SPSRandomGenerator::GenRandPosTheta(), G4SPSRandomGenerator::GenRandTheta(), G4SPSRandomGenerator::GenRandX(), G4SPSRandomGenerator::GenRandY(), G4SPSRandomGenerator::GenRandZ(), G4GDMLReadMaterials::PropertyRead(), G4GDMLReadSolids::PropertyRead(), G4SPSRandomGenerator::SetEnergyBias(), G4SPSRandomGenerator::SetPhiBias(), G4SPSRandomGenerator::SetPosPhiBias(), G4SPSRandomGenerator::SetPosThetaBias(), G4SPSRandomGenerator::SetThetaBias(), G4SPSRandomGenerator::SetXBias(), G4SPSRandomGenerator::SetYBias(), G4SPSRandomGenerator::SetZBias(), G4SPSAngDistribution::UserDefAngPhi(), G4SPSAngDistribution::UserDefAngTheta(), and G4SPSEneDistribution::UserEnergyHisto().

PutValue()

void G4PhysicsFreeVector::PutValue ( const std::size_t index, const G4double e, const G4double value )

inline

Definition at line 97 of file G4PhysicsFreeVector.hh.

{
  PutValues(index, e, value);
}

Referenced by G4PartialWidthTable::AddWidths(), G4DiffuseElastic::BuildAngleTable(), G4NuclNuclDiffuseElastic::BuildAngleTable(), G4hhElastic::BuildTableT(), G4hhElastic::BuildTableTest(), G4XDeltaDeltastarTable::CrossSectionTable(), G4XDeltaDeltaTable::CrossSectionTable(), G4XDeltaNstarTable::CrossSectionTable(), G4XNDeltastarTable::CrossSectionTable(), G4XNDeltaTable::CrossSectionTable(), G4XNNstarTable::CrossSectionTable(), G4EMDissociationCrossSection::GetCrossSectionForProjectile(), G4NeutronElectronElModel::Initialise(), G4BaryonPartialWidth::MassDependentWidth(), G4BaryonWidth::MassDependentWidth(), G4DiffuseElastic::TestAngleTable(), and G4NuclNuclDiffuseElastic::TestAngleTable().

PutValues()

void G4PhysicsFreeVector::PutValues	(	const std::size_t	index,
		const G4double	energy,
		const G4double	value )

Definition at line 112 of file G4PhysicsFreeVector.cc.

{
  if (index >= numberOfNodes)
  {
    PrintPutValueError(index, value, "G4PhysicsFreeVector::PutValues ");
    return;
  }
  binVector[index]  = e;
  dataVector[index] = value;
  if(index == 0)
  {
    edgeMin = e;
  }
  else if(numberOfNodes == index + 1)
  {
    edgeMax = e;
  }
}

Referenced by G4PenelopeCrossSection::AddCrossSectionPoint(), G4PenelopeCrossSection::AddShellCrossSectionPoint(), G4PenelopeBremsstrahlungFS::BuildScaledXSTable(), G4PenelopeCrossSection::NormalizeShellCrossSections(), G4PenelopeBremsstrahlungAngular::PrepareTables(), PutValue(), and G4PenelopeBremsstrahlungFS::SampleGammaEnergy().

---

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

• G4PhysicsFreeVector.hh
• G4PhysicsFreeVector.cc