Search for Very-High-Energy (E $>$ 100 GeV) Emission from Geminga Supernova by VERITAS

A. U. Abeysekara (for the VERITAS collaboration)

REVIEW 3 major objections 4 minor 18 references

Search for Very-High-Energy (E $>$ 100 GeV) Emission from Geminga Supernova by VERITAS

T0 review · 3 major / 4 minor · reviewed 2026-08-14 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read This paper claims that the Matched Runs Method (MRM) lets VERITAS estimate the gamma-ray background over the entire field of view from archival point-source observations, validated on source-free dwarf spheroidal fields.

desk verdict A plausible, honestly framed methods-in-progress paper: the Matched Runs Method is a sensible automated reuse of archival VERITAS runs for extended sources, but the abstract claim that it is 'demonstrated effective' outruns null tests in two empty fields. read the letter →

arxiv 1908.05369 v1 pith:HHW2BJXW submitted 2019-08-14 astro-ph.HE

A. U. Abeysekara (for the VERITAS Collaboration) This is my paper · ORCID

classification astro-ph.HE

keywords very-high-energygammaraysGemingapulsarextendedsourceanalysisbackgroundestimationimagingatmosphericCherenkovtelescopecosmic-raypositronexcessMatchedRunsMethodVERITAS

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

VERITAS can detect gamma rays from 100 GeV to above 30 TeV, but its standard background estimators need empty patches in the field of view, which fails when a source fills the field. This paper argues that the Matched Runs Method (MRM) cures that by searching archived observations of unrelated point sources for runs whose elevation, azimuth, date, and event rate match the target field, then using those runs to model the background. If MRM works, VERITAS can measure the extended TeV halo around Geminga, filling the energy gap between Fermi-LAT and HAWC and testing whether Geminga still explains the cosmic-ray positron excess. The paper's validation is statistical: MRM maps of the dwarf spheroidals Segue 1 and Ursa Minor show significance distributions consistent with pure noise.

What carries the argument

The central object is the Matched Runs Method (MRM), a background-estimation algorithm that reuses archival observations of point-like sources as OFF runs for an extended-source ON field. Its matching function minimizes $\Delta N_{\rm CR}=\sum_{\rm runs}|N_{\rm CR}^{\rm ON}-N_{\rm CR}^{\rm OFF}|$, the summed difference in non-gamma-like event counts between candidate ON and OFF runs, subject to cuts $\Delta {\rm Az}\le 10^\circ$, $\Delta {\rm El}\le 10^\circ$, and $\Delta T\le 360$ days. The method then forms the background from ring or wobble-offset regions within those matched OFF fields, so no dedicated blank-field observation is required. The machinery converts the problem of estimating background under an extended source into a database search for geometrically and temporally similar runs.

What would settle it

Run MRM on simulated data with an injected extended source matching the HAWC Geminga halo (roughly 2-degree radius with known flux and morphology) and compare the reconstructed sky map to the input; if the significance or surface brightness comes out biased, or the source is smeared or suppressed, the central claim fails. Alternatively, if the MRM analysis of the 93-hour Geminga dataset does not recover the HAWC-detected excess, the method's validity as an estimator for real extended sources is contradicted.

Watch

Extended reading notes

Core claim

The central claim is that a gamma-ray source with spatial extent comparable to the camera field of view can be analyzed with VERITAS by borrowing background from archival runs on point-like sources, instead of taking separate blank-field OFF runs. The Matched Runs Method (MRM) automates the selection of matched runs using four criteria—elevation difference, azimuth difference, time gap, and difference in non-gamma-like event counts—and then uses ring or wobble-offset regions in those OFF observations to model the background across the whole field. Applied to 26 hours of Segue 1 and 5 hours of Ursa Minor, MRM produces sky maps with no excess and significance distributions consistent with a standard normal ($\mu=0$, $\sigma=1$), which the paper takes as evidence that the background is correctly estimated everywhere in the field of view. The intended application is the 93-hour Geminga dataset, where MRM would fill the 500 GeV to 1 TeV gap between Fermi-LAT and HAWC and reduce the surface-brightness uncertainties that currently hamper diffusion-model tests.

Load-bearing premise

The method is validated only on fields known to contain no gamma-ray source, where a null significance map shows the background estimate is not grossly biased; the paper does not yet show that MRM correctly reproduces a real, extended source's signal.

Editorial extensions

If this is right

If MRM is valid, VERITAS can produce background-modeled sky maps of sources as large as the field of view using archival data, without doubling the observation time.
The method gives a path to measure the Geminga TeV halo across 100 GeV to 30 TeV, bridging the energy gap between Fermi-LAT and HAWC.
MRM enables deep observations of faint extended sources that would be impractical with the ON/OFF technique because of its factor-of-two exposure penalty.
Because matched runs can come from any point-like calibration source such as the Crab or 1ES 0229+200, archival VERITAS data becomes reusable as background for many extended-source analyses.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

The null-map validation on source-free dwarf spheroidals does not by itself prove that MRM preserves the signal of a real extended source; a natural next test is injecting a simulated Geminga-like halo into a source-free field and checking that the reconstructed flux and morphology match.
If MRM succeeds on Geminga, the same procedure could be applied to other HAWC-detected pulsar halos mentioned in the paper, such as Monogem, to map their energy-dependent morphology at VERITAS's finer angular resolution.
The paper's rationale connects to the debate over radially varying positron diffusion: with sub-degree surface-brightness profiles from MRM, one could test whether the diffusion constant inferred near the pulsar differs from the value outside the halo, which the paper only notes as an unresolved question.
A residual systematic risk is that matched runs from different sky positions could imprint camera-response gradients; the Gaussian significance test would catch large-scale biases but not necessarily small distortions in a faint extended source.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

MRM is a practical extension of the ON/OFF idea: instead of dedicating new OFF runs, it automatically matches archival runs to the ON pointing using elevation, azimuth, time gap, and a non-gamma-ray rate difference. That is genuinely useful for VERITAS, where the standard ring and reflected-region methods lose power when a source fills the field of view. The automated matching objective is new, and the two null tests on Segue 1 and Ursa Minor are real evidence that the method does not obviously bias empty fields. The maps look featureless, the significance distributions are reported as consistent with a standard normal, and no part of the result is circular: the validation is not fit to the target signal. The prose is clear about the intended application and honest about the fact that Geminga analysis is still pending.

The soft spot is the gap between what is tested and what is claimed. The abstract says MRM 'has been demonstrated to be an effective technique,' and the conclusion says it is 'capable of estimating the background of the entire field-of-view.' What the paper actually shows is that two source-free dwarf-spheroidal fields analyzed with MRM look consistent with pure background. That does not test the load-bearing case: a bright extended source like Geminga filling much of the field, where background subtraction and exposure normalization could partially absorb the signal. It also does not test whether MRM recovers the correct flux or morphology of a known extended source. The global Gaussian shape is a weak probe, since a localized few-sigma excess in a multi-thousand-bin map would barely move the histogram. And the paper gives no quantitative fit statistics for the null distributions, just an eyeball 'consistent.' The thresholds in the matching algorithm are hand-set, which is fine as a starting point, but there is no study of how sensitive the results are to them.

These are real limitations, but they are not fatal for what this paper actually is. It is an ICRC proceedings contribution, not a refereed method paper, and its own Section 5 says the procedure will be applied to Geminga only 'when the validation tests are completed.' The abstract simply gets ahead of the evidence. The right fix is to temper the wording and, for the next version, add a known extended source with a mild extension, like IC 443, or an injected-source Monte Carlo, to show that the method preserves signal while subtracting background.

Who should read this? Anyone working on IACT analysis of extended sources, especially within VERITAS, and anyone trying to bridge the Fermi-LAT/HAWC energy gap. It is a credible methodological contribution, not a detection. I would cite it as a reference for the method, and I would send a revised version to a serious referee, because the method is plausible and the energy gap it targets is scientifically important. But the 'demonstrated effective' language should be revised before publication.

Referee Report

3 major / 4 minor

Summary. This paper presents the Matched Runs Method (MRM), a background estimation technique for VERITAS observations of gamma-ray sources extended on the scale of the camera field of view. The authors describe the method, in which OFF runs are selected from archival data using matching criteria in elevation, azimuth, time, and non-gamma-ray event rate, and present validation tests on two dwarf spheroidal galaxies (Segue 1 and Ursa Minor) showing that significance maps are consistent with a standard normal distribution. The paper also summarizes 93 hours of VERITAS observations of Geminga but states that the MRM analysis awaits completion of validation tests.

Significance. If fully validated, MRM would be a useful addition to IACT background estimation, potentially enabling searches for sources extended on FOV scales without dedicated blank-field observations. The paper's strength is that it describes a concrete, automatable algorithm and provides an empirical null test on archival data. However, the current evidence only demonstrates that MRM produces globally background-like maps in source-free fields; it does not yet establish that the method preserves the signal of an extended source, which is the intended application. The claims in the abstract and conclusion go beyond what the presented tests support.

major comments (3)

[Section 4 and Section 5] Section 4 validates MRM only on source-free fields (Segue 1 and Ursa Minor), and Section 5 concludes that the method is 'capable of estimating the background of the entire field-of-view.' This conclusion does not follow from the evidence because a blank-field null test cannot rule out that MRM partially subtracts or distorts a real extended source such as the ~2° Geminga halo. The paper should include a validation with a known or simulated extended source, demonstrating that the source flux and morphology are recovered, before claiming that MRM is demonstrated for the Geminga search.
[Section 4, Figures 1b/2b] The global significance distribution check (μ=0, σ=1) is insensitive to localized, spatially correlated residuals. A single 4–5σ excess in a multi-thousand-bin significance map would barely change the global histogram, so this statistic does not validate the method for source searching. The authors should report a quantitative goodness-of-fit test (e.g., a Kolmogorov-Smirnov or Anderson-Darling test with the p-value, or a χ² over the histogram) and an analysis of the distribution of the maximum significance or spatial autocorrelation of residuals.
[Abstract and Section 5] The abstract states that MRM 'has been demonstrated to be an effective technique by applying it to archival VERITAS data,' but Section 5 states that 'The MRM procedure will be applied to the Geminga data set when the validation tests are completed.' This internal inconsistency overstates the current status. The claims should be revised to match the evidence: at present, the method has passed null tests, and extended-source validation is still in progress.

minor comments (4)

[Abstract] There is a typo in the abstract: 'spin-down time sca le' should be 'spin-down time scale.'
[Section 3] The definition of ΔNCR in equation (3.1) is unclear: it is written as a sum over 'All Runs,' but the text says the algorithm 'minimizes the ΔNCR to a list of ON and OFF runs.' Please clarify whether this is a sum over pairs of matched ON/OFF runs or a global matching criterion.
[Title] The paper title promises a 'Search for Very-High-Energy Emission from Geminga,' but the paper contains no Geminga results; consider changing the title to reflect that this is a methods paper, or adding a short statement that the search will be presented separately.
[References] Reference [18] cites 'ICRC proceeding 2007' but the arXiv number 1708.07447 corresponds to ICRC 2017; please correct the year.

Circularity Check

0 steps flagged · score 0.0 of 10

No circularity: the MRM background estimate is validated against source-free fields, and the Geminga measurement is explicitly deferred; no result is derived from its own target quantity.

full rationale

The paper does not derive any claimed result from the same quantity it is supposed to measure. The Matched Runs Method selects OFF runs using hardware-level matching parameters (ΔAz ≤ 10°, ΔEl ≤ 10°, ΔT ≤ 360 days, minimized ΔNCR in Eq. 3.1) and then validates the resulting significance maps of Segue 1 and Ursa Minor against a standard normal distribution. Those null fields provide an external empirical check; they are not fitted inputs used to force the conclusion. The matching thresholds are specified before the validation and are not tuned to make the blank-field maps appear null. The strongest claims, 'MRM has been demonstrated to be an effective technique' (abstract) and 'Validation tests with Segue 1 and Ursa Minor show that MRM is capable of estimating the background of the entire field-of-view' (§5), rest on the null tests, which are legitimate though limited in scope. The limitation that no test yet shows the method preserves a bright extended source like Geminga is a correctness/scope concern, not a circularity concern, and the paper itself acknowledges the Geminga application is pending ('The MRM procedure will be applied to the Geminga data set when the validation tests are completed'). Citations to prior VERITAS/HAWC papers are contextual or provide independent archival data (e.g., the dSph non-detection result [18]), and none of them is used to define the method's output in terms of its input. Therefore no circular step can be exhibited, and the appropriate score is 0.

Assumptions & free parameters 1 free parameters · 3 assumptions · 0 invented entities

The central method relies on the domain assumptions above; no independent calibration against a known extended source is provided. The only free parameters are the hand-chosen matching thresholds, which are not fitted to Geminga data.

free parameters (1)

Matched-run selection thresholds (Delta El, Delta Az, Delta T) = Delta El <= 10 deg, Delta Az <= 10 deg, Delta T <= 360 days
Hand-chosen thresholds used to select OFF runs from the archive; not fitted to Geminga data.

assumptions (3)

domain assumption A match in elevation, azimuth, time gap, and non-gamma-like event rate guarantees that the OFF-run background matches the ON-run background in the ROI.
Central assumption of MRM, introduced in Section 3; not directly tested by blank-field null maps.
domain assumption A significance distribution consistent with a standard normal in a source-free field proves unbiased background estimation across the full field of view.
Assumed in Section 4; it only establishes the null behavior, not background subtraction under a real extended source.
standard math Gamma-ray event counts follow Poisson statistics and the Li-Ma significance formula.
Implicit in the significance distribution test in Section 4; standard in IACT analyses.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Search for Very-High-Energy (E $>$ 100 GeV) Emission from Geminga Supernova by VERITAS." pith.science (2026). https://pith.science/paper/HHW2BJXW

@misc{pith2026190805369,
  author       = {Pith},
  title        = {Pith review of: Search for Very-High-Energy (E $>$ 100 GeV) Emission from Geminga Supernova by VERITAS},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/HHW2BJXW}},
  note         = {Machine review of arXiv:1908.05369}
}

read the original abstract

Geminga is a nearby (250 pc) middle-aged (spin-down time scale ~12,000 years) pulsar associated with a supernova remnant. Geminga has been a prime candidate for the origin of the unexpectedly high flux of cosmic-ray positrons above 10 GeV detected at Earth. Extended TeV gamma-ray emission from a 2-degree region around the Geminga pulsar was detected by the HAWC observatory, thus suggesting efficient, high-energy leptonic acceleration. Fermi-LAT observations show that the density of GeV leptons in the TeV nebula is lower than predicted by single zone and two zone diffusion models constrained with the HAWC measurements. However, the energy gap between Fermi-LAT and HAWC (~500 GeV to ~1 TeV) remains under-examined. The VERITAS gamma-ray observatory is sensitive in the energy range from 100 GeV to greater than 30 TeV, filling the gap between Fermi-LAT and HAWC. Therefore, VERITAS measurements potentially provide missing information. VERITAS has observed Geminga for 93 hours since 2009 including 28 hours in the 2018/2019 season. However, the standard VERITAS data analysis techniques have insufficient sensitivity to sources extended at the scale of the HAWC detection, due to difficulties with background estimation. We developed the Matched Runs Method (MRM) for VERITAS analysis of spatially extended sources. MRM has been demonstrated to be an effective technique by applying it to archival VERITAS data, and we are currently applying it to the Geminga observations. Here we present the summary of the MRM.

Figures

Figures reproduced from arXiv: 1908.05369 by the authors.

**Figure 2.** Sky map and the significance distribution of the ent [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

18 extracted references · 17 canonical work pages

[1]

2013, PRL, 111, 081102

Adriani, O., et al. 2013, PRL, 111, 081102

work page 2013
[2]

1996, ApJ, 457, L103

Golden, R.L., et al. 1996, ApJ, 457, L103

work page 1996
[3]

2001, ApJ, 559, 296

DuV ernois, M.A., et al. 2001, ApJ, 559, 296

work page 2001
[4]

2000, ApJ, 532, 653

Boezio, M., et al. 2000, ApJ, 532, 653

work page 2000
[5]

2007, Phys

Aguilar, M., et al. 2007, Phys. Lett. B, 646, 145

work page 2007
[6]

2004, PRL, 93, 241102

Beaty, J.J., et al. 2004, PRL, 93, 241102

work page 2004
[7]

2012, PRL, 108, 011103

Ackermann, M., et al. 2012, PRL, 108, 011103

work page 2012
[8]

Aguilar, M., et al., 2013, PRL, 110, 141102

work page 2013

Show all 18 references

[9]

E., et al

Fichtel, C. E., et al. 1975, ApJ, 198, 163

1975
[10]

1977, A & A, 56, 469

Bennett, K., et al. 1977, A & A, 56, 469

1977
[11]

2009, arXiv:0810.2784

Y uksel, H., et al. 2009, arXiv:0810.2784

2009 arXiv
[12]

2010, Astron

Delahaye, T., et al. 2010, Astron. Astrophys., 524, A51

2010
[13]

2009, ApJ, 700, L127

Abdo, A.A, et al. 2009, ApJ, 700, L127

2009
[14]

U., et al

Abeysekara, A. U., et al. 2017, ApJ, 843, 40

2017
[15]

Profumo, S. et al. 2018, PRD , 97, 123008

2018
[16]

U., et al

Abeysekara, A. U., et al. 2017, Science, 358, 911

2017
[17]

Berge, S

D. Berge, S. Funk, and J. Hinton, Background modelling i n very-high-energy γ -ray astronomy, A&A 466 (May, 2007) 1219-1229

2007
[18]

Benjamin Zitzer for the VERITAS Collaboration, ICRC pr oceeding 2007, arXiv:1708.07447 4

2007 arXiv

Pith tools

Reviewed August 14, 2026 · model on record in the stance chip above.

[1] [1]

2013, PRL, 111, 081102

Adriani, O., et al. 2013, PRL, 111, 081102

work page 2013

[2] [2]

1996, ApJ, 457, L103

Golden, R.L., et al. 1996, ApJ, 457, L103

work page 1996

[3] [3]

2001, ApJ, 559, 296

DuV ernois, M.A., et al. 2001, ApJ, 559, 296

work page 2001

[4] [4]

2000, ApJ, 532, 653

Boezio, M., et al. 2000, ApJ, 532, 653

work page 2000

[5] [5]

2007, Phys

Aguilar, M., et al. 2007, Phys. Lett. B, 646, 145

work page 2007

[6] [6]

2004, PRL, 93, 241102

Beaty, J.J., et al. 2004, PRL, 93, 241102

work page 2004

[7] [7]

2012, PRL, 108, 011103

Ackermann, M., et al. 2012, PRL, 108, 011103

work page 2012

[8] [8]

Aguilar, M., et al., 2013, PRL, 110, 141102

work page 2013