Hw7Selector.h

// -*- C++ -*-
//
// Hw7Selector.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_Hw7Selector_H
#define HERWIG_Hw7Selector_H
//
// This is the declaration of the Hw7Selector class.
//
 
#include "HadronSelector.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class Hw7Selector: public HadronSelector {
 
public:
 
  Hw7Selector() : HadronSelector(1),
                   _pwtDquark( 1.0 ),_pwtUquark( 1.0 ),_pwtSquark( 1.0 ),_pwtCquark( 0.0 ),
                   _pwtBquark( 0.0 ),_pwtDIquarkS0( 1.0 ),_pwtDIquarkS1( 1.0 ),
                   _sngWt( 1.0 ), _decWt( 1.0 ),
                   _mode(1), _enhanceSProb(0), _m0Decay(1.*GeV),
                   _scHadronWtFactor(1.), _sbHadronWtFactor(1.)
  {}
 
  virtual PDPtr makeDiquark(tcPDPtr par1, tcPDPtr par2);
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
  
  virtual double baryonWeight(long id) const;
 
  std::tuple<bool,bool,bool> selectBaryon(const Energy cluMass, tcPDPtr par1, tcPDPtr par2) const;
 
  virtual double strangeWeight(const Energy cluMass, tcPDPtr par1, tcPDPtr par2) const;
 
  virtual void insertOneHalf(HadronInfo a, int flav1, int flav2);
 
  virtual void insertThreeHalf(HadronInfo a, int flav1, int flav2);
  
  tcPDPair lightestBaryonPair(tcPDPtr ptr1, tcPDPtr ptr2, int pspin) const;
 
protected:
 
   virtual IBPtr clone() const;
 
   virtual IBPtr fullclone() const;
 
protected:
 
  virtual void doinit();
 
private:
 
  Hw7Selector & operator=(const Hw7Selector &) = delete;
 
private:
 
  double _pwtDquark;
 
  double _pwtUquark;
 
  double _pwtSquark;
 
  double _pwtCquark;
 
  double _pwtBquark;
 
  double _pwtDIquarkS0;
 
  double _pwtDIquarkS1;
 
  double _sngWt;
 
  double _decWt;
  
private:
 
  unsigned int _mode;
 
  int _enhanceSProb;
 
  Energy _m0Decay;
 
  int _massMeasure;
 
  const double _maxScale = 20.;
 
  double _scHadronWtFactor;
 
  double _sbHadronWtFactor;
 
  map<pair<long,long>,tcPDPair> lightestBaryonsS0_,lightestBaryonsS1_;
};
 
}
 
#endif /* HERWIG_Hw7Selector_H */