rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

DELPHI_1999_I499183

Event Shapes at 133, 161, 172 and 183 GeV
Experiment: DELPHI (LEP)
Inspire ID: 499183
Status: VALIDATED
Authors:
  • Peter Richardson
References:
  • Phys.Lett. B456 (1999) 322-340, 1999
Beams: e- e+
Beam energies: (66.5, 66.5); (80.5, 80.5); (86.0, 86.0); (91.5, 91.5) GeV
Run details:
  • e+e- to hadrons.

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}