Rivet analyses referenceDELPHI_1999_I499183Event Shapes at 133, 161, 172 and 183 GeVExperiment: DELPHI (LEP) Inspire ID: 499183 Status: VALIDATED Authors:
Beam energies: (66.5, 66.5); (80.5, 80.5); (86.0, 86.0); (91.5, 91.5) GeV Run details:
A wide range of event shapes in $e^+e^-$ collisions for centre-of-mass energies 133, 161, 172 and 183 GeV. Source code: DELPHI_1999_I499183.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/Beam.hh"
4#include "Rivet/Projections/FinalState.hh"
5#include "Rivet/Projections/FastJets.hh"
6#include "Rivet/Projections/Thrust.hh"
7#include "Rivet/Projections/Sphericity.hh"
8#include "Rivet/Projections/Hemispheres.hh"
9#include "Rivet/Projections/ParisiTensor.hh"
10
11namespace Rivet {
12
13
14 /// @brief event shapes at 133, 161 172 and 183 GeV
15 class DELPHI_1999_I499183 : public Analysis {
16 public:
17
18 /// Constructor
19 RIVET_DEFAULT_ANALYSIS_CTOR(DELPHI_1999_I499183);
20
21
22 /// @name Analysis methods
23 /// @{
24
25 /// Book histograms and initialise projections before the run
26 void init() {
27
28 // Initialise and register projections
29 declare(Beam(), "Beams");
30 const FinalState fs;
31 declare(fs, "FS");
32 const Thrust thrust(fs);
33 declare(thrust, "Thrust");
34 declare(Sphericity(fs), "Sphericity");
35 declare(ParisiTensor(fs), "Parisi");
36 declare(Hemispheres(thrust), "Hemispheres");
37
38 // Book histograms
39 size_t ie=0;
40 for (double eVal : allowedEnergies()) {
41
42 const string en = toString(int(eVal));
43 if (isCompatibleWithSqrtS(eVal)) _sqs = en;
44
45 size_t offset , offset2;
46 if (ie==3) { offset = 1; offset2 = 1; }
47 else if (ie==2) { offset = 0; offset2 = 3; }
48 else if (ie==1) { offset = 0; offset2 = 2; }
49 else { offset = 0; offset2 = 1; }
50
51 book(_h[en+"thrust"] , 13+offset, 1, offset2);
52 book(_h[en+"major"] , 15+offset, 1, offset2);
53 book(_h[en+"minor"] , 17+offset, 1, offset2);
54 book(_h[en+"oblateness"] , 19+offset, 1, offset2);
55 book(_h[en+"sphericity"] , 21+offset, 1, offset2);
56 book(_h[en+"planarity"] , 23+offset, 1, offset2);
57 book(_h[en+"aplanarity"] , 25+offset, 1, offset2);
58 book(_h[en+"heavy_jet_mass"] , 27+offset, 1, offset2);
59 book(_h[en+"light_jet_mass"] , 29+offset, 1, offset2);
60 book(_h[en+"diff_jet_mass"] , 31+offset, 1, offset2);
61 book(_h[en+"wide_broading"] , 33+offset, 1, offset2);
62 book(_h[en+"narrow_broading"], 35+offset, 1, offset2);
63 book(_h[en+"total_broading"] , 37+offset, 1, offset2);
64 book(_h[en+"diff_broading"] , 39+offset, 1, offset2);
65 book(_h[en+"CParam"] , 41+offset, 1, offset2);
66 book(_h[en+"DParam"] , 43+offset, 1, offset2);
67
68 if (ie < 3) {
69 book(_p["thrust"][ie], 1,1,1+ie);
70 book(_p["major"][ie], 2,1,1+ie);
71 book(_p["minor"][ie], 3,1,1+ie);
72 book(_p["obl"][ie], 4,1,1+ie);
73 book(_p["heavy"][ie], 5,1,1+ie);
74 book(_p["light"][ie], 6,1,1+ie);
75 book(_p["diff"][ie], 7,1,1+ie);
76 book(_p["bmax"][ie], 8,1,1+ie);
77 book(_p["bmin"][ie], 9,1,1+ie);
78 book(_p["bsum"][ie], 10,1,1+ie);
79 book(_p["bdiff"][ie], 11,1,1+ie);
80 book(_p["C"][ie], 12,1,1+ie);
81 }
82 ++ie;
83 }
84 if (_sqs == "" && !merging()) {
85 throw BeamError("Invalid beam energy for " + name() + "\n");
86 }
87 }
88
89 void fillMoment(array<BinnedProfilePtr<int>,3>& mom, double val) {
90 double tmp=val;
91 for (size_t ix=0; ix<3; ++ix) {
92 mom[ix]->fill(std::stoi(_sqs), tmp);
93 tmp *= val;
94 }
95 }
96
97 /// Perform the per-event analysis
98 void analyze(const Event& event) {
99
100 // Get beams and average beam momentum
101 const ParticlePair& beams = apply<Beam>(event, "Beams").beams();
102 const double meanBeamMom = 0.5*(beams.first.p3().mod() + beams.second.p3().mod());
103 MSG_DEBUG("Avg beam momentum = " << meanBeamMom);
104
105 const Thrust& thrust = apply<Thrust>(event, "Thrust");
106 // thrust related observables
107 _h[_sqs+"thrust"] ->fill(1.-thrust.thrust() );
108 _h[_sqs+"major"] ->fill(thrust.thrustMajor());
109 _h[_sqs+"minor"] ->fill(thrust.thrustMinor());
110 _h[_sqs+"oblateness"]->fill(thrust.oblateness() );
111 fillMoment(_p["thrust"],1.-thrust.thrust() );
112 fillMoment(_p["major"] ,thrust.thrustMajor());
113 fillMoment(_p["minor"] ,thrust.thrustMinor());
114 fillMoment(_p["obl"] ,thrust.oblateness() );
115 // sphericity related
116 const Sphericity& sphericity = apply<Sphericity>(event, "Sphericity");
117 _h[_sqs+"sphericity"]->fill(sphericity.sphericity());
118 _h[_sqs+"planarity"] ->fill(sphericity.planarity() );
119 _h[_sqs+"aplanarity"]->fill(sphericity.aplanarity());
120 // hemisphere related
121 const Hemispheres& hemi = apply<Hemispheres>(event, "Hemispheres");
122 // standard jet masses
123 _h[_sqs+"heavy_jet_mass"]->fill(hemi.scaledM2high());
124 _h[_sqs+"light_jet_mass"]->fill(hemi.scaledM2low() );
125 _h[_sqs+"diff_jet_mass"] ->fill(hemi.scaledM2diff());
126 fillMoment(_p["heavy"],hemi.scaledM2high());
127 fillMoment(_p["light"],hemi.scaledM2low() );
128 fillMoment(_p["diff"] ,hemi.scaledM2diff());
129 // jet broadening
130 _h[_sqs+"wide_broading"] ->fill(hemi.Bmax() );
131 _h[_sqs+"narrow_broading"]->fill(hemi.Bmin() );
132 _h[_sqs+"total_broading"] ->fill(hemi.Bsum() );
133 _h[_sqs+"diff_broading"] ->fill(hemi.Bdiff());
134 fillMoment(_p["bmax"],hemi.Bmax() );
135 fillMoment(_p["bmin"],hemi.Bmin() );
136 fillMoment(_p["bsum"],hemi.Bsum() );
137 fillMoment(_p["bdiff"],hemi.Bdiff());
138 MSG_DEBUG("Calculating Parisi params");
139 const ParisiTensor& parisi = apply<ParisiTensor>(event, "Parisi");
140 _h[_sqs+"CParam"]->fill(parisi.C());
141 _h[_sqs+"DParam"]->fill(parisi.D());
142 fillMoment(_p["C"],parisi.C());
143 }
144
145
146 /// Normalise histograms etc., after the run
147 void finalize() {
148 for (double eVal : allowedEnergies()) {
149
150 const string en = toString(int(eVal));
151
152 for (auto& item : _h) {
153 if (item.first.substr(0,3) != en) continue;
154 normalize(item.second);
155 }
156
157 }
158 }
159
160 /// @}
161
162
163 /// @name Histograms
164 /// @{
165 map<string,Histo1DPtr> _h;
166
167 map<string,array<BinnedProfilePtr<int>,3>> _p;
168
169 string _sqs = "";
170 /// @}
171
172
173 };
174
175
176 RIVET_DECLARE_PLUGIN(DELPHI_1999_I499183);
177
178
179}
|