SMTopPOWHEGDecayer.h

// -*- C++ -*-
#ifndef Herwig_SMTopPOWHEGDecayer_H
#define Herwig_SMTopPOWHEGDecayer_H
//
// This is the declaration of the SMTopPOWHEGDecayer class.
//

#include "SMTopDecayer.h"

namespace Herwig {

using namespace ThePEG;

class SMTopPOWHEGDecayer: public SMTopDecayer {

public:

  SMTopPOWHEGDecayer();

  virtual POWHEGType hasPOWHEGCorrection() {return FSR;}

  virtual HardTreePtr generateHardest(ShowerTreePtr);

  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

  static void Init();

protected:

  virtual IBPtr clone() const;

  virtual IBPtr fullclone() const;


protected:
  bool deadZoneCheck(double xw, double xg);

protected:

  double matrixElementRatio(vector<Lorentz5Momentum> particleMomenta);

protected:

  bool calcMomenta(int j, Energy pT, double y, double phi, double& xg, 
                   double& xw, double& xb, double& xb_z, 
                   vector<Lorentz5Momentum>& particleMomenta);

protected:

  bool psCheck(double xg, double xw);

protected:

  vector<Lorentz5Momentum> hardMomenta();

private:

  SMTopPOWHEGDecayer & operator=(const SMTopPOWHEGDecayer &);

private:

  Energy mt_;

  double w_;

  double b_;

  double w2_;

  double b2_;

  Energy pTmin_;

  Energy pT_;

};

}

#endif /* Herwig_SMTopPOWHEGDecayer_H */