InvariantMassCut.h

// -*- C++ -*-
//
// InvariantMassCut.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_InvariantMassCut_H
#define Herwig_InvariantMassCut_H
//
// This is the declaration of the InvariantMassCut class.
//
 
#include "ThePEG/Cuts/TwoCutBase.h"
#include "ThePEG/PDT/MatcherBase.h"
 
namespace Herwig {
 
using namespace ThePEG;
 
class InvariantMassCut: public TwoCutBase {
 
public:
 
  InvariantMassCut() : 
    theMinMass(0*GeV), theMaxMass(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 minMass() const { return theMinMass; }
 
  Energy maxMass() const { return theMaxMass; }
 
  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 theMinMass;
 
  Energy theMaxMass;
 
  bool theSameFlavourOnly;
 
  bool theOppositeSignOnly;
 
  Ptr<MatcherBase>::ptr theFirstMatcher;
  Ptr<MatcherBase>::ptr theSecondMatcher;
 
private:
 
  InvariantMassCut & operator=(const InvariantMassCut &) = delete;
 
};
 
}
 
#endif /* Herwig_InvariantMassCut_H */