IILightInvertedTildeKinematics.h

// -*- C++ -*-
//
// IILightInvertedTildeKinematics.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_IILightInvertedTildeKinematics_H
#define HERWIG_IILightInvertedTildeKinematics_H
//
// This is the declaration of the IILightInvertedTildeKinematics class.
//
 
#include "Herwig/MatrixElement/Matchbox/Phasespace/InvertedTildeKinematics.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class IILightInvertedTildeKinematics: public Herwig::InvertedTildeKinematics {
 
public:
 
  virtual bool doMap(const double *);
 
  virtual Energy lastPt() const;
 
  virtual double lastZ() const;
 
  virtual Energy ptMax() const;
 
  virtual pair<double,double> zBounds(Energy pt, Energy hardPt = ZERO) const;
 
  virtual bool doesTransform() const { return true; }
 
  virtual Lorentz5Momentum transform(const Lorentz5Momentum& p) const {
    return p-(2.*(KplusKtilde*p)/KplusKtilde2)*KplusKtilde+(2.*(Ktilde*p)/K2)*K;
  }
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
 
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
 
private:
 
  Lorentz5Momentum K;
  Energy2 K2;
  Lorentz5Momentum Ktilde;
  Lorentz5Momentum KplusKtilde;
  Energy2 KplusKtilde2;
 
private:
 
  IILightInvertedTildeKinematics & operator=(const IILightInvertedTildeKinematics &) = delete;
 
};
 
}
 
#endif /* HERWIG_IILightInvertedTildeKinematics_H */