deformednucleus.h

Go to the documentation of this file.

/*
 *    Copyright (c) 2014-2022,2024-2025
 *      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,
      SpinInteractionType spin_interaction_type = SpinInteractionType::Off);
  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);
 
  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 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;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_DEFORMEDNUCLEUS_H_