1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/ZFinder.hh"
 5
 6namespace Rivet {
 7
 8
 9  /// @brief Measurement of Z(->muon muon) pT differential cross-section
10  ///
11  /// @author Flavia Dias
12  class D0_2010_S8671338 : public Analysis {
13  public:
14
15    RIVET_DEFAULT_ANALYSIS_CTOR(D0_2010_S8671338);
16
17
18    ///@name Analysis methods
19    //@{
20
21    /// Add projections and book histograms
22    void init() {
23      Cut cut = Cuts::abseta < 1.7 && Cuts::pT > 15*GeV;
24      ZFinder zfinder(FinalState(), cut, PID::MUON, 65*GeV, 115*GeV, 0.2, ZFinder::ClusterPhotons::NONE, ZFinder::AddPhotons::YES);
25      declare(zfinder, "ZFinder");
26
27      book(_h_Z_pT_normalised ,1, 1, 1);
28      book(_h_Z_pT_xs ,2, 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        double ZpT = zfinder.bosons()[0].pT()/GeV;
37        _h_Z_pT_normalised->fill(ZpT);
38        _h_Z_pT_xs->fill(ZpT);
39      }
40    }
41
42
43    /// Finalize
44    void finalize() {
45      normalize(_h_Z_pT_normalised);
46      scale(_h_Z_pT_xs, crossSection()/sumOfWeights());
47    }
48
49    //@}
50
51
52  private:
53
54    /// @name Histogram
55    Histo1DPtr _h_Z_pT_normalised;
56    Histo1DPtr _h_Z_pT_xs;
57
58  };
59
60
61
62  RIVET_DECLARE_ALIASED_PLUGIN(D0_2010_S8671338, D0_2010_I856972);
63
64}