`find.binding.positions` return empty results for narrow peak call

[liuyang2006]

Dear author,

I follow the tutorial (https://hbctraining.github.io/Intro-to-ChIPseq/lessons/peak_calling_spp.html) to call narrow peaks from CUT & TAG data, however, the returned results is always empty (0 peaks detected).

When I call broad peaks, there are lots of results.

Would you guide me how to find the issues?

Here is the codes snippet:

# determine binding positions
bp <- find.binding.positions(signal.data = chip_filt,
                             control.data = input_filt,
                             fdr = use_fdr,
                             whs = detection.window.halfsize)

print(paste("detected", sum(unlist(lapply(bp$npl, function(d) length(d$x)))), "peaks"))

Best, Yang