StoSFFDecayer.h

// -*- C++ -*-
#ifndef THEPEG_StoSFFDecayer_H
#define THEPEG_StoSFFDecayer_H
//
// This is the declaration of the StoSFFDecayer class.
//
 
#include "GeneralThreeBodyDecayer.h"
#include "ThePEG/Helicity/Vertex/AbstractSSSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractRFSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractVSSVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFVVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractSSTVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFTVertex.h"
#include "ThePEG/Helicity/Vertex/AbstractFFSSVertex.h"
 
namespace Herwig {
  using namespace ThePEG;
 
class StoSFFDecayer: public GeneralThreeBodyDecayer {
 
public:
  
  double me2(const int ichan,const Particle & part,
             const tPDVector & outgoing,
             const vector<Lorentz5Momentum> & momenta,
             MEOption meopt) const;
 
  virtual void constructSpinInfo(const Particle & part,
                                 ParticleVector outgoing) const;
  
  virtual WidthCalculatorBasePtr threeBodyMEIntegrator(const DecayMode & dm) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
protected:
 
  virtual void setupDiagrams(bool checkKinematics);
 
private:
 
  StoSFFDecayer & operator=(const StoSFFDecayer &) = delete;
 
private:
  
  vector<pair<AbstractSSSVertexPtr, AbstractFFSVertexPtr> > sca_;
 
  vector<pair<AbstractFFSVertexPtr, AbstractFFSVertexPtr> > fer_;
 
  vector<pair<AbstractRFSVertexPtr, AbstractRFSVertexPtr> > RSfer_;
 
  vector<pair<AbstractVSSVertexPtr, AbstractFFVVertexPtr> > vec_;
 
  vector<pair<AbstractSSTVertexPtr, AbstractFFTVertexPtr> > ten_;
 
  vector<AbstractFFSSVertexPtr> four_;
 
  mutable RhoDMatrix rho_;
 
  mutable ScalarWaveFunction swave_;
 
  mutable pair<vector<SpinorWaveFunction>,vector<SpinorBarWaveFunction> > outspin_[3];
};
 
}
 
#endif /* THEPEG_StoSFFDecayer_H */