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 ()
G4double	GetDlog ()
G4bool	IsScatProjToProjCase ()
void	Write (G4String file_name)
void	Read (G4String file_name)

Detailed Description

Definition at line 60 of file G4AdjointCSMatrix.hh.

Constructor & Destructor Documentation

G4AdjointCSMatrix::G4AdjointCSMatrix

(

G4bool

aBool

)

Definition at line 38 of file G4AdjointCSMatrix.cc.

                                                {
    theLogPrimEnergyVector.clear();
    theLogCrossSectionVector.clear();
    theLogSecondEnergyMatrix.clear();
    theLogProbMatrix.clear();
    theLogProbMatrixIndex.clear();
    log0Vector.clear();
    nb_of_PrimEnergy=0;
    is_scat_proj_to_proj_case  =aBool;
    dlog =0;
}

G4AdjointCSMatrix::~G4AdjointCSMatrix

(

)

Definition at line 51 of file G4AdjointCSMatrix.cc.

                                     {
    theLogPrimEnergyVector.clear();
    theLogCrossSectionVector.clear();
    theLogSecondEnergyMatrix.clear();
    theLogProbMatrix.clear();
}

Member Function Documentation

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

Definition at line 71 of file G4AdjointCSMatrix.cc.

References G4AdjointInterpolator::FindPosition(), G4AdjointInterpolator::GetInstance(), int(), and G4INCL::Math::min().

                                                                                              {
    
    G4AdjointInterpolator* theInterpolator=G4AdjointInterpolator::GetInstance();
    
    //At this time we consider that the energy is increasing monotically
    theLogPrimEnergyVector.push_back(aLogPrimEnergy);
    theLogCrossSectionVector.push_back(aLogCS);
    theLogSecondEnergyMatrix.push_back(aLogSecondEnergyVector);
    theLogProbMatrix.push_back(aLogProbVector);
    
    std::vector< size_t>* aLogProbVectorIndex = 0;
    dlog =0;
    
    if (n_pro_decade > 0 && aLogProbVector->size()>0) {
        aLogProbVectorIndex = new std::vector< size_t>();
        dlog=std::log(10.)/n_pro_decade;
        G4double log_val = int(std::min((*aLogProbVector)[0],aLogProbVector->back())/dlog)*dlog;
        log0Vector.push_back(log_val);
        
        while(log_val<0.) {
            aLogProbVectorIndex->push_back(theInterpolator->FindPosition(log_val,(*aLogProbVector)));
            log_val+=dlog;
        } 
    }
    else {
        log0Vector.push_back(0.);
    }
    theLogProbMatrixIndex.push_back(aLogProbVectorIndex);
    
    
    nb_of_PrimEnergy++;
    
    
}

void G4AdjointCSMatrix::Clear

(

)

Definition at line 59 of file G4AdjointCSMatrix.cc.

{
    theLogPrimEnergyVector.clear();
    theLogCrossSectionVector.clear();
    theLogSecondEnergyMatrix.clear();
    theLogProbMatrix.clear();
    theLogProbMatrixIndex.clear();
    log0Vector.clear();
    nb_of_PrimEnergy=0;
}

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 108 of file G4AdjointCSMatrix.cc.

Referenced by G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix().

{   if (i>= nb_of_PrimEnergy) return false;
    //G4cout<<"Test Get Data "<<G4endl;
    aLogPrimEnergy = theLogPrimEnergyVector[i];
    aLogCS = theLogCrossSectionVector[i];
    aLogSecondEnergyVector = theLogSecondEnergyMatrix[i];
    aLogProbVector = theLogProbMatrix[i];
    aLogProbVectorIndex = theLogProbMatrixIndex[i];
    log0=log0Vector[i];
    return true;
    
}

G4double G4AdjointCSMatrix::GetDlog ( )

inline

Definition at line 82 of file G4AdjointCSMatrix.hh.

{return dlog;}  

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

inline

Definition at line 81 of file G4AdjointCSMatrix.hh.

{return &theLogCrossSectionVector;}

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

inline

Definition at line 80 of file G4AdjointCSMatrix.hh.

Referenced by G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix().

{return &theLogPrimEnergyVector;}

G4bool G4AdjointCSMatrix::IsScatProjToProjCase ( )

inline

Definition at line 83 of file G4AdjointCSMatrix.hh.

Referenced by G4VEmAdjointModel::SampleAdjSecEnergyFromCSMatrix().

{return is_scat_proj_to_proj_case;} 

void G4AdjointCSMatrix::Read

(

G4String

file_name

)

Definition at line 161 of file G4AdjointCSMatrix.cc.

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

void G4AdjointCSMatrix::Write

(

G4String

file_name

)

Definition at line 123 of file G4AdjointCSMatrix.cc.

References G4endl, and python.hepunit::MeV.

{   std::fstream FileOutput(file_name, std::ios::out);            
    FileOutput<<std::setiosflags(std::ios::scientific);
    FileOutput<<std::setprecision(6);
    FileOutput<<theLogPrimEnergyVector.size()<<G4endl;
    for (size_t i=0;i<theLogPrimEnergyVector.size();i++){
        FileOutput<<std::exp(theLogPrimEnergyVector[i])/MeV<<'\t'<<std::exp(theLogCrossSectionVector[i])<<G4endl;
        size_t j1=0;
        FileOutput<<theLogSecondEnergyMatrix[i]->size()<<G4endl;
        for (size_t j=0;j<theLogSecondEnergyMatrix[i]->size();j++){
            FileOutput<<std::exp((*theLogSecondEnergyMatrix[i])[j]);
            j1++;
            if (j1<10) FileOutput<<'\t';
            else  {
                FileOutput<<G4endl;
                j1=0;
            }   
        }
        if (j1>0) FileOutput<<G4endl;
        j1=0;
        FileOutput<<theLogProbMatrix[i]->size()<<G4endl;
        for (size_t j=0;j<theLogProbMatrix[i]->size();j++){
            FileOutput<<std::exp((*theLogProbMatrix[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