InvertedTildeKinematics.h

// -*- C++ -*-
//
// InvertedTildeKinematics.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_InvertedTildeKinematics_H
#define HERWIG_InvertedTildeKinematics_H
//
// This is the declaration of the InvertedTildeKinematics class.
//
 
#include "ThePEG/Handlers/HandlerBase.h"
#include "ThePEG/Handlers/StandardXComb.h"
#include "ThePEG/Repository/EventGenerator.h"
#include "Herwig/MatrixElement/Matchbox/Dipoles/SubtractionDipole.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class InvertedTildeKinematics: public HandlerBase {
 
public:
 
  InvertedTildeKinematics();
 
public:
 
  const Lorentz5Momentum& realEmitterMomentum() const { return theRealEmitterMomentum; }
 
  const Lorentz5Momentum& realEmissionMomentum() const { return theRealEmissionMomentum; }
 
  const Lorentz5Momentum& realSpectatorMomentum() const { return theRealSpectatorMomentum; }
 
  const Lorentz5Momentum& bornEmitterMomentum() const { 
    return theBornXComb->meMomenta()[theDipole->bornEmitter()];
  }
 
  const Lorentz5Momentum& bornSpectatorMomentum() const { 
    return theBornXComb->meMomenta()[theDipole->bornSpectator()];
  }
 
  double emitterX() const { 
    return 
      theDipole->bornEmitter() == 0 ?
      theBornXComb->lastX1() :
      theBornXComb->lastX2();
  }
 
  double spectatorX() const { 
    return 
      theDipole->bornSpectator() == 0 ?
      theBornXComb->lastX1() :
      theBornXComb->lastX2();
  }
 
  const vector<double>& subtractionParameters() const { return theDipole->subtractionParameters(); }
 
  virtual bool doesTransform() const { return false; }
 
  virtual Lorentz5Momentum transform(const Lorentz5Momentum& p) const { return p; }
 
  Energy2 sHat() const { return theBornXComb->lastSHat(); }
 
public:
 
  Ptr<InvertedTildeKinematics>::ptr cloneMe() const {
    return dynamic_ptr_cast<Ptr<InvertedTildeKinematics>::ptr>(clone());
  }
 
  void prepare(tcStdXCombPtr newRealXComb,
               tcStdXCombPtr newBornXComb) {
    theRealXComb = newRealXComb; theBornXComb = newBornXComb;
  }
 
  tcStdXCombPtr realXComb() const { return theRealXComb; }
 
  tcStdXCombPtr bornXComb() const { return theBornXComb; }
 
  void dipole(Ptr<SubtractionDipole>::tptr dip) { theDipole = dip; }
 
  Ptr<SubtractionDipole>::tptr dipole() { return theDipole; }
 
  Ptr<SubtractionDipole>::tcptr dipole() const { return theDipole; }
 
  virtual int nDimRadiation() const { return 3; }
 
  virtual bool doMap(const double *) = 0;
 
  void ptCut(Energy pt) { thePtCut = pt; }
  
  Energy ptCut() const { return thePtCut; }
 
  virtual int evolutionVariable() const { return 0; }
 
  virtual double evolutionCutoff() const { return 0.0; }
 
  virtual Energy lastPt() const = 0;
 
  virtual double lastZ() const = 0;
 
  virtual Energy lastScale() const;
 
  virtual Energy ptMax() const = 0;
 
  virtual pair<double,double> zBounds(Energy pt, Energy hardPt = ZERO) const = 0;
 
  virtual pair<Energy,double> generatePtZ(double& jac, const double * r,
                                          double power=1., vector<double>* values = NULL) const;
 
  double jacobian() const { return theJacobian; }
 
  cPDPtr realEmitterData() const { 
    return 
      (theDipole && theRealXComb) ? 
      theRealXComb->mePartonData()[theDipole->realEmitter()] :
      cPDPtr();
  }
 
  cPDPtr realEmissionData() const { 
    return 
      (theDipole && theRealXComb) ? 
      theRealXComb->mePartonData()[theDipole->realEmission()] :
      cPDPtr();
  }
 
  cPDPtr realSpectatorData() const { 
    return 
      (theDipole && theRealXComb) ? 
      theRealXComb->mePartonData()[theDipole->realSpectator()] :
      cPDPtr();
  }
 
  cPDPtr bornEmitterData() const { 
    return 
      (theDipole && theBornXComb) ? 
      theBornXComb->mePartonData()[theDipole->bornEmitter()] :
      cPDPtr();
  }
 
  cPDPtr bornSpectatorData() const { 
    return 
      (theDipole && theBornXComb) ? 
      theBornXComb->mePartonData()[theDipole->bornSpectator()] :
      cPDPtr();
  }
 
protected:
 
  Lorentz5Momentum& realEmitterMomentum() { return theRealEmitterMomentum; }
 
  Lorentz5Momentum& realEmissionMomentum() { return theRealEmissionMomentum; }
 
  Lorentz5Momentum& realSpectatorMomentum() { return theRealSpectatorMomentum; }
 
  vector<double>& subtractionParameters() { return theDipole->subtractionParameters(); }
 
  void jacobian(double w) { theJacobian = w; }
 
  Lorentz5Momentum getKt(const Lorentz5Momentum& p1,
                         const Lorentz5Momentum& p2,
                         Energy pt,
                         double phi,
                         bool spacelike = false) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
 
// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).
 
protected:
 
  virtual void rebind(const TranslationMap & trans);
 
  virtual IVector getReferences();
 
 
private:
 
  Ptr<SubtractionDipole>::tptr theDipole;
 
  tcStdXCombPtr theRealXComb;
 
  tcStdXCombPtr theBornXComb;
 
  Lorentz5Momentum theRealEmitterMomentum;
 
  Lorentz5Momentum theRealEmissionMomentum;
 
  Lorentz5Momentum theRealSpectatorMomentum;
 
  double theJacobian;
 
  Energy thePtCut;
 
private:
 
  InvertedTildeKinematics & operator=(const InvertedTildeKinematics &) = delete;
 
};
 
}
 
#endif /* HERWIG_InvertedTildeKinematics_H */