Kinematics.h

Go to the documentation of this file.

#ifndef AMPGEN_KINEMATICS_H
#define AMPGEN_KINEMATICS_H
#include <stddef.h>
#include <array>
#include <tuple>
#include <vector>
 
#include <TLorentzVector.h>
#include <TVector3.h>
 
namespace AmpGen {
  class Event;


  class HelicityCosine {
  public:
    HelicityCosine(const std::vector<unsigned> &p1, const std::vector<unsigned> &p2, const std::vector<unsigned> &pR);
 
    HelicityCosine(const unsigned &i, const unsigned &j, const std::vector<unsigned> &pR);
 
    double operator()(std::vector<Event>::iterator evt) const;
    double operator()(const Event &evt) const;
 
  private:
    std::vector<unsigned> _i, _j, _pR;
  };

  class MomentumTransfer {
  public:
    MomentumTransfer(const std::vector<unsigned> &_p1, const std::vector<unsigned> &_p2);
    double operator()(const Event &evt) const;
 
  private:
    double Q2(const double &s, const double &s1, const double &s2) const;
    std::vector<unsigned> p1;
    std::vector<unsigned> p2;
    std::vector<unsigned> s;
  };

  double acoplanarity(const Event &evt);
 
  double PHI(const Event &evt);
  double phi(const Event &evt, int i, int j, int k, int w);
 
  std::vector<double> rotate(const std::vector<double> &input, const std::vector<double> &n, const double &v);
  void boost(Event &evt, const std::tuple<double, double, double> &n, const double &v);
  void rotate(Event &evt, const std::tuple<double, double, double> &n, const double &v);
 
  void rotateBasis(Event &evt, const TVector3 &p1, const TVector3 &p2, const TVector3 &p3);

  double dotProduct(const TLorentzVector &p1, const TLorentzVector &p2, const TLorentzVector &pX);

  TLorentzVector pFromEvent(const Event &evt, const unsigned &ref);

  TLorentzVector pFromEvent(const Event &evt, const std::vector<unsigned> &ref);
 
} // namespace AmpGen
#endif