SimpleBaryonRemnantDecayer.h

// -*- C++ -*-
//
// SimpleBaryonRemnantDecayer.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_SimpleBaryonRemnantDecayer_H
#define THEPEG_SimpleBaryonRemnantDecayer_H
//
// This is the declaration of the SimpleBaryonRemnantDecayer class.
//

#include "ThePEG/PDT/RemnantDecayer.h"
#include "ThePEG/Handlers/ZGenerator.h"
#include "ThePEG/Handlers/FlavourGenerator.h"
#include "ThePEG/Utilities/VSelector.h"

namespace ThePEG {

class SimpleBaryonRemnantDecayer: public RemnantDecayer {

public:

  typedef Ptr<ZGenerator>::pointer ZGPtr;
  typedef Ptr<FlavourGenerator>::pointer FlGPtr;

public:

  SimpleBaryonRemnantDecayer()
    : theMargin(1.0*GeV), useSpecialValence(false) {}

  SimpleBaryonRemnantDecayer(const SimpleBaryonRemnantDecayer & x)
    : RemnantDecayer(x), theZGenerator(x.theZGenerator),
      theFlavourGenerator(x.theFlavourGenerator), theMargin(x.theMargin),
      useSpecialValence(x.useSpecialValence) {}

  virtual ~SimpleBaryonRemnantDecayer();

public:

  virtual bool accept(const DecayMode & dm) const;

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

  virtual bool canHandle(tcPDPtr parent, tcPDPtr extracted) const;

  virtual bool checkExtract(tcPPtr parent, tcPPtr extracted,
                            const LorentzMomentum & pnew) const;

protected:

  virtual void doinit();

  virtual bool preInitialize() const;

public:

  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

  static void Init();

public:

  struct NoDISRespect: public Exception {};

  struct DecayFailed: public Exception {};

public:

  ZGenerator & zGenerator() const { return *theZGenerator; }

  FlavourGenerator & flavourGenerator() const { return *theFlavourGenerator; }

  Energy margin() const { return theMargin; }

  bool specialValence() const { return useSpecialValence; }

protected:

  virtual IBPtr clone() const;

  virtual IBPtr fullclone() const;

public:

  struct BaryonContent {
    vector<int> flav;
    VSelector< pair<int,int> > flavsel;
    int sign;
  };

  const BaryonContent & getBaryonInfo(tcPDPtr baryon) const;    

private:

  ZGPtr theZGenerator;

  FlGPtr theFlavourGenerator;

  Energy theMargin;

  bool useSpecialValence;

  mutable map<tcPDPtr,BaryonContent> baryonmap;

private:

  static ClassDescription<SimpleBaryonRemnantDecayer>
  initSimpleBaryonRemnantDecayer;

  SimpleBaryonRemnantDecayer & operator=(const SimpleBaryonRemnantDecayer &) = delete;

};

}

#include "ThePEG/Utilities/ClassTraits.h"

namespace ThePEG {

template <>
struct BaseClassTrait<SimpleBaryonRemnantDecayer,1> {
  typedef RemnantDecayer NthBase;
};

template <>
struct ClassTraits<SimpleBaryonRemnantDecayer>
  : public ClassTraitsBase<SimpleBaryonRemnantDecayer> {
  static string className() { return "ThePEG::SimpleBaryonRemnantDecayer"; }
};

}

#endif /* THEPEG_SimpleBaryonRemnantDecayer_H */