Rivet analyses referenceJADE_1984_I203145Spectra of $\rho^0$ and $K^{*\pm}$ in $e^+e^-$ collisions at 35 GeVExperiment: JADE (PETRA) Inspire ID: 203145 Status: VALIDATED Authors:
Beam energies: (17.5, 17.5) GeV Run details:
Measurement of the $\rho^0$ and $K^{*\pm}$ spectra in $e^+e^-$ collisions for centre-of-mass energy 35 GeV by the JADE experiment at Petra. Source code: JADE_1984_I203145.cc 1// -*- C++ -*-
2#include "Rivet/Analysis.hh"
3#include "Rivet/Projections/UnstableParticles.hh"
4#include "Rivet/Projections/Beam.hh"
5
6namespace Rivet {
7
8
9 /// @brief rho0 and K*+/- spectra at 35 GeV
10 class JADE_1984_I203145 : public Analysis {
11 public:
12
13 /// Constructor
14 RIVET_DEFAULT_ANALYSIS_CTOR(JADE_1984_I203145);
15
16
17 /// @name Analysis methods
18 //@{
19
20 /// Book histograms and initialise projections before the run
21 void init() {
22
23 // Initialise and register projections
24 declare(Beam(), "Beams");
25 declare(UnstableParticles(), "UFS");
26
27 // Book histograms
28 book(_h_rho , 2, 1, 1);
29 book(_h_kstar, 3, 1, 1);
30
31 }
32
33
34 /// Perform the per-event analysis
35 void analyze(const Event& event) {
36
37 // Get beams and average beam momentum
38 const ParticlePair& beams = apply<Beam>(event, "Beams").beams();
39 const double meanBeamMom = ( beams.first.p3().mod() +
40 beams.second.p3().mod() ) / 2.0;
41 MSG_DEBUG("Avg beam momentum = " << meanBeamMom);
42
43 for (const Particle& p : apply<UnstableParticles>(event, "UFS").
44 particles(Cuts::abspid==323 or Cuts::pid==113)) {
45 double xE = p.E()/meanBeamMom;
46 double modp = p.p3().mod();
47 double beta = modp/p.E();
48 if(p.pid()==113) {
49 _h_rho->fill(xE,1./beta);
50 }
51 else {
52 _h_kstar->fill(xE,1./beta);
53 }
54 }
55 }
56
57
58 /// Normalise histograms etc., after the run
59 void finalize() {
60
61 scale(_h_rho , sqr(sqrtS())*crossSection()/nanobarn/sumOfWeights());
62 scale(_h_kstar, sqr(sqrtS())*crossSection()/nanobarn/sumOfWeights());
63
64 }
65
66 //@}
67
68
69 /// @name Histograms
70 //@{
71 Histo1DPtr _h_rho, _h_kstar;
72 //@}
73
74
75 };
76
77
78 // The hook for the plugin system
79 RIVET_DECLARE_PLUGIN(JADE_1984_I203145);
80
81
82}
|