G4PhysicsFreeVector Class Reference

#include <G4PhysicsFreeVector.hh>

Inheritance diagram for G4PhysicsFreeVector:

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

Protected Member Functions
virtual void	Initialise ()
void	PrintPutValueError (std::size_t index, G4double value, const G4String &text)

Protected Attributes
std::vector< G4double >	binVector
std::vector< G4double >	dataVector
G4double	edgeMax = 0.0
G4double	edgeMin = 0.0
G4int	idxmax = 0
G4double	invdBin = 0.0
G4double	logemin = 0.0
std::size_t	numberOfNodes = 0
std::vector< G4double >	secDerivative
G4PhysicsVectorType	type = T_G4PhysicsFreeVector
G4int	verboseLevel = 0

Private Member Functions
void	ComputeSecDerivative0 ()
void	ComputeSecDerivative1 ()
void	ComputeSecDerivative2 (const G4double firstPointDerivative, const G4double endPointDerivative)
std::size_t	GetBin (const G4double energy) const
G4double	Interpolation (const std::size_t idx, const G4double energy) const

Private Attributes
G4bool	useSpline = false

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 39 of file G4PhysicsFreeVector.cc.

  : G4PhysicsVector(spline)
{}

G4PhysicsFreeVector() [2/6]

G4PhysicsFreeVector::G4PhysicsFreeVector ( G4int  length )

explicit

Definition at line 44 of file G4PhysicsFreeVector.cc.

  : G4PhysicsFreeVector((std::size_t)length, false)
{}

G4PhysicsFreeVector() [3/6]

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

explicit

Definition at line 49 of file G4PhysicsFreeVector.cc.

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

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, G4PhysicsVector::Initialise(), and G4PhysicsVector::numberOfNodes.

G4PhysicsFreeVector() [4/6]

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

explicit

Definition at line 62 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 68 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();
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, FatalException, G4Exception(), G4PhysicsVector::Initialise(), and G4PhysicsVector::numberOfNodes.

G4PhysicsFreeVector() [6/6]

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

explicit

Definition at line 88 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();
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, G4PhysicsVector::Initialise(), and G4PhysicsVector::numberOfNodes.

~G4PhysicsFreeVector()

virtual G4PhysicsFreeVector::~G4PhysicsFreeVector ( )

virtualdefault

Member Function Documentation

ComputeLogVectorBin()

std::size_t G4PhysicsVector::ComputeLogVectorBin ( const G4double  logenergy ) const

inlineinherited

ComputeSecDerivative0()

void G4PhysicsVector::ComputeSecDerivative0 ( )

privateinherited

Definition at line 272 of file G4PhysicsVector.cc.

{
  std::size_t n = numberOfNodes - 1;
 
  for(std::size_t i = 1; i < n; ++i)
  {
    secDerivative[i] = 3.0 *
      ((dataVector[i + 1] - dataVector[i]) / (binVector[i + 1] - binVector[i]) -
       (dataVector[i] - dataVector[i - 1]) /
         (binVector[i] - binVector[i - 1])) /
      (binVector[i + 1] - binVector[i - 1]);
  }
  secDerivative[n] = secDerivative[n - 1];
  secDerivative[0] = secDerivative[1];
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, CLHEP::detail::n, G4PhysicsVector::numberOfNodes, and G4PhysicsVector::secDerivative.

Referenced by G4PhysicsVector::FillSecondDerivatives().

ComputeSecDerivative1()

void G4PhysicsVector::ComputeSecDerivative1 ( )

privateinherited

Definition at line 290 of file G4PhysicsVector.cc.

{
  G4int n = numberOfNodes - 1;
  G4double* u = new G4double[n];
  G4double p, sig;
 
  u[1] = ((dataVector[2] - dataVector[1]) / (binVector[2] - binVector[1]) -
          (dataVector[1] - dataVector[0]) / (binVector[1] - binVector[0]));
  u[1] = 6.0 * u[1] * (binVector[2] - binVector[1]) /
         ((binVector[2] - binVector[0]) * (binVector[2] - binVector[0]));
 
  // Decomposition loop for tridiagonal algorithm. secDerivative[i]
  // and u[i] are used for temporary storage of the decomposed factors.
 
  secDerivative[1] = (2.0 * binVector[1] - binVector[0] - binVector[2]) /
                     (2.0 * binVector[2] - binVector[0] - binVector[1]);
 
  for(G4int i = 2; i < n - 1; ++i)
  {
    sig =
      (binVector[i] - binVector[i - 1]) / (binVector[i + 1] - binVector[i - 1]);
    p                = sig * secDerivative[i - 1] + 2.0;
    secDerivative[i] = (sig - 1.0) / p;
    u[i] =
      (dataVector[i + 1] - dataVector[i]) / (binVector[i + 1] - binVector[i]) -
      (dataVector[i] - dataVector[i - 1]) / (binVector[i] - binVector[i - 1]);
    u[i] =
      (6.0 * u[i] / (binVector[i + 1] - binVector[i - 1])) - sig * u[i - 1] / p;
  }
 
  sig =
    (binVector[n - 1] - binVector[n - 2]) / (binVector[n] - binVector[n - 2]);
  p = sig * secDerivative[n - 3] + 2.0;
  u[n - 1] =
    (dataVector[n] - dataVector[n - 1]) / (binVector[n] - binVector[n - 1]) -
    (dataVector[n - 1] - dataVector[n - 2]) /
      (binVector[n - 1] - binVector[n - 2]);
  u[n - 1] = 6.0 * sig * u[n - 1] / (binVector[n] - binVector[n - 2]) -
             (2.0 * sig - 1.0) * u[n - 2] / p;
 
  p = (1.0 + sig) + (2.0 * sig - 1.0) * secDerivative[n - 2];
  secDerivative[n - 1] = u[n - 1] / p;
 
  // The back-substitution loop for the triagonal algorithm of solving
  // a linear system of equations.
 
  for(G4int k = n - 2; k > 1; --k)
  {
    secDerivative[k] *=
      (secDerivative[k + 1] - u[k] * (binVector[k + 1] - binVector[k - 1]) /
                                (binVector[k + 1] - binVector[k]));
  }
  secDerivative[n] =
    (secDerivative[n - 1] - (1.0 - sig) * secDerivative[n - 2]) / sig;
  sig = 1.0 - ((binVector[2] - binVector[1]) / (binVector[2] - binVector[0]));
  secDerivative[1] *= (secDerivative[2] - u[1] / (1.0 - sig));
  secDerivative[0] = (secDerivative[1] - sig * secDerivative[2]) / (1.0 - sig);
 
  delete[] u;
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, CLHEP::detail::n, G4PhysicsVector::numberOfNodes, and G4PhysicsVector::secDerivative.

Referenced by G4PhysicsVector::FillSecondDerivatives().

ComputeSecDerivative2()

void G4PhysicsVector::ComputeSecDerivative2 ( const G4double  firstPointDerivative, const G4double  endPointDerivative  )

privateinherited

Definition at line 355 of file G4PhysicsVector.cc.

{
  G4int n = numberOfNodes - 1;
  G4double* u = new G4double[n];
  G4double p, sig, un;
 
  u[0] = (6.0 / (binVector[1] - binVector[0])) *
         ((dataVector[1] - dataVector[0]) / (binVector[1] - binVector[0]) -
          firstPointDerivative);
 
  secDerivative[0] = -0.5;
 
  // Decomposition loop for tridiagonal algorithm. secDerivative[i]
  // and u[i] are used for temporary storage of the decomposed factors.
 
  for(G4int i = 1; i < n; ++i)
  {
    sig =
      (binVector[i] - binVector[i - 1]) / (binVector[i + 1] - binVector[i - 1]);
    p                = sig * (secDerivative[i - 1]) + 2.0;
    secDerivative[i] = (sig - 1.0) / p;
    u[i] =
      (dataVector[i + 1] - dataVector[i]) / (binVector[i + 1] - binVector[i]) -
      (dataVector[i] - dataVector[i - 1]) / (binVector[i] - binVector[i - 1]);
    u[i] =
      6.0 * u[i] / (binVector[i + 1] - binVector[i - 1]) - sig * u[i - 1] / p;
  }
 
  sig =
    (binVector[n - 1] - binVector[n - 2]) / (binVector[n] - binVector[n - 2]);
  p  = sig * secDerivative[n - 2] + 2.0;
  un = (6.0 / (binVector[n] - binVector[n - 1])) *
         (endPointDerivative - (dataVector[n] - dataVector[n - 1]) /
                                 (binVector[n] - binVector[n - 1])) -
       u[n - 1] / p;
  secDerivative[n] = un / (secDerivative[n - 1] + 2.0);
 
  // The back-substitution loop for the triagonal algorithm of solving
  // a linear system of equations.
 
  for(G4int k = n - 1; k > 0; --k)
  {
    secDerivative[k] *=
      (secDerivative[k + 1] - u[k] * (binVector[k + 1] - binVector[k - 1]) /
                                (binVector[k + 1] - binVector[k]));
  }
  secDerivative[0] = 0.5 * (u[0] - secDerivative[1]);
 
  delete[] u;
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, CLHEP::detail::n, G4PhysicsVector::numberOfNodes, and G4PhysicsVector::secDerivative.

Referenced by G4PhysicsVector::FillSecondDerivatives().

DumpValues()

void G4PhysicsVector::DumpValues ( G4double  unitE = 1.0, G4double  unitV = 1.0  ) const

inherited

Definition at line 163 of file G4PhysicsVector.cc.

{
  for(std::size_t i = 0; i < numberOfNodes; ++i)
  {
    G4cout << binVector[i] / unitE << "   " << dataVector[i] / unitV 
           << G4endl;
  }
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, G4cout, G4endl, and G4PhysicsVector::numberOfNodes.

Referenced by G4OpWLS::DumpPhysicsTable(), G4OpWLS2::DumpPhysicsTable(), G4PhysicsVector::FillSecondDerivatives(), G4SPSEneDistribution::GenUserHistEnergies(), and G4SPSEneDistribution::LinearInterpolation().

Energy()

G4double G4PhysicsVector::Energy ( const std::size_t  index ) const

inlineinherited

Referenced by G4SPSEneDistribution::ArbInterpolate(), G4EmCorrections::BarkasCorrection(), G4DiffuseElasticV2::BuildAngleTable(), G4PolarizedCompton::BuildAsymmetryTable(), G4PolarizedAnnihilation::BuildAsymmetryTables(), G4PolarizedIonisation::BuildAsymmetryTables(), G4EmCorrections::BuildCorrectionVector(), G4KokoulinMuonNuclearXS::BuildCrossSectionTable(), G4LossTableBuilder::BuildInverseRangeTable(), G4VLEPTSModel::BuildMeanFreePathTable(), G4Cerenkov::BuildPhysicsTable(), G4Scintillation::BuildPhysicsTable(), G4OpWLS::BuildPhysicsTable(), G4OpWLS2::BuildPhysicsTable(), G4LossTableBuilder::BuildRangeTable(), G4PenelopeBremsstrahlungFS::BuildScaledXSTable(), G4LossTableBuilder::BuildTableForModel(), G4AdjointCSManager::BuildTotalSigmaTables(), G4MaterialPropertiesTable::CalculateGROUPVEL(), G4OpRayleigh::CalculateRayleighMeanFreePaths(), G4PenelopeRayleighModelMI::CalculateThetaAndAngFun(), G4eeToHadronsModel::ComputeCMCrossSectionPerElectron(), G4JAEAElasticScatteringModel::ComputeCrossSectionPerAtom(), G4JAEAPolarizedElasticScatteringModel::ComputeCrossSectionPerAtom(), G4LivermoreComptonModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedComptonModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedRayleighModel::ComputeCrossSectionPerAtom(), G4LivermoreRayleighModel::ComputeCrossSectionPerAtom(), G4LowEPComptonModel::ComputeCrossSectionPerAtom(), G4LowEPPolarizedComptonModel::ComputeCrossSectionPerAtom(), G4TablesForExtrapolator::ComputeElectronDEDX(), G4TablesForExtrapolator::ComputeMuonDEDX(), G4TablesForExtrapolator::ComputeProtonDEDX(), G4TablesForExtrapolator::ComputeTrasportXS(), G4PenelopeRayleighModelMI::CrossSectionPerVolume(), G4PAIModelData::CrossSectionPerVolume(), G4PAIPhotData::CrossSectionPerVolume(), G4PAIModelData::DEDXPerVolume(), G4PAIPhotData::DEDXPerVolume(), G4EmModelManager::DumpModelList(), G4EmModelManager::FillDEDXVector(), G4EmModelManager::FillLambdaVector(), G4HadElementSelector::G4HadElementSelector(), G4LivermorePhotoElectricModel::GetBindingEnergy(), G4ICRU90StoppingData::GetDEDX(), G4PAIPhotData::GetEnergyPhotonTransfer(), G4PAIPhotData::GetEnergyPlasmonTransfer(), G4PAIModelData::GetEnergyTransfer(), G4PAIPhotData::GetEnergyTransfer(), G4PAIPhotData::GetPlasmonRatio(), G4PAIPhotData::Initialise(), G4PAIModelData::Initialise(), G4eeToHadronsModel::Initialise(), G4eplusTo2GammaOKVIModel::Initialise(), G4WentzelVIModel::Initialise(), G4DNABornExcitationModel2::Initialise(), G4HadronXSDataTable::Initialise(), G4NeutronCaptureXS::IsoCrossSection(), G4Cerenkov::PostStepDoIt(), G4GDMLWriteMaterials::PropertyVectorWrite(), G4IonICRU73Data::RetrieveVector(), G4PAIPhotData::SampleAlongStepPhotonTransfer(), G4PAIPhotData::SampleAlongStepPlasmonTransfer(), G4PAIPhotData::SampleAlongStepTransfer(), G4PAIModelData::SampleAlongStepTransfer(), G4PAIPhotData::SamplePostStepPhotonTransfer(), G4PAIPhotData::SamplePostStepPlasmonTransfer(), G4PAIPhotData::SamplePostStepTransfer(), G4PAIModelData::SamplePostStepTransfer(), G4DiffuseElasticV2::SampleTableThetaCMS(), G4VEnergyLossProcess::ScaledKinEnergyForLoss(), G4VEnergyLossProcess::SetLambdaTable(), and G4VEmProcess::StreamInfo().

FillSecondDerivatives()

void G4PhysicsVector::FillSecondDerivatives ( const G4SplineType  stype = G4SplineType::Base, const G4double  dir1 = 0.0, const G4double  dir2 = 0.0  )

inherited

Definition at line 205 of file G4PhysicsVector.cc.

{
  if(!useSpline) { return; }
  // cannot compute derivatives for less than 5 points
  const std::size_t nmin = (stype == G4SplineType::Base) ? 5 : 4;
  if(nmin > numberOfNodes) 
  {
    if(0 < verboseLevel)
    { 
      G4cout << "### G4PhysicsVector: spline cannot be used for "
         << numberOfNodes << " points - spline disabled" 
         << G4endl;
      DumpValues();
    }
    useSpline = false;
    return;
  }
  // check energies of free vector
  if(type == T_G4PhysicsFreeVector)
  {
    for(G4int i=0; i<=idxmax; ++i) 
    {
      if(binVector[i + 1] <= binVector[i])
      {
        if(0 < verboseLevel) 
        {
      G4cout << "### G4PhysicsVector: spline cannot be used, because "
         << " E[" << i << "]=" << binVector[i]
         << " >= E[" << i+1 << "]=" << binVector[i + 1]
         << G4endl;
      DumpValues();
        }
        useSpline = false;
        return;
      }
    }
  }
 
  // spline is possible
  Initialise();
  secDerivative.resize(numberOfNodes);
 
  if(1 < verboseLevel)
  {
    G4cout << "### G4PhysicsVector:: FillSecondDerivatives N=" 
           << numberOfNodes << G4endl;
    DumpValues();
  }
 
  switch(stype) 
  {
    case G4SplineType::Base:
      ComputeSecDerivative1();
      break;
 
    case G4SplineType::FixedEdges:
      ComputeSecDerivative2(dir1, dir2);
      break;
 
    default:
      ComputeSecDerivative0();
  }
}

References Base, G4PhysicsVector::binVector, G4PhysicsVector::ComputeSecDerivative0(), G4PhysicsVector::ComputeSecDerivative1(), G4PhysicsVector::ComputeSecDerivative2(), G4PhysicsVector::DumpValues(), FixedEdges, G4cout, G4endl, G4PhysicsVector::idxmax, G4PhysicsVector::Initialise(), G4PhysicsVector::numberOfNodes, G4PhysicsVector::secDerivative, T_G4PhysicsFreeVector, G4PhysicsVector::type, G4PhysicsVector::useSpline, and G4PhysicsVector::verboseLevel.

Referenced by G4WaterStopping::AddData(), G4PSTARStopping::AddData(), G4ASTARStopping::AddData(), G4ICRU90StoppingData::AddData(), G4UPiNuclearCrossSection::AddDataSet(), G4MaterialPropertiesTable::AddProperty(), G4PolarizedCompton::BuildAsymmetryTable(), G4PolarizedAnnihilation::BuildAsymmetryTables(), G4VEnergyLossProcess::BuildDEDXTable(), G4LossTableBuilder::BuildDEDXTable(), G4PenelopeRayleighModel::BuildFormFactorTable(), G4PenelopeRayleighModelMI::BuildFormFactorTable(), G4LossTableBuilder::BuildInverseRangeTable(), G4VEmProcess::BuildLambdaTable(), G4VEnergyLossProcess::BuildLambdaTable(), G4LossTableBuilder::BuildRangeTable(), G4LossTableBuilder::BuildTableForModel(), G4TablesForExtrapolator::ComputeElectronDEDX(), G4TablesForExtrapolator::ComputeMuonDEDX(), G4TablesForExtrapolator::ComputeProtonDEDX(), G4TablesForExtrapolator::ComputeTrasportXS(), G4eplusTo2GammaOKVIModel::Initialise(), G4WentzelVIModel::Initialise(), G4HadronXSDataTable::Initialise(), G4LindhardSorensenData::InitialiseData(), G4PenelopeBremsstrahlungAngular::PrepareTables(), G4ComponentSAIDTotalXS::ReadData(), G4LivermorePhotoElectricModel::ReadData(), G4JAEAElasticScatteringModel::ReadData(), G4JAEAPolarizedElasticScatteringModel::ReadData(), G4LivermoreNuclearGammaConversionModel::ReadData(), G4LivermorePolarizedGammaConversionModel::ReadData(), G4BoldyshevTripletModel::ReadData(), G4IonICRU73Data::ReadElementData(), and G4IonICRU73Data::RetrieveVector().

FindBin()

std::size_t G4PhysicsVector::FindBin ( const G4double  energy, std::size_t  idx  ) const

inherited

Definition at line 173 of file G4PhysicsVector.cc.

{
  if(idx + 1 < numberOfNodes && 
     energy >= binVector[idx] && energy <= binVector[idx])
  {
    return idx;
  } 
  else if(energy <= binVector[1])
  {
    return 0;
  }
  else if(energy >= binVector[idxmax])
  {
    return idxmax;
  }
  return GetBin(energy); 
}

References G4PhysicsVector::binVector, G4INCL::KinematicsUtils::energy(), G4PhysicsVector::GetBin(), G4PhysicsVector::idxmax, and G4PhysicsVector::numberOfNodes.

Referenced by G4PAIModelData::CrossSectionPerVolume(), G4PAIPhotData::CrossSectionPerVolume(), G4PAIPhotData::DEDXPerVolume(), G4SPSEneDistribution::GetArbEneWeight(), G4PAIPhotData::GetPlasmonRatio(), G4PAIPhotData::SampleAlongStepPhotonTransfer(), G4PAIPhotData::SampleAlongStepPlasmonTransfer(), G4PAIPhotData::SampleAlongStepTransfer(), G4PAIModelData::SampleAlongStepTransfer(), G4PAIPhotData::SamplePostStepPhotonTransfer(), G4PAIPhotData::SamplePostStepPlasmonTransfer(), G4PAIPhotData::SamplePostStepTransfer(), G4PAIModelData::SamplePostStepTransfer(), and G4DiffuseElasticV2::SampleTableThetaCMS().

FindLinearEnergy()

G4double G4PhysicsVector::FindLinearEnergy ( const G4double  rand ) const

inlineinherited

Referenced by G4AdjointPrimaryGenerator::SampleDistanceAlongBackRayAndComputeWeightCorrection().

GetBin()

std::size_t G4PhysicsVector::GetBin ( const G4double  energy ) const

inlineprivateinherited

Referenced by G4PhysicsVector::FindBin().

GetEnergy()

G4double G4PhysicsVector::GetEnergy ( const G4double  value ) const

inherited

Definition at line 431 of file G4PhysicsVector.cc.

{
  if(0 == numberOfNodes)
  {
    return 0.0;
  }
  if(1 == numberOfNodes || val <= dataVector[0])
  {
    return edgeMin;
  }
  if(val >= dataVector[numberOfNodes - 1])
  {
    return edgeMax;
  }
  std::size_t bin = 
    std::lower_bound(dataVector.begin(), dataVector.end(), val) -
    dataVector.begin() - 1;
  if(static_cast<G4int>(bin) > idxmax) { bin = idxmax; } 
  G4double res = binVector[bin];
  G4double del = dataVector[bin + 1] - dataVector[bin];
  if(del > 0.0)
  {
    res += (val - dataVector[bin]) * (binVector[bin + 1] - res) / del;
  }
  return res;
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, G4PhysicsVector::edgeMax, G4PhysicsVector::edgeMin, G4PhysicsVector::idxmax, and G4PhysicsVector::numberOfNodes.

Referenced by G4SPSEneDistribution::GenEpnHistEnergies(), G4SPSAngDistribution::GenerateUserDefPhi(), G4SPSAngDistribution::GenerateUserDefTheta(), G4SPSRandomGenerator::GenRandEnergy(), G4SPSRandomGenerator::GenRandPhi(), G4SPSRandomGenerator::GenRandPosPhi(), G4SPSRandomGenerator::GenRandPosTheta(), G4SPSRandomGenerator::GenRandTheta(), G4SPSRandomGenerator::GenRandX(), G4SPSRandomGenerator::GenRandY(), G4SPSRandomGenerator::GenRandZ(), G4SPSEneDistribution::GenUserHistEnergies(), G4Scintillation::PostStepDoIt(), G4OpWLS::PostStepDoIt(), and G4OpWLS2::PostStepDoIt().

GetLowEdgeEnergy()

G4double G4PhysicsVector::GetLowEdgeEnergy ( const std::size_t  index ) const

inlineinherited

Referenced by G4EMDissociation::ApplyYourself(), G4VXTRenergyLoss::BuildAngleForEnergyBank(), G4VXTRenergyLoss::BuildAngleTable(), G4DiffuseElastic::BuildAngleTable(), G4NuclNuclDiffuseElastic::BuildAngleTable(), G4hRDEnergyLoss::BuildDEDXTable(), G4PenelopeIonisationXSHandler::BuildDeltaTable(), G4VXTRenergyLoss::BuildEnergyTable(), G4VXTRenergyLoss::BuildGlobalAngleTable(), G4hRDEnergyLoss::BuildInverseRangeTable(), G4hRDEnergyLoss::BuildLabTimeVector(), G4hImpactIonisation::BuildLambdaTable(), G4hImpactIonisation::BuildLossTable(), G4hRDEnergyLoss::BuildProperTimeVector(), G4hRDEnergyLoss::BuildRangeVector(), G4hhElastic::BuildTableT(), G4ForwardXrayTR::BuildXrayTRtables(), G4PenelopeIonisationXSHandler::BuildXSTable(), G4PenelopeBremsstrahlungModel::BuildXSTable(), G4SPSEneDistribution::ConvertEPNToEnergy(), G4PenelopeRayleighModel::DumpFormFactorTable(), G4PenelopeRayleighModelMI::DumpFormFactorTable(), G4SPSEneDistribution::ExpInterpolation(), G4SPSEneDistribution::GenArbPointEnergies(), G4SPSEneDistribution::GenEpnHistEnergies(), G4SPSAngDistribution::GenerateUserDefPhi(), G4SPSAngDistribution::GenerateUserDefTheta(), G4SPSRandomGenerator::GenRandEnergy(), G4SPSRandomGenerator::GenRandPhi(), G4SPSRandomGenerator::GenRandPosPhi(), G4SPSRandomGenerator::GenRandPosTheta(), G4SPSRandomGenerator::GenRandTheta(), G4SPSRandomGenerator::GenRandX(), G4SPSRandomGenerator::GenRandY(), G4SPSRandomGenerator::GenRandZ(), G4SPSEneDistribution::GenUserHistEnergies(), G4ForwardXrayTR::GetEnergyTR(), G4VXTRenergyLoss::GetMeanFreePath(), G4VXTRenergyLoss::GetNumberOfPhotons(), G4VXTRenergyLoss::GetXTRrandomEnergy(), G4NeutronElectronElXsc::Initialise(), G4NeutronElectronElModel::Initialise(), G4InitXscPAI::IntegralCherenkov(), G4InitXscPAI::IntegralPAIdEdx(), G4InitXscPAI::IntegralPAIxSection(), G4InitXscPAI::IntegralPlasmon(), G4hRDEnergyLoss::InvertRangeVector(), G4SPSEneDistribution::LinearInterpolation(), G4SPSEneDistribution::LogInterpolation(), G4PenelopeCrossSection::NormalizeShellCrossSections(), G4ForwardXrayTR::PostStepDoIt(), G4XNNElasticLowE::Print(), G4XnpElasticLowE::Print(), G4XnpTotalLowE::Print(), G4DiffuseElastic::SampleTableThetaCMS(), G4NuclNuclDiffuseElastic::SampleTableThetaCMS(), and G4SPSEneDistribution::SplineInterpolation().

GetMaxEnergy()

G4double G4PhysicsVector::GetMaxEnergy ( ) const

inlineinherited

Referenced by G4SPSEneDistribution::ArbInterpolate(), G4EmModelManager::FillLambdaVector(), G4PAIModelData::GetEnergyTransfer(), G4GammaNuclearXS::GetIsoCrossSection(), G4Scintillation::GetScintillationYieldByParticleType(), G4NeutronElasticXS::Initialise(), G4NeutronInelasticXS::Initialise(), G4ParticleInelasticXS::Initialise(), G4Cerenkov::PostStepDoIt(), G4PAIModelData::SampleAlongStepTransfer(), G4PAIPhotData::SamplePostStepPhotonTransfer(), G4PAIPhotData::SamplePostStepPlasmonTransfer(), G4PAIPhotData::SamplePostStepTransfer(), G4PAIModelData::SamplePostStepTransfer(), and G4VEmProcess::StreamInfo().

GetMaxValue()

G4double G4PhysicsVector::GetMaxValue ( ) const

inlineinherited

Referenced by G4Scintillation::GetScintillationYieldByParticleType(), G4Cerenkov::PostStepDoIt(), G4Scintillation::PostStepDoIt(), G4OpWLS::PostStepDoIt(), G4OpWLS2::PostStepDoIt(), and G4Cerenkov::PostStepGetPhysicalInteractionLength().

GetMinEnergy()

G4double G4PhysicsVector::GetMinEnergy ( ) const

inlineinherited

GetMinValue()

G4double G4PhysicsVector::GetMinValue ( ) const

inlineinherited

GetSpline()

G4bool G4PhysicsVector::GetSpline ( ) const

inlineinherited

Referenced by G4PolarizedAnnihilation::BuildAsymmetryTables().

GetType()

G4PhysicsVectorType G4PhysicsVector::GetType ( ) const

inlineinherited

GetValue()

G4double G4PhysicsVector::GetValue ( const G4double  energy, G4bool &  isOutRange  ) const

inlineinherited

Referenced by G4hRDEnergyLoss::BuildInverseRangeTable(), G4hRDEnergyLoss::BuildLabTimeVector(), G4hRDEnergyLoss::BuildProperTimeVector(), G4hRDEnergyLoss::BuildRangeCoeffATable(), G4hRDEnergyLoss::BuildRangeCoeffBTable(), G4hRDEnergyLoss::BuildRangeCoeffCTable(), G4hRDEnergyLoss::BuildRangeVector(), G4XNNElasticLowE::CrossSection(), G4XnpElasticLowE::CrossSection(), G4XnpTotalLowE::CrossSection(), G4XResonance::CrossSection(), G4PartialWidthTable::Dump(), G4ChargeExchangeProcess::GetElementCrossSection(), G4ParticleHPJENDLHEData::getXSfromThisIsotope(), G4hRDEnergyLoss::LabTimeIntLog(), G4XNNElasticLowE::Print(), G4XnpElasticLowE::Print(), G4XnpTotalLowE::Print(), G4hRDEnergyLoss::ProperTimeIntLog(), G4hRDEnergyLoss::RangeIntLin(), G4hRDEnergyLoss::RangeIntLog(), G4XAnnihilationChannel::VariablePartialWidth(), and G4XAnnihilationChannel::VariableWidth().

GetVectorLength()

std::size_t G4PhysicsVector::GetVectorLength ( ) const

inlineinherited

Referenced by G4PolarizedAnnihilation::BuildAsymmetryTables(), G4PolarizedIonisation::BuildAsymmetryTables(), G4LossTableBuilder::BuildDEDXTable(), G4PenelopeIonisationXSHandler::BuildDeltaTable(), G4hRDEnergyLoss::BuildInverseRangeTable(), G4LossTableBuilder::BuildInverseRangeTable(), G4Cerenkov::BuildPhysicsTable(), G4Scintillation::BuildPhysicsTable(), G4OpWLS::BuildPhysicsTable(), G4OpWLS2::BuildPhysicsTable(), G4LossTableBuilder::BuildRangeTable(), G4PenelopeIonisationXSHandler::BuildXSTable(), G4PenelopeBremsstrahlungModel::BuildXSTable(), G4MaterialPropertiesTable::CalculateGROUPVEL(), G4OpRayleigh::CalculateRayleighMeanFreePaths(), G4JAEAElasticScatteringModel::ComputeCrossSectionPerAtom(), G4JAEAPolarizedElasticScatteringModel::ComputeCrossSectionPerAtom(), G4LivermoreComptonModel::ComputeCrossSectionPerAtom(), G4LivermoreNuclearGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedComptonModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedRayleighModel::ComputeCrossSectionPerAtom(), G4LivermoreRayleighModel::ComputeCrossSectionPerAtom(), G4LowEPComptonModel::ComputeCrossSectionPerAtom(), G4LowEPPolarizedComptonModel::ComputeCrossSectionPerAtom(), G4SPSEneDistribution::ConvertEPNToEnergy(), G4PAIModelData::CrossSectionPerVolume(), G4PAIPhotData::CrossSectionPerVolume(), G4PAIModelData::DEDXPerVolume(), G4PAIPhotData::DEDXPerVolume(), G4PenelopeRayleighModel::DumpFormFactorTable(), G4PenelopeRayleighModelMI::DumpFormFactorTable(), G4EmModelManager::DumpModelList(), G4SPSEneDistribution::ExpInterpolation(), G4EmModelManager::FillDEDXVector(), G4EmModelManager::FillLambdaVector(), G4SPSEneDistribution::GenArbPointEnergies(), G4SPSEneDistribution::GenEpnHistEnergies(), G4SPSAngDistribution::GenerateUserDefPhi(), G4SPSAngDistribution::GenerateUserDefTheta(), G4SPSRandomGenerator::GenRandEnergy(), G4SPSRandomGenerator::GenRandPhi(), G4SPSRandomGenerator::GenRandPosPhi(), G4SPSRandomGenerator::GenRandPosTheta(), G4SPSRandomGenerator::GenRandTheta(), G4SPSRandomGenerator::GenRandX(), G4SPSRandomGenerator::GenRandY(), G4SPSRandomGenerator::GenRandZ(), G4SPSEneDistribution::GenUserHistEnergies(), G4SPSEneDistribution::GetArbEneWeight(), G4Cerenkov::GetAverageNumberOfPhotons(), G4PAIPhotData::GetEnergyPhotonTransfer(), G4PAIPhotData::GetEnergyPlasmonTransfer(), G4PAIModelData::GetEnergyTransfer(), G4PAIPhotData::GetEnergyTransfer(), G4PenelopeCrossSection::GetHardCrossSection(), G4PenelopeCrossSection::GetNormalizedShellCrossSection(), G4PAIPhotData::GetPlasmonRatio(), G4PenelopeCrossSection::GetShellCrossSection(), G4PenelopeCrossSection::GetSoftStoppingPower(), G4PenelopeCrossSection::GetTotalCrossSection(), G4eeToHadronsModel::Initialise(), G4NeutronElasticXS::Initialise(), G4NeutronInelasticXS::Initialise(), G4ParticleInelasticXS::Initialise(), G4SPSEneDistribution::LinearInterpolation(), G4SPSEneDistribution::LogInterpolation(), G4GDMLWriteMaterials::PropertyVectorWrite(), G4PAIPhotData::SampleAlongStepPhotonTransfer(), G4PAIPhotData::SampleAlongStepPlasmonTransfer(), G4PAIPhotData::SampleAlongStepTransfer(), G4PAIModelData::SampleAlongStepTransfer(), G4PAIPhotData::SamplePostStepPhotonTransfer(), G4PAIPhotData::SamplePostStepPlasmonTransfer(), G4PAIPhotData::SamplePostStepTransfer(), G4PAIModelData::SamplePostStepTransfer(), G4VEnergyLossProcess::SetLambdaTable(), G4SPSEneDistribution::SplineInterpolation(), and G4VEmProcess::StreamInfo().

Initialise()

void G4PhysicsVector::Initialise ( )

protectedvirtualinherited

Reimplemented in G4PhysicsLinearVector, and G4PhysicsLogVector.

Definition at line 44 of file G4PhysicsVector.cc.

{
  idxmax = numberOfNodes - 2;
  if(0 < numberOfNodes)
  {
    edgeMin = binVector[0];
    edgeMax = binVector[numberOfNodes - 1];
  }
}

References G4PhysicsVector::binVector, G4PhysicsVector::edgeMax, G4PhysicsVector::edgeMin, G4PhysicsVector::idxmax, and G4PhysicsVector::numberOfNodes.

Referenced by G4PhysicsVector::FillSecondDerivatives(), G4PhysicsFreeVector(), InsertValues(), G4PhysicsVector::Retrieve(), and G4PhysicsVector::ScaleVector().

InsertValues()

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

Definition at line 133 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();
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, G4INCL::KinematicsUtils::energy(), G4PhysicsVector::Initialise(), and G4PhysicsVector::numberOfNodes.

Referenced by G4MaterialPropertiesTable::AddEntry(), G4SPSEneDistribution::ArbEnergyHisto(), G4SPSEneDistribution::ArbEnergyHistoFile(), G4Cerenkov::BuildPhysicsTable(), G4Scintillation::BuildPhysicsTable(), G4OpWLS::BuildPhysicsTable(), G4OpWLS2::BuildPhysicsTable(), G4MaterialPropertiesTable::CalculateGROUPVEL(), G4OpRayleigh::CalculateRayleighMeanFreePaths(), G4AdjointPrimaryGenerator::ComputeAccumulatedDepthVectorAlongBackRay(), G4SPSEneDistribution::ConvertEPNToEnergy(), G4SPSEneDistribution::EpnEnergyHisto(), G4SPSEneDistribution::ExpInterpolation(), G4SPSEneDistribution::GenEpnHistEnergies(), G4SPSAngDistribution::GenerateUserDefPhi(), G4SPSAngDistribution::GenerateUserDefTheta(), G4SPSRandomGenerator::GenRandEnergy(), G4SPSRandomGenerator::GenRandPhi(), G4SPSRandomGenerator::GenRandPosPhi(), G4SPSRandomGenerator::GenRandPosTheta(), G4SPSRandomGenerator::GenRandTheta(), G4SPSRandomGenerator::GenRandX(), G4SPSRandomGenerator::GenRandY(), G4SPSRandomGenerator::GenRandZ(), G4SPSEneDistribution::GenUserHistEnergies(), G4SPSEneDistribution::LinearInterpolation(), G4SPSEneDistribution::LogInterpolation(), G4GDMLReadMaterials::PropertyRead(), G4GDMLReadSolids::PropertyRead(), G4SPSRandomGenerator::SetEnergyBias(), G4SPSRandomGenerator::SetPhiBias(), G4SPSRandomGenerator::SetPosPhiBias(), G4SPSRandomGenerator::SetPosThetaBias(), G4SPSRandomGenerator::SetThetaBias(), G4SPSRandomGenerator::SetXBias(), G4SPSRandomGenerator::SetYBias(), G4SPSRandomGenerator::SetZBias(), G4SPSEneDistribution::SplineInterpolation(), G4SPSAngDistribution::UserDefAngPhi(), G4SPSAngDistribution::UserDefAngTheta(), and G4SPSEneDistribution::UserEnergyHisto().

Interpolation()

G4double G4PhysicsVector::Interpolation ( const std::size_t  idx, const G4double  energy  ) const

inlineprivateinherited

LogVectorValue()

G4double G4PhysicsVector::LogVectorValue ( const G4double  energy, const G4double  theLogEnergy  ) const

inlineinherited

Referenced by G4GammaGeneralProcess::ComputeGeneralLambda(), G4IonICRU73Data::GetDEDX(), G4GammaGeneralProcess::GetProbability(), G4ParticleInelasticXS::IsoCrossSection(), G4NeutronCaptureXS::IsoCrossSection(), and G4NeutronInelasticXS::IsoCrossSection().

operator!=()

G4bool G4PhysicsVector::operator!= ( const G4PhysicsVector &  right ) const

deleteinherited

operator()()

G4double G4PhysicsVector::operator() ( const std::size_t  index ) const

inlineinherited

operator==()

G4bool G4PhysicsVector::operator== ( const G4PhysicsVector &  right ) const

deleteinherited

operator[]()

G4double G4PhysicsVector::operator[] ( const std::size_t  index ) const

inlineinherited

PrintPutValueError()

void G4PhysicsVector::PrintPutValueError ( std::size_t  index, G4double  value, const G4String &  text  )

protectedinherited

Definition at line 459 of file G4PhysicsVector.cc.

{
  G4ExceptionDescription ed;
  ed << "Vector type: " << type << " length= " << numberOfNodes
     << "; an attempt to put data at index= " << index
     << " value= " << val << " in " << text;
  G4Exception("G4PhysicsVector:", "gl0005", 
              FatalException, ed, "Wrong operation");
}

References FatalException, G4Exception(), G4PhysicsVector::numberOfNodes, and G4PhysicsVector::type.

Referenced by PutValues().

PutValue() [1/2]

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

inline

Definition at line 95 of file G4PhysicsFreeVector.hh.

{
  PutValues(index, e, value);
}

References PutValues().

Referenced by G4PartialWidthTable::AddWidths(), G4DiffuseElastic::BuildAngleTable(), G4NuclNuclDiffuseElastic::BuildAngleTable(), G4PenelopeIonisationXSHandler::BuildDeltaTable(), G4PenelopeRayleighModel::BuildFormFactorTable(), G4PenelopeRayleighModelMI::BuildFormFactorTable(), G4hhElastic::BuildTableT(), G4hhElastic::BuildTableTest(), G4PenelopeRayleighModelMI::CalculateThetaAndAngFun(), G4XDeltaDeltaTable::CrossSectionTable(), G4XNDeltaTable::CrossSectionTable(), G4XDeltaDeltastarTable::CrossSectionTable(), G4XDeltaNstarTable::CrossSectionTable(), G4XNDeltastarTable::CrossSectionTable(), G4XNNstarTable::CrossSectionTable(), G4VXTRenergyLoss::GetAngleVector(), G4EMDissociationCrossSection::GetCrossSectionForProjectile(), G4PenelopeRayleighModel::GetPMaxTable(), G4PenelopeRayleighModelMI::GetPMaxTable(), G4NeutronElectronElModel::Initialise(), G4PAIPhotData::Initialise(), G4PAIModelData::Initialise(), G4BaryonPartialWidth::MassDependentWidth(), G4BaryonWidth::MassDependentWidth(), G4PenelopeGammaConversionModel::ReadDataFile(), G4PenelopePhotoElectricModel::ReadDataFile(), G4PenelopeRayleighModel::ReadDataFile(), G4PenelopeRayleighModelMI::ReadDataFile(), G4PenelopeRayleighModelMI::ReadMolInterferenceData(), G4DiffuseElastic::TestAngleTable(), and G4NuclNuclDiffuseElastic::TestAngleTable().

PutValue() [2/2]

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

inlineinherited

Referenced by G4VXTRenergyLoss::BuildAngleForEnergyBank(), G4PolarizedCompton::BuildAsymmetryTable(), G4PolarizedAnnihilation::BuildAsymmetryTables(), G4PolarizedIonisation::BuildAsymmetryTables(), G4EmCorrections::BuildCorrectionVector(), G4KokoulinMuonNuclearXS::BuildCrossSectionTable(), G4hRDEnergyLoss::BuildDEDXTable(), G4LossTableBuilder::BuildDEDXTable(), G4VXTRenergyLoss::BuildEnergyTable(), G4VXTRenergyLoss::BuildGlobalAngleTable(), G4hRDEnergyLoss::BuildInverseRangeTable(), G4hRDEnergyLoss::BuildLabTimeVector(), G4hImpactIonisation::BuildLambdaTable(), G4hImpactIonisation::BuildLossTable(), G4VLEPTSModel::BuildMeanFreePathTable(), G4VLEPTSModel::BuildPhysicsTable(), G4ChargeExchangeProcess::BuildPhysicsTable(), G4hRDEnergyLoss::BuildProperTimeVector(), G4hRDEnergyLoss::BuildRangeCoeffATable(), G4hRDEnergyLoss::BuildRangeCoeffBTable(), G4hRDEnergyLoss::BuildRangeCoeffCTable(), G4LossTableBuilder::BuildRangeTable(), G4hRDEnergyLoss::BuildRangeVector(), G4LossTableBuilder::BuildTableForModel(), G4AdjointCSManager::BuildTotalSigmaTables(), G4ForwardXrayTR::BuildXrayTRtables(), G4eeToHadronsModel::ComputeCMCrossSectionPerElectron(), G4TablesForExtrapolator::ComputeElectronDEDX(), G4TablesForExtrapolator::ComputeMuonDEDX(), G4TablesForExtrapolator::ComputeProtonDEDX(), G4TablesForExtrapolator::ComputeTrasportXS(), G4EmModelManager::FillDEDXVector(), G4EmModelManager::FillLambdaVector(), G4NeutronElectronElXsc::G4NeutronElectronElXsc(), G4XNNElasticLowE::G4XNNElasticLowE(), G4XnpElasticLowE::G4XnpElasticLowE(), G4XnpTotalLowE::G4XnpTotalLowE(), G4NeutronElectronElXsc::Initialise(), G4PAIPhotData::Initialise(), G4PAIModelData::Initialise(), G4eeToHadronsModel::Initialise(), G4eplusTo2GammaOKVIModel::Initialise(), G4WentzelVIModel::Initialise(), G4DNABornExcitationModel2::Initialise(), G4HadronXSDataTable::Initialise(), G4LindhardSorensenData::InitialiseData(), G4InitXscPAI::IntegralCherenkov(), G4InitXscPAI::IntegralPAIdEdx(), G4InitXscPAI::IntegralPAIxSection(), G4InitXscPAI::IntegralPlasmon(), G4hRDEnergyLoss::InvertRangeVector(), G4IonICRU73Data::ReadElementData(), G4ChannelingECHARM::ReadFromECHARM(), G4IonICRU73Data::RetrieveVector(), and G4ChannelingMaterialData::SetBR().

PutValues()

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

Definition at line 111 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;
  }
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, G4PhysicsVector::edgeMax, G4PhysicsVector::edgeMin, G4PhysicsVector::numberOfNodes, and G4PhysicsVector::PrintPutValueError().

Referenced by G4PenelopeCrossSection::AddCrossSectionPoint(), G4WaterStopping::AddData(), G4PSTARStopping::AddData(), G4ASTARStopping::AddData(), G4ICRU90StoppingData::AddData(), G4UPiNuclearCrossSection::AddDataSet(), G4PenelopeCrossSection::AddShellCrossSectionPoint(), G4LossTableBuilder::BuildInverseRangeTable(), G4PenelopeBremsstrahlungFS::BuildScaledXSTable(), G4ParticleHPData::DoPhysicsVector(), G4EmCorrections::Initialise(), G4PenelopeBremsstrahlungFS::InitializeEnergySampling(), G4PenelopeCrossSection::NormalizeShellCrossSections(), G4PenelopeBremsstrahlungAngular::PrepareTables(), PutValue(), G4ParticleHPJENDLHEData::readAFile(), G4LivermorePhotoElectricModel::ReadData(), and G4PenelopeBremsstrahlungFS::SampleGammaEnergy().

Retrieve()

G4bool G4PhysicsVector::Retrieve ( std::ifstream &  fIn, G4bool  ascii = false  )

inherited

Definition at line 87 of file G4PhysicsVector.cc.

{
  // clear properties;
  dataVector.clear();
  binVector.clear();
  secDerivative.clear();
 
  // retrieve in ascii mode
  if(ascii)
  {
    // binning
    fIn >> edgeMin >> edgeMax >> numberOfNodes;
    if(fIn.fail() || numberOfNodes < 2)
    {
      return false;
    }
    // contents
    G4int siz = 0;
    fIn >> siz;
    if(fIn.fail() || siz != G4int(numberOfNodes))
    {
      return false;
    }
 
    binVector.reserve(siz);
    dataVector.reserve(siz);
    G4double vBin, vData;
 
    for(G4int i = 0; i < siz; ++i)
    {
      vBin  = 0.;
      vData = 0.;
      fIn >> vBin >> vData;
      if(fIn.fail())
      {
        return false;
      }
      binVector.push_back(vBin);
      dataVector.push_back(vData);
    }
    Initialise();
    return true;
  }
 
  // retrieve in binary mode
  // binning
  fIn.read((char*) (&edgeMin), sizeof edgeMin);
  fIn.read((char*) (&edgeMax), sizeof edgeMax);
  fIn.read((char*) (&numberOfNodes), sizeof numberOfNodes);
 
  // contents
  std::size_t size;
  fIn.read((char*) (&size), sizeof size);
 
  G4double* value = new G4double[2 * size];
  fIn.read((char*) (value), 2 * size * (sizeof(G4double)));
  if(G4int(fIn.gcount()) != G4int(2 * size * (sizeof(G4double))))
  {
    delete[] value;
    return false;
  }
 
  binVector.reserve(size);
  dataVector.reserve(size);
  for(std::size_t i = 0; i < size; ++i)
  {
    binVector.push_back(value[2 * i]);
    dataVector.push_back(value[2 * i + 1]);
  }
  delete[] value;
 
  Initialise();
  return true;
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, G4PhysicsVector::edgeMax, G4PhysicsVector::edgeMin, G4PhysicsVector::Initialise(), G4PhysicsVector::numberOfNodes, and G4PhysicsVector::secDerivative.

Referenced by G4ComponentSAIDTotalXS::ReadData(), G4LivermorePhotoElectricModel::ReadData(), G4LivermoreComptonModel::ReadData(), G4LivermoreNuclearGammaConversionModel::ReadData(), G4LivermorePolarizedComptonModel::ReadData(), G4LivermorePolarizedGammaConversionModel::ReadData(), G4LivermoreRayleighModel::ReadData(), G4LowEPComptonModel::ReadData(), G4LowEPPolarizedComptonModel::ReadData(), G4BoldyshevTripletModel::ReadData(), G4LivermoreGammaConversion5DModel::ReadData(), G4LivermoreGammaConversionModel::ReadData(), G4LivermorePolarizedRayleighModel::ReadData(), G4PhysicsTable::RetrievePhysicsTable(), G4GammaNuclearXS::RetrieveVector(), G4IonICRU73Data::RetrieveVector(), G4NeutronCaptureXS::RetrieveVector(), G4NeutronInelasticXS::RetrieveVector(), and G4ParticleInelasticXS::RetrieveVector().

ScaleVector()

void G4PhysicsVector::ScaleVector ( const G4double  factorE, const G4double  factorV  )

inherited

Definition at line 193 of file G4PhysicsVector.cc.

{
  for(std::size_t i = 0; i < numberOfNodes; ++i)
  {
    binVector[i] *= factorE;
    dataVector[i] *= factorV;
  }
  Initialise();
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, G4PhysicsVector::Initialise(), and G4PhysicsVector::numberOfNodes.

Referenced by G4SPSEneDistribution::ExpInterpolation(), G4SPSEneDistribution::LinearInterpolation(), G4SPSEneDistribution::LogInterpolation(), G4ComponentSAIDTotalXS::ReadData(), G4LivermorePhotoElectricModel::ReadData(), G4LivermoreComptonModel::ReadData(), G4LivermorePolarizedComptonModel::ReadData(), G4LowEPComptonModel::ReadData(), G4LowEPPolarizedComptonModel::ReadData(), G4IonICRU73Data::RetrieveVector(), and G4SPSEneDistribution::SplineInterpolation().

SetVerboseLevel()

void G4PhysicsVector::SetVerboseLevel ( G4int  value )

inlineinherited

Referenced by G4MaterialPropertiesTable::AddProperty(), and G4MaterialPropertiesTable::CalculateGROUPVEL().

Store()

G4bool G4PhysicsVector::Store ( std::ofstream &  fOut, G4bool  ascii = false  ) const

inherited

Definition at line 55 of file G4PhysicsVector.cc.

{
  // Ascii mode
  if(ascii)
  {
    fOut << *this;
    return true;
  }
  // Binary Mode
 
  // binning
  fOut.write((char*) (&edgeMin), sizeof edgeMin);
  fOut.write((char*) (&edgeMax), sizeof edgeMax);
  fOut.write((char*) (&numberOfNodes), sizeof numberOfNodes);
 
  // contents
  std::size_t size = dataVector.size();
  fOut.write((char*) (&size), sizeof size);
 
  G4double* value = new G4double[2 * size];
  for(std::size_t i = 0; i < size; ++i)
  {
    value[2 * i]     = binVector[i];
    value[2 * i + 1] = dataVector[i];
  }
  fOut.write((char*) (value), 2 * size * (sizeof(G4double)));
  delete[] value;
 
  return true;
}

References G4PhysicsVector::binVector, G4PhysicsVector::dataVector, G4PhysicsVector::edgeMax, G4PhysicsVector::edgeMin, and G4PhysicsVector::numberOfNodes.

Value() [1/2]

G4double G4PhysicsVector::Value ( const G4double  energy ) const

inlineinherited

Value() [2/2]

G4double G4PhysicsVector::Value ( const G4double  energy, std::size_t &  lastidx  ) const

inlineinherited

Referenced by G4EmCorrections::BarkasCorrection(), G4EmCorrections::BuildCorrectionVector(), G4PenelopeRayleighModelMI::BuildFormFactorTable(), G4LossTableBuilder::BuildRangeTable(), G4OpBoundaryProcess::CalculateReflectivity(), G4Track::CalculateVelocityForOpticalPhoton(), G4eeToHadronsModel::ComputeCMCrossSectionPerElectron(), G4JAEAElasticScatteringModel::ComputeCrossSectionPerAtom(), G4JAEAPolarizedElasticScatteringModel::ComputeCrossSectionPerAtom(), G4LivermoreComptonModel::ComputeCrossSectionPerAtom(), G4LivermoreNuclearGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedComptonModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedRayleighModel::ComputeCrossSectionPerAtom(), G4LivermoreRayleighModel::ComputeCrossSectionPerAtom(), G4LowEPComptonModel::ComputeCrossSectionPerAtom(), G4LowEPPolarizedComptonModel::ComputeCrossSectionPerAtom(), G4PenelopeGammaConversionModel::ComputeCrossSectionPerAtom(), G4PenelopePhotoElectricModel::ComputeCrossSectionPerAtom(), G4PenelopeRayleighModel::ComputeCrossSectionPerAtom(), G4PenelopeRayleighModelMI::ComputeCrossSectionPerAtom(), G4BoldyshevTripletModel::ComputeCrossSectionPerAtom(), G4eplusTo2GammaOKVIModel::ComputeCrossSectionPerAtom(), G4LivermoreGammaConversion5DModel::ComputeCrossSectionPerAtom(), G4LivermoreGammaConversionModel::ComputeCrossSectionPerAtom(), G4eeToHadronsModel::ComputeCrossSectionPerElectron(), G4PolarizedAnnihilation::ComputeSaturationFactor(), G4PolarizedCompton::ComputeSaturationFactor(), G4PolarizedIonisation::ComputeSaturationFactor(), G4eplusTo2GammaOKVIModel::CrossSectionPerVolume(), G4DNABornExcitationModel2::CrossSectionPerVolume(), G4EmCorrections::EffectiveChargeCorrection(), G4eeToHadronsModel::GenerateCMPhoton(), G4LivermorePolarizedRayleighModel::GenerateCosTheta(), G4Cerenkov::GetAverageNumberOfPhotons(), G4ChannelingMaterialData::GetBR(), G4ICRU90StoppingData::GetDEDX(), G4PenelopeIonisationXSHandler::GetDensityCorrection(), G4ChannelingECHARM::GetEC(), G4NeutronElectronElXsc::GetElementCrossSection(), G4KokoulinMuonNuclearXS::GetElementCrossSection(), G4PAIPhotData::GetEnergyPhotonTransfer(), G4PAIPhotData::GetEnergyPlasmonTransfer(), G4PAIModelData::GetEnergyTransfer(), G4PAIPhotData::GetEnergyTransfer(), G4PenelopeRayleighModel::GetFSquared(), G4PenelopeRayleighModelMI::GetFSquared(), G4PenelopeCrossSection::GetHardCrossSection(), G4OpAbsorption::GetMeanFreePath(), G4OpMieHG::GetMeanFreePath(), G4OpRayleigh::GetMeanFreePath(), G4OpWLS::GetMeanFreePath(), G4OpWLS2::GetMeanFreePath(), G4PenelopeCrossSection::GetNormalizedShellCrossSection(), G4SPSEneDistribution::GetProbability(), G4OpBoundaryProcess::GetReflectivity(), G4Scintillation::GetScintillationYieldByParticleType(), G4PenelopePhotoElectricModel::GetShellCrossSection(), G4PenelopeCrossSection::GetShellCrossSection(), G4PenelopeCrossSection::GetSoftStoppingPower(), G4PenelopeCrossSection::GetTotalCrossSection(), G4ElementData::GetValueForElement(), G4NeutronElectronElXsc::Initialise(), G4PAIPhotData::Initialise(), G4eeToHadronsModel::Initialise(), G4EmCorrections::KShellCorrection(), G4EmCorrections::LShellCorrection(), G4Cerenkov::PostStepDoIt(), G4OpBoundaryProcess::PostStepDoIt(), G4PenelopeBremsstrahlungAngular::PrepareTables(), G4DNABornExcitationModel2::RandomSelect(), G4IonICRU73Data::RetrieveVector(), G4PAIModelData::SampleAlongStepTransfer(), G4PenelopeBremsstrahlungAngular::SampleDirection(), G4PAIModelData::SamplePostStepTransfer(), G4PenelopeRayleighModel::SampleSecondaries(), G4PenelopeRayleighModelMI::SampleSecondaries(), G4eplusTo2GammaOKVIModel::SampleSecondaries(), G4VEnergyLossProcess::ScaledKinEnergyForLoss(), G4PenelopePhotoElectricModel::SelectRandomShell(), and G4EmCorrections::ShellCorrection().

Field Documentation

binVector

std::vector<G4double> G4PhysicsVector::binVector

protectedinherited

Definition at line 214 of file G4PhysicsVector.hh.

Referenced by G4PhysicsVector::ComputeSecDerivative0(), G4PhysicsVector::ComputeSecDerivative1(), G4PhysicsVector::ComputeSecDerivative2(), G4PhysicsVector::DumpValues(), G4PhysicsVector::FillSecondDerivatives(), G4PhysicsVector::FindBin(), G4PhysicsFreeVector(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsVector::GetEnergy(), G4PhysicsVector::Initialise(), G4PhysicsLinearVector::Initialise(), G4PhysicsLogVector::Initialise(), InsertValues(), PutValues(), G4PhysicsVector::Retrieve(), G4PhysicsVector::ScaleVector(), and G4PhysicsVector::Store().

dataVector

std::vector<G4double> G4PhysicsVector::dataVector

protectedinherited

Definition at line 215 of file G4PhysicsVector.hh.

Referenced by G4PhysicsVector::ComputeSecDerivative0(), G4PhysicsVector::ComputeSecDerivative1(), G4PhysicsVector::ComputeSecDerivative2(), G4PhysicsVector::DumpValues(), G4PhysicsFreeVector(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsVector::GetEnergy(), InsertValues(), PutValues(), G4PhysicsVector::Retrieve(), G4PhysicsVector::ScaleVector(), and G4PhysicsVector::Store().

edgeMax

G4double G4PhysicsVector::edgeMax = 0.0

protectedinherited

Definition at line 202 of file G4PhysicsVector.hh.

Referenced by G4PhysicsVector::GetEnergy(), G4PhysicsVector::Initialise(), G4PhysicsLinearVector::Initialise(), G4PhysicsLogVector::Initialise(), PutValues(), G4PhysicsVector::Retrieve(), and G4PhysicsVector::Store().

edgeMin

G4double G4PhysicsVector::edgeMin = 0.0

protectedinherited

Definition at line 201 of file G4PhysicsVector.hh.

Referenced by G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsVector::GetEnergy(), G4PhysicsVector::Initialise(), G4PhysicsLinearVector::Initialise(), G4PhysicsLogVector::Initialise(), PutValues(), G4PhysicsVector::Retrieve(), and G4PhysicsVector::Store().

idxmax

G4int G4PhysicsVector::idxmax = 0

protectedinherited

Definition at line 208 of file G4PhysicsVector.hh.

Referenced by G4PhysicsVector::FillSecondDerivatives(), G4PhysicsVector::FindBin(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsVector::GetEnergy(), G4PhysicsVector::Initialise(), G4PhysicsLinearVector::Initialise(), and G4PhysicsLogVector::Initialise().

invdBin

G4double G4PhysicsVector::invdBin = 0.0

protectedinherited

Definition at line 204 of file G4PhysicsVector.hh.

Referenced by G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsLinearVector::Initialise(), and G4PhysicsLogVector::Initialise().

logemin

G4double G4PhysicsVector::logemin = 0.0

protectedinherited

Definition at line 205 of file G4PhysicsVector.hh.

Referenced by G4PhysicsLogVector::Initialise().

numberOfNodes

std::size_t G4PhysicsVector::numberOfNodes = 0

protectedinherited

Definition at line 209 of file G4PhysicsVector.hh.

Referenced by G4PhysicsVector::ComputeSecDerivative0(), G4PhysicsVector::ComputeSecDerivative1(), G4PhysicsVector::ComputeSecDerivative2(), G4PhysicsVector::DumpValues(), G4PhysicsVector::FillSecondDerivatives(), G4PhysicsVector::FindBin(), G4PhysicsFreeVector(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsVector::GetEnergy(), G4PhysicsVector::Initialise(), G4PhysicsLinearVector::Initialise(), G4PhysicsLogVector::Initialise(), InsertValues(), G4PhysicsVector::PrintPutValueError(), PutValues(), G4PhysicsVector::Retrieve(), G4PhysicsVector::ScaleVector(), and G4PhysicsVector::Store().

secDerivative

std::vector<G4double> G4PhysicsVector::secDerivative

protectedinherited

Definition at line 216 of file G4PhysicsVector.hh.

Referenced by G4PhysicsVector::ComputeSecDerivative0(), G4PhysicsVector::ComputeSecDerivative1(), G4PhysicsVector::ComputeSecDerivative2(), G4PhysicsVector::FillSecondDerivatives(), and G4PhysicsVector::Retrieve().

type

G4PhysicsVectorType G4PhysicsVector::type = T_G4PhysicsFreeVector

protectedinherited

Definition at line 211 of file G4PhysicsVector.hh.

Referenced by G4PhysicsVector::FillSecondDerivatives(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLogVector::G4PhysicsLogVector(), and G4PhysicsVector::PrintPutValueError().

useSpline

G4bool G4PhysicsVector::useSpline = false

privateinherited

Definition at line 220 of file G4PhysicsVector.hh.

Referenced by G4PhysicsVector::FillSecondDerivatives().

verboseLevel

G4int G4PhysicsVector::verboseLevel = 0

protectedinherited

Definition at line 207 of file G4PhysicsVector.hh.

Referenced by G4PhysicsVector::FillSecondDerivatives().

---

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

• source/global/management/include/G4PhysicsFreeVector.hh
• source/global/management/src/G4PhysicsFreeVector.cc