G4PhysicsVector Class Reference

#include <G4PhysicsVector.hh>

Inheritance diagram for G4PhysicsVector:

Public Member Functions
	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

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

Protected Attributes
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

Friends
std::ostream &	operator<< (std::ostream &out, const G4PhysicsVector &pv)

Detailed Description

Definition at line 54 of file G4PhysicsVector.hh.

Constructor & Destructor Documentation

G4PhysicsVector() [1/3]

G4PhysicsVector::G4PhysicsVector ( G4bool spline = false )

explicit

Definition at line 39 of file G4PhysicsVector.cc.

  : useSpline(val)
{}

G4PhysicsVector() [2/3]

G4PhysicsVector::G4PhysicsVector ( const G4PhysicsVector & )

default

G4PhysicsVector() [3/3]

G4PhysicsVector::G4PhysicsVector ( const G4PhysicsVector && )

delete

~G4PhysicsVector()

virtual G4PhysicsVector::~G4PhysicsVector ( )

virtualdefault

Member Function Documentation

ComputeLogVectorBin()

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

inline

DumpValues()

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

Definition at line 165 of file G4PhysicsVector.cc.

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

Referenced by G4OpWLS2::DumpPhysicsTable(), G4OpWLS::DumpPhysicsTable(), and FillSecondDerivatives().

Energy()

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

inline

Referenced by G4MaterialPropertiesTable::AddEntry(), G4MaterialPropertiesTable::AddProperty(), G4SPSEneDistribution::ArbInterpolate(), G4DiffuseElasticV2::BuildAngleTable(), G4KokoulinMuonNuclearXS::BuildCrossSectionTable(), G4LossTableBuilder::BuildInverseRangeTable(), G4Cerenkov::BuildPhysicsTable(), G4OpWLS2::BuildPhysicsTable(), G4OpWLS::BuildPhysicsTable(), G4Scintillation::BuildPhysicsTable(), G4LossTableBuilder::BuildRangeTable(), G4PenelopeBremsstrahlungFS::BuildScaledXSTable(), G4LossTableBuilder::BuildTableForModel(), G4AdjointCSManager::BuildTotalSigmaTables(), G4JAEAElasticScatteringModel::ComputeCrossSectionPerAtom(), G4JAEAPolarizedElasticScatteringModel::ComputeCrossSectionPerAtom(), G4LivermoreComptonModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedComptonModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedRayleighModel::ComputeCrossSectionPerAtom(), G4LivermoreRayleighModel::ComputeCrossSectionPerAtom(), G4LowEPComptonModel::ComputeCrossSectionPerAtom(), G4LowEPPolarizedComptonModel::ComputeCrossSectionPerAtom(), G4PAIModelData::CrossSectionPerVolume(), G4PAIPhotData::CrossSectionPerVolume(), G4PenelopeRayleighModelMI::CrossSectionPerVolume(), G4PAIModelData::DEDXPerVolume(), G4PAIPhotData::DEDXPerVolume(), G4EmModelManager::DumpModelList(), G4EmModelManager::FillDEDXVector(), G4EmModelManager::FillLambdaVector(), G4EmUtility::FillPeaksStructure(), G4EmUtility::FindCrossSectionMax(), G4HadElementSelector::G4HadElementSelector(), G4LivermorePhotoElectricModel::GetBindingEnergy(), G4PAIPhotData::GetPlasmonRatio(), G4DNABornExcitationModel2::Initialise(), G4eeToHadronsModel::Initialise(), G4eplusTo2GammaOKVIModel::Initialise(), G4HadronXSDataTable::Initialise(), G4PAIModelData::Initialise(), G4PAIPhotData::Initialise(), G4WentzelVIModel::Initialise(), G4NeutronCaptureXS::IsoCrossSection(), G4Cerenkov::PostStepDoIt(), G4GDMLWriteMaterials::PropertyVectorWrite(), G4PAIPhotData::SampleAlongStepPhotonTransfer(), G4PAIPhotData::SampleAlongStepPlasmonTransfer(), G4PAIModelData::SampleAlongStepTransfer(), G4PAIPhotData::SampleAlongStepTransfer(), G4PAIPhotData::SamplePostStepPhotonTransfer(), G4PAIPhotData::SamplePostStepPlasmonTransfer(), G4PAIModelData::SamplePostStepTransfer(), G4PAIPhotData::SamplePostStepTransfer(), and G4DiffuseElasticV2::SampleTableThetaCMS().

FillSecondDerivatives()

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

Definition at line 207 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 (std::size_t 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();
  }
}

Referenced by G4EmTableUtil::BuildDEDXTable(), G4EmTableUtil::BuildLambdaTable(), G4EmTableUtil::BuildLambdaTable(), G4LossTableBuilder::BuildRangeTable(), G4LossTableBuilder::BuildTableForModel(), G4eplusTo2GammaOKVIModel::Initialise(), G4HadronXSDataTable::Initialise(), G4WentzelVIModel::Initialise(), G4PenelopeBremsstrahlungAngular::PrepareTables(), and G4ExtDEDXTable::RetrievePhysicsTable().

FindBin()

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

Definition at line 175 of file G4PhysicsVector.cc.

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

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

FindLinearEnergy()

G4double G4PhysicsVector::FindLinearEnergy ( const G4double rand ) const

inline

Referenced by G4AdjointPrimaryGenerator::SampleDistanceAlongBackRayAndComputeWeightCorrection().

GetEnergy()

G4double G4PhysicsVector::GetEnergy

(

const G4double

value

)

const

Definition at line 433 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.cbegin(), dataVector.cend(), val)
                  - dataVector.cbegin() - 1;
  if (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;
}

Referenced by G4SPSRandomGenerator::GenRandEnergy(), G4SPSRandomGenerator::GenRandPhi(), G4SPSRandomGenerator::GenRandPosPhi(), G4SPSRandomGenerator::GenRandPosTheta(), G4SPSRandomGenerator::GenRandTheta(), G4SPSRandomGenerator::GenRandX(), G4SPSRandomGenerator::GenRandY(), G4SPSRandomGenerator::GenRandZ(), G4OpWLS2::PostStepDoIt(), G4OpWLS::PostStepDoIt(), and G4Scintillation::PostStepDoIt().

GetLowEdgeEnergy()

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

inline

Referenced by G4EMDissociation::ApplyYourself(), G4VXTRenergyLoss::BuildAngleForEnergyBank(), G4DiffuseElastic::BuildAngleTable(), G4NuclNuclDiffuseElastic::BuildAngleTable(), G4VXTRenergyLoss::BuildAngleTable(), G4hRDEnergyLoss::BuildDEDXTable(), G4VXTRenergyLoss::BuildEnergyTable(), G4VXTRenergyLoss::BuildGlobalAngleTable(), G4hhElastic::BuildTableT(), G4ForwardXrayTR::BuildXrayTRtables(), G4PenelopeIonisationXSHandler::BuildXSTable(), G4PenelopeRayleighModel::DumpFormFactorTable(), G4PenelopeRayleighModelMI::DumpFormFactorTable(), G4SPSRandomGenerator::GenRandEnergy(), G4SPSRandomGenerator::GenRandPhi(), G4SPSRandomGenerator::GenRandPosPhi(), G4SPSRandomGenerator::GenRandPosTheta(), G4SPSRandomGenerator::GenRandTheta(), G4SPSRandomGenerator::GenRandX(), G4SPSRandomGenerator::GenRandY(), G4SPSRandomGenerator::GenRandZ(), G4ForwardXrayTR::GetEnergyTR(), G4VXTRenergyLoss::GetMeanFreePath(), G4VXTRenergyLoss::GetNumberOfPhotons(), G4VXTRenergyLoss::GetXTRrandomEnergy(), G4NeutronElectronElModel::Initialise(), G4NeutronElectronElXsc::Initialise(), G4InitXscPAI::IntegralCherenkov(), G4InitXscPAI::IntegralPAIdEdx(), G4InitXscPAI::IntegralPAIxSection(), G4InitXscPAI::IntegralPlasmon(), G4PenelopeCrossSection::NormalizeShellCrossSections(), G4ForwardXrayTR::PostStepDoIt(), G4XNNElasticLowE::Print(), G4XnpElasticLowE::Print(), G4XnpTotalLowE::Print(), G4DiffuseElastic::SampleTableThetaCMS(), and G4NuclNuclDiffuseElastic::SampleTableThetaCMS().

GetMaxEnergy()

G4double G4PhysicsVector::GetMaxEnergy ( ) const

inline

Referenced by G4SPSEneDistribution::ArbInterpolate(), G4EmModelManager::FillLambdaVector(), G4GammaNuclearXS::GetElementCrossSection(), G4Scintillation::GetScintillationYieldByParticleType(), G4Cerenkov::PostStepDoIt(), G4PAIModelData::SampleAlongStepTransfer(), G4PAIPhotData::SamplePostStepPhotonTransfer(), G4PAIPhotData::SamplePostStepPlasmonTransfer(), G4PAIModelData::SamplePostStepTransfer(), and G4PAIPhotData::SamplePostStepTransfer().

GetMaxValue()

G4double G4PhysicsVector::GetMaxValue ( ) const

inline

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

GetMinEnergy()

G4double G4PhysicsVector::GetMinEnergy ( ) const

inline

GetMinValue()

G4double G4PhysicsVector::GetMinValue ( ) const

inline

GetSpline()

G4bool G4PhysicsVector::GetSpline ( ) const

inline

GetType()

G4PhysicsVectorType G4PhysicsVector::GetType ( ) const

inline

Referenced by G4ExtDEDXTable::StorePhysicsTable().

GetValue()

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

inline

Referenced by G4XNNElasticLowE::CrossSection(), G4XnpElasticLowE::CrossSection(), G4XnpTotalLowE::CrossSection(), G4XResonance::CrossSection(), G4PartialWidthTable::Dump(), G4ChargeExchangeProcess::GetElementCrossSection(), G4XNNElasticLowE::Print(), G4XnpElasticLowE::Print(), and G4XnpTotalLowE::Print().

GetVectorLength()

std::size_t G4PhysicsVector::GetVectorLength ( ) const

inline

Referenced by G4MaterialPropertiesTable::AddEntry(), G4MaterialPropertiesTable::AddProperty(), G4LossTableBuilder::BuildDEDXTable(), G4LossTableBuilder::BuildInverseRangeTable(), G4Cerenkov::BuildPhysicsTable(), G4OpWLS2::BuildPhysicsTable(), G4OpWLS::BuildPhysicsTable(), G4Scintillation::BuildPhysicsTable(), G4LossTableBuilder::BuildRangeTable(), G4PenelopeIonisationXSHandler::BuildXSTable(), G4JAEAElasticScatteringModel::ComputeCrossSectionPerAtom(), G4JAEAPolarizedElasticScatteringModel::ComputeCrossSectionPerAtom(), G4LivermoreComptonModel::ComputeCrossSectionPerAtom(), G4LivermoreNuclearGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedComptonModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedRayleighModel::ComputeCrossSectionPerAtom(), G4LivermoreRayleighModel::ComputeCrossSectionPerAtom(), G4LowEPComptonModel::ComputeCrossSectionPerAtom(), G4LowEPPolarizedComptonModel::ComputeCrossSectionPerAtom(), G4PAIModelData::CrossSectionPerVolume(), G4PAIPhotData::CrossSectionPerVolume(), G4PAIModelData::DEDXPerVolume(), G4PAIPhotData::DEDXPerVolume(), G4PenelopeRayleighModel::DumpFormFactorTable(), G4PenelopeRayleighModelMI::DumpFormFactorTable(), G4EmModelManager::DumpModelList(), G4EmModelManager::FillDEDXVector(), G4EmModelManager::FillLambdaVector(), G4EmUtility::FillPeaksStructure(), G4EmUtility::FindCrossSectionMax(), G4SPSRandomGenerator::GenRandEnergy(), G4SPSRandomGenerator::GenRandPhi(), G4SPSRandomGenerator::GenRandPosPhi(), G4SPSRandomGenerator::GenRandPosTheta(), G4SPSRandomGenerator::GenRandTheta(), G4SPSRandomGenerator::GenRandX(), G4SPSRandomGenerator::GenRandY(), G4SPSRandomGenerator::GenRandZ(), G4SPSEneDistribution::GetArbEneWeight(), G4Cerenkov::GetAverageNumberOfPhotons(), G4PenelopeCrossSection::GetHardCrossSection(), G4PenelopeCrossSection::GetNormalizedShellCrossSection(), G4PAIPhotData::GetPlasmonRatio(), G4PenelopeCrossSection::GetShellCrossSection(), G4PenelopeCrossSection::GetSoftStoppingPower(), G4PenelopeCrossSection::GetTotalCrossSection(), G4eeToHadronsModel::Initialise(), G4GDMLWriteMaterials::PropertyVectorWrite(), G4PAIPhotData::SampleAlongStepPhotonTransfer(), G4PAIPhotData::SampleAlongStepPlasmonTransfer(), G4PAIModelData::SampleAlongStepTransfer(), G4PAIPhotData::SampleAlongStepTransfer(), G4PAIPhotData::SamplePostStepPhotonTransfer(), G4PAIPhotData::SamplePostStepPlasmonTransfer(), G4PAIModelData::SamplePostStepTransfer(), and G4PAIPhotData::SamplePostStepTransfer().

Initialise()

void G4PhysicsVector::Initialise ( )

protectedvirtual

Reimplemented in G4PhysicsLinearVector, and G4PhysicsLogVector.

Definition at line 44 of file G4PhysicsVector.cc.

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

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

LogFreeVectorValue()

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

inline

LogVectorValue()

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

inline

Referenced by G4GammaGeneralProcess::ComputeGeneralLambda(), G4NeutronGeneralProcess::ComputeGeneralLambda(), G4IonICRU73Data::GetDEDX(), G4KokoulinMuonNuclearXS::GetElementCrossSection(), G4GammaGeneralProcess::GetProbability(), G4NeutronGeneralProcess::GetProbability(), and G4NeutronCaptureXS::IsoCrossSection().

operator!=()

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

delete

operator()()

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

inline

operator=() [1/2]

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

delete

operator=() [2/2]

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

default

operator==()

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

delete

operator[]()

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

inline

PrintPutValueError()

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

protected

Definition at line 460 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");
}

Referenced by G4PhysicsFreeVector::PutValues().

PutValue()

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

inline

Referenced by G4KokoulinMuonNuclearXS::BuildCrossSectionTable(), G4hRDEnergyLoss::BuildDEDXTable(), G4ChargeExchangeProcess::BuildPhysicsTable(), G4LossTableBuilder::BuildRangeTable(), G4LossTableBuilder::BuildTableForModel(), G4AdjointCSManager::BuildTotalSigmaTables(), G4EmModelManager::FillDEDXVector(), G4EmModelManager::FillLambdaVector(), G4NeutronElectronElXsc::G4NeutronElectronElXsc(), G4XNNElasticLowE::G4XNNElasticLowE(), G4XnpElasticLowE::G4XnpElasticLowE(), G4XnpTotalLowE::G4XnpTotalLowE(), G4DNABornExcitationModel2::Initialise(), G4eeToHadronsModel::Initialise(), G4eplusTo2GammaOKVIModel::Initialise(), G4HadronXSDataTable::Initialise(), G4NeutronElectronElXsc::Initialise(), G4WentzelVIModel::Initialise(), G4InitXscPAI::IntegralCherenkov(), G4InitXscPAI::IntegralPAIdEdx(), G4InitXscPAI::IntegralPAIxSection(), G4InitXscPAI::IntegralPlasmon(), G4ChannelingECHARM::ReadFromECHARM(), G4ChannelingMaterialData::SetBR(), and G4ChannelingMaterialData::SetBR().

Retrieve()

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

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 siz0 = 0;
    fIn >> siz0;
    if (siz0 < 2) { return false; }
    std::size_t siz = static_cast<std::size_t>(siz0);
    if (fIn.fail() || siz != numberOfNodes)
    {
      return false;
    }
 
    binVector.reserve(siz);
    dataVector.reserve(siz);
    G4double vBin, vData;
 
    for (std::size_t 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 (static_cast<G4int>(fIn.gcount()) != static_cast<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;
}

Referenced by G4ExtDEDXTable::RetrievePhysicsTable(), and G4PhysicsTable::RetrievePhysicsTable().

ScaleVector()

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

Definition at line 195 of file G4PhysicsVector.cc.

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

SetVerboseLevel()

void G4PhysicsVector::SetVerboseLevel ( G4int value )

inline

Store()

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

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;
}

Referenced by G4ExtDEDXTable::StorePhysicsTable().

Value() [1/2]

G4double G4PhysicsVector::Value ( const G4double energy ) const

inline

Value() [2/2]

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

inline

Referenced by G4EmCorrections::BarkasCorrection(), G4Track::CalculateVelocityForOpticalPhoton(), G4BoldyshevTripletModel::ComputeCrossSectionPerAtom(), G4eplusTo2GammaOKVIModel::ComputeCrossSectionPerAtom(), G4JAEAElasticScatteringModel::ComputeCrossSectionPerAtom(), G4JAEAPolarizedElasticScatteringModel::ComputeCrossSectionPerAtom(), G4LivermoreComptonModel::ComputeCrossSectionPerAtom(), G4LivermoreGammaConversion5DModel::ComputeCrossSectionPerAtom(), G4LivermoreGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermoreNuclearGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedComptonModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermorePolarizedRayleighModel::ComputeCrossSectionPerAtom(), G4LivermoreRayleighModel::ComputeCrossSectionPerAtom(), G4LowEPComptonModel::ComputeCrossSectionPerAtom(), G4LowEPPolarizedComptonModel::ComputeCrossSectionPerAtom(), G4PenelopeGammaConversionModel::ComputeCrossSectionPerAtom(), G4PenelopePhotoElectricModel::ComputeCrossSectionPerAtom(), G4PenelopeRayleighModel::ComputeCrossSectionPerAtom(), G4PenelopeRayleighModelMI::ComputeCrossSectionPerAtom(), G4eeToHadronsModel::ComputeCrossSectionPerElectron(), G4DNABornExcitationModel2::CrossSectionPerVolume(), G4eplusTo2GammaOKVIModel::CrossSectionPerVolume(), G4EmCorrections::EffectiveChargeCorrection(), G4eeToHadronsModel::GenerateCMPhoton(), G4Cerenkov::GetAverageNumberOfPhotons(), G4ChannelingMaterialData::GetBR(), G4PenelopeIonisationXSHandler::GetDensityCorrection(), G4ChannelingECHARM::GetEC(), G4NeutronElectronElXsc::GetElementCrossSection(), G4PenelopeCrossSection::GetHardCrossSection(), G4OpAbsorption::GetMeanFreePath(), G4OpMieHG::GetMeanFreePath(), G4OpWLS2::GetMeanFreePath(), G4OpWLS::GetMeanFreePath(), G4PenelopeCrossSection::GetNormalizedShellCrossSection(), G4SPSEneDistribution::GetProbability(), G4Scintillation::GetScintillationYieldByParticleType(), G4PenelopeCrossSection::GetShellCrossSection(), G4PenelopePhotoElectricModel::GetShellCrossSection(), G4PenelopeCrossSection::GetSoftStoppingPower(), G4PenelopeCrossSection::GetTotalCrossSection(), G4eeToHadronsModel::Initialise(), G4NeutronElectronElXsc::Initialise(), G4EmCorrections::KShellCorrection(), G4GammaNuclearXS::LowEnergyCrossSection(), G4EmCorrections::LShellCorrection(), G4Cerenkov::PostStepDoIt(), G4OpBoundaryProcess::PostStepDoIt(), G4PenelopeBremsstrahlungAngular::PrepareTables(), G4XrayReflection::Reflectivity(), G4PAIModelData::SampleAlongStepTransfer(), G4PenelopeBremsstrahlungAngular::SampleDirection(), G4PAIModelData::SamplePostStepTransfer(), G4eplusTo2GammaOKVIModel::SampleSecondaries(), G4LivermorePhotoElectricModel::SampleSecondaries(), G4PenelopeRayleighModel::SampleSecondaries(), G4PenelopeRayleighModelMI::SampleSecondaries(), and G4EmCorrections::ShellCorrection().

Friends And Related Symbol Documentation

operator<<

std::ostream & operator<< ( std::ostream & out, const G4PhysicsVector & pv )

friend

Definition at line 414 of file G4PhysicsVector.cc.

{
  // binning
  G4long prec = out.precision();
  out << std::setprecision(12) << pv.edgeMin << " " << pv.edgeMax << " "
      << pv.numberOfNodes << G4endl;
 
  // contents
  out << pv.dataVector.size() << G4endl;
  for (std::size_t i = 0; i < pv.dataVector.size(); ++i)
  {
    out << pv.binVector[i] << "  " << pv.dataVector[i] << G4endl;
  }
  out.precision(prec);
 
  return out;
}

Member Data Documentation

binVector

std::vector<G4double> G4PhysicsVector::binVector

protected

Definition at line 226 of file G4PhysicsVector.hh.

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

dataVector

std::vector<G4double> G4PhysicsVector::dataVector

protected

Definition at line 227 of file G4PhysicsVector.hh.

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

edgeMax

G4double G4PhysicsVector::edgeMax = 0.0

protected

Definition at line 209 of file G4PhysicsVector.hh.

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

edgeMin

G4double G4PhysicsVector::edgeMin = 0.0

protected

Definition at line 208 of file G4PhysicsVector.hh.

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

iBin1

G4double G4PhysicsVector::iBin1 = 0.0

protected

Definition at line 214 of file G4PhysicsVector.hh.

Referenced by G4PhysicsFreeVector::EnableLogBinSearch().

idxmax

std::size_t G4PhysicsVector::idxmax = 0

protected

Definition at line 218 of file G4PhysicsVector.hh.

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

imax1

std::size_t G4PhysicsVector::imax1 = 0

protected

Definition at line 219 of file G4PhysicsVector.hh.

Referenced by G4PhysicsFreeVector::EnableLogBinSearch().

invdBin

G4double G4PhysicsVector::invdBin = 0.0

protected

Definition at line 211 of file G4PhysicsVector.hh.

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

lmin1

G4double G4PhysicsVector::lmin1 = 0.0

protected

Definition at line 215 of file G4PhysicsVector.hh.

Referenced by G4PhysicsFreeVector::EnableLogBinSearch().

logemin

G4double G4PhysicsVector::logemin = 0.0

protected

Definition at line 212 of file G4PhysicsVector.hh.

Referenced by G4PhysicsLogVector::Initialise().

nLogNodes

std::size_t G4PhysicsVector::nLogNodes = 0

protected

Definition at line 221 of file G4PhysicsVector.hh.

Referenced by G4PhysicsFreeVector::EnableLogBinSearch().

numberOfNodes

std::size_t G4PhysicsVector::numberOfNodes = 0

protected

Definition at line 220 of file G4PhysicsVector.hh.

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

scale

std::vector<std::size_t> G4PhysicsVector::scale

protected

Definition at line 229 of file G4PhysicsVector.hh.

Referenced by G4PhysicsFreeVector::EnableLogBinSearch().

secDerivative

std::vector<G4double> G4PhysicsVector::secDerivative

protected

Definition at line 228 of file G4PhysicsVector.hh.

Referenced by FillSecondDerivatives(), and Retrieve().

type

G4PhysicsVectorType G4PhysicsVector::type = T_G4PhysicsFreeVector

protected

Definition at line 223 of file G4PhysicsVector.hh.

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

verboseLevel

G4int G4PhysicsVector::verboseLevel = 0

protected

Definition at line 217 of file G4PhysicsVector.hh.

Referenced by FillSecondDerivatives().

---

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

• G4PhysicsVector.hh
• G4PhysicsVector.cc