TNucleus.h

Go to the documentation of this file.

#ifndef TNUCLEUS_H
#define TNUCLEUS_H

#include <iostream>
#include <math.h>
#include <string.h>
#include <iomanip>
#include <stdlib.h>
#include <fstream>
#include <string>


#include "TTransition.h"

#include "TObject.h"
#include "TNamed.h"
#include "TList.h"

class TNucleus : public TNamed{

 private:
  static const char *grsipath;
  static std::string massfile; //The massfile to be used, which includes Z, N, atomic symbol, and mass excess
  //static const char *massfile; //The massfile to be used, which includes Z, N, atomic symbol, and mass excess
  //static std::string masspath;

 public:
  TNucleus() { };                                                         // Should not be use, here so we can write things to a root file.
  TNucleus(const char* symbol);                                           // Creates a nucleus based on symbol and sets all parameters from mass.dat
  TNucleus(int Z, int N, double mass, const char* symbol);                // Creates a nucleus with Z, N, mass, and symbol
  TNucleus(int Z, int N, const char* MassFile = 0);                      // Creates a nucleus with Z, N using mass table (default MassFile = "mass.dat")

  virtual ~TNucleus();
        
  //static void SetMassFile(const char *tmp = NULL);// {massfile = tmp;} //Sets the mass file to be used

  static const char* SortName(const char *name);
  //void SetName(const char *name) { fName = name; }
  void SetZ(int);                                       // Sets the Z (# of protons) of the nucleus
  void SetN(int);                                       // Sets the N (# of neutrons) of the nucleus
  void SetMassExcess(double);                           // Sets the mass excess of the nucleus (in MeV)
  void SetMass(double);                                 // Sets the mass manually (in MeV)
  void SetMass();                                       // Sets the mass based on the A and mass excess of nucleus (in MeV)
  void SetSymbol(const char*);                          // Sets the atomic symbol for the nucleus


  // void AddTransition(Double_t energy, Double_t intensity, Double_t energy_uncertainty = 0.0, Double_t intensity_uncertainty = 0.0);
  // void AddTransition(TTransition *tran);
  //Bool_t RemoveTransition(Int_t idx);
  //TGRSITransition *GetTransition(Int_t idx);

  //const char* GetName()        { return fName.c_str(); }
  int GetZ() const             { return fZ; }                   // Gets the Z (# of protons) of the nucleus
  int GetN() const             { return fN; }                   // Gets the N (# of neutrons) of the nucleus
  int GetA() const             { return fN+fZ; }                // Gets the A (Z + N) of the nucleus
  double GetMassExcess() const { return fMassExcess; }          // Gets the mass excess of the nucleus (in MeV)
  double GetMass() const       { return fMass; }                // Gets the mass of the nucleus (in MeV)
  const char* GetSymbol() const{ return fSymbol.c_str(); }      // Gets the atomic symbol of the nucleus


  void AddTransition(Double_t energy, Double_t intensity, Double_t energy_uncertainty = 0.0, Double_t intensity_uncertainty = 0.0);
  void AddTransition(TTransition *tran);
  //Bool_t RemoveTransition(Int_t idx);
  TTransition *GetTransition(Int_t idx);
  
  Int_t NTransitions() const { return TransitionList.GetSize();};
  double GetRadius() const;                                     
  int GetZfromSymbol(char*);                            

  TList TransitionList;
  //bool SetSourceData();

  void Print(Option_t *opt = "") const;
  void WriteSourceFile(std::string outfilename = "");



private:
  int fA;                                               // Number of nucleons (Z + N) 
  int fN;                                               // Number of neutrons (N)
  int fZ;                                               // Number of protons (Z)
  double fMass;                                         // Mass (in MeV)
  double fMassExcess;                                   // Mass excess (in MeV)
  std::string fSymbol;                                  // Atomic symbol (ex. Ba, C, O, N)
  //std::string fName;                                    // Name, whatever user inputs (ex. 9C, 9c);

  //JAB
  bool LoadTransitionFile();

        ClassDef(TNucleus,1);                           // Creates a nucleus with corresponding nuclear information
};



#endif