|
|
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 | // -*- C++ -*-
#include "Rivet/Analysis.hh"
#include "Rivet/Projections/FinalState.hh"
namespace Rivet {
/// @brief forward neutron production cross-section at 13 TeV
class LHCF_2018_I1692008 : public Analysis {
public:
/// Constructor
RIVET_DEFAULT_ANALYSIS_CTOR(LHCF_2018_I1692008);
/// @name Analysis methods
//@{
/// Book histograms and initialise projections before the run
void init() {
// Initialise and register projections
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) {
// Select photons above threshold
const FinalState &fs = apply<FinalState> (event, "FS");
Particles fs_particles = fs.particles(Cuts::abspid==PID::NEUTRON && Cuts::E>=500/GeV && Cuts::abseta>8.812347);
for (const Particle& p : fs_particles) {
// Double analysis efficiency with a two-sided LHCf --- NOTE: taken care of in finalize division by 2
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.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
RIVET_DECLARE_PLUGIN(LHCF_2018_I1692008);
}
|
|