EventShapes.hh

Go to the documentation of this file.

// -*- C++ -*-
//
// EventShapes.h is a part of Herwig++ - A multi-purpose Monte Carlo
// event generator Copyright (C) 2002-2011 The Herwig Collaboration
//
// Herwig++ is licenced under version 2 of the GPL, see COPYING for
// details. Please respect the MCnet academic guidelines, see
// GUIDELINES for details.
//

// Jan Hajer - replaced ThePEG classes with BoCA classes

#pragma once

#include "boca/math/Matrix3.hh"
#include "boca/fastjet/Jet.hh"
#include "boca/event/HemisphereMasses.hh"

namespace boca
{

class EventShapes
{

public:

    EventShapes();

    EventShapes(std::vector<Jet> const& jets);

    EventShapes(std::vector<LorentzVector<Momentum>> const& lorentz_vectors);

    double Thrust() const;

    double ThrustMajor() const;

    double ThrustMinor() const;

    double Oblateness() const;

    Vector3<double> ThrustAxis() const;

    Vector3<double> MajorAxis() const;

    Vector3<double> MinorAxis() const;

    double CParameter() const;

    double DParameter() const;

    std::vector<double> LinTenEigenValues() const;

    std::vector<Vector3<double>> LinTenEigenVectors() const;

    double Sphericity() const;

    double Aplanarity() const;

    double Planarity() const;

    Vector3<double> SphericityAxis() const;

    std::vector<double> SphericityEigenValues() const;

    std::vector<Vector3<double>> SphericityEigenVectors() const;

    boca::HemisphereMasses HemisphereMasses() const;

    Momentum PtInT(LorentzVector<Momentum> const& lorentz_vector) const;

    Momentum PtOutT(LorentzVector<Momentum> const& lorentz_vector) const;

    Angle RapidityT(LorentzVector<Momentum> const& lorentz_vector) const;

    Momentum PtInS(LorentzVector<Momentum> const& lorentz_vector) const;

    Momentum PtOutS(LorentzVector<Momentum> const& lorentz_vector) const;

    Angle RapidityS(LorentzVector<Momentum> const& lorentz_vector) const;

private:

    Momentum ScalarMomentum() const;

    Array3<GradedVector3<double>> Thrusts() const;

    Array3<GradedVector3<double>> LinearTensors() const;

    Array3<GradedVector3<double>> SphericalTensors() const;

    Array3<GradedVector3<double>> Thrusts1() const;

    Array3<GradedVector3<double>> Thrusts2() const;

    Array3<GradedVector3<double>> Thrusts3() const;

    Array3<GradedVector3<double>> Thrusts4() const;

    GradedVector3<Momentum> Thrust(std::vector<Vector3<Momentum>> const& vectors) const;

    GradedVector3<Momentum> Major(std::vector<Vector3<Momentum>> const& vectors) const;

    std::vector<LorentzVector<Momentum>> LorentzVectors() const;

    std::vector<Vector3<Momentum>> Vectors() const;

    std::vector<LorentzVector<Momentum>> lorentz_vectors_;

    Mutable<std::vector<Vector3<Momentum>>> vectors_;

    Mutable<Momentum> scalar_momentum_;

    Mutable<boca::HemisphereMasses> hemishpere_masses_;

    Mutable<Array3<GradedVector3<double>>> thrusts_;

    Mutable<Array3<GradedVector3<double>>> spherical_tensors_;

    Mutable<Array3<GradedVector3<double>>> linear_tensors_;

};

}