1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/FinalState.hh"
 4#include "Rivet/Projections/FastJets.hh"
 5#include "Rivet/Projections/UnstableParticles.hh"
 6
 7namespace Rivet {
 8
 9
10  /// @brief D_s spectrum at 34.7
11  class TASSO_1984_I194774 : public Analysis {
12  public:
13
14    /// Constructor
15    RIVET_DEFAULT_ANALYSIS_CTOR(TASSO_1984_I194774);
16
17
18    /// @name Analysis methods
19    ///@{
20
21    /// Book histograms and initialise projections before the run
22    void init() {
23      declare(UnstableParticles(), "UFS");
24      book(_h_Ds, 1, 1, 1);
25    }
26
27
28    /// Perform the per-event analysis
29    void analyze(const Event& event) {
30      UnstableParticles ufs = apply<UnstableParticles>(event,"UFS");
31      for(const Particle & p : ufs.particles(Cuts::abspid==431)) {
32      	double xE = 2.*p.E()/sqrtS();
33      	Vector3 mom3 = p.p3();
34        const double energy = p.E();
35      	double modp = mom3.mod();
36      	double beta = modp/energy;
37        _h_Ds ->fill(xE,1./beta);
38      }
39    }
40
41
42    /// Normalise histograms etc., after the run
43    void finalize() {
44      scale( _h_Ds , sqr(sqrtS())*crossSection()/nanobarn/sumOfWeights());
45    }
46
47    ///@}
48
49
50    /// @name Histograms
51    ///@{
52    Histo1DPtr _h_Ds;
53    ///@}
54
55
56  };
57
58
59  RIVET_DECLARE_PLUGIN(TASSO_1984_I194774);
60
61}