tally.h

#ifndef PYNE_IQ4M73STINHJDPRV6KWUZZXOYE
#define PYNE_IQ4M73STINHJDPRV6KWUZZXOYE

#include <iostream>
#include <fstream>
#include <string>
#include <map>
#include <set>
#include <stdio.h>
#include <stdlib.h>

#ifndef PYNE_IS_AMALGAMATED
  #include "h5wrap.h"
  #include "utils.h"
#endif


namespace pyne
{
  class Tally
  {
  public:
    std::map<std::string, std::string> rx2fluka;
    std::map<std::string, std::string> rx2mcnp5;
    std::map<std::string, std::string> rx2mcnp6;

    Tally (); 

    Tally(std::string type, std::string particle_name, int entity,
      std::string entity_type, std::string entity_name,
      std::string tally_name = "", double entity_size = 0.0, 
      double normalization = 1.0);

    ~Tally (); 


    // Create hdf5 datatable for tallies
    hid_t create_dataspace(hid_t file, std::string datapath);

    // create filetype of data
    hid_t create_filetype();

    // create memory type for tally struct
    hid_t create_memtype();
    
    void from_hdf5(char * filename, char *datapath, int row = -1);

    void from_hdf5(std::string filename, std::string datapath, int row = -1);

    void write_hdf5( char * filename, char * datapath);

    void write_hdf5(std::string filename, std::string datapath);

    // mcnp tally
    std::string mcnp(int tally_index = 1, std::string mcnp_version = "mcnp5" );
    
    // fluka tally 
    std::string fluka(std::string unit_number = "-21");

    std::string entity_type; 
    std::string entity_name; 
    std::string particle_name; 
    std::string tally_type; 
    std::string tally_name; 
    int entity_id; 
    double entity_size; 
    double normalization; 
  };

  std::ostream& operator<< (std::ostream& os, Tally tally);


  typedef struct tally_struct {
    int entity_id;
    int entity_type;
    int tally_type;
    const char * particle_name;
    const char * entity_name;
    const char * tally_name;
    double entity_size;
    double normalization;
  } tally_struct;
  
// End pyne namespace
}

#endif