HwppSelector.h

// -*- C++ -*-
//
// HwppSelector.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_HwppSelector_H
#define HERWIG_HwppSelector_H
//
// This is the declaration of the HwppSelector class.
//
 
#include "HadronSelector.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class HwppSelector: public HadronSelector {
 
public:
 
  HwppSelector() : HadronSelector(1),
                   _pwtDquark( 1.0 ),_pwtUquark( 1.0 ),_pwtSquark( 1.0 ),_pwtCquark( 0.0 ),
                   _pwtBquark( 0.0 ),_pwtDIquark(1.0 ),
                   _sngWt( 1.0 ), _decWt( 1.0 ),
                   _mode(1), _enhanceSProb(0), _m0Decay(1.*GeV),
                   _scHadronWtFactor(1.), _sbHadronWtFactor(1.)
  {}
 
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;
 
  inline Energy massLightestBaryonPair(tcPDPtr ptr1, tcPDPtr ptr2) const {
    map<pair<long,long>,tcPDPair>::const_iterator lightest =
      lightestBaryons_.find(make_pair(abs(ptr1->id()),abs(ptr2->id())));
    assert(lightest!=lightestBaryons_.end());
    return lightest->second.first->mass()+lightest->second.second->mass();
  }
  
  tcPDPair lightestBaryonPair(tcPDPtr ptr1, tcPDPtr ptr2) const;
 
protected:
 
   virtual IBPtr clone() const;
 
   virtual IBPtr fullclone() const;
 
protected:
 
  virtual void doinit();
 
private:
 
  HwppSelector & operator=(const HwppSelector &) = delete;
 
private:
 
  double _pwtDquark;
 
  double _pwtUquark;
 
  double _pwtSquark;
 
  double _pwtCquark;
 
  double _pwtBquark;
 
  double _pwtDIquark;
 
  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> lightestBaryons_;
 
};
 
}
 
#endif /* HERWIG_HwppSelector_H */