ACEIsotope.h

Go to the documentation of this file.

#ifndef _ACE_ISOTOPE_
#define _ACE_ISOTOPE_
#include "RIPLData.h"
// c++ includes
#include<iostream>
#include<iomanip>
#include<fstream>
#include<sstream>
#include<cmath>
#include<string>
#include<vector>
#include<map>
#include<algorithm>
using namespace std;
 
class ACEIsotope
{
public:
  ACEIsotope();
  ACEIsotope(SampleACE* a_samp);
  int loadACEData(string a_ZAID, int a_verbose=0);
private:
  void getNuBarData(const vector<double>& rawData);
  void getDelayedNData(const vector<double>& rawData);
  void getDelayedPreNData(const vector<double>& rawData);
  void readOutgoingSpec(const vector<double>& rawData, string a_type);
  int getLawLength(const vector<double>& data, int a_JED, int a_IDAT, int a_LAW);
 
public:
  int drawXSec(int a_MT, bool a_writeToFile=false);
 
  angDist getAngDists(int a_MT);
  secondarySpectra getOutgoingSpec(int a_MT);
 
  void reset();
 
  int sampleMT(double a_energyIn);
  vector<outgoingData> sampleOutgoing(double a_energyIn, 
                                      int a_MT,
                                      bool a_sampleGamma=true);
 
  map<int, ACECrossSection> getMTXSecs();
 
  void getDrawableMTs();
 
  double getTotalXSec(double a_energyIn);
private:
 
  int m_verbose=0;
 
  map<string, int> m_nxsFlags;
  map<string, int> m_jxsFlags;
 
  map<int, ACECrossSection> m_MTxsec;
  vector<int> m_nMTsProd;
  vector<int> m_gMTsProd;
 
  map<int, int> m_MTMult;
 
  RIPL m_riplData;
 
  vector<double> m_energy;
  vector<double> m_totalXSec;
  vector<double> m_totalRemoval;
  vector<double> m_totalElas;
  vector<double> m_totalHeat;
  vector<double> m_nu_coeffsTotal;
  vector<double> m_nu_coeffsPrompt;
  map<int, tabData> m_pt_nuBar;
  tabData m_d_nuBar;
  map<int, pair<double, tabData>> m_delayNeutronPrecusor;
  map<int, pair<int,vector<double>>> m_sigsMT;
  map< int, angDist> m_angDistsMT;
  map<int, secondarySpectra> m_secondSpecNMT;
 
  SampleACE* m_samp;
 
};
 
#endif