JetFinder.h

// -*- C++ -*-
//
// JetFinder.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad
// Copyright (C) 2009-2019 Simon Platzer
//
// 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_JetFinder_H
#define THEPEG_JetFinder_H
//
// This is the declaration of the JetFinder class.
//

#include "ThePEG/Interface/Interfaced.h"
#include "ThePEG/PDT/MatcherBase.h"
#include "Cuts.fh"

namespace ThePEG {

class JetFinder: public Interfaced {

public:

  JetFinder();

  virtual ~JetFinder();

public:

  virtual bool cluster(tcPDVector & ptype, vector<LorentzMomentum> & p,
                       tcCutsPtr parent, tcPDPtr t1 = tcPDPtr(),
                       tcPDPtr t2 = tcPDPtr()) const = 0;

  Ptr<MatcherBase>::tptr unresolvedMatcher() const { return theUnresolvedMatcher; }

  void minOutgoing(unsigned int n) { theMinOutgoing = n; }

  unsigned int minOutgoing() const { return theMinOutgoing; }

  bool restrictConsitutents() const { return theRestrictConstituents; }

  void restrictConsitutents(bool on) { theRestrictConstituents = on; }

  const pair<double,double>& constituentRapidityRange() const { return theConstituentRapidityRange; }

  void constituentRapidityRange(const pair<double,double>& r) { theConstituentRapidityRange = r; }

  virtual void describe() const = 0;

public:

  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

  static void Init();


// If needed, insert declarations of virtual function defined in the
// InterfacedBase class here (using ThePEG-interfaced-decl in Emacs).


private:

  string doYRange(string);

  Ptr<MatcherBase>::ptr theUnresolvedMatcher;

  unsigned int theMinOutgoing;

  bool theRestrictConstituents;

  pair<double,double> theConstituentRapidityRange;

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

};

}

#endif /* THEPEG_JetFinder_H */