source.h

#ifndef PYNE_IQ4M73STINHJDPRV6KWUZZXOYE
#define PYNE_IQ4M73STINHJDPRV6KWUZZXOYE

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

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

namespace pyne {

class Source {
 public:
  Source();  

  ~Source();  

  // mcnp tally
  virtual std::string mcnp(int version=5) const {return "";}

  // fluka tally
  virtual std::string fluka() {}
};

class PointSource : public Source {
 public:
  PointSource(double _x = 0, double _y = 0, double _z = 0, double _u = 0,
              double _v = 0, double _w = 0, double _E = 14, std::string _particle = "Neutron",
              double _weight = 1);  

  ~PointSource();  

  // mcnp tally
  virtual std::string mcnp(int version=5) const;

  // fluka tally
  std::string fluka();



  double x;
  double y;
  double z;

  double u;
  double v;
  double w;

  double E;
  double weight;
  std::string particle;
};

std::ostream& operator<<(std::ostream& os, Source source);

// End pyne namespace
}  // namespace pyne

#endif