QuarksToHadronsDecayer.h

// -*- C++ -*-
//
// QuarksToHadronsDecayer.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2007 Leif Lonnblad
//
// ThePEG is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef THEPEG_QuarksToHadronsDecayer_H
#define THEPEG_QuarksToHadronsDecayer_H
// This is the declaration of the QuarksToHadronsDecayer class.

#include "ThePEG/PDT/Decayer.h"
#include "ThePEG/Handlers/FlavourGenerator.h"

namespace ThePEG {

ThePEG_DECLARE_CLASS_POINTERS(FlavourGenerator, FlavGenPtr);

class QuarksToHadronsDecayer: public Decayer {

public:

  QuarksToHadronsDecayer()
    : theFixedN(0), theMinN(2), theC1(4.5), theC2(0.7*GeV), theC3(0.0) {}

  virtual ~QuarksToHadronsDecayer();

public:

  virtual bool accept(const DecayMode & dm) const;

  virtual ParticleVector decay(const DecayMode & dm, const Particle & p) const;

  virtual int getN(Energy m0, Energy summq, int Nq) const;

  virtual PVector getHadrons(int Nh, tcPDVector quarks) const;

  virtual void distribute(const Particle & parent, PVector & children) const;

  virtual double reweight(const Particle & parent,
                          const PVector & children) const;

public:

  int fixedN() const { return theFixedN; }

  int minN() const { return theMinN; }

  double c1() const { return theC1; }

  Energy c2() const { return theC2; }

  double c3() const { return theC3; }

  tcFlavGenPtr flavourGenerator() const { return theFlavourGenerator; }

public:


  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

  static void Init();

protected:


protected:

  virtual IBPtr clone() const;

  virtual IBPtr fullclone() const;

private:

  int theFixedN;

  int theMinN;

  double theC1;
  Energy theC2;

  double theC3;

  FlavGenPtr theFlavourGenerator;

private:

  static ClassDescription<QuarksToHadronsDecayer> initQuarksToHadronsDecayer;

  QuarksToHadronsDecayer & operator=(const QuarksToHadronsDecayer &);

};

}


namespace ThePEG {

template <>
struct BaseClassTrait<QuarksToHadronsDecayer,1>: public ClassTraitsType {
  typedef Decayer NthBase;
};

template <>
struct ClassTraits<QuarksToHadronsDecayer>
  : public ClassTraitsBase<QuarksToHadronsDecayer> {
  static string className() { return "ThePEG::QuarksToHadronsDecayer"; }
  static string library() { return "QuarksToHadronsDecayer.so"; }
};

}

#endif /* THEPEG_QuarksToHadronsDecayer_H */

---