FoxWolframMoments.hh

// -*- C++ -*-
#ifndef RIVET_FoxWolframMoments_HH
#define RIVET_FoxWolframMoments_HH

#include "Rivet/Config/RivetCommon.hh"
#include "Rivet/Projection.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Projections/VetoedFinalState.hh"
#include "Rivet/Projections/VisibleFinalState.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Event.hh"

#include <gsl/gsl_sf_legendre.h>

#define MAXMOMENT 5

namespace Rivet {


  class FoxWolframMoments : public Projection {
  public:

    FoxWolframMoments(const FinalState& fsp) {
      setName("FoxWolframMoments");
      addProjection(fsp, "FS");

      VisibleFinalState vfs(fsp);
      addProjection(vfs, "VFS");

      // Initialize moments vector
      for (int i = 0; i < MAXMOMENT ; ++i) {
        _fwmoments.push_back(0.0);
      }
    }

    DEFAULT_RIVET_PROJ_CLONE(FoxWolframMoments);


    double getFoxWolframMoment(unsigned int l) const {
      if (l < MAXMOMENT) {
        return _fwmoments[l];
      }
      return -666.0;
    }


  protected:

    void project(const Event& e);

    int compare(const Projection& p) const;


  private:

    vector<double> _fwmoments;

  };


}


#endif