On geometrical interpretation of alignment phenomenon

The observed alignment of spots in the X-ray films in cosmic ray emulsion experiments is analyzed and interpreted in the framework of geometrical approach. It is shown that the high degree of alignment can appear partly due to the selection procedure of most energetic particles itself and the threshold on the energy deposition together with the transverse momentum conservation.


Introduction
The intricate phenomenon of the strong collinearity of the most energetic cores of γ-ray-hadron families, closely related to coplanar scattering of secondary particles in the interaction, has been observed long time ago in mountainbased [1,2,3] and stratospheric [4] x-ray-emulsion chamber experiments.So far there is no simple satisfactory explanation of these cosmic ray observations in spite of numerous attempts to find it (see, for instance, [3,5,6,7,8,9,10] and references therein).The collider experiments do not show any evidence for exotic types of interactions although the collision energy at the LHC is already higher than the estimated threshold √ s eff > ∼ 4 TeV after which the alignment or some its collider analog should be clearly seen.The ridge effect [11] observed by the CMS collaboration in pp-collisions at the LHC stimulated the search of any manifestations [12,13] of alignment phenomenon at the LHC conditions.However, it is very difficult to establish some kind of the connection between these striking phenomena since they are observed in the quite different intervals of the rapidity practically without overlapping and in the different reference frames.Moreover after the intensive investigations the ridge effect is found a natural explanation in the framework of known standard interactions: for instance, as a simple interplay between the elliptic and triangular flows [14].This strengthens a rather widespread point of view (questioned in [7,15]) that the alignment phenomenon is not more than a tail in a distribution caused fluctuations.The main purpose of the present paper is just to add more argumentations in favour of this opinion and a great significance of the selection procedure of the most energetic particles with the threshold on energy deposition in the framework of the geometrical modeling.In Sect. 2 we formulate the problem on the whole.Section 3 describes the results of numerical simula-tion made under conditions close to emulsion experiments in the framework of the geometrical approach.A summary can be found in Sect. 4.

Alignment phenomenon and kinematics
In order to be clear let us recall thatin the Pamir experiment [1,3] the families with the total energy of the γ-quanta larger than a certain threshold and at least one hadron present were selected and analyzed.The alignment becomes apparent considerably at E γ > 0.5 PeV (that corresponds to interaction energies √ s eff > ∼ 4 TeV).The families are produced, mostly, by a proton with energy > ∼ 10 4 TeV interacting at a height h of several hundred meters to several kilometers in the atmosphere above the chamber [1,3].The collision products are observed within a radial distance r max up to several centimeters in the emulsion where the spot separation r min is of the order of 1 mm.Let us consider the kinematics in detail, since the kinematical relations are important in comparison the colliding beams experiment results with one with the fixed target.It is convenient to parameterize the 4-momentum of each particle i under consideration with its transverse momentum p T i (with respect to the collision axis z), rapidity η i and azimuthal angle φ i in the center-of-mass system: In this case the transformation from the center-of-mass system to the laboratory frame amounts to the rapidity shift: ζ i = η 0 + η i , where η 0 , ζ i are the rapidities of the center-of-mass system and the particle i respectively in the laboratory reference frame.Neglecting further the interaction of particles in the atmosphere (this would yield an estimate of the maximum alignment), one can easily determine the particle positions in the (xy)-plane in the film: where v zi and v ri are the longitudinal and radial components of particle velocity respectively.Since the size of the observation area is about several centimeters, these distances must obey the following restriction: The condition (3) means that the spots do not merge with the center formed by the particles emitted along the collision axis (predominantly region of incident-hadron fragmentation).The distinguishability of spots in the film imposes yet another constraint on the distance between particles: It must exceed r min : Otherwise the particles are combined into a cluster until there remain only particles and/or particle clusters with mutual distances larger than r min .The coordinates of the new cluster formed by two particles are determined in just the same way as the center-of-mass coordinates of two bodies in classical mechanics: Among clusters that satisfy the restrictions (3), ( 4) and ( 6) one selects the 2, ..., 7 clusters/particles which are most energetic.After that one calculates the alignment λ N using the common definition introduced by A. Borisov [2]: and taking into account the central cluster, i.e.N − 1 = 2, ..., 7.Here φ ijk is the angle between the two vectors (r k − r j ) and (r k − r i ) (for the central spot r = 0).This parameter, which changes from −1/(N − 1) to 1, characterizes precisely the disposition of N points just along the straight line.The combinatorial normalization factor (N (N −1)(N −2)) takes into account the number of choices to select the three different points among N points.In the case of N = 3 the first point can be selected by the three ways (the vertices of a triangle), then the second point can be selected by the two ways only.It means that every angle of a triangle is taken twice into account.By way of example, λ 3 = −0.5 in the case of the symmetrical configuration of three points in a plane (the equilateral triangle).Such a normalization allows the parameter λ N to be equal to 1 independently of the number of points under consideration if they lie exactly along the straight line.The parameter λ N describes the degree of alignment better than the possible other parameters of asymmetry like the eccentricity or the thrust.For instance, λ 4 will be equal to 1 if all four points lie strictly on the same straight line, but it will be considerably less than 1 if these points form four vertices of a long rectangle.
The degree of alignment P N is defined as the fraction of the events for which λ N > 0.8 [3] among the total number of events in which the number of cores not less than N .

Geometrical simulation of alignment
An absence of any collider experiment evidences for exotic types of interactions encourages us to investigate more intently the influence of the selection procedure of most energetic particles itself and the threshold on the energy deposition in appearing of alignment phenomenon.These studies are based on the geometrical and kinematical considerations being not influenced by the specific dynamics and therefore having more firm conclusions.We try to understand how the selection procedure and its restrictions influence on the degree of alignment P N for chaotically located spots in the x-ray film.
Let us consider simply the position of points in the restricted area of the (xy)-plane as an event.Their coordinates are randomly generated in the square centered in an origin.The size of this square can be fixed that the radius of the inscribed circle is equal to 1 as a characteristic radial scale.After that all distances needed can be easily recalculated in the traditional units.In this geometrical approach we can calculate the alignment for any configuration of points in the (xy)-plane including an origin into account.The results for the degree of alignment P N are the following: The degree of alignment P 3 (for the two chaotically located points and the central point in an origin) can be easily estimated having in mind that the probability to find a point in some area is simply proportional to this area square.Let us draw a circle through the most remote point A 1 with the center in an origin as it is shown in Fig. 1.The location of point B on the chord A 1 A 2 is determined from the condition that the triangle OA 1 B is similar to the triangle OA 1 A 2 .The position of point B on the chord A 1 A 2 is determined the same way for the triangles OA 1 B and OA 1 A 2 .At small enough angle α the alignment λ for the three points OA 1 A 2 and OA 1 B is easily calculated: The expansion cos(x) 1 − x 2 /2 being used.At any location of point A 2 inside the triangles OA 1 B and OA 1 B the alignment for the three points OA 1 A 2 will be larger        than 1 − 4α 2 since the both acute angles of the triangle OA 1 A 2 will be less than α in this case.At moving point A 2 along the line A 2 O to an origin the alignment for the three points OA 1 A 2 reduces to since the one acute angle of the triangle OA 1 A 2 will be close to zero and the second acute angle will be close to 2α in this case.At any location of point A 2 inside the sector OA 2 A 2 the alignment for the three points OA 1 A 2 will be larger than this minimum value (11) Thus we can use this minimum value of alignment for all points from the area painted over in Fig. 1 to obtain the simple restriction from below on the degree of alignment P min 3 (λ 0 ).This minimum value is simply the ratio of the square of the area painted over in Fig. 1 to the square of a circle: At λ 0 = 0.8 one obtains P min 3 (λ 0 ) 0.16 as a minimum estimation that is close to the real value 0.2 of the numerical simulation.Thus, to have the high degree of alignment the second randomly generated point should be located in the restricted specific region with respect to the first point.In fact, this real specific region is wider than one used in our simple estimations.Besides the characteristic angle is evaluated from the condition to provide the high degree of alignment larger than a some given value λ 0 in the whole area.The relative square of this specific region determines just the degree of alignment for three points: P 3 = 0.2 at λ 3 = 0.8.One can suppose that each additional point decreases the degree of alignment by this geometrical suppression factor therefore in this geometrical picture one can expect: that is in a good agreement with the numerical calculations (9) at P 3 = 0.2.We also verify that the degree of alignment P 3 (λ 0 ) is indeed proportional to √ 1 − λ 0 in accordance with the geometrical argumentations above.It worth noting that these findings are independent of the size of the observation region.Now we examine how the additional restrictions, imitating the experimental one in some sense, influence on the degree of alignment.Naturally, the experimental conditions and restrictions cannot be reproduced fully nevertheless the root of things can be adequately caught even in our simple geometrical procedure.The particle positions are determined by their transverse momenta and rapidities in the (xy)-plane in the film.The radius vector r i directs along the particle transverse momentum and its value r i = |r i | increases with the growth of transverse momentum and decreases with the rapidity growth.An existence of total energy threshold means that the threshold on the radial distances should be also.In our generation procedure we imitate this energetic threshold in the form: where R has meaning as the average "energetic threshold" on one particle.Further the total transverse momentum of all particles should be zero.The influence of transverse momentum conservation we take into account in the form of missing transverse momentum: The smaller the value ∆ the better the transverse momentum of the selected particles (or the total radius-vector of generated points in our language) should be balanced and vice versa the larger the value ∆ the greater the disbalance of the transverse momentum (or the radius-vector) is permissible for the selected particles (points).The results of our geometrical modeling with the restrictions ( 15) and ( 16) are shown in Fig. 2 for the three points.These restrictions together increase considerably the degree of alignment for three points.There is a reasonable and wide interval of the values of the threshold R and the disbalance ∆ in which the 100% degree of alignment is observed.In accordance with Fig. 1 the maximum disbalance of the radius-vector if the second randomly generated point is located inside the sector OA 2 A 2 .At ∆ < ∆ 0 for all generated points the degree of alignment P 3 is equal to 1 and at ∆ > ∆ 0 this degree begins to decrease.The numerical simulation shown in Fig. 2 confirms the existence of a such crucial value of the disbalance in accordance with simple geometrical reasons.
Since the energetic threshold ( 15) has the azimuthal symmetry, while the region in which the points are generated is a square with other type of symmetry, then at the large enough value of the threshold R (close to the maximum possible value R max ) the degree of alignment P N can be sensitive to the region form.To test this sensitivity we calculate also the degree of alignment as a function of the disbalance for regions with other symmetry types.Now the position of points are also randomly generated inside the ellipse centered in an origin.The large half axis of an ellipse is fixed to be equal 1 again to do not complicate our generation procedure.This is not crucial at all since the relative size values are meaningful in the geometrical approach only.Figure 2(b,c,d) shows the results for the ellipse eccentricities e = 0, 0.2, 0.5 respectively.In the case of e = 0 (i.e. of a circle) the approximate formula (17) is in a good agreement even for R close to the maximum possible value (R max = 1 in this case) since its derivation assumes the azimuthal symmetry.In the case of a square (Fig. 2(a)) at R close to the maximum possible value (R max = 1.414 in this case) the points are mainly located in the opposite corners of a square (the back-toback configurations), and the 100% degree of alignment is observed for the larger interval of ∆ in comparison with the estimation (17).The degree of alignment begins to decrease when the condition ( 16) allows the point positions to be in all corners of a square and the configurations with the small λ parameter (not only back-to-back) is admissible.The similar behaviour takes place for the ellipse with the large enough eccentricity e = 0.5 (Fig. 2(d)).At the moderate R the behaviour of the degree of alignment is insensitive to the region form as it can be expected.
With the restrictions ( 15) and ( 16) the degree of alignment for the four points increases also noticeably by a factor of the order of 3, in comparison with the value P 4 = 0.04 in Eq. ( 9) without these restrictions, at the relatively small ∆ and the moderate R as it is shown in Fig. 3.These results are practically insensitive to the region form.At the threshold R close to the maximum possible value the unexpected behaviour is observed at first glance: very small alignment at the small ∆.In fact, this is the effect of the odd number of the points randomly generated.At the large enough R the disbalance for these points is of the order R, if they lie close to the same line passing through an origin.The small value of ∆ does not allow such configurations to be generated and the degree of alignment is close to zero.This effect takes place for all types of the region symmetry considered.At the moderate R the small disbalance is achievable for the configurations above with the large enough value of λ, because the different radial distances are possible and the small radial distances are not completely excluded by the restriction (15) in this case.
At R close to R max and at the large enough ∆, the nonregular behaviour with the large fluctuations takes place in the case of a square (Fig. 3(a)) and in the case of an ellipse with the large enough eccentricity (Fig. 3(d)).In the case of a circle (Fig. 3(b)) such fluctuations are absent since the type of circle symmetry is the same as for the restriction (15) unlike the two cases above.A square has singled out directions along it diagonals and for an ellipse this is its large axis.
The degree of alignment for the five points increases also noticeably by a factor of the order of 3, in comparison with the value P 5 = 0.008 in Eq. ( 9) without restrictions ( 15) and ( 16), at the relatively small ∆ and the moderate R as it is shown in Fig. 4.These results are again practically insensitive to the region form.At R close to R max the fluctuations with a large magnitude are observed in the case of a square (Fig. 5(a)) and in the case of an ellipse with the large enough eccentricity (Fig. 5(b)) and are not observed in the case of a circle (Fig. 4(b)).Such different behaviour reflects the existence of singled out directions in the case of a square and an ellipse.
Thus our modeling demonstrates clearly that the conservation law of the transverse momentum in the form of the disbalance of the total radius-vector allows one to select automatically more line up configurations.The energetic threshold in the form of the threshold on the radial distances allows one to use the quite reasonable values of this disbalance to obtain the high degree of alignment.Our modeling is also independent of the absolute size scale of the region, but is sensitive to its form and boundary.
It worth noting that the degree of alignment P N is not enough large to match the central values of experimental measurements for N ≥ 4. For comparison the Pamir measurements are P exp 4 = 0.67 ± 0.33 and P exp 5 = 0.33 ± 0.23 [6], while our "best" values are P 4 P 5 0.1 taking into account the magnitude of fluctuations for the five points.Nevertheless these values are in the limits of the two standard deviations from the experimentally measured values as a minimum.We seem to exhaust all possibilities of the geometrical approach to explain the alignment phenomenon.In principle, to increase the degree of alignment for the four and five points in our modeling one can probe to generate points, for instance, with the inhomogeneous density introducing a some singled out direction.However such possible improvement would be no special meaning except as the simple phenomenological fitting, but could violate a coherent logical structure.

Discussion and conclusions
Let us discuss physical consequences which can be extracted from our geometrical modeling of alignment.First of all one should note that our results (9) without any restrictions coincide exactly with the results [8,9] obtained for P N in the framework of the Monte Carlo generator PYTHIA [16] for minimum bias events.It looks unexpectedly at first glance only.In fact the azimuthal angle distribution of particles is isotropic in the PYTHIA generator as the same as the angular distribution of points in our procedure.It means that the angular distribution is more meaningful than the radial distribution in appearing of alignment.The effect of the conservation law of transverse momentum manifests itself [8,9] at the relatively small height of the primary interaction in comparison with the height [1,3] estimated by the Pamir collaboration.In this case particles from both hard jets (with rapidities near η i 0 in the center-of-mass system) hit on the observation region and the high degree of alignment is possible.The introduction of threshold, E thr Σ ∼ E lab /2, on the total energy of all (N c − 1) selected clusters (without taking into account the energy deposition in the centre around r = 0), allows one to select the events with hard jets only in a "natural" physical way.If the process hardness is close to maximum for the given energy √ s, the estimated degree [8,9] of alignment is already comparable with the experimentally observed one.Since the most energetic particles are selected they should satisfy the conservation law of transverse momentum with some missing transverse momentum.The larger the total particle energy the better the transverse momentum of the selected particles should be balanced.In our geometrical modeling the effect of the conservation law of transverse momentum imitates over the disbalance of the total radius-vector and manifests itself in the noticeable increasing of the degree of alignment, especially for the three points, for which 100%alignment is feasible.Our simple geometrical modeling matches the azimuthally isotropic decay of a some fireball moving with the large collective velocity along the z-axis.
For instance, it can be the products of the quark-gluon plasma hadronization from the narrow forward rapidity interval.In this case the radial distances are determined by the particle transverse momenta only and our procedure above has a good physical foundation.Thus our investigations demonstrates that the high degree of alignment can appear partly due to the selection procedure of most energetic particles itself and the threshold on the energy deposition together with the transverse momentum conservation.These studies are based on the geometrical and kinematical considerations being not influenced by the specific dynamics and can explain an absence of any collider experiment evidences for exotic types of interactions.However one should note that the degree of alignment P N in our modeling is not enough large to match the central values of experimental measurements for N ≥ 4 living some room for other explanations.

Fig. 1 .
Fig.1.The painted over region is the permissible location of the second randomly generated point A2 with respect to the first one A1 where the high degree of alignment is observed for the three points.The central point O in an origin is fixed.

Fig. 2 .
Fig. 2. The degree of alignment P3 for the three points as a function of the disbalance ∆ at the different values of the threshold R for a square (a) and for an ellipse with the different eccentricities e = 0 (b), e = 0.2 (c) and e = 0.5 (d).

Fig. 4 .
Fig. 4. The degree of alignment P5 for the five points as a function of the disbalance ∆ at the different values of the threshold R for a square (a) and for an ellipse with the different eccentricities e = 0 (b), e = 0.2 (c) and e = 0.5 (d).

Fig. 5 .
Fig. 5.The degree of alignment P5 for the five points as a function of the disbalance ∆ at the the threshold R close to to the maximum possible value for a square with Rmax = 1.414(a) and for an ellipse with the large enough e = 0.5 and with Rmax = 1 (b).