Rivet  1.8.3
SVertex.hh
1 // -*- C++ -*-
2 #ifndef RIVET_SVertex_HH
3 #define RIVET_SVertex_HH
4 
5 #include "Rivet/Rivet.hh"
6 #include "Rivet/Projection.hh"
7 #include "Rivet/Projections/PVertex.hh"
8 #include "Rivet/Projections/ChargedFinalState.hh"
9 #include "Rivet/Event.hh"
10 
11 namespace Rivet {
12 
13 
40  class SVertex : public Projection {
41  public:
42 
44 
45 
46 
48  const vector<FourMomentum>& jetaxes, double deltaR,
49  double detEta, double IPres, double DLS, double DLSres=0.0)
50  : _jetaxes(jetaxes), _deltaR(deltaR),
51  _detEta(detEta), _IPres(IPres), _DLS(DLS),
52  _DLSres(DLSres)
53  {
54  setName("SVertex");
55  addProjection(PVertex(), "PV");
56  addProjection(chfs, "FS");
57  if (_DLSres == 0.0) {
58  _DLSres = _IPres;
59  }
60  }
61 
63  virtual const Projection* clone() const {
64  return new SVertex(*this);
65  }
67 
68 
69  public:
71  const vector<FourMomentum>& getTaggedJets() const {
72  return _taggedjets;
73  }
74 
75  protected:
76 
78  void project(const Event& e);
79 
81  int compare(const Projection& p) const;
82 
83  private:
84 
86  const vector<FourMomentum>& _jetaxes;
87 
89  double _deltaR;
90 
93  //bool (*_applyVtxTrackCuts) (const ParticleVector&, const Vector3&, FourMomentum);
94  bool _applyVtxTrackCuts(const ParticleVector&, const Vector3&, FourMomentum);
95 
97  double _detEta;
98 
100  double _IPres;
101 
103  double _DLS;
104 
106  double _DLSres;
107 
109  vector<FourMomentum> _taggedjets;
110  };
111 
112 }
113 
114 #endif