Rank3TensorSpinInfo.h

// -*- C++ -*-
//
// Rank3TensorSpinInfo.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 2003-2019 Peter Richardson, 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_Rank3TensorSpinInfo_H
#define THEPEG_Rank3TensorSpinInfo_H
// This is the declaration of the Rank3TensorSpinInfo class.
 
#include "ThePEG/EventRecord/SpinInfo.h"
#include "ThePEG/Helicity/LorentzRank3Tensor.h"
#include "Rank3TensorSpinInfo.fh"
// #include "Rank3TensorSpinInfo.xh"
#include <array>
 
namespace ThePEG {
namespace Helicity {
 
class Rank3TensorSpinInfo: public SpinInfo {
 
public:
 
  Rank3TensorSpinInfo() : SpinInfo(PDT::Spin3), _decaycalc(false) {}
  
  Rank3TensorSpinInfo(const Lorentz5Momentum & p, bool time)
    : SpinInfo(PDT::Spin3, p, time), _decaycalc(false) {}
 
public:
 
  void setBasisState(unsigned int hel, LorentzRank3Tensor<double> in) const {
    assert(hel<7);
    _productionstates[hel]=in;
    _currentstates   [hel]=in;
  }
 
  void setDecayState(unsigned int hel, LorentzRank3Tensor<double> in) const {
    assert(hel<7);
    _decaycalc = true;
    _decaystates[hel] = in;
  }
 
  const LorentzRank3Tensor<double> & getProductionBasisState(unsigned int hel) const {
    assert(hel<7);
    return _productionstates[hel];
  }
  
  const LorentzRank3Tensor<double> & getCurrentBasisState(unsigned int hel) const {
    assert(hel<7);
    return _currentstates[hel];
  }
 
  const LorentzRank3Tensor<double> & getDecayBasisState(unsigned int hel) const {
    assert(hel<7);
    if(!_decaycalc) {
      for(unsigned int ix=0;ix<7;++ix)
        _decaystates[ix]=_currentstates[ix].conjugate();
      _decaycalc=true;
    }
    return _decaystates[hel];
  }
 
  virtual void transform(const LorentzMomentum &,const LorentzRotation &);
 
  virtual void undecay() const {
    _decaycalc=false;
    SpinInfo::undecay();
  }
 
  virtual void reset() {
    undecay();
    _currentstates = _productionstates;
    SpinInfo::reset();
  }
 
public:
 
  static void Init();
 
  virtual EIPtr clone() const;
 
private:
 
  Rank3TensorSpinInfo & operator=(const Rank3TensorSpinInfo &) = delete;
 
private:
 
  mutable std::array<LorentzRank3Tensor<double>,7> _productionstates;
 
  mutable std::array<LorentzRank3Tensor<double>,7> _decaystates;
 
  mutable std::array<LorentzRank3Tensor<double>,7> _currentstates;
 
  mutable bool _decaycalc;
 
};
 
}
}
 
 
namespace ThePEG {
 
}
#endif /* THEPEG_Rank3TensorSpinInfo_H */