rivet is hosted by Hepforge, IPPP Durham

Rivet analyses reference

LHCF_2018_I1692008

Measurement of forward neutron production cross-section in proton-proton collisions at $\mathrm{\sqrt{s}=13\,TeV}$ with the LHCf detector
Experiment: LHCF (LHC)
Inspire ID: 1692008
Status: VALIDATED
Authors:
  • Eugenio Berti
  • LHCf collaboration
References: Beams: p+ p+
Beam energies: (6500.0, 6500.0) GeV
Run details:
  • Inclusive differential production cross section of forward neutrons (and antineutrons) produced in p-p collisions at $\mathrm{\sqrt{s} = 13\,TeV$}. Cross section is expressed as a function of energy in three different pseudorapidity regions, namely $\eta > 10.76$, $8.99 < \eta < 9.22$ and $8.81 < \eta < 8.99$. The measurements refer to all neutrons (and antineutrons) directly produced in the collisions or resulting from the decay of short life particles ($\mathrm{c \tau<1\,cm}$).

In this analysis, we report the measurement relative to the production of forward neutrons (and antineutrons) in proton-proton collisions at $\mathrm{\sqrt{s} = 13\,TeV}$ obtained using the LHCf Arm2 detector at the Large Hadron Collider. The results for the inclusive differential production cross section are presented as a function of energy in three different pseudorapidity regions: $\eta > 10.76$, $8.99 < \eta < 9.22$ and $8.81 < \eta < 8.99$. The analysis was performed using a data set acquired in June 2015 that corresponds to an integrated luminosity of $\mathrm{0.194\,nb^{-1}}$. The measurements were compared with the predictions of several hadronic interaction models used to simulate air showers generated by Ultra High Energy Cosmic Rays. None of these generators showed good agreement with the data for all pseudorapidity intervals. For $\eta > 10.76$, no model is able to reproduce the observed peak structure at around $\mathrm{5\,TeV}$ and all models underestimate the total production cross section: among them, QGSJET II-04 shows the smallest deficit with respect to data for the whole energy range. For $8.99 < \eta < 9.22$ and $8.81 < \eta < 8.99$, the models having the best overall agreement with data are SIBYLL 2.3 and EPOS-LHC, respectively: in particular, in both regions SIBYLL 2.3 is able to reproduce the observed peak structure at around $\mathrm{1.5-2.5\,TeV}$.

Source code: LHCF_2018_I1692008.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
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
// -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/Beam.hh"
#include "Rivet/Tools/BinnedHistogram.hh"
#include "Rivet/Projections/UnstableParticles.hh"

namespace Rivet {

  /// @brief forward neutron production cross-section at 13 TeV
  class LHCF_2018_I1692008 : public Analysis {
  public:

    /// Constructor
    DEFAULT_RIVET_ANALYSIS_CTOR(LHCF_2018_I1692008);

    /// @name Analysis methods
    //@{

    /// Book histograms and initialise projections before the run
    void init() {

      // Initialise and register projections
      declare(Beam(), "Beam");
      declare(FinalState(), "FS");

      // Book histograms
      book(_h_n_en_eta1, 1, 1, 1);
      book(_h_n_en_eta2, 2, 1, 1);
      book(_h_n_en_eta3, 3, 1, 1);
    }

    /// Perform the per-event analysis
    void analyze(const Event& event) {

      const FinalState &fs = apply<FinalState> (event, "FS");
      Particles fs_particles = fs.particles();

      for (Particle& p : fs_particles) {

        // select photons above threshold
        if (p.abspid() != 2112) continue; //Select neutrons and antineutrons only
        if (p.E()/GeV < 500.)   continue; //Select only energy above 500 GeV
        //if(p.pT()/GeV == 0.0) continue; //Artificially remove QGSJet II-04 wrong events

        // Double analysis efficiency with a two-sided LHCf
        const double eta = abs(p.eta());
        const double energy = p.E()/GeV;

        if ( eta > 10.758267 ) {
          _h_n_en_eta1->fill( energy );
        } 
        else if (eta > 8.994669 && eta < 9.217812) {
          _h_n_en_eta2->fill( energy );
        } 
        else if (eta > 8.812347 && eta < 8.994669) {
          _h_n_en_eta3->fill( energy );
        }
      }

    }

    /// Normalise histograms etc., after the run
    void finalize() {
      //Scale considering the LHCf Arm2 side
      scale(_h_n_en_eta1, crossSection()/millibarn/sumOfWeights()/2.); // norm to cross section
      scale(_h_n_en_eta2, crossSection()/millibarn/sumOfWeights()/2.); // norm to cross section
      scale(_h_n_en_eta3, crossSection()/millibarn/sumOfWeights()/2.); // norm to cross section
    }
    //@}

  private:

    /// @name Histograms
    //@{
    Histo1DPtr _h_n_en_eta1;
    Histo1DPtr _h_n_en_eta2;
    Histo1DPtr _h_n_en_eta3;
    //@}
  };

  // The hook for the plugin system
  DECLARE_RIVET_PLUGIN(LHCF_2018_I1692008);

}