rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

CELLO_1990_I283026

Spectra for strange hadron production at $E_{\text{CMS}}=35$ GeV
Experiment: CELLO (Petra)
Inspire ID: 283026
Status: VALIDATED
Authors:
  • Peter Richardson
References:
  • Z.Phys. C46 (1990) 397-404, 1990
Beams: e+ e-
Beam energies: (17.5, 17.5) GeV
Run details:
  • e+e- to hadrons at 35 GeV

Spectra for $K^0$, $K^{*+}$ and $\Lambda$ production in $e^+e^-$ collisions at a centre-of-mass energy of 35 GeV measured by the CELLO experiment at Petra. The data in HEPData only contains the centroids of the bins, not the bin edges which have therefore been inferred.

Source code: CELLO_1990_I283026.cc
 1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/ChargedFinalState.hh"
 5#include "Rivet/Projections/UnstableParticles.hh"
 6#include "Rivet/Projections/Beam.hh"
 7
 8namespace Rivet {
 9
10
11  /// @brief Strange hadrons at 35GeV
12  class CELLO_1990_I283026 : public Analysis {
13  public:
14
15    /// Constructor
16    RIVET_DEFAULT_ANALYSIS_CTOR(CELLO_1990_I283026);
17
18
19    /// @name Analysis methods
20    //@{
21
22    /// Book histograms and initialise projections before the run
23    void init() {
24      declare(Beam(), "Beams");
25      declare(ChargedFinalState(), "FS");
26      declare(UnstableParticles(), "UFS");
27
28      // Book histograms
29      book(_h_K0    , 1, 1, 1);
30      book(_h_Kstar , 2, 1, 1);
31      book(_h_Lambda, 3, 1, 1);
32
33    }
34
35
36    /// Perform the per-event analysis
37    void analyze(const Event& event) {
38
39      // First, veto on leptonic events by requiring at least 4 charged FS particles
40      const FinalState& fs = apply<FinalState>(event, "FS");
41      const size_t numParticles = fs.particles().size();
42      
43      // Even if we only generate hadronic events, we still need a cut on numCharged >= 2.
44      if (numParticles < 2) {
45        MSG_DEBUG("Failed leptonic event cut");
46        vetoEvent;
47      }
48      MSG_DEBUG("Passed leptonic event cut");
49      
50      // Get beams and average beam momentum
51      const ParticlePair& beams = apply<Beam>(event, "Beams").beams();
52      const double meanBeamMom = ( beams.first.p3().mod() +
53                                   beams.second.p3().mod() ) / 2.0;
54      MSG_DEBUG("Avg beam momentum = " << meanBeamMom);
55      
56      // Final state of unstable particles to get particle spectra
57      const UnstableParticles& ufs = apply<UnstableParticles>(event, "UFS");
58      
59      for (const Particle& p : ufs.particles()) {
60        const int id = p.abspid();
61        if (id == PID::K0S || id == PID::K0L) {
62          _h_K0->fill(p.E()/meanBeamMom);
63        }
64	else if(abs(id)==323) {
65          _h_Kstar->fill(p.E()/meanBeamMom);
66	}
67	else if(abs(id)==3122) {
68          _h_Lambda->fill(p.E()/meanBeamMom);
69	}
70      }
71    }
72
73
74    /// Normalise histograms etc., after the run
75    void finalize() {
76      scale(_h_K0    , 1./sumOfWeights());
77      scale(_h_Kstar , 1./sumOfWeights());
78      scale(_h_Lambda, 1./sumOfWeights());
79    }
80
81    //@}
82
83
84    /// @name Histograms
85    //@{
86    Histo1DPtr _h_K0, _h_Kstar, _h_Lambda;
87    //@}
88
89
90  };
91
92
93  // The hook for the plugin system
94  RIVET_DECLARE_PLUGIN(CELLO_1990_I283026);
95
96
97}