ClusterDecayer.h

// -*- C++ -*-
//
// ClusterDecayer.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_ClusterDecayer_H
#define HERWIG_ClusterDecayer_H
 
#include <ThePEG/Interface/Interfaced.h>
#include <ThePEG/EventRecord/Step.h>
#include "CluHadConfig.h"
#include "HadronSelector.h"
#include "ClusterDecayer.fh"
 
namespace Herwig {
using namespace ThePEG;
 
class ClusterDecayer: public Interfaced {
 
public:
 
  ClusterDecayer();
 
  void decay(const ClusterVector & clusters, tPVector & finalhadrons)
   ;
 
public:
 
  void persistentOutput(PersistentOStream &) const;
 
  void persistentInput(PersistentIStream &, int);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
private:
 
  ClusterDecayer & operator=(const ClusterDecayer &) = delete;
 
public:
 
  pair<PPtr,PPtr> decayIntoTwoHadrons(tClusterPtr ptr);
 
private:
 
  void calculatePositions( const Lorentz5Momentum &, const LorentzPoint &,
                           const Lorentz5Momentum &, const Lorentz5Momentum &,
                           LorentzPoint &, LorentzPoint &) const;
 
  Ptr<HadronSelector>::pointer _hadronsSelector;
 
 
  bool _clDirLight;
  bool _clDirBottom;
  bool _clDirCharm;
  bool _clDirExotic;
 
  double _clSmrLight;
  double _clSmrBottom;
  double _clSmrCharm;
  double _clSmrExotic;
 
  bool _onshell;
 
  unsigned int _masstry;
 
 
};
 
 
}
 
#endif /* HERWIG_ClusterDecayer_H */