PartonicDecayerBase.h

// -*- C++ -*-
//
// PartonicDecayerBase.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_PartonicDecayerBase_H
#define HERWIG_PartonicDecayerBase_H
//
// This is the declaration of the PartonicDecayerBase class.
//
 
#include "Herwig/Decay/HwDecayerBase.h"
#include "Herwig/Hadronization/PartonSplitter.h"
#include "Herwig/Hadronization/ClusterFinder.h"
#include "Herwig/Hadronization/ClusterFissioner.h"
#include "Herwig/Hadronization/LightClusterDecayer.h"
#include "Herwig/Hadronization/ClusterDecayer.h"
#include "Herwig/Hadronization/Cluster.h"
#include "Herwig/Shower/QTilde/Base/PartnerFinder.h"
#include "Herwig/Shower/QTilde/SplittingFunctions/SplittingGenerator.h"
#include "Herwig/Shower/QTilde/Kinematics/KinematicsReconstructor.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
ThePEG_DECLARE_MULTISET(tcPDPtr,cParticleMSet);
 
class PartonicDecayerBase: public HwDecayerBase {
 
public:
 
  PartonicDecayerBase();
 
  virtual bool accept(tcPDPtr parent, const tPDVector & children) const = 0;
  
  virtual ParticleVector decay(const Particle & parent,
                               const tPDVector & children) const = 0;
 
  virtual ParticleVector decay(const DecayMode & dm, const Particle & p) const;
 
  virtual void dataBaseOutput(ofstream & os,bool header) const;
 
protected:
 
  PVector shower(const Particle & parent, ParticleVector & partons) const;
  
  bool timeLikeShower(tShowerParticlePtr particle,
                      Branching fb, bool first, PVector & output) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual void doinit();
 
protected:
 
  bool duplicateMode(const Particle & parent,
                     const vector<tPPtr> & hadrons) const;
 
private:
 
  PartonicDecayerBase & operator=(const PartonicDecayerBase &) = delete;
 
private:
 
  PartonSplitterPtr      _partonSplitter;
 
  ClusterFinderPtr       _clusterFinder;
 
  ClusterFissionerPtr    _clusterFissioner;
 
  LightClusterDecayerPtr _lightClusterDecayer;
 
  ClusterDecayerPtr      _clusterDecayer; 
 
  bool _exclusive;
 
  unsigned int _partontries;
 
  bool _inter;
 
  // @name Parton Shower related variables
 
  bool _shower;
  
  SplittingGeneratorPtr _splittingGenerator;
 
  PartnerFinderPtr _partnerFinder;
  
  KinematicsReconstructorPtr _reconstructor;
 
};
 
}
 
#endif /* HERWIG_PartonicDecayerBase_H */