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) {}
 
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 */