Systematic KMTNet Planetary Anomaly Search. V. Complete Sample of 2018 Prime-Field

We complete the analysis of all 2018 prime-field microlensing planets identified by the KMTNet AnomalyFinder. Among the 10 previously unpublished events with clear planetary solutions, 8 are either unambiguously planetary or are very likely to be planetary in nature: OGLE-2018-BLG-1126, KMT-2018-BLG-2004, OGLE-2018-BLG-1647, OGLE-2018-BLG-1367, OGLE-2018-BLG-1544, OGLE-2018-BLG-0932, OGLE-2018-BLG-1212, and KMT-2018-BLG-2718. Combined with the 4 previously published new AnomalyFinder events and 12 previously published (or in preparation) planets that were discovered by eye, thismakes a total of 24 2018 prime-field planets discovered or recovered by AnomalyFinder. Together with a paper in preparation on 2018 sub-prime planets, this work lays the basis for the first statistical analysis of the planet mass-ratio function based on planets identified in KMTNet data. By systematically applying the heuristic analysis of Hwang et al. (2022) to each event, we identify the small modification in their formalism that is needed to unify the so-called close/wide and inner/outer degeneracies, as conjectured by


Introduction
From its inception, and even conception, the Korea Microlensing Telescope Network (KMTNet, Kim et al. 2016) had as its major aim the construction and analysis of a large-scale statistical sample of microlensing planets. Nevertheless, during its first five years of full operations (2016)(2017)(2018)(2019)(2020), the overwhelming focus was on the detection and analysis of individual events of high scientific interest. In part, this focus reflected the new possibilities opened by KMTNet's continuous wide field coverage from three continents. For example, KMTNet played a major or decisive role in the detections of all three of the planets with mass ratios q < 3 × 10 −5 that were known by 2020 (Gould et al. 2020;Yee et al. 2021;Zang et al. 2021a) During this period, substantial work was carried out that would ultimately lay the basis for large-scale statistical studies. This included the development of a tiered observing strategy covering 97 deg 2 of the Galactic bulge ( Figure 12 of Kim et al. 2018a), as well as robust methods of identifying of order 3000 microlensing events per year using the EventFinder and AlertFinder systems (Kim et al. 2018a,b).
However, a number of practical, technical, and scientific challenges impeded the inauguration of large-scale statistical studies. At the most basic level, the online photometry remained of mixed quality until 2019. This did not prevent high-precision analysis of individual events because, from the beginning, KMT-Net had a tender-loving care (TLC) system of data re-reduction based on pySIS (Albrow et al. 2009), which returned high-A&A proofs: manuscript no. ms quality photometry on an event-by-event basis. However, it did mean that automated planet-searches of the KMTNet database would have yielded difficult-to-interpret results. In 2019, a new end-of-season pipeline was put into place that produced goodquality photometry for the great majority of events. This enabled the first KMTNet statistical study, a search for free-floating planet (FFP) candidates in the 2019 database (Kim et al. 2021). The same pipeline was gradually applied to the three previous seasons, but this labor-intensive work was only completed in November 2021.
Making use of these improved databases, Zang et al. (2021b) developed a new AnomalyFinder algorithm that was adapted to the characteristics of KMTNet, i.e., combining unprecedented quantities of microlensing data from three sites operating under very different conditions. The key innovation was to fit for "anomalies" in the residuals rather than for planets in the original light curves, which permitted a reduction of the search from three to two dimensions and also vastly simplified the modeling. This dimensional reduction is adapted from the KMT-Net EventFinder algorithm (Gould 1996;Kim et al. 2018a), and like EventFinder, it results in many false positives for each true anomaly, which must then be rejected by human review. However, in contrast to EventFinder, which annually results in O(5×10 5 ) false positives on O(5×10 8 ) catalog stars, the Anoma-lyFinder yields O(1 × 10 4 ) false positives on O(3 × 10 3 ) microlensing events. That is, while the specific false-positive rate is larger by 3.5 orders of magnitude, the total number of false positives is smaller by a factor 50, making human review much more tractable. In particular, it is quite feasible for several people to independently conduct this review as a cross check.
The specific false-positive rate is larger because the search is much more aggressive, i.e., attempting to discover all planetary anomalies down to a very low threshold. In particular, for AnomalyFinder, the operator may be shown dozens of potential anomalies, whereas for EventFinder only the highest-χ 2 candidate event is shown. Stated otherwise, the search can be much more aggressive because the number of light curves has been reduced from 5 × 10 8 to 3 × 10 3 , i.e., by 10 5 .
Another practical obstacle was the large human effort required for TLC reductions, which often took of order one day of work for each event. Again, this is not a major problem if one is publishing of order a dozen events per year. However, a statistical analysis requires not only the accurate parameter characterization of all "interesting" planets, but of all planets, and more dauntingly, all anomalous (or potentially anomalous) events that might plausibly be planetary. We estimate that this will be of order 200 TLC reductions for 2016-2019. Motivated by this challenging situation, H. Yang et al. (2022, in preparation) developed a quasi-automated TLC system that reduces the average reduction time to about one hour.
Our immediate goal is to prepare a complete sample of AnomalyFinder events from 2018 that can be compared to planet detection efficiency calculator (Jung et al. 2022, in preparation). This will be the first step toward the analysis of the 2016-2019 sample.
All of the remaining 85 events were fitted using online data, i.e., pipeline photometry. Of these, four new planets have already been published, including three by Hwang et al. (2022) in a study of low-q planets, and one by Wang et al. (2022) as part of a study of wide-orbit planets. Of the remaining 81, 56 were found to have q > 0.06, one (OGLE-2018-BLG-2014 was determined to be in the range (0.05 < q < 0.06), and 24 required TLC reductions, either because they were potentially planetary, q online < 0.05, or because the light curve could not be reliably characterized without TLC reductions. Of these 24, the 7 that have planetary solutions are analyzed here. Note that the 28 events that required TLC (24 analyzed here, and four previously published planets), were distributed among the five classification categories (planet, planet/binary, binary/planet, binary, finite source) as (9,11,4,3,1) of which (8,1,0,0,1) ultimately proved to have unambiguous planetary solutions and (1,1,0,0,0) ultimately proved to have planetary solutions that were ambiguous. This shows that great majority of events that ultimately prove to have planetary solutions are classified at the first stage as "planets" and that the great majority of events so classified prove to be planetary. We also analyze 3 of the 4 such events that were listed as "in preparation" in Table 11 of (Hwang et al. 2022) (namely, OGLE-2018-BLG-0932, OGLE-2018-BLG-1554, and OGLE-2018-BLG-1647 ), for a total of 10 events with planetary solutions 1 . These 10 include 8 that are clearly or very likely planetary in nature (q < 0.03) and 2 others that have an ambiguous nature. One additional event (OGLE-2018-BLG-0856/KMT-2018-BLG-2392) was found to have a solution with 0.03 < q < 0.05 based on TLC data. We follow a policy of noting such events but excluding them from the sample.
In sum, based on previous analyses and the current work, the prime-field sample has a total 26 planets or possible planets, of which 23 have unambiguous mass-ratio determinations, making them potentially suitable for a statistical analysis. Note that these must still be vetted for various effects, for which we provide detailed guidance in the text. The 3 others are clearly unsuitable because they are subject to multiple interpretations in q.
Note that the fraction of events that were subjected to AnomalyFinder analysis that were initially classified as planet and/or binary (110/843 = 13%) and those finally determined to be planetary or possibly planetary (26/843 = 3.1%) are both very similar the corresponding ratios in the first highcadence 24/7 microlensing survey that was carried out by Shvartzvald et al. (2016). They found 26/244 = 13% anomalous events and 9/244 = 3.8% planets, i.e., both identical within Poisson errors.
For 2018, AnomalyFinder2.0 has already been run on the 21 KMTNet fields with lower cadence, Γ < 2 hr −1 , covering 84 deg 2 and yielding a total of 173 anomalous events, which are distributed among the five classifications as (17,4,19,126,7). These include nine published planets and three in preparation. However, among the nine published planets, three have ambigu- ous or larger-error mass-ratio measurements and so, are not suitable for studying the planet-host mass-ratio function. Therefore, we may expect that after lower-cadence AnomalyFinder output is fully studied, there will be a total of 30-35 planets in the 2018 sample. This will be comparable in size to the largest previous study (Suzuki et al. 2016), which included 22 planets from six years of the Microlensing Observations in Astrophysics II (MOA-II) survey, augmented by 8 planets found in two earlier survey/followup studies (Gould et al. 2010;Cassan et al. 2012). However, the AnomalyFinder sample will subsequently be expanded to cover at least four years, as we continue to publish all planets 2016-2019.

Observations
As described in Section 1, all of the planetary (or potentially planetary) events that are presented in this paper were initially identified by applying the AnomalyFinder2.0 algorithm  to the 843 events that were originally found by the KMTNet EventFinder and AlertFinder systems in the prime fields during 2018. As described by Hwang et al. (2022), when available, we use data from independent alerts from the Optical Gravitational Lensing Experiment (OGLE) and MOA to vet the anomalies for systematics (otherwise, we study these anomalies at the image level). We also include OGLE and MOA data in the analysis of the events. These were taken using the OGLE 1.3m telescope with 1.4 deg 2 field of view at Las Campanas Observatory in Chile, and the MOA 1.8m telescope with 2.2 deg 2 field of view at Mt. John Observatory in New Zealand. The OGLE and MOA data analyzed here are in the I band and a broad, customized, R-I filter, respectively. Table 1 gives basic observational information about each event. Column 1 gives the event names in the order of discovery (if discovered by multiple teams), which enables cross identification. However, in most of the rest of the paper, we refer to events only by the name given by the group who made the first discovery. The nominal cadences are given in column 2, and column 3 shows the first discovery date. The remaining four columns show the event coordinates in the equatorial and galactic systems. Events with OGLE names were originally discovered by the OGLE Early Warning System (Udalski et al. 1994;Udalski 2003).
To the best of our knowledge, there were no ground-based follow-up observations of any of these events. One event, OGLE-2018-BLG-0932, lies in the field of the UKIRT Microlensing Survey (Shvartzvald et al. 2017), and we make use of these data to determine its source color. This survey employs a 3.8m telescope in Hawaii, with an effective field of view of 0.2 deg 2 , to observe in the H and K bands. OGLE-2018-BLG-0932 was also observed by the Spitzer space telescope, but the analysis of the resulting data is beyond the scope of the present work and will be presented elsewhere.

Preamble
We begin by describing the light-curve analysis methods and notations that are common to all events. All of the events in this paper appear, to a first approximation as simple 1L1S light curves, which can be described by three Paczyński (1986) parameters, (t 0 , u 0 , t E ), i.e., the time of lens-source closest approach, the impact parameter in units of θ E and the Einstein timescale, where M is the lens mass, π rel and µ rel are the lens-source relative parallax and proper-motion, respectively, and µ rel ≡ |µ rel |.
Here nLmS means n lenses and m sources. In addition, to these three non-linear parameters, two flux parameters, ( f S , f B ), are required for each observatory, representing the source flux and the A&A proofs: manuscript no. ms blended flux that does not participate in the event. Note, however, that these are linear parameters, which can be determined by regression after the model is specified by the non-linear parameters. We then search for "static" 2L1S solutions, which require four additional parameters (s, q, α, ρ), i.e., the planet-host separation in units of θ E , the planet-host mass ratio, the angle of the source trajectory relative to the binary axis, and the angular source size normalized to θ E , i.e., ρ = θ * /θ E .
We conduct this search in two phases. In the first phase, we search on a 2-dimensional (2-D) grid. For each (s, q) pair, we construct a magnification map following Dong et al. (2009b). We then conduct a downhill search using the Monte Carlo Markov chain (MCMC) technique. We seed this search with the 1L1S solution for the Paczyński (1986) parameters, (t 0 , u 0 , t E ). We use the approach of Gaudi et al. (2002) to find the seed for ρ. For α we seed at a grid of values around the unit circle. This procedure yields a χ 2 map on the (s, q) plane, which we use to identify one or several local minima.
In the second phase, we refine the best solution at each local minimum by allowing all seven parameters to vary in the MCMC.
In this analysis, we often make use of a modified version of the heuristic analysis introduced by Hwang et al. (2022). If a brief anomaly at t anom is assumed to be generated by the source crossing the planet-host axis, then Hwang et al. (2022) suggested analytic estimates for (s, α) of where u 2 anom = τ 2 anom + u 2 0 , τ anom = (t anom − t 0 )/t E , and where ∆s (i.e., half the difference between the two solutions) generally cannot be determined from by-eye inspection. In the great majority of cases, s † + > 1 corresponds to anomalous bumps and s † − < 1 corresponds to anomalous dips. This formalism was designed to reflect the "inner/outer" degeneracy (Gaudi & Gould 1997) whereby the source passes the planetary caustic(s) on the side closer to (or farther from) the central caustic. However, following the work of Herrera-Martin et al. (2020) and Yee et al. (2021), it was already recognized to have somewhat wider application.
In the course of the present investigation, in which we systematically applied Equation (2) to all 10 events, we encountered OGLE-2018-BLG-1647, which proved to be the "Rosetta Stone" that unified the "inner/outer" degeneracy for planetary caustics (Gaudi & Gould 1997) with the "close/wide" degeneracy for central caustics (Griest & Safizadeh 1998), as conjectured by Yee et al. (2021). For this event, the formula for s † + in Equation (2) proved to be a better approximation to the geometric mean of the two empirically derived solutions, s ± , i.e., s † = √ s + s − , rather than the arithmetic mean, i.e., s † = (s + + s − )/2. This fact immediately led to several realizations. First, this reformulation did not contradict any of the four cases examined by Hwang et al. (2022), nor the many other cases examined in the current work, because for these ∆ ln s ≡ (1/2) ln(s + /s − ) was always small, ∆ ln s ≪ 1. In this limit, for which Equation (2) worked quite well, the arithmetic and geometric means differ by only ∼ (∆ ln s) 2 /2, which is generally too small to notice. Second, the mathematical representation of this reformulation, is equivalent to the usual expression for the "close/wide" degeneracy, s − = 1/s + , provided that s † → 1. Moreover, be-cause Griest & Safizadeh (1998) derived this relation in the limit of central caustics, i.e., high-magnification events for which u anom ≪ 1, the limit s † → 1 does indeed apply to this case. Third, what made OGLE-2018-BLG-1647 a "Rosetta Stone" is that the geometric mean of Equation (3) applied, even though s † 1 (contrary to the "close/wide" limit). Fourth, the several historical examples that inspired Yee et al. (2021) to suggest unification were all "inner/outer" degeneracies in which one of the two solutions had the source passing between the central and planetary caustics, while the other had it passing outside the planetary wing of a resonant caustic. That is, one solution appeared more closely related to the "inner/outer" degeneracy and the other to the "close/wide" degeneracy 2 . The pair of solutions were dubbed "inner/outer" primarily because both solutions had the same logarithmic sign, (ln s + )(ln s − ) > 0. This had already indicated a continuous degeneracy to Yee et al. (2021). However, in the course of this (and other) work, we noted additional cases with similar topologies, but for which (ln s + )(ln s − ) < 0 (as in the "close/wide" limit), but for which Equation (3) remained a better approximation than the s − = 1/s + prediction of Griest & Safizadeh (1998). We regarded this as further evidence for a continuum of (s − , s + ) degeneracies from inner/outer (s † < 1, minor-image caustics), through close/wide (s † ≃ 1, central and resonant caustics), to outer/inner (s † > 1, major-image caustics).
Subsequently, Ryu et al. (2022) have provided uniform notation for this formalism in their Equations (2)-(7). We follow their conventions here. In particular s † ± (with "±" subscript) denotes the theoretical prediction of Equation (2), while s † (without subscript) denotes the geometric mean of the two empirical solutions, whose offset is characterized by ∆ ln s, Hwang et al. (2022) also introduced an estimate of the mass ratio q for dip-type anomalies, which is ultimately based on the theoretical analysis by Han (2006): where ∆t dip is the full duration of the dip. Ryu et al. (2022) noted that this expression can be rewritten in terms of "direct observables": where they pointed out that ∆t dip and δt anom can be read directly off the light curve, while t eff ≡ u 0 t E ≃ FWHM/ √ 12 for even moderately high magnification events, A max > ∼ 5. Indeed, Yee et al. (2012) had already pointed out that t q = qt E is also an invariant for high-magnification events.
In some cases, we investigate whether the microlens parallax vector (Gould 1992(Gould , 2000(Gould , 2004 can be constrained by the data. Note that if this quantity can be measured, then by combining Equations (1) and (7) one can infer the lens and mass and distance, where π S is the parallax of the source, which usually is approximately known. However, even if π E cannot be measured (e.g., it is consistent with zero at 1 σ), it can significantly constrain (M, D L ) after imposing priors from a Galactic model, provided that the error ellipse on π E is sufficiently small, at least in one dimension (see the Appendix in Han et al. 2016).
To model the parallax effects due to Earth's orbital motion, we add two parameters (π E,N , π E,E ), which are the components of π E in equatorial coordinates. Because these effects can be mimicked by those due to lens orbital motion (Batista et al. 2011;Skowron et al. 2011), we always add (at least initially) two parameters γ = [(ds/dt)/s, dα/dt], where sγ are the first derivatives of projected lens orbital motion at t 0 , i.e., parallel and perpendicular to the projected separation of the planet at that time, respectively. In order to eliminate unphysical solutions, we impose the constraint on the ratio of the transverse kinetic to potential energy Dong et al. 2009a), Note that while orbits are only unbound if β > 1, we impose a slightly stronger constraint because it is extremely rare for planets to be in such high-eccentricity orbits and observed at the right orientation and epoch to yield β > 0.8. It often happens that γ is neither significantly constrained nor significantly correlated with π E . In these cases, we suppress these two degrees of freedom.
Very frequently, including several cases in this paper, the parallax contours in the π E plane take the form of elongated ellipses (Gould et al. 1994) with the orientation angle of short axis, ψ, being approximately aligned with the projected position angle of the Sun, ψ ⊙ , at the peak of the event, t 0 . That is, ψ ≃ ψ ⊙ . This is because, for events with t E ≪ 1 yr, Earth's acceleration is approximately constant, under which condition lens-source motion along the direction of acceleration gives rise to much more pronounced effects than does the transverse motion (Smith et al. 2003). When this occurs, it can be substantially more informative to characterize π E = (π E, , π E,⊥ ) in terms of these two components. For example, unless ψ is closely aligned with one of the cardinal directions, σ(π E, ) can be much smaller than either σ(π E,N ) or σ(π E,E ). For reference, we note that the (Gaussian) likelihood associated with the parallax measurement can be expressed as, where π E,0 is the best fit, b ≡ c −1 , c is the covariance matrix, and where we have written the determinant of this matrix explicitly in terms of its eigenvectors in order to make contact with future applications. As pointed out by Gaudi (1998), 1L2S events can mimic 2L1S events, particularly if there are no sharp caustic-crossing features in the light curve. If ∆χ 2 = χ 2 (1L2S) − χ 2 (2L1S) is strongly negative, then we conclude that the event is 1L2S, and we eliminate it from consideration. If we test for 1L2S and find The anomaly is a dip that is centered at 8298.7, which is detected at ∆χ 2 = χ 2 (1L1S) − χ 2 (2L1S) = 69. As in all 10 light-curve figures in this paper, we show the full light curve and anomaly region in separate panels, we show the caustic topologies in one or more insets, we show residual panels for indicated models, and we color the data points by observatory, as indicated in the legend. that ∆χ 2 is strongly positive, we remark that such solutions are ruled out. If 1L2S and 2L1S have either competitive or roughly comparable χ 2 we report both solutions. The former class of events are ambiguous in nature and cannot be included in planetary catalogs, nor certainly in mass-ratio function studies. However, we report such events because it may be possible in the future to resolve the degeneracy for some of them using auxiliary data.
We carry out 1L2S modeling by adding at least 3 parameters (t 0,2 , u 0,2 , q F ) to the 3 Paczyński (1986) parameters. These are the time of closest approach and impact parameter of the second source and the ratio of the second to the first source flux in the I-band (Hwang et al. 2013). If either lens-source approach can be interpreted as exhibiting finite source effects, then we must add one or two further parameters, i.e., ρ 1 and/or ρ 2 . And, if the two sources are projected closely enough on the sky, one must also consider source orbital motion (e.g., Hwang et al. 2018b).

OGLE-2018-BLG-1126
The KMTC data exhibit a systematic decline relative to the 1L1S model centered on 8298.7. See Figure 1. The formal significance of this deviation is modest: ∆χ 2 = χ 2 (1L1S) − χ 2 (2L1S) = 69. Moreover, because the coverage of the anomaly is incomplete, one must be concerned that this deviation is due to some systematic effect. The main potential cause of such an effect would be the Moon, which was full when it passed through the bulge (about 11 • north of the event) roughly 36 hours before the anomaly. There is a well-known mechanism for the Moon to induce a spurious excess (though not deficit) in the tabulated flux, which generates many false alerts of short timescale events A&A proofs: manuscript no. ms by the EventFinder (Kim et al. 2018a): the higher background pushes a bright star above the pixel well depth, causing charge to bleed into a column and so pollute the photometry of fainter stars that are downstream in the same column. These bleeds are often invisible on normal displays of the original images because the stretch is generally too weak to detect them. However, they are easily visible on difference images, for which the stretch can be made much stronger. We carefully examine the difference images throughout the night and find no such signatures. Another possibility is that the Moon caused excess flux on the previous night when it resulted in much higher background (13000 versus 4000), thus affecting the overall light-curve model, thereby giving the appearance of an anomaly on the following night. However, we see no evidence for bleeds on the previous night. Thus we conclude that the anomaly is real. Adopting Paczyński (1986) parameters (t 0 , u 0 , t E ) = (8298.17, 0.0083, 53.3 day) and light curve features (t dip , ∆t dip ) = (8298.8, 1.2day), the heuristic approach outlined in Section 3.1 yields τ anom = +0.0118, α = 35 • , s † − = 0.993, and q = 7 × 10 −4 . The grid search returns two local minima. After refining these as described in Section 3.1, we find that they generally agree with heuristic prediction. See Table 2. The main discrepancy is in α (29 • versus 35 • ), which is mainly due to the difficulty of judging the center of dip from the incomplete light curve. Of particular note is the striking agreement of s † ≡ √ s + s − = 0.992 (compared to the prediction s † − = 0.993). Thus, although this degeneracy would normally be considered as a classic example of the "close/wide degeneracy" for central and resonant caustics because s close s wide ≃ 1, the prediction of the s † formalism (derived in the limit of planetary caustics) is actually 10 times more accurate 3 . Note that there is essentially no constraint on ρ for this planet.
Due to the faintness of the source, we do not attempt a parallax analysis.
While we have concluded that the planet is real, it may not be suitable for mass-ratio function studies. From Table 2, we see that the 1 σ error in log q is 0.28 dex, which corresponds to a factor of ∼ 1.9. The goal of the present paper is not to impose a boundary for this parameter, but rather to present a comprehensive account of all planets that meet much broader criteria in order to provide a basis for such choices in future analyses of 3 Note that in the geometric-mean formalism of Equation (3), s † = √ s + s − < 1, which conforms to the minor-image caustic morphology of the light curve. However, this is not true of the arithmetic-mean prediction s † = (s + + s − )/2 > 1, This was a significant puzzle for us when we initially analyzed this event, but was resolved after analyzing the "Rosetta Stone" event OGLE-2018-BLG-1647. See Section 3.1 the mass-ratio function. However, we remark that it is at least questionable whether this planet will enter such studies.
We note that although this planet meets the q < 2 × 10 −4 selection criterion of Hwang et al. (2022), it was not included in their sample. This is because it was detected by AnomalyFinder2.0 ), but not Anoma-lyFinder1.0 (Zang et al. 2021b), which was the basis of the Hwang et al. (2022) study.

KMT-2018-BLG-2004
The anomaly in Figure 2 consists of a short bump, which is traced by both KMTS and KMTC data, centered on t anom = 8242.7, when the Moon was about 10 • north of the event. While in this case, the Moon was 4 days past full (so the background at passage was 5000, compared to 13000 for OGLE-2018-BLG-1126), it is far more plausible that the Moon would cause a bump in the light curve, rather than a dip. Indeed, given that the bump is continuous across two observatories separated by 8000 km, it is difficult to conceive of any other source of systematics. However, we again carefully examine the subtracted images and find no evidence of bleeding columns. Hence, we again conclude that the anomaly is due to microlensing.
Using the above t anom , combined with the 1L1S parameters (t 0 , u 0 , t E ) = (8239.17, 0.23, 31 day), the heuristic formalism (see Equation (3)) predicts s † + = 1.14 and α = 244 • . The grid search returns only two solutions, which after refinement agree quite well with these predictions. See Table 3. In particular, s † = √ s inner s outer = 1.14. The anomaly is detected at Given that the anomaly is a featureless bump, it is essential to check whether it can be explained by a binary source (1L2S) model. From Table 3, we see that such models are disfavored   by ∆χ 2 = 14.8, which is substantial, though not overwhelming, evidence in favor of 2L1S. In the 1L2S model, the best fit value of the flux ratio is q F = 2.2 × 10 −3 , corresponding to a magnitude difference of ∆I = −2.5 log(q F ) = 6.6 magnitudes. We will show in Section 4.2 that the source lies about 3.6 mag below the clump. Hence, the putative source companion would have an absolute magnitude of M I,comp ∼ 10. Such stars are common, so the 1L2S solution cannot be regarded as implausible on these grounds.
The 1L2S model makes the definite prediction that the "bump" should be basically invisible in the V band. That is, the source companion should have (V − I) comp,0 ∼ 3.3 whereas (as we show in Section 4.2), (V − I) S,0 ∼ 0.7. Thus, the relative amplitude of the bump should be 10 0.4(3.3−0.7) = 11 times smaller in V than I. This implies that the V-band light curve should follow the I-band light curve for 2L1S but should follow the 1L1S curve for 1L2S. See Hwang et al. (2018b). Unfortunately, the V data are not good enough to test this prediction. Of the four potential data sets, (KMTC & KMTS)×(BLG01 & BLG41), only KMTS BLG01 provides useful information. This has only one V-band point over the bump. The point lies almost exactly on the 2L1S curve. However, it is only 0.5 σ from the 1L1S curve, due to the relative large V-band error bars.
Thus, the only strong argument against the 1L2S solution is that ∆χ 2 = 14.8. (If we incorporate the V-band test, just mentioned, this becomes ∆χ 2 = 15.1). We consider that the planet solution is strongly preferred, but we cannot rule out the binarysource solution unconditionally.
The event is moderately long and has good photometry, so we attempt to fit it for parallax. Figure 3 shows the parallax contours for two of the four cases, namely the "inner" solution with u 0 > 0, and u 0 < 0.
The parallax fit reveals interesting information. The basic form is of a so-called 1-dimensional (1-D) parallax measurement, which occurs because Earth's acceleration toward the projected position of the Sun (ψ ⊙ = 96.7 • north though east) is roughly constant over the relatively short duration of the event. See Section 3.1. Formally the error ellipses have an access ratio of ∼ 12. The two "lobes" toward the north and south imply that the measurement is subject to the so-called jerk-parallax degeneracy (Gould 2004;Park et al. 2004). While these are striking to the eye, in part because of their large values, π E ∼ 2, they are favored by only ∆χ 2 ∼ 4, which would have marginal statistical significance even if the errors could be treated as Gaussian. That is, even in this case, their weight would be overwhelmed by the Galactic priors in a Bayesian analysis, which heavily disfavors For both events, these contours have very large axis ratios that are characteristic of so-called 1-D parallax measurements. We argue in the text that only the short-axis information in these contours is reliable and reduce them to truly 1-D constraints. See Equations (10) and (11) and Sections 3.3 and 3.5.
such large parallax values. Moreover, in addition to having larger statistical errors along the long axis of the ellipse, the result is also more subject to systematic errors because the information is coming primarily from the wings of the light curve (Smith et al. 2003;Gould 2004).
The actual information in these contours comes from their small width, not their best-fit values. In principle, if these narrow contours were displaced from the origin, as in the first microlensing planet with such features, OGLE-2005-BLG-071 (Dong et al. 2009a), then they would be strong evidence for a minimum value π E ≥ π E, , even if the exact value was not determined. However, in the present case, the contours pass through the origin, so the result has less discriminatory value.
Nevertheless, we proceed to extract the essence of the parallax information, while suppressing possible systematic effects, by retaining the short-axis information σ(π E, ), while setting σ(π E,⊥ ) → ∞, and using the fact that the contours pass through the origin. Noting that the contours "bend" at the origin, we adopt for the four cases (sgn(u 0 ) = ±; sgn(π E,N ) = ±)), and we set π E,0 = 0. Because this is a 1-D constraint (albeit on a 2-D space), we substitute 2πσ(π E, )πσ(π E,⊥ ) → √ 2πσ(π E, ). Note that, by construction, b is a degenerate matrix.

OGLE-2018-BLG-1647
Figure 4 shows a pronounced bump ∆τ = −0.083 before the peak. The grid search returns two local minima, whose refinements are shown in Table 4. Traditionally, this would be interpreted as the close/wide degeneracy in which the source passes similar-looking central caustics (Figure 4), for which we would expect the geometric mean to be unity, compared √ s close × s wide = 1.07, for these two reported solutions. On the other hand, adopting u 0 = 0.105, the heuristic analysis of Section 3.1 yields α = −52 • and s † + = 1.07, i.e., essentially identical to the geometric mean. Hence, this event is much closer mathematically to the inner/outer degeneracy (derived in the limit of planetary caustics) than it is to the close/wide degeneracy (derived in the limit of central and planetary caustics).
Note that the arithmetic mean of Equation (2) would yield (s + +s − )/2 = 1.11. As we discussed in some detail in Section 3.1, it was the fact that the geometric mean worked better than the arithmetic mean that led us to adopt Equation (3) to unify the inner/outer and close/wide degeneracies.
Because the wide/inner model is preferred by ∆χ 2 = 17, we adopt it over the close/outer model. In any case, the two models have essentially identical mass ratios, q ≃ 0.010. We also search for 1L2S models, but find that they are disfavored by ∆χ 2 = 28. See Table 4. Hence, they are decisively rejected.
Due to the faintness of the source, we do not attempt a parallax analysis.
OGLE-2018-BLG-1647 is one of three previously known planets that are listed by Hwang et al. (2022) as "in preparation" but are analyzed here for the first time.

OGLE-2018-BLG-1367
Figure 5 shows a flattened, or perhaps slightly depressed peak. A natural way to produce a flattened peak is a 1L1S geometry with finite source effects as the lens transits the face of the source, i.e., so-called finite-source/point-lens (FSPL) events. We search for such a model, but it produces a relatively poor fit, χ 2 (FSPL) − χ 2 (2L1S) = 65. In addition, the FSPL fit parameters (t E , ρ) = (22.0 day, 0.048), would imply an extraordinarily long source self-crossing time (t * = 1.1 day), given that the source is a turnoff star (see Section 4.4). Hence, the Einstein radius would be θ E ≃ 16 µas, while the proper motion would be an extraordinarily slow µ rel ≃ 0.27 mas yr −1 , with prior probability p = 8 × 10 −5 . See Equation (14), below. That is, we expect only about one event with such a low proper motion during the five  years of KMT normal observations, and this one event would have only a few percent chance of giving rise to finite-source effects (thus enabling its low µ rel to come to our attention). By contrast, the 2L1S models (Table 5) fit the data quite well and do not require exceptional physical parameters. The flattening (or depression) near the peak is then explained by the source passing roughly perpendicular to the planet-host axis on the opposite side of the planet, a region that is characterized by a negative magnification deviation relative to 1L1S.
For perpendicular trajectories, s † − = ( 4 + u 2 0 − u 0 )/2 → 0.987. Hence, the geometric mean of the two solutions (0.981) is slightly closer to this value than it is to unity (the close/wide prediction). This tends to confirm our conjecture that Equation (4) is the correct generalization of the s † formalism, even though the event is qualitatively well described by the "close/wide" degeneracy This is another massive planet, q ≃ 3.4 × 10 −3 , i.e., 3.5 times larger than the Jupiter/Sun ratio.
Because the source is relatively bright and the photometry is good, we attempt to measure π E . Figure 3 shows the parallax contours for one of the four solutions, namely the close solution for u 0 > 0. As in the case of KMT-2018-BLG-2004, the contours are highly elongated (1-D parallax) with two lobes, indicating that the event is subject to the jerk-parallax degeneracy. However, contrary to that case, the contours do not pass through the origin, but rather cross the π E,N axis at π E,E ≃ 0.165, which is 4 times larger than the error. Hence, this parallax measurement contains significant information.
To extract this information, we follow similar procedures to those of Section 3.3, but with some difference. First, contrary to the previous case, there is essentially no bend between the positive and negative π E,N regimes. Second, the contours are essentially identical for positive and negative u 0 . Third, as mentioned above, the contours do not pass through the origin. The first two of these imply that there is one regime: (σ(π E, ), ψ)) = (0.0396, 87.30 • ). To implement the third within the framework of Equation (10), we rotate the measured π E, ,0 = 0.165 to Equatorial coordinates. π E,0 (N, E) = π E, ,0 (cos ψ, sin ψ) = (0.008, 0.165).
3.6. OGLE-2018-BLG-1544 Figure 6 shows a dip starting near the peak, followed by a bump centered at t bump = 8352.7. If the latter is attributed to the source crossing the planet-host axis on the planet side, then the heuristic formalism gives α = 208 • and s † + = 1.03. The angle, in particular, implies that the dip is generated by passage along one of the Fig. 6. Light curve and model for OGLE-2018-BLG-1544, The anomaly is a long dip near the peak followed by a shorter bump. The heuristic analysis is anchored in the latter, which implies a shallow source trajectory α = 208 • . The dip is then understood as the lateral passage of one wall of a central caustic. See inset. long sides of the central caustic due to a low-mass (but not necessarily planetary) companion. In principle, there might be other geometries. However, the grid search finds only two local minima, which correspond to the close and wide versions of the one anticipated above, with q = 0.019 and q = 0.016, respectively, the former being favored by ∆χ 2 = 3. See Table 6. Hence, this is another very massive planet (under the planet definition q < 0.03).
Due to the faintness of the source, we do not attempt a parallax analysis.
Because this event has a major-image "bump generating" caustic topology, and despite the fact that it does not exhibit the classical "isolated bump" morphology that would normally induce concerns about a possible binary-source interpretation, we fit for 1L2S models. We find that ∆χ 2 = χ 2 (1L2S) − χ 2 (2L1S) = A&A proofs: manuscript no. ms Fig. 7. Light curve and model for OGLE-2018-BLG-0932. The anomaly is a bump centered at 8273.5. Unlike most smooth, isolated bumps, this one is due to a source passage over a minor-image caustic, with the smoothness due to the fact that source is very large compared to the caustic. See inset. Among the 10 events analyzed here, this is the only one for which the source-size parameter ρ = θ * /θ E is precisely measured. Table 6. Hence, while the planetary interpretation is favored, there is a significant possibility that the anomaly is actually due to a binary source.

See
OGLE-2018-BLG-1544 is one of three previously known planets that are listed by Hwang et al. (2022) as "in preparation" but are analyzed here for the first time.

OGLE-2018-BLG-0932
OGLE-2018-BLG-0932 is a good example of a case for which the heuristic formalism gives relatively imprecise guidance. The 1L1S approximation has (t 0 , u 0 , t E ) ≃ (8301.1, 0.85, 27 day), and t anom ≃ 8273.5, i.e., τ anom = −1.02. These imply s † ± = 1.20±0.66 and α + = 320 • (or α − = 140 • ). The fact that the anomaly is a "bump" rather than a "dip" leads one to expect that this is major image perturbation, so s † ∼ 1.86, α = 320 • . In fact, however, a full grid search shows that there is only one solution, for which the bump is due to the source transiting a triangular caustic from a minor-image perturbation and for which the heuristic prediction is s † − ∼ 0.54, α = 140 • . Comparison to Table 7 shows that s ≃ s † − , as expected for cases with no inner/outer degeneracy. However, α differs from the prediction by 5 • , which is much larger than any of the other cases examined here or the 11 cases to which the heuristic analysis was systematically applied by Hwang et al. (2022) and Ryu et al. (2022). The reason is that the heuristic analysis implicitly assumes that the anomaly is centered on the planet-host axis. This basically holds for majorimage planetary perturbations, for dip-like minor-image planetary perturbations, and even for minor-image caustic crossings for the cases of very small q (because the caustics are then very close to the minor-image axis). However, for the present case, q ∼ 10 −3 , the caustic is 0.1 Einstein radii from the axis (see Figure 7, i.e., at an angle sin −1 (0.1/u anom ) = 4 • relative to this axis, which accounts for the "error" in the heuristic prediction. The results shown in Table 7 have blending fixed to zero, specifically using the baseline source flux as determined by OGLE. A free fit to blending gives f B / f base = −0.30 ± 0.09, with an improvement ∆χ 2 = 6.7. For such a bright source, such large negative blending cannot be the result of unmodeled fluxes from unresolved stars. In principle, it could be a statistical fluctuation (Gaussian probability p = 4%), but is more likely due to lowlevel systematics or source variability, or possibly to unmodeled physical effects, such as parallax.
From the present perspective, we simply impose zero blending, while noting that the parameters do not change much for the negative blending solutions. For example, the value of q rises from 1.19×10 −3 to 1.26×10 −3 . We do not investigate parallax solutions here because this event has Spitzer parallax observations under a large program that was outlined by Yee et al. (2015). These will be analyzed elsewhere.
We searched for 1L2S solutions, but find that these are ruled out by ∆χ 2 = 564.
OGLE-2018-BLG-0932 is one of three previously known planets that are listed by Hwang et al. (2022) as "in preparation" but are analyzed here for the first time.

OGLE-2018-BLG-1212
The light curve for this event shows a strong asymmetry due to parallax, even when the anomaly is removed. Hence, contrary to our usual procedures, we fit for parallax prior to searching for 2L1S solutions. Both the 1L1S and 2L1S models in Figure 8 include parallax. We can still carry out a heuristic analysis using the 1L1S parallax-model parameters (t 0 , u 0 , t E ) = (8393.76, 0.014, 51 day), together with the midpoint and width of the dip: t anom = 8394.1 and ∆t dip = 0.75 day. These yield s † − = 0.992, α = 64 • , and q = 7 × 10 −4 . These should be compared with the results from the full grid search shown in Table 8, i.e., s † = √ s close × s wide = 0.993, α = 63 • , and q = 12 × 10 −4 . For the record, we note that in our initial 2L1S fit, we obtained a very well-localized solution at π E (N, E) = (0.534, 0.550). However, we found that the jerk-parallax degeneracy formalism (Equations (7)-(9) of Park et al. 2004) predicts another solution 4 at π E (N, E) = (−0.404, 0.550), and numerical 4 The event peaked almost at quadrature, i.e., ψ = 101 • in this formalism. Moreover the projected position of the Sun at this time is only 0.2 • from due west, implying that π E,⊥ ≃ −π E,N . Thus, to an excellent approximation, Park et al. (2004) Equation (9)   . Fig. 8. Light curve and model for OGLE-2018-BLG-1212. The anomaly is a dip centered at 8394.1, which is traced by both KMTA and MOA data. The event has a very strong parallax signal and large parallax parameter, π E = 0.767 ± 0.019, almost certainly implying a nearby lens. See Sections 3.8, 4.7, and 5.7. tion (7) becomes π ′ E,N = −π E,N + (4/3)(au/v ⊕ t E )) sin(−β) = −0.534 + 0.126 = −0.408, very close to the more exact calculation. investigation then showed that this was recovered to high precision. See Table 8, While this second set of solutions is disfavored by ∆χ 2 ∼ 11, we keep track of its potential implications because the π E,⊥ (≃ π E,N ) parameter is among the most sensitive to subtle systematic errors.
The wide solution is favored by ∆χ 2 = 3, which is far below the level that would be required to distinguish between the two solutions. However, the parameters (apart from s) of the two solutions are essentially identical.
This value tends to favor lens distances D L ∼ 1-2 kpc. That is, ignoring the peculiar motions of the lens relative to the disk and of the source relative to the bulge,ṽ lsr (l, b) ≃ [(D S /D L −1) −1 v rot , 0] for a flat rotation curve with rotation speed v rot = 235 km s −1 . This would imply D L ∼ (1 + v rot /ṽ lsr ) −1 D S → 1.4 kpc. Because the lens and source peculiar motions cannot truly be ignored, and because there is more phase space at larger distances, this argument is only suggestive. Nevertheless, we discuss its potential implications in Sections 4.7 and 5.7.

KMT-2018-BLG-2718
From Figure 9, this event does not, at first sight, appear to be planetary in nature. The anomaly is a dip near the peak of the event, which is of very long duration t dip ∼ 20 days. Estimating Article number, page 11 of 22 A&A proofs: manuscript no. ms  9. Light curve and model for KMT-2018-BLG-2718. The anomaly is a dip near the peak, which is flanked by two bumps. This morphology is the classic signature of the planet/binary degeneracy identified by Han & Gaudi (2008). See insets. In this case, the planetary interpretation is favored by ∆χ 2 = 12.7. The invariant parameter t q = qt E = 3.16 ± 0.16 days would imply non-planetary mass ratios (by our definition, q > 0.03), unless t E 100 days. In fact, the fits imply much longer timescales. See Table 9.
t eff ∼ 10 days and ∆t anom ≪ t eff (so s † − ∼ 1), we can expect 5 from Equation (6) that t q ≃ 2.5 days, so that this event would only meet our planet definition q < q max = 0.03 provided that t E > ∼ t q /q max ∼ 83 days. Nevertheless, the morphology of this very faint (I peak ∼ 18.7) event does suggest such a long duration. This emphasizes the importance of carefully reviewing all detections of the AnomalyFinder even if they do not look planetary at first sight.
The grid search indeed returns a wide/close pair of planetary solutions with q = 0.020 and q = 0.014 that are in accord with the above heuristic analysis, i.e., with timescales t E ∼ 160 days and 230 days, respectively. However, it also returns a pair of binary solutions with q > ∼ 0.6. See Table 9. The planetary solutions are favored by ∆χ 2 = 12.7. If the statistics could be assumed to  (2008) planet binary degeneracy (see insets), but contrary to that case, the planetary interpretation is not decisively favored (see Table 10). Therefore, the lens companion cannot be claimed as a planet.
be Gaussian, then this would decisively resolve the planet/binary ambiguity. However, given the quality of the data and the general inapplicability of Gaussian statistics to microlensing data, we would rather regard this as "basically resolved".
Due to the faintness of the source, we do not attempt a parallax analysis.

KMT-2018-BLG-2164
Figure 10 shows a dip near the overall peak, flanked by roughly equal bumps. In principle, this could be caused by the source passing roughly perpendicular to the planet-star axis on the opposite side of the planet, similarly to OGLE-2018-BLG-1367. The grid search indeed returns a close/wide pair that corresponds to this geometry. But it also finds a second pair of minima, in which the source passes diagonally outside a Chang-Refsdal caustic. Refinement of these minima indicate a planet-versusbinary degeneracy, i.e., q ∼ 0.001 versus q ∼ 0.15, which was predicted by Han & Gaudi (2008). The planetary solution is favored by ∆χ 2 = 3.5, but this is far below the level what would be required to confidently claim a planet. See Table 10. This object is presented here because our protocols demand that we include all companions that are consistent with being planetary, even if this designation cannot be confirmed. In this case, the planetary and binary solutions predict similar source fluxes and there are no proper-motion estimates (because there is no ρ measurement). Hence, future adaptic optics (AO) observations cannot distinguish between the solutions. This could only be done using RV follow-up observations on extremely large telescopes (ELTs), or possibly even larger telescopes that will operate in the more distant future. Note, however, that even if this proves to be a planet, the uncertainty in log q is 0.2 dex, corresponding to a factor 1.6. This large uncertainty is related to the fact that the improvement relative to 1L1S is only ∆χ 2 = 89.
Due to the faintness of the source, we do not attempt a parallax analysis.

OGLE-2018-BLG-1554
As shown in Figure 11, the light curve exhibits a long-term deviation over the peak, which is relatively small, but nonetheless we find to be statistically significant at ∆χ 2 = 413. The grid search returns two pairs of solutions, one being a planetary pair with q ∼ 0.025 and the other being a binary pair with q ∼ 0.075. In addition to these four solutions, we find a 1L2S solution. All three classes have a member that lies within the overall minimum at ∆χ 2 < 1.4, so all three are "equally good" in this sense. See Table 11.
Only the planetary solutions have a ρ measurement, ρ ∼ 0.03, corresponding to t * ≡ ρt E ∼ 0.4 days. In Section 4.10, we will show that θ * ≃ 0.93 µas. Hence, if the planetary solution were correct, then µ rel = θ * /t * ∼ 0.8 mas yr −1 . As we explain just below in Section 4, the fraction of events with such low proper motions is p < (µ rel /6.4 mas yr −1 ) 3 ≃ 2×10 −3 . Thus, we consider the planetary solution to be extremely unlikely.
In any case, given that the planetary solution cannot (at present) be distinguished from the binary-lens and 1L2S solutions, this event cannot be included in (present-day) mass-ratio function studies.
For completeness, we remark that if future AO followup observations confirm the very low µ rel < ∼ 1 mas yr −1 predicted by the planetary solutions, this would constitute strong evidence (though not proof) that it was correct. However, such confir- Fig. 11. Light curve and model for OGLE-2018-BLG-1554. The anomaly is characterized by weak deviations both before and after the peak. Like the previous two events, this one is subject to the Han & Gaudi (2008) planet/binary degeneracy (see insets), but even more severely (see Table 11). In addition, there is a severe 1L2S/2L1S degeneracy. See Table 11. Therefore, it is not established that the lens has a companion, and even if it does, this companion cannot be claimed as a planet. mation would face extreme observational challenges, even with next-generation 30m class telescopes.
The first point is that if the planetary solution is correct, then θ E ∼ 30 µas, and so π rel ∼ 0.11 µas (M/M ⊙ ) −1 . That is, the lens will be invisibly faint unless the lens and source are within D LS ≡ D S − D L ≃ D 2 S π rel /au < ∼ 100 pc, which is itself highly improbable. Moreover, it means that the "correction" from the measured geocentric to the relevant heliocentric proper motion, ∆µ = µ rel,hel − µ rel = v ⊕⊥ π rel /au will be extremely small. Here v ⊕⊥ (N, E) = (−3.7, +13.7) mas yr −1 is the projected velocity of Earth at t 0 . That is, |∆µ| ∼ 0.02(M/0.075 M ⊙ ) −1 mas yr −1 , so that µ rel,hel ≃ µ rel = 0.8mas yr −1 . Given the faintness of the lens, this A&A proofs: manuscript no. ms would require waiting of order 3 decades even with ELTs. Thus, even in the unlikely case that the planetary solution is correct, the prospects for confirming it are distant at best.

Source Properties
For a substantial majority of planetary microlensing events that have been reported in the past, ρ was measured. Hence, if the angular source size, θ * , could be determined, it yielded θ E and µ rel : Then, if π E could also be measured, one could directly infer the lens mass and distance via Equation (8). However, even if π E could not be measured, the combination of (t E , θ E ) [so, also, µ rel ] provided more powerful constraints on the Bayesian mass and distances estimates using Galactic-model priors than is possible from the t E constraint alone. Moreover, the determination of µ rel allows one to accurately estimate how long one must wait in order to separately resolve the lens and source in high-resolution follow-up observations using, e.g., AO on large telescopes or telescopes in space (e.g., Batista et al. 2015;Bennett et al. 2020Bennett et al. , 2015. For this reason, virtually all papers on planetary microlensing events make a serious effort to measure θ * . We follow this general practice here, but we note in advance that, with the exception of two events (OGLE-2018-BLG-1647 and OGLE-2018-BLG-0932), the value of doing so is likely to be minimal. This is because, for all of the other events analyzed here, there are only weak upper limits on ρ, or in some cases no limits at all.
The limit on ρ can be characterized as "weak" if it leads to a "weak" lower limit on the proper motion µ lim = θ * /t * ,lim , where t * ≡ ρt E and t * ,lim ≡ ρ lim t E . In turn, µ lim is "weak" if it does not exclude a significant fraction of the parameter space.
We quantify this as follows. Following the Appendix of Gould et al. (2021), we note that for events with bulge lenses and bulge sources, the fraction of events with µ rel < µ lim ≪ σ is where we have modeled the bulge proper-motion distribution as an isotropic Gaussian with dispersion σ = 2.9 mas yr −1 . One can show that in this low µ lim regime, the probability for disk-bulge lensing is even lower. Thus, for example, if µ lim < ∼ 0.5 mas yr −1 (as in most of our events) then fewer than p < ∼ 10 −3 of simulated events will be eliminated by imposing this limit, implying negligible impact on the Bayesian estimate.
Nevertheless, while θ * is itself of little use in these cases, the measurements of the source color and magnitude, which are needed to determine θ * , can be important for the interpretation of future AO observations. Together, they will enable prediction of the source flux in the observed band (e.g., H or J), and so allow one to determine which of the two stars is the source, with the other being the lens, whose properties will be the main subject of interest. These observations will, by themselves, yield µ rel (from the observed separation and elapsed time), and so θ E = µ rel t E . Together with the lens flux, this will enable good estimates of M and D L . Thus, even though these θ * measurements are likely to be of little use, either now or in the future, they are a small additional step relative to the actually necessary color and magnitude measurements. Hence, we report them as well.
Our general approach (with a few exceptions that are explicitly noted) will be to obtain pyDIA (Albrow 2017) reductions of KMT data at one (or possibly several) observatory/field combinations. These yield the microlensing light curve and field-star photometry on the same system. We then determine the source color by regression of the V-band light curve on the I-band light curve, and the source magnitude by regression of the I-band light curve on the best model. We then transform the instrumental KMT photometry to calibrated OGLE photometry, usually OGLE-III (Szymański et al. 2011), but in two cases, OGLE-II (Szymański 2005;Kubiak & Szymański 1997;Udalski et al. 2002). If there is inadequate V-band signal in a single observatory/field, we repeat the procedure for several, check for consistency, and then combine them. In two cases, we are not able to measure (V − I) from the light curve. In one of these cases, we infer the color by combining OGLE-IV I-band observations with H-band observations from the UKIRT microlensing project (Shvartzvald et al. 2017). In the other, we make use of a deep, high-resolution color-magnitude diagram (CMD) based on archival Hubble Space Telescope (HST) data (Holtzman et al. 1998). Figures 12 and 13 show the resulting CMD for each event, with the position of the source and the centroid of the red clump indicated in blue and red respectively. Table 12 lists these values   Figure 12 for the remaining 4 (out of 10) of the events analyzed in this paper. and also shows the steps leading to the calculation of θ * for each event.
For this, we follow the method of Yoo et al. (2004). We adopt the intrinsic color of the clump (V − I) 0,cl = 1.06 from Bensby et al. (2013) and its intrinsic magnitude from Table 1  , I] cl . We convert from V/I to V/K using the VIK color-color relations of Bessell & Brett (1988) and then derive θ * from the color/surface-brightness relations of Kervella et al. (2004). After propagating errors, we add 5% in quadrature to account for errors induced by the overall method.
Where relevant, we report the offset of source from the baseline object. In all cases, this is found by comparing the difference image near peak to the baseline object position in the template.
Comments on individual events follow.

OGLE-2018-BLG-1126
The CMD is shown in Figure 12. There are no useful constraints on ρ. We note that the baseline object has [(V − I), I] base = (2.14, 18.70), implying that the blend has [(V − I), I] B = (2.13, 18.98), i.e., similar in color but about 9 times brighter than the source. We find that it is displaced from the event by 260 mas, meaning that it is almost certainly unrelated to the event. Most likely, it is a bulge sub-giant. Its brightness and proximity prevent any useful constraints on the lens flux. On the positive side, it is unlikely to interfere with future AO observations.

KMT-2018-BLG-2004
The CMD is shown in Figure 13. The constraints on ρ have practically no impact. The baseline object (I base = 18.88) is offset from the source by about 600 mas, meaning that the blend has I B ≃ 19.8 and is almost certainly unrelated to the event. Moreover, the blend color is very poorly determined. Hence, we do not display it in the CMD. We adopt I L > 19.6, which corresponds to I L,0 > 18.1 for bulge lenses (and other lenses that are behind essentially all the dust). This will have a minor effect. See Section 5.2. The magnitude listed in Table 12 is for the planetary solution with the lower χ 2 , as will always be the case except when otherwise specified. In this case, the other solution would have a larger θ * by 1.4%, i.e., a small difference compared to the error bars.
This event is not in the OGLE-III footprint, but fortunately it is in the OGLE-II footprint (Szymański 2005;Kubiak & Szymański 1997;Udalski et al. 2002). As indicated in Figure 13, we therefore calibrate the photometry using OGLE-II.

OGLE-2018-BLG-1647
The CMD is shown in Figure 12. In this case, there are ρ measurements for both solutions. Because the wide solution is favored by ∆χ 2 = 17, we do not further consider the close solution. While the fractional error in ρ is fairly large (20%), we note that very low values are strongly excluded. For example, ρ > 0.0023 at 2.5 σ, which is very similar to the naive extrapolation from the 1 σ error bar. This corresponds to θ E < 0.20 mas and µ rel < 1.4 mas yr −1 at the same significance. Hence, this is likely to be a low-mass lens in the bulge.
OGLE-III photometry, which resolves out a nearby neighbor at about 600 mas thereby showing a baseline magnitude I base = 19.96, implies an estimated blend magnitude I B = 20.57. We set a more conservative limit on the lens brightness I L > 20.30. Given the extinction toward this field, A I = 1.43, this corresponds to I L,0 > 18.87 for lenses that that are behind essentially all the dust. Hence, given that the θ E measurement already favors a low-mass bulge host, the flux constraint plays a limited A&A proofs: manuscript no. ms role. Because we do not have a color determination for the baseline object (hence, also for the blend), we do not display it on the CMD.

OGLE-2018-BLG-1367
The CMD is shown in Figure 12. Again, the limit on ρ is very weak, corresponding to θ E > 0.048 mas and µ rel > 0.77 mas yr −1 , which are hardly constraining.
OGLE-III shows a baseline magnitude I base = 18.57, leaving an estimated blend magnitude I B = 19.92. We set a more conservative limit on the lens brightness I L > 19.70, which corresponds to I L,0 > 18.79 for lenses behind essentially all the dust. This is a very similar, mildly constraining limit as in the case of OGLE-2018-BLG-1647. Again, we do not display the blend on the CMD due to poor color determination.

OGLE-2018-BLG-1544
The CMD is shown in Figure 12. The source is blended with a clump giant [(V − I), I] base = (2.88, 16.74), which is separated by 600 mas. Hence, the blended light cannot be constrained. Following the logic that was applied to OGLE-2018-BLG-1647, the limit, ρ < 0.012, implies µ rel > 0.45 mas yr −1 , which is not useful.

OGLE-2018-BLG-0932
The CMD is shown in Figure 13. We are not able to accurately measure the V-band source flux in spite of the source being in or near the clump, for two reasons: the source is heavily reddened and the peak magnification is low (A max = 1.47). Fortunately, the event lies in the UKIRT microlensing footprint (Shvartzvald et al. 2017), which allows us to determine the source color on an [(I − H), I] CMD. To this end, we match OGLE-IV I and UKIRT H data, which are shown in Figure 13. We find that the source is ∆(I − H) = −0.016 ± 0.054 bluer than the clump, from which we infer that it is ∆(V − I) = −0.01±0.03, which is the basis of our color determination in Table 12.
Note that for this field, I OGLE−III − I OGLE−IV = 0.04. We do not correct for this offset from standard magnitudes in Table 12. This makes no difference for our estimate of θ * , which depends only on relative photometry. However, it should be noted in the unlikely event that there is future, high-precision, I-band photometry that could probe this level of difference.
Of the 10 events analyzed in this paper, OGLE-2018-BLG-0932 is the only one with a precise ρ measurement and one of only two with any ρ measurement. Combining this with our determination of θ * , we find, θ E = 0.458 ± 0.033 mas µ rel = 6.22 ± 0.44 mas yr −1 .
As discussed in Section 3.7, the blending is consistent with zero, but is not well measured. we have set f B = 0 in the fit, but (given that the source is a clump giant), we cannot set any useful limits on the lens flux.

OGLE-2018-BLG-1212
Before evaluating the CMD information for this event, it is important to recall that there is a very precise, and fairly large, parallax measurement π E = 0.767 ± 0.019. As discussed in Section 3.8, this result strongly favors (but does not prove) that the lens is relatively nearby, i.e., only a few kpc from the Sun. In light of this, is is notable that the Gaia measurement of the "baseline object", π base =3.41 ± 0.93 µ base (N, E) =(−4.73, −7.32) ± (0.98, 1.16) mas yr −1 , suggests that the baseline object may be a very nearby star, or possibly a blend of a nearby object with a more distant star. In particular, one scenario is that this "object" is comprised of a bulge source and a very nearby disk lens (or a companion to the lens). If, for example, the lens contributed half the light (and if the Gaia measurement were not itself corrupted, see below), then the lens should have π L ∼ 6 mas. In this case, π rel ≃ π L , in which case the relative proper motion would be µ rel = π rel /π E t E ∼ (57 mas yr −1 )(π rel /6 mas).
Such a high lens-source relative proper motion would have two consequences that are not confirmed. First, the Gaia proper motion itself would be fractionally affected at the same level as the parallax (in this example, by 50%), whereas the actual Gaia proper motion is just ∼ 2 σ from the mean of the bulge distribution. Second, such a high-motion star would be separately resolved in OGLE-II images (from 1999) and would be recognizable either as a "new star" (at position angle φ ∼ 224 • ) compared to the OGLE-IV finding chart (from 2010) or as being displaced from the corresponding OGLE-IV object in the same direction. We find no such high proper motion stars in the OGLE-II images.
Thus, while the large Gaia parallax may be suggestive of the presence of a nearby star in the Gaia aperture (whether related to the event of not), it is difficult to infer anything about the lens from this measurement. In addition, we note that Gaia reports a RUWE value of 1.89, probably indicating some form of contamination of the measurement.
Interestingly, OGLE-2018-BLG-1212 was the subject of a Gaia alert 6 , on 2018-10-05 08:38:24, as being a transient of unknown origin. Two of the Gaia points, at HJD ′ = 8396.86 and 8396.93, were just 3 days after the anomaly. However, there are only four significantly magnified Gaia points in total. Hence, these data do not help constrain the event.
The CMD is shown in Figure (1.59, 19.16). From its position on the CMD, the blend light could very well be dominated by a companion to the source.
Of more direct interest, the blend light cannot be dominated by the lens. For example, given the parallax measurement π E ≃ 0.767, an M = 0.25 M ⊙ lens would lie at D L ∼ 0.76 kpc, and so would have roughly I L ∼ 19, thus approximately accounting for the I B light. However, after accounting for E(V − I) L ∼ 0.4 of reddening, it would have (V − I) L ∼ 3.4, implying V ∼ 22.4, which is almost 2 magnitudes redder than the blend. On the other hand, if the lens were at D L ∼ 1.5 kpc (as crudely estimated in Section 3.8 based on kinematic arguments), then M ∼ 0.11 M ⊙ . In this case, the lens would not contribute significantly to the blended light, thereby avoiding all photometric constraints. In principle, the lens could be farther and so have yet lower mass, but these distances are disfavored by both the declining mass function and the kinematic arguments. These will automatically be taken into account when we carry out a Bayesian analysis in Section 5.7.
Thus, in spite of the several intriguing facts about the blend, in the end, its only implication for the analysis is that it places an upper limit on the lens light, for which we adopt I L > 19.0. However, as we discuss in Section 5.7, even this role has a relatively modest practical effect.
Finally we note that the proper motion can be expressed µ rel = θ E /t E = κMπ E /t E , implying, µ rel,hel = 45(M/M ⊙ )mas yr −1 . Hence, if the lens is luminous (M > ∼ 0.075 M ⊙ ), then its proper motion is > ∼ 3.3 mas yr −1 . Therefore, it will be separated from the source by at least 40 mas by 2030, a plausible first light for AO on ELTs. Note that even if the lens were a white dwarf (WD), it would almost certainly be visible in AO follow-up. For example, at M = 0.6 M ⊙ , a relatively dim WD with M K = 14, would be at D L ∼ 0.33 kpc and so K ∼ 21.6, which would be visible in ELT observations. In this case, the proper motion would be µ rel,hel = 27 mas yr −1 , so that the separation in 2030 would be ∼ 300 mas. Hence, a second epoch would be required for confirmation. Nevertheless, this does mean that a non-detection in ELT AO follow-up would imply that the host is a brown dwarf.

KMT-2018-BLG-2718
The CMD is shown in Figure 13. Due to the small variation in the V-band light curve, our standard procedure for determining the source color yields a very imprecise result: (V − I) 0,S = 1.54 ± 0.33. We therefore estimate the color from the I-band offset between the source and the clump, which yields (V − I) 0,S = 1.37±0.14, using the Galactic bulge CMD derived from HST observation by Holtzman et al. (1998). (As usual, all aspects of this evaluation are based on the lowest-χ 2 solution, i.e., the planetary solution with s > 1.) For the four solutions, the limits on ρ shown in Table 9 correspond to t * = (1.6, 1.1, 2.4, 3.1) days. For the second of these, i.e., the best fit, this corresponds to µ rel > 0.12 mas yr −1 . The excluded region contains a fraction p < (µ rel,lim /2.9 mas yr −1 ) 3 /6 √ π = 7 × 10 −6 . That is, this limit is completely unconstraining. For the other three cases, the limit is even weaker.
It is unlikely that the ambiguity between planetary and binary solutions can be decisively resolved until RV observations become feasible for this very faint host. Because the planetary solution is formally favored by ∆χ 2 = 12.7, the event can plausibly be included in mass-ratio function studies. However, this will require a specific decision.

KMT-2018-BLG-2164
The CMD is shown in Figure 12. There are no useful constraints on ρ. The OGLE-III baseline object has I base = 20.54, yielding I B = 20.70. We adopt I L > 20.40, corresponding to I L,0 > 18.85 for lenses lying behind essentially all the dust, which is mildly constraining. We remind the reader that there is a factor ∼ 200 ambiguity in q for the two classes of solutions that we presented in Section 3.10, which cannot be resolved except by RV observations in the far future. Hence, we believe that this event is unlikely to attract interest for AO follow-up observations.

OGLE-2018-BLG-1554
The CMD is shown in Figure 13. As we discussed in Section 3.11, there is a ρ measurement only for the planetary solution. We argued that its high value, ρ ≃ 0.03, rendered the planetary solution highly unlikely.
The OGLE-III baseline object has [(V − I), I] base = (2.01, 18.99), which is very similar to the source values from Ta-ble 12, [(V − I), I] S = (2.02, 19.10), implying that the source is almost unblended. We adopt an upper limit on lens light I L > 22.80, corresponding to I L,0 > 21.17 for lenses lying behind most or all of the dust. This would be a significant constraint. However, because the event is not clearly planetary, this constraint has no practical impact. See Section 5.10.

Physical Parameters
None of the 10 events reported in this paper have both θ E and π E measurements. Hence, as is customary for a substantial majority of microlensing planets, we make Bayesian estimates of the physical parameters of the system by incorporating priors from a Galactic model. In the subsections below, we summarize the constraints that are derived from the light-curve analysis and CMD analysis, as reported in Sections 3 and 4. Our general approach is to simulate events based on a Galactic model and then assign each event a weight (possibly zero) depending on how well it matches these constraints. For example, if (as is true of several events), the only constraint is the measurement of the Einstein timescale t E ± σ(t E ), then the weight of the sim- The Galactic model is summarized in Section 5 of Han et al. (2021b).
In Table 13, we present the resulting Bayesian estimates of the host mass M host , the planet mass M planet , the distance to the lens system D L , and the planet-host projected separation a ⊥ . For the majority of events, there are two or more competing solutions. For these cases we show the results of the Bayesian analysis for each solution separately, and we then show the "adopted" values below these. For M host , M planet , and D L , these are simply the weighted averages of the separate solutions, where the weights are the product of the two factors at the right side of each row. The first factor is simply the total weight from the Bayesian analysis. The second is exp(−∆χ 2 /2) where ∆χ 2 is the χ 2 difference relative to the best solution. See Ryu et al. (2022). For a ⊥ , we follow a similar approach provided that either the individual solutions are strongly overlapping or that one solution is strongly dominant. However, if neither condition is met, we enter "bi-modal" instead.
We present Bayesian analyses for 8 of the 10 events, but not for KMT-2018-BLG-2164 and OGLE-2018-BLG-1554. See Sections 5.9 and 5.10. Figures 14 and 15 show histograms for M host and D L for these 8 events.

OGLE-2018-BLG-1126
The only constraint is the measurement of t E . As a result the histograms of host mass and distance are extremely broad. See Figures 14 and 15. The planet has a similarly broad distribution, but is generally in the Neptune-class range. We recall from Section 3.2 that the planet is detected by only ∆χ 2 = 69.

KMT-2018-BLG-2004
This event has three constraints in addition to the t E measurement. First, there is the 1-D parallax measurement, π E, = 0 ± σ(π E, ), where the error bar and orientation ψ of the π E, measurement take on 4 pairs of values that depend on the signs of u 0 and π E,N , as given just above Equation (11). In addition, there are limits on ρ (< 0.021 or < 0.024) and on lens light, I L > 19.60. The 1-D parallax measurement is incorporated via Equations (10) and (11) as described in Section 3.3. The ρ con-A&A proofs: manuscript no. ms −0.31 and hence it plays virtually no role. The main information comes from the π E, measurement. Because this measurement is consistent with π E ∼ 0, bulge lenses are permitted. Of course, the contours extend up into the north-east quadrant of the π E diagram, which is preferred by disk lenses, so these are also permitted. However, because the parallax constraint has constant width, it is more restrictive of disk lenses (which have higher π E ) than bulge lenses. Hence, bulge lenses, which are already favored by higher phase-space density, receive a further boost. Within this context the flux constraint plays a modest secondary role by eliminating some bulge lenses at the very top of the main sequence. The planet has a Saturn-class mass, and the system is very likely in, or at least close to, the bulge.

OGLE-2018-BLG-1647
The wide solution is favored by ∆χ 2 = 17, so we consider the close/wide degeneracy to be resolved, and so we only show one solution in Table 13. Both t E and ρ are measured from the light curve, and so t E and θ E = θ * /ρ enter as constraints (Tables 4 and  12). For the latter we adopt θ E = 91 ± 18 µas. Although the error in this measurement is large, θ E is nevertheless constrained to be much smaller than in typical events, which strongly favors a low mass M host ∼ 0.1 M ⊙ host in or near the Galactic bulge. See Figures 14 and 15. Hence, despite its high mass ratio, q ≃ 10 −2 , the planet is likely to be of Jovian mass. We also incorporate the limit on lens light, I L > 20.30 from Section 4.3. However, this plays only a small role because the θ E measurement already heavily disfavors lenses that are this bright.

OGLE-2018-BLG-1367
Like KMT-2018-BLG-2004, this event has three constraints in addition to the t E measurement. There is 1-D parallax measurement, π E, = 0.165 ± 0.040, as well as limits on ρ < 0.016 and on lens light I L > 19.70. The 1-D parallax measurement is incorporated via Equations (10) and (11) with ψ = 87.30 • , as described in Section 3.5. The ρ constraint implies θ E > 48 µas, corresponding to µ rel > 0.8 mas yr −1 , and hence it plays almost no role. The main information comes from the π E, measurement. First, it implies π E ≥ π E, ≃ 0.165, so if the lens is in the bulge (π rel < ∼ 10 µas), then M = π rel /κπ 2 E < ∼ 0.1 M ⊙ , which greatly reduces the phase space accessible to bugle lenses. Second, the smallest values of π E are in the north-east quadrant of the π E diagram, which is the preferred location of disk lenses. Hence, the lens distance distribution broadly peaks in the disk at D L ∼ 5 kpc (i.e., π rel ∼ 120 µas) and so at masses M = π rel /κπ 2 E < ∼ 0.5 M ⊙ .  Table 13.
Where there are several solutions, we show the distribution for the one with the lowest χ 2 . However, as can be assessed from the Table 13, the other solutions hardly differ.
The flux constraint therefore plays a relatively minor role because lenses that would violate it are already heavily disfavored. The planet is again Jovian class.

OGLE-2018-BLG-1544
Nominally, this event has two constraints, a t E measurement and an upper limit on ρ. However, the latter leads to a very weak proper-motion constraint µ rel > ∼ 0.4 mas yr −1 , which therefore plays virtually no role. As with OGLE-2018-BLG-1126 (which has only a t E measurement), the posterior Bayesian distributions of mass and distance are extremely broad. However, because t E is smaller in the present case by a factor ∼ 0.65, these distributions are shifted to somewhat lower mass and distances. See Figures 14 and 15. Because of the event's high mass ratio, q > ∼ 0.01, the planet mass estimate is centered near the planet-BD boundary, but with a wide dispersion.

OGLE-2018-BLG-0932
In addition to the t E measurement, this event has two constraints, a measurement of ρ (leading to measurements of θ E = 0.458 ± 0033 mas and µ rel = 6.22 ± 0.44 mas yr −1 ), and a Gaia measurement of the source proper motion µ S (N, E) = (−7.53, −8.81) ± (0.17, 0.26) mas yr −1 . There are also Spitzer microlensing data for this event, which should ultimately yield a π E measurement. However, the analysis of these data is beyond the scope of the present work and will be presented elsewhere.  Table 13.
Where there are several solutions, we show the distribution for the one with the lowest χ 2 . However, as can be assessed from the Table 13, the other solutions hardly differ.
We note that in Galactic coordinates, the source proper motion is µ S (l, b) = (−10.96, +3.78) masyr −1 , which is ∼ 6.2 mas yr −1 from the bulge mean, i.e., slightly more than 2 σ and tending in the direction of anti-rotation. This means that a bulge lens would be expected to generate µ rel ∼ 7 mas yr −1 (quite consistent with what is observed), while disk lenses would be expected to generate µ rel ∼ 11 mas yr −1 . Thus, the Gaia measurement increases the likelihood of bulge lenses, which are already strongly favored by phase-space considerations. The net result can be judged from Figures 14 and 15. The planet is in inferred to have Jovian mass.
It will be interesting to compare the host mass estimate in Table 13 to the results of from the future Spitzer analysis. Roughly speaking M ≃ 0.72 ± 0.27 M ⊙ corresponds to π E = θ E /κM = 0.078 ± 0.029.

OGLE-2018-BLG-1212
For this event, there are two constraints in addition to the t E measurement. First, there is a well very well-localized parallax measurement, π E = 0.767 ± 0.019, whose direction (in the LSR frame) is closely aligned to Galactic rotation, i.e., 2 • ± 1 • . Second, there is a limit on lens flux, I L > 19.0. Table 13 shows that the Bayesian mass and distance estimates are in good accord with the expectations outlined in Section 3.8, which was based purely on kinematic arguments. For comparison, we conducted a separate Bayesian analysis in which the flux constraint was ignored. This test showed that inclusion A&A proofs: manuscript no. ms of the flux constraint drove the distance down from 1.62 kpc to 1.55 kpc, i.e., a small effect.
We also show in Table 13 the results from the alternate jerk-parallax solutions. These have almost no formal statistical weight, as indicated by the last two columns. The main point is to show that the principal implications for the host and planet masses, the system distance, and the projected separation are not qualitatively different. In particular, the planet is a Saturn-class object that is at 1-2 kpc, orbiting a late M dwarf. As discussed in Section 4.7, these predictions can be tested at first AO light on ELTs.

KMT-2018-BLG-2718
Beyond the t E measurement that is common to all events, there is only a weak constraint on the normalized source size, ρ < 0.0068, which leads to an exceedingly weak limit on the proper motion, µ rel > ∼ 0.1 mas yr −1 . Thus, the only real information from the photometric light curve is that the Einstein timescale is exceptionally long. Lenses essentially anywhere along the line of sight can generate such long timescale events by virtue of the rare chance that the source and lens proper motions are very similar. At any distance, large masses M ∝ t 2 E are favored, and these general remarks are well reflected in the distributions shown in Figures 14 and 15.

KMT-2018-BLG-2164
KMT-2018-BLG-2164 is neither unambiguously planetary in nature nor is the planetary interpretation significantly preferred. That is, it has only ∆χ 2 = 4.7 relative to the binary interpretation. Even if Gaussian statistics applied, the binary probability would be ∼ 10%. Therefore, it should not be "registered as a planet" in community databases, and we therefore refrain from trying to characterize it using Bayesian estimates. It is included in the present study only for completeness, i.e., to identify all events with viable planetary solutions, regardless of whether these are unique. .

OGLE-2018-BLG-1554
The case for a planetary interpretation for OGLE-2018-BLG-1554 is even weaker than for KMT-2018-BLG-2164. First, the 1L2S solution is slightly preferred by χ 2 . Second, there is a competing binary solution at ∆χ 2 ≃ 0. Third, as we remarked in Section 4.10, the measured θ E and µ rel for the planetary (but not binary or 1L2S) solution are highly unlikely a priori. Again, this event is only included in this study for completeness. We again counsel against its "registration" as a planet in community databases, and so we refrain from a Bayesian characterization. .

Conclusions
The goal of this paper was to complete the analysis of all events from 2018 with viable planetary solutions that were identified by the KMTNet AnomalyFinder system and that lie in one or more of the 6 KMT prime fields. Because the main motivation was to prepare a complete sample for statistical analysis, we pushed the boundaries of this sample beyond what will ultimately be used in such studies, and we provide sufficient information to permit future workers to set their own detailed boundaries. In Notes. Event names are abbreviations for, e.g., OGLE-2018-BLG-1185 and KMT-2018-BLG-1024. a: s degeneracy. b: Nearly factor 2 q degeneracy. c: 1L2S/2L1S degeneracy. d: Two-planet system. e: planet in binary system. f: planet/binary degeneracy. g: large q degeneracy. particular, we report on all events with viable solutions with mass ratios q < 0.06, and we provide detailed analysis of all events that have viable solutions with q < 0.03, even for cases that would not normally be published due to ambiguity with binary-lens (q > 0.03) and/or binary-source (1L2S) solutions. Indeed, two of the 10 events that we have analyzed are in one or both of the last two categories and would not normally be pub- There is one additional AlertFinder recovery, OGLE-2018-BLG-0100, that remains in preparation, but this has an ambiguous mass ratio q at the factor 100 level. Table 14 shows the 26 events with viable planetary solutions that were recovered or discovered by AnomalyFinder from the 2018 KMT prime-field events. The 4 previously published discoveries are from Hwang et al. (2022) and Wang et al. (2022). References are given for the 11 previously published recoveries. Note that among these, OGLE-2018-BLG-1700 is marked as a planet in a binary system because the statistical properties of AnomalyFinder discoveries/recoveries may differ for such systems. The 10 entries marked "This Work" include 7 discoveries and 3 recoveries, while one previously known planetary solution remains "in preparation". We consider that the 3 entries below the double line are unlikely to enter a mass-ratio function analysis, while 5 others (OGLE-2018-BLG-1126, KMT-2018-BLG-1025, KMT-2018-BLG-2004, OGLE-2018-BLG-1700, and OGLE-2018-BLG-1544 will require detailed assessment. Here, we provide only the information necessary for these assessments but not the assessments themselves. Through the course of our systematic study of the 10 events published here, we noticed that the "s † " formalism that was introduced by Hwang et al. (2022) for heuristic analysis should be slightly modified, from using the arithmetic to the geometric mean of the two solutions. In this form, it unifies the so-called close/wide degeneracy of Griest & Safizadeh (1998) for central and resonant caustics with the so-called inner/outer degeneracy of Gaudi & Gould (1997) for planetary caustics, a unification that was previously conjectured by Yee et al. (2021).