This repository contains the fully reproducible likelihood and parameter-scan pipeline used in the phenomenological analysis of baseline- and energy-dependent phase corrections in reactor neutrino oscillations. The pipeline implements: • Daya Bay near–far spectral-ratio likelihood fits • KamLAND covariance-matrix spectral fits • One-parameter and multi-parameter phase-correction scans • Profile-likelihood construction with nuisance-parameter profiling • Global likelihood combination • Automated generation of figures and tables used in the manuscript The Daya Bay and KamLAND inputs are based on publicly released experimental spectral data products. The JUNO-like configuration included in the pipeline is a sensitivity-style simulation with multi-baseline averaging and detector-response smearing. It is not based on real JUNO data. All numerical results, χ² profiles, confidence intervals, tables, and figures are produced directly by this pipeline without manual post-processing. This software is intended for reproducible research and independent verification of the analysis results reported in the associated paper.