Clipping test gives false negative (algorithm false to detect clipping)
Some traces are passing the clipping detection test when they should not be. I assume the histogram algorithm should detect this clipping.
Example data: event.json KO.KIZT.xml KO.KIZT..HHE.mseed KO.KIZT..HHN.mseed KO.KIZT..HHZ.mseed
Code to run clipping test for example data:
#!/usr/bin/env python3
import pathlib
import obspy
import json
from gmprocess.core.stationstream import StationStream
from gmprocess.waveform_processing.clipping.clipping_check import check_clipping
from gmprocess.utils.event import get_event_object
EVENT = "us6000jlqa"
STATION = "KO.KIZT"
traces = obspy.read(f"{STATION}*.mseed")
inventory = obspy.read_inventory(f"{STATION}.xml")
station_stream = StationStream(traces, inventory)
with open("event.json") as fin:
ev_info = json.load(fin)
event = get_event_object(ev_info)
import pdb; pdb.set_trace()
check_clipping(station_stream, event)