/usr/include/Rivet/Projections/FoxWolframMoments.hh is in librivet-dev 1.8.3-2build1.
This file is owned by root:root, with mode 0o644.
The actual contents of the file can be viewed below.
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 | // -*- C++ -*-
#ifndef RIVET_FoxWolframMoments_HH
#define RIVET_FoxWolframMoments_HH
#include "Rivet/Rivet.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 {
/// @brief Calculate Fox-Wolfram moments
class FoxWolframMoments : public Projection {
public:
/// Constructor.
FoxWolframMoments(const FinalState& fsp)
{
setName("FoxWolframMoments");
addProjection(fsp, "FS");
/// @todo Let the user supply any projection they like?
VisibleFinalState vfs(fsp);
addProjection(vfs, "VFS");
// Initialize moments vector
for (int i = 0; i < MAXMOMENT ; ++i) {
_fwmoments.push_back(0.0);
}
}
/// Clone on the heap.
virtual const Projection* clone() const {
return new FoxWolframMoments(*this);
}
public:
/// The projected Fox-Wolfram Moment of order l
double getFoxWolframMoment(unsigned int l) const {
if (l < MAXMOMENT) {
return _fwmoments[l];
}
/// @todo What?!?
return -666.0;
}
protected:
/// Apply the projection to the event.
void project(const Event& e);
/// Compare projections.
int compare(const Projection& p) const;
private:
vector<double> _fwmoments;
};
}
#endif
|