ParticleDecay.h

Go to the documentation of this file.

/*
 * Thermal-FIST package
 * 
 * Copyright (c) 2019-2019 Volodymyr Vovchenko
 *
 * GNU General Public License (GPLv3 or later)
 */
#ifndef PARTICLEDECAY_H
#define PARTICLEDECAY_H

#include <string>
#include <vector>
#include <cmath>

#include "HRGBase/ThermalModelParameters.h"
#include "HRGBase/IdealGasFunctions.h"
#include "HRGBase/xMath.h"

namespace thermalfist {
  struct Feeddown {
    enum Type { 
      Primordial = 0,      
      StabilityFlag = 1,   
      Weak = 2,            
      Electromagnetic = 3, 
      Strong = 4           
      };
    static const int NumberOfTypes = 5;
  };

  struct ParticleDecayType
  {
    enum DecayType { 
      Stable = 0,           
      Default = 1,          
      Weak = 2,             
      Electromagnetic = 3,  
      Strong = 4            
      };
    static const int NumberOfDecayTypes = 5;
  };

  struct ParticleDecayChannel {
    double mBratio;                      
    std::vector<long long> mDaughters;   
    double mM0;                          
    double mPole;
    double mL;                      
    std::vector<double> mBratioVsM; 
    double mBratioAverage;          
    std::string mChannelName;

    ParticleDecayChannel(double bratio = 0., const std::vector<long long> &daughters = std::vector<long long>(0)) :
      mBratio(bratio), mDaughters(daughters), mM0(0.), mPole(0.), mL(0.),
      mBratioVsM(std::vector<double>(0)), mBratioAverage(bratio), mChannelName("decay1") {
    }

    double ModifiedWidth(double m) const;

    bool operator==(const ParticleDecayChannel &rhs) const; // TODO: improve
    bool operator!=(const ParticleDecayChannel &rhs) const { return !(*this == rhs); }
  };  

} // namespace thermalfist

#endif