#include <G4AdjointCSMatrix.hh>

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

void	Clear ()

	G4AdjointCSMatrix (G4bool aBool)

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 > *	GetLogCrossSectionvector ()

std::vector< double > *	GetLogPrimEnergyVector ()

G4bool	IsScatProjToProj ()

void	Read (G4String file_name)

void	Write (G4String file_name)

	~G4AdjointCSMatrix ()

Private Attributes
std::vector< double >	fLog0Vector

std::vector< double >	fLogCrossSectionVector

std::vector< double >	fLogPrimEnergyVector

std::vector< std::vector< double > * >	fLogProbMatrix

std::vector< std::vector< size_t > * >	fLogProbMatrixIndex

std::vector< std::vector< double > * >	fLogSecondEnergyMatrix

size_t	fNbPrimEnergy = 0

G4bool	fScatProjToProj

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.

36{ fScatProjToProj = aBool; }

G4AdjointCSMatrix::fScatProjToProj

G4bool fScatProjToProj

Definition: G4AdjointCSMatrix.hh:98

References fScatProjToProj.

◆ ~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();
}

References fLogCrossSectionVector, fLogPrimEnergyVector, fLogProbMatrix, fLogProbMatrixIndex, and fLogSecondEnergyMatrix.

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

References G4AdjointInterpolator::FindPosition(), fLog0Vector, fLogCrossSectionVector, fLogPrimEnergyVector, fLogProbMatrix, fLogProbMatrixIndex, fLogSecondEnergyMatrix, fNbPrimEnergy, G4AdjointInterpolator::GetInstance(), and G4INCL::Math::min().

Referenced by G4AdjointCSManager::BuildCrossSectionsModelAndElement(), and G4AdjointCSManager::BuildCrossSectionsModelAndMaterial().

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

References fLog0Vector, fLogCrossSectionVector, fLogPrimEnergyVector, fLogProbMatrix, fLogProbMatrixIndex, fLogSecondEnergyMatrix, and fNbPrimEnergy.

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

References fLog0Vector, fLogCrossSectionVector, fLogPrimEnergyVector, fLogProbMatrix, fLogProbMatrixIndex, fLogSecondEnergyMatrix, and fNbPrimEnergy.

Referenced by G4AdjointCSManager::ComputeAdjointCS(), and G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix().

◆ GetLogCrossSectionvector()

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

inline

Definition at line 71 of file G4AdjointCSMatrix.hh.

  {
    return &fLogCrossSectionVector;
  }

References fLogCrossSectionVector.

◆ GetLogPrimEnergyVector()

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

inline

Definition at line 66 of file G4AdjointCSMatrix.hh.

  {
    return &fLogPrimEnergyVector;
  }

References fLogPrimEnergyVector.

Referenced by G4AdjointCSManager::ComputeAdjointCS(), and G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix().

◆ IsScatProjToProj()

G4bool G4AdjointCSMatrix::IsScatProjToProj ( )

inline

Definition at line 76 of file G4AdjointCSMatrix.hh.

76{ return fScatProjToProj; }

References fScatProjToProj.

Referenced by G4AdjointCSManager::ComputeAdjointCS(), and 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);
    }
  }
}

References fLogCrossSectionVector, fLogPrimEnergyVector, fLogProbMatrix, and fLogSecondEnergyMatrix.

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

References fLogCrossSectionVector, fLogPrimEnergyVector, fLogProbMatrix, fLogSecondEnergyMatrix, G4endl, and MeV.

Field Documentation

◆ fLog0Vector

std::vector<double> G4AdjointCSMatrix::fLog0Vector

private

Definition at line 94 of file G4AdjointCSMatrix.hh.

Referenced by AddData(), Clear(), and GetData().

◆ fLogCrossSectionVector

std::vector<double> G4AdjointCSMatrix::fLogCrossSectionVector

private

Definition at line 85 of file G4AdjointCSMatrix.hh.

Referenced by AddData(), Clear(), GetData(), GetLogCrossSectionvector(), Read(), Write(), and ~G4AdjointCSMatrix().

◆ fLogPrimEnergyVector

std::vector<double> G4AdjointCSMatrix::fLogPrimEnergyVector

private

Definition at line 83 of file G4AdjointCSMatrix.hh.

Referenced by AddData(), Clear(), GetData(), GetLogPrimEnergyVector(), Read(), Write(), and ~G4AdjointCSMatrix().

◆ fLogProbMatrix

std::vector<std::vector<double>*> G4AdjointCSMatrix::fLogProbMatrix

private

Definition at line 88 of file G4AdjointCSMatrix.hh.

Referenced by AddData(), Clear(), GetData(), Read(), Write(), and ~G4AdjointCSMatrix().

◆ fLogProbMatrixIndex

std::vector<std::vector<size_t>*> G4AdjointCSMatrix::fLogProbMatrixIndex

private

Definition at line 93 of file G4AdjointCSMatrix.hh.

Referenced by AddData(), Clear(), GetData(), and ~G4AdjointCSMatrix().

◆ fLogSecondEnergyMatrix

std::vector<std::vector<double>*> G4AdjointCSMatrix::fLogSecondEnergyMatrix

private

Definition at line 87 of file G4AdjointCSMatrix.hh.

Referenced by AddData(), Clear(), GetData(), Read(), Write(), and ~G4AdjointCSMatrix().

◆ fNbPrimEnergy

size_t G4AdjointCSMatrix::fNbPrimEnergy = 0

private

Definition at line 96 of file G4AdjointCSMatrix.hh.

Referenced by AddData(), Clear(), and GetData().

◆ fScatProjToProj

G4bool G4AdjointCSMatrix::fScatProjToProj

private

Definition at line 98 of file G4AdjointCSMatrix.hh.

Referenced by G4AdjointCSMatrix(), and IsScatProjToProj().

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

source/processes/electromagnetic/adjoint/include/G4AdjointCSMatrix.hh
source/processes/electromagnetic/adjoint/src/G4AdjointCSMatrix.cc

Public Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ G4AdjointCSMatrix()

◆ ~G4AdjointCSMatrix()

Member Function Documentation

◆ AddData()

◆ Clear()

◆ GetData()

◆ GetLogCrossSectionvector()

◆ GetLogPrimEnergyVector()

◆ IsScatProjToProj()

◆ Read()

◆ Write()

Field Documentation

◆ fLog0Vector

◆ fLogCrossSectionVector

◆ fLogPrimEnergyVector

◆ fLogProbMatrix

◆ fLogProbMatrixIndex

◆ fLogSecondEnergyMatrix

◆ fNbPrimEnergy

◆ fScatProjToProj