BetaSpectra.hxx

Go to the documentation of this file.

/*
  This file is part of MURE,
  Copyright (C) 2007-2021 MURE developers.
 
  MURE is free software: you can redistribute it and/or modify
  it under the terms of the GNU Lesser General Public License as published by
  the Free Software Foundation, either version 3 of the License, or
  (at your option) any later version.
 
  MURE is distributed in the hope that it will be useful,
  but WITHOUT ANY WARRANTY; without even the implied warranty of
  MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
  GNU Lesser General Public License for more details.
 
  You should have received a copy of the GNU Lesser General Public License
  along with MURE.  If not, see <http://www.gnu.org/licenses/>.
*/
 
#ifndef _BetaSpectra_
#define _BetaSpectra_
 
// GetBetaSpectrum(35, 90, 0, "JENDL", 100);
#include <vector>
#include <iostream>
#include <string>
#include <map>
 
#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