suep_shower.h

Go to the documentation of this file.

/*
 * This file is taken from the available public code at: 
 * https://gitlab.com/simonknapen/suep_generator 
 * by Simon Knapen.
 */
#ifndef SUEPSHOWER_h
#define SUEPSHOWER_h
 
#include <vector>
#include <math.h>
#include <functional> // for bind() 
#include <boost/math/tools/roots.hpp>
#include "Pythia8/PhaseSpace.h"
#include "Pythia8/Basics.h"
 
//** Auxiliary class for tolerance checks */
class tolerance {
 public:
  tolerance(double eps) :
  m_eps(eps) {
  }
  bool operator()(double a, double b) {
    return (std::abs(b - a) <= m_eps);
  }
 private:
  double m_eps;
};

class Suep_shower 
{ 
 public:
  
  Suep_shower(double mass, double temperature, double energy, Pythia8::Rndm *rndm = 0);
  
  std::vector< Pythia8::Vec4 > generate_shower();
  
 protected: 

  Pythia8::Rndm *m_rndmEngine;
  
  double m_m;
  
  double m_Temp;
  
  double m_Etot; 
  
 protected:
  
  double f(double p);
  
  double fp(double p);
  
  double test_fun(double p);
  
  Pythia8::Vec4 generateFourVector(); 
  
  double reballance_func(double a, const std::vector< Pythia8::Vec4 > &event);
  
 private:
  double m_A;
  double m_p_m;
  double m_p_plus, m_p_minus;
  double m_lambda_plus,m_lambda_minus;
  double m_q_plus, m_q_minus, m_q_m;
  
}; 
 
#endif