HwDecayerBase.h

// -*- C++ -*-
//
// HwDecayerBase.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_HwDecayerBase_H
#define HERWIG_HwDecayerBase_H
//
// This is the declaration of the HwDecayerBase class.
//
 
#include "ThePEG/PDT/Decayer.h"
#include "Herwig/Shower/RealEmissionProcess.fh"
#include "HwDecayerBase.fh"
 
namespace Herwig {
 
struct Branching;
 
using namespace ThePEG;
 
class HwDecayerBase: public Decayer {
 
public:
 
  HwDecayerBase() : _initialize(false), _dbOutput(false) {}
 
  virtual bool accept(const DecayMode & dm) const;
 
  virtual ParticleVector decay(const DecayMode & dm, const Particle & p) const;
 
public:
 
  
  enum POWHEGType {No, ISR, FSR, Both};
 
 
  virtual POWHEGType hasPOWHEGCorrection() {return No;}
 
 
  virtual bool hasMECorrection() {return false;}
 
  virtual void initializeMECorrection(RealEmissionProcessPtr , double & ,
                                      double & );
 
  virtual RealEmissionProcessPtr applyHardMatrixElementCorrection(RealEmissionProcessPtr);
 
  virtual bool softMatrixElementVeto(PPtr parent,
                                     PPtr progenitor,
                                     const bool & fs,
                                     const Energy & highestpT,
                                     const vector<tcPDPtr> & ids,
                                     const double & z,
                                     const Energy & scale,
                                     const Energy & pT);
 
  virtual RealEmissionProcessPtr generateHardest(RealEmissionProcessPtr);
 
protected:
 
  virtual bool accept(tcPDPtr parent, const tPDVector & children) const = 0;
  
  virtual ParticleVector decay(const Particle & parent,
                               const tPDVector & children) const = 0;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
public:
 
  virtual void dataBaseOutput(ofstream & os,bool header) const = 0;
 
  bool initialize() const {return _initialize;}
 
  bool databaseOutput() const {return _dbOutput;}
 
protected:
 
  void fixRho(RhoDMatrix &) const;
  
protected:
 
  virtual void dofinish();
 
private:
 
  HwDecayerBase & operator=(const HwDecayerBase &) = delete;
 
private:
 
  bool _initialize;
 
  bool _dbOutput;
};
 
}
 
#endif /* HERWIG_HwDecayerBase_H */