EvtRaritaSchwinger.hh

Go to the documentation of this file.

//--------------------------------------------------------------------------
//
// Environment:
//      This software is part of the EvtGen package developed jointly
//      for the BaBar and CLEO collaborations.  If you use all or part
//      of it, please give an appropriate acknowledgement.
//
// Copyright Information: See EvtGen/COPYRIGHT
//      Copyright (C) 2000      Caltech, UCSB
//
// Module: EvtGen/EvtRaritaSchwing.hh
//
// Description:Class to handle spin 3/2
//
// Modification history:
//
//    RYD     April 23, 2000         Module created
//
//------------------------------------------------------------------------
 
#ifndef EVTRARITASCHWINGER_HH
#define EVTRARITASCHWINGER_HH
 
#include "EvtComplex.hh"
#include "EvtDiracSpinor.hh"
#include "EvtVector4C.hh"
#include "EvtVector4R.hh"
 
class EvtRaritaSchwinger;
EvtRaritaSchwinger rotateEuler( const EvtRaritaSchwinger& rs, double alpha, double beta,
                                double gamma );
EvtRaritaSchwinger boostTo( const EvtRaritaSchwinger& rs, const EvtVector4R p4 );
EvtRaritaSchwinger boostTo( const EvtRaritaSchwinger& rs, const EvtVector3R boost );
EvtRaritaSchwinger dirProd( EvtVector4R v, EvtDiracSpinor u );
EvtRaritaSchwinger dirProd( EvtVector4C v, EvtDiracSpinor u );
EvtRaritaSchwinger conj( const EvtRaritaSchwinger& u2 ); // pangcy add
EvtRaritaSchwinger operator+( const EvtRaritaSchwinger& u1, const EvtRaritaSchwinger& u2 );
EvtRaritaSchwinger operator-( const EvtRaritaSchwinger& u1, const EvtRaritaSchwinger& u2 );
EvtComplex         operator*( const EvtRaritaSchwinger& u1, const EvtRaritaSchwinger& u2 );
 
//========================
class EvtRaritaSchwinger {
 
  friend EvtRaritaSchwinger rotateEuler( const EvtRaritaSchwinger& rs, double alpha,
                                         double beta, double gamma );
  friend EvtRaritaSchwinger boostTo( const EvtRaritaSchwinger& rs, const EvtVector4R p4 );
  friend EvtRaritaSchwinger boostTo( const EvtRaritaSchwinger& rs, const EvtVector3R boost );
 
  friend EvtRaritaSchwinger dirProd( EvtVector4R v, EvtDiracSpinor u );
  friend EvtRaritaSchwinger dirProd( EvtVector4C v, EvtDiracSpinor u );
  friend EvtRaritaSchwinger conj( const EvtRaritaSchwinger& u2 ); // pangcy add
  friend EvtRaritaSchwinger operator+( const EvtRaritaSchwinger& u1,
                                       const EvtRaritaSchwinger& u2 );
  friend EvtRaritaSchwinger operator-( const EvtRaritaSchwinger& u1,
                                       const EvtRaritaSchwinger& u2 );
 
  friend EvtComplex operator*( const EvtRaritaSchwinger& u1, const EvtRaritaSchwinger& u2 );
 
public:
  inline EvtRaritaSchwinger();
  virtual ~EvtRaritaSchwinger();
  inline EvtRaritaSchwinger( const EvtRaritaSchwinger& rs );
  inline EvtRaritaSchwinger& operator=( const EvtRaritaSchwinger& rs );
 
  void set( int i, int j, const EvtComplex& sp );
  //  inline const EvtComplex& get(int i, int j) const;  //pangcy add
 
  void               applyRotateEuler( double alpha, double beta, double gamma );
  void               applyBoostTo( const EvtVector4R p4 );
  void               applyBoostTo( const EvtVector3R boost );
  EvtRaritaSchwinger conj() const; // pangcy add
 
  EvtRaritaSchwinger& operator+=( const EvtRaritaSchwinger& u2 );
  EvtRaritaSchwinger& operator-=( const EvtRaritaSchwinger& u2 );
 
  EvtComplex           get( int i, int j ) const;
  friend std::ostream& operator<<( std::ostream& s, const EvtRaritaSchwinger& rs );
 
  EvtVector4C    getVector( int i ) const;
  EvtDiracSpinor getSpinor( int i ) const;
 
  void setVector( int i, const EvtVector4C& v );
  void setSpinor( int i, const EvtDiracSpinor& sp );
 
private:
  // First index in spinor index, second is Lorentz index.
  EvtComplex _rs[4][4];
};
 
EvtRaritaSchwinger::EvtRaritaSchwinger() {
 
  int i, j;
  for ( i = 0; i < 4; i++ )
  {
    for ( j = 0; j < 4; j++ ) { _rs[i][j] = 0.0; }
  }
}
 
EvtRaritaSchwinger::EvtRaritaSchwinger( const EvtRaritaSchwinger& rs ) {
 
  int i, j;
  for ( i = 0; i < 4; i++ )
  {
    for ( j = 0; j < 4; j++ ) { _rs[i][j] = rs._rs[i][j]; }
  }
}
 
EvtRaritaSchwinger& EvtRaritaSchwinger::operator=( const EvtRaritaSchwinger& rs ) {
 
  int i, j;
  for ( i = 0; i < 4; i++ )
  {
    for ( j = 0; j < 4; j++ ) { _rs[i][j] = rs._rs[i][j]; }
  }
 
  return *this;
}
 
#endif