decaytype.h

Go to the documentation of this file.

/*
 *    Copyright (c) 2015-2018
 *      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_