MECuts.h

// -*- C++ -*-
//
// MECuts.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_MECuts_H
#define ThePEG_MECuts_H
// This is the declaration of the MECuts class.
 
#include "ThePEG/Interface/Interfaced.h"
#include "ThePEG/Vectors/LorentzRotation.fh"
#include "ThePEG/Utilities/Triplet.h"
#include "ThePEG/PDT/StandardMatchers.h"
 
namespace ThePEG {
 
class MECuts: public Interfaced {
 
public:
 
  MECuts();
 
public:
 
  void cut(const SubProcess &) const;
 
public:
 
  Energy mHatMin() const { return theMHatMin; }
  Energy mHatMax() const { 
    theMHatMax > mHatMin()? theMHatMax: Constants::MaxEnergy; 
  }
  Energy sHatMin() const { return sqr(mHatMin()); }
  Energy sHatMax() const { return sqr(mHatMax()); }
 
  Energy pTHatMin() const { return thePTHatMin; }
  Energy pTHatMax() const { 
    return thePTHatMax > pTHatMin()? thePTHatMax: Constants::MaxEnergy; 
  }
 
  Energy pTHatSingularMin() const { return thePTHatSingularMin; }
  Energy singularMassMax() const { return theSingularMassMax; }
 
  double cTHMin() const { return theCTHMin; }
  double cTHMax() const { return theCTHMax; }
 
  Energy2 tHatMin() const { return theTHatMin; }
  Energy2 tHatMax() const { 
    return theTHatMax > tHatMin()? theTHatMax: Constants::MaxEnergy2; 
  }
 
  Energy2 uHatMin() const { return theUHatMin; }
  Energy2 uHatMax() const { 
    return theUHatMax > uHatMin()? theUHatMax: Constants::MaxEnergy2; 
  }
 
  Energy2 scaleMin() const { return theScaleMin; }
  Energy2 scaleMax() const {
    return theScaleMax > scaleMin()? theScaleMax: Constants::MaxEnergy2; 
  }
 
public:
 
  void persistentOutput(PersistentOStream &) const;
  void persistentInput(PersistentIStream &, int);
 
  static void Init();
 
protected:
 
  virtual void newcut(const SubProcess &) const;
 
protected:
 
  virtual void doupdate();
 
  virtual IBPtr clone() const;
 
private:
 
  Energy theMHatMin;
  Energy theMHatMax;
 
  Energy thePTHatMin;
  Energy thePTHatMax;
 
  Energy thePTHatSingularMin;
  Energy theSingularMassMax;
 
  double theCTHMin;
  double theCTHMax;
 
  Energy2 theTHatMin;
  Energy2 theTHatMax;
 
  Energy2 theUHatMin;
  Energy2 theUHatMax;
 
  Energy2 theScaleMin;
  Energy2 theScaleMax;
 
private:
 
  static ClassDescription<MECuts> initMECuts;
 
  MECuts & operator=(const MECuts &) = delete;
 
};
 
template <>
struct BaseClassTrait<MECuts,1>: public ClassTraitsType {
  typedef Interfaced NthBase;
};
 
template <>
struct ClassTraits<MECuts>:
    public ClassTraitsBase<MECuts> {
  static string className() { return "ThePEG::MECuts"; }
};
 
}
 
#endif /* ThePEG_MECuts_H */