deformednucleus.h

Go to the documentation of this file.

/*
 *    Copyright (c) 2014-2022
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 */
#ifndef SRC_INCLUDE_SMASH_DEFORMEDNUCLEUS_H_
#define SRC_INCLUDE_SMASH_DEFORMEDNUCLEUS_H_
 
#include <map>
 
#include "angles.h"
#include "configuration.h"
#include "forwarddeclarations.h"
#include "nucleus.h"
#include "threevector.h"
 
namespace smash {
 
double y_l_m(int l, int m, double cosx, double phi);
 
class DeformedNucleus : public Nucleus {
 public:
  DeformedNucleus(const std::map<PdgCode, int> &particle_list, int nTest);
  DeformedNucleus(Configuration &config, int nTest, bool auto_deformation);
  ThreeVector distribute_nucleon() override;
 
  void set_deformation_parameters_automatic();
 
  void set_deformation_parameters_from_config(Configuration &config);
 
  void set_orientation_from_config(Configuration &orientation_config);
  void rotate() override;
 
  inline double get_saturation_density() const { return saturation_density_; }
  double nucleon_density(double r, double cosx, double phi) const override;
  double nucleon_density_unnormalized(double r, double cosx,
                                      double phi) const override;
  double integrant_nucleon_density_phi(double r, double cosx) const;
  double calculate_saturation_density() const override;
  inline void set_beta_2(double b2) { beta2_ = b2; }
  inline void set_gamma(double ga) { gamma_ = ga; }
  inline void set_beta_3(double b3) { beta3_ = b3; }
  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);
  }
  inline void set_angle_psi(double psi) { nuclear_orientation_.set_psi(psi); }
  inline double get_beta2() { return beta2_; }
  inline double get_beta3() { return beta3_; }
  inline double get_beta4() { return beta4_; }
 
 private:
  double beta2_ = 0.0;
  double gamma_ = 0.0;
  double beta3_ = 0.0;
  double beta4_ = 0.0;
  Angles nuclear_orientation_;
  bool random_rotation_ = false;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_DEFORMEDNUCLEUS_H_