deformednucleus.h

Go to the documentation of this file.

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

#include <map>

#include "angles.h"
#include "configuration.h"
#include "forwarddeclarations.h"
#include "nucleus.h"
#include "threevector.h"

namespace smash {

class DeformedNucleus : public Nucleus {
 public:
  DeformedNucleus(const std::map<PdgCode, int> &particle_list, int nTest);
  DeformedNucleus(Configuration &config, int nTest);

  double deformed_woods_saxon(double r, double cosx) const;

  ThreeVector distribute_nucleon() const override;

  void set_deformation_parameters_automatic();

  void set_deformation_parameters_from_config(Configuration &config);

  void rotate() override;

  void generate_fermi_momenta() override;

  double y_l_0(int l, double cosx) const;

  inline void set_beta_2(double b2) { beta2_ = b2; }
  inline void set_beta_4(double b4) { beta4_ = b4; }
  inline void set_polar_angle(double theta) {
    nuclear_orientation_.set_theta(theta);
  }
  inline void set_azimuthal_angle(double phi) {
    nuclear_orientation_.set_phi(phi);
  }

 private:
  double beta2_ = 0.0;
  double beta4_ = 0.0;
  Angles nuclear_orientation_;
};

}  // namespace smash

#endif  // SRC_INCLUDE_DEFORMEDNUCLEUS_H_