FermionSpinInfo.h

// -*- C++ -*-
//
// FermionSpinInfo.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 2003-2007 Peter Richardson, 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_FermionSpinInfo_H
#define ThePEG_FermionSpinInfo_H
// This is the declaration of the FermionSpinInfo class.

#include "SpinInfo.h"
#include "ThePEG/Helicity/LorentzSpinor.h"
#include "FermionSpinInfo.fh"
// #include "FermionSpinInfo.xh"

namespace ThePEG {
namespace Helicity {

class FermionSpinInfo: public SpinInfo {

public:

  FermionSpinInfo() : SpinInfo(PDT::Spin1Half), _productionstates(2), 
                      _currentstates(2), _decaystates(2), 
                      _decaycalc(false) {}

  FermionSpinInfo(const Lorentz5Momentum & p, bool time)
    : SpinInfo(PDT::Spin1Half, p, time),
      _productionstates(2), _currentstates(2), _decaystates(2),
      _decaycalc(false) {}

public:

  void setBasisState(unsigned int hel,
                     const LorentzSpinor<SqrtEnergy> & in) const {
    assert(hel<2);
    _productionstates[hel] = in;
    _currentstates   [hel] = in;
  }

  void setDecayState(unsigned int hel,
                     const LorentzSpinor<SqrtEnergy> & in) const {
    assert(hel<2);
    _decaycalc = true;
    _decaystates[hel] = in;
  }

  const LorentzSpinor<SqrtEnergy> & getProductionBasisState(unsigned int hel) const {
    assert(hel<2);
    return _productionstates[hel];
  }

  const LorentzSpinor<SqrtEnergy> & getCurrentBasisState(unsigned int hel) const {
    assert(hel<2);
    return _currentstates[hel];
  }

  const LorentzSpinor<SqrtEnergy> & getDecayBasisState(unsigned int hel) const {
    assert(hel<2);
    if(!_decaycalc) {
      for(unsigned int ix=0;ix<2;++ix) _decaystates[ix]=_currentstates[ix];
      _decaycalc=true;
    }
    return _decaystates[hel];
  }

  virtual void transform(const LorentzMomentum &,const LorentzRotation &);

public:

  static void Init();

  virtual EIPtr clone() const;

private:

  static NoPIOClassDescription<FermionSpinInfo> initFermionSpinInfo;

  FermionSpinInfo & operator=(const FermionSpinInfo &);

private:

  mutable vector<LorentzSpinor<SqrtEnergy> > _productionstates;

  mutable vector<LorentzSpinor<SqrtEnergy> > _currentstates;

  mutable vector<LorentzSpinor<SqrtEnergy> > _decaystates;

  mutable bool _decaycalc;

};

}
}

namespace ThePEG {

template <>
struct BaseClassTrait<ThePEG::Helicity::FermionSpinInfo,1>
  : public ClassTraitsType {
  typedef ThePEG::Helicity::SpinInfo NthBase;
};

template <>
struct ClassTraits<ThePEG::Helicity::FermionSpinInfo>
  : public ClassTraitsBase<ThePEG::Helicity::FermionSpinInfo> {
  static string className() { return "ThePEG::Helicity::FermionSpinInfo"; }
};

}

#endif /* ThePEG_FermionSpinInfo_H */

---