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