G4ExtendedPhysicsVector Class Reference

#include <G4ExtendedPhysicsVector.hh>

Public Member Functions
	G4ExtendedPhysicsVector (G4PhysicsVector *, G4int nxsec=0)
virtual	~G4ExtendedPhysicsVector ()
	G4ExtendedPhysicsVector (const G4ExtendedPhysicsVector &)=default
G4ExtendedPhysicsVector &	operator= (const G4ExtendedPhysicsVector &)=default
	G4ExtendedPhysicsVector (const G4ExtendedPhysicsVector &&)=delete
G4ExtendedPhysicsVector &	operator= (const G4ExtendedPhysicsVector &&)=delete
G4bool	operator== (const G4ExtendedPhysicsVector &right) const =delete
const G4PhysicsVector *	GetPhysicsVector () const
G4double	Value (const G4double energy, std::size_t &lastidx) const
G4double	LogLogValue (const G4double energy, std::size_t &lastidx) const
void	SetDataLength (G4int dlength)
void	PutPartialXSData (const std::size_t idx, const G4double *y)
G4int	SampleReactionChannel (const G4double energy, const G4double rand, std::size_t &lastidx) const
G4int	SampleReactionChannelLogLog (const G4double energy, const G4double rand, std::size_t &lastidx) const
void	DumpValues (G4double unitE=1.0, G4double unitV=1.0) const

Detailed Description

Definition at line 44 of file G4ExtendedPhysicsVector.hh.

Constructor & Destructor Documentation

G4ExtendedPhysicsVector() [1/3]

G4ExtendedPhysicsVector::G4ExtendedPhysicsVector ( G4PhysicsVector * ptr, G4int nxsec = 0 )

explicit

Definition at line 39 of file G4ExtendedPhysicsVector.cc.

  : totalData(ptr)
{
  if (nullptr == ptr)
  {
    if (0 < verboseLevel)
    {
      G4cout << "### G4ExtendedPhysicsVector with undefined G4PhysicsVector " << G4endl;
    }
    totalData = new G4PhysicsFreeVector(false);
  }
  if (nxsec > 1)
  {
    nPartialXS = nxsec - 1;
    dataPartialXS = new std::vector<std::vector<G4float>* >;
    dataPartialXS->resize((std::size_t)nPartialXS, nullptr);
  }
}

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

~G4ExtendedPhysicsVector()

G4ExtendedPhysicsVector::~G4ExtendedPhysicsVector ( )

virtual

Definition at line 59 of file G4ExtendedPhysicsVector.cc.

{
  if (nPartialXS > 0)
  {
    for (auto const & p : *dataPartialXS) { delete p; }
    delete dataPartialXS;
  }
  delete totalData;
}

G4ExtendedPhysicsVector() [2/3]

G4ExtendedPhysicsVector::G4ExtendedPhysicsVector ( const G4ExtendedPhysicsVector & )

default

G4ExtendedPhysicsVector() [3/3]

G4ExtendedPhysicsVector::G4ExtendedPhysicsVector ( const G4ExtendedPhysicsVector && )

delete

Member Function Documentation

DumpValues()

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

Definition at line 213 of file G4ExtendedPhysicsVector.cc.

{
  G4cout << "====== Data length " << numberOfNodes << " =====" << G4endl;
  // e, partial cumulative normalized x-sections, total x-section
  for (std::size_t i = 0; i < numberOfNodes; ++i)
  {
    G4cout << i << ". " << totalData->Energy(i) / unitE;
    for (G4int j = 0; j < nPartialXS; ++j)
    {
      G4cout << "  " << (*((*dataPartialXS)[i]))[j];
    }
    G4cout << "  " << (*totalData)[i] / unitV << G4endl;
  }
}

GetPhysicsVector()

const G4PhysicsVector * G4ExtendedPhysicsVector::GetPhysicsVector ( ) const

inline

Definition at line 109 of file G4ExtendedPhysicsVector.hh.

{
  return totalData;
}

LogLogValue()

G4double G4ExtendedPhysicsVector::LogLogValue	(	const G4double	energy,
		std::size_t &	lastidx ) const

Definition at line 99 of file G4ExtendedPhysicsVector.cc.

{
  G4bool interpolation = totalData->CheckIndex(e, idx);
 
  G4double y1 = (*totalData)[idx];
  if (interpolation)
  {
    G4double e1 = totalData->Energy(idx);
    G4double e2 = totalData->Energy(idx + 1);
    G4double y2 = (*totalData)[idx + 1];
    if (e1 > 0.0 && e2 > e1 && y1 > 0.0 && y2 > 0.0)
    {
      y1 *= G4Exp(G4Log(e/e1) * G4Log(y2/y1) / G4Log(e2/e1));
    }
  }
  return y1;
}

operator=() [1/2]

G4ExtendedPhysicsVector & G4ExtendedPhysicsVector::operator= ( const G4ExtendedPhysicsVector && )

delete

operator=() [2/2]

G4ExtendedPhysicsVector & G4ExtendedPhysicsVector::operator= ( const G4ExtendedPhysicsVector & )

default

operator==()

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

delete

PutPartialXSData()

void G4ExtendedPhysicsVector::PutPartialXSData	(	const std::size_t	idx,
		const G4double *	y )

Definition at line 118 of file G4ExtendedPhysicsVector.cc.

{
  if (idx >= numberOfNodes)
  {
    if (0 < verboseLevel)
    {
      G4cout << "### G4ExtendedPhysicsVector::PutPartialXSData(..) idx=" << idx
             << " is out of range " << numberOfNodes 
             << G4endl;
    }
    return;
  }
 
  // prepare data vector, the last vector is not created
  // because of normalisation of the sum to 1.
  G4float sum = 0.f;
  for (G4int i = 0; i < nPartialXS; ++i)
  {
    sum += (G4float)y[i];
    (*((*dataPartialXS)[i]))[idx] = sum;
  }
  sum += (G4float)y[nPartialXS];
  if (sum > 0.f)
  {
    sum = 1.f/sum;
    for (G4int i = 0; i < nPartialXS; ++i)
    {
      (*((*dataPartialXS)[i]))[idx] *= sum;
    }
  }
}

Referenced by G4MatUtils::BuildExtendedVector().

SampleReactionChannel()

G4int G4ExtendedPhysicsVector::SampleReactionChannel	(	const G4double	energy,
		const G4double	rand,
		std::size_t &	lastidx ) const

Definition at line 151 of file G4ExtendedPhysicsVector.cc.

{
  if (nPartialXS <= 1) { return 0; }
  G4bool interpolation = totalData->CheckIndex(e, idx);
 
  G4double e1 = totalData->Energy(idx);
  G4double e2 = e1;
  if (interpolation)
  {
    e2 = totalData->Energy(idx + 1) - e1;
    if (e2 <= 0.0) { interpolation = false; }
  }
  for (G4int i=0; i < nPartialXS; ++i)
  {
    G4double xs = (G4double)(*((*dataPartialXS)[i]))[idx];
    if (interpolation)
    {
      G4double xs2 = (G4double)((*((*dataPartialXS)[i]))[idx + 1]) - xs;
      xs += (e - e1) * xs2 / e2;
    }
    if (xs >= rand) { return i; }  
  }
  return nPartialXS;
}

SampleReactionChannelLogLog()

G4int G4ExtendedPhysicsVector::SampleReactionChannelLogLog	(	const G4double	energy,
		const G4double	rand,
		std::size_t &	lastidx ) const

Definition at line 180 of file G4ExtendedPhysicsVector.cc.

{
  if (nPartialXS <= 1) { return 0; }
  G4bool interpolation = totalData->CheckIndex(e, idx);
 
  G4double e1 = totalData->Energy(idx);
  if (e1 <= 0.0) { interpolation = false; }
  G4double e2 = e1;
  if (interpolation)
  {
    e2 = totalData->Energy(idx + 1);
    if (e2 <= e1) { interpolation = false; }
    e2 = G4Log(e2 / e1);
  }
  for (G4int i = 0; i < nPartialXS; ++i)
  {
    G4double xs = (G4double)(*((*dataPartialXS)[i]))[idx];
    if (interpolation && xs > 0.0)
    {
      G4double xs2 = (G4double)((*((*dataPartialXS)[i]))[idx + 1]) - xs;
      if (xs2 > 0.0)
      {
        xs *= G4Exp(G4Log(e / e1) * G4Log(xs2 / xs) / e2);
      }
    }
    if (xs >= rand) { return i; }
  }
  return nPartialXS;
}

SetDataLength()

void G4ExtendedPhysicsVector::SetDataLength

(

G4int

dlength

)

Definition at line 70 of file G4ExtendedPhysicsVector.cc.

{
  // this method may be applied for empty vector only
  totalData->SetDataLength(dlength);
  numberOfNodes = totalData->GetVectorLength();
  if (1 < numberOfNodes)
  {
    idxmax = numberOfNodes - 2;
  }
  else
  {
   if (0 < verboseLevel)
    {
      G4cout << "### G4ExtendedPhysicsVector::SetDataLength numberOfNodes="
             << numberOfNodes << " data structure left empty." << G4endl;
    }
    return;
  }
 
  if (nPartialXS > 0)
  {
    for (G4int i = 0; i < nPartialXS; ++i)
    {
      (*dataPartialXS)[i] = new std::vector<G4float>(numberOfNodes, 0.f); 
    }
  }
}

Referenced by G4MatUtils::BuildExtendedVector().

Value()

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

inline

Definition at line 114 of file G4ExtendedPhysicsVector.hh.

{
  return totalData->Value(e, idx);
}

---

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

• G4ExtendedPhysicsVector.hh
• G4ExtendedPhysicsVector.cc