GeneralCurrentDecayer.h

// -*- C++ -*-
//
// GeneralCurrentDecayer.h is a part of Herwig - A multi-purpose Monte Carlo event generator
// Copyright (C) 2002-2019 The Herwig Collaboration
//
// Herwig is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef HERWIG_GeneralCurrentDecayer_H
#define HERWIG_GeneralCurrentDecayer_H
//
// This is the declaration of the GeneralCurrentDecayer class.
//
 
#include "Herwig/Decay/DecayIntegrator.h"
#include "Herwig/Decay/WeakCurrents/WeakCurrent.h"
#include "Herwig/Decay/PhaseSpaceMode.h"
#include "ThePEG/Helicity/Vertex/VertexBase.h"
#include "GeneralCurrentDecayer.fh"
 
namespace Herwig {
using namespace ThePEG;
using Helicity::VertexBasePtr;
 
class GeneralCurrentDecayer: public DecayIntegrator {
 
public:
 
  GeneralCurrentDecayer() : 
    maxmass_(5.*GeV), wgtmax_(0.) {}
 
  virtual int modeNumber(bool & cc, tcPDPtr parent,const tPDVector & children) const;
  
  virtual Energy partialWidth(tPDPtr inpart, tPDPtr outa,
                              vector<tPDPtr> currout) = 0;
 
  void setDecayInfo(PDPtr in, PDPtr out, const vector<tPDPtr> & outCurrent,
                    VertexBasePtr vertex, WeakCurrentPtr current,
                    Energy maxmass);
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual void doinit();
 
  virtual void doinitrun();
 
protected:
 
  int modeNumber(bool & cc, vector<long> id) const;
 
  unsigned int mode() const { return mode_; }
 
  WeakCurrentPtr weakCurrent() const { return current_; }
 
  VertexBasePtr vertex() const { return theVertex_; }
 
private:
 
  GeneralCurrentDecayer & operator=(const GeneralCurrentDecayer &) = delete;
 
private:
  
  VertexBasePtr theVertex_;
  
  PDPtr inpart_;
 
  PDPtr outpart_;
 
  vector<tPDPtr> currentOut_; 
 
  WeakCurrentPtr current_;
 
  Energy maxmass_;
 
  unsigned int mode_;
 
  int wgtloc_;
 
  double wgtmax_;
 
  vector<double> weights_;
};
 
}
 
#endif /* HERWIG_GeneralCurrentDecayer_H */