SMHiggsWidthGenerator.h

// -*- C++ -*-
//
// SMHiggsWidthGenerator.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_SMHiggsWidthGenerator_H
#define HERWIG_SMHiggsWidthGenerator_H
//
// This is the declaration of the SMHiggsWidthGenerator class.
//
 
#include "GenericWidthGenerator.h"
 
namespace Herwig {
using namespace ThePEG;
 
class SMHiggsWidthGenerator: public GenericWidthGenerator {
 
public:
 
  SMHiggsWidthGenerator() 
    : widthopt_(2), offshell_(10.), mw_(ZERO), mz_(ZERO), gamw_(ZERO),
      gamz_(ZERO), qmass_(7,ZERO), lmass_(3,ZERO),
      sw2_(0.), ca_(0.), cf_(0.), qlast_(ZERO)
  {}
 
private:
 
 
  virtual Energy width(const ParticleData &, Energy m) const;
 
  virtual DecayMap rate(const ParticleData &) const;
 
  virtual DecayMap rate(const Particle &);
 
  virtual pair<Energy,Energy> width(Energy, const ParticleData &) const;
 
  Energy partialWidth(int imode, Energy mH) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
protected:
 
  double HwDoubleBW(double x, double y) const;
 
  Complex HwW2(double tau) const;
 
  DecayMap branching(Energy scale, const ParticleData &) const;
 
protected:
 
  virtual void doinit();
 
 
private:
 
  SMHiggsWidthGenerator & operator=(const SMHiggsWidthGenerator &) = delete;
 
private:
 
  unsigned int widthopt_;
 
  double offshell_;
 
  Energy mw_;
 
  Energy mz_;
 
  Energy gamw_;
 
  Energy gamz_;
 
  vector<Energy> qmass_;
 
  vector<Energy> lmass_;
 
  double sw2_;
 
  double ca_;
 
  double cf_;
 
  mutable Energy qlast_;
 
  mutable Energy lambdaQCD_;
 
  mutable double alphaEM_;
 
  mutable double alphaS_;
 
  mutable double cd_;
 
  mutable Energy2 gfermiinv_;
 
  mutable double gam0_;
 
  mutable double beta0_;
 
  map<int,int> locMap_;
};
 
}
 
#endif /* HERWIG_SMHiggsWidthGenerator_H */