SpinHalfLorentzRotation.h

// -*- C++ -*-
//
// SpinHalfLorentzRotation.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef THEPEG_SpinHalfLorentzRotation_H
#define THEPEG_SpinHalfLorentzRotation_H
//
// This is the declaration of the SpinHalfLorentzRotation class.
//
#include "ThePEG/Helicity/HelicityDefinitions.h"
#include "ThreeVector.h"
 
namespace ThePEG {
 
class SpinHalfLorentzRotation {
 
  friend SpinHalfLorentzRotation inverseOf ( const SpinHalfLorentzRotation & lt );
 
public:
 
 
  SpinHalfLorentzRotation();
 
  SpinHalfLorentzRotation (double bx, double by, double bz, double gamma=-1.);
 
  SpinHalfLorentzRotation (const Boost & b,double gamma=-1.);
 
  bool isIdentity() const;
 
  SpinHalfLorentzRotation inverse() const;
 
  SpinHalfLorentzRotation & invert() { return *this = inverse(); }
 
  std::ostream & print( std::ostream & os ) const;
 
 
  SpinHalfLorentzRotation & setBoost (double bx, double by, double bz,double gamma=-1.);
 
  SpinHalfLorentzRotation & setBoost (const Boost & b,double gamma=-1.);
 
  SpinHalfLorentzRotation & setBoostX (double & boost);
 
  SpinHalfLorentzRotation & setBoostY (double & boost);
 
  SpinHalfLorentzRotation & setBoostZ (double & boost);
 
  SpinHalfLorentzRotation & setRotate(double delta, const Axis & axis);
 
  SpinHalfLorentzRotation & setRotateX (double & angle);
 
  SpinHalfLorentzRotation & setRotateY (double & angle);
 
  SpinHalfLorentzRotation & setRotateZ (double & angle);
  
 
 
  Complex s1s1() const { return _mx[0][0]; }
 
  Complex s1s2() const { return _mx[0][1]; }
 
  Complex s1s3() const { return _mx[0][2]; }
 
  Complex s1s4() const { return _mx[0][3]; }
 
  Complex s2s1() const { return _mx[1][0]; }
 
  Complex s2s2() const { return _mx[1][1]; }
 
  Complex s2s3() const { return _mx[1][2]; }
 
  Complex s2s4() const { return _mx[1][3]; }
 
  Complex s3s1() const { return _mx[2][0]; }
 
  Complex s3s2() const { return _mx[2][1]; }
 
  Complex s3s3() const { return _mx[2][2]; }
 
  Complex s3s4() const { return _mx[2][3]; }
 
  Complex s4s1() const { return _mx[3][0]; }
 
  Complex s4s2() const { return _mx[3][1]; }
 
  Complex s4s3() const { return _mx[3][2]; }
 
  Complex s4s4() const { return _mx[3][3]; }
 
  Complex operator()(unsigned int i, unsigned int j) const {
    assert(i<=3 && j<=3);
    return _mx[i][j];
  }
 
 
 
  SpinHalfLorentzRotation operator * (const SpinHalfLorentzRotation & lt) const;
 
   SpinHalfLorentzRotation & operator *= (const SpinHalfLorentzRotation & );
 
   SpinHalfLorentzRotation & transform   (const SpinHalfLorentzRotation & );
 
  SpinHalfLorentzRotation & rotateX(double delta);
 
  SpinHalfLorentzRotation & rotateY(double delta);
 
  SpinHalfLorentzRotation & rotateZ(double delta);
  
  SpinHalfLorentzRotation & rotate(double delta, const Axis & axis);
 
  SpinHalfLorentzRotation & boostX(double beta);
 
  SpinHalfLorentzRotation & boostY(double beta);
 
  SpinHalfLorentzRotation & boostZ(double beta);
 
  SpinHalfLorentzRotation & boost(double bx, double by, double bz, double gamma=-1.);
 
  SpinHalfLorentzRotation & boost(const Boost & bv, double gamma=-1.);
 
protected:
 
  SpinHalfLorentzRotation(Complex s1s1,Complex s1s2,Complex s1s3,Complex s1s4,
                                 Complex s2s1,Complex s2s2,Complex s2s3,Complex s2s4,
                                 Complex s3s1,Complex s3s2,Complex s3s3,Complex s3s4,
                                 Complex s4s1,Complex s4s2,Complex s4s3,Complex s4s4);
 
private:
 
  using MatrixT = array<array<Complex,4>,4>;
 
  SpinHalfLorentzRotation(const MatrixT & m) : _mx(m) {}
 
  MatrixT _mx;
};
 
inline SpinHalfLorentzRotation inverseOf ( const SpinHalfLorentzRotation & lt ) {
  return lt.inverse();
}
 
inline std::ostream & operator<< ( std::ostream & os,
                                   const  SpinHalfLorentzRotation& lt ) {
  return lt.print(os);
}
 
}
 
#endif /* THEPEG_SpinHalfLorentzRotation_H */