G4PhysicsVector Class Reference

#include <G4PhysicsVector.hh>

Inheritance diagram for G4PhysicsVector:

Public Member Functions
	G4PhysicsVector (G4bool spline=false)
	G4PhysicsVector (const G4PhysicsVector &)
G4PhysicsVector &	operator= (const G4PhysicsVector &)
virtual	~G4PhysicsVector ()
void *	operator new (size_t)
void	operator delete (void *)
G4double	Value (G4double theEnergy)
G4double	GetValue (G4double theEnergy, G4bool &isOutRange)
G4int	operator== (const G4PhysicsVector &right) const
G4int	operator!= (const G4PhysicsVector &right) const
G4double	operator[] (const size_t binNumber) const
G4double	operator() (const size_t binNumber) const
void	PutValue (size_t index, G4double theValue)
virtual void	ScaleVector (G4double factorE, G4double factorV)
G4double	Energy (size_t index) const
G4double	GetMaxEnergy () const
virtual G4double	GetLowEdgeEnergy (size_t binNumber) const
size_t	GetVectorLength () const
void	FillSecondDerivatives ()
void	ComputeSecDerivatives ()
void	ComputeSecondDerivatives (G4double firstPointDerivative, G4double endPointDerivative)
G4bool	IsFilledVectorExist () const
G4PhysicsVectorType	GetType () const
void	SetSpline (G4bool)
virtual G4bool	Store (std::ofstream &fOut, G4bool ascii=false)
virtual G4bool	Retrieve (std::ifstream &fIn, G4bool ascii=false)
G4double	GetLastEnergy () const
G4double	GetLastValue () const
size_t	GetLastBin () const
void	SetVerboseLevel (G4int value)
G4int	GetVerboseLevel (G4int)

Protected Member Functions
virtual size_t	FindBinLocation (G4double theEnergy) const =0
void	DeleteData ()
void	CopyData (const G4PhysicsVector &vec)

Protected Attributes
G4PhysicsVectorType	type
G4double	edgeMin
G4double	edgeMax
size_t	numberOfNodes
G4PhysicsVectorCache *	cache
G4PVDataVector	dataVector
G4PVDataVector	binVector
G4PVDataVector	secDerivative
G4double	dBin
G4double	baseBin
G4int	verboseLevel

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

Detailed Description

Definition at line 76 of file G4PhysicsVector.hh.

Constructor & Destructor Documentation

G4PhysicsVector() [1/2]

G4PhysicsVector::G4PhysicsVector

(

G4bool

spline = false

)

Definition at line 63 of file G4PhysicsVector.cc.

 : type(T_G4PhysicsVector),
   edgeMin(0.), edgeMax(0.), numberOfNodes(0),
   useSpline(spline), 
   dBin(0.), baseBin(0.),
   verboseLevel(0)
{
  cache      = new G4PhysicsVectorCache();
}

G4PhysicsVector() [2/2]

G4PhysicsVector::G4PhysicsVector

(

const G4PhysicsVector &

right

)

Definition at line 82 of file G4PhysicsVector.cc.

{
  cache        = new G4PhysicsVectorCache();
  dBin         = right.dBin;
  baseBin      = right.baseBin;
  verboseLevel = right.verboseLevel;
 
  DeleteData();
  CopyData(right);
}

~G4PhysicsVector()

G4PhysicsVector::~G4PhysicsVector ( )

virtual

Definition at line 75 of file G4PhysicsVector.cc.

{
  delete cache; cache =0;
}

Member Function Documentation

ComputeSecDerivatives()

void G4PhysicsVector::ComputeSecDerivatives

(

)

Definition at line 441 of file G4PhysicsVector.cc.

{
  if(!SplinePossible())  { return; }
 
  if(3 > numberOfNodes)  // cannot compute derivatives for less than 4 bins
  {
    useSpline = false;
    return;
  }
 
  size_t n = numberOfNodes-1;
 
  for(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];
}

Referenced by ComputeSecondDerivatives(), and FillSecondDerivatives().

ComputeSecondDerivatives()

void G4PhysicsVector::ComputeSecondDerivatives	(	G4double	firstPointDerivative,
		G4double	endPointDerivative
	)

Definition at line 308 of file G4PhysicsVector.cc.

{
  if(4 > numberOfNodes)   // cannot compute derivatives for less than 4 bins
  {
    ComputeSecDerivatives();
    return;
  }
 
  if(!SplinePossible()) { return; }
 
  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;
}

CopyData()

void G4PhysicsVector::CopyData ( const G4PhysicsVector &  vec )

protected

Definition at line 130 of file G4PhysicsVector.cc.

{
  type = vec.type;
  edgeMin = vec.edgeMin;
  edgeMax = vec.edgeMax;
  numberOfNodes = vec.numberOfNodes;
  cache->lastEnergy = vec.GetLastEnergy();
  cache->lastValue = vec.GetLastValue();
  cache->lastBin = vec.GetLastBin();
  useSpline = vec.useSpline;
 
  size_t i;
  dataVector.clear();
  for(i=0; i<(vec.dataVector).size(); i++){ 
    dataVector.push_back( (vec.dataVector)[i] );
  }
  binVector.clear();
  for(i=0; i<(vec.binVector).size(); i++){ 
    binVector.push_back( (vec.binVector)[i] );
  }
  secDerivative.clear();
  for(i=0; i<(vec.secDerivative).size(); i++){ 
    secDerivative.push_back( (vec.secDerivative)[i] );
  }
}

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

DeleteData()

void G4PhysicsVector::DeleteData ( )

protected

Definition at line 123 of file G4PhysicsVector.cc.

{
  secDerivative.clear();
}

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

Energy()

G4double G4PhysicsVector::Energy ( size_t  index ) const

inline

Referenced by G4EmCorrections::BarkasCorrection(), G4KokoulinMuonNuclearXS::BuildCrossSectionTable(), G4LossTableBuilder::BuildInverseRangeTable(), G4ePolarizedIonisation::BuildPhysicsTable(), G4LossTableBuilder::BuildRangeTable(), G4LossTableBuilder::BuildTableForModel(), G4Scintillation::BuildThePhysicsTable(), G4LivermoreRayleighModel::ComputeCrossSectionPerAtom(), G4EmModelManager::DumpModelList(), G4EmModelManager::FillDEDXVector(), G4EmModelManager::FillLambdaVector(), G4NeutronCaptureXS::GetElementCrossSection(), G4NeutronElasticXS::GetElementCrossSection(), G4NeutronInelasticXS::GetElementCrossSection(), G4NeutronCaptureXS::GetIsoCrossSection(), G4GDMLWriteMaterials::PropertyVectorWrite(), and G4VEnergyLossProcess::SetLambdaTable().

FillSecondDerivatives()

void G4PhysicsVector::FillSecondDerivatives

(

)

Definition at line 371 of file G4PhysicsVector.cc.

{  
  if(5 > numberOfNodes)  // cannot compute derivatives for less than 4 points
  {
    ComputeSecDerivatives();
    return;
  }
 
  if(!SplinePossible()) { return; }
 
  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;
}

Referenced by G4VEnergyLossProcess::BuildDEDXTable(), G4LossTableBuilder::BuildDEDXTable(), G4LossTableBuilder::BuildInverseRangeTable(), G4VEnergyLossProcess::BuildLambdaTable(), G4LossTableBuilder::BuildRangeTable(), and G4LossTableBuilder::BuildTableForModel().

FindBinLocation()

virtual size_t G4PhysicsVector::FindBinLocation ( G4double  theEnergy ) const

protectedpure virtual

Implemented in G4PhysicsFreeVector, G4PhysicsLinearVector, G4PhysicsLnVector, and G4PhysicsLogVector.

GetLastBin()

size_t G4PhysicsVector::GetLastBin ( ) const

inline

Referenced by CopyData().

GetLastEnergy()

G4double G4PhysicsVector::GetLastEnergy ( ) const

inline

Referenced by CopyData().

GetLastValue()

G4double G4PhysicsVector::GetLastValue ( ) const

inline

Referenced by CopyData().

GetLowEdgeEnergy()

G4double G4PhysicsVector::GetLowEdgeEnergy ( size_t  binNumber ) const

virtual

Reimplemented in G4PhysicsOrderedFreeVector.

Definition at line 158 of file G4PhysicsVector.cc.

{
  return binVector[binNumber];
}

Referenced by G4EMDissociation::ApplyYourself(), G4VXTRenergyLoss::BuildAngleForEnergyBank(), G4VXTRenergyLoss::BuildAngleTable(), G4DiffuseElastic::BuildAngleTable(), G4NuclNuclDiffuseElastic::BuildAngleTable(), G4eplusPolarizedAnnihilation::BuildAsymmetryTable(), G4PolarizedCompton::BuildAsymmetryTable(), G4hRDEnergyLoss::BuildDEDXTable(), G4VXTRenergyLoss::BuildEnergyTable(), G4VXTRenergyLoss::BuildGlobalAngleTable(), G4VeLowEnergyLoss::BuildInverseRangeTable(), G4VRangeToEnergyConverter::BuildLossTable(), G4PAIModel::BuildPAIonisationTable(), G4PAIPhotonModel::BuildPAIonisationTable(), G4VRangeToEnergyConverter::BuildRangeVector(), G4AdjointCSManager::BuildTotalSigmaTables(), G4ForwardXrayTR::BuildXrayTRtables(), G4VRangeToEnergyConverter::ConvertCutToKineticEnergy(), G4PenelopeRayleighModel::DumpFormFactorTable(), G4PAIPhotonModel::GetAlongStepTransfer(), G4ForwardXrayTR::GetEnergyTR(), G4VXTRenergyLoss::GetMeanFreePath(), G4PAIModel::GetPostStepTransfer(), G4PAIPhotonModel::GetPostStepTransfer(), G4VXTRenergyLoss::GetXTRrandomEnergy(), G4eeToHadronsModel::Initialise(), G4InitXscPAI::IntegralCherenkov(), G4InitXscPAI::IntegralPAIdEdx(), G4InitXscPAI::IntegralPAIxSection(), G4InitXscPAI::IntegralPlasmon(), G4ForwardXrayTR::PostStepDoIt(), G4XNNElasticLowE::Print(), G4XnpElasticLowE::Print(), G4XnpTotalLowE::Print(), G4PAIModel::SampleFluctuations(), G4DiffuseElastic::SampleTableThetaCMS(), and G4NuclNuclDiffuseElastic::SampleTableThetaCMS().

GetMaxEnergy()

G4double G4PhysicsVector::GetMaxEnergy ( ) const

inline

Referenced by G4VEmProcess::PrintInfoDefinition().

GetType()

G4PhysicsVectorType G4PhysicsVector::GetType ( ) const

inline

GetValue()

G4double G4PhysicsVector::GetValue ( G4double  theEnergy, G4bool &  isOutRange  )

inline

Referenced by G4VeLowEnergyLoss::BuildInverseRangeTable(), G4VeLowEnergyLoss::BuildRangeCoeffATable(), G4VEnergyLoss::BuildRangeCoeffATable(), G4VeLowEnergyLoss::BuildRangeCoeffBTable(), G4VEnergyLoss::BuildRangeCoeffBTable(), G4VeLowEnergyLoss::BuildRangeCoeffCTable(), G4VEnergyLoss::BuildRangeCoeffCTable(), G4eeToHadronsModel::ComputeCrossSectionPerElectron(), G4XNNElasticLowE::CrossSection(), G4XnpElasticLowE::CrossSection(), G4XnpTotalLowE::CrossSection(), G4XResonance::CrossSection(), G4PartialWidthTable::Dump(), G4eeToHadronsModel::GenerateCMPhoton(), G4ChargeExchangeProcess::GetElementCrossSection(), G4PolarizedCompton::GetMeanFreePath(), G4eeToHadronsModel::Initialise(), G4PolarizedCompton::PostStepGetPhysicalInteractionLength(), G4XNNElasticLowE::Print(), G4XnpElasticLowE::Print(), and G4XnpTotalLowE::Print().

GetVectorLength()

size_t G4PhysicsVector::GetVectorLength ( ) const

inline

Referenced by G4LossTableBuilder::BuildDEDXTable(), G4LossTableBuilder::BuildInverseRangeTable(), G4VeLowEnergyLoss::BuildInverseRangeTable(), G4ePolarizedIonisation::BuildPhysicsTable(), G4LossTableBuilder::BuildRangeTable(), G4Scintillation::BuildThePhysicsTable(), G4AdjointCSManager::BuildTotalSigmaTables(), G4LivermoreGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermoreRayleighModel::ComputeCrossSectionPerAtom(), G4PenelopeRayleighModel::DumpFormFactorTable(), G4EmModelManager::DumpModelList(), G4EmModelManager::FillDEDXVector(), G4EmModelManager::FillLambdaVector(), G4SPSRandomGenerator::GenRandEnergy(), G4SPSRandomGenerator::GenRandPhi(), G4SPSRandomGenerator::GenRandPosPhi(), G4SPSRandomGenerator::GenRandPosTheta(), G4SPSRandomGenerator::GenRandTheta(), G4SPSRandomGenerator::GenRandX(), G4SPSRandomGenerator::GenRandY(), G4SPSRandomGenerator::GenRandZ(), G4NeutronElasticXS::GetElementCrossSection(), G4NeutronInelasticXS::GetElementCrossSection(), G4PenelopeCrossSection::GetHardCrossSection(), G4PenelopeCrossSection::GetNormalizedShellCrossSection(), G4PenelopeCrossSection::GetShellCrossSection(), G4PenelopeCrossSection::GetSoftStoppingPower(), G4PenelopeCrossSection::GetTotalCrossSection(), G4eeToHadronsModel::Initialise(), G4VEmProcess::PrintInfoDefinition(), G4GDMLWriteMaterials::PropertyVectorWrite(), and G4VEnergyLossProcess::SetLambdaTable().

GetVerboseLevel()

G4int G4PhysicsVector::GetVerboseLevel ( G4int  )

inline

IsFilledVectorExist()

G4bool G4PhysicsVector::IsFilledVectorExist ( ) const

inline

operator delete()

void G4PhysicsVector::operator delete ( void *  )

inline

operator new()

void * G4PhysicsVector::operator new ( size_t  )

inline

operator!=()

G4int G4PhysicsVector::operator!=

(

const G4PhysicsVector &

right

)

const

Definition at line 116 of file G4PhysicsVector.cc.

{
  return (this != &right);
}

operator()()

G4double G4PhysicsVector::operator() ( const size_t  binNumber ) const

inline

operator=()

G4PhysicsVector & G4PhysicsVector::operator=

(

const G4PhysicsVector &

right

)

Definition at line 95 of file G4PhysicsVector.cc.

{
  if (&right==this)  { return *this; }
  dBin         = right.dBin;
  baseBin      = right.baseBin;
  verboseLevel = right.verboseLevel;
 
  DeleteData();
  CopyData(right);
  return *this;
}

operator==()

G4int G4PhysicsVector::operator==

(

const G4PhysicsVector &

right

)

const

Definition at line 109 of file G4PhysicsVector.cc.

{
  return (this == &right);
}

operator[]()

G4double G4PhysicsVector::operator[] ( const size_t  binNumber ) const

inline

PutValue()

void G4PhysicsVector::PutValue ( size_t  index, G4double  theValue  )

inline

Referenced by G4RToEConvForGamma::BuildAbsorptionLengthVector(), G4VXTRenergyLoss::BuildAngleForEnergyBank(), G4eplusPolarizedAnnihilation::BuildAsymmetryTable(), G4PolarizedCompton::BuildAsymmetryTable(), G4KokoulinMuonNuclearXS::BuildCrossSectionTable(), G4hRDEnergyLoss::BuildDEDXTable(), G4LossTableBuilder::BuildDEDXTable(), G4VXTRenergyLoss::BuildEnergyTable(), G4VXTRenergyLoss::BuildGlobalAngleTable(), G4VeLowEnergyLoss::BuildInverseRangeTable(), G4PAIModel::BuildLambdaVector(), G4PAIPhotonModel::BuildLambdaVector(), G4VRangeToEnergyConverter::BuildLossTable(), G4PAIModel::BuildPAIonisationTable(), G4PAIPhotonModel::BuildPAIonisationTable(), G4ePolarizedIonisation::BuildPhysicsTable(), G4ChargeExchangeProcess::BuildPhysicsTable(), G4VeLowEnergyLoss::BuildRangeCoeffATable(), G4VEnergyLoss::BuildRangeCoeffATable(), G4VeLowEnergyLoss::BuildRangeCoeffBTable(), G4VEnergyLoss::BuildRangeCoeffBTable(), G4VeLowEnergyLoss::BuildRangeCoeffCTable(), G4VEnergyLoss::BuildRangeCoeffCTable(), G4LossTableBuilder::BuildRangeTable(), G4VRangeToEnergyConverter::BuildRangeVector(), G4LossTableBuilder::BuildTableForModel(), G4AdjointCSManager::BuildTotalSigmaTables(), G4ForwardXrayTR::BuildXrayTRtables(), G4EmModelManager::FillDEDXVector(), G4EmModelManager::FillLambdaVector(), G4XNNElasticLowE::G4XNNElasticLowE(), G4XnpElasticLowE::G4XnpElasticLowE(), G4XnpTotalLowE::G4XnpTotalLowE(), G4eeToHadronsModel::Initialise(), G4InitXscPAI::IntegralCherenkov(), G4InitXscPAI::IntegralPAIdEdx(), G4InitXscPAI::IntegralPAIxSection(), G4InitXscPAI::IntegralPlasmon(), and G4VRangeToEnergyConverter::operator=().

Retrieve()

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

virtual

Reimplemented in G4PhysicsLinearVector, G4PhysicsLnVector, and G4PhysicsLogVector.

Definition at line 198 of file G4PhysicsVector.cc.

{
  // clear properties;
  cache->lastEnergy=-DBL_MAX;
  cache->lastValue =0.;
  cache->lastBin   =0;
  dataVector.clear();
  binVector.clear();
  secDerivative.clear();
 
  // retrieve in ascii mode
  if (ascii){
    // binning
    fIn >> edgeMin >> edgeMax >> numberOfNodes; 
    if (fIn.fail())  { return false; }
    // contents
    G4int siz=0;
    fIn >> siz;
    if (fIn.fail())  { return false; }
    if (siz<=0)
    {
#ifdef G4VERBOSE  
      G4cerr << "G4PhysicsVector::Retrieve():";
      G4cerr << " Invalid vector size: " << siz << G4endl;
#endif
      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);
    }
 
    // to remove any inconsistency 
    numberOfNodes = siz;
    edgeMin = binVector[0];
    edgeMax = binVector[numberOfNodes-1];
    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
  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(size_t i = 0; i < size; ++i)
  {
    binVector.push_back(value[2*i]);
    dataVector.push_back(value[2*i+1]);
  }
  delete [] value;
 
  // to remove any inconsistency 
  numberOfNodes = size;
  edgeMin = binVector[0];
  edgeMax = binVector[numberOfNodes-1];
 
  return true;
}

Referenced by G4PhysicsLinearVector::Retrieve(), G4PhysicsLnVector::Retrieve(), G4PhysicsLogVector::Retrieve(), and G4PhysicsTable::RetrievePhysicsTable().

ScaleVector()

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

virtual

Reimplemented in G4PhysicsLinearVector, G4PhysicsLnVector, and G4PhysicsLogVector.

Definition at line 285 of file G4PhysicsVector.cc.

{
  size_t n = dataVector.size();
  size_t i;
  if(n > 0) { 
    for(i=0; i<n; ++i) {
      binVector[i]  *= factorE;
      dataVector[i] *= factorV;
    } 
  }
  //  n = secDerivative.size();
  // if(n > 0) { for(i=0; i<n; ++i) { secDerivative[i] *= factorV; } }
  secDerivative.clear();
 
  edgeMin *= factorE;
  edgeMax *= factorE;
  cache->lastEnergy = factorE*(cache->lastEnergy);
  cache->lastValue  = factorV*(cache->lastValue);
}

Referenced by G4PhysicsLinearVector::ScaleVector(), G4PhysicsLnVector::ScaleVector(), and G4PhysicsLogVector::ScaleVector().

SetSpline()

void G4PhysicsVector::SetSpline ( G4bool  )

inline

Referenced by G4VEnergyLossProcess::BuildDEDXTable(), G4LossTableBuilder::BuildDEDXTable(), G4LossTableBuilder::BuildInverseRangeTable(), G4VEnergyLossProcess::BuildLambdaTable(), G4LossTableBuilder::BuildRangeTable(), G4LossTableBuilder::BuildTableForModel(), G4EmElementSelector::G4EmElementSelector(), G4HadronNucleonXsc::InitialiseKaonNucleonTotXsc(), G4VEmProcess::LambdaPhysicsVector(), G4VEnergyLossProcess::LambdaPhysicsVector(), G4ee2KChargedModel::PhysicsVector(), G4ee2KNeutralModel::PhysicsVector(), G4eeTo3PiModel::PhysicsVector(), G4eeToPGammaModel::PhysicsVector(), and G4eeToTwoPiModel::PhysicsVector().

SetVerboseLevel()

void G4PhysicsVector::SetVerboseLevel ( G4int  value )

inline

Store()

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

virtual

Definition at line 165 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
  size_t size = dataVector.size(); 
  fOut.write((char*)(&size), sizeof size);
 
  G4double* value = new G4double[2*size];
  for(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;
}

Value()

G4double G4PhysicsVector::Value ( G4double  theEnergy )

inline

Referenced by G4EmCorrections::BarkasCorrection(), G4LossTableBuilder::BuildRangeTable(), G4Track::CalculateVelocityForOpticalPhoton(), G4LivermoreGammaConversionModel::ComputeCrossSectionPerAtom(), G4LivermoreRayleighModel::ComputeCrossSectionPerAtom(), G4PenelopeGammaConversionModel::ComputeCrossSectionPerAtom(), G4PenelopePhotoElectricModel::ComputeCrossSectionPerAtom(), G4PenelopeRayleighModel::ComputeCrossSectionPerAtom(), G4VRangeToEnergyConverter::ConvertCutToKineticEnergy(), G4EmCorrections::EffectiveChargeCorrection(), G4PenelopeIonisationXSHandler::GetDensityCorrection(), G4NeutronCaptureXS::GetElementCrossSection(), G4NeutronElasticXS::GetElementCrossSection(), G4NeutronInelasticXS::GetElementCrossSection(), G4PenelopeCrossSection::GetHardCrossSection(), G4NeutronCaptureXS::GetIsoCrossSection(), G4HadronNucleonXsc::GetKmNeutronTotXscVector(), G4HadronNucleonXsc::GetKmProtonTotXscVector(), G4HadronNucleonXsc::GetKpNeutronTotXscVector(), G4HadronNucleonXsc::GetKpProtonTotXscVector(), G4PenelopeCrossSection::GetNormalizedShellCrossSection(), G4SPSEneDistribution::GetProbability(), G4PenelopePhotoElectricModel::GetShellCrossSection(), G4PenelopeCrossSection::GetShellCrossSection(), G4PenelopeCrossSection::GetSoftStoppingPower(), G4PenelopeCrossSection::GetTotalCrossSection(), G4ElementData::GetValueForElement(), G4EmCorrections::KShellCorrection(), G4EmCorrections::LShellCorrection(), G4Cerenkov::PostStepDoIt(), G4Scintillation::PostStepDoIt(), G4OpBoundaryProcess::PostStepDoIt(), G4PenelopeBremsstrahlungAngular::SampleDirection(), G4PenelopeRayleighModel::SampleSecondaries(), G4NeutronCaptureXS::SelectIsotope(), G4VEnergyLossProcess::SetCSDARangeTable(), G4VEnergyLossProcess::SetDEDXTable(), and G4EmCorrections::ShellCorrection().

Friends And Related Function Documentation

operator<<

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

friend

Definition at line 487 of file G4PhysicsVector.cc.

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

Member Data Documentation

baseBin

G4double G4PhysicsVector::baseBin

protected

Definition at line 232 of file G4PhysicsVector.hh.

Referenced by G4PhysicsLinearVector::FindBinLocation(), G4PhysicsLnVector::FindBinLocation(), G4PhysicsLogVector::FindBinLocation(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLnVector::G4PhysicsLnVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsVector(), operator=(), G4PhysicsLinearVector::Retrieve(), G4PhysicsLnVector::Retrieve(), G4PhysicsLogVector::Retrieve(), G4PhysicsLinearVector::ScaleVector(), G4PhysicsLnVector::ScaleVector(), and G4PhysicsLogVector::ScaleVector().

binVector

G4PVDataVector G4PhysicsVector::binVector

protected

Definition at line 212 of file G4PhysicsVector.hh.

Referenced by ComputeSecDerivatives(), ComputeSecondDerivatives(), CopyData(), G4LPhysicsFreeVector::DumpValues(), FillSecondDerivatives(), G4PhysicsFreeVector::FindBinLocation(), G4LPhysicsFreeVector::G4LPhysicsFreeVector(), G4PhysicsFreeVector::G4PhysicsFreeVector(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLnVector::G4PhysicsLnVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsOrderedFreeVector::GetLowEdgeEnergy(), GetLowEdgeEnergy(), G4PhysicsOrderedFreeVector::InsertValues(), G4PhysicsFreeVector::PutValue(), G4PhysicsLinearVector::Retrieve(), G4PhysicsLnVector::Retrieve(), G4PhysicsLogVector::Retrieve(), Retrieve(), G4PhysicsLinearVector::ScaleVector(), G4PhysicsLnVector::ScaleVector(), G4PhysicsLogVector::ScaleVector(), ScaleVector(), and Store().

cache

G4PhysicsVectorCache* G4PhysicsVector::cache

protected

Definition at line 209 of file G4PhysicsVector.hh.

Referenced by CopyData(), G4PhysicsVector(), Retrieve(), ScaleVector(), and ~G4PhysicsVector().

dataVector

G4PVDataVector G4PhysicsVector::dataVector

protected

Definition at line 211 of file G4PhysicsVector.hh.

Referenced by ComputeSecDerivatives(), ComputeSecondDerivatives(), CopyData(), G4LPhysicsFreeVector::DumpValues(), FillSecondDerivatives(), G4LPhysicsFreeVector::G4LPhysicsFreeVector(), G4PhysicsFreeVector::G4PhysicsFreeVector(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLnVector::G4PhysicsLnVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsOrderedFreeVector::InsertValues(), G4PhysicsFreeVector::PutValue(), Retrieve(), ScaleVector(), and Store().

dBin

G4double G4PhysicsVector::dBin

protected

Definition at line 231 of file G4PhysicsVector.hh.

Referenced by G4PhysicsLinearVector::FindBinLocation(), G4PhysicsLnVector::FindBinLocation(), G4PhysicsLogVector::FindBinLocation(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLnVector::G4PhysicsLnVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsVector(), operator=(), G4PhysicsLinearVector::Retrieve(), G4PhysicsLnVector::Retrieve(), G4PhysicsLogVector::Retrieve(), G4PhysicsLinearVector::ScaleVector(), G4PhysicsLnVector::ScaleVector(), and G4PhysicsLogVector::ScaleVector().

edgeMax

G4double G4PhysicsVector::edgeMax

protected

Definition at line 205 of file G4PhysicsVector.hh.

Referenced by CopyData(), G4LPhysicsFreeVector::G4LPhysicsFreeVector(), G4PhysicsFreeVector::G4PhysicsFreeVector(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLnVector::G4PhysicsLnVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsOrderedFreeVector::InsertValues(), G4PhysicsFreeVector::PutValue(), Retrieve(), ScaleVector(), and Store().

edgeMin

G4double G4PhysicsVector::edgeMin

protected

Definition at line 204 of file G4PhysicsVector.hh.

Referenced by CopyData(), G4LPhysicsFreeVector::G4LPhysicsFreeVector(), G4PhysicsFreeVector::G4PhysicsFreeVector(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLnVector::G4PhysicsLnVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsOrderedFreeVector::InsertValues(), G4PhysicsFreeVector::PutValue(), Retrieve(), ScaleVector(), and Store().

numberOfNodes

size_t G4PhysicsVector::numberOfNodes

protected

Definition at line 207 of file G4PhysicsVector.hh.

Referenced by ComputeSecDerivatives(), ComputeSecondDerivatives(), CopyData(), G4LPhysicsFreeVector::DumpValues(), FillSecondDerivatives(), G4PhysicsFreeVector::FindBinLocation(), G4LPhysicsFreeVector::G4LPhysicsFreeVector(), G4PhysicsFreeVector::G4PhysicsFreeVector(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLnVector::G4PhysicsLnVector(), G4PhysicsLogVector::G4PhysicsLogVector(), G4PhysicsOrderedFreeVector::InsertValues(), G4PhysicsFreeVector::PutValue(), Retrieve(), and Store().

secDerivative

G4PVDataVector G4PhysicsVector::secDerivative

protected

Definition at line 213 of file G4PhysicsVector.hh.

Referenced by ComputeSecDerivatives(), ComputeSecondDerivatives(), CopyData(), DeleteData(), FillSecondDerivatives(), Retrieve(), and ScaleVector().

type

G4PhysicsVectorType G4PhysicsVector::type

protected

Definition at line 202 of file G4PhysicsVector.hh.

Referenced by CopyData(), G4LPhysicsFreeVector::G4LPhysicsFreeVector(), G4PhysicsFreeVector::G4PhysicsFreeVector(), G4PhysicsLinearVector::G4PhysicsLinearVector(), G4PhysicsLnVector::G4PhysicsLnVector(), G4PhysicsLogVector::G4PhysicsLogVector(), and G4PhysicsOrderedFreeVector::G4PhysicsOrderedFreeVector().

verboseLevel

G4int G4PhysicsVector::verboseLevel

protected

Definition at line 234 of file G4PhysicsVector.hh.

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

---

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

• G4PhysicsVector.hh
• G4PhysicsVector.cc