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 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 // -*- C++ -*- #include "Rivet/Analysis.hh" #include "Rivet/Projections/UnstableFinalState.hh" #include "Rivet/Tools/BinnedHistogram.hh" namespace Rivet { class LHCF_2012_I1115479 : public Analysis { public: LHCF_2012_I1115479() : Analysis("LHCF_2012_I1115479") { } public: void init() { declare(UnstableFinalState(),"UFS"); _binnedHistos_y_pT.addHistogram( 8.9, 9.0, bookHisto1D(1, 1, 1)); _binnedHistos_y_pT.addHistogram( 9.0, 9.2, bookHisto1D(2, 1, 1)); _binnedHistos_y_pT.addHistogram( 9.2, 9.4, bookHisto1D(3, 1, 1)); _binnedHistos_y_pT.addHistogram( 9.4, 9.6, bookHisto1D(4, 1, 1)); _binnedHistos_y_pT.addHistogram( 9.6, 10.0, bookHisto1D(5, 1, 1)); _binnedHistos_y_pT.addHistogram(10.0, 11.0, bookHisto1D(6, 1, 1)); } void analyze(const Event& event) { const UnstableFinalState& ufs = apply(event, "UFS"); const double weight = event.weight(); const double dphi = TWOPI; foreach (const Particle& p, ufs.particles()) { if (p.pid() == 111) { double pT = p.pT(); double y = p.rapidity(); if (pT > 0.6*GeV) continue; const double scaled_weight = weight/(dphi*pT/GeV); _binnedHistos_y_pT.fill(y, pT/GeV, scaled_weight); } } } void finalize() { _binnedHistos_y_pT.scale( 1./sumOfWeights() , this); } private: BinnedHistogram _binnedHistos_y_pT; }; // The hook for the plugin system DECLARE_RIVET_PLUGIN(LHCF_2012_I1115479); }