RIPLData.h

Go to the documentation of this file.

#ifndef _RIPL_DATA_
#define _RIPL_DATA_
 
#include <vector>
#include <map>
#include <string>
#include <cstdint>
class riplLevel
{
public:
  riplLevel() {;}
  void prepareLevel();
  uint16_t getDecayLevel(double a_randNum);
  void setEnergy(double a_energy) {m_energy = a_energy;}
  void setENSDFLevel(int a_ensdfLevel) {m_ensdfLevel = a_ensdfLevel;}
  int getENSDFLevel() {return m_ensdfLevel;}
  double getEnergy() {return m_energy;};
  // void addLevelID(uint16_t a_levelIDs) {m_nextLevelID.push_back(a_levelIDs);}
  void setLevelID(std::vector<uint16_t> a_levelIDS) {m_nextLevelID = a_levelIDS;}
  std::vector<uint16_t> getLevelIDs() {return m_nextLevelID;}
  // void addBR(double a_BR) {m_branchRatio.push_back(a_BR);}
  void setBR(std::vector<double> a_BR) {m_branchRatio = a_BR;}
  void setSpin(float a_spin) {m_spin = a_spin;}
  float getSpin() {return m_spin;}
  void setParity(int a_parity) {m_parity = a_parity;}
  int getParity() {return m_parity;}
  int getNumDecayPaths();
private:
 
  int m_ensdfLevel = 0;
  double m_energy=0;
  float m_spin =-1;
  int m_parity = 0;
  std::vector<uint16_t> m_nextLevelID;
  std::vector<double> m_branchRatio;
};
 
class RIPL
{
public:
  RIPL() {;}
  int readRIPL(int a_ZID, 
                int a_AID,
                double BR_cutoff);
  double getLevelEnergy(uint16_t a_level);
  uint16_t getLevelNumber(int a_level);
  float getLevelSpin(uint16_t a_level);
  int getLevelParity(uint16_t a_level);
  void setNumLevels(int a_numLevels) {m_numLevels = a_numLevels;}
  std::vector<double> doCascade(uint16_t a_initLevel);
 
  std::map<uint16_t, riplLevel> getLevelMap();
  std::vector<uint16_t> getTransitionLevelIDs(uint16_t a_initLevel);
  int getNumLevels();
private:
  std::map<uint16_t, riplLevel> m_data;
  int m_numLevels =0;
};
 
#endif