Skip to content

Commit

Permalink
Merge pull request #407 from ahnitz/ifar_limits
Browse files Browse the repository at this point in the history
fix limits on pycbc_page_ifar search for timeslides
  • Loading branch information
spxiwh committed Oct 1, 2015
2 parents 756732f + 52355e5 commit 04a221b
Showing 1 changed file with 5 additions and 2 deletions.
7 changes: 5 additions & 2 deletions bin/hdfcoinc/pycbc_page_ifar
Original file line number Diff line number Diff line change
Expand Up @@ -97,9 +97,12 @@ xs, ys = pylab.poly_between(expected_ifar, error_minus, error_plus)
pylab.fill(xs, ys, facecolor='y', alpha=0.2, label='$2N^{1/2}$ Errors')

# get a unique list of timeslide_ids and loop over them
id_limit = (len(back_tsid)/2) // opts.decimation_factor
tsids = [x for x in range(-id_limit*opts.decimation_factor, (id_limit+1)*opts.decimation_factor, opts.decimation_factor) if x != 0]
min_tsid = numpy.ceil(back_tsid.min() / float(opts.decimation_factor))
max_tsid = int(back_tsid.max() / opts.decimation_factor)
tsids = numpy.arange(min_tsid, max_tsid, 1).astype(numpy.int64) * opts.decimation_factor
for tsid in tsids:
if tsid == 0:
continue

# find all triggers in this time slide
ts_indx = numpy.where(back_tsid == tsid)
Expand Down

0 comments on commit 04a221b

Please sign in to comment.