BetaSpectrum.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 _BetaSpectrum_
#define _BetaSpectrum_
 
#include <complex>
#include <cmath>
#include <vector>
#include <iostream>
#include <string>
 
using namespace std;
 
#include "Spectrum.hxx"
 
class Material;
 
 
class BetaSpectrum : public Spectrum
{
  public:
 
 
    BetaSpectrum(double Emin, double Emax, int NumberOfBins, bool log = false);
 
 
    BetaSpectrum(vector < double > LowerEnergy);
 
 
    BetaSpectrum(double *LowerEnergy, int NumberOfBins);
 
 
    BetaSpectrum(string BinningFileName);
 
    BetaSpectrum(int NumberOfBins, double a, double b, double c);
 
    BetaSpectrum(const BetaSpectrum  &spect);
 
 
    BetaSpectrum *ReadENSDF(int Z, int A, int I, double NAtoms = 1);
 
  protected:
 
    void ExtractDataFromLine(string Line, Material *TheParent, double NAtoms);
 
    double CalculateIntensity(int Zdaughter, double EndPointEnergy, double Energy);
    void GetSpinParity(string sSpinParity, double  &Spin, char  &Parity);
 
 
    complex < double > cgamma(const complex < double > &z, bool LOG = true) ;
 
    bool fAllowToAdd;           
    double fDecayConst;         
    double fRelativeIntensity;  
    double fBranchingRatio;     
    double fNB;         
    double fEofDecayLevel;      
    double fEnergyLevel;        
    double fEndPointEnergy;     
    double fQvalue;         
    double fIntensity;      
    bool fIsomStateIsOk;        
    char fParityParent;     
    double fSpinParent;     
    char fParityDaughter;       
    double fSpinDaughter;       
};
#endif