Commit d73c5cf5 authored by Reinhard Prix's avatar Reinhard Prix
Browse files

HierarchSearchGCT: fixed potential problem in SignalOnly handling

- number of IFOs derived from total catalog instead of per-segment
  catalog, which could lead to errors in CreateFstatInput() when
  number of IFOs in a segment is less than for others
- refs #1531
Original: c7d1655f6a882591aa8b214a8a64ac62c776991a
parent cabc3353
......@@ -2041,7 +2041,8 @@ void SetUpSFTs( LALStatus *status, /**< pointer to LALStatus structure */
/* if single-only flag is given, assume a PSD with sqrt(S) = 1.0 */
MultiNoiseFloor assumeSqrtSX, *p_assumeSqrtSX;
if ( in->SignalOnly ) {
assumeSqrtSX.length = XLALCountIFOsInCatalog(catalog);
const SFTCatalog *catalog_k = &(catalogSeq.data[k]);
assumeSqrtSX.length = XLALCountIFOsInCatalog ( catalog_k );
for (UINT4 X = 0; X < assumeSqrtSX.length; ++X) {
assumeSqrtSX.sqrtSn[X] = 1.0;
}
......
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment