ParticleDecaysMC.h

Go to the documentation of this file.

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

#include <cmath>
#include <vector>
#include <map>

#include "HRGEventGenerator/SimpleParticle.h"

namespace thermalfist {

  namespace ParticleDecaysMC {
    double ThreeBodym12F2(double m12, double M, double m1, double m2, double m3);

    double TernaryThreeBodym12Maximum(double M, double m1, double m2, double m3);

    extern int threebodysucc, threebodytot;

    double GetRandomThreeBodym12(double M, double m1, double m2, double m3, double fm12max);

    SimpleParticle LorentzBoost(const SimpleParticle &part, double vx, double vy, double vz);

    std::vector<SimpleParticle> TwoBodyDecay(const SimpleParticle & Mother, double m1, long long pdg1, double m2, long long pdg2);
    
    std::vector<SimpleParticle> ManyBodyDecay(const SimpleParticle & Mother, std::vector<double> masses, std::vector<long long> pdgs); // TODO: proper implementation for 4+ - body decays


    void ShuffleDecayProducts(std::vector<double> &masses, std::vector<long long> &pdgs); // TODO: proper implementation for 4+ - body decays
  }

} // namespace thermalfist

#endif