clear all; close all; clc;
start_time = [162 23 10 0];
window = 5000;
das = '9FCC';

[seisTREX,hdr] = segycut(start_time,window,das,'3');

tvec=linspace(0,window,window*hdr.sps);

conv_factor_1INCH_DM_V_V2PA=  1/800;
vect_signalPA=-seisTREX*conv_factor_1INCH_DM_V_V2PA*hdr.atod;

npoles = 2

high = .25

[B2,A2] = butter(npoles,[high]/(hdr.sps/2),'high')

highpassV = filter(B2,A2,vect_signalPA);

figure(2)
plot(tvec,highpassV);
axis tight;

set(gcf,'Position',[0 0 700 150])