5 #define FRLIST_NAME "input/H1_LDAS_C02_L2.frl" 6 #define CHANNEL_NAME "H1:LDAS-STRAIN" 18 plot =
new watplot(const_cast<char*>(
"plot"),200,20,800,500);
21 plot->
goptions(const_cast<char*>(
"alp logx logy"), 1, x.
start(), x.
start()+100.,
true, 10., 0.,
true, 8.);
26 const char *mkf_argv[11] = {
"mkf",
"-Bu",
"-Hp",
"-o",
"6",
"-a",
"200",
"-s",
"16384",
"-l",NULL};
27 Filter
filter(mkf_argc, mkf_argv);
31 cout <<
"nfiles : " << nfiles << endl;
35 cout << frName << endl;
48 cout << plName << endl;
void gtitle(TString title="", TString xtitle="", TString ytitle="")
void setFrName(TString frName)
cout<< endl;cout<< "ts size = "<< ts.size()<< " ts rate = "<< ts.rate()<< endl;tf.Forward(ts, wdm);int levels=tf.getLevel();cout<< "tf size = "<< tf.size()<< endl;double dF=tf.resolution();double dT=1./(2 *dF);cout<< "rate(hz) : "<< RATE<< "\ layers : "<< nLAYERS<< "\ dF(hz) : "<< dF<< "\ dT(ms) : "<< dT *1000.<< endl;int itime=TIME_PIXEL_INDEX;int ifreq=FREQ_PIXEL_INDEX;int index=(levels+1) *itime+ifreq;double time=itime *dT;double freq=(ifreq >0) ? ifreq *dF :dF/4;cout<< endl;cout<< "PIXEL TIME = "<< time<< " sec "<< endl;cout<< "PIXEL FREQ = "<< freq<< " Hz "<< endl;cout<< endl;wavearray< double > x
virtual void start(double s)
void goptions(char *opt=NULL, int col=1, double t1=0., double t2=0., bool fft=false, float f1=0., float f2=0., bool psd=false, float t3=0., bool oneside=false)
virtual void stop(double s)
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)