Hw64Decayer.h

// -*- C++ -*-
//
// Hw64Decayer.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_Hw64Decayer_H
#define HERWIG_Hw64Decayer_H
//
// This is the declaration of the Hw64Decayer class.
//
#include <ThePEG/Config/ThePEG.h>
#include <HwDecayerBase.h>
#include <ThePEG/Interface/Interfaced.h>
#include <ThePEG/PDT/DecayMode.h>
#include <ThePEG/Repository/Strategy.fh>
#include <fstream>
 
namespace Herwig {
 
using namespace ThePEG;
 
class Hw64Decayer: public HwDecayerBase {
 
public:
 
  Hw64Decayer() : MECode(0),_masstry(50) {} 
 
  virtual bool accept(tcPDPtr parent, const tPDVector & children) const;
 
  virtual ParticleVector decay(const Particle & parent,
                               const tPDVector & children) const;
 
  virtual void dataBaseOutput(ofstream & os,bool header) const;
 
public:
 
  static void Init();
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
protected:
 
  virtual IBPtr clone() const {return new_ptr(*this);}
 
  virtual IBPtr fullclone() const {return new_ptr(*this);}
 
private:
 
  static double VAWt(Energy2 t0, Energy2 t1, Energy2 t2, InvEnergy4 t3);
 
  void setParticleMomentum(ParticleVector & out, const cPDVector & particles, 
                           const vector<Lorentz5Momentum> & moms) const {
    unsigned int numProds = particles.size();
    for(unsigned int ix=0;ix<numProds;++ix)
      out.push_back(particles[ix]->produceParticle(moms[ix]));
  }
 
private:
 
  const Hw64Decayer & operator=(const Hw64Decayer &) = delete;
 
private:
 
  int MECode;
 
  unsigned int _masstry;
};
 
}
 
#endif /* HERWIG_Hw64Decayer_H */