Does superhigh-energy cosmic-ray coplanarity contradict LHC data?

One of phenomena observed in cosmic-ray experiments at E0≳1016\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$E_0 \gtrsim 10^{16}$$\end{document} eV (s≳5\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\sqrt{s} \gtrsim 5$$\end{document} TeV) is a so-called coplanarity of most energetic subcores of γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}-ray–hadron families. This effect requires appearance of some coplanar generation of most energetic particles in hadron-nucleus interactions. To improve research for this phenomenon, a few versions of the new phenomenological FANSY 2.0 model are designed. This model makes it possible to simulate interactions via both traditional and coplanar hadron generation modes. The model reproduces ALICE, ATLAS, CMS, TOTEM data on low-rapidity particle generation, LHCf’s high-xF\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$x_\text {F}$$\end{document} and high-η\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\eta $$\end{document} data on γ\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\gamma $$\end{document}-rays and neutron-like hadrons as well as low-energy data on generation of stable and resonance hadrons, light-quark, strange and charmed particles. Features of coplanar particle generation of hadrons at superhigh energies are discussed. A test experiment at the LHC is proposed.


Introduction
A number of models of high-energy hadron-nucleus interactions is concurrently applied in cosmic-ray experiments. However, none of the models can reproduce the complete set of experimental extensive air shower (EAS) data with some sufficient accuracy (see [1], e.g.).
Besides, some results obtained in mountain-based and stratospheric X-ray-emulsion chamber (XREC) experiments are not explained yet. The present paper considers one of these phenomena, namely, a strong azimuthal effect that manifests itself in the form of a tendency to a coplanarity of most energetic subcores of so-called γ -ray-hadron families, i.e., groups of high-energy (E n · 1 TeV) particles (γ -rays, e ± , hadrons) in cores of relatively young extensive air showers, i.e., showers which do not reach the maximum of their a e-mail: rauf_m@mail.ru development, and the number of interactions of most energetic particles in these showers is much less than the average one. These showers are initiated mainly by protons and light nuclei of the primary cosmic radiation (PCR) with energies E 0 10 16 eV ( √ s 5 TeV). The effect was found by the Pamir Collaboration [2][3][4][5][6] and observed later in other experiments [7][8][9][10][11].
There is a point of view that this effect is not more than a giant fluctuation. This opinion is mainly based, on the one hand, on the fact that collider experiments do not inform us on such interactions. On the other hand, it is caused by desultory and fragmentary awareness based, as a rule, on separated results presented by different experimental groups with different XREC types, namely, so-called "carbon", "lead", "iron" (C-, Pb-, Fe-XREC) with a lateral resolution of the order of 100 µm.
However, to analyze this problem, it is necessary to take into account the following circumstances.
First, LHC experiments cannot give non-inclusive data on high-x F multiple hadron generation.
Secondly, five independent sets of experimental data were accumulated, namely, data of Pamir Collaboration's C-and Pb-XRECs on γ -ray families with energies E γ ≥ 700 TeV and data of Mt.Canbala Collaboration's Fe-XRECs E γ ≥ 500 TeV . Besides, only two stratospheric events (the Strana and JF2af2) with E γ > 1 PeV have been detected at very high altitudes, and both of these events demonstrate very high coplanarity of most energetic particles (MEP). It should be noted that the features of the stratospheric events are much more sensitive to the first interactions of PCR particles (due to a much smaller thickness of the atmospheric layer above the XREC detectors in each of these experiments), whereas events at mountain levels are strongly contaminated by secondary interactions of cascade hadrons.
All these data are related to hadron-nucleus interaction energies E 0 10 16 eV. The parameter was proposed by the Pamir Collaboration [2] to calculate the degree of a so-called alignment of N point-like tracks on a target plane that have appeared due to falling of MEPs to the XREC. Here ϕ k i j is the angle between two vectors, connecting the k-th point with the i-th and j-th points. The summation is performed over all possible combinations. The parameter λ N is equal to unity for N points aligned along a straight line and λ N ≈ −1/(N − 1) in the isotropic cases. Table 1 shows types of experimental data sets and used XRECs; criteria used to select experimental coplanar (or aligned) events, i.e., values of λ N and energy E γ of γray families; number of aligned families, N copl exp , and total family number N tot ; experimental fraction of aligned families, F copl exp , observed in each experimental data set; corresponding references (Columns 1−6, respectively). Background fractions of aligned γ -ray families, F copl sim , simulated using the QGSM-based FANSY 1.0 model [14] are given in Column 7.
It is obvious that although the statistical security of the data of each of the ground-based experiments is not too high, each of the data sets demonstrates some excess over the corresponding simulated value. Both stratospheric events are particularly unique in terms of both energy and extreme alignment of MEP tracks.
This phenomenon under consideration cannot be studied in EAS experiments due to a poor coordinate and energy resolution of detectors aimed at measurements of the total lowenergy charged-particle number, but not the narrow highenergy structure of the central EAS core.
The phenomenon is related to energies which were unreachable with the use of previous generation of accelerators. However, the LHC makes it possible to study this problem.
The probability for the total set of these experimental results to be produced by EAS fluctuations is rather negligible less than ∼ 10 −10 [12][13][14].
The phenomenon was initially interpreted as a manifestation of coplanar generation of most energetic particles in the fragmentation range (x Lab = E/E 0 0.01−0.05) characterized by relatively large transverse momenta, p copl t > 1 GeV/c. Besides, it was shown with using a nuclear-emulsion detector that p copl t can rise up to ∼ n × 10 GeV/c as early as at E 0 ∼ 10 16 eV [15,16]. However, it should be emphasized that this result should be interpreted with caution because of some methodological uncertainties of this experiment.
As regards γ -ray families observed in ground-based experiments, a more accurate measurement or, at least, an estimation of the MEPs' average p copl t in aligned events was not carried out in reality at the time of realization of the above-mentioned experiments due to difficulties in simulation of XREC response.
However, one can say that the average ratio of the MEPs' transverse momentum components lying in the coplanarity plane, p copl t , to the p t components directed normally to the coplanarity plane, p norm t , is large, i.e. p copl t / p norm t ∼ 10 [4].
In general, the coplanarity phenomenon is [12] • associated with most energetic secondary particles; • not described by traditional QGS models and QCD approach; • characterized by a relatively large cross section that is comparable with traditional inelastic proton-air highenergy interaction cross section estimations, σ p−air prod .
Just a few theoretical ideas were proposed to explain this phenomenon, namely, • conservation of the angular momentum of a relativistic fast-rotating quark-gluon string (QGS) stretched between colliding hadrons [17]; • semihard double inelastic diffraction (SHDID) [18], which assumes the coplanarity to be a result of some QGS tension in the diffraction cluster between a semihardly scattered constituent quark and other spectator quarks of the projectile hadron; • projectile's diquark breaking [19]; • appearance of very-high-spin leading systems [20,21]; • the most extraordinary and speculative assumption of a connection between this phenomenon and the hypothesis of "crystal world" with latticed and anisotropic spatial dimensions [22].
At least, it cannot be excluded that this phenomenon is a signature of some coherent Cherenkov-like hadron radiation at ultrahigh-energies [23].
It should be noted that the first four hypotheses include large transverse momenta as an almost mandatory element while the "crystal world" idea implies a gradual reduction of space dimensionality from three to two dimensions at sufficiently high energies. It entails the p t localization in some plane and a reduction of p t components directed perpendicularly to this plane.
Unfortunately, all of these ideas are not quantitatively developed, especially, the last one.
It looks like that in this situation we can use both largeand normalp t approaches to study this phenomenon.
While considering the angular-momentum conservation approach with respect to generation of secondary particles, Fermi [24] many years ago has neglected it "…to simplify the mathematics". Hagedorn [25] noted later that "this whole question, though of practical importance, seems to be still not understood. …it seems most reasonable to disregard angular momentum at all". Unfortunately, this tendency dominates until now, although energies under consideration increased since Fermi and Hagedorn's time by many orders of magnitude. As a result, to analyze the observed coplanarity phenomenon, we can construct only quite arbitrary models without any theoretical basis.
The coplanar generation of most energetic particles could be phenomenologically reproduced in the framework of two concepts, namely, (a) initial generation of quark-gluon strings with the following multiparticle production; (b) initial generation of an energetic leading system (like diffraction clusters, e.g.) with a lower multiplicity and higher average energies of particles. In this work, the first approach is only considered.
It should be noted that the second way must be considered separately as the contribution of diffraction processes into phenomena observed in cosmic-ray experiments is more significant than their formal contribution into the inelastic cross section.
The energy range corresponding to the coplanar generation of MEPs is studied at the LHC that makes it possible to test this phenomenon and, maybe, to make this problem more clear.
It would be very interesting to search for a possible relation of this phenomenon with the "ridge" effect found by the CMS Collaboration [35]. The "ridge" effect shows a two-particle long-range azimuthal-angle correlation at φ ≈ 0 and Δη 4. At first glance, these two phenomena could be somehow related. However, collider and cosmic-ray experiments are carried out with completely different event-selection criteria (very high particle multiplicity in the central pseudorapidity range in the case of the "ridge" effect and a few high-x F MEPs in coplanar γ -ray families). As a result, at present one can neither exclude nor confirm a relationship between these effects.
To study this problem, the phenomenological model FANSY 1.0 [13,14] has been developed a dozen years ago (before the launch of the LHC). It described main characteristics of γ -ray families (for instance, their intensity and mean transversal size) and helped to understand general features of the coplanar experimental events. However, FANSY 1.0 cannot properly describe some observed features of γ -ray families, namely, their average transversal size [36]. Besides, FANSY 1.0 cannot properly describe such LHC data as, for example, LHCf results on high-x F γ -rays ( [38], e.g.) and neutron-like hadrons (n,n, neutral kaons, Λ 0 ,Λ 0 etc.) which are really important for cosmic ray experiments.
To improve this situation, a new phenomenological model FANSY 2.0 [37] is designed. The model makes it possible to simulate QGSM-based traditional and coplanar particle generation (QGSJ and QGSCPG versions, respectively) with a higher accuracy of description of experimental data as compared with FANSY 1.0. Both the FANSY 2.0 versions reproduce high-x F and high-η data on γ -rays and neutron-like hadrons accumulated by the LHCf experiment, as well as ALICE, ATLAS, CMS, TOTEM data on generation of stable particles and resonances, mesons and baryons, light-quarks, strange and charmed particles. In addition, FANSY 2.0 reproduces low-energy data on generation of stable particles and resonances, strange and charmed particles and so on.
One can summarize the following key findings.
• a coplanarity of most energetic subcores of γ -ray-hadron families is observed in XREC experiments; • the phenomena is not described by traditional QGS models and QCD approach as well as fluctuations of EAS development; • the effect is characterized by a relatively large cross section that is comparable with σ p−air prod ; • a few qualitative theoretical ideas were proposed to explain this phenomenon; • a phenomenological model could be only used to study main features of this phenomenon.
This paper is organized as follows. Section 2 presents a short discussion on simulation problems related to the primary concept of coplanarity origin and describes shortly the FANSY 2.0 QGSCPG version with coplanar particle generation (CPG). Section 3 compares some experimental and simulated results on interactions features. Problems related to the proposed LHC experiment are considered in Sect. 4. Final comments are given in Conclusion.
More details of comparison of experimental data with results of simulation with the basic FANSY 2.0 QGSJ version are given in Ref. [37].

Coplanarity concepts
To analyze the coplanarity problem, first of all, the hypothesis of conservation of the angular momentum of relativistic fast-rotating quark-gluon string (QGS) stretched between colliding hadrons [17] is applied to develop the FANSY 2.0 QGSCPG version which simulates coplanarization of MEPs generated in nondiffraction processes.
Simulations using FANSY 2.0 QGSCPG tentative versions shows the following fundamental problem.
The original concept exploited by FANSY 1.0 [13,14] qualitatively explains the observed coplanarity of longitudinal momenta of MEPs with assuming their high transverse momenta, p copl t , in a coplanarity plane. However, in this case a significant p t growth suppresses dσ/dy and dσ/dη cross sections of hadron generation at highest |y| and |η| values and creates robust peaks at 2 |η| 4 which are contrary to LHC data (Fig. 1). This was an unsolvable problem for all the considered QGSCPG versions based on the initial highp copl t coplanarity concept. However, the coplanarity is actually observed in cosmicray experiments. Is it possible to reconcile this result with the LHC data?
Simulation shows that some agreement of LHC data and idea of coplanar generation becomes real only using a new concept of coplanarity origin, namely, assuming some decrease of particles' transverse momentum components directed normally to the coplanarity plane. This decrease While introducing azimuth effects, the algorithm keeps transverse momenta of each of MEPs constant. In this paper, some turning of particle's p t with no change of its absolute value is exploited. As a result, corresponding p t distributions are unchanged. Figure 2 shows three versions of observed tracks of particles, generated in the same imaginary interaction and arriving on a target plane placed at some distance from the interaction point. Most energetic particles are shown with large black symbols. One can see cases of traditional QGSM-like interaction (left), primary-concept coplanar interaction with increased p copl t (middle), and new-concept coplanar interaction with traditional average p copl t magnitudes (right). The geometric scale is given in arbitrary units.

Coplanarity simulation
As it is mentioned in Sect. 2.1, the hypothesis of conservation of the angular momentum of a relativistic fast-rotating quarkgluon string (QGS) stretched between colliding hadrons [17] is applied to develop the QGSCPG version which simulates coplanarization of the main bulk of secondary particles.
All interaction characteristics (excluding azimuthal ones) simulated with FANSY 2.0 QGSJ and QGSCPG versions, are similar. With decreasing energy, both the QGSJ and QGSCPG versions gradually become identical at √ s 2 TeV.
Simulation of CPG processes is realized as follows. All characteristics of particles in the QGSCPG version are primarily simulated with the traditional way, i.e., applying the QGSJ option. After this, hadron longitudinal momenta are listed in decreasing order, i.e., p max > p 1 > p 2 > p 3 > p 4 . . . p n > − p max . The corresponding rapidity values of Fig. 2 Tracks of particles, generated in the same interaction, on a target plane in cases of traditional particle generation (left), coplanar particle generation within the primary concept of increased p t (middle) and coplanar particle generation within the new concept of traditional p t (right). Most energetic particles are shown with large filled circles  these hadrons are denoted by y 1 , y 2 , y 3 , y 4 . . . y n . In pp interactions, p max = s/4 − m 2 p . After this, the so-called "coplanarization" procedure is implemented for MEPs in the coplanarity range (i.e., |y| > y copl thr ) of the hadron rapidity distribution. Figure 3 shows a schematic drawing of the traditional and coplanarity ranges of dn/dy distributions.
While beginning from the most energetic particle, transversal momenta of high-rapidity particles with |y| > y copl thr are coplanarized with a probability w copl (s) = σ CPG inel /σ inel , if the summary energy of secondary particles, s e f f , is higher than a fixed value, s CPG thr (although real observable effects appear at significantly higher s e f f values).
The position of the coplanarity plane is determined by the longitudinal momenta of the interacting projectile hadrons and the summary transverse momentum of two most energetic hadrons with momenta p 1 and p n after the interaction.
An algorithm rotates the transverse-momentum vector of all the hadrons with |y| > y copl thr towards the coplanarity plane along the shortest path and stops this procedure at |y| < y copl thr . The azimuth-angle distribution of transversal momenta of most energetic particles with rapidity y relative to the coplanarity plane obeys a Gaussian distribution with a standard deviation of σ copl Δφ (y) = σ copl 0 ·(y 2 /|y|). As a result, σ copl Δφ (y) is minimum for large-rapidity hadrons and increases with decreasing rapidity, i.e., the coplanarity is suppressed with decreasing hadron rapidities.
The value of y copl thr varies from ∼ 2 to ∼ 5 for hadrons with p z > 0 even at a constant Δy copl due to fluctuations of s e f f and hadron rapidity values. Similar procedure is carried out for hadrons with p z < 0. If the "ridge" effect [35] is associated with the coplanar generation of MEPs, then y copl thr seems to be of ∼ 2 − 3. Figure 4 shows the energy dependence of cross section of the coplanarization process simulation, σ CPG inel (s), as well as inelastic cross section, σ inel (s), in pp interactions. Observable results are determined by values of s CPG thr , Δy copl , σ copl 0 , and interaction energy. One can see that σ copl 0 is characterized with an energy threshold of √ s ≈ 1.4 TeV and fast growth up to 0.89 σ inel at √ s 10 TeV.

LHC data and FANSY 2.0 QGSCPG results
Experimental LHC data are mainly related to the central low-rapidity range. However, the coplanarity phenomenon is associated with most energetic particles in the cosmic-rayinitiated EAS core at the initial development stage, i.e., close to first interaction of the primary particle.
More detailed comparison of LHC data and results of simulation of pp interactions with the FANSY 2.0 QGSJ in a wide energy range from E 0 = 158 GeV to √ s = 13 TeV is given in Ref. [37]. Figure 5 shows dn ch /dη distributions based on data of the CMS+TOTEM experiment [41] for events selected using different criteria. Figure 5a shows so-called "NSD-enhanced" dn ch /dη distribution derived for selected events with n ch > 0 in the ranges −6.5 < η < −5.3 and 5.3 < η < 6.5 in inelastic pp interactions at √ s = 8 TeV. Figure 5b shows "more-forward" dn ch /dη distribution for events selected with n ch > 0 in the ranges −6.5 < η < −5.3 or 5.3 < η < 6.5. Figure 5c shows "SD-enhanced" dn ch /dη data derived for selected events with n ch > 0 only in the ranges −6.5 < η < −5.3 or only 5.3 < η < 6.5 at √ s = 8 TeV. The left part of this symmetric distribution is cut off. Figure 5d shows "more-forward" dn ch /dη data derived with displaced interaction points for selected events with n ch > 0 in the ranges −7.0 < η < −6.0 or 3.7 < η < 4.8.

CMS+TOTEM data
Experimental statistical and systematical errors are shown by shorter and longer limits, respectively. Corresponding FANSY 2.0 QGSCPG results are also shown in Fig. 5. Experimental and simulated distributions are obviously consistent within the systematical errors.
One can see a good agreement between the FANSY 2.0 results and collider data.

Search for coplanarity at the LHC
The LHC energy range is appropriate to study the CPG process. However, to test the above-described model in LHC experiments, it is necessary to solve a number of difficulties caused by collider's specific kinematics of experiments and a high beam luminosity. First of all, this effect can be studied only by event-by-event way, i.e., at lower collider luminosity. Second, the effect relates mainly to a large-x F range.
The LHCf detector with characteristic sizes of ∼ 10 cm, aimed at the study of high x F and η values and placed at L = 140 m from interaction points, is capable of detecting neutral particles with relatively low p t . However, this experiment detects, as a rule, only one particle to a few tens of interactions that makes it impossible to study azimuthal correlations of particles.
First of all, one should emphasize that the coplanarity effects are most clear and visible immediately after the CPG interaction. The following cascade development rapidly erode coplanarity effects with increasing depth of observation in the atmosphere. As a result, it can be observed only at high altitudes (in mountain or stratospheric experiments). It makes no sense to look for this phenomenon at the sea level. At mountain levels, it is necessary to separate energy flows in γ -families, not single particles only, with using a so-called decascading algorithm [42,43] for the following analysis.  [26,27], DPMJET 3.04 ( ) [28], SIBYLL 2.1 ( ) [29][30][31], EPOS 1.99 (×) [32,33], PYTHIA 8.145 ( ) [34] models are shown as well  Second, the proton-initiated γ -ray families dominate among observable events at any, more or less, reliable PCR mass composition due to XREC's high-energy threshold of detection and, as a result, lower efficiency of nuclei in formation of families, which decreases quickly with increasing mass of primary nuclei. On the other hand, each nucleusinitiated γ -ray family is a result of superposition of separated nucleon-initiated cascades. If so, such families are characterized by a lower coplanarity in the framework of any CPG-type model.
Therefore, to analyze the coplanarity problem in the first approximation, it is advisable to consider exactly the protonproton interactions. Table 2 shows fractions, F(λ 4 > 0.8), of interactions with four aligned tracks of stable MEPs on an imaginary target plane calculated with FANSY 2.0's traditional (QGSJ) and coplanar generation (QGSCPG) versions at √ s = 7 TeV. Definition of λ 4 is determined by Eq. (1). One can see a significant difference in results.
The CASTOR experiment [44] seems to be promising. Its operation pseudorapidity range is 5.25 < η < 6.5, while some CPG effects could manifest themselves in this range (at least, in the framework of FANSY 2.0).
Let us simplistically estimate potentialities of the CAS-TOR taking into account that much more detailed simulation is undoubtedly required for serious and quantitative analysis of this problem. Figure 9 shows a simplified cross section scheme of CAS-TOR and an example of detection of one coplanar interaction. The detector consists of 16 segments and is divided in the middle by a vertical gap. Particles are assumed to be detected if their pseudorapidity values are in the range of 5.3 < η < 6.5 and they do not fall into the vertical gap.
Black circles in Fig. 9 show tracks of particles. The larger the circle size, the higher the energy of the particle. Highenergy particles show a tendency to some coplanarity. Lowenergy particles form a more or less azimuthally symmetric halo.
To analyze observed interactions, energy values "measured" in each of the segments, E i , are used. Events with total energy E i ≥ 1 TeV are only analyzed, in which the number of segments, N s , with measured energy E i > E i min = 100 GeV is equal to or greater than the minimum Fig. 9 CASTOR's simplified cross section scheme with an example of detection of one coplanar interaction. The larger the circle size, the higher the energy of the particle Fig. 10 Probability dω/dε copl distributions predicted by FANSY 2.0 QGSJ and QGSCPG versions for the CASTOR detector value, N s N s min = 2. The first number is assigned to a segment with the maximum release of energy, E max . Here and below, E 1 = E max , E copl = E 1 + E 9 ; E tr = E 5 + E 13 , i.e., E tr is the energy measured in 9th and 13th segments, perpendicular to the first segment. A simple parameter is applied, namely, ε copl = E copl /(E copl + E tr ), which characterizes the event coplanarity degree. The degree of event coplanarity is maximum at ε copl = 1. Figure 10 shows dω/dε copl probability distributions (normalized to unity) reproduced with FANSY 2.0 QGSJ and QGSCPG versions. Obviously, different FANSY 2.0 versions predict different probability magnitudes to observe coplanar interactions at ε copl → 1. Of course, it is necessary to pay our attention on fluctuations, which form some peak near ε copl → 1 even in the case of the traditional interactions simulated in the framework of the FANSY 2.0 QGSJ version.

Conclusion
The FANSY 2.0 model is developed to study superhighenergy cosmic-ray "forward physics" interactions in general and the coplanar particle generation (CPG) phenomenon in particular.
At present, FANSY 2.0 includes two versions, i.e., the traditional QGSM-based QGSJ version and unconventional QGSCPG one including a CPG mode.
Comparison of results of these versions with experimental data demonstrates that the CPG process simulated with FANSY 2.0 QGSCPG version does not contradict LHC data. However, to reconcile experimental and simulated data, it is necessary to move from the primary concept of growth of transverse momenta of MEPs in the coplanarity plane to a new concept of reduction of transverse momentum component directed normally to the coplanarity plane so that the average p t values remain traditional. This new concept could be most naturally related to the speculative "crystal word" hypothesis only.
Some increase of the average transverse momentum of MEPs, p copl t (y > y copl thr ) , is not excluded. A model assuming that coplanar generation process takes place in double-diffraction interactions only [18] requires a separate consideration.
It would be very desirable to test the above-described model in LHC experiments.

Data Availability Statement
This manuscript has no associated data or the data will not be deposited. [Authors' comment: FANSY 2.0 data will be released after the completion of the model design, that is, the development of meson-proton and hadron-nucleus simulation algorithms and publication of relevant results.] Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecomm ons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. Funded by SCOAP 3 .