EtaPiPiPiDecayer.h

// -*- C++ -*-
//
// EtaPiPiPiDecayer.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_EtaPiPiPiDecayer_H
#define HERWIG_EtaPiPiPiDecayer_H
// This is the declaration of the EtaPiPiPiDecayer class.
 
#include "Herwig/Decay/DecayIntegrator.h"
#include "Herwig/Decay/PhaseSpaceMode.h"
 
namespace Herwig {
using namespace ThePEG;
 
class EtaPiPiPiDecayer: public DecayIntegrator {
 
public:
 
  EtaPiPiPiDecayer() {
    // no intermediates
    generateIntermediates(false);
  }
 
  
  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 WidthCalculatorBasePtr threeBodyMEIntegrator(const DecayMode & dm) const;
  
  virtual InvEnergy threeBodydGammads(const int imode, const Energy2 q2, const  Energy2 s,
                                   const Energy m1, const Energy m2, 
                                   const Energy m3) const;
 
  virtual void dataBaseOutput(ofstream & os,bool header) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
protected:
 
  virtual void doinit();
 
  virtual void doinitrun();
 
private:
 
  EtaPiPiPiDecayer & operator=(const EtaPiPiPiDecayer &) = delete;
 
public:
 
  string setUpDecayMode(string arg);
 
private:
 
  vector<int> incoming_;
 
  vector<int> outgoing_;
 
  vector<bool> charged_;
 
  vector<double> prefactor_;
 
  //*{
  vector<double> a_;
 
  vector<double> b_;
 
  vector<double> c_;
 
  vector<double> maxWeight_;
 
  mutable RhoDMatrix rho_;
};
 
}
 
 
#endif /* HERWIG_EtaPiPiPiDecayer_H */