Bright Star Subtraction Pipeline for LSST: Phase one report

We present the phase one report of the Bright Star Subtraction (BSS) pipeline for the Vera C. Rubin Observatory's Legacy Survey of Space and Time (LSST). This pipeline is designed to create an extended PSF model by utilizing observed stars, followed by subtracting this model from the bright stars present in LSST data. Running the pipeline on Hyper Suprime-Cam (HSC) data shows a correlation between the shape of the extended PSF model and the position of the detector within the camera's focal plane. Specifically, detectors positioned closer to the focal plane's edge exhibit reduced circular symmetry in the extended PSF model. To mitigate this effect, we present an algorithm that enables users to account for the location dependency of the model. Our analysis also indicates that the choice of normalization annulus is crucial for modeling the extended PSF. Smaller annuli can exclude stars due to overlap with saturated regions, while larger annuli may compromise data quality because of lower signal-to-noise ratios. This makes finding the optimal annulus size a challenging but essential task for the BSS pipeline. Applying the BSS pipeline to HSC exposures allows for the subtraction of, on average, 100 to 700 stars brighter than 12th magnitude measured in g-band across a full exposure, with a full HSC exposure comprising ~100 detectors.