72 std::vector< double>* aLogProbVector,
size_t n_pro_decade){
82 std::vector< size_t>* aLogProbVectorIndex = 0;
85 if (n_pro_decade > 0 && aLogProbVector->size()>0) {
86 aLogProbVectorIndex =
new std::vector< size_t>();
87 dlog=std::log(10.)/n_pro_decade;
93 aLogProbVectorIndex->push_back(theInterpolator->
FindPosition(log_val,(*aLogProbVector)));
110 std::vector< double>*& aLogProbVector, std::vector< size_t>*& aLogProbVectorIndex)
125 { std::fstream FileOutput(file_name, std::ios::out);
126 FileOutput<<std::setiosflags(std::ios::scientific);
127 FileOutput<<std::setprecision(6);
136 if (j1<10) FileOutput<<
'\t';
142 if (j1>0) FileOutput<<
G4endl;
148 if (j1<10) FileOutput<<
'\t';
154 if (j1>0) FileOutput<<
G4endl;
172 for (
size_t i=0; i<n1;i++){
181 for (
size_t j=0; j<n2;j++){
187 for (
size_t j=0; j<n2;j++){
std::vector< std::vector< size_t > * > theLogProbMatrixIndex
static constexpr double MeV
G4AdjointCSMatrix(G4bool aBool)
std::vector< std::vector< double > * > theLogSecondEnergyMatrix
void AddData(G4double aPrimEnergy, G4double aCS, std::vector< double > *aLogSecondEnergyVector, std::vector< double > *aLogProbVector, size_t n_pro_decade=0)
size_t FindPosition(G4double &x, std::vector< G4double > &x_vec, size_t ind_min=0, size_t ind_max=0)
std::vector< double > theLogPrimEnergyVector
void Write(G4String file_name)
typedef int(XMLCALL *XML_NotStandaloneHandler)(void *userData)
std::vector< double > theLogCrossSectionVector
static G4AdjointInterpolator * GetInstance()
std::vector< std::vector< double > * > theLogProbMatrix
unsigned int nb_of_PrimEnergy
G4bool GetData(unsigned int i, G4double &aPrimEnergy, G4double &aCS, G4double &log0, std::vector< double > *&aLogSecondEnergyVector, std::vector< double > *&aLogProbVector, std::vector< size_t > *&aLogProbVectorIndex)
void Read(G4String file_name)
G4bool is_scat_proj_to_proj_case
std::vector< double > log0Vector
T min(const T t1, const T t2)
brief Return the smallest of the two arguments