customnucleus.h

Go to the documentation of this file.

/*
 *    Copyright (c) 2014-2020
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 */
#ifndef SRC_INCLUDE_SMASH_CUSTOMNUCLEUS_H_
#define SRC_INCLUDE_SMASH_CUSTOMNUCLEUS_H_
 
#include <fstream>
#include <map>
#include <memory>
#include <string>
#include <vector>
 
#include "nucleus.h"
#include "pdgcode.h"
#include "threevector.h"
namespace smash {
 
struct Nucleoncustom {
  double x;
  double y;
  double z;
  bool spinprojection;
  bool isospin;
};
 
class CustomNucleus : public Nucleus {
 public:
  CustomNucleus(Configuration& config, int testparticles, bool same_file);
  void fill_from_list(const std::vector<Nucleoncustom>& vec);
  ThreeVector distribute_nucleon() override;
  void arrange_nucleons() override;
  std::vector<Nucleoncustom> readfile(std::ifstream& infile) const;
  std::string file_path(const std::string& file_directory,
                        const std::string& file_name);
  void generate_fermi_momenta() override;
 
 private:
  /*
   * The unique_ptr is only required to work around a bug in GCC 4.8, because it
   * seems to be trying to use the non-existing copy-constructor of
   * `std::ifstream`. Newer compilers don't require this unneccessary
   * indirection.
   */
  static std::unique_ptr<std::ifstream> filestream_shared_;
  std::unique_ptr<std::ifstream> filestream_;
  std::unique_ptr<std::ifstream>* used_filestream_;
  int number_of_nucleons_ = 0;
  int number_of_protons_ = 0;
  int number_of_neutrons_ = 0;
  std::vector<Nucleoncustom> custom_nucleus_;
  size_t index_ = 0;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_CUSTOMNUCLEUS_H_