Rivet analyses referenceTASSO_1989_I279165Jet masses in e+e− between 14 and 44 GeVExperiment: TASSO (Petra) Inspire ID: 279165 Status: VALIDATED Authors:
Beam energies: (7.0, 7.0); (11.0, 11.0); (17.4, 17.4); (21.8, 21.8) GeV Run details:
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*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}
|
