1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/FastJets.hh"
 5
 6namespace Rivet {
 7
 8
 9  class ATLAS_2011_I889546 : public Analysis {
10  public:
11
12    /// Constructor
13    RIVET_DEFAULT_ANALYSIS_CTOR(ATLAS_2011_I889546);
14
15
16    /// @name Analysis methods
17    /// @{
18
19    /// Book histograms and initialise projections before the run
20    void init() {
21      /// Initialise and register projections
22      declare(FastJets(FinalState(), JetAlg::ANTIKT, 0.6), "AntiKtJets06");
23
24      /// Book histograms
25      book(_h_deltaPhi, {110., 160., 210., 260., 310., 400., 500., 600., 800., 10000.});
26      for (auto& b : _h_deltaPhi->bins()) {
27        book(b, 1, 1, b.index());
28      }
29    }
30
31
32    /// Perform the per-event analysis
33    void analyze(const Event& event) {
34
35      Jets jets06 = apply<FastJets>(event, "AntiKtJets06").jetsByPt(Cuts::absrap < 2.8 && Cuts::pT > 100*GeV);
36      if (jets06.size()>1){
37        if (jets06[0].absrap()<0.8 && jets06[1].absrap()<0.8) {
38          double observable = mapAngle0ToPi(jets06[0].phi()-jets06[1].phi()) / M_PI;
39          _h_deltaPhi->fill(jets06[0].pT(), observable);
40        }
41      }
42    }
43
44
45    /// Normalise histograms etc., after the run
46    void finalize() {
47      normalize(_h_deltaPhi, 1.0/M_PI);
48    }
49
50    /// @}
51
52
53  private:
54
55    /// Histograms
56    Histo1DGroupPtr _h_deltaPhi;
57
58  };
59
60
61
62  RIVET_DECLARE_ALIASED_PLUGIN(ATLAS_2011_I889546, ATLAS_2011_S8971293);
63
64}