1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/Beam.hh"
 4#include "Rivet/Projections/Thrust.hh"
 5#include "Rivet/Projections/Hemispheres.hh"
 6#include "Rivet/Projections/ChargedFinalState.hh"
 7
 8namespace Rivet {
 9
10
11  /// @brief jet masses for 14, 22, 34.8 and 43,5 GeV
12  class TASSO_1989_I279165 : public Analysis {
13  public:
14
15    /// Constructor
16    RIVET_DEFAULT_ANALYSIS_CTOR(TASSO_1989_I279165);
17
18
19    /// @name Analysis methods
20    /// @{
21
22    /// Book histograms and initialise projections before the run
23    void init() {
24
25      const ChargedFinalState cfs;
26      declare(cfs, "CFS");
27      // Thrust
28      const Thrust thrust(cfs);
29      declare(thrust, "Thrust");
30      declare(Hemispheres(thrust), "Hemispheres");
31
32      int offset = 0;
33      if (isCompatibleWithSqrtS(14.0*GeV))
34	offset=1;
35      else if(isCompatibleWithSqrtS(22.0*GeV))
36	offset=2;
37      else if(isCompatibleWithSqrtS(34.8*GeV))
38	offset=3;
39      else if(isCompatibleWithSqrtS(43.5*GeV))
40	offset=4;
41      else
42	MSG_ERROR("Beam energy " << sqrtS() << " not supported!");
43      // Book histograms
44      book(_h_diff , 1, 1, offset);
45      book(_h_heavy, 2, 1, offset);
46      book(_h_light, 3, 1, offset);
47
48    }
49
50
51    /// Perform the per-event analysis
52    void analyze(const Event& event) {
53      const ChargedFinalState& cfs = apply<ChargedFinalState>(event, "CFS");
54
55      if (cfs.particles().size() < 3 ) vetoEvent;
56
57      const Hemispheres& hemi = apply<Hemispheres>(event, "Hemispheres");
58      _h_heavy->fill(hemi.scaledM2high());
59      _h_light->fill(hemi.scaledM2low() );
60      _h_diff ->fill(hemi.scaledM2diff());
61    }
62
63
64    /// Normalise histograms etc., after the run
65    void finalize() {
66      normalize(_h_diff );
67      normalize(_h_heavy);
68      normalize(_h_light);
69    }
70    /// @}
71
72
73    /// @name Histograms
74    /// @{
75    Histo1DPtr _h_diff,_h_heavy,_h_light;
76    /// @}
77
78
79  };
80
81
82  RIVET_DECLARE_PLUGIN(TASSO_1989_I279165);
83
84
85}