O2AlphaS.h

// -*- C++ -*-
//
// O2AlphaS.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_O2AlphaS_H
#define HERWIG_O2AlphaS_H
//
// This is the declaration of the O2AlphaS class.
//
 
#include "ThePEG/StandardModel/AlphaSBase.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class O2AlphaS: public AlphaSBase {
 
public:
 
  O2AlphaS() : _lambdaQCD(180.*MeV), _bcoeff(6,0.0), _ccoeff(6,0.0),
               _lambdas(7), _threshold(6), _match(6,0.0), _copt(0) {}
 
  virtual double value(Energy2 scale, const StandardModelBase & sm) const;
 
  virtual vector<Energy2> flavourThresholds() const;
 
  virtual vector<Energy> LambdaQCDs() const;
 
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();
 
 
private:
 
  O2AlphaS & operator=(const O2AlphaS &) = delete;
 
private:
 
  Energy _lambdaQCD;
 
  vector<double> _bcoeff;
 
  vector<double> _ccoeff;
 
  vector<Energy> _lambdas;
 
  vector<Energy> _threshold;
  
  vector<double> _match;
  unsigned int _copt;
};
 
}
 
#endif /* HERWIG_O2AlphaS_H */