ClusterHadronizationHandler.h

// -*- C++ -*-
//
// ClusterHadronizationHandler.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_ClusterHadronizationHandler_H
#define HERWIG_ClusterHadronizationHandler_H
 
#include <ThePEG/Handlers/HadronizationHandler.h>
#include "PartonSplitter.h"
#include "ClusterFinder.h"
#include "ColourReconnector.h"
#include "ClusterFissioner.h"
#include "LightClusterDecayer.h"
#include "ClusterDecayer.h"
#include "ClusterHadronizationHandler.fh"
#include "Herwig/Utilities/Reshuffler.h"
#include "GluonMassGenerator.h"
 
namespace Herwig {
using namespace ThePEG;
 
 
class ClusterHadronizationHandler: 
    public HadronizationHandler, public Reshuffler {
 
public:
 
  virtual void handle(EventHandler & ch, const tPVector & tagged,
                      const Hint & hint);
 
  Energy2 minVirtuality2() const
  { return _minVirtuality2; }
 
  Length maxDisplacement() const
  { return _maxDisplacement; }
 
  bool isSoftUnderlyingEventON() const
  { return _underlyingEventHandler; }
  
  static const ClusterHadronizationHandler * currentHandler() {
    if(!currentHandler_){
      cerr<< " \nCreating new ClusterHadronizationHandler without input from infiles.";
      cerr<< " \nWhen using for example the string model ";
      cerr<< " hadronic decays are still treated by the Cluster model\n";
      currentHandler_=new ClusterHadronizationHandler();;
    }
    return currentHandler_;
  }
 
  Ptr<GluonMassGenerator>::tptr gluonMassGenerator() const {
    return gluonMassGenerator_;
  }
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
private:
 
  ClusterHadronizationHandler & operator=(const ClusterHadronizationHandler &) = delete;
 
  PartonSplitterPtr      _partonSplitter;
 
  ClusterFinderPtr       _clusterFinder;
 
  ColourReconnectorPtr   _colourReconnector;
 
  ClusterFissionerPtr    _clusterFissioner;
 
  LightClusterDecayerPtr _lightClusterDecayer;
 
  ClusterDecayerPtr      _clusterDecayer; 
 
  bool reshuffle_ = false;
  
  int reshuffleMode_ = 0;
  
  Ptr<GluonMassGenerator>::ptr gluonMassGenerator_;
 
  Energy2 _minVirtuality2 = 0.1_GeV2;
 
  Length _maxDisplacement = 1.0e-10_mm;
 
  StepHdlPtr _underlyingEventHandler;
 
  bool _reduceToTwoComponents = true;
 
  void _setChildren(const ClusterVector & clusters) const;
 
   
   void splitIntoColourSinglets(PVector thelist,
                                vector<PVector>& reshufflelists);
 
  
  static ClusterHadronizationHandler * currentHandler_;
 
 
};
 
 
}
 
#endif /* HERWIG_ClusterHadronizationHandler_H */