1 #define FRLIST_TARGET "input/L1_LDAS_C02_L2.frl" 2 #define FRLIST_WITNESS "input/L-R-9425.frl" 27 "L1:SUS-MMT1_SENSOR_LL",
28 "L1:SUS-MMT1_SENSOR_LR",
29 "L1:SUS-MMT1_SENSOR_UR",
30 "L1:SUS-MMT1_SENSOR_UL",
31 "L1:SUS-MMT2_SENSOR_LL",
32 "L1:SUS-MMT2_SENSOR_LR",
33 "L1:SUS-MMT2_SENSOR_UR",
34 "L1:SUS-MMT2_SENSOR_UL",
35 "L1:SUS-MMT3_SENSOR_LL",
36 "L1:SUS-MMT3_SENSOR_LR",
37 "L1:SUS-MMT3_SENSOR_UR",
38 "L1:SUS-MMT3_SENSOR_UL",
39 "L1:SUS-ITMX_SENSOR_LL",
40 "L1:SUS-ITMX_SENSOR_LR",
41 "L1:SUS-ITMX_SENSOR_UR",
42 "L1:SUS-ITMX_SENSOR_UL",
43 "L1:SUS-ITMY_SENSOR_LL",
44 "L1:SUS-ITMY_SENSOR_LR",
45 "L1:SUS-ITMY_SENSOR_UR",
46 "L1:SUS-ITMY_SENSOR_UL" 53 double fLow=110;
double fHigh=130;
56 double fLow=170;
double fHigh=190;
60 plot =
new watplot(const_cast<char*>(
"plot"),200,20,800,500);
62 plot->
goptions(const_cast<char*>(
"alp logy"), 1, 0., 0.,
true, fLow, fHigh,
true, 64);
81 cout <<
"yy : " << n <<
" " << yy[
n]->
size() <<
" " << yy[
n]->rate() << endl;
93 int level=xx.rate()/(2*
dFREQ);
100 rr.
add(*yy[0],const_cast<char*>(
"witness0"),59.,61.);
101 for(
int n=1;
n<
nWITNESS;
n++) rr.
add(*yy[
n],const_cast<char*>(
"witnessn"),0,3.);
102 for(
int n=1; n<
nWITNESS; n++) rr.
add(1,n+1,const_cast<char*>(
"bi-witness"));
133 plot->
gtitle(
"L1 regression @ 60Hz");
134 *plot >>
const_cast<char*
>(
"L1_regression_60Hz.png");
137 plot->gtitle(
"L1 regression @ 120Hz");
138 *plot >>
const_cast<char*
>(
"L1_regression_1200Hz.png");
141 plot->gtitle(
"L1 regression @ 180Hz");
142 *plot >>
const_cast<char*
>(
"L1_regression_180Hz.png");
void gtitle(TString title="", TString xtitle="", TString ytitle="")
void setMatrix(double edge=0., double f=1.)
wavearray< double > target
wavearray< double > getNoise()
frame frt(FRLIST_NAME, HCHANNEL)
size_t add(WSeries< double > &target, char *name, double fL=0., double fH=0.)
virtual void start(double s)
void setChName(TString chName)
virtual size_t size() const
void apply(double threshold=0., char c='a')
int getLayer(wavearray< DataType_t > &w, double n)
param: n - layer number
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)
void solve(double th, int nE=0, char c='s')
wavearray< double > getClean()
virtual void stop(double s)
void Forward(int n=-1)
param: wavelet - n is number of steps (-1 means full decomposition)
void mask(int n, double flow=0., double fhigh=0.)