SMHiggsGGHiggsPPDecayer.h

// -*- C++ -*-
//
// SMHiggsGGHiggsPPDecayer.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_SMHiggsGGHiggsPPDecayer_H
#define HERWIG_SMHiggsGGHiggsPPDecayer_H
//
// This is the declaration of the SMHiggsGGHiggsPPDecayer class.
//
 
#include "Herwig/Decay/PerturbativeDecayer.h"
#include "Herwig/Decay/PhaseSpaceMode.h"
#include "Herwig/Models/StandardModel/StandardModel.h"
#include "ThePEG/Helicity/Vertex/AbstractVVSVertex.h"
 
namespace Herwig {
using namespace ThePEG;
using namespace ThePEG::Helicity;
  
class SMHiggsGGHiggsPPDecayer: public PerturbativeDecayer {
  
public:
 
  SMHiggsGGHiggsPPDecayer() : _h0wgt(3,1.), _minloop(6), _maxloop(6), _massopt(0)
  {}
  
  double me2(const int ichan,const Particle & part,
             const tPDVector & outgoing,
             const vector<Lorentz5Momentum> & momenta,
             MEOption meopt) const;
 
  virtual void constructSpinInfo(const Particle & part,
                                 ParticleVector outgoing) const;
  
  virtual bool accept(tcPDPtr parent, const tPDVector & children) const;
  
  virtual int modeNumber(bool &, tcPDPtr, const tPDVector & ) const {return -1;}
  
  virtual ParticleVector decay(const Particle & parent,const tPDVector & children) const;
 
  virtual void dataBaseOutput(ofstream & os,bool header) const;
 
  virtual double matrixElementRatio(const Particle & inpart, const ParticleVector & decay2,
                                    const ParticleVector & decay3, MEOption meopt,
                                    ShowerInteraction inter);
 
  virtual POWHEGType hasPOWHEGCorrection() {return FSR;}
  
public:
  
  void persistentOutput(PersistentOStream & os) const;
  
  void persistentInput(PersistentIStream & is, int version);
  
  static void Init();
  
protected:
  
  virtual IBPtr clone() const {return new_ptr(*this);}
  
  virtual IBPtr fullclone() const {return new_ptr(*this);}
  
protected:
  
  virtual void doinit();
  
  virtual void doinitrun();
 
protected:
 
  double realME(const vector<cPDPtr> & partons, 
                const vector<Lorentz5Momentum> & momenta) const;
  
  Energy2 loME(Energy mh) const;
  
private:
  
  SMHiggsGGHiggsPPDecayer & operator=(const SMHiggsGGHiggsPPDecayer &) = delete;
  
  AbstractVVSVertexPtr _hggvertex;
  
  AbstractVVSVertexPtr _hppvertex;
  
  AbstractVVSVertexPtr _hzpvertex;
  
  vector<double> _h0wgt;
  
  mutable RhoDMatrix _rho;
 
  mutable ScalarWaveFunction _swave;
 
  mutable vector<VectorWaveFunction> _vwave[2];
 
private:
 
  int _minloop;
 
  int _maxloop;
 
  unsigned int _massopt;
};
  
}
 
#endif /* HERWIG_SMHiggsGGHiggsPPDecayer_H */