VVKinematics.h

// -*- C++ -*-
#ifndef HERWIG_VVKinematics_H
#define HERWIG_VVKinematics_H
//
// This is the declaration of the VVKinematics class.
//
 
#include "ThePEG/Vectors/Lorentz5Vector.h"
 
namespace Herwig {
using namespace ThePEG;
 
class bornVVKinematics {
 
public:
 
  bornVVKinematics();
 
  bornVVKinematics(vector<Lorentz5Momentum> Momenta, double x1, double x2);
 
public:
 
  double x1b()   const { return x1b_; }
 
  double eta1b() const { return eta1b_; } 
 
  double x2b()   const { return x2b_; }
 
  double eta2b() const { return eta2b_; }
 
  Lorentz5Momentum p1b() const { return p1b_; }
  
  Lorentz5Momentum p2b() const { return p2b_; }
  
  Lorentz5Momentum k1b() const { return k1b_; }
  
  Lorentz5Momentum k2b() const { return k2b_; }
 
  Energy2 sb() const { return sb_; }
 
  Energy2 tb() const { return tb_; }
 
  Energy2 ub() const { return ub_; }
 
  double Yb() const { return Yb_; }
 
  Energy2 k12b() const { return k12b_; }
 
  Energy2 k22b() const { return k22b_; }
 
  double theta1b() const { return theta1b_; }
 
  void sanityCheck() const;
 
 private:
 
  double x1b_;
 
  double eta1b_;
 
  double x2b_;
 
  double eta2b_;
 
  Lorentz5Momentum p1b_;
  
  Lorentz5Momentum p2b_;
  
  Lorentz5Momentum k1b_;
  
  Lorentz5Momentum k2b_;
 
  Energy2 sb_;
 
  Energy2 tb_;
 
  Energy2 ub_;
 
  double Yb_;
 
  Energy2 k12b_;
 
  Energy2 k22b_;
 
  double theta1b_;     
};
 
 
 
class realVVKinematics {
 
public:
 
  realVVKinematics();
 
  realVVKinematics(bornVVKinematics bornVariables,double xt, double y, double theta2);
 
  void sanityCheck() const;
     
public:
 
  bornVVKinematics bornVariables() const { return bornVariables_; }
 
  double xbar() const { return xbar_; }
 
  double xt() const { return xt_; }
 
  double y() const { return y_; }
 
  double xr() const { return xr_; }
 
  double x1r() const { return x1r_; }
 
  double x2r() const { return x2r_; }
 
  Energy2 s2r() const { return s2r_; }
 
  Energy2 k12r() const { return k12r_; }
 
  Energy2 k22r() const { return k22r_; }
 
  double  theta1r() const { return theta1r_; }
 
  double  theta2r() const { return theta2r_; }
 
  Energy2 pT2_in_lab() const { return tkr_*ukr_/sr_; }
 
  Energy2 sr() const { return sr_; }
 
  Energy2 tkr() const { return tkr_; }   
 
  Energy2 ukr() const { return ukr_; }
 
  double  cpsir() const { return cpsir_; }
 
  double  cpsipr() const { return cpsiprr_; }
 
  double  betaxr() const { return betaxr_; }
 
  double  v1r() const { return v1r_; }
 
  double  v2r() const { return v2r_; }
 
  Energy2 q1r() const { return q1r_; }
 
  Energy2 q2r() const { return q2r_; } 
 
  Energy2 q1hatr() const { return q1hatr_; }
 
  Energy2 q2hatr() const { return q2hatr_; }
 
  Energy2 w1r() const { return w1r_; }
 
  Energy2 w2r() const { return w2r_; }
 
  Lorentz5Momentum p1r() const { return p1r_; }
 
  Lorentz5Momentum p2r() const { return p2r_; }
 
  Lorentz5Momentum kr()  const { return kr_ ; }
 
  Lorentz5Momentum k1r() const { return k1r_; }
 
  Lorentz5Momentum k2r() const { return k2r_; }
 
  void p1r(Lorentz5Momentum p1r) { p1r_ = p1r; }
 
  void p2r(Lorentz5Momentum p2r) { p2r_ = p2r; }
 
  void kr (Lorentz5Momentum kr ) { kr_  = kr ; }
 
  void k1r(Lorentz5Momentum k1r) { k1r_ = k1r; }
 
  void k2r(Lorentz5Momentum k2r) { k2r_ = k2r; }
 
private:
 
  bornVVKinematics bornVariables_;
 
  double xbar_;
 
  // The `raw' radiative variables.
  double xt_;
 
  double y_;
 
  double xr_;
 
  double x1r_;
 
  double x2r_;
 
  Energy2 s2r_;
 
  Energy2 k12r_;
 
  Energy2 k22r_;
 
  double  theta1r_;
 
  double  theta2r_;
 
  Energy2 sr_;
 
  Energy2 tkr_;
 
  Energy2 ukr_;
 
  double  cpsir_;
 
  double  cpsiprr_;
 
  double  betaxr_;
 
  double  v1r_;
 
  double  v2r_;
 
  Energy2 q1r_;
 
  Energy2 q2r_;
 
  Energy2 q1hatr_;
 
  Energy2 q2hatr_;
 
  Energy2 w1r_;
 
  Energy2 w2r_;
 
  Lorentz5Momentum p1r_;
 
  Lorentz5Momentum p2r_;
 
  Lorentz5Momentum kr_;
 
  Lorentz5Momentum k1r_;
 
  Lorentz5Momentum k2r_;
 
};
 
}
 
#endif /* HERWIG_VVKinematics_H */