isoparticletype.h

Go to the documentation of this file.

/*
 *    Copyright (c) 2015-2018
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 */

#ifndef SRC_INCLUDE_ISOPARTICLETYPE_H_
#define SRC_INCLUDE_ISOPARTICLETYPE_H_

#include <string>
#include <unordered_map>

#include "particletype.h"
#include "tabulation.h"

namespace smash {

class IsoParticleType {
 public:
  IsoParticleType(const std::string &n, double m, double w, unsigned int s,
                  Parity p);

  IsoParticleType(const IsoParticleType &) = delete;
  IsoParticleType &operator=(const IsoParticleType &) = delete;

  IsoParticleType(IsoParticleType &&) = default;
  IsoParticleType &operator=(IsoParticleType &&) = default;

  bool operator==(const IsoParticleType &rhs) const {
    return states_[0]->pdgcode() == rhs.states_[0]->pdgcode();
  }

  const std::string &name() const { return name_; }

  double mass() const { return mass_; }

  double width() const { return width_; }

  int isospin() const { return states_.size() - 1; }

  unsigned int spin() const { return spin_; }

  Parity parity() const { return parity_; }

  bool is_hadron() const { return states_[0]->is_hadron(); }

  ParticleTypePtrList get_states() const { return states_; }

  void add_state(const ParticleType &type);

  bool has_anti_multiplet() const;

  static const IsoParticleTypeList &list_all();

  static const IsoParticleType *try_find(const std::string &name);

  static const IsoParticleType &find(const std::string &name);

  static IsoParticleType *find(const ParticleType &type);

  struct ParticleNotFoundFailure : public std::runtime_error {
    using std::runtime_error::runtime_error;
  };

  static bool exists(const std::string &name);

  static const ParticleTypePtr find_state(const std::string &name);

  static void create_multiplet(const ParticleType &type);

  double get_integral_NR(double sqrts);

  double get_integral_RR(const ParticleType &type_res_2, double sqrts);

  TabulationPtr integrate_RR(ParticleTypePtr &type_res_2);

  double get_integral_RK(double sqrts);

  double get_integral_piR(double sqrts);

 private:
  std::string name_;
  double mass_;
  double width_;
  unsigned int spin_;
  Parity parity_;
  ParticleTypePtrList states_;

  TabulationPtr XS_piR_tabulation_;
  TabulationPtr XS_RK_tabulation_;
  TabulationPtr XS_NR_tabulation_;
  TabulationPtr XS_DR_tabulation_;
  std::unordered_map<IsoParticleType *, TabulationPtr> XS_RR_tabulations;

  static IsoParticleType &find_private(const std::string &name);
};

}  // namespace smash

#endif  // SRC_INCLUDE_ISOPARTICLETYPE_H_