10 #define WF_HP "Wolfy_hplus.txt" // Wolfy waveform 13 #define USE_WDM_00_90 // use 00 + 90 components 46 double rate = (time.
size()-1)/(time[time.
size()-1]-time[0]);
48 cout <<
"rate time series " << x.
rate() << endl;
49 cout <<
"size time series " << x.
size() << endl;
50 cout <<
"start time series " << x.
start()<< endl;
51 cout <<
"stop time series " << x.
stop() << endl;
62 cout <<
"E_TIME = " << E_TIME << endl;
70 double dF = R/(2*level);
73 cout<<
"R="<< R <<
" dF= " << dF <<
" dT= " << dT << endl;
79 cout <<
"size WDM " << w.
size() << endl;
80 cout <<
"start WDM " << w.
start() << endl;
81 cout <<
"stop WDM " << w.
stop() << endl;
88 float dt2 = 1./(2*df2);
104 geWDM->
data[
i-1] += en;
108 cout <<
"slices "<< slices <<
" layers"<< layers << endl;
109 cout <<
"E_WDM = " << E_WDM << endl;
110 cout <<
"E_WDM/E_TIME = " << E_WDM/E_TIME << endl;
124 plot->
gtitle(gtitle,
"time(sec)",
"amplitude");
134 WTS =
new watplot(const_cast<char*>(
"wtswrc"));
140 WTS->
plot(&w, 2, start, stop,const_cast<char*>(
"COLZ"));
142 WTS->
plot(&w, 4, start, stop,const_cast<char*>(
"COLZ"));
144 WTS->
hist2D->GetYaxis()->SetRangeUser(flow, fhigh);
146 WTS->
canvas->SaveAs(
"Wolfy_WDM.png");
157 for(
int i=0;
i<greT->
size()-1;
i++) {
172 gfile=
"Wolfy_Time_Freq_Comp.png";
size_t append(const wavearray< DataType_t > &)
void gtitle(TString title="", TString xtitle="", TString ytitle="")
virtual void rate(double r)
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
gwavearray< double > * geT
plot gtitle(gtitle,"frequency (Hz)","strain/#sqrt{Hz}")
virtual void start(double s)
void plot(wavearray< double > &, char *=NULL, int=1, double=0., double=0., bool=false, float=0., float=0., bool=false, float=0., bool=false)
gwavearray< double > * geWDM
virtual size_t size() const
gwavearray< double > * gx
double GetTimeRange(double &tMin, double &tMax, double efraction=0.9999999)
gwavearray< double > * greT
WDM< double > wdtf(lev, 2 *lev, 6, 10)
virtual void stop(double s)
void resample(const wavearray< DataType_t > &, double, int=6)
void Forward(int n=-1)
param: wavelet - n is number of steps (-1 means full decomposition)
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)
DataType_t getSample(int n, double m)
void Draw(GWAT_DRAW type=GWAT_TIME, TString options="ALP", Color_t color=kBlack)