decaytype.h

Go to the documentation of this file.

/*
 *    Copyright (c) 2015-2020
 *      SMASH Team
 *
 *    GNU General Public License (GPLv3 or later)
 */
 
#ifndef SRC_INCLUDE_DECAYTYPE_H_
#define SRC_INCLUDE_DECAYTYPE_H_
 
#include <memory>
#include <vector>
 
#include "forwarddeclarations.h"
#include "particletype.h"
#include "tabulation.h"
 
namespace smash {
 
class DecayType {
 public:
  DecayType(ParticleTypePtrList part_types, int l)
      : particle_types_(part_types), L_(l) {}
  virtual ~DecayType() = default;
  virtual unsigned int particle_number() const = 0;
  virtual bool has_particles(ParticleTypePtrList list) const = 0;
  virtual bool has_mother(ParticleTypePtr mother) const {
    SMASH_UNUSED(mother);
    return true;
  }
  const ParticleTypePtrList &particle_types() const { return particle_types_; }
  inline int angular_momentum() const { return L_; }
  virtual double width(double m0, double G0, double m) const = 0;
  virtual double in_width(double m0, double G0, double m, double m1,
                          double m2) const = 0;
 
 protected:
  ParticleTypePtrList particle_types_;
  int L_;
};
 
class TwoBodyDecay : public DecayType {
 public:
  TwoBodyDecay(ParticleTypePtrList part_types, int l);
  unsigned int particle_number() const override;
  bool has_particles(ParticleTypePtrList list) const override;
  double threshold() const {
    return particle_types_[0]->min_mass_spectral() +
           particle_types_[1]->min_mass_spectral();
  }
 
 protected:
  virtual double rho(double mass) const {
    SMASH_UNUSED(mass);
    return 1.;
  }
};
 
class TwoBodyDecayStable : public TwoBodyDecay {
 public:
  TwoBodyDecayStable(ParticleTypePtrList part_types, int l);
 
  double width(double m0, double G0, double m) const override;
 
  double in_width(double m0, double G0, double m, double m1,
                  double m2) const override;
 
 protected:
  double rho(double m) const override;
};
 
class TwoBodyDecaySemistable : public TwoBodyDecay {
 public:
  TwoBodyDecaySemistable(ParticleTypePtrList part_types, int l);
 
  double width(double m0, double G0, double m) const override;
 
  double in_width(double m0, double G0, double m, double m1,
                  double m2) const override;
 
 protected:
  double rho(double m) const override;
 
  double get_Lambda();
 
  double Lambda_;
 
  mutable std::unique_ptr<Tabulation> tabulation_;
};
 
class TwoBodyDecayUnstable : public TwoBodyDecay {
 public:
  TwoBodyDecayUnstable(ParticleTypePtrList part_types, int l);
  double width(double m0, double G0, double m) const override;
  double in_width(double m0, double G0, double m, double m1,
                  double m2) const override;
 
 protected:
  double rho(double m) const override;
  double get_Lambda();
 
  double Lambda_;
 
  mutable std::unique_ptr<Tabulation> tabulation_;
};
 
class TwoBodyDecayDilepton : public TwoBodyDecayStable {
 public:
  TwoBodyDecayDilepton(ParticleTypePtrList part_types, int l);
 
  double width(double m0, double G0, double m) const override;
};
 
class ThreeBodyDecay : public DecayType {
 public:
  ThreeBodyDecay(ParticleTypePtrList part_types, int l);
 
  unsigned int particle_number() const override;
  bool has_particles(ParticleTypePtrList list) const override;
  double width(double m0, double G0, double m) const override;
  double in_width(double m0, double G0, double m, double m1,
                  double m2) const override;
};
 
class ThreeBodyDecayDilepton : public ThreeBodyDecay {
 public:
  ThreeBodyDecayDilepton(ParticleTypePtr mother, ParticleTypePtrList part_types,
                         int l);
 
  bool has_mother(ParticleTypePtr mother) const override;
 
  static double diff_width(double m_par, double m_l, double m_dil,
                           double m_other, ParticleTypePtr other,
                           ParticleTypePtr t);
  double width(double m0, double G0, double m) const override;
 
 protected:
  mutable std::unique_ptr<Tabulation> tabulation_;
 
  ParticleTypePtr mother_;
};
 
}  // namespace smash
 
#endif  // SRC_INCLUDE_DECAYTYPE_H_