rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

TASSO_1989_I279165

Jet masses in $e^+e^-$ between 14 and 44 GeV
Experiment: TASSO (Petra)
Inspire ID: 279165
Status: VALIDATED
Authors:
  • Peter Richardson
References:
  • Z.Phys. C45 (1989) 11
Beams: e- e+
Beam energies: (7.0, 7.0); (11.0, 11.0); (17.4, 17.4); (21.8, 21.8) GeV
Run details:
  • $e^+ e^- \to$ hadrons. Beam energy must be specified as analysis option "ENERGY" when rivet-merging samples.

Light, heavy jet masses and their difference at 14, 22, 34.8, and 43.5 measured by the TASSO experiment. Beam energy must be specified as analysis option "ENERGY" when rivet-merging samples.

Source code: TASSO_1989_I279165.cc
 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))
34	offset=1;
35      else if(isCompatibleWithSqrtS(22.0))
36	offset=2;
37      else if(isCompatibleWithSqrtS(34.8))
38	offset=3;
39      else if(isCompatibleWithSqrtS(43.5))
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  // The hook for the plugin system
83  RIVET_DECLARE_PLUGIN(TASSO_1989_I279165);
84
85
86}