ClusterSequence.hh

Go to the documentation of this file.

#pragma once

#include "fastjet/ClusterSequence.hh"

#include "boca/units/ElectronVolt.hh"

namespace fastjet
{
class ClusterSequence;
class JetDefinition;
}

namespace boca
{

class Jet;

class ClusterSequence
{

public :

    ClusterSequence(std::vector<fastjet::PseudoJet> const& jets, fastjet::JetDefinition const& jet_definition);

    ClusterSequence(std::vector<Jet> const& jets, fastjet::JetDefinition const& jet_definition);

    ClusterSequence(ClusterSequence const& cluster_sequence);

    ClusterSequence(ClusterSequence && cluster_sequence) noexcept;

    ~ClusterSequence();

    ClusterSequence& operator=(ClusterSequence const& cluster_sequence);

    ClusterSequence& operator=(ClusterSequence && cluster_sequence) noexcept;

    std::vector<Jet> ExclusiveJets(int jet_number) const;

    std::vector<Jet> ExclusiveJetsUpTo(int sub_jet_number) const;

    std::vector<Jet> InclusiveJets(Momentum const& min_pt = at_rest) const;

    fastjet::ClusterSequence& Get();

    fastjet::ClusterSequence const& Get() const;

    void AddConstituents(fastjet::PseudoJet const& jet, std::vector<fastjet::PseudoJet>& subjet_vector) const;

    void NoLongerNeeded();

private :

    fastjet::ClusterSequence* cluster_sequence_;

    mutable bool success_ = true;

};

}