/usr/include/Rivet/Projections/LossyFinalState.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 | // -*- C++ -*-
#ifndef RIVET_LossyFinalState_HH
#define RIVET_LossyFinalState_HH
#include "Rivet/Tools/Logging.hh"
#include "Rivet/Rivet.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Event.hh"
#include "Rivet/Projection.hh"
#include "Rivet/Projections/FinalState.hh"
namespace Rivet {
/// @brief Templated FS projection which can lose some of the supplied particles.
template <typename FILTER>
class LossyFinalState : public FinalState {
public:
/// @name Constructors
//@{
/// Constructor from FinalState.
LossyFinalState(const FinalState& fsp, FILTER filter)
: _filter(filter)
{
setName("LossyFinalState");
addProjection(fsp, "FS");
}
/// Stand-alone constructor. Initialises the base FinalState projection.
LossyFinalState(FILTER filter,
double mineta = -MAXRAPIDITY,
double maxeta = MAXRAPIDITY,
double minpt = 0.0)
: _filter(filter)
{
setName("LossyFinalState");
addProjection(FinalState(mineta, maxeta, minpt), "FS");
}
/// Virtual destructor, to allow subclassing
virtual ~LossyFinalState() { }
/// Clone on the heap.
virtual const Projection* clone() const {
return new LossyFinalState<FILTER>(*this);
}
//@}
protected:
/// Apply the projection on the supplied event.
void project(const Event& e) {
const FinalState& fs = applyProjection<FinalState>(e, "FS");
getLog() << Log::DEBUG << "Pre-loss number of FS particles = " << fs.particles().size() << endl;
_theParticles.clear();
std::remove_copy_if(fs.particles().begin(), fs.particles().end(),
std::back_inserter(_theParticles), _filter);
getLog() << Log::DEBUG << "Filtered number of FS particles = " << _theParticles.size() << endl;
}
/// Compare projections.
int compare(const Projection& p) const {
const LossyFinalState<FILTER>& other = pcast< LossyFinalState<FILTER> >(p);
const int fscmp = mkNamedPCmp(other, "FS");
if (fscmp) return fscmp;
return _filter.compare(other._filter);
}
protected:
/// Filtering object: must support operator(const Particle&) and compare(const Filter&)
FILTER _filter;
};
}
#endif
|