3 #pragma GCC system_header
10 #include "TObjArray.h"
11 #include "TObjString.h"
21 cout <<
"-----> CWB_Plugin_WDM_freqCuts.C" << endl;
22 cout <<
"ifo " << ifo.Data() << endl;
23 cout <<
"type " << type << endl;
28 cout <<
"CWB_Plugin_freqCuts.C is implmenented only for 2G analysis" << endl;
34 double freqCuts[100][2];
35 freqCuts[0][0]=58; freqCuts[0][1]=62;
36 freqCuts[1][0]=78; freqCuts[1][1]=81;
37 freqCuts[2][0]=118; freqCuts[2][1]=122;
38 freqCuts[3][0]=176; freqCuts[3][1]=184;
39 freqCuts[4][0]=188; freqCuts[4][1]=191;
40 freqCuts[5][0]=198; freqCuts[5][1]=201;
44 for(
int level=cfg->
l_high; level>=cfg->
l_low; level--) {
58 double layer =
j+0.01;
59 for(
int k=0;
k<nfreqCuts;
k++) {
61 if((freq>freqCuts[
k][0])&&(freq<freqCuts[k][1])) {
monster wdmMRA
list of pixel pointers for MRA
void CWB_Plugin(TFile *jfile, CWB::config *cfg, network *net, WSeries< double > *x, TString ifo, int type)
COHERENCE.
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<< "\t layers : "<< nLAYERS<< "\t dF(hz) : "<< dF<< "\t 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
network ** net
NOISE_MDC_SIMULATION.
int getLayer(wavearray< DataType_t > &w, double n)
param: n - layer number
void Forward(int n=-1)
param: wavelet - n is number of steps (-1 means full decomposition)
void Inverse(int n=-1)
param: n - number of steps (-1 means full reconstruction)
void putLayer(wavearray< DataType_t > &, double n)
param: n - layer number