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);

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

    SimpleParticle LorentzBoostMomentaAndCoordinates(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
  
    
    double ParticleDistanceSquared(const SimpleParticle& part1, const SimpleParticle& part2);

    double ComputeDCA(const SimpleParticle& part);
  }
 
} // namespace thermalfist
 
#endif