#include <deformednucleus.h>

DeformedNucleus: Child of nucleus for deformed nuclei.

All options from the nucleus will still apply. The deformed nucleus adds new or updated features which are outlined below.

Definition at line 38 of file deformednucleus.h.

Inheritance diagram for smash::DeformedNucleus:

Public Member Functions
	DeformedNucleus (const std::map< PdgCode, int > &particle_list, int nTest)
	Constructor for DeformedNucles which takes a particle list and the number of testparticles. More...

	DeformedNucleus (Configuration &config, int nTest, bool auto_deformation)
	Constructor for DeformedNucleus, that needs the configuration parameters from the inputfile and the number of testparticles. More...

ThreeVector	distribute_nucleon () override
	Deformed Woods-Saxon sampling routine. More...

void	set_deformation_parameters_automatic ()
	Sets the deformation parameters of the radius according to the current mass number. More...

void	set_deformation_parameters_from_config (Configuration &config)
	Set parameters for spherical deformation of the nucleus from the values specified in the configuration file. More...

void	set_orientation_from_config (Configuration &orientation_config)
	Set angles for orientation of nucleus from config file. More...

void	rotate () override
	Rotates the nucleus according to members nucleus_polar_angle_ and nucleus_azimuthal_angle_ and updates nucleon positions. More...

double	get_saturation_density () const

double	nucleon_density (double r, double cosx, double phi) const override
	Return the deformed Woods-Saxon probability density for the given position. More...

double	nucleon_density_unnormalized (double r, double cosx, double phi) const override
	Return the unnormalized deformed Woods-Saxon distribution for the given position. More...

double	integrant_nucleon_density_phi (double r, double cosx) const
	Return the integral over the azimuthal angle phi. More...

double	calculate_saturation_density () const override

void	set_beta_2 (double b2)
	Set deformation coefficient for Y_2_0. More...

void	set_gamma (double ga)
	Set the triaxiality coefficient gamma for Y_2_0 and Y_2_2. More...

void	set_beta_3 (double b3)
	Set deformation coefficient for Y_3_0. More...

void	set_beta_4 (double b4)
	Set deformation coefficient for Y_4_0. More...

void	set_polar_angle (double theta)
	Set the nucleus polar angle. More...

void	set_azimuthal_angle (double phi)
	Set the nucleus azimuthal angle. More...

void	set_angle_psi (double psi)
	Set the angle psi. More...

double	get_beta2 ()
	return the beta2 value. More...

double	get_beta3 ()
	return the beta3 value. More...

double	get_beta4 ()
	return the beta4 value. More...

Public Member Functions inherited from smash::Nucleus
	Nucleus ()=default
	default constructor More...

	Nucleus (Configuration &config, int nTest)
	Constructor for Nucleus, that needs the configuration parameters from the inputfile and the number of testparticles. More...

	Nucleus (const std::map< PdgCode, int > &particle_list, int nTest)
	Constructor which directly initializes the Nucleus with particles and respective counts. More...

virtual	~Nucleus ()=default

double	mass () const

double	woods_saxon (double x)
	Woods-Saxon distribution. More...

virtual void	arrange_nucleons ()
	Sets the positions of the nucleons inside a nucleus. More...

virtual void	set_parameters_automatic ()
	Sets the deformation parameters of the Woods-Saxon distribution according to the current mass number. More...

virtual void	set_parameters_from_config (Configuration &config)
	Sets the parameters of the Woods-Saxon according to manually added values in the configuration file. More...

virtual void	generate_fermi_momenta ()
	Generates momenta according to Fermi motion for the nucleons. More...

void	boost (double beta_scalar)
	Boosts the nuclei into the computational frame, such that the nucleons have the appropriate momentum and the nuclei are lorentz-contracted. More...

void	fill_from_list (const std::map< PdgCode, int > &particle_list, int testparticles)
	Adds particles from a map PDG code => Number_of_particles_with_that_PDG_code to the nucleus. More...

void	shift (double z_offset, double x_offset, double simulation_time)
	Shifts the nucleus to correct impact parameter and z displacement. More...

void	copy_particles (Particles *particles)
	Copies the particles from this nucleus into the particle list. More...

size_t	size () const
	Number of numerical (=test-)particles in the nucleus: More...

size_t	number_of_particles () const
	Number of physical particles in the nucleus: More...

size_t	number_of_protons () const
	Number of physical protons in the nucleus: More...

FourVector	center () const
	Calculate geometrical center of the nucleus. More...

void	set_label (BelongsTo label)
	Sets target / projectile labels on nucleons. More...

void	align_center ()
	Shifts the nucleus so that its center is at (0,0,0) More...

virtual void	set_saturation_density (double density)
	Sets the saturation density of the nucleus. More...

std::vector< ParticleData >::iterator	begin ()
	For iterators over the particle list: More...

std::vector< ParticleData >::iterator	end ()
	For iterators over the particle list: More...

std::vector< ParticleData >::const_iterator	cbegin () const
	For const iterators over the particle list: More...

std::vector< ParticleData >::const_iterator	cend () const
	For const iterators over the particle list: More...

void	set_diffusiveness (double diffuse)
	Sets the diffusiveness of the nucleus. More...

double	get_diffusiveness () const

double	get_saturation_density () const

double	default_nuclear_radius ()
	Default nuclear radius calculated as: More...

void	set_nuclear_radius (double rad)
	Sets the nuclear radius. More...

double	get_nuclear_radius () const

Private Attributes
double	beta2_ = 0.0
	Deformation parameter for angular momentum l=2. More...

double	gamma_ = 0.0
	Triaxiality parameter for angular momentum l=2. More...

double	beta3_ = 0.0
	Deformation parameter for angular momentum l=3. More...

double	beta4_ = 0.0
	Deformation parameter for angular momentum l=4. More...

Angles	nuclear_orientation_
	Nucleus orientation (initial profile in xz plane) in terms of a pair of angles (theta, phi) More...

bool	random_rotation_ = false
	Whether the nuclei should be rotated randomly. More...

Additional Inherited Members
Protected Member Functions inherited from smash::Nucleus
void	random_euler_angles ()
	Randomly generate Euler angles. More...

Protected Attributes inherited from smash::Nucleus
std::vector< ParticleData >	particles_
	Particles associated with this nucleus. More...

double	saturation_density_ = nuclear_density
	Saturation density of this nucleus. More...

double	euler_phi_
	Euler angel phi. More...

double	euler_theta_
	Euler angel theta. More...

double	euler_psi_
	Euler angel psi. More...

Constructor & Destructor Documentation

◆ DeformedNucleus() [1/2]

smash::DeformedNucleus::DeformedNucleus	(	const std::map< PdgCode, int > &	particle_list,
		int	nTest
	)

Constructor for DeformedNucles which takes a particle list and the number of testparticles.

This constructor is only used for testing purposes.

Parameters

[in]	particle_list	Map with PDGCode and number of particles which make up the nucleus
[in]	nTest	number of testparticles

Definition at line 21 of file deformednucleus.cc.

23 : Nucleus(particle_list, nTest) {}

smash::Nucleus::Nucleus

Nucleus()=default

default constructor

◆ DeformedNucleus() [2/2]

smash::DeformedNucleus::DeformedNucleus	(	Configuration &	config,
		int	nTest,
		bool	auto_deformation
	)

Constructor for DeformedNucleus, that needs the configuration parameters from the inputfile and the number of testparticles.

Parameters

[in]	config	contains the parameters from the inputfile on the numbers of particles with a certain PDG code
[in]	nTest	number of testparticles
[in]	auto_deformation	whether or not deformation parameters should be set automatically

Definition at line 25 of file deformednucleus.cc.

     : Nucleus(config, nTest) {
   if (auto_deformation) {
     set_deformation_parameters_automatic();
     set_saturation_density(calculate_saturation_density());
   } else {
     set_deformation_parameters_from_config(config);
   }
   if (config.has_value({"Deformed", "Orientation"})) {
     Configuration sub_conf =
         config.extract_sub_configuration({"Deformed", "Orientation"});
     set_orientation_from_config(sub_conf);
   }
 }

Member Function Documentation

◆ distribute_nucleon()

ThreeVector smash::DeformedNucleus::distribute_nucleon ( )

overridevirtual

Deformed Woods-Saxon sampling routine.

Returns: Spatial position from uniformly sampling the deformed woods-saxon distribution

Reimplemented from smash::Nucleus.

Definition at line 41 of file deformednucleus.cc.

                                                 {
   double a_radius;
   Angles a_direction;
   // Set a sensible maximum bound for radial sampling.
   double radius_max =
       Nucleus::get_nuclear_radius() / Nucleus::get_diffusiveness() +
       Nucleus::get_nuclear_radius() * Nucleus::get_diffusiveness();
  
   // Sample the distribution.
   do {
     a_direction.distribute_isotropically();
     // sample r**2 dr
     a_radius = radius_max * std::cbrt(random::canonical());
   } while (random::canonical() > nucleon_density(a_radius,
                                                  a_direction.costheta(),
                                                  a_direction.phi()) /
                                      Nucleus::get_saturation_density());
  
   // Update (x, y, z) positions.
   return a_direction.threevec() * a_radius;
 }

◆ set_deformation_parameters_automatic()

void smash::DeformedNucleus::set_deformation_parameters_automatic ( )

Sets the deformation parameters of the radius according to the current mass number.

The deformation parameters are taken from Moller:1993ed [38]. Corrections to the deformation parameter beta2 in Uranium come from Kuhlman:2005ts [32]. For finite nucleon size corrections to the nuclear density and radius for copper and gold, see Hirano:2009ah [25], and Hirano:2010jg [27] for uranium.

Definition at line 63 of file deformednucleus.cc.

                                                            {
   // Set the deformation parameters
   // reference for U, Pb, Au, Cu: \iref{Moller:1993ed}
   // reference for Zr and Ru: \iref{Schenke:2019ruo}
   // reference for Xe: \iref{Moller:2015fba}
   bool listed = 0;
   const std::map<int, std::string> A_map = {{238, "Uranium"},
                                             {208, "Lead"},
                                             {197, "Gold"},
                                             {63, "Copper"},
                                             {129, "Xenon"}};
   const std::map<std::string, std::string> Z_map = {{"Uranium", "92"},
                                                     {"Lead", "82"},
                                                     {"Gold", "79"},
                                                     {"Copper", "29"},
                                                     {"Xenon", "54"}};
   int A = Nucleus::number_of_particles();
   int Z = Nucleus::number_of_protons();
   switch (A) {
     case 238:  // Uranium
       if (Z == 92) {
         set_beta_2(0.28);
         set_beta_4(0.093);
       } else {
         listed = true;
       }
       break;
     case 208:  // Lead
       if (Z == 82) {
         set_beta_2(0.0);
         set_beta_4(0.0);
       } else {
         listed = true;
       }
       break;
     case 197:  // Gold
       if (Z == 79) {
         set_beta_2(-0.131);
         set_beta_4(-0.031);
       } else {
         listed = true;
       }
       break;
     case 129:  // Xenon
       if (Z == 54) {
         set_beta_2(0.162);
         set_beta_4(-0.003);
       } else {
         listed = true;
       }
       break;
     case 63:  // Copper
       if (Z == 29) {
         set_beta_2(0.162);
         set_beta_4(-0.006);
       } else {
         listed = true;
       }
       break;
     case 96:
       if (Z == 40) {  // Zirconium
         set_beta_2(0.0);
         set_beta_4(0.0);
       } else if (Z == 44) {  // Ruthenium
         set_beta_2(0.158);
         set_beta_4(0.0);
       } else {
         throw std::domain_error(
             "Number of protons for nuclei with mass number A = 96 does not "
             "match that of Zirconium or Ruthenium. The deformation parameters "
             "for additional isobars are currently not implemented."
             " Please specify at least \"Beta_2\" and \"Beta_4\" "
             "manually and set \"Automatic: False.\" ");
       }
       break;
     default:
       throw std::domain_error(
           "Mass number not listed for automatically setting deformation "
           "parameters. Please specify at least \"Beta_2\" and \"Beta_4\" "
           "manually and set \"Automatic: False.\" ");
   }
   if (listed) {
     throw std::domain_error("Mass number is listed under " + A_map.at(A) +
                             " but the proton "
                             "number of " +
                             std::to_string(Z) +
                             " does not match "
                             "its " +
                             Z_map.at(A_map.at(A)) +
                             " protons."
                             "Please specify at least \"Beta_2\" and \"Beta_4\" "
                             "manually and set \"Automatic: False.\" ");
   }
 }

◆ set_deformation_parameters_from_config()

void smash::DeformedNucleus::set_deformation_parameters_from_config ( Configuration & config )

Set parameters for spherical deformation of the nucleus from the values specified in the configuration file.

Parameters

config The configuration for the deformation of this nucleus (projectile or target).

Definition at line 158 of file deformednucleus.cc.

                            {
   // Deformation parameters.
   if (config.has_value({"Deformed", "Beta_2"})) {
     set_beta_2(static_cast<double>(config.take({"Deformed", "Beta_2"})));
   }
   if (config.has_value({"Deformed", "Gamma"})) {
     set_gamma(static_cast<double>(config.take({"Deformed", "Gamma"})));
   }
   if (config.has_value({"Deformed", "Beta_3"})) {
     set_beta_3(static_cast<double>(config.take({"Deformed", "Beta_3"})));
   }
   if (config.has_value({"Deformed", "Beta_4"})) {
     set_beta_4(static_cast<double>(config.take({"Deformed", "Beta_4"})));
   }
 }

◆ set_orientation_from_config()

void smash::DeformedNucleus::set_orientation_from_config ( Configuration & orientation_config )

Set angles for orientation of nucleus from config file.

Parameters

[in] orientation_config The configuration for the deformation of this nucleus (projectile or target).

Definition at line 175 of file deformednucleus.cc.

                                        {
   // Read in orientation if provided, otherwise, the defaults are
   // theta = pi/2, phi = 0, as declared in the angles class
  
   if (orientation_config.has_value({"Theta"})) {
     if (orientation_config.has_value({"Random_Rotation"}) &&
         orientation_config.take({"Random_Rotation"})) {
       throw std::domain_error(
           "Random rotation of nuclei is activated although"
           " theta is provided. Please specify only either of them. ");
     } else {
       set_polar_angle(static_cast<double>(orientation_config.take({"Theta"})));
     }
   }
  
   if (orientation_config.has_value({"Phi"})) {
     if (orientation_config.has_value({"Random_Rotation"}) &&
         orientation_config.take({"Random_Rotation"})) {
       throw std::domain_error(
           "Random rotation of nuclei is activated although"
           " phi is provided. Please specify only either of them. ");
     } else {
       set_azimuthal_angle(
           static_cast<double>(orientation_config.take({"Phi"})));
     }
   }
  
   if (orientation_config.has_value({"Psi"})) {
     if (orientation_config.has_value({"Random_Rotation"}) &&
         orientation_config.take({"Random_Rotation"})) {
       throw std::domain_error(
           "Random rotation of nuclei is activated although"
           " psi is provided. Please specify only either of them. ");
     } else {
       set_angle_psi(static_cast<double>(orientation_config.take({"Psi"})));
     }
   }
  
   if (orientation_config.take({"Random_Rotation"}, false)) {
     random_rotation_ = true;
   }
 }

◆ rotate()

void smash::DeformedNucleus::rotate ( )

overridevirtual

Rotates the nucleus according to members nucleus_polar_angle_ and nucleus_azimuthal_angle_ and updates nucleon positions.

Reimplemented from smash::Nucleus.

Definition at line 219 of file deformednucleus.cc.

                              {
   if (random_rotation_) {
     // Randomly generate euler angles for theta and phi. Psi needs not be
     // assigned, as the nucleus objects are symmetric with respect to psi.
     Nucleus::random_euler_angles();
     set_azimuthal_angle(euler_phi_);
     set_polar_angle(euler_theta_);
     set_angle_psi(euler_psi_);
   }
   for (auto &particle : *this) {
     /* Rotate every vector by the nuclear azimuth phi, polar angle
      * theta and psi (the Euler angles). This means applying the matrix for a
      * rotation of phi about z, followed by the matrix for a rotation
      * theta about the rotated x axis. If the triaxiality coefficient is not
      * zero, one has to include the third rotation around psi as the nucleus is
      * not symmetric under rotation of any axis.*/
     ThreeVector three_pos = particle.position().threevec();
     three_pos.rotate(nuclear_orientation_.phi(), nuclear_orientation_.theta(),
                      nuclear_orientation_.psi());
     particle.set_3position(three_pos);
   }
 }

◆ get_saturation_density()

double smash::DeformedNucleus::get_saturation_density ( ) const

inline

Returns: the saturation density of the deformed_nucleus

See also: saturation_density_

Definition at line 102 of file deformednucleus.h.

102 { return saturation_density_; }

smash::Nucleus::saturation_density_

double saturation_density_

Saturation density of this nucleus.

Definition: nucleus.h:292

◆ nucleon_density()

double smash::DeformedNucleus::nucleon_density	(	double	r,
		double	cosx,
		double	phi
	)		const

overridevirtual

Return the deformed Woods-Saxon probability density for the given position.

This corresponds to the nuclear density at the very same position.

Parameters

[in]	r	The radius at which to sample
[in]	cosx	The cosine of the polar angle at which to sample
[in]	phi	The azimuthal angle at which to sample

Returns: The Woods-Saxon density

Reimplemented from smash::Nucleus.

Definition at line 260 of file deformednucleus.cc.

                                                           {
   return Nucleus::get_saturation_density() /
          (1 + std::exp((r - Nucleus::get_nuclear_radius() *
                                 (1 +
                                  beta2_ * (std::cos(gamma_) *
                                                y_l_m(2, 0, cosx, phi) +
                                            std::sqrt(2) * std::sin(gamma_) *
                                                y_l_m(2, 2, cosx, phi)) +
                                  beta3_ * y_l_m(3, 0, cosx, phi) +
                                  beta4_ * y_l_m(4, 0, cosx, phi))) /
                        Nucleus::get_diffusiveness()));
 }

◆ nucleon_density_unnormalized()

double smash::DeformedNucleus::nucleon_density_unnormalized	(	double	r,
		double	cosx,
		double	phi
	)		const

overridevirtual

Return the unnormalized deformed Woods-Saxon distribution for the given position.

Parameters

[in]	r	The radius
[in]	cosx	The cosine of the polar angle
[in]	phi	The azimuthal angle

Returns: The unnormalized Woods-Saxon distribution

Reimplemented from smash::Nucleus.

Definition at line 274 of file deformednucleus.cc.

                                                                        {
   return 1.0 /
          (1 + std::exp((r - Nucleus::get_nuclear_radius() *
                                 (1 +
                                  beta2_ * (std::cos(gamma_) *
                                                y_l_m(2, 0, cosx, phi) +
                                            std::sqrt(2) * std::sin(gamma_) *
                                                y_l_m(2, 2, cosx, phi)) +
                                  beta3_ * y_l_m(3, 0, cosx, phi) +
                                  beta4_ * y_l_m(4, 0, cosx, phi))) /
                        Nucleus::get_diffusiveness()));
 }

◆ integrant_nucleon_density_phi()

double smash::DeformedNucleus::integrant_nucleon_density_phi	(	double	r,
		double	cosx
	)		const

Return the integral over the azimuthal angle phi.

Parameters

[in]	r	The radius
[in]	cosx	The cosine of the polar angle

Returns: The unnormalized Woods-Saxon distribution integrated over dphi

Definition at line 288 of file deformednucleus.cc.

                                                                          {
   Integrator integrate;
   // Perform the phi integration. This is needed if the triaxiality coefficient
   // gamma is included, which includes a dependency around the phi axis.
   // Unfortunately the Integrator class does not support 3d integration which is
   // why this intermediate integral is needed. It has been checked that the
   // integral factorizes.
   const auto result = integrate(0.0, 2.0 * M_PI, [&](double phi) {
     return nucleon_density_unnormalized(r, cosx, phi);
   });
   return result.value();
 }

◆ calculate_saturation_density()

double smash::DeformedNucleus::calculate_saturation_density ( ) const

overridevirtual

Returns: the normalized ground state density for the corresponding Woods-Saxon parameter. This is done by integrating the Woods-Saxon distribution and setting the normalization such that the integral of the Woods-Saxon distribution yields the number of particles in the nucleus \(\int\rho(r)d^3r = N_{particles}\).

Reimplemented from smash::Nucleus.

Definition at line 302 of file deformednucleus.cc.

                                                            {
   Integrator2d integrate;
   // Transform integral from (0, oo) to (0, 1) via r = (1 - t) / t.
   // To prevent overflow, the integration is only performed to t = 0.01 which
   // corresponds to r = 99fm. Additionally the precision settings in the
   // Integrator2d scheme are equally important. However both these point affect
   // the result only after the seventh digit which should not be relevant here.
   if (gamma_ == 0.0) {
     const auto result = integrate(0.01, 1, -1, 1, [&](double t, double cosx) {
       const double r = (1 - t) / t;
       return twopi * std::pow(r, 2.0) *
              nucleon_density_unnormalized(r, cosx, 0.0) / std::pow(t, 2.0);
     });
     const auto rho0 = number_of_particles() / result.value();
     return rho0;
   } else {
     const auto result = integrate(0.01, 1, -1, 1, [&](double t, double cosx) {
       const double r = (1 - t) / t;
       return std::pow(r, 2.0) * integrant_nucleon_density_phi(r, cosx) /
              std::pow(t, 2.0);
     });
     const auto rho0 = number_of_particles() / result.value();
     return rho0;
   }
 }

◆ set_beta_2()

void smash::DeformedNucleus::set_beta_2 ( double b2 )

inline

Set deformation coefficient for Y_2_0.

Parameters

[in] b2 deformation coefficient for l=2

Definition at line 145 of file deformednucleus.h.

145 { beta2_ = b2; }

◆ set_gamma()

void smash::DeformedNucleus::set_gamma ( double ga )

inline

Set the triaxiality coefficient gamma for Y_2_0 and Y_2_2.

Parameters

[in] ga triaxiality coefficient for l=2

Definition at line 150 of file deformednucleus.h.

150 { gamma_ = ga; }

◆ set_beta_3()

void smash::DeformedNucleus::set_beta_3 ( double b3 )

inline

Set deformation coefficient for Y_3_0.

Parameters

[in] b3 deformation coefficient for l=3

Definition at line 155 of file deformednucleus.h.

155 { beta3_ = b3; }

◆ set_beta_4()

void smash::DeformedNucleus::set_beta_4 ( double b4 )

inline

Set deformation coefficient for Y_4_0.

Parameters

[in] b4 deformation coefficient for l=4

Definition at line 160 of file deformednucleus.h.

160 { beta4_ = b4; }

◆ set_polar_angle()

void smash::DeformedNucleus::set_polar_angle ( double theta )

inline

Set the nucleus polar angle.

Parameters

[in] theta Polar angle of position inside nucleus

Definition at line 165 of file deformednucleus.h.

                                             {
     nuclear_orientation_.set_theta(theta);
   }

◆ set_azimuthal_angle()

void smash::DeformedNucleus::set_azimuthal_angle ( double phi )

inline

Set the nucleus azimuthal angle.

Parameters

[in] phi Azimuthal angle of position inside nucleus

Definition at line 172 of file deformednucleus.h.

                                               {
     nuclear_orientation_.set_phi(phi);
   }

◆ set_angle_psi()

void smash::DeformedNucleus::set_angle_psi ( double psi )

inline

Set the angle psi.

Parameters

[in] psi Angle psi for properly rotating nucleus

Definition at line 179 of file deformednucleus.h.

179 { nuclear_orientation_.set_psi(psi); }

smash::Angles::set_psi

void set_psi(const double psi)

Sets the euler angle psi.

Definition: angles.h:213

◆ get_beta2()

double smash::DeformedNucleus::get_beta2 ( )

inline

return the beta2 value.

Definition at line 183 of file deformednucleus.h.

183 { return beta2_; }

◆ get_beta3()

double smash::DeformedNucleus::get_beta3 ( )

inline

return the beta3 value.

Definition at line 187 of file deformednucleus.h.

187 { return beta3_; }

◆ get_beta4()

double smash::DeformedNucleus::get_beta4 ( )

inline

return the beta4 value.

Definition at line 191 of file deformednucleus.h.

191 { return beta4_; }

Member Data Documentation

◆ beta2_

double smash::DeformedNucleus::beta2_ = 0.0

private

Deformation parameter for angular momentum l=2.

Definition at line 195 of file deformednucleus.h.

◆ gamma_

double smash::DeformedNucleus::gamma_ = 0.0

private

Triaxiality parameter for angular momentum l=2.

Definition at line 197 of file deformednucleus.h.

◆ beta3_

double smash::DeformedNucleus::beta3_ = 0.0

private

Deformation parameter for angular momentum l=3.

Definition at line 199 of file deformednucleus.h.

◆ beta4_

double smash::DeformedNucleus::beta4_ = 0.0

private

Deformation parameter for angular momentum l=4.

Definition at line 201 of file deformednucleus.h.

◆ nuclear_orientation_

Angles smash::DeformedNucleus::nuclear_orientation_

private

Nucleus orientation (initial profile in xz plane) in terms of a pair of angles (theta, phi)

Definition at line 206 of file deformednucleus.h.

◆ random_rotation_

bool smash::DeformedNucleus::random_rotation_ = false

private

Whether the nuclei should be rotated randomly.

Definition at line 210 of file deformednucleus.h.

The documentation for this class was generated from the following files:

src/include/smash/deformednucleus.h
src/deformednucleus.cc

Public Member Functions

Private Attributes

Additional Inherited Members

Constructor & Destructor Documentation

◆ DeformedNucleus() [1/2]

◆ DeformedNucleus() [2/2]

Member Function Documentation

◆ distribute_nucleon()

◆ set_deformation_parameters_automatic()

◆ set_deformation_parameters_from_config()

◆ set_orientation_from_config()

◆ rotate()

◆ get_saturation_density()

◆ nucleon_density()

◆ nucleon_density_unnormalized()

◆ integrant_nucleon_density_phi()

◆ calculate_saturation_density()

◆ set_beta_2()

◆ set_gamma()

◆ set_beta_3()

◆ set_beta_4()

◆ set_polar_angle()

◆ set_azimuthal_angle()

◆ set_angle_psi()

◆ get_beta2()

◆ get_beta3()

◆ get_beta4()

Member Data Documentation

◆ beta2_

◆ gamma_

◆ beta3_

◆ beta4_

◆ nuclear_orientation_

◆ random_rotation_