pith. sign in

arxiv: 2606.18517 · v1 · pith:NYUXC65Nnew · submitted 2026-06-16 · ✦ hep-ex

Electromagnetic Shower Reconstruction and Identification in FASER's Emulsion Detector for LHC Forward Neutrino Measurements

FASER Collaboration: Roshan Mammen Abraham , Xiaocong Ai , Saul Alonso Monsalve , John Anders , Emma Kate Anderson , Akitaka Ariga , Tomoko Ariga , Jeremy Atkinson
show 126 more authors
Florian U. Bernlochner Jianming Bian Tobias Boeckh Eliot Bornand Jamie Boyd Lydia Brenner Angela Burger Franck Cadoux Roberto Cardella David W. Casper Charlotte Cavanagh Shiyang Chen Xin Chen Xing Cheng Dhruv Chouhan Andrea Coccaro Fabio Cufino Stephane D\'ebieux Ansh Desai Sergey Dmitrievsky Radu Dobre Monica D'Onofrio Sinead Eley Yannick Favre Jonathan L. Feng Carlo Alberto Fenoglio Didier Ferrere Max Fieg Wissal Filali Elena Firu Haruhi Fujimori Edward Galantay Stephen Gibson Sergio Gonzalez-Sevilla Yuri Gornushkin Yotam Granov Jinjing Gu Carl Gwilliam Elie Hammou Daiki Hayakawa Michael Holzbock Shih-Chieh Hsu Zhen Hu Giuseppe Iacobucci Tomohiro Inada Luca Iodice Sune Jakobsen Cesar Jesus-Valls Arash Jofrehei Hans Joos Enrique Kajomovitz Alex Keyken Felix Kling Daniela K\"ock Pantelis Kontaxakis Jelle Koorn Umut Kose Peter Krack Susanne Kuehn Thanushan Kugathasan Sebastian Laudage Lorne Levinson Botao Li Jiaxi Liu Jinfeng Liu Yi Liu Margaret S. Lutz Joern Mahlstedt Toni M\"akel\"a Yasuhiro Maruya Anna Mascellani Lawson McCoy Josh McFayden Andrea Pizarro Medina Hiroaki Menjo Th\'eo Moretti Toshiyuki Nakano Laurie Nevay Yuma Ohara Ken Ohashi Hidetoshi Otono Lorenzo Paolozzi Annabelle Parry Pawan Pawan Brian Petersen Titi Preda Markus Prim Junkai Qin Michaela Queitsch-Maitland Juan Rojo Hiroki Rokujo Andr\'e Rubbia Osamu Sato Paola Scampoli Kristof Schmieden Matthias Schott Cristiano Sebastiani Anna Sfyrla Davide Sgalaberna Mansoora Shamim Yosuke Takubo Kakeru Tanaka Noshin Tarannum Simon Thor Eric Torrence Serhan Tufanli Oscar Ivan Valdes Martinez Svetlana Vasina Emanuele Villa Benedikt Vormwald Chi Wang Yuxiao Wang Eli Welch Aaron White Monika Wielers Benjamin James Wilson Zhongyi Wu Yue Xu Heng Yang Lekai Yao Daichi Yoshikawa Stefano Zambito Shunliang Zhang Yuxuan Zhang Xingyu Zhao Zijian Zhao
This is my paper

Pith reviewed 2026-06-26 21:32 UTC · model grok-4.3

classification ✦ hep-ex
keywords electromagnetic shower reconstructionemulsion detectorneutrino identificationFASERnubackground rejectionBDT classifiertest-beam validationenergy resolution
0
0 comments X

The pith

A multi-level identification chain reconstructs and identifies electromagnetic showers in the FASERnu emulsion detector with 99.99% background rejection at 100 GeV.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper develops methods to reconstruct electromagnetic showers from electrons in the FASERnu emulsion detector and to separate them from background. It uses 100 GeV and 200 GeV electron test-beam data collected at the CERN SPS H4 beamline to validate a clustering algorithm for shower axis finding that requires no energy-dependent tuning. A three-stage identification process of track pre-selection, cut-based selection, and a BDT classifier is shown to deliver combined background rejection of 99.99% at 100 GeV and 99.94% at 200 GeV while retaining total reconstruction-plus-identification efficiencies of 58.9% and 70.8%. Energy is estimated from the total number of reconstructed segments, producing relative biases of +0.6% and -0.8% with resolutions of 25.4% and 22.6%. These results supply a validated procedure for identifying electron neutrinos in the forward LHC neutrino program.

Core claim

The paper establishes that a clustering-based shower-axis algorithm without energy-dependent tuning, followed by a multi-level identification chain of track pre-selection, cut-based selection, and a BDT classifier, reaches background rejection rates of 99.99% (100 GeV) and 99.94% (200 GeV) together with total efficiencies of 58.9% (100 GeV) and 70.8% (200 GeV) when evaluated on simulated samples; the same chain yields energy reconstruction with relative biases of +0.6% and -0.8% and resolutions of 25.4% and 22.6% using the total number of reconstructed segments as the calorimetric estimator, with systematics dominated by emulsion film detection efficiency variations of roughly 10%.

What carries the argument

The multi-level identification chain of track pre-selection, cut-based selection, and BDT classifier, paired with a clustering-based shower-axis algorithm that requires no energy-dependent tuning.

If this is right

  • The validated chain supplies a framework for electron neutrino identification with the FASERnu detector at the LHC.
  • Energy reconstruction from total reconstructed segments carries relative biases below 1% and resolutions of 22-25%.
  • Systematic uncertainties on energy scale are dominated by variations in emulsion film detection efficiency, reaching approximately +10%/-8% at 100 GeV.
  • The reconstruction algorithm operates without energy-dependent tuning across the tested range.
  • Combined performance meets the requirements stated for forward neutrino measurements.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The same chain could be applied to measure the electron neutrino component of the forward flux and thereby test predictions of neutrino production in the LHC forward region.
  • Adaptation of the segment-counting energy estimator to other energies or detector configurations would require only re-calibration of the BDT rather than redesign of the axis finder.
  • Pairing the emulsion shower tag with timing or muon-veto information from the rest of the FASER apparatus could further reduce residual neutral-hadron background.
  • Running the identical selection on simulated samples at additional beam energies would test whether the quoted efficiencies and resolutions scale as expected.

Load-bearing premise

The simulated samples used to evaluate efficiencies and the test-beam data at the CERN SPS H4 beamline accurately represent the shower development, backgrounds, and detector response expected from electron neutrinos in the FASERnu detector during LHC operation.

What would settle it

Observation of background rejection or efficiency values in actual LHC neutrino data that lie well outside the 99.94-99.99% rejection and 58.9-70.8% efficiency ranges reported from the test-beam validation would falsify readiness of the method for physics use.

Figures

Figures reproduced from arXiv: 2606.18517 by Aaron White, Akitaka Ariga, Alex Keyken, Andrea Coccaro, Andrea Pizarro Medina, Andr\'e Rubbia, Angela Burger, Annabelle Parry, Anna Mascellani, Anna Sfyrla, Ansh Desai, Arash Jofrehei, Benedikt Vormwald, Benjamin James Wilson, Botao Li, Brian Petersen, Carl Gwilliam, Carlo Alberto Fenoglio, Cesar Jesus-Valls, Charlotte Cavanagh, Chi Wang, Cristiano Sebastiani, Daichi Yoshikawa, Daiki Hayakawa, Daniela K\"ock, Davide Sgalaberna, David W. Casper, Dhruv Chouhan, Didier Ferrere, Edward Galantay, Elena Firu, Elie Hammou, Eliot Bornand, Eli Welch, Emanuele Villa, Emma Kate Anderson, Enrique Kajomovitz, Eric Torrence, Fabio Cufino, FASER Collaboration: Roshan Mammen Abraham, Felix Kling, Florian U. Bernlochner, Franck Cadoux, Giuseppe Iacobucci, Hans Joos, Haruhi Fujimori, Heng Yang, Hidetoshi Otono, Hiroaki Menjo, Hiroki Rokujo, Jamie Boyd, Jelle Koorn, Jeremy Atkinson, Jianming Bian, Jiaxi Liu, Jinfeng Liu, Jinjing Gu, Joern Mahlstedt, John Anders, Jonathan L. Feng, Josh McFayden, Juan Rojo, Junkai Qin, Kakeru Tanaka, Ken Ohashi, Kristof Schmieden, Laurie Nevay, Lawson McCoy, Lekai Yao, Lorenzo Paolozzi, Lorne Levinson, Luca Iodice, Lydia Brenner, Mansoora Shamim, Margaret S. Lutz, Markus Prim, Matthias Schott, Max Fieg, Michaela Queitsch-Maitland, Michael Holzbock, Monica D'Onofrio, Monika Wielers, Noshin Tarannum, Osamu Sato, Oscar Ivan Valdes Martinez, Pantelis Kontaxakis, Paola Scampoli, Pawan Pawan, Peter Krack, Radu Dobre, Roberto Cardella, Saul Alonso Monsalve, Sebastian Laudage, Sergey Dmitrievsky, Sergio Gonzalez-Sevilla, Serhan Tufanli, Shih-Chieh Hsu, Shiyang Chen, Shunliang Zhang, Simon Thor, Sinead Eley, Stefano Zambito, Stephane D\'ebieux, Stephen Gibson, Sune Jakobsen, Susanne Kuehn, Svetlana Vasina, Thanushan Kugathasan, Th\'eo Moretti, Titi Preda, Tobias Boeckh, Tomohiro Inada, Tomoko Ariga, Toni M\"akel\"a, Toshiyuki Nakano, Umut Kose, Wissal Filali, Xiaocong Ai, Xin Chen, Xing Cheng, Xingyu Zhao, Yannick Favre, Yasuhiro Maruya, Yi Liu, Yosuke Takubo, Yotam Granov, Yue Xu, Yuma Ohara, Yuri Gornushkin, Yuxiao Wang, Yuxuan Zhang, Zhen Hu, Zhongyi Wu, Zijian Zhao.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic view of the FASER [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Experimental setup at the H4 beamline of the CERN SPS. The emulsion module and electromagnetic [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Overview of the EM shower reconstruction and identification pipeline, showing the three main [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Example of DBSCAN clustering on a 200 GeV electron MC shower. The identified shower core cluster [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Complete reconstruction chain for a representative 200 GeV test-beam shower candidate in Z-X [PITH_FULL_IMAGE:figures/full_fig_p009_5.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. Representative reconstructed events from the 200 GeV test beam module: a muon track (left) and [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗
Figure 7
Figure 7. Figure 7: FIG. 7. Distributions of the ten BDT input variables for EM shower signal (blue, from 100 GeV electron [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. BDT classifier performance: ROC curve with AUC = 0.997 (left), BDT score distributions for [PITH_FULL_IMAGE:figures/full_fig_p013_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. Distributions of shower observables at 100 GeV (top row) and 200 GeV (bottom row): Shower Max [PITH_FULL_IMAGE:figures/full_fig_p014_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Averaged longitudinal shower profiles from MC (red) and test beam data (blue) for 100 GeV (left) [PITH_FULL_IMAGE:figures/full_fig_p015_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Energy calibration using the total number of reconstructed segments as the estimator. Left: [PITH_FULL_IMAGE:figures/full_fig_p016_11.png] view at source ↗
read the original abstract

We present methods for electromagnetic shower reconstruction and identification in the FASERnu emulsion detector using 100 GeV and 200 GeV electron test-beam data from the CERN SPS H4 beamline. The reconstruction employs a clustering-based algorithm without energy-dependent tuning to determine shower axes. A multi-level identification chain comprising track pre-selection, a cut-based selection, and a BDT classifier achieves combined background rejection rates of 99.99% (100 GeV) and 99.94% (200 GeV). The method reaches total reconstruction and identification efficiencies of 58.9% (100 GeV) and 70.8% (200 GeV) evaluated from simulated samples. Energy reconstruction using the total number of reconstructed segments as the calorimetric estimator yields relative biases of +0.6% (100 GeV) and -0.8% (200 GeV), with resolutions of 25.4% and 22.6%, respectively. Systematic uncertainties on the energy reconstruction are dominated by variations in emulsion film detection efficiency, contributing (+10.9%/-8.2%) at 100 GeV and (+10.3%/-6.9%) at 200 GeV. The methodology provides a validated framework for electron neutrino identification with the FASERnu detector at the LHC.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript presents a multi-level electromagnetic shower reconstruction and identification method for the FASERnu emulsion detector, developed and tested with 100 GeV and 200 GeV electron test-beam data at CERN SPS H4. It employs a clustering algorithm for shower axis determination without energy-dependent tuning, followed by track pre-selection, cut-based selection, and a BDT classifier. Background rejection reaches 99.99% (100 GeV) and 99.94% (200 GeV). Total reconstruction and identification efficiencies of 58.9% (100 GeV) and 70.8% (200 GeV) are reported from simulated samples. Energy reconstruction via total reconstructed segments shows relative biases of +0.6% and -0.8% with resolutions 25.4% and 22.6%, respectively. Systematic uncertainties are dominated by emulsion film detection efficiency variations. The work claims to provide a validated framework for electron neutrino identification in FASERnu at the LHC.

Significance. If the reported performance holds under LHC conditions, the method supplies a practical, high-rejection tool for isolating electron neutrino interactions in the forward region, directly supporting FASERnu's goals of measuring neutrino cross sections and searching for beyond-Standard-Model physics. The test-beam validation of the reconstruction chain is a concrete strength, and the parameter-free clustering approach plus calorimetric energy estimator offer reproducible elements that could be adopted by similar emulsion-based experiments.

major comments (2)
  1. [Abstract and performance evaluation sections] Abstract and performance evaluation sections: The headline efficiencies (58.9% at 100 GeV, 70.8% at 200 GeV) and background rejection rates (99.99%/99.94%) are evaluated exclusively on simulated samples, yet the manuscript provides no quantitative comparison of simulated shower development, segment multiplicity, or background composition against the test-beam data or any other independent benchmark. This assumption is load-bearing for the central claim that the method constitutes a validated framework for LHC neutrino measurements.
  2. [BDT classifier description (likely §4)] BDT classifier description (likely §4): The manuscript does not specify BDT training and validation procedures, feature list, hyperparameter choices, or data exclusion criteria. Without these, the robustness of the combined rejection rates cannot be assessed, and the quoted performance numbers remain difficult to reproduce or extrapolate.
minor comments (1)
  1. [Systematic uncertainties paragraph] The systematic uncertainty estimation for emulsion film detection efficiency (contributing +10.9%/-8.2% at 100 GeV) is stated to dominate but lacks a description of how the variations were sampled or propagated through the reconstruction chain.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript. We address the major comments point by point below, agreeing that additional details and comparisons are required to strengthen the presentation.

read point-by-point responses
  1. Referee: [Abstract and performance evaluation sections] Abstract and performance evaluation sections: The headline efficiencies (58.9% at 100 GeV, 70.8% at 200 GeV) and background rejection rates (99.99%/99.94%) are evaluated exclusively on simulated samples, yet the manuscript provides no quantitative comparison of simulated shower development, segment multiplicity, or background composition against the test-beam data or any other independent benchmark. This assumption is load-bearing for the central claim that the method constitutes a validated framework for LHC neutrino measurements.

    Authors: We thank the referee for identifying this key point. The test-beam data provide pure electron showers and were used to develop and optimize the clustering algorithm, pre-selection, cuts, and BDT. Efficiencies and rejection rates were evaluated on simulation to incorporate the mixed backgrounds expected in neutrino interactions at the LHC, which are absent from the test-beam sample. We agree that a direct quantitative comparison is needed to support the validation claim. In the revised manuscript we will add a dedicated subsection comparing data and simulation for shower axis resolution, segment multiplicity distributions, and background composition, thereby demonstrating simulation fidelity and strengthening the applicability to FASERnu. revision: yes

  2. Referee: [BDT classifier description (likely §4)] BDT classifier description (likely §4): The manuscript does not specify BDT training and validation procedures, feature list, hyperparameter choices, or data exclusion criteria. Without these, the robustness of the combined rejection rates cannot be assessed, and the quoted performance numbers remain difficult to reproduce or extrapolate.

    Authors: We acknowledge that the BDT description is incomplete for reproducibility. The revised manuscript will expand the relevant section to include the complete list of input features, the composition of the training sample and the train/validation split, the validation procedure (including any cross-validation), the hyperparameter optimization method and final values, and the criteria used for data exclusion during training. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper develops a clustering-based reconstruction algorithm, cut-based selection, and BDT classifier using 100/200 GeV test-beam electron data, then reports efficiencies and background rejection evaluated on independent simulated samples. No step reduces by construction to its own inputs: the segment-count energy estimator is a direct observable, not a fitted parameter renamed as a prediction; no load-bearing self-citations or uniqueness theorems are invoked; and the multi-level identification chain is described without self-definitional loops. Standard separation of development data from evaluation samples keeps the chain non-circular.

Axiom & Free-Parameter Ledger

2 free parameters · 1 axioms · 0 invented entities

The central claim rests on the transferability of test-beam validation to LHC conditions and the fidelity of simulations for efficiency evaluation; limited information available from abstract alone.

free parameters (2)
  • BDT classifier thresholds and training parameters
    The BDT and selection cuts are optimized to achieve the reported rejection and efficiency; values chosen on simulated data.
  • Clustering algorithm parameters
    Parameters for the clustering-based shower axis determination without energy-dependent tuning.
axioms (1)
  • domain assumption Test beam conditions at CERN SPS H4 with 100 and 200 GeV electrons accurately mimic electromagnetic showers from electron neutrinos in the FASERnu detector at the LHC.
    Used to validate the method for the physics application at the LHC.

pith-pipeline@v0.9.1-grok · 6394 in / 1426 out tokens · 39834 ms · 2026-06-26T21:32:21.681124+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

17 extracted references · 5 canonical work pages · 1 internal anchor

  1. [1]

    Located 480 m downstream of the ATLAS interaction point along the beam collision axis, FASER accesses a kinematic region not covered by the main LHC detectors

    INTRODUCTION The Forward Search ExpeRiment [1, 2] (FASER) is designed to search for light, weakly inter- acting particles and to study high-energy neutrino interactions at the LHC [3]. Located 480 m downstream of the ATLAS interaction point along the beam collision axis, FASER accesses a kinematic region not covered by the main LHC detectors. The FASERνem...

  2. [2]

    EXPERIMENTAL SETUP AND DATA PROCESSING 2.1. Test Beam Configuration and Data Quality Emulsion modules were exposed to electron and muon beams at CERN SPS H4 beamline in 2023, providing mono-energetic electron samples at 100 GeV and 200 GeV. The experimental setup is shown in Figure 2. FIG. 2. Experimental setup at the H4 beamline of the CERN SPS. The emul...

  3. [3]

    SHOWER RECONSTRUCTION AND IDENTIFICATION Figure 3 provides an overview of the complete analysis chain. The pipeline proceeds in three stages: pre-selection applies track-level quality cuts to reject obvious muon backgrounds; shower reconstruction uses adaptive DBSCAN (Density-Based Spatial Clustering of Applications with Noise) clustering [20] and Huber r...

  4. [4]

    PERFORMANCE VALIDATION Tables 3 and 4 summarize the step-by-step efficiency of the full reconstruction and identification chain, where each entry is the fraction of events passing that step relative to all events entering it. Background rejection rates are measured directly from test beam data: since the sample is>99% muons, the rejection rate is well def...

  5. [5]

    This observable captures the full shower development and provides more robust energy estimation than localized estimators such asN 7-core

    ENERGY RECONSTRUCTION The total number of segments across all plates in the reconstructed shower,N total, is used as the calorimetric energy estimator. This observable captures the full shower development and provides more robust energy estimation than localized estimators such asN 7-core. SinceN total depends directly on the single-film reconstruction ef...

  6. [6]

    The beam energy spreadat the SPS H4 beamline is ∆p/p≈1.4% [25], which propagates directly to a∼1.4% uncertainty on the reconstructed energy

    SYSTEMATIC UNCERTAINTIES The systematic uncertainties in this analysis are dominated by three sources: the beam energy spread, background contamination, and single-film efficiency fluctuations. The beam energy spreadat the SPS H4 beamline is ∆p/p≈1.4% [25], which propagates directly to a∼1.4% uncertainty on the reconstructed energy. The beam has a transve...

  7. [7]

    Detecting and studying high-energy collider neutrinos with FASER at the LHC,

    CONCLUSION We have presented and validated methods for electromagnetic shower reconstruction and identi- fication in the FASERνemulsion detector using 100 GeV and 200 GeV electron test-beam data from the CERN SPS H4 beamline. The reconstruction algorithm combines adaptive DBSCAN cluster- ing, sequential centroid validation with a 60µm deviation threshold,...

  8. [8]

    LHC machine,

    L. Evans and P. Bryant, “LHC machine,”JINST3(2008) S08001. [4]FASERCollaboration, “The FASER detector,”JINST19(2024) P05066,arXiv:2207.11427. [5]FASERCollaboration, “Technical proposal of FASERνneutrino detector,” Tech. Rep. CERN-LHCC-2019-017, CERN, 2019.arXiv:2001.03073. [6]OPERACollaboration, “The OPERA experiment in the CERN to Gran Sasso neutrino bea...

  9. [9]

    Nuclear emulsions,

    A. Ariga, T. Ariga, G. D. Lellis, A. Ereditato, and K. Niwa, “Nuclear emulsions,”Particle Physics Reference Library: Volume 2: Detectors for Particles and Radiation(2020) 383–438

  10. [10]

    Hyper track selector nuclear emulsion readout system aimed at scanning an area of one thousand square meters,

    M. Yoshimoto, T. Nakano, R. Komatani, and H. Kawahara, “Hyper track selector nuclear emulsion readout system aimed at scanning an area of one thousand square meters,”Prog. Theor. Exp. Phys. 2017(2017) 103H01. [9]FASERCollaboration, “Neutrino rate predictions for FASER,”Phys. Rev. D110(2024) 012009, arXiv:2402.13318

  11. [11]

    Review of particle physics,

    Particle Data Group, “Review of particle physics,”Prog. Theor. Exp. Phys.2022(2022) 083C01. [11]FASERCollaboration, “First neutrino interaction candidates at the LHC,”Phys. Rev. D104(2021) L091101,arXiv:2105.06197. [12]FASERCollaboration, “First direct observation of collider neutrinos with FASER at the LHC,” Phys. Rev. Lett.131(2023) 031801,arXiv:2303.14...

  12. [12]

    A density-based algorithm for discovering clusters in large spatial databases with noise,

    M. Ester, H.-P. Kriegel, J. Sander, and X. Xu, “A density-based algorithm for discovering clusters in large spatial databases with noise,” inProceedings of the 2nd International Conference on Knowledge Discovery and Data Mining (KDD’96), pp. 226–231. AAAI Press, 1996

  13. [13]

    A K-means clustering algorithm,

    J. A. Hartigan and M. A. Wong, “A K-means clustering algorithm,”J. R. Stat. Soc. C28(1979) no. 1, 100–108

  14. [14]

    Robust estimation of a location parameter,

    P. J. Huber, “Robust estimation of a location parameter,” inBreakthroughs in Statistics, pp. 492–518. Springer, 1992

  15. [15]

    Scikit-learn: Machine learning in python,

    F. Pedregosa, G. Varoquaux, A. Gramfort, V. Michel, B. Thirion, O. Grisel, M. Blondel, P. Prettenhofer, R. Weiss, V. Dubourg, J. Vanderplas, A. Passos, D. Cournapeau, M. Brucher, M. Perrot, and ´Edouard Duchesnay, “Scikit-learn: Machine learning in python,”Journal of Machine Learning Research12(2011) no. 85, 2825–2830.http://jmlr.org/papers/v12/pedregosa11a.html

  16. [16]

    XGBoost: A Scalable Tree Boosting System

    T. Chen and C. Guestrin, “XGBoost: A scalable tree boosting system,”arXiv(2016) , arXiv:1603.02754

  17. [17]

    H4 beam line

    “H4 beam line.”https://be-dep-ea.web.cern.ch/experimental-areas/north-area/h4