particledata.h

Go to the documentation of this file.

/*
 *    Copyright (c) 2012-2020
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 */
#ifndef SRC_INCLUDE_SMASH_PARTICLEDATA_H_
#define SRC_INCLUDE_SMASH_PARTICLEDATA_H_
 
#include <limits>
 
#include "forwarddeclarations.h"
#include "fourvector.h"
#include "particletype.h"
#include "pdgcode.h"
#include "processbranch.h"
 
namespace smash {
 
struct HistoryData {
  int32_t collisions_per_particle = 0;
  int32_t id_process = 0;
  ProcessType process_type = ProcessType::None;
  double time_last_collision = 0.0;
  PdgCode p1 = 0x0;
  PdgCode p2 = 0x0;
};
 
class ParticleData {
 public:
  explicit ParticleData(const ParticleType &particle_type, int unique_id = -1)
      : id_(unique_id), type_(&particle_type) {}
 
  int32_t id() const { return id_; }
  void set_id(int i) { id_ = i; }
 
  PdgCode pdgcode() const { return type_->pdgcode(); }
 
  // Convenience accessors to PdgCode:
  bool is_hadron() const { return type_->is_hadron(); }
 
  bool is_baryon() const { return pdgcode().is_baryon(); }
 
  bool is_nucleus() const { return pdgcode().is_nucleus(); }
 
  bool is_rho() const { return type_->is_rho(); }
 
  bool is_proton() const { return pdgcode().is_proton(); }
 
  bool is_neutron() const { return pdgcode().is_neutron(); }
 
  bool is_pion() const { return pdgcode().is_pion(); }
 
  double pole_mass() const { return type_->mass(); }
  double effective_mass() const;
  const ParticleType &type() const { return *type_; }
 
  uint32_t id_process() const { return history_.id_process; }
  HistoryData get_history() const { return history_; }
  void set_history(int ncoll, uint32_t pid, ProcessType pt, double time_of_or,
                   const ParticleList &plist);
 
  const FourVector &momentum() const { return momentum_; }
 
  void set_4momentum(const FourVector &momentum_vector) {
    momentum_ = momentum_vector;
  }
 
  void set_4momentum(double mass, const ThreeVector &mom) {
    momentum_ = FourVector(std::sqrt(mass * mass + mom * mom), mom);
  }
 
  void set_4momentum(double mass, double px, double py, double pz) {
    momentum_ = FourVector(std::sqrt(mass * mass + px * px + py * py + pz * pz),
                           px, py, pz);
  }
  void set_3momentum(const ThreeVector &mom) {
    momentum_ = FourVector(momentum_.x0(), mom);
  }
 
  const FourVector &position() const { return position_; }
  void set_4position(const FourVector &pos) { position_ = pos; }
  void set_3position(const ThreeVector &pos) {
    position_ = FourVector(position_.x0(), pos);
  }
 
  ParticleData translated(const ThreeVector &delta) const {
    ParticleData p = *this;
    p.position_[1] += delta[0];
    p.position_[2] += delta[1];
    p.position_[3] += delta[2];
    return p;
  }
 
  double formation_time() const { return formation_time_; }
  double begin_formation_time() const { return begin_formation_time_; }
 
  void set_formation_time(const double &form_time) {
    formation_time_ = form_time;
    // cross section scaling factor will be a step function in time
    begin_formation_time_ = form_time;
  }
  void set_slow_formation_times(double begin_form_time, double form_time) {
    begin_formation_time_ = begin_form_time;
    formation_time_ = form_time;
  }
 
  const double &initial_xsec_scaling_factor() const {
    return initial_xsec_scaling_factor_;
  }
  void set_cross_section_scaling_factor(const double &xsec_scal) {
    initial_xsec_scaling_factor_ = xsec_scal;
  }
 
  ThreeVector velocity() const { return momentum_.velocity(); }
 
  double inverse_gamma() const {
    return std::sqrt(1. - momentum_.sqr3() / (momentum_.x0() * momentum_.x0()));
  }
 
  void boost(const ThreeVector &v) {
    set_4momentum(momentum_.lorentz_boost(v));
    set_4position(position_.lorentz_boost(v));
  }
 
  void boost_momentum(const ThreeVector &v) {
    set_4momentum(momentum_.lorentz_boost(v));
  }
 
  bool operator==(const ParticleData &a) const { return this->id_ == a.id_; }
  bool operator<(const ParticleData &a) const { return this->id_ < a.id_; }
 
  bool operator==(int id_a) const { return this->id_ == id_a; }
  bool operator<(int id_a) const { return this->id_ < id_a; }
 
  ParticleData(const ParticleType &ptype, int uid, int index)
      : id_(uid), index_(index), type_(&ptype) {}
 
  double xsec_scaling_factor(double delta_time = 0.) const;
 
  static double formation_power_;
 
 private:
  friend class Particles;
  ParticleData() = default;
 
  void copy_to(ParticleData &dst) const {
    dst.history_ = history_;
    dst.momentum_ = momentum_;
    dst.position_ = position_;
    dst.formation_time_ = formation_time_;
    dst.initial_xsec_scaling_factor_ = initial_xsec_scaling_factor_;
    dst.begin_formation_time_ = begin_formation_time_;
  }
 
  int32_t id_ = -1;
 
  unsigned index_ = std::numeric_limits<unsigned>::max();
 
  ParticleTypePtr type_;
 
  // this leaves us two Bytes padding to use for "free"
  static_assert(sizeof(ParticleTypePtr) == 2, "");
  // make sure we don't exceed that space
  static_assert(sizeof(bool) <= 2, "");
  bool hole_ = false;
 
  FourVector momentum_;
  FourVector position_;
  double formation_time_ = 0.0;
  double begin_formation_time_ = 0.0;
  double initial_xsec_scaling_factor_ = 1.0;
  HistoryData history_;
};
 
std::ostream &operator<<(std::ostream &s, const ParticleData &p);
 
std::ostream &operator<<(std::ostream &out, const ParticleList &particle_list);
 
struct PrintParticleListDetailed {
  const ParticleList &list;
};
inline PrintParticleListDetailed detailed(const ParticleList &list) {
  return {list};
}
 
std::ostream &operator<<(std::ostream &out,
                         const PrintParticleListDetailed &particle_list);
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_SMASH_PARTICLEDATA_H_