BetaSpectra.hxx

Go to the documentation of this file.

#ifndef _BetaSpectra_
#define _BetaSpectra_ 

// GetBetaSpectrum(35,90,0,"JENDL",100);
#include <math.h>
#include <vector>
#include <iostream>
#include <string>

#include "MureHeaders.hxx"
using namespace std;


enum BaseNumber{UNDEFINED=-1, BDB, BDBN, ENDF, JEFF, JENDL, RUDSTAM};

class BetaSpectra
{
public:
    BetaSpectra(); //constructor
    ~BetaSpectra(); //  destructor
    void ReadAllBetaData();
    void ReadBDB();
    void ReadBDBN();
    void ReadENDFB68();
    void ReadJEFF31();
    void ReadJENDL();
    void ReadRUDSTAM();
    vector <double> GetBetaSpectrum(int Z,int A,int I=0, BaseNumber BN=UNDEFINED); // don't forget methods to read end points
    //vector <double> GetBetaSpectrum(int Z,int A,int I=0, BaseNumber BN, int BinSize=100) {RebinSpectrum(GetBetaSpectrum(Z,A,I,BN), BinSize);}
    vector <double> RebinSpectrum(vector <double> Spectrum, int BinSize);
    //vector<double> AddSpectra(vector <double> Spectrum1,vector <double> Spectrum2);
    //vector<double> AddSpectra(BetaCollector* BC); ????

// il faut coupler ca au nombre de noyaux maintenant
vector<double> GetNbFP(int Z, int A, int I=0);
void ComputeNbFP();

//vector<double> AddAllSpectra(); 

//on fait 2 possibilites ? celle d'avoir les spectres construits dans et en dehors de MURE ?



protected:
    bool fDataIsRead;
    string fBDBname;
    string fBDBNname;
    string fENDFB68name;
    string fJEFF31name;
    string fJENDLname;
    string fRUDSTAMname;
    map<int ,map<int ,map<int ,map<int,map<int, double > > > > > fBetaSpectra; 
    // Z         A         I    Basenum   SpecVal
    map<int ,map<int ,map<int ,map<int, double > > > > fBetaEndPoint; 
    map<int, int> fBinSize; 

        map<int ,map<int ,map<int ,map<int, double > > > > fNbFP; 
    // Z         A         I    time   NbVal


};
#endif