SimpleFlavour.h

// -*- C++ -*-
//
// SimpleFlavour.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef THEPEG_SimpleFlavour_H
#define THEPEG_SimpleFlavour_H
// This is the declaration of the SimpleFlavour class.
 
#include "ThePEG/Handlers/FlavourGenerator.h"
#include "ThePEG/Utilities/VSelector.h"
// #include "SimpleFlavour.fh"
// #include "SimpleFlavour.xh"
 
namespace ThePEG {
 
class SimpleFlavour: public FlavourGenerator {
 
public:
 
  typedef map<long, VSelector< pair<long,long> > > ProbabilityMap;
 
public:
 
  SimpleFlavour();
 
  virtual ~SimpleFlavour();
 
public:
 
  virtual tcPDPair generateHadron(tcPDPtr quark) const;
 
  virtual tcPDPtr getHadron(long iq1, long iq2) const;
  using FlavourGenerator::getHadron;
 
  virtual tcPDPtr getBaryon(long q1, long q2, long q3) const;
  using FlavourGenerator::getBaryon;
 
  virtual long selectQuark() const;
 
  virtual long selectFlavour() const;
 
public:
 
  double sSup() const { return theSSup; }
 
  double diSup() const { return theDiSup; }
 
  double di1Sup() const { return theDi1Sup; }
 
  double diSSup() const  { return theDiSSup; }
 
  double etaSup() const  { return theEtaSup; }
 
  double etaPSup() const  { return theEtaPSup; }
 
  double baryon10Sup() const { return theBaryon10Sup; }
 
  double pSpin1() const { return thePSpin1; }
 
  double pSpinS1() const { return thePSpinS1; }
 
  double pSpinC1() const { return thePSpinC1; }
 
protected:
 
  virtual void setProbabilities(long iq) const;
 
  virtual double vectorMesonProbability(long iq1, long iq2) const;
 
  virtual double baryonDecupletProbability(long iq1, long iq2) const;
 
  virtual tcPDPtr pseudoScalarMeson(long iq, long iqbar) const;
 
  virtual tcPDPtr vectorMeson(long iq, long iqbar) const;
 
  virtual tcPDPtr baryonOctet(long iq, long idq) const;
 
  virtual tcPDPtr baryonDecuplet(long iq, long idq) const;
 
  virtual long pseudoScalarId(long iqh, long iql) const;
 
  virtual long vectorId(long iqh, long iql) const;
 
  virtual long baryonOctetId(long iqa, long iqb, long iqc,
                             long iq, bool dqs1) const;
 
  virtual long baryonDecupletId(long iqa, long iqb, long iqc) const;
 
  virtual vector< pair<long,double> >
  pseudoScalarIds(long iqh, long iql) const;
 
  virtual vector< pair<long,double> > vectorIds(long iqh, long iql) const;
 
  virtual vector< pair<long,double> >
  baryonOctetIds(long iqa, long iqb, long iqc,
                 long iq, bool dqs1) const;
 
  virtual vector< pair<long,double> >
  baryonDecupletIds(long iqa, long iqb, long iqc) const;
 
  void clear();
 
  static double weightSU6QDiQSpin(long iq, long idq, int spin);
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
protected:
 
  virtual void doinit();
 
  virtual void doinitrun();
 
private:
 
  double theSSup;
 
  double theDiSup;
 
  double theDi1Sup;
 
  double theDiSSup;
 
  double theEtaSup;
 
  double theEtaPSup;
 
  double theBaryon10Sup;
 
  double thePSpin1;
 
  double thePSpinS1;
 
  double thePSpinC1;
 
  mutable VSelector<long> theFlavourSelector;
 
  mutable ProbabilityMap theProbabilities;
 
 
private:
 
  static ClassDescription<SimpleFlavour> initSimpleFlavour;
 
  SimpleFlavour & operator=(const SimpleFlavour &) = delete;
 
};
 
}
 
 
namespace ThePEG {
 
template <>
struct BaseClassTrait<SimpleFlavour,1>: public ClassTraitsType {
  typedef FlavourGenerator NthBase;
};
 
template <>
struct ClassTraits<SimpleFlavour>
  : public ClassTraitsBase<SimpleFlavour> {
  static string className() { return "ThePEG::SimpleFlavour"; }
  static string library() { return "SimpleFlavour.so"; }
 
};
 
}
 
#endif /* THEPEG_SimpleFlavour_H */