27 double*
tmp =
new double[pWDM->
nWWS/2];
28 for(
int j=1*32;
j<19*32; ++
j)
29 for(
int i=0;
i<=32; ++
i)
33 double* TFMap = pWDM->
pWWS;
34 for(
int j=32;
j<19*32; ++
j)
35 for(
int i=0;
i<=32; ++
i)
36 TFMap[
j*33+
i] = tmp[
j*33+
i];
45 TH1F*
h =
new TH1F(
"h",
"Time delay filter relative error", 100, -0.01, 0.01);
46 for(
int j=1024*2;
j<18*1024; ++
j)
47 h->Fill(ts2[
j]/
ts[
j-13] - 1);
virtual void rate(double r)
WDM< double > wdm(nLAYERS, nLAYERS, 6, 10)
void setTDFilter(int nCoeffs, int L=1)
double getPixelAmplitude(int, int, int, bool=false)
wavearray< double > ts(N)
int getLayer(wavearray< DataType_t > &w, double n)
param: n - layer number
unsigned long nWWS
pointer to wavelet work space
void Forward(int n=-1)
param: wavelet - n is number of steps (-1 means full decomposition)
WaveDWT< DataType_t > * pWavelet
void Inverse(int n=-1)
param: n - number of steps (-1 means full reconstruction)