FastJetFinder.h

// -*- C++ -*-
//
// FastJetFinder.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_FastJetFinder_H
#define THEPEG_FastJetFinder_H
//
// This is the declaration of the FastJetFinder class.
//

#include "ThePEG/Cuts/JetFinder.h"

namespace ThePEG {

class FastJetFinder: public JetFinder {

public:

  FastJetFinder();

  virtual ~FastJetFinder();

public:

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

  virtual void describe() const;

public:
  virtual void doinit();

  virtual void doinitrun();

  void persistentOutput(PersistentOStream & os) const;

  void persistentInput(PersistentIStream & is, int version);

  static void Init();

protected:

  virtual IBPtr clone() const;

  virtual IBPtr fullclone() const;


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

private:

  Energy2 theDCut;

  double theConeRadius;

  enum variants {
    kt = 1,
    CA = 2,
    antiKt = 3,
    sphericalKt = 4,
    sphericalCA = 5,
    sphericalAntiKt = 6
  };

  int theVariant;

  enum modes {
    inclusive = 1,
    exclusive = 2
  };

  int theMode;

  enum recombinations {
    recoPt = 1,
    recoE = 2
  };

  int theRecombination;

private:

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

};

}

#endif /* THEPEG_FastJetFinder_H */