particle.h

// defines the primary particle types that are allowed by most monte carlo codes
// some monte carlo codes allow us to score so called "heavy ions", in fact we
// define heavy ions to be particles with more than one neutron or proton

#ifndef PYNE_UWPOP4EE6BEB5CZK4BQQHYTCEI
#define PYNE_UWPOP4EE6BEB5CZK4BQQHYTCEI

#include <string>
#include <map>
#include <set>

#ifndef PYNE_IS_AMALGAMATED
  #include "nucname.h"
#endif

#define NUM_PARTICLES 32

namespace pyne
{
namespace particle
{
  extern int _pdcids[NUM_PARTICLES]; 
  extern std::string _docs[NUM_PARTICLES];
  extern std::string _names[NUM_PARTICLES];
  extern std::set<std::string> names;
  extern std::set<int> pdc_nums;
  extern std::map<std::string,int> name_id;
  extern std::map<int,std::string> id_name;
  extern std::map<std::string,std::string> docs;
  extern std::map<std::string,int> altnames;
  extern std::map<std::string,std::string> part_to_mcnp;
  extern std::map<std::string,std::string> part_to_mcnp6;
  extern std::map<std::string,std::string> part_to_fluka;
  extern std::map<std::string,std::string> part_to_geant4;


  bool is_hydrogen(int n);
  bool is_hydrogen(char *s);
  bool is_hydrogen(std::string s);

  bool is_heavy_ion(int s);
  bool is_heavy_ion(char *s);
  bool is_heavy_ion(std::string s);

  bool is_valid(int n);
  bool is_valid(char *s);
  bool is_valid(std::string s);

  int id(int n);
  int id(char *s);
  int id(std::string s);

  std::string name(int s);
  std::string name(char *s);
  std::string name(std::string s);

  std::string mcnp(int s);
  std::string mcnp(char *s);
  std::string mcnp(std::string s);

  std::string mcnp6(int s);
  std::string mcnp6(char *s);
  std::string mcnp6(std::string s);

  std::string fluka(int s);
  std::string fluka(char *s);
  std::string fluka(std::string s);

  std::string geant4(int s);
  std::string geant4(char *s);
  std::string geant4(std::string s);


  std::string describe(int s);
  std::string describe(char *s);
  std::string describe(std::string s);

  void * _fill_maps();
  extern void * filler;  


  class NotAParticle : public std::exception
  {
  public:
    NotAParticle () {};

    ~NotAParticle () throw () {};

    NotAParticle(std::string particle_name)
    {
       part_name = particle_name;
    }

    virtual const char* what() const throw()
    {
      std::string pname ("Not a valid particle name ");
      if(!part_name.empty())
        pname += part_name;
      const char *pname_rtn = pname.c_str();
      return pname_rtn;
    }

    private:
       std::string part_name;  


  };
}
}

#endif