Rivet  1.8.0
ClusteredPhotons.hh
1 // -*- C++ -*-
2 #ifndef RIVET_ClusteredPhotons_HH
3 #define RIVET_ClusteredPhotons_HH
4 
5 #include "Rivet/Tools/Logging.hh"
6 #include "Rivet/Rivet.hh"
7 #include "Rivet/Particle.hh"
8 #include "Rivet/Event.hh"
9 #include "Rivet/Projection.hh"
10 #include "Rivet/Projections/FinalState.hh"
11 #include "Rivet/Projections/IdentifiedFinalState.hh"
12 
13 namespace Rivet {
14 
15 
17  class ClusteredPhotons : public FinalState {
18  public:
19 
21 
22 
23 
24  ClusteredPhotons(const FinalState& fs, const FinalState& signal, double dRmax)
25  : _dRmax(dRmax)
26  {
27  setName("ClusteredPhotons");
28  IdentifiedFinalState photonfs(fs);
29  photonfs.acceptId(PHOTON);
30  addProjection(photonfs, "Photons");
31  addProjection(signal, "Signal");
32  }
33 
34 
36  virtual const Projection* clone() const {
37  return new ClusteredPhotons(*this);
38  }
40 
41 
42  public:
43 
44  protected:
45 
47  void project(const Event& e);
48 
50  int compare(const Projection& p) const;
51 
52 
53  private:
54 
56  double _dRmax;
57 
58  };
59 
60 
61 }
62 
63 
64 #endif