MECuts.h

// -*- C++ -*-
//
// MECuts.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2007 Leif Lonnblad
//
// ThePEG is licenced under version 2 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 &);

};

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

---