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