1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/ZFinder.hh"
 4
 5namespace Rivet {
 6
 7
 8  /// @brief Measurement of D0 Run II Z \f$ p_\perp \f$ diff cross-section shape
 9  ///
10  /// @author Andy Buckley
11  /// @author Gavin Hesketh
12  /// @author Frank Siegert
13  class D0_2007_S7075677 : public Analysis {
14  public:
15
16    RIVET_DEFAULT_ANALYSIS_CTOR(D0_2007_S7075677);
17
18
19    /// @name Analysis methods
20    /// @{
21
22    /// Book histograms
23    void init() {
24      ZFinder zfinder(FinalState(), Cuts::open(), PID::ELECTRON,
25                      71*GeV, 111*GeV, 0.2, ZFinder::ClusterPhotons::NODECAY, ZFinder::AddPhotons::YES);
26      declare(zfinder, "ZFinder");
27
28      book(_h_yZ ,1, 1, 1);
29    }
30
31
32    /// Do the analysis
33    void analyze(const Event & e) {
34      const ZFinder& zfinder = apply<ZFinder>(e, "ZFinder");
35      if (zfinder.bosons().size() == 1) {
36        const Particles& el(zfinder.constituents());
37        if (el[0].pT() > 25*GeV || el[1].pT() > 25*GeV) {
38          _h_yZ->fill(fabs(zfinder.bosons()[0].rapidity()));
39        }
40      } else {
41        MSG_DEBUG("No unique lepton pair found.");
42      }
43    }
44
45
46    // Finalize
47    void finalize() {
48      // Data seems to have been normalized for the avg of the two sides
49      // (+ve & -ve rapidity) rather than the sum, hence the 0.5:
50      normalize(_h_yZ, 0.5);
51    }
52
53    /// @}
54
55
56  private:
57
58    /// @name Histograms
59    /// @{
60    Histo1DPtr _h_yZ;
61    /// @}
62
63  };
64
65
66
67  RIVET_DECLARE_ALIASED_PLUGIN(D0_2007_S7075677, D0_2007_I744624);
68
69}