FPDistribution.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 _FPDistribution_
#define _FPDistribution_
 
#include<fstream>
#include<string>
#include<vector>
#include<map>
 
#include "ZAI.hxx"
#include "FissionProduct.hxx"
 
class FissionProduct;
class ZAI;
 
using namespace std;
 
 
class FPRecord
{
  public:
 
    FPRecord(unsigned Rec, short NE);
 
    ~FPRecord();
    void Add(double E, unsigned Rec);
 
    int GetNE()
    {
        return fNE;   
    }
    double GetE(int i)
    {
        return fE[i];   
    }
    unsigned GetERecord(int i)
    {
        return fEStartRec[i];   
    }
    unsigned GetStartRecord()
    {
        return fStartRecord;   
    }
 
  private:
    unsigned            fStartRecord; 
    short int           fNE; 
    vector<double>      fE; 
    vector<unsigned>    fEStartRec; 
};
 
 
class FPDistribution
{
  public:
 
    FPDistribution(int Z, int A, int I = 0);
    ~FPDistribution(); 
 
    void    Add(int Z, int A, int I, double Yield);
    void    Add(FissionProduct FP); 
    vector<double> GetEnergyGroups();
    double  Yield(int i)
    {
        return fYield[i];   
    }
    ZAI    *GetZAI(int i)
    {
        return fZAI[i];   
    }
    int     GetNFP()
    {
        return fNFP;   
    }
    int     GetNZAI()
    {
        return fZAI.size();   
    }
    void ReadYield(double E); 
    void ReadYield(int Ebin);
    void FindClosestZA(); 
    bool FindZA(int Z, int A, int I, int SF); 
    void SetSpontaneousFission(int SF)
    {
        fSF = SF;   
    }
    void SetASCIIFileName(string name = "FPavailable.dat")
    {
        fASCIIFileName = name;   
    }
    void SetBinaryFileName(string name = "FPyield.bin")
    {
        fBinaryFileName = name;   
    }
    void BuildMap(); 
 
  private:
    void ReadYieldFromRecord(unsigned rec);
    int FindZAIinVector(int aZ, int aA, int aI);
    void DeleteYields();
    int                         fZ; 
    int                         fA; 
    int                         fI; 
    vector<ZAI *>               fZAI; 
    vector<double>              fYield; 
    vector<double>              fE; 
    int                         fSF; 
    int                         fNFP; 
    string                      fASCIIFileName; 
    string                      fBinaryFileName; 
    map<int, map<int, map<int, map<int, FPRecord *> > > > fAvaillableFP; 
    bool                        fIsMapBuild;
 
};
#endif