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 */