/usr/include/Rivet/Projections/FParameter.hh is in librivet-dev 1.8.3-1.1.
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 79 80 81 82 83 84 85 86 87 88 89 | // -*- C++ -*-
#ifndef RIVET_Sphericity_HH
#define RIVET_FParameter_HH
#include "Rivet/Projection.hh"
#include "Rivet/Projections/FinalState.hh"
#include "Rivet/Event.hh"
namespace Rivet {
class FParameter : public Projection {
public:
/// @name Constructors etc.
//@{
/// Constructor
FParameter(const FinalState& fsp);
/// Clone on the heap.
virtual const Projection* clone() const {
return new FParameter(*this);
}
//@}
protected:
/// Perform the projection on the Event
void project(const Event& e);
/// Compare with other projections
//int compare(const Projection& p) const;
// Taken from Thrust.hh
int compare(const Projection& p) const {
return mkNamedPCmp(p, "FS");
}
public:
/// Reset the projection
void clear();
/// @name Access the event shapes by name
/// @{
/// F-Parametr
double F() const { return lambda1() >= lambda2() ? lambda2()/lambda1() : lambda1()/lambda2(); }
/// @}
/// @name Access the linearised transverse momentum tensor eigenvalues
/// @{
double lambda1() const { return _lambdas[0]; }
double lambda2() const { return _lambdas[1]; }
/// @}
/// @name Direct methods
/// Ways to do the calculation directly, without engaging the caching system
//@{
/// Manually calculate the sphericity, without engaging the caching system
void calc(const FinalState& fs);
/// Manually calculate the sphericity, without engaging the caching system
void calc(const vector<Particle>& fsparticles);
/// Manually calculate the sphericity, without engaging the caching system
void calc(const vector<FourMomentum>& fsmomenta);
/// Manually calculate the sphericity, without engaging the caching system
void calc(const vector<Vector3>& fsmomenta);
//@}
private:
/// Eigenvalues.
vector<double> _lambdas;
private:
/// Actually do the calculation
void _calcFParameter(const vector<Vector3>& fsmomenta);
};
}
#endif
|