Search for highly-ionizing particles in pp collisions at the LHC’s Run-1 using the prototype MoEDAL detector

A search for highly electrically charged objects (HECOs) and magnetic monopoles is presented using 2.2 fb-1\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$\hbox {fb}{^{-1}}$$\end{document} of p-p\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$p-p$$\end{document} collision data taken at a centre of mass energy (ECM\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$_{CM}$$\end{document}) of 8 TeV by the MoEDAL detector during LHC’s Run-1. The data were collected using MoEDAL’s prototype Nuclear Track Detectord array and the Trapping Detector array. The results are interpreted in terms of Drell–Yan pair production of stable HECO and monopole pairs with three spin hypotheses (0, 1/2 and 1). The search provides constraints on the direct production of magnetic monopoles carrying one to four Dirac magnetic charges and with mass limits ranging from 590 GeV/c2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{2}$$\end{document} to 1 TeV/c2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{2}$$\end{document}. Additionally, mass limits are placed on HECOs with charge in the range 10e to 180e, where e is the charge of an electron, for masses between 30 GeV/c2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{2}$$\end{document} and 1 TeV/c2\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{mathrsfs} \usepackage{upgreek} \setlength{\oddsidemargin}{-69pt} \begin{document}$$^{2}$$\end{document}.

according to the Bethe-Bloch formula [20], massive singly charged particles traversing matter can be highly ionizing at low velocity, β (the particle velocity expressed as a fraction of the speed of light, c).The physics program of the MoEDAL experiment described [4] a number of such massive Singly Electrically Charged Object (SECO) HIP scenarios.Further studies have examined MoEDAL's sensitivity to SUSY SECOs [21,22] and doubly charged [22] SUSY particles.
In 1931 Dirac formulated a consistent description of a magnetic monopole [23] within the framework of quantum physics.This monopole is associated with a line of singularity called a Dirac string.Dirac derived his Quantization Condition (DQC) in order that this string has no physical effect: µ 0 e n SI units of Ampere-metres (1) where e is the electric charge of the particle probe, h is Planck's constant divided by 2π, g D is the magnetic charge, µ 0 is the permeability of free space and n is an integer.
The DQC indicates that if magnetic charge exists then the electric charge is quantized in units of e = 2πh/(µ 0 g D ).The value of g D is approximately 68.5e.Dirac's theory did not constrain the mass or the spin of the monopole.Further, the Dirac quantization condition indicates a coupling strength much bigger than one: α m = µ 0 g 2 D /(4πhc) ≈ 34.Thus, perturbation theory cannot be applied and cross-section calculations based on perturbation theory are not physically valid, although useful as a benchmark.
In 1974 't Hooft [24] and Polyakov [25] discovered monopole solutions of the non-Abelian Georgi-Glashow model [26].This model has only one gauge symmetry, SO(3), with a three component Higgs field.The mass of the 't Hooft-Polyakov MM was predicted to be around 100 GeV/c 2 .However, MMs with such a low mass were ruled out by experiment.Subsequently, Georgi and Glashow combined their electroweak theory with a theoretical description of strong nuclear forces to form a Grand Unified Theory (GUT) [27] using the single non-Abelian gauge symmetry, SU (5).In this GUT theory the MM would have a mass of ∼10 15 GeV/c 2 which is far too heavy to be directly produced at any foreseeable terrestrial collider.
The SM has an SU (2) × U (1) group structure that does not allow a finite-energy monopole.However, Cho and co-workers have modified its structure to admit the possibility of an "electroweak" monopole [28,29] with a magnetic charge of 2g D .Based on this work, Cho, Kim and Yoon (CKY) [30] have more recently presented an adaptation of the SM -including a non-minimal coupling of its Higgs field to the square of its U (1) gauge coupling strength -that permits the possibility of a finite energy dyon [31].
In another extension of the SM there exists a topologically stable, finite energy magnetic monopole with a mass estimated to lie in range 900 GeV/c 2 to 3 TeV/c 2 [32,33].This extension retains the same gauge group as the SM but possesses an extended fermion and Higgs sector where right-handed neutrinos are non-sterile.
The question of whether it is possible to create generalizations of the CKY model that are consistent with the SM was considered by Ellis, Mavromatos and You (EMY) [38].EMY concluded that there is a possibility that an "electroweak" monopole, consistent with the current constraints on the SM, may exist and be detectable at the LHC.The existence of a MM is such a theoretically well predicated and revolutionary possibility that the search for a MM has been carried out as each new energy frontier is broached.
We consider here only those models that admit a magnetic charge quantized in units of Dirac charge, g D , or a multiple of the Dirac charge.As g D = 68.5e, a relativistic monopole with a single Dirac charge will ionize ∼4700 times more than a relativistic proton.It is thus a prime example of a HIP.
The first searches for MMs and/or HECOs at the LHC were performed by the ATLAS and MoEDAL Collaborations in 8 TeV p − p collisions [5,6,9].At this stage, the ATLAS monopole search was sensitive to singly magnetically charged (1g D ) monopoles, whereas the MoEDAL search was sensitive to single and multiply charged monopoles.ATLAS and MoEDAL continued the quest for HIPs at Run-2.
In the case of MMs, the ATLAS and MoEDAL searches were complementary, in the sense that ATLAS utilized the MMs highly ionizing signature [8,15] whereas, until now, the MoEDAL experiment only exploited the induction technique to directly detect the magnetic charge [10][11][12].Extensive accelerator searches for HIPs at the LHC have also been undertaken [5,7,8,15].The latest result from the LHC is from an ATLAS experiment search for HECOs and monopoles using data taken during LHC's Run-2 at a centre-of-mass energy of 13 TeV [16].
In this paper we report the first use of the prototype MoEDAL Nuclear Track Detector (NTD) System, which relies on an ionization signal to detect HIPs in conjunction with the prototype MoEDAL trapping detector system that utilizes a Superconducting Quantum Interference Device (SQUID) to detect the presence of trapped magnetic charge.The complete prototype detector is shown in Fig. 2. A total of 2.2 fb −1 of p − p collision data was obtained during LHC's Run-1 at intersection point IP8 on the LHC ring using this detector and analyzed for evidence of HECOs.The precision of the luminosity measurement at IP8 during Run-1 is estimated to be 1.16% [45].
A DY mechanism provides a simple model for HIP pair production.Monopole pair production and spin-0 and spin-1 HECO pair production cross sections are computed using the Feynman-like diagram shown in Fig. 1(top).In the case of spin-1/2 HECOs DY production can take place via virtual photon or Z exchange [46], as depicted in Fig. 1(bottom).In the case of Drell-Yan processes of magnetic monopole production, the coupling of the magnetic charge to the Z boson is usually assumed to be absent.In specific models of such monopoles this is proven explicitly [32,33], since any Z-flux that could exist in the monopole solution would be concentrated inside the monopole core.In the case of electrically charged dyons this issue is model dependent.The spin 1 magnetic monopole [34,35] can have a non-zero magnetic moment, characterized by the parameter κ [36].In this analysis the value κ=1 is used as it is the only one that respects unitarity [37].
It should be noted that the large monopole-photon coupling places such calculations in the non-perturbative regime.In the case of HECOs, which are characterized by large electrical charges, the Drell-Yan (DY) diagram shown in Fig. 1 should undergo appropriate resumma-tion, to account for potential non perturbative quantum corrections, see e.g.[47], [48].Such techniques are beyond the scope of this paper, and will be the topic of a future investigation.

THE MOEDAL DETECTOR
MoEDAL's detector technology is radically different from the general-purpose LHC experiments, ATLAS and CMS.The MoEDAL detector, deployed alongside LHCb's VELO (VErtex LOcator) detector at IP8, employs two unconventional passive detection methodologies tuned to the discovery of HIPs.The first of these is a plastic NTD stack array to detect the ionization trail of HIPs.The second is a detector system comprised of aluminium absorber elements.This detector system is called the MMT (Magnetic Monopole Trapper) since it was used to trap HIPs with magnetic charge, that slow down and stop within its sensitive volume, for further laboratory analysis.Both of these detector systems are passive, requiring neither a trigger or readout electronics.The MoEDAL detector is described in more detail below.
The MoEDAL detector is exemplified by its ability to retain a permanent record, and even capture new particles for further study.The NTDs provide a tried-andtested and cost effective method to accurately measure the track of a HIP and its effective charge.Importantly, the NTD response was directly calibrated using heavyion beams at the CERN SPS.The second detector system, the MMT, ensures that a small but significant fraction of the HIPs produced are slowed down, stopped and trapped for further study in the laboratory.There are no SM particles that can produce such distinct signatures.Thus, even the detection in MoEDAL of few HIP messengers of new physics would herald a discovery.

Energy Loss of HIPs in MoEDAL
In the MoEDAL detector HIPs lose energy by ionization.The energy loss by ionization in the MMT detector is computed using Bethe-Bloch formula.For NTDs, the relevant quantity is the Restricted Energy Loss (REL) [49].For β <10 −2 , the REL is equal to the particle's total energy loss in the medium.At larger velocities, REL is the fraction of the electronic energy loss leading to the formation of δ-rays with energies lower than a cut-off energy T cut .The REL can be computed from the Bethe-Bloch formula restricted to energy transfers T<T cut with T cut a constant characteristic of the medium.For Makfrofol, which is the MoEDAL NTD used for the analysis reported in this paper, T cut ≤ 350 eV.The RELs for MMs and for HECOs in Makrofol are shown in Fig. 3 and Fig. 4, respectively.

The MMT Detector
The prototype MMT detector deployed for LHC's Run-1 consisted of 198 aluminium rods weighing a total of 163 kg.These rods were housed in an enclosure placed just underneath the beampipe at the upstream end of LHCb's VELO detector as shown in Fig. 5.After exposure the MMT's aluminium volumes are sent to the ETH Zurich Laboratory for Natural Magnetism where they are passed through a SQUID magnetometer to scan for the presence of trapped magnetic charge.A monopole will stop in the MMT detector when its speed falls below β ≤ 10 −3 .It then binds due to the interaction between the  The anomalously large magnetic moment of an aluminium nucleus gives rise to a monopole-nucleus binding energy of 0.5 -2.5 MeV [50], comparable to the shell model splittings.In any case, it is reasonable to assume that the very strong magnetic field of the monopole will rearrange the nucleus, permitting it to bind strongly to the nucleus.As reported in Ref. [50] monopoles bound in such a way would be trapped indefinitely.It would require fields well in excess of several Tesla for the lifetime of the trapped monopole state to compromise its detection by the MoEDAL trapping detector.We note that the MoEDAL detector is only subject to fields lower than ∼10 mT.

Calibration of the MMT Detector
A magnetic monopole captured in an MMT volume is tagged and measured as a persistent current in the SQUID coil encircling the samples' transport axis that passes through the SQUID magnetometer.The calibration of the magnetometer response is achieved using two independent techniques.In brief, the magnetometer calibration was obtained using a convolution method applied to a dipole sample, and validated using long thin solenoids that simulate a monopole of well-known magnetic charge.For more details see Ref. [54].These calibration methods agree to within 10%, which is taken as the pole strength calibration uncertainty.The magnetometer response has been determined by measurement to be charge-symmetric and linear in a range of magnetic charge 0.3 -300 g D .

The Nuclear Track Detector System
The MoEDAL NTD is arranged in modules deployed around IP8 in the VELO cavern.A prototype NTD array of 125 × 25 cm × 25 cm stacks was installed for Run-1 as shown in Fig. 2. Each module comprises three layers of 1.5 mm thick CR39 ® polymer, three layers of Makrofol DE ® and three layers of Lexan ® 0.5 and 0.25 mm thick, respectively, inside Aluminium bags (Fig. 6).Currently the Lexan ® foils serve as protective layers and are not analyzed.
In this analysis only the Makrofol NTDs are utilized.This is due to roughly a factor ten higher detection threshold in Makrofol than CR39 which results in substantially less "visual noise" in the etched plastic large due to spallation products arising from beam backgrounds.Beam background particles are generated by the interaction of the LHC beam with LHC machine elements and the machine environment.Thus, the analysis of the CR39 NTDs is considerably more time intensive.HIPs produced via the DY mechanism in LHC collisions during Run-1 are sufficiently highly ionizing that they can easily be detected with the Makrofol NTDs, obviating the need to scan the CR39 in the first pass.In the event of the observation of a candidate event in the Makrofol all 6 NTD sheets in the stack would have been analyzed.

The Etching Procedure
In plastic track-etch detectors, the passage of a heavily ionizing particle can produce a permanent damage of polymeric bonds in a cylindrical region ("latent track") extending few tens of nanometers around the particle trajectory (Fig. 7).By subsequent chemical etching the latent track is "amplified" and can be made visible under FIG.6. NTD module composition an optical microscope.In the etching process, the bulk of the material is removed at a rate v B and at a higher rate v T along the latent track.The damage zone is revealed under an optical microscope as a pair of cone shaped etch-pits, one on each face of the NTD sheet.Etch-pits surface openings have a circular shape for normally incident particles, otherwise they are elliptical.A single well measured etch pit is called a "track" candidate.If another etch pit is measured on the NTD sheet that is consistent with being the twin then we have a confirmed track candidate.
A sketch of an etch-pit at different etching times is shown in Fig. 7 for a normally incident particle crossing the detector with a constant energy loss.Two etching conditions were applied (Table I).The first is the socalled "strong" etching condition, allowing faster etching and yielding larger etch-pits that were easier to detect under visual scanning.Strong etching was applied to the first, most upstream, Makrofol foil in each module.The second, "soft etching," condition results in a slower etching process.This allows the etching process to proceed in several steps in order to follow the formation of etchpits.Soft etching is applied to subsequent Makrofol foils in the stack, if a candidate track is found in the first layer.In Fig. 8 are shown microphotographs of relativistic Pb 82+ tracks in Makrofol foils etched in (left) "strong conditions"; (right) "soft conditions."

Calibration of the NTD Detector
The response of the NTD is measured by the etching rate ratio, also called the reduced etch-rate, p = v T /v B , as a function of the particle's REL.Heavy ion beams are used to determine the detector response over a large range of energy losses, as discussed in ref. [55].The Makrofol was calibrated with ion beams of 158 A GeV Pb 82+ and 13 A GeV Xe 54+ energy per nucleon, at CERN's SPS.The calibration set-up included a stack of Makrofol foils placed upstream and downstream of an  After etching, the size of the etch pits was measured with an automatic scanning system providing the cone base area, and the coordinates of the center of the etch pits.Etch pits diameters typically range from 10 µm to 100 µm, with a modal value in the range 30 µm to 40 µm.The base area distributions of incoming ions and of their fragments is shown in Fig. 9.The projectile fragments have the same velocity and approximately the same direction as the incident ions.From the base area spectrum, the charge corresponding to each nuclear fragment peak can be identified, and the corresponding REL determined.A detailed description of the calibration procedure can be found in [55].
For each identified peak the reduced etch rate p, the Z/β and eventually the restricted energy loss are computed (Fig. 9).Calibration data thus obtained are shown in Fig. 10.For Makrofol, the mimimum detectable relativistic charge is Z/β ∼ 50, for both soft and strong etching.The REL corresponding to this detector threshold is ∼ 2700 MeV g −1 cm 2 .

Etching and Scanning of MoEDAL NTD
After exposure in the LHC IP8 region, the MoEDAL NTD stacks were taken to the INFN etching and scanning laboratory in Bologna.A global module reference system is created by drilling three reference holes -2 mm diameter -on each detector module.This coordinate system provides an accuracy of 100 µm on the determination of the position of a particle track over the detector surface.The stacks are then unpacked, the detectors foils labelled and their thickness measured on a grid of points uniformly distributed over the foil surface.
For the search reported in this paper only Makrofol foils were analysed.In each exposed stack, the most upstream Makrofol layer was etched in 6 N KOH + 20% ethyl alcohol at 65 • C.After 6 hours etching, etch-pits as small as 10µm would be detected under 20 × magnification.An efficiency of ∼ 99% was estimated by scanning foils exposed to ions as described below.Each Makrofol layer examined was manually scanned.Every detected surface structure was further observed under higher magnification and classified either as material defects or particle's track.If a pair of etch-pits is detected on the front and back sides of the foil, it was observed at larger (100 -500×) magnification.From the etch-pit size on the "front" and "back" surfaces, and the bulk etching rate, the incidence angle on each surface is estimated.It takes around 2.5 hours to scan one side of a Makrofol layer when using the microscope at 32× total magnification.
A pair of collinear incident and exiting etch-pits, consistent with pointing to the IP, is defined as a potential candidate "track".In particular, if a candidate is found in the first layer of a module the downstream Makrofol foils would be etched in 6 N KOH + 20% ethyl alcohol at 50 • C. The vast majority of spallation products arising from beam backgrounds have a very limited range in the NTD sheet, typically tens of microns, and only give rise to a single pit when the NTD sheet is etched.
An accurate scan of the downstream Makrofol sheets would then be performed using an optical microscope with high magnification (100 -200 ×) in a square region of about 1 cm 2 around the candidate expected position.If collinear etch-pits are found in all three Makrofol sheets the CR39 would then be scanned for etch-pits collinear with those in the Makrofol layers.A HIP candidate track requires collinear etch-pits in all six NTD sheets in the stack that points to the IP.Additionally, the REL estimated from the etch-pits dimensions (surface area, etchpit length) has to be consistent with the HIP hypothesis.However, no candidate "track" candidate was found.

The Detection Threshold for Makrofol
For the HIP to be detected its REL must be greater than the detection threshold of the Makrofol.The detection threshold will vary with the etching conditions.It will also vary with the angle of incidence (δ) of the HIP on the NTD.The connection between the threshold and the maximum angle of incidence (δ M ax ) to the normal to the NTD that the HIP can make and still be detected, is expressed by the relationship: p = 1 cos(δ M ax ) , where p is reduced etch-rate described above.The lowest threshold is obtained for a HIP impinging normally to the NTD.The curve obtained from a empirically based parameterization of the relation between δ M ax and the REL is shown in Fig. 11.This parameterization is used in the determination of the acceptance for HIPs incident on the MoEDAL's NTD stacks.

Efficiency and False Positives in the NTD Detectors
As described above the signal for the passage of a HIP messenger of new physics through a MoEDAL NTD stack would be a string of etch-pits in the stack, where an etchpit pair is due to the ingress and egress of the HIP passing though an NTD sheet.No such signal has ever been seen in this search, or observed by any other HIP search employing NTDs [2].Indeed, no candidates were seen in the 125 stacks examined (corresponding to 7.8 m 2 ), where only the most upstream sheet of the NTD stacks FIG.11.The maximum angle to the normal of the NTD plane within which the HIP will be detected.
were examined.
The absence of false positives using the NTD technique was also a feature of the astroparticle physics experiments MACRO [60] and SLIM [61] , which deployed a surface area of 1263 m 2 and 427 m 2 , respectively.Neither experiment observed a single HIP candidate.It should be noted that the NTD technique employed by these experiments are essentially identical to those employed at colliders.
The lack of false positives in the NTD technique at colliders or in astroparticle physics experiment raises the question of the false negatives or detector efficiency, where a signal exists but is not seen.This can be evaluated using the heavy ion beams that are used to calibrate NTD detectors.In the absence of beam backgrounds, the detection, or scanning, efficiency for the etch-pits due to heavy-ion HIPs with ionizing power above the NTD threshold was measured to be in excess of 99% as described directly below.
In order to estimate the detection efficiency of NTDs for HIPs in the presence of beam backgrounds we utilized NTD calibration stacks exposed to a relativistic lead-ion beam as described above.The stacks were comprised of sheets of Makrofol NTDs exposed to the beam backgrounds (LHC-exposed sheets) in the VELO cavern at the LHC for a year of data taking, interleaved with unexposed Makrofol NTD sheets (LHC-unexposed sheets).Plastic from the same production batch was used in calibration and standard data taking.
The NTDs sheets comprising the calibration stacks were then etched in the same way as the standard NTD stacks deployed for data taking during Run-1.The individual sheets were scanned using the same manually controlled optical scanning microscope technology employed to examine all MoEDAL NTD stacks.
The relativistic lead-ion calibration beam particles penetrate the whole stack allowing the signal etch-pits seen in the LHC unexposed sheets -where the signal can clearly be observed with a 100% efficiency -to serve as a map.The identification of etch-pits in the LHCunexposed sheets is measured to be 100% by making independent comparison scans of the other LHC-unexposed sheets in the stack which, of course, have the identical etch-pit number and pattern.
Using the LHC-unexposed sheets in the stack as a oneto-one same-scale map for the hits in the adjacent LHCexposed sheets the scanning efficiency for LHC-exposed sheets can be measured.Such measurements indicate that the overall scanning efficiency for detection above threshold was in excess of 99%.This number was found by scanning the LHC-exposed sheets and then comparing the etch-pits found with the etch-pits identified in the adjacent LHC-unexposed stacks.Each sheet has exactly the same number and pattern of " signal" etch-pits, resulting from the calibration beam -since the beam passes through the complete stack.
The MoEDAL detector's acceptance is defined to be the fraction of the number of events in which at least one HIP of the DY produced pair was detected in MoEDAL in either the NTD detector or the MMT detector.The acceptance for DY production of HECOs and magnetic monopoles is described by an interplay of the geometrical disposition of MoEDAL NTD modules and MMT detectors, energy loss in the detectors, mass of the particle and the spin-dependent kinematics of the interaction products.In the case of the HECOs, MoEDAL's NTD system provides the only means of detection.
For a given HIP mass and charge, the pair-production model determines the kinematics and the overall trapping acceptance obtained.The uncertainty in the acceptance is dominated by uncertainties in the material description [9][10][11].This contribution is estimated by performing simulations with hypothetical material conservatively added and removed from the nominal geometry model.An example, showing the MoEDAL NTD acceptance curves for spin-1/2, spin-0, spin-1 HECOs with charge 125e -produced by a DY process via virtual photon exchange -is shown in Fig. 12.Note that HECOs can only be detected by the NTDs since HECOs trapped in the MMT detector do not have a magnetic charge and hence cannot be detected in MoEDAL's SQUID detector.These curves are determined by the kinematics of the produced particles convoluted with the VELO material immediately surrounding IP8 and the distribution of the MoEDAL NTD detectors.The prototype MoEDAL detector provides only a partial, non-uniform, coverage FIG.12. Acceptance for spin-1, spin-0 and spin-1/2 HECOs with charge 125e, produced via a DY process with virtual photon exchange only.The dashed lines represent 4th order polynomial fits to the data.FIG. 13.Acceptance for monopole pair production with magnetic charge 2gD.The dashed lines represent 4th order polynomial fits to the data.
of the available solid angle.
The acceptance curves for spin-0, spin-1/2 and spin-1 monopoles, found using the NTD and MMT detectors combined, are shown in Fig. 13.The acceptances shown in Fig. 12 and Fig. 13 refer to the prototype detector deployed for LHC's Run-1.The acceptance for the Full Run-2 detector is somewhat larger.The difference between the acceptances for HECOs with spins with spin-1/2, spin-0 and spin-1, is mainly due to how well the disposition of the detector elements in theta (θ) matches the corresponding theta distribution of the DY produced HECO's.To illustrate this point we have compared the kinematic quantities (momentum ( − → p ), theta (θ) and phi (φ) ) of the uncut signal to the corresponding distributions for the events that pass the selection criteria.The angles theta (θ) and phi (φ) are defined in Figure 14.

FIG. 14. The coordinate system used in the analysis
In order to make the comparison clearer, given the large disparity between the number of generated and selected events, each plot was normalized to the the highest amplitude bins as selected by the Freedman -Diaconis rule [62].As an example we plotted the kinematic variables for DY production of HECOs with mass of 100 GeV/c 2 and electric charge 50e are shown in Figure 15, Figure 16 and Figure 17.The corresponding plots for DY production via photon/Z 0 * exchange, which affects only the spin-1/2 case, are shown in Figure 18.

ANALYSIS RESULTS
The first Makrofol sheet of each of MoEDAL's 125 NTD stacks, exposed during LHC's Run-1 were etched and scanned, as described above, for evidence of the passage through the sheet of a highly ionizing object such as a HECO or a magnetic monopole.The total area of plastic analyzed was 7.8 m 2 .No candidate events were observed.In addition, no monopole candidates were observed to be trapped in the MMT detector.This is the first time that the data from the full MoEDAL prototype detector, deployed during Run-1, has been presented.
Both the MMTs and the NTDs can be used to detect magnetic monopoles.Consequently, both the NTD and MMT detectors are incorporated into the analysis together.The Run-1 geometric acceptance of MoEDAL's NTD and MMT detectors did not overlap.Thus, the procedure to combine the monopole signal detected in the NTD and MMT detectors is simple addition.However, only the NTDs can be utilized for the HECO analysis since we have no way to detect electrically charged particles trapped in the MMT detectors.
FIG. 15.The momentum distribution for HECOs of mass 100 GeV/c 2 produced via the DY process with virtual photon exchange, with electric charge 50e, for each spin assignment.In each case, 60K events were originally generated.The plots are normalized to the maximum amplitude.The blue histogram represents generated events and the underlying light brown histogram shows the distribution of selected events.FIG.16.The theta angle distribution for HECOs of mass 100 GeV/c 2 with electric charge 50e produced via the DY process vias virtual photon exchange, for each spin assignment.In each case, 60K events were originally generated.The plots are normalized to the maximum amplitude.The blue histogram represents generated events and the underlying light brown histogram shows the distribution of selected events.
The dominant source of systematic error in this analysis arises from the imperfect knowledge of the amount of material between the interaction point and the MoEDAL NTD modules, due to LHCb's VELO detector.The VELO vacuum vessel and the various elements of the VELO detector within LHCb's physics acceptance are simulated with great precision in the LHCb geometry.However, detailed technical drawings of other elements of VELO outside of the physics acceptance such as cables, in-situ electronics, cooling pipes, various flanges, a vacuum pump and a vacuum manifold, are not available.
Nominally, this intervening material is between 0.1 and 8.0 radiation lengths (X 0 ) in thickness and on average around 1.4 X 0 [56] thick.The main contribution to the systematic uncertainty in this analysis arises from the estimate of the material in the Geant4 geometry description.The uncertainty in the material map is modelled by two geometries which represent an excess and a deficit of material, using conservative estimates of uncertainties on material thicknesses and densities, compared to the best assessment of the material budget that is compatible with direct measurement and existing drawings.
FIG. 17.The phi angle distribution for HECOs of mass 100 GeV/c 2 with electric charge 50e, produced via the DY process with virtual photon exchange, for each spin assignment.In each case, 60K events were originally generated.The plots are normalized to the maximum amplitude.The blue histogram represents generated events and the underlying light brown histogram shows the distribution of selected events.FIG.18.The momentum, theta and phi angle distribution for spin-1/2 HECOs of mass 100 GeV/c 2 and electric charge 50e, produced via the DY process via virtual photon and Z 0 exchange.In each case, 60K events were originally generated.The plots are normalized to the maximum amplitude.The blue histogram represents generated events and the underlying light brown histogram shows the distribution of selected events.This systematic uncertainty in the material map gives rise to uncertainties in the DY acceptances.For singly charge monopoles (|g| = g D ) the resulting relative uncertainty is of the order of 10% [9].This uncertainty increases with electric and magnetic charge.For a doubly charged monopoles (|g| = 2g D ) it is of the order of 10 -20% for intermediate masses, around 1 TeV.
Other sources of systematic error are an uncertainty due to a conservative estimate of 1 cm uncertainty in the trapping detector position.Simulations show this error lies in the range 1-17% [9].Another source of systematics is the uncertainty in dE/dx as a function of β, resulting in a 1-10% relative uncertainty in the acceptance [9].
In the case of monopoles and HECOs a systematic error on the variable p, due to the NTD etching and calibration process is given in Fig. 10 (bottom).This error on p can give rise to an error on the threshold value for detection of the plastic as well as an error on the variation of efficiency with angle of the NTD.However, these uncertainties are negligible compared to the error on the material map discussed above.All of the above sources of systematic error were added in quadrature and included FIG. 19.95% C.L. upper limits to the cross-section using as a measure a DY model with single virtual photon exchange for spin of 0, 1/2 and 1 HECO production.Upper limits for the production of spin-1/2 HECOS via DY production including virtual photon and Z 0 exchange is also included The solid lines denote the DY cross-sections for each case considered.
in the final limit calculation.
We calculated the 95% C.L. upper limits to the cross-section using as a measure a DY model for HECO and magnetic monopole production assuming a β−independent monopole coupling and that the monopole can have a spin of 0, 1/2 and 1.The limit curves obtained are shown in Fig. 19 for HECOs.The change in trend of the limit curves as the HECO charge 95% C.L. upper limits to the cross-section using as a measure a DY model for spin of 0, 1/2 and 1 magnetic monopole production assuming a β−independent monopole coupling.The solid lines denote the DY cross-sections for each case considered.
goes above 50e is due to the NTD threshold for detection of highly relativistic HECOs that lies at roughly 50e.For monopoles the cross-section upper limits versus mass are given in Fig. 20 for spin 0, 1/2 and 1.The values of the corresponding 95% C.L. mass limits are listed in Table II and Table III, for HECOs and magnetic monopoles, respectively.In the case of spin-1/2 HECOs we have included exclusive DY production limits, both from combined photon/Z 0 exchange and from simply photon exchange.This allows our result to be be compared with the best published charge limit on HECO production prior to this, from the ATLAS Collaboration [16], that considered only DY production of HECOs via photon exchange.
Note that in the spin-1/2 case and at the mass scale explored in this analysis, the cross-section of DY production of HECOs via photon exchange is slightly larger than the DY HECO production cross-section via photon/Z 0 exchange, due to destructive interference effects.The convolution of this lower cross-section with the harder momentum spectrum and better acceptance of HECOs produced by DY production via photon/Z 0 exchange, results in similar spin-1/2 mass limits to those obtained from HECOs produced by DY production via photon exchange only, over most of the charge range.CONCLUSIONS Both MoEDAL's prototype NTD system and aluminium elements of the MoEDAL MMT detector, were exposed to 8 TeV LHC collisions during LHC's Run-1.At the end of Run-1 both detector systems were examined for the presence of magnetic monopoles and/or HECOs.The NTDs were etched and scanned to reveal evidence for the passage of a magnetic monopole or a HECO using semi-automatic and manual optical microscopes.In the case of the MMT a SQUID-based magnetometer was also utilized to search for the presence of trapped magnetic charge.This is the first time that search results utilizing the NTD detectors are presented.
In previous MoEDAL searches [9] only MoEDAL's MMT detectors were utilized.Consequently, the HIP search was limited to magnetic monopoles.In this search the use of the NTDs allows the highly ionizing signature of the HIP to be registered.This permits both magnetically charged and electrically charged HIPs (HECOs) to be detected.
No magnetic monopole candidates were found.Consequently, limits on the DY production of magnetic monopole pair with cross-section in the range of approximately 40 fb to 5 pb were set for magnetic charges up to 4g D and mass as high as 1.2 TeV/c 2 .These limits are not competitive with recent Run-2 collider limits [12,16]    despite the use of the NTDs as well as the MMT subdetectors.This is due to a combination of: the limited acceptance of MoEDAL's MMT and NTD Run-1 prototype detectors compared to Run-2; the smaller E CM and DY cross-section at Run-1; and, the smaller luminosity of Run-1 compared to Run-2.
No evidence was found for DY produced HECO pairs.Thus, limits were placed on the DY production of HECO pairs with cross-sections from around 30 fb to 70 pb, for electric charges in the range 15e to 175e and masses from 110 GeV/c 2 to 1020 GeV/c 2 .The limits on the DY production of HECOs are the strongest to date, in terms of charge reach, at any collider experiment.

FIG. 4 .
FIG. 4. Restricted Energy Loss in Makrofol for HECOs of different electric charge.The horizontal dashed line indicates the Makrofol detection threshold.

C 3 .
FIG.7.Illustration of the track-etch technique: a) latent track forming along the trajectory of a high ionizing particle impinging perpendicularly on the NTD surface ; b) development of conical pits during the etching process; c) etch-pits joining after a prolonged etching, forming a hole in the detector.

FIG. 8 .
FIG. 8. Microphotographs of relativistic Pb 82+ tracks and of nuclear fragments (Z<82) in two consecutive foils of Makrofol.Each image frame measures 0.64 mm x 0.80 mm.Etch-pits are from the same ions crossing the detector foils: (left) Makrofol foil etched in "strong conditions"; (right) Makrofol foil etched in "soft conditions" .Note that the microphotographs also show two clearly differentiated fragmentation products of Pb: La (Z = 57); and, Pm (Z=61).

FIG. 10 .
FIG. 10.Reduced etch-rate (p) versus REL for Makrofol exposed to relativistic Lead and Xenon ion beams: (top) detectors etched in soft conditions;(bottom) detectors etched in strong conditions.The upper and low curves are drawn through the ±1σ value of the error on each p value, where the error bars represent a convolution of the statistical and systematic errors on each point.

FIG. 20 .
FIG. 20.95% C.L. upper limits to the cross-section using as a measure a DY model for spin of 0, 1/2 and 1 magnetic monopole production assuming a β−independent monopole coupling.The solid lines denote the DY cross-sections for each case considered.

TABLE I .
Etching Conditions of Makrofol

TABLE II .
95% CL mass limits for the HECO search.

TABLE III .
95% CL mass limits for the magnetic monopole search.