AlphaSBase.h

// -*- C++ -*-
//
// AlphaSBase.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_AlphaSBase_H
#define ThePEG_AlphaSBase_H
// This is the declaration of the AlphaSBase class.
 
#include "RunningCoupling.h"
 
namespace ThePEG {
 
class AlphaSBase: public RunningCoupling {
 
public:
 
  virtual vector<Energy2> flavourThresholds() const = 0;
 
  virtual vector<Energy> LambdaQCDs() const = 0;
 
  unsigned int Nf(Energy2 scale) const {
    unsigned int i = 0;
    while ( i < theFlavourThresholds.size() &&
            theFlavourThresholds[i] < scale ) ++i;
    return i;
  }
 
  Energy LambdaQCD(unsigned int nflav) const  { return theLambdaQCDs[nflav]; }
 
  Energy LambdaQCD(Energy2 scale) const { return theLambdaQCDs[Nf(scale)]; }
 
  const vector<Energy>& quarkMasses() const { return theQuarkMasses; }
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
 
  virtual void doinit();
 
 
private:
 
  vector<Energy> theQuarkMasses;
 
  vector<Energy2> theFlavourThresholds;
 
  vector<Energy> theLambdaQCDs;
 
private:
 
  static AbstractClassDescription<AlphaSBase> initAlphaSBase;
 
  AlphaSBase & operator=(const AlphaSBase &) = delete;
 
};
 
template <>
struct BaseClassTrait<AlphaSBase,1>: public ClassTraitsType {
  typedef RunningCoupling NthBase;
};
 
template <>
struct ClassTraits<AlphaSBase>: public ClassTraitsBase<AlphaSBase> {
  static string className() { return "ThePEG::AlphaSBase"; }
};
 
}
 
#endif /* ThePEG_AlphaSBase_H */