customnucleus.h

Go to the documentation of this file.

/*
 *    Copyright (c) 2014-2019
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 */
#ifndef SRC_INCLUDE_CUSTOMNUCLEUS_H_
#define SRC_INCLUDE_CUSTOMNUCLEUS_H_

#include <fstream>
#include <map>
#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,
                                      int particle_number) const;
  std::string file_path(const std::string& file_directory,
                        const std::string& file_name);
  void generate_fermi_momenta() override;
  int number_of_nucleons_ = 0;
  std::vector<Nucleoncustom> custom_nucleus_;
  size_t index = 0;

 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_;
};

}  // namespace smash

#endif  // SRC_INCLUDE_CUSTOMNUCLEUS_H_