5 #define HCHANNEL "H1:LSC-STRAIN"
6 #define FRLIST_NAME "S5_H1_RDS_C03_L2.frl"
8 #define FRLIST_NAME2 "H-R-8159.frl"
10 #define START 815974196
19 TString auch[]={
"H0:PEM-COIL_MAGZ",
"H0:PEM-BSC10_MAGX",
"H0:PEM-BSC10_MAGY",
"H0:PEM-BSC10_MAGZ"};
24 "H1:SUS-ETMX_COIL_UR",
"H1:SUS-ETMY_COIL_UR",
25 "H1:SUS-ETMY_SENSOR_UL",
"H1:SUS-ETMY_SENSOR_LL",
26 "H1:SUS-ETMY_SENSOR_UR",
"H1:SUS-ETMY_SENSOR_LR",
27 "H1:SUS-ETMX_SENSOR_UL",
"H1:SUS-ETMX_SENSOR_LL",
28 "H1:SUS-ETMX_SENSOR_UR",
"H1:SUS-ETMX_SENSOR_LR"};
54 r.
add(tfmap,
"hchannel");
68 r.
add(x,auch[i].Data());
69 cout << auch[
i].Data() << endl;
71 for (
int j=0;
j<
NB;
j++)
76 r.
add(x,channels[j].Data(),0,10);
77 cout << channels[
j].Data() << endl;
80 for (
int i=0;
i<
NA;
i++)
for (
int j=0; j<
NB; j++) r.
add(
i+1,NA+1+j,channels[j].Data());
83 for (
int j=0; j<
NB; j++) r.
mask(NA+1+j);
96 plot >>
"Regression_H1_bic.png";
101 for (
int i=0;
i<freq.
size();
i++) cout << freq.
data[
i] <<
" " << rank.
data[
i] << endl;
virtual size_t size() const
void setMatrix(double edge=0., double f=1.)
frame frt(FRLIST_NAME, HCHANNEL)
WDM< double > wdtf(lev, 2 *lev, 6, 10)
size_t add(WSeries< double > &target, char *name, double fL=0., double fH=0.)
virtual void start(double s)
void setChName(TString chName)
void apply(double threshold=0., char c='a')
wavearray< double > getRank(int n)
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 > vfreq
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 unmask(int n, double flow=0., double fhigh=0.)
void mask(int n, double flow=0., double fhigh=0.)