Qjets.h

Go to the documentation of this file.

#ifndef RecoJets_JetAlgorithms_QJets_h
#define RecoJets_JetAlgorithms_QJets_h
#include <queue>
#include <vector>
#include <list>
#include <algorithm>
#include "fastjet/JetDefinition.hh"
#include "fastjet/PseudoJet.hh"
#include "fastjet/ClusterSequence.hh"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/isFinite.h"
#include "FWCore/Utilities/interface/RandomNumberGenerator.h"
#include "FWCore/ServiceRegistry/interface/Service.h"
#include "CLHEP/Random/RandomEngine.h"

struct JetDistance {
  double dij;
  int j1;
  int j2;
};

class JetDistanceCompare {
public:
  JetDistanceCompare(){};
  bool operator()(const JetDistance& lhs, const JetDistance& rhs) const { return lhs.dij > rhs.dij; };
};

class Qjets {
private:
  double omega;
  bool _rand_seed_set;
  unsigned int _seed;
  double _zcut, _dcut, _dcut_fctr, _exp_min, _exp_max, _rigidity, _truncation_fctr;
  std::map<int, bool> _merged_jets;
  std::priority_queue<JetDistance, std::vector<JetDistance>, JetDistanceCompare> _distances;
  CLHEP::HepRandomEngine* _rnEngine;

  double d_ij(const fastjet::PseudoJet& v1, const fastjet::PseudoJet& v2) const;
  void computeDCut(fastjet::ClusterSequence& cs);

  double Rand();
  bool Prune(JetDistance& jd, fastjet::ClusterSequence& cs);
  bool JetsUnmerged(const JetDistance& jd) const;
  bool JetUnmerged(int num) const;
  void ComputeNewDistanceMeasures(fastjet::ClusterSequence& cs, unsigned int new_jet);
  void ComputeAllDistances(const std::vector<fastjet::PseudoJet>& inp);
  double ComputeMinimumDistance();
  double ComputeNormalization(double dmin);
  JetDistance GetNextDistance();
  bool Same(const JetDistance& lhs, const JetDistance& rhs);

public:
  Qjets(double zcut,
        double dcut_fctr,
        double exp_min,
        double exp_max,
        double rigidity,
        double truncation_fctr,
        CLHEP::HepRandomEngine* rnEngine)
      : _rand_seed_set(false),
        _zcut(zcut),
        _dcut(-1.),
        _dcut_fctr(dcut_fctr),
        _exp_min(exp_min),
        _exp_max(exp_max),
        _rigidity(rigidity),
        _truncation_fctr(truncation_fctr),
        _rnEngine(rnEngine){};

  void Cluster(fastjet::ClusterSequence& cs);
  void SetRandSeed(unsigned int seed); /* In case you want reproducible behavior */
};

class QjetsBaseExtras : public fastjet::ClusterSequence::Extras {
public:
  QjetsBaseExtras() : _wij(-1.) {}
  ~QjetsBaseExtras() override {}
  virtual double weight() const { return _wij; }
  friend class Qjets;

protected:
  double _wij;
};

#endif