G4AdjointCSMatrix Class Reference

#include <G4AdjointCSMatrix.hh>

Public Member Functions
	G4AdjointCSMatrix (G4bool aBool)
	~G4AdjointCSMatrix ()
void	Clear ()
void	AddData (G4double aPrimEnergy, G4double aCS, std::vector< double > *aLogSecondEnergyVector, std::vector< double > *aLogProbVector, size_t n_pro_decade=0)
G4bool	GetData (unsigned int i, G4double &aPrimEnergy, G4double &aCS, G4double &log0, std::vector< double > *&aLogSecondEnergyVector, std::vector< double > *&aLogProbVector, std::vector< size_t > *&aLogProbVectorIndex)
std::vector< double > *	GetLogPrimEnergyVector ()
std::vector< double > *	GetLogCrossSectionvector ()
G4bool	IsScatProjToProj ()
void	Write (G4String file_name)
void	Read (G4String file_name)

Detailed Description

Definition at line 49 of file G4AdjointCSMatrix.hh.

Constructor & Destructor Documentation

G4AdjointCSMatrix()

G4AdjointCSMatrix::G4AdjointCSMatrix

(

G4bool

aBool

)

Definition at line 36 of file G4AdjointCSMatrix.cc.

{ fScatProjToProj = aBool; }

~G4AdjointCSMatrix()

G4AdjointCSMatrix::~G4AdjointCSMatrix

(

)

Definition at line 39 of file G4AdjointCSMatrix.cc.

{
  fLogPrimEnergyVector.clear();
  fLogCrossSectionVector.clear();
 
  for (auto p : fLogSecondEnergyMatrix) {
    p->clear();
    delete p;
    p = nullptr;
  }
  fLogSecondEnergyMatrix.clear();
 
  for (auto p : fLogProbMatrix) {
    p->clear();
    delete p;
    p = nullptr;
  }
  fLogProbMatrix.clear();
 
  for (auto p : fLogProbMatrixIndex) {
    if (p) {
      p->clear();
      delete p;
      p = nullptr;
    }
  }
  fLogProbMatrixIndex.clear();
}

Member Function Documentation

AddData()

void G4AdjointCSMatrix::AddData	(	G4double	aPrimEnergy,
		G4double	aCS,
		std::vector< double > *	aLogSecondEnergyVector,
		std::vector< double > *	aLogProbVector,
		size_t	n_pro_decade = 0
	)

Definition at line 81 of file G4AdjointCSMatrix.cc.

{
  G4AdjointInterpolator* theInterpolator = G4AdjointInterpolator::GetInstance();
 
  // At this time we consider that the energy is increasing monotically
  fLogPrimEnergyVector.push_back(aLogPrimEnergy);
  fLogCrossSectionVector.push_back(aLogCS);
  fLogSecondEnergyMatrix.push_back(aLogSecondEnergyVector);
  fLogProbMatrix.push_back(aLogProbVector);
 
  std::vector<size_t>* aLogProbVectorIndex = nullptr;
 
  if(n_pro_decade > 0 && !aLogProbVector->empty())
  {
    aLogProbVectorIndex = new std::vector<size_t>();
    G4double dlog       = std::log(10.) / n_pro_decade;
    G4double log_val =
      int(std::min((*aLogProbVector)[0], aLogProbVector->back()) / dlog) * dlog;
    fLog0Vector.push_back(log_val);
 
    // Loop checking, 07-Aug-2015, Vladimir Ivanchenko
    while(log_val < 0.)
    {
      aLogProbVectorIndex->push_back(
        theInterpolator->FindPosition(log_val, (*aLogProbVector)));
      log_val += dlog;
    }
  }
  else
  {
    fLog0Vector.push_back(0.);
  }
  fLogProbMatrixIndex.push_back(aLogProbVectorIndex);
 
  ++fNbPrimEnergy;
}

Clear()

void G4AdjointCSMatrix::Clear

(

)

Definition at line 69 of file G4AdjointCSMatrix.cc.

{
  fLogPrimEnergyVector.clear();
  fLogCrossSectionVector.clear();
  fLogSecondEnergyMatrix.clear();
  fLogProbMatrix.clear();
  fLogProbMatrixIndex.clear();
  fLog0Vector.clear();
  fNbPrimEnergy = 0;
}

GetData()

G4bool G4AdjointCSMatrix::GetData	(	unsigned int	i,
		G4double &	aPrimEnergy,
		G4double &	aCS,
		G4double &	log0,
		std::vector< double > *&	aLogSecondEnergyVector,
		std::vector< double > *&	aLogProbVector,
		std::vector< size_t > *&	aLogProbVectorIndex
	)

Definition at line 122 of file G4AdjointCSMatrix.cc.

{
  if(i >= fNbPrimEnergy)
    return false;
  aLogPrimEnergy         = fLogPrimEnergyVector[i];
  aLogCS                 = fLogCrossSectionVector[i];
  aLogSecondEnergyVector = fLogSecondEnergyMatrix[i];
  aLogProbVector         = fLogProbMatrix[i];
  aLogProbVectorIndex    = fLogProbMatrixIndex[i];
  log0                   = fLog0Vector[i];
  return true;
}

Referenced by G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix().

GetLogCrossSectionvector()

std::vector< double > * G4AdjointCSMatrix::GetLogCrossSectionvector ( )

inline

Definition at line 71 of file G4AdjointCSMatrix.hh.

  {
    return &fLogCrossSectionVector;
  }

GetLogPrimEnergyVector()

std::vector< double > * G4AdjointCSMatrix::GetLogPrimEnergyVector ( )

inline

Definition at line 66 of file G4AdjointCSMatrix.hh.

  {
    return &fLogPrimEnergyVector;
  }

Referenced by G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix().

IsScatProjToProj()

G4bool G4AdjointCSMatrix::IsScatProjToProj ( )

inline

Definition at line 76 of file G4AdjointCSMatrix.hh.

{ return fScatProjToProj; }

Referenced by G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix().

Read()

void G4AdjointCSMatrix::Read

(

G4String

file_name

)

Definition at line 186 of file G4AdjointCSMatrix.cc.

{
  std::fstream FileOutput(file_name, std::ios::in);
  size_t n1, n2;
 
  fLogPrimEnergyVector.clear();
  fLogCrossSectionVector.clear();
  fLogSecondEnergyMatrix.clear();
  fLogProbMatrix.clear();
  FileOutput >> n1;
  for(size_t i = 0; i < n1; ++i)
  {
    G4double E, CS;
    FileOutput >> E >> CS;
    fLogPrimEnergyVector.push_back(E);
    fLogCrossSectionVector.push_back(CS);
    FileOutput >> n2;
    fLogSecondEnergyMatrix.push_back(new std::vector<G4double>());
    fLogProbMatrix.push_back(new std::vector<G4double>());
 
    for(size_t j = 0; j < n2; ++j)
    {
      G4double E1;
      FileOutput >> E1;
      fLogSecondEnergyMatrix[i]->push_back(E1);
    }
    FileOutput >> n2;
    for(size_t j = 0; j < n2; ++j)
    {
      G4double prob;
      FileOutput >> prob;
      fLogProbMatrix[i]->push_back(prob);
    }
  }
}

Write()

void G4AdjointCSMatrix::Write

(

G4String

file_name

)

Definition at line 140 of file G4AdjointCSMatrix.cc.

{
  std::fstream FileOutput(file_name, std::ios::out);
  FileOutput << std::setiosflags(std::ios::scientific);
  FileOutput << std::setprecision(6);
  FileOutput << fLogPrimEnergyVector.size() << G4endl;
  for(size_t i = 0; i < fLogPrimEnergyVector.size(); ++i)
  {
    FileOutput << std::exp(fLogPrimEnergyVector[i]) / MeV << '\t'
               << std::exp(fLogCrossSectionVector[i]) << G4endl;
    size_t j1 = 0;
    FileOutput << fLogSecondEnergyMatrix[i]->size() << G4endl;
    for(size_t j = 0; j < fLogSecondEnergyMatrix[i]->size(); ++j)
    {
      FileOutput << std::exp((*fLogSecondEnergyMatrix[i])[j]);
      ++j1;
      if(j1 < 10)
        FileOutput << '\t';
      else
      {
        FileOutput << G4endl;
        j1 = 0;
      }
    }
    if(j1 > 0)
      FileOutput << G4endl;
    j1 = 0;
    FileOutput << fLogProbMatrix[i]->size() << G4endl;
    for(size_t j = 0; j < fLogProbMatrix[i]->size(); ++j)
    {
      FileOutput << std::exp((*fLogProbMatrix[i])[j]);
      ++j1;
      if(j1 < 10)
        FileOutput << '\t';
      else
      {
        FileOutput << G4endl;
        j1 = 0;
      }
    }
    if(j1 > 0)
      FileOutput << G4endl;
  }
}

---

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

• G4AdjointCSMatrix.hh
• G4AdjointCSMatrix.cc