Rivet analyses referenceLHCF_2012_I1115479Measurement of forward neutral pion transverse momentum spectra for $\sqrt{s}$ = 7 TeV proton-proton collisions at LHCExperiment: LHCF (LHC) Inspire ID: 1115479 Status: VALIDATED Authors:
Beam energies: (3500.0, 3500.0) GeV Run details:
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}
|