PairPtCut.h

// -*- C++ -*-
//
// PairPtCut.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_PairPtCut_H
#define Herwig_PairPtCut_H
//
// This is the declaration of the PairPtCut class.
//
 
#include "ThePEG/Cuts/TwoCutBase.h"
#include "ThePEG/PDT/MatcherBase.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class PairPtCut: public TwoCutBase {
 
public:
 
  PairPtCut() : 
    theMinPt(0*GeV), theMaxPt(Constants::MaxEnergy),
    theSameFlavourOnly(false), theOppositeSignOnly(false) {}
 
public:
 
  virtual bool passCuts(tcCutsPtr parent, tcPDPtr pitype, tcPDPtr pjtype,
                        LorentzMomentum pi, LorentzMomentum pj,
                        bool inci = false, bool incj = false) const;
 
  virtual Energy2 minSij(tcPDPtr , tcPDPtr ) const { return ZERO; }
 
  virtual Energy2 minTij(tcPDPtr , tcPDPtr ) const { return ZERO; }
 
  virtual double minDeltaR(tcPDPtr , tcPDPtr ) const { return ZERO; }
 
  virtual Energy minKTClus(tcPDPtr , tcPDPtr ) const { return ZERO; }
 
  virtual double minDurham(tcPDPtr , tcPDPtr ) const { return ZERO; }
 
 
  virtual void describe() const;
 
public:
 
  Energy minPt() const { return theMinPt; }
 
  Energy maxPt() const { return theMaxPt; }
 
  bool sameFlavourOnly() const { return theSameFlavourOnly; }
 
  bool oppositeSignOnly() const { return theOppositeSignOnly; }
 
  Ptr<MatcherBase>::tptr firstMatcher() const { return theFirstMatcher; }
  Ptr<MatcherBase>::tptr secondMatcher() const { return theSecondMatcher; }
 
protected:
 
  int family(long id) const;
 
public:
 
  void persistentOutput(PersistentOStream & os) const;
 
  void persistentInput(PersistentIStream & is, int version);
 
  static void Init();
 
protected:
 
  virtual IBPtr clone() const;
 
  virtual IBPtr fullclone() const;
 
private:
 
  Energy theMinPt;
 
  Energy theMaxPt;
 
  bool theSameFlavourOnly;
 
  bool theOppositeSignOnly;
 
  Ptr<MatcherBase>::ptr theFirstMatcher;
  Ptr<MatcherBase>::ptr theSecondMatcher;
 
private:
 
  PairPtCut & operator=(const PairPtCut &) = delete;
 
};
 
}
 
#endif /* Herwig_PairPtCut_H */