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  RIVET_DECLARE_PLUGIN(JADE_1984_I203145);
79
80
81}