PairRapidityCut.h

// -*- C++ -*-
//
// PairRapidityCut.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_PairRapidityCut_H
#define Herwig_PairRapidityCut_H
//
// This is the declaration of the PairRapidityCut class.
//
 
#include "ThePEG/Cuts/TwoCutBase.h"
#include "ThePEG/PDT/MatcherBase.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class PairRapidityCut: public TwoCutBase {
 
public:
 
  PairRapidityCut();
 
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:
 
//   /**
//    * Return the minimal allowed rapidity
//    */
//   double minRapidity() const { return theMinRapidity; }
// 
//   /**
//    * Return the maximal allowed rapidity
//    */
//   double maxRapidity() const { return theMaxRapidity; }
 
  const vector<pair<double,double> >& yRanges() const { return theYRanges; }
 
  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:
 
  string doYRange(string);
 
  vector<pair<double,double> > theYRanges;
 
  bool thePseudo;
 
  bool theSameFlavourOnly;
 
  bool theOppositeSignOnly;
 
  Ptr<MatcherBase>::ptr theFirstMatcher;
  Ptr<MatcherBase>::ptr theSecondMatcher;
 
private:
 
  PairRapidityCut & operator=(const PairRapidityCut &) = delete;
 
};
 
}
 
#endif /* Herwig_PairRapidityCut_H */