f1FourPiDecayer.h

// -*- C++ -*-
#ifndef Herwig_f1FourPiDecayer_H
#define Herwig_f1FourPiDecayer_H
//
// This is the declaration of the f1FourPiDecayer class.
//
 
#include "DecayIntegrator.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class f1FourPiDecayer: public DecayIntegrator {
 
public:
 
  f1FourPiDecayer();
 
  virtual int modeNumber(bool & cc, tcPDPtr parent, 
                         const tPDVector & children) const;
 
  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 void dataBaseOutput(ofstream & os,bool header) const;
  
  virtual ParticleVector decay(const Particle & parent,
                               const tPDVector & children) 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 doinit();
 
  virtual void doinitrun();
 
private:
 
  f1FourPiDecayer & operator=(const f1FourPiDecayer &) = delete;
 
private:
 
  double gRhoPiPi_;
 
  Energy ga1RhoPi_;
 
  InvEnergy gf1a1Pi_;
 
  Energy ma1_;
 
  Energy ga1_;
 
  Energy mrho_;
 
  Energy grho_;
  
 
  vector<double> maxWeight_;
 
  mutable RhoDMatrix rho_;
 
  mutable vector<Helicity::LorentzPolarizationVector> vector_;
  
};
 
}
 
#endif /* Herwig_f1FourPiDecayer_H */