--- title: Band-Sampled-Data Directed Pipeline url: https://www.emergentmind.com/topics/band-sampled-data-directed-pipeline type: topic --- # Band-Sampled-Data Directed Pipeline Searching arXiv for the exact topic and closely related papers. Band-Sampled-Data Directed Pipeline most specifically denotes the **BSD-based semi-coherent continuous-wave search architecture** used for the Galactic Center search in Advanced LIGO O2 data, where **Band-Sampled-Data (BSD)** files provide reduced, band-organized detector data and a **directed** search fixes the sky position while scanning broad frequency and first spin-down ranges [1910.05097]. In that canonical usage, the pipeline combines BSD heterodyne-style Doppler correction, long coherent FFTs, peakmap construction, the **FrequencyHough** transform, inter-detector coincidence, and upper-limit setting for continuous gravitational waves from isolated asymmetric rotating neutron stars near **Sgr A\*** [1910.05097]. A broader use of the term is suggested by later sampled-data work in system identification, control, speech reconstruction, wideband conversion, and wireless sensing, where band-organized or nonuniform sampled data are routed through explicitly staged, direction-preserving processing chains; that broader reading is interpretive rather than standardized. ## 1. Defining meaning and scope In its primary and explicit sense, the Band-Sampled-Data directed pipeline is the search framework introduced for a **directed semi-coherent search** for continuous gravitational waves from the **Galactic Center** using Advanced LIGO **O2** data [1910.05097]. The search is **directed** because the sky location is assumed known—taken to be the position of **Sgr A\***—while the intrinsic signal frequency and first spin-down are unknown and searched over broad ranges [1910.05097]. It is therefore neither a fully targeted search, because the source rotational parameters are not known, nor an all-sky search, because only one sky position/bin is searched [1910.05097]. The target source model is the standard continuous-wave model of an **isolated spinning neutron star** that is non-axisymmetric about its rotation axis [1910.05097]. The intrinsic strain amplitude is written as \[ h_0=\frac{4\pi^2G}{c^4}\frac{I_{zz}f_{GW}^2}{d}\epsilon, \] with the usual dependence on moment of inertia, gravitational-wave frequency, source distance, and ellipticity [1910.05097]. Operationally, the observed signal is treated as nearly monochromatic but modulated by intrinsic spin-down, detector motion relative to the source, and antenna-pattern amplitude modulation [1910.05097]. A broader *Editor's term* use can be inferred from later literature: namely, a pipeline that ingests **band-limited, band-organized, or irregularly sampled data**, applies stage-wise harmonization or correction, and then performs a directed estimation, reconstruction, or detection task. This broader interpretation is suggested by sampled-data continuous-time identification under band-limited intersample assumptions [2409.09299], sampled-data observer design with predictor-observer decomposition [1911.07601], and multirate wideband DSP pipelines in hardware [2312.14392]. ## 2. BSD data architecture and processing chain The **BSD framework** is the reduced-data architecture on which the directed Galactic Center pipeline is built [1910.05097]. Each BSD file contains the **reprocessed strain time series** \(h(t)\), represented as a **complex time series**, and **down-sampled to 10 Hz** from the original 16 kHz data [1910.05097]. BSD files are organized by **10 Hz frequency bands** and by run sub-periods of about one month, so the data can be manipulated bandwise without repeatedly processing the full-rate detector stream [1910.05097]. For each BSD file covering one 10 Hz band and one sub-period, the single-detector pipeline follows a fixed sequence: **partial Doppler correction** at the assumed sky position using a modified BSD heterodyne in **1 Hz sub-bands**; **longer coherent FFTs** and peakmap construction; **FrequencyHough transform** from the peakmap to a \((f,\dot f)\) map; **summation of FrequencyHough maps** from all BSD files spanning the same frequency/spin-down region; **candidate selection** by standard FrequencyHough ranking; **inter-detector coincidence** between Hanford and Livingston candidate sets; **post-processing vetoes** including a significance threshold and known-line exclusion; **candidate inspection**; and finally **upper-limit setting** by Monte Carlo injections and sensitivity-depth extrapolation [1910.05097]. The BSD-specific contribution is the reduced data format, the complex low-rate band-limited data products, BSD heterodyne-style corrections, the new **partial Doppler correction in 1 Hz sub-bands** (“multi-Doppler”), and efficient handling in 10 Hz chunks and monthly pieces [1910.05097]. The FrequencyHough-specific contribution is the use of peakmaps, mapping of time-frequency peaks into the \((f,\dot f)\) plane, Hough-map construction and summation, candidate ranking by Hough number count and critical ratio, and coincidence logic in Hough-grid coordinates [1910.05097]. The overall architecture is therefore a BSD front-end wrapped around a FrequencyHough semi-coherent back-end. The multi-Doppler correction is implemented by multiplying each 1 Hz sub-band by \[ \exp\left(i\frac{2\pi}{c}p_{\vec n}f_i\right), \] where \(p_{\vec n}\) is the detector position projected along the source direction and \(f_i\) is the sub-band central frequency [1910.05097]. The paper reports that, within a **1 Hz** sub-band, the approximation is valid with a **maximum error of 5% with respect to the source frequency**, in the sense that the corrected signal remains in the expected frequency bin [1910.05097]. ## 3. Search model, resolution, and computational strategy The Galactic Center search fixes the sky position to **Sgr A\***, with RA(J2000) \(=17^\mathrm{h}\,45^\mathrm{m}\,40.04^\mathrm{s}\), Dec(J2000) \(=-29^\circ\,00'\,28.1''\), or ecliptic coordinates \((266.8517,-5.6077)^\circ\) [1910.05097]. Only this sky position is used in the search grid, so the total number of templates is simply the product of the number of frequency bins and spin-down bins [1910.05097]. The searched parameter space is \[ f \in [10,710]\ \mathrm{Hz}, \qquad \dot f \in [-1.8\times10^{-9},\,3.7\times10^{-11}] \ \mathrm{Hz/s}, \] including both spin-down and a small spin-up range [1910.05097]. After sky-position Doppler correction, the intrinsic frequency evolution is treated as effectively linear over the observing span, \[ f(t)\simeq f_0+\dot f(t-t_{\rm ref}), \] with higher derivatives not searched [1910.05097]. The FrequencyHough grid resolutions are \[ \delta f_{FH}=\frac{1}{T_{coh}K_f}, \qquad \delta \dot f_{FH}=\frac{1}{T_{coh}T_{obs}K_{\dot f}}, \] with \(K_f=10\) and \(K_{\dot f}=2\) in this search [1910.05097]. Coincidence between Hanford and Livingston candidates is defined by the normalized distance \[ d=\sqrt{\left(\frac{\Delta f}{\delta f_{FH}}\right)^2+\left(\frac{\Delta \dot f}{\delta \dot f_{FH}}\right)^2}, \] with threshold \(d_{thr}=4\) [1910.05097]. Even though only one sky bin is searched, the corresponding angular resolution is discussed through Doppler resolution. The number of frequency bins affected by Doppler modulation is reported as \(N_D=273\) at the lowest frequencies and \(N_D=1623\) at the highest frequencies, corresponding at 8 kpc to a Galactic-Center-centered patch radius from about **150 pc** at low frequency to **25 pc** at high frequency [1910.05097]. The central computational advantage is that partial Doppler correction before peakmap generation allows **longer coherence times** at fixed cost [1910.05097]. The coherence time scales as \(T_{coh}\propto 1/\sqrt{f_{max}}\), with reported values of **64208 s** for 10–20 Hz and **10776 s** for 700–710 Hz; crucially, the partial correction allows a coherence time **4 times longer** than without correction [1910.05097]. The effective Hough-grid frequency resolution ranges from \(1.6\times10^{-6}\,\mathrm{Hz}\) to \(9.3\times10^{-6}\,\mathrm{Hz}\), while the spin-down natural resolution ranges from \(3.3\times10^{-13}\) to \(2.0\times10^{-12}\,\mathrm{Hz/s}\) for Hanford and from \(3.8\times10^{-13}\) to \(2.3\times10^{-12}\,\mathrm{Hz/s}\) for Livingston [1910.05097]. The total template counts are reported as **\(2.4\times10^{11}\)** for Livingston and **\(2.7\times10^{11}\)** for Hanford, yet the search required only **207 jobs per detector**, about **30 min/job**, and about **200 core hours** total, excluding BSD production [1910.05097]. ## 4. Galactic Center O2 deployment and empirical performance The search used **Advanced LIGO O2 open data** from **Hanford** and **Livingston** [1910.05097]. O2 ran from **2016-11-30** to **2017-08-25**; only science segments of the latest calibrated data were used, poor-quality periods were excluded, Livingston data before **2017-01-04** were discarded, and **35 days** of Hanford data from mid-March to mid-April were excluded [1910.05097]. Virgo was not used because of shorter observing time and lower sensitivity [1910.05097]. The final search employed **1120 BSD files** spanning **70 frequency bands** of width 10 Hz across 10–710 Hz [1910.05097]. Peakmap construction uses the standard FrequencyHough threshold with noise-peak probability \[ p_0=0.0755 \] [1910.05097]. Candidate ranking retains approximately **1000 candidates per job**, producing **203961 candidates** for Livingston and **202556 candidates** for Hanford [1910.05097]. Coincidence with \(d_{thr}=4\) yields **237 coincident candidates**; a second-stage veto uses the **Critical Ratio** \[ CR=\frac{n-Np_0}{\sqrt{Np_0(1-p_0)}}, \] with a threshold chosen so that, on average, only **one false candidate over the full parameter space** is expected [1910.05097]. The threshold range is \(CR_{thr}\in[6.00,6.55]\) for Hanford and \(CR_{thr}\in[5.98,6.53]\) for Livingston [1910.05097]. After the critical-ratio cut, only **9** coincident candidates survive [1910.05097]. Of these, **4** are due to known instrumental lines, **1** is due to the hardware injection **Pulsar_10**, and the remaining **4** were judged non-astrophysical because they were associated with transient lines, especially in Livingston before 14 March 2017 [1910.05097]. No astrophysical continuous-wave candidate was found [1910.05097]. Sensitivity is reported through upper limits and sensitivity depth. Using injection campaigns in 26 clean 1 Hz bands, the mean depths are about **44.30 \(1/\sqrt{\mathrm{Hz}}\)** for Hanford and **52.44 \(1/\sqrt{\mathrm{Hz}}\)** for Livingston [1910.05097]. The detection efficiency is fit with \[ D(x)=K\left(1-e^{-A_1(x-x_{min})^{A_2}}\right), \] and the \(95\%\) upper limit is defined by \(D=0.95\) [1910.05097]. The most stringent \(95\%\) confidence upper limits are **\(\sim 1.4\times10^{-25}\)** near **161 Hz** for Livingston and **\(\sim 1.6\times10^{-25}\)** near **195 Hz** for Hanford; these do **not** include calibration uncertainty [1910.05097]. With \(d=8\) kpc and fiducial \(I_{zz}=10^{38}\,\mathrm{kg\,m^2}\), the strongest ellipticity constraint is about **\(\sim 4\times10^{-6}\)** at the highest frequencies for Livingston [1910.05097]. The paper characterizes the result as the **most sensitive directed search** for Galactic Center continuous waves to date and the **first such search using O2 data** [1910.05097]. ## 5. Broader sampled-data formulations in estimation and control Outside gravitational-wave astronomy, the phrase can be generalized only cautiously. A plausible broader reading is a staged architecture that starts from sampled data whose intersample, spectral, or timing structure matters, and then performs explicitly directed estimation or control. In **continuous-time system identification from sampled data**, the decisive step is to specify both the **intersample behavior** and the **past behavior** of the continuous-time input [2409.09299]. For **band-limited** input under the Nyquist frequency and **periodically appended** past behavior, the input can be reconstructed exactly from the samples by a finite Fourier series, and the kernel-based regularization estimator of the continuous-time impulse response acquires a **closed form** [2409.09299]. The estimator is \[ \hat g(\tau)=\pmb{\Sigma}_{gy_0}(\tau)(\pmb{\Sigma}_{y_0}+\gamma I_N)^{-1}\pmb y, \] and the paper reports from Monte Carlo experiments that the broader sampled-data KRM framework is **more robust** than SRIVC and PEM and **more accurate when the sample size is small**, although the experiments focus on the ZOH case rather than the band-limited case directly [2409.09299]. This suggests a band-sampled-data pipeline in which sampled inputs are first lifted into a continuous-time model class before identification. In **sampled-data observer design**, the central architecture is a **predictor-observer decomposition**: a continuous-time observer is driven not by the unavailable continuous output \(y(t)\), but by an auxiliary predictor state \(w(t)\) that is reset at sample times and propagated between samples [1911.07601]. The general predictor is \[ \dot w(t)=\nabla h(z(t))f(z(t),u(t)) - K(z(t),w(t),u(t))(w(t)-h(z(t))), \] with reset \(w(t_k)=y(t_k)\), while the observer evolves as \[ \dot z(t)=f(z(t),u(t))+g(z(t),w(t),u(t))(w(t)-h(z(t))) . \] Under an IOS-type assumption on the underlying continuous-time observer and the sampling constraint \[ 2\gamma L\int_0^T \exp(2qs)\,ds<1, \] the sampled-data observer inherits exponential convergence in the noiseless case and robustness to measurement noise [1911.07601]. Here the “directed” character lies in the explicit information flow sample \(\to\) predictor \(\to\) observer \(\to\) state estimate. A related directed sampled-data architecture appears in **formation control with local measurements** under directed graphs [1909.04819]. The sampled-data controller is implemented with synchronous zero-order hold, requires only local-frame relative measurements, and admits local exponential stability for a static target when \[ 0