1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4
 5namespace Rivet {
 6
 7
 8  /// @brief e+e- -> hadrons
 9  class FENICE_1996_I426675 : public Analysis {
10  public:
11
12    /// Constructor
13    RIVET_DEFAULT_ANALYSIS_CTOR(FENICE_1996_I426675);
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(FinalState(), "FS");
23
24      // Book histograms
25      book(_c_hadrons, 1, 1, 1);
26      for (const string& en : _c_hadrons.binning().edges<0>()) {
27        const double end = sqrt(std::stod(en)*GeV);
28        if (isCompatibleWithSqrtS(end)) {
29          _s = en;
30          break;
31        }
32      }
33      if (_s.empty())
34        MSG_ERROR("Beam energy incompatible with analysis.");
35    }
36
37
38    /// Perform the per-event analysis
39    void analyze(const Event& event) {
40      const FinalState& fs = apply<FinalState>(event, "FS");
41
42      map<long,int> nCount;
43      int ntotal(0);
44      for (const Particle& p : fs.particles()) {
45	nCount[p.pid()] += 1;
46	++ntotal;
47      }
48      // mu+mu- + photons
49      if(nCount[-13]==1 and nCount[13]==1 &&
50	 ntotal==2+nCount[22])
51	vetoEvent;
52      _c_hadrons->fill(_s);
53    }
54
55
56    /// Normalise histograms etc., after the run
57    void finalize() {
58      scale(_c_hadrons,crossSection()/ sumOfWeights() /nanobarn);
59    }
60
61    /// @}
62
63
64    /// @name Histograms
65    /// @{
66    BinnedHistoPtr<string> _c_hadrons;
67    string _s;
68    /// @}
69
70
71  };
72
73
74  RIVET_DECLARE_PLUGIN(FENICE_1996_I426675);
75
76
77}