#include <G4PenelopeSamplingData.hh>

Public Member Functions
void	AddPoint (G4double x0, G4double pac0, G4double a0, G4double b0, size_t ITTL0, size_t ITTU0)

void	Clear ()

void	DumpTable ()

	G4PenelopeSamplingData (const G4PenelopeSamplingData &)=delete

	G4PenelopeSamplingData (G4int npoints=150)

G4double	GetA (size_t index)

G4double	GetB (size_t index)

size_t	GetNumberOfStoredPoints ()

G4double	GetPAC (size_t index)

G4double	GetX (size_t index)

G4PenelopeSamplingData &	operator= (const G4PenelopeSamplingData &right)=delete

G4double	SampleValue (G4double rndm)

	~G4PenelopeSamplingData ()

Private Attributes
G4DataVector *	fA

G4DataVector *	fB

std::vector< size_t > *	fITTL

std::vector< size_t > *	fITTU

G4int	fNP

G4DataVector *	fPAC

G4DataVector *	fX

Detailed Description

Definition at line 48 of file G4PenelopeSamplingData.hh.

Constructor & Destructor Documentation

◆ G4PenelopeSamplingData() [1/2]

G4PenelopeSamplingData::G4PenelopeSamplingData ( G4int npoints = 150 )

explicit

Definition at line 36 of file G4PenelopeSamplingData.cc.

                                                            : 
  fNP(nPoints)
{
  //create vectors
  fX = new G4DataVector();
  fPAC = new G4DataVector();
  fA = new G4DataVector();
  fB = new G4DataVector();
  fITTL = new std::vector<size_t>;
  fITTU = new std::vector<size_t>;
}

References fA, fB, fITTL, fITTU, fPAC, and fX.

◆ ~G4PenelopeSamplingData()

G4PenelopeSamplingData::~G4PenelopeSamplingData ( )

Definition at line 49 of file G4PenelopeSamplingData.cc.

{
  if (fX) delete fX;
  if (fPAC) delete fPAC;
  if (fA) delete fA;
  if (fB) delete fB;
  if (fITTL) delete fITTL;
  if (fITTU) delete fITTU;
}

References fA, fB, fITTL, fITTU, fPAC, and fX.

◆ G4PenelopeSamplingData() [2/2]

G4PenelopeSamplingData::G4PenelopeSamplingData ( const G4PenelopeSamplingData & )

delete

Member Function Documentation

◆ AddPoint()

void G4PenelopeSamplingData::AddPoint	(	G4double	x0,
		G4double	pac0,
		G4double	a0,
		G4double	b0,
		size_t	ITTL0,
		size_t	ITTU0
	)

Definition at line 97 of file G4PenelopeSamplingData.cc.

{
  fX->push_back(x0);
  fPAC->push_back(pac0);
  fA->push_back(a0);
  fB->push_back(b0);
  fITTL->push_back(ITTL0);
  fITTU->push_back(ITTU0);
 
  //check how many points we do have now
  size_t nOfPoints = GetNumberOfStoredPoints();
 
  if (nOfPoints > ((size_t)fNP))
    {
      G4cout << "G4PenelopeSamplingData::AddPoint() " << G4endl;
      G4cout << "WARNING: Up to now there are " << nOfPoints << " points in the table" << G4endl;
      G4cout << "while the anticipated (declared) number is " << fNP << G4endl;
    }
    return;
}

References a0, fA, fB, fITTL, fITTU, fNP, fPAC, fX, G4cout, G4endl, and GetNumberOfStoredPoints().

Referenced by G4PenelopeRayleighModel::InitializeSamplingAlgorithm(), and G4PenelopeRayleighModelMI::InitializeSamplingAlgorithm().

◆ Clear()

void G4PenelopeSamplingData::Clear ( )

Definition at line 79 of file G4PenelopeSamplingData.cc.

{
  if (fX) delete fX;
  if (fPAC) delete fPAC;
  if (fA) delete fA;
  if (fB) delete fB;
  if (fITTL) delete fITTL;
  if (fITTU) delete fITTU;
  //create vectors
  fX = new G4DataVector();
  fPAC = new G4DataVector();
  fA = new G4DataVector();
  fB = new G4DataVector();
  fITTL = new std::vector<size_t>;
  fITTU = new std::vector<size_t>;
}

References fA, fB, fITTL, fITTU, fPAC, and fX.

◆ DumpTable()

void G4PenelopeSamplingData::DumpTable ( )

Definition at line 120 of file G4PenelopeSamplingData.cc.

{
  
  G4cout << "*************************************************************************" << G4endl;
  G4cout << GetNumberOfStoredPoints() << " points" << G4endl;
  G4cout << "*************************************************************************" << G4endl;
  for (size_t i=0;i<GetNumberOfStoredPoints();i++)
    {
      G4cout << i << " " << (*fX)[i] << " " << (*fPAC)[i] << " " << (*fA)[i] << " " << 
    (*fB)[i] << " " << (*fITTL)[i] << " " << (*fITTU)[i] << G4endl;
    }
  G4cout << "*************************************************************************" << G4endl;
}

References G4cout, G4endl, and GetNumberOfStoredPoints().

Referenced by G4PenelopeRayleighModel::InitializeSamplingAlgorithm(), and G4PenelopeRayleighModelMI::InitializeSamplingAlgorithm().

◆ GetA()

G4double G4PenelopeSamplingData::GetA ( size_t index )

Definition at line 153 of file G4PenelopeSamplingData.cc.

{
  if (index < fA->size())
    return (*fA)[index];
  else
    return 0;
}

References fA.

Referenced by G4PenelopeRayleighModel::GetPMaxTable(), and G4PenelopeRayleighModelMI::GetPMaxTable().

◆ GetB()

G4double G4PenelopeSamplingData::GetB ( size_t index )

Definition at line 162 of file G4PenelopeSamplingData.cc.

{
  if (index < fB->size())
    return (*fB)[index];
  else
    return 0;
}

References fB.

Referenced by G4PenelopeRayleighModel::GetPMaxTable(), and G4PenelopeRayleighModelMI::GetPMaxTable().

◆ GetNumberOfStoredPoints()

size_t G4PenelopeSamplingData::GetNumberOfStoredPoints ( )

Definition at line 61 of file G4PenelopeSamplingData.cc.

{
  size_t points = fX->size();
 
  //check everything is all right
  if (fPAC->size() != points || fA->size() != points || 
      fB->size() != points || fITTL->size() != points ||
      fITTU->size() != points)
    {
      G4ExceptionDescription ed;
      ed << "Data vectors look to have different dimensions !" << G4endl;
      G4Exception("G4PenelopeSamplingData::GetNumberOfStoredPoints()","em2040",
          FatalException,ed);      
    }
  return points;
}

References fA, FatalException, fB, fITTL, fITTU, fPAC, fX, G4endl, and G4Exception().

Referenced by AddPoint(), DumpTable(), G4PenelopeRayleighModel::GetPMaxTable(), G4PenelopeRayleighModelMI::GetPMaxTable(), G4PenelopeRayleighModel::SampleSecondaries(), G4PenelopeRayleighModelMI::SampleSecondaries(), and SampleValue().

◆ GetPAC()

G4double G4PenelopeSamplingData::GetPAC ( size_t index )

Definition at line 144 of file G4PenelopeSamplingData.cc.

{
  if (index < fPAC->size())
    return (*fPAC)[index];
  else
    return 0;
}

References fPAC.

Referenced by G4PenelopeRayleighModel::GetPMaxTable(), and G4PenelopeRayleighModelMI::GetPMaxTable().

◆ GetX()

G4double G4PenelopeSamplingData::GetX ( size_t index )

Definition at line 135 of file G4PenelopeSamplingData.cc.

{
  if (index < fX->size())
    return (*fX)[index];
  else
    return 0;
}

References fX.

Referenced by G4PenelopeRayleighModel::GetPMaxTable(), G4PenelopeRayleighModelMI::GetPMaxTable(), G4PenelopeRayleighModel::SampleSecondaries(), and G4PenelopeRayleighModelMI::SampleSecondaries().

◆ operator=()

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

delete

◆ SampleValue()

G4double G4PenelopeSamplingData::SampleValue ( G4double rndm )

Definition at line 171 of file G4PenelopeSamplingData.cc.

{
  //One passes here a random number in (0,1).
  //Notice: it possible that is between (0,b) with b<1
  size_t points = GetNumberOfStoredPoints();
 
  size_t itn = (size_t) (maxRand*(points-1)); 
  size_t i = (*fITTL)[itn];
  size_t j = (*fITTU)[itn];
 
  while ((j-i) > 1)
    {
      size_t k = (i+j)/2;
      if (maxRand > (*fPAC)[k])
    i = k;
      else
    j = k;
    }
 
  //Sampling from the rational inverse cumulative distribution
  G4double result = 0;
 
  G4double rr = maxRand - (*fPAC)[i];
  if (rr > 1e-16)
    {
      G4double d = (*fPAC)[i+1]-(*fPAC)[i];
      result = (*fX)[i]+
    ((1.0+(*fA)[i]+(*fB)[i])*d*rr/
     (d*d+((*fA)[i]*d+(*fB)[i]*rr)*rr))*((*fX)[i+1]-(*fX)[i]);      
    }
  else
    result = (*fX)[i]; 
  
  return result;
}