rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

LHCF_2012_I1115479

Measurement of forward neutral pion transverse momentum spectra for $\sqrt{s}$ = 7 TeV proton-proton collisions at LHC
Experiment: LHCF (LHC)
Inspire ID: 1115479
Status: VALIDATED
Authors:
  • Sercan Sen
References:
  • Phys.Rev. D86 (2012) 092001
  • arXiv: 1205.4578
Beams: p+ p+
Beam energies: (3500.0, 3500.0) GeV
Run details:
  • Inelastic events (ND+SD+DD) at $\sqrt{s}$ = 7 TeV.

The inclusive production rate of neutral pions has been measured by LHCf experiment during $\sqrt{s}=7$ TeV pp collision operation in early 2010. In order to ensure good event reconstruction efficiency, the range of the $\pi^0$ rapidity and $p_\perp$ are limited to $8.9 < y < 11.0$ and $p_\perp < 0.6$ GeV, respectively.

Source code: LHCF_2012_I1115479.cc
 1// -*- C++ -*-
 2#include "Rivet/Analysis.hh"
 3#include "Rivet/Projections/UnstableParticles.hh"
 4
 5namespace Rivet {
 6
 7
 8  class LHCF_2012_I1115479 : public Analysis {
 9  public:
10
11    RIVET_DEFAULT_ANALYSIS_CTOR(LHCF_2012_I1115479);
12
13  public:
14
15    void init() {
16      declare(UnstableParticles(),"UFS");
17
18      book(_g, {8.9, 9., 9.2, 9.4, 9.6, 10., 11.});
19      for (auto& b : _g->bins()) {
20        book(b, b.index(), 1, 1);
21      }
22    }
23
24
25    void analyze(const Event& event) {
26      const UnstableParticles& ufs = apply<UnstableParticles>(event, "UFS");
27
28      const double dphi = TWOPI;
29
30      for (const Particle& p : ufs.particles()) {
31        if (p.pid() == 111) {
32          double pT = p.pT();
33          double y  = p.rapidity();
34
35          if (pT > 0.6*GeV) continue;
36
37          const double scaled_weight = 1.0/(dphi*pT/GeV);
38          _g->fill(y, pT/GeV, scaled_weight);
39        }
40      }
41    }
42
43
44    void finalize() {
45      scale(_g, 1./sumOfWeights());
46      divByGroupWidth(_g);
47    }
48
49  private:
50
51    Histo1DGroupPtr _g;
52
53  };
54
55
56  RIVET_DECLARE_PLUGIN(LHCF_2012_I1115479);
57
58}