1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/GammaGammaFinalState.hh"
 4#include "Rivet/Projections/FastJets.hh"
 5
 6namespace Rivet {
 7
 8
 9  /// @brief Jet production in photon-photon collisions at 198 GeV
10  class L3_2004_I661114 : public Analysis {
11  public:
12
13    /// Constructor
14    RIVET_DEFAULT_ANALYSIS_CTOR(L3_2004_I661114);
15
16
17    /// @name Analysis methods
18    /// @{
19
20    /// Book histograms and initialise projections before the run
21    void init() {
22      // get the hadronic final state
23      const FinalState & fs = declare(GammaGammaFinalState(), "FS");
24      declare(FastJets(fs, JetAlg::KT,1.),"Jets");
25
26      // Book histograms
27      book(_h_y, 1, 1, 1);
28
29    }
30
31
32    /// Perform the per-event analysis
33    void analyze(const Event& event) {
34      Jets jets = apply<FastJets>(event, "Jets").jetsByPt(Cuts::pT > 3*GeV and Cuts::abseta < 1.0);
35      if(jets.empty()) vetoEvent;
36      for(const Jet & jet : jets) {
37      	_h_y->fill(jet.pT());
38      }
39    }
40
41
42    /// Normalise histograms etc., after the run
43    void finalize() {
44
45      scale(_h_y, crossSection()/picobarn/sumOfWeights());
46
47    }
48
49    /// @}
50
51
52    /// @name Histograms
53    /// @{
54    Histo1DPtr _h_y;
55    /// @}
56
57
58  };
59
60
61  RIVET_DECLARE_PLUGIN(L3_2004_I661114);
62
63
64}