SpinHalfLorentzRotation.h

// -*- C++ -*-
//
// SpinHalfLorentzRotation.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2007 Leif Lonnblad
//
// ThePEG is licenced under version 2 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:


  vector<vector<Complex> > _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 */

---