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 */