pycbc icon indicating copy to clipboard operation
pycbc copied to clipboard

`pycbc_single_template --psd-output` output .hdf file error

Open Fergus-Hayes opened this issue 6 years ago • 7 comments

Calling the command:

$ pycbc_single_template --segment-length 528 --segment-start-pad 144 --segment-end-pad 16 --psd-estimation median --psd-segment-length 16 --psd-segment-stride 8 --psd-inverse-length 16 --approximant "SPAtmplt:mtotal<4" "SEOBNRv4_ROM:else" --order -1 --psd-num-segments 63 --taper-data 1 --allow-zero-padding --autogating-threshold 100 --autogating-cluster 0.5 --autogating-width 0.25 --autogating-taper 0.25 --autogating-pad 16 --minimum-chisq-bins 3 --strain-high-pass 15. --pad-data 8 --template-start-frequency 27 --chisq-bins '0.72*get_freq("fSEOBNRv4Peak",params.mass1,params.mass2,params.spin1z,params.spin2z)**0.7' --sample-rate 2048 --mass1 19.19753 --mass2 1.0981159 --spin1z -0.33484802 --spin2z 0.011133107 --low-frequency-cutoff 25.8 --gps-start-time 1169299588 --gps-end-time 1169300612 --trigger-time 1169300100.8 --window 1.48876953125 --frame-type H1_HOFT_C00 --channel-name H1:GDS-CALIB_STRAIN --psd-output ./PSD_1169300100_H1.hdf --output-file ./SNRTS_1169300100_H1.hdf

brings up the error:

  return self._data.__rdiv__(other)
Traceback (most recent call last):
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/bin/pycbc_single_template", line 4, in <module>
    __import__('pkg_resources').run_script('PyCBC==1.9.4', 'pycbc_single_template')
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/local/lib/python2.7/site-packages/pkg_resources/__init__.py", line 658, in run_script
    self.require(requires)[0].run_script(script_name, ns)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/local/lib/python2.7/site-packages/pkg_resources/__init__.py", line 1438, in run_script
    exec(code, namespace, namespace)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/deploy/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/lib/python2.7/site-packages/PyCBC-1.9.4-py2.7.egg/EGG-INFO/scripts/pycbc_single_template", line 238, in <module>
    precision='single')
  File "/cvmfs/oasis.opensciencegrid.org/ligo/deploy/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/lib/python2.7/site-packages/PyCBC-1.9.4-py2.7.egg/pycbc/psd/__init__.py", line 464, in associate_psds_to_segments
    precision=precision)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/deploy/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/lib/python2.7/site-packages/PyCBC-1.9.4-py2.7.egg/pycbc/psd/__init__.py", line 425, in generate_overlapping_psds
    dyn_range_factor=dyn_range_factor, precision=precision)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/deploy/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/lib/python2.7/site-packages/PyCBC-1.9.4-py2.7.egg/pycbc/psd/__init__.py", line 113, in from_cli
    (psd.astype(float64) / (dyn_range_factor ** 2)).save(opt.psd_output)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/deploy/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/lib/python2.7/site-packages/PyCBC-1.9.4-py2.7.egg/pycbc/types/frequencyseries.py", line 418, in save
    compression_opts=9, shuffle=True)
  File "/cvmfs/oasis.opensciencegrid.org/ligo/deploy/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/lib/python2.7/site-packages/h5py/_hl/group.py", line 109, in create_dataset
    self[name] = dset
  File "/cvmfs/oasis.opensciencegrid.org/ligo/deploy/sw/pycbc/x86_64_deb_8/virtualenv/pycbc-v1.9.4/lib/python2.7/site-packages/h5py/_hl/group.py", line 277, in __setitem__
    h5o.link(obj.id, self.id, name, lcpl=lcpl, lapl=self._lapl)
  File "h5py/_objects.pyx", line 54, in h5py._objects.with_phil.wrapper
  File "h5py/_objects.pyx", line 55, in h5py._objects.with_phil.wrapper
  File "h5py/h5o.pyx", line 202, in h5py.h5o.link
RuntimeError: Unable to create link (name already exists)

while removing the --psd-output option to .hdf will produce the snr time series fine.

Fergus-Hayes avatar May 16 '18 14:05 Fergus-Hayes