SimpleParticle.h

Go to the documentation of this file.

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

#include <cmath>

namespace thermalfist {
  struct SimpleParticle {
    double px, py, pz; 
    double m;          
    double p0;         
    long long PDGID;         
    long long MotherPDGID;   
    int epoch;               
    bool processed;          

    SimpleParticle(double inPx = 0., double inPy = 0., double inPz = 0., double inM = 0., long long inPDGID = 0., long long inMotherPDGID = 0) :
      px(inPx),
      py(inPy),
      pz(inPz),
      m(inM),
      p0(sqrt(m*m + px * px + py * py + pz * pz)),
      PDGID(inPDGID), MotherPDGID(inMotherPDGID), epoch(0), processed(false) { }

    double GetP() const {
      return sqrt(p0*p0 - m * m);
    }

    double GetPt() const {
      return sqrt(px*px + py * py);
    }

    double GetMt() const {
      return sqrt(m*m + px * px + py * py);
    }

    double GetY() const {
      return 0.5*log((p0 + pz) / (p0 - pz));
    }

    double GetEta() const {
      //return atanh(pz/GetP());
      return 0.5*log((GetP() + pz) / (GetP() - pz));
    }

    void RapidityBoost(double dY) {
      pz = GetMt() * sinh(GetY() + dY);
      p0 = sqrt(m*m + px * px + py * py + pz * pz);
    }
  };

} // namespace thermalfist

#endif