Rivet  3.1.2
LossyFinalState.hh
1 // -*- C++ -*-
2 #ifndef RIVET_LossyFinalState_HH
3 #define RIVET_LossyFinalState_HH
4 
5 #include "Rivet/Tools/Logging.hh"
6 #include "Rivet/Config/RivetCommon.hh"
7 #include "Rivet/Particle.hh"
8 #include "Rivet/Event.hh"
9 #include "Rivet/Projection.hh"
10 #include "Rivet/Projections/FinalState.hh"
11 
12 namespace Rivet {
13 
14 
16  template <typename FILTER>
17  class LossyFinalState : public FinalState {
18  public:
19 
21 
22 
24  LossyFinalState(const FinalState& fsp, FILTER filter)
25  : _filter(filter)
26  {
27  setName("LossyFinalState");
28  declare(fsp, "FS");
29  }
30 
32  LossyFinalState(FILTER filter, const Cut& c=Cuts::open())
33  : _filter(filter)
34  {
35  setName("LossyFinalState");
36  declare(FinalState(c), "FS");
37  }
38 
40  virtual ~LossyFinalState() { }
41 
44 
46 
47 
49  void project(const Event& e) {
50  const FinalState& fs = applyProjection<FinalState>(e, "FS");
51  getLog() << Log::DEBUG << "Pre-loss number of FS particles = " << fs.particles().size() << '\n';
52  _theParticles.clear();
53  std::remove_copy_if(fs.particles().begin(), fs.particles().end(),
54  std::back_inserter(_theParticles), _filter);
55  getLog() << Log::DEBUG << "Filtered number of FS particles = " << _theParticles.size() << '\n';
56  }
57 
58 
60  CmpState compare(const Projection& p) const {
61  const LossyFinalState<FILTER>& other = pcast< LossyFinalState<FILTER> >(p);
62  const CmpState fscmp = mkNamedPCmp(other, "FS");
63  if (fscmp != CmpState::EQ) return fscmp;
64  return _filter.compare(other._filter);
65  }
66 
67 
68  protected:
69 
71  FILTER _filter;
72 
73  };
74 
75 
76 }
77 
78 #endif
void setName(const std::string &name)
Used by derived classes to set their name.
Definition: Projection.hh:142
Definition: MC_Cent_pPb.hh:10
FinalState(const Cut &c=Cuts::open())
Construction using Cuts object.
virtual ~LossyFinalState()
Virtual destructor, to allow subclassing.
Definition: LossyFinalState.hh:40
Templated FS projection which can lose some of the supplied particles.
Definition: LossyFinalState.hh:17
virtual const Particles & particles() const
Get the particles in no particular order, with no cuts.
Definition: ParticleFinder.hh:49
const Cut & open()
Fully open cut singleton, accepts everything.
void project(const Event &e)
Apply the projection on the supplied event.
Definition: LossyFinalState.hh:49
LossyFinalState(FILTER filter, const Cut &c=Cuts::open())
Stand-alone constructor. Initialises the base FinalState projection.
Definition: LossyFinalState.hh:32
CmpState compare(const Projection &p) const
Compare projections.
Definition: LossyFinalState.hh:60
Representation of a HepMC event, and enabler of Projection caching.
Definition: Event.hh:22
LossyFinalState(const FinalState &fsp, FILTER filter)
Constructor from FinalState.
Definition: LossyFinalState.hh:24
Project out all final-state particles in an event. Probably the most important projection in Rivet! ...
Definition: FinalState.hh:12
const PROJ & declare(const PROJ &proj, const std::string &name)
Register a contained projection (user-facing version)
Definition: ProjectionApplier.hh:170
Cmp< Projection > mkNamedPCmp(const Projection &otherparent, const std::string &pname) const
Base class for all Rivet projections.
Definition: Projection.hh:29
Log & getLog() const
Get a Log object based on the getName() property of the calling projection object.
Definition: Projection.hh:136
DEFAULT_RIVET_PROJ_CLONE(LossyFinalState)
Clone on the heap.