/cvmfs/atlas.cern.ch/repo/sw/ASG/AnalysisBase/2.4.31/JetRec/JetRec/JetFinder.h

// JetFinder.h

#ifndef JetFinder_H
#define JetFinder_H

// David Adams
// January 2014
//
// Tool to find jets using fastjet.
// The output collection is of type xAOD::JetContainer.

#include "AsgTools/AsgTool.h"
#include "JetInterface/IJetFinder.h"
#include "JetInterface/IJetFromPseudojet.h"
#include "fastjet/PseudoJet.hh"
#include "fastjet/JetDefinition.hh"
#include "xAODJet/JetContainer.h"
#include "AsgTools/ToolHandle.h"

namespace fastjet {
  class ClusterSequence;
}



class JetFinder
: public asg::AsgTool,
  virtual public IJetFinder {
ASG_TOOL_CLASS(JetFinder, IJetFinder)

public:

  // Ctor.
  JetFinder(std::string name);

  // Initialization.
  StatusCode initialize();

  // Find jets and put them in a container.
  int find(const PseudoJetVector& inps, xAOD::JetContainer& jets,
           xAOD::JetInput::Type contype, const NameList& ghostlabs) const;

  // Save a cluster sequence in the event store.
  // So it lives as long as this event.
  void save(fastjet::ClusterSequence* pcs) const;

  // Is variable-finding done?
  bool isVariableR() const;

  // Dump to log.
  void print() const;

private:  //data

  // Job options.
  std::string m_jetalg;                 // Kt, AntiKt, CamKt
  float m_jetrad;                       // Jet size parameter.
  float m_minrad;                       // Variable-R min radius
  float m_massscale;                    // Variable-R mass scale
  float m_ptmin;                        // pT min in MeV
  float m_ghostarea;                    // Area for ghosts. 0==>no ghosts.
  int m_ranopt;                         // Rand option: 0=fj default, 1=run/event
  ToolHandle<IJetFromPseudojet> m_bld;  // Tool to build jets.

  // Data
  fastjet::JetAlgorithm m_fjalg;
  bool m_isVariableR;

};

#endif

---