collidermodus.h

Go to the documentation of this file.

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

#include <cstring>
#include <memory>
#include <utility>

#include "deformednucleus.h"
#include "forwarddeclarations.h"
#include "fourvector.h"
#include "interpolation.h"
#include "modusdefault.h"
#include "nucleus.h"
#include "pdgcode.h"

namespace smash {

struct ExperimentParameters;

class ColliderModus : public ModusDefault {
 public:
  explicit ColliderModus(Configuration modus_config,
                         const ExperimentParameters &parameters);

  double initial_conditions(Particles *particles,
                            const ExperimentParameters &parameters);
  int total_N_number() const { return target_->size() + projectile_->size(); }
  int proj_N_number() const { return projectile_->size(); }
  double velocity_projectile() const { return velocity_projectile_; }
  double velocity_target() const { return velocity_target_; }
  bool cll_in_nucleus() { return cll_in_nucleus_; }
  FermiMotion fermi_motion() { return fermi_motion_; }
  bool is_collider() const { return true; }
  double impact_parameter() const { return impact_; }
  struct ColliderEmpty : public ModusDefault::BadInput {
    using ModusDefault::BadInput::BadInput;
  };

 private:
  std::unique_ptr<Nucleus> projectile_;
  std::unique_ptr<Nucleus> target_;
  double total_s_;
  double sqrt_s_NN_;
  static std::unique_ptr<DeformedNucleus> create_deformed_nucleus(
      Configuration &nucleus_cfg, const int ntest,
      const std::string &nucleus_type);
  double impact_ = 0.;
  Sampling sampling_ = Sampling::Quadratic;
  double imp_min_ = 0.0;
  double imp_max_ = 0.0;
  double yield_max_ = 0.0;
  std::unique_ptr<InterpolateDataLinear<double>> impact_interpolation_ =
      nullptr;
  void sample_impact();
  double initial_z_displacement_ = 2.0;
  CalculationFrame frame_ = CalculationFrame::CenterOfVelocity;
  FermiMotion fermi_motion_ = FermiMotion::Off;
  bool cll_in_nucleus_ = false;
  double velocity_projectile_ = 0.0;
  double velocity_target_ = 0.0;
  std::pair<double, double> get_velocities(double mandelstam_s, double m_a,
                                           double m_b);

  friend std::ostream &operator<<(std::ostream &, const ColliderModus &);
};

}  // namespace smash

#endif  // SRC_INCLUDE_COLLIDERMODUS_H_