SimpleDISCut.h

// -*- C++ -*-
//
// SimpleDISCut.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_SimpleDISCut_H
#define THEPEG_SimpleDISCut_H
//
// This is the declaration of the SimpleDISCut class.
//

#include "ThePEG/Cuts/TwoCutBase.h"

namespace ThePEG {

class SimpleDISCut: public TwoCutBase {

public:

  SimpleDISCut()
    : theMinQ2(1.0*GeV2), theMaxQ2(100.0*GeV2),theMiny(0.0), theMaxy(1.0),
      theMinW2(100.0*GeV2), theMaxW2(1000000.0*GeV2), chargedCurrent(false) {}

public:

  virtual Energy2 minSij(tcPDPtr pi, tcPDPtr pj) const;

  virtual Energy2 minTij(tcPDPtr pi, tcPDPtr po) const;

  virtual double minDeltaR(tcPDPtr pi, tcPDPtr pj) const;

  virtual Energy minKTClus(tcPDPtr pi, tcPDPtr pj) const;

  virtual double minDurham(tcPDPtr pi, tcPDPtr pj) const;

  virtual bool passCuts(tcCutsPtr parent, tcPDPtr pitype, tcPDPtr pjtype,
                        LorentzMomentum pi, LorentzMomentum pj,
                        bool inci = false, bool incj = false) const;

  virtual void describe() const;

protected:

  bool check(long idi, long ido) 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:

  double maxMiny() const;

  double minMaxy() const;

  Energy2 maxMinQ2() const;

  Energy2 minMaxQ2() const;

  Energy2 maxMinW2() const;

  Energy2 minMaxW2() const;

private:

  Energy2 theMinQ2;

  Energy2 theMaxQ2;

  double theMiny;

  double theMaxy;



  Energy2 theMinW2;

  Energy2 theMaxW2;

  bool chargedCurrent;

private:

  static ClassDescription<SimpleDISCut> initSimpleDISCut;

  SimpleDISCut & operator=(const SimpleDISCut &) = delete;

};

}

#include "ThePEG/Utilities/ClassTraits.h"

namespace ThePEG {

template <>
struct BaseClassTrait<SimpleDISCut,1> {
  typedef TwoCutBase NthBase;
};

template <>
struct ClassTraits<SimpleDISCut>
  : public ClassTraitsBase<SimpleDISCut> {
  static string className() { return "ThePEG::SimpleDISCut"; }
  static string library() { return "SimpleDISCut.so"; }
};

}

#endif /* THEPEG_SimpleDISCut_H */