PseudoJet.hh

Go to the documentation of this file.

#pragma once

#include <boost/operators.hpp>

#include "fastjet/PseudoJet.hh"
#include "boca/units/Prefixes.hh"

class TLorentzVector;

namespace boca
{

template<typename Value_>
class LorentzVector;

template<typename Value_>
class Vector3;

template<typename Value_>
class Vector2;

class PseudoJet : fastjet::PseudoJet
    , boost::totally_ordered<PseudoJet>
{

public:

    PseudoJet();

    PseudoJet(fastjet::PseudoJet const &pseudo_jet);

    PseudoJet(TLorentzVector const &vector);

    PseudoJet(boca::LorentzVector<Momentum> const &vector);

    PseudoJet(boca::Vector3<boca::Momentum> const& spatial, const boca::Energy &e);

    PseudoJet(Vector2<boca::Momentum> const& transverse, Momentum const& z, boca::Energy const& e);

    PseudoJet(Momentum const& x, Momentum const& y, Momentum const& z, boca::Energy const& e);

    PseudoJet(PseudoJet const&) = default;

    PseudoJet(PseudoJet &&) = default;


    void Reset(PseudoJet const &pseudo_jet);

    void ScaleMomentum(double factor);

    fastjet::PseudoJet &FastJet();

    fastjet::PseudoJet const &FastJet() const;


    virtual UserInfoBase const &Info() const;

    virtual UserInfoBase &Info();

    bool HasInfo();


    Momentum Pt() const;

    boca::Energy Energy() const;

    boca::Energy E() const;

    boca::Energy Et2() const;

    boca::Energy Et() const;

    boca::Mass Mass() const;

    boca::MassSquare MassSquare() const;

    Momentum Px() const;

    Momentum Py() const;

    Momentum Pz() const;

    MomentumSquare ModP2() const;

    Momentum ModP() const;


    Angle Rap() const;

    Angle Phi() const;

    Angle DeltaPhiTo(PseudoJet const &jet) const;

    Angle DeltaRapTo(PseudoJet const &jet) const;

    Angle DeltaRTo(PseudoJet const &jet) const;


    Vector2<Angle> DeltaTo(PseudoJet const &jet) const;

    Vector2<Angle> Angles(bool wrap_phi = false) const;

    Vector2<Angle> AnglesMinTo(PseudoJet const &jet) const;

    boca::LorentzVector< Momentum > LorentzVector() const;

    boca::Vector3<Momentum> Spatial() const;

    Vector2<Momentum> Transverse() const;


    bool operator<(PseudoJet const &jet) const;

    bool operator==(const PseudoJet &pseudo_jet) const;

    PseudoJet &operator=(PseudoJet const&) & = default;

    PseudoJet &operator=(PseudoJet &&) & = default;

    friend std::ostream & operator<<(std::ostream & stream, PseudoJet const& jet);


    virtual ~PseudoJet() {}

protected:


};

}