  We describe a novel, very fast and robust, directed search incoherent method
for periodic gravitational waves (GWs) from neutron stars in binary systems. As
directed search, we assume the source sky position to be known with enough
accuracy, but all other parameters are supposed to be unknown. We exploit the
frequency-modulation due to source orbital motion to unveil the signal
signature by commencing from a collection of time and frequency peaks. We
validate our pipeline adding 131 artificial continuous GW signals from pulsars
in binary systems to simulated detector Gaussian noise, characterized by a
power spectral density Sh = 4x10^-24 Hz^-1/2 in the frequency interval [70,
200] Hz, which is overall commensurate with the advanced detector design
sensitivities. The pipeline detected 128 signals, and the weakest signal
injected and detected has a GW strain amplitude of ~10^-24, assuming one month
of gapless data collected by a single advanced detector. We also provide
sensitivity estimations, which show that, for a single- detector data covering
one month of observation time, depending on the source orbital Doppler
modulation, we can detect signals with an amplitude of ~7x10^-25. By using
three detectors, and one year of data, we would easily gain more than a factor
3 in sensitivity, translating into being able to detect weaker signals. We also
discuss the parameter estimate proficiency of our method, as well as
computational budget, which is extremely cheap. In fact, sifting one month of
single-detector data and 131 Hz-wide frequency range takes roughly 2.4 CPU
hours. Due to the high computational speed, the current procedure can be
readily applied in ally-sky schemes, sieving in parallel as many sky positions
as permitted by the available computational power.
