Definition at line 27 of file wseries.hh.
Public Member Functions | |
WSeries () | |
WSeries (const Wavelet &w) | |
WSeries (const wavearray< DataType_t > &value, const Wavelet &w) | |
param: value - data to initialize the WSeries object More... | |
WSeries (const WSeries< DataType_t > &value) | |
param: value - object to copy from More... | |
virtual | ~WSeries () |
void | bandpass (wavearray< DataType_t > &ts, double flow, double fhigh, int n=-1) |
void | bandpass (double flow, double fhigh, double a=0.) |
virtual void | Browse (TBrowser *b) |
virtual WSeries< double > | calibrate (size_t, double, d_complex *, d_complex *, wavearray< double > &, wavearray< double > &, size_t ch=0) |
param: number of samples in calibration arrays R & C param: frequency resolution param: pointer to response function R in Fourier domain param: pointer to sensing function C in Fourier domain param: time dependent calibration coefficient alpha param: time dependent calibration coefficient gamma param: 0/1 - AS_Q/DARM_ERR calibration, by default is 0 return array with calibration constants for each wavelet layer More... | |
virtual double | coincidence (WSeries< DataType_t > &, int=0, int=0, double=0.) |
param: WSeries object used for coincidence param: coincidence window in seconds return pixel occupancy More... | |
virtual double | Coincidence (WSeries< DataType_t > &, double=0., double=0.) |
param: WSeries object used for coincidence param: coincidence window in seconds param: threshold on significance return pixel occupancy More... | |
virtual void | Dump (const char *, int=0) |
virtual wavearray< double > | filter (size_t) |
param: n - number of decomposition steps algorithm: 1) do forward wavelet transform with n decomposition steps 2) whiten wavelet layers and calculate noise rms as 1/Sum(1/var) 3) do inverse wavelet transform with n reconstruction steps More... | |
void | Forward (int n=-1) |
param: wavelet - n is number of steps (-1 means full decomposition) More... | |
void | Forward (wavearray< DataType_t > &, int n=-1) |
void | Forward (wavearray< DataType_t > &, Wavelet &, int n=-1) |
virtual double | fraction (double=0., double=0., int=0) |
param: t - sub interval duration. If can not divide on integer More... | |
double | frequency (int l) |
double | Gamma2Gauss (TH1F *=NULL) |
double | getbpp () const |
double | gethigh () const |
int | getLayer (wavearray< DataType_t > &w, double n) |
param: n - layer number More... | |
int | getLevel () |
double | getlow () const |
int | getMaxLevel () |
DataType_t | getSample (int n, double m) |
std::slice | getSlice (double n) |
virtual double | gSignificance (double, double=1., double=0.) |
param: T - sliding window duration in seconds param: f - black pixel fraction param: t - sliding step in seconds options: f = 0 - returns black pixel occupancy options: t = 0 - sliding step = wavelet time resolution. More... | |
void | Inverse (int n=-1) |
param: n - number of steps (-1 means full reconstruction) More... | |
bool | isWDM () |
int | layer (double f) |
virtual void | lprFilter (double, int=0, double=0., double=0.) |
double | maxEnergy (wavearray< DataType_t > &ts, Wavelet &w, double=0, int=1, int=0, TH1F *=NULL) |
size_t | maxIndex () |
int | maxLayer () |
virtual void | median (double t, bool norm=false) |
void | mul (WSeries< DataType_t > &) |
virtual WSeries< DataType_t > & | operator*= (WSeries< DataType_t > &) |
virtual WSeries< DataType_t > & | operator*= (wavearray< DataType_t > &) |
virtual WSeries< DataType_t > & | operator*= (const DataType_t) |
virtual WSeries< DataType_t > & | operator+= (WSeries< DataType_t > &) |
virtual WSeries< DataType_t > & | operator+= (wavearray< DataType_t > &) |
virtual WSeries< DataType_t > & | operator+= (const DataType_t) |
virtual WSeries< DataType_t > & | operator-= (WSeries< DataType_t > &) |
virtual WSeries< DataType_t > & | operator-= (wavearray< DataType_t > &) |
virtual WSeries< DataType_t > & | operator-= (const DataType_t) |
WSeries< DataType_t > & | operator= (const wavearray< DataType_t > &) |
WSeries< DataType_t > & | operator= (const WSeries< DataType_t > &) |
WSeries< DataType_t > & | operator= (const DataType_t) |
virtual WSeries< DataType_t > & | operator[] (const std::slice &) |
virtual double | percentile (double=0., int=0, WSeries< DataType_t > *=NULL) |
param: f - black pixel fraction param: m - mode options: f = 0 - returns black pixel occupancy m = 1 - set threshold f, returns percentile amplitudes m =-1 - set threshold f, returns wavelet amplitudes m > 1 - random policy,returns percentile amplitudes m <-1 - random policy,returns wavelet amplitudes m = 0 - random pixel selection if m<0 return wavelet amplitudes instead of the percentile amplitude More... | |
virtual double | pixclean (double=0.) |
param: S - threshold on pixel significance return pixel occupancy. More... | |
void | print () |
param: int n if n<0, zero pixels defined in mask (regression) if n>=0, zero all pixels except ones defined in the mask param: bool - if true, set WSeries data to be positive if pMask.size()=0, mask(0,true) is equivalent to abs(data) return core pixel occupancy More... | |
void | putLayer (wavearray< DataType_t > &, double n) |
param: n - layer number More... | |
void | putSample (DataType_t a, int n, double m) |
virtual void | resample (double, int=6) |
virtual void | resize (unsigned int) |
double | resolution (int=0) |
virtual double | rsignificance (size_t=0, double=1.) |
param: n - sub-interval duration in domain units param: f - black pixel fraction options: f = 0 - returns black pixel occupancy More... | |
virtual double | rSignificance (double, double=1., double=0.) |
param: T - sliding window duration in seconds param: f - black pixel fraction param: t - sliding step in seconds options: f = 0 - returns black pixel occupancy options: t = 0 - sliding step = wavelet time resolution. More... | |
void | setbpp (double f) |
void | sethigh (double f) |
void | setLevel (size_t n) |
void | setlow (double f) |
void | setWavelet (const Wavelet &w) |
virtual double | significance (double, double=1.) |
param: n - sub-interval duration in seconds param: f - black pixel fraction options: f = 0 - returns black pixel occupancy More... | |
size_t | sizeZero () |
virtual WSeries< float > | variability (double=0., double=-1., double=-1.) |
param: first - time window to calculate normalization constants second - low frequency boundary for correction third - high frequency boundary for correction algorithm: 1) sort wavelet amplitudes with the same time stamp 2) calculate left(p) and right(p) amplitudes put (right(p)-left(p))/2 into output array 3) if first parameter >0 - devide WSeries by average variability More... | |
void | wavescan (WSeries< DataType_t > **, int, TH1F *=NULL) |
double | wdmPacket (int pattern, char opt='L', TH1F *=NULL) |
virtual WSeries< double > | white (double, int, double=0., double=0.) |
what it does: each wavelet layer is devided into k intervals. More... | |
virtual bool | white (WSeries< double > ws, int mode=0) |
void | wrate (double r) |
double | wrate () const |
size_t | xsize () |
![]() | |
wavearray (int) | |
wavearray () | |
wavearray (const wavearray< DataType_t > &) | |
template<class T > | |
wavearray (const T *, unsigned int, double=0.) | |
virtual | ~wavearray () |
void | add (const wavearray< DataType_t > &, int=0, int=0, int=0) |
size_t | append (const wavearray< DataType_t > &) |
size_t | append (DataType_t) |
void | cpf (const wavearray< DataType_t > &, int=0, int=0, int=0) |
virtual void | delay (double T) |
virtual void | DumpBinary (const char *, int=0) |
virtual void | DumpObject (const char *) |
virtual void | DumpShort (const char *, int=0) |
virtual void | edge (double s) |
virtual double | edge () const |
virtual void | exponential (double) |
virtual void | FFT (int=1) |
virtual void | FFTW (int=1) |
DataType_t | get (size_t i) |
DataType_t | get (double t) |
virtual wavearray< double > | getLPRFilter (size_t, size_t=0) |
virtual int | getSampleRank (size_t n, size_t l, size_t r) const |
virtual int | getSampleRankE (size_t n, size_t l, size_t r) const |
virtual std::slice | getSlice () const |
double | getStatistics (double &mean, double &rms) const |
void | hann (void) |
virtual size_t | limit () const |
virtual size_t | limit (const std::slice &) const |
virtual size_t | limit (const wavearray< DataType_t > &) const |
virtual void | lprFilter (wavearray< double > &) |
virtual DataType_t | max () const |
virtual void | max (wavearray< DataType_t > &) |
virtual double | mean () const |
virtual double | mean (double f) |
virtual double | mean (const std::slice &) |
virtual void | mean (double t, wavearray< DataType_t > *in, bool fl=false, size_t n=1) |
virtual double | median (size_t=0, size_t=0) const |
virtual void | median (double t, wavearray< DataType_t > *in, bool fl=false, size_t n=1) |
virtual DataType_t | min () const |
virtual wavearray< DataType_t > & | operator<< (wavearray< DataType_t > &) |
wavearray< DataType_t > & | operator= (const wavearray< DataType_t > &) |
wavearray< DataType_t > & | operator= (const DataType_t) |
virtual char * | operator>> (char *) |
virtual DataType_t & | operator[] (const unsigned int) |
void | print () |
long | rand48 (long k=1024) |
DataType_t | rank (double=0.5) const |
virtual void | rate (double r) |
virtual double | rate () const |
virtual void | ReadBinary (const char *, int=0) |
virtual void | ReadShort (const char *) |
void | Resample (const wavearray< DataType_t > &, double, int=6) |
void | resample (const wavearray< DataType_t > &, double, int=6) |
virtual void | Resample (double) |
virtual void | resetFFTW () |
virtual double | rms () |
virtual double | rms (const std::slice &) |
virtual void | rms (double t, wavearray< DataType_t > *in, bool fl=false, size_t n=1) |
virtual void | setSlice (const std::slice &s) |
virtual size_t | size () const |
virtual void | spesla (double, double, double=0.) |
virtual void | SQRT () |
double | Stack (const wavearray< DataType_t > &, int) |
double | Stack (const wavearray< DataType_t > &, int, int) |
double | Stack (const wavearray< DataType_t > &, double) |
virtual void | start (double s) |
virtual double | start () const |
virtual void | stop (double s) |
virtual double | stop () const |
void | sub (const wavearray< DataType_t > &, int=0, int=0, int=0) |
long | uniform () |
size_t | wavecount (double x, int n=0) |
virtual void | waveSort (DataType_t **pp, size_t l=0, size_t r=0) const |
virtual void | waveSort (size_t l=0, size_t r=0) |
virtual void | waveSplit (DataType_t **pp, size_t l, size_t r, size_t m) const |
virtual DataType_t | waveSplit (size_t l, size_t r, size_t m) |
virtual void | waveSplit (size_t m) |
virtual wavearray< double > | white (double, int=0, double=0., double=0.) const |
Public Attributes | |
double | bpp |
double | f_high |
double | f_low |
WaveDWT< DataType_t > * | pWavelet |
size_t | w_mode |
double | wRate |
![]() | |
DataType_t * | data |
double | Edge |
TFFTRealComplex * | fftw |
TFFTComplexReal * | ifftw |
pointer to direct fftw object More... | |
double | Rate |
size_t | Size |
data array More... | |
std::slice | Slice |
double | Start |
double | Stop |
Additional Inherited Members | |
![]() | |
static int | compare (const void *x, const void *y) |
pointer to inverse fftw object More... | |
#include <wseries.hh>
Definition at line 23 of file wseries.cc.
Definition at line 35 of file wseries.cc.
|
explicit |
param: value - data to initialize the WSeries object
Definition at line 48 of file wseries.cc.
param: value - object to copy from
Definition at line 61 of file wseries.cc.
Definition at line 76 of file wseries.cc.
void WSeries< DataType_t >::bandpass | ( | wavearray< DataType_t > & | ts, |
double | flow, | ||
double | fhigh, | ||
int | n = -1 |
||
) |
Definition at line 295 of file wseries.cc.
void WSeries< DataType_t >::bandpass | ( | double | flow, |
double | fhigh, | ||
double | a = 0. |
||
) |
Definition at line 322 of file wseries.cc.
|
inlinevirtual |
Reimplemented from wavearray< DataType_t >.
Reimplemented in gWSeries< DataType_t >.
Definition at line 433 of file wseries.hh.
|
virtual |
param: number of samples in calibration arrays R & C param: frequency resolution param: pointer to response function R in Fourier domain param: pointer to sensing function C in Fourier domain param: time dependent calibration coefficient alpha param: time dependent calibration coefficient gamma param: 0/1 - AS_Q/DARM_ERR calibration, by default is 0 return array with calibration constants for each wavelet layer
Definition at line 2189 of file wseries.cc.
|
virtual |
param: WSeries object used for coincidence param: coincidence window in seconds return pixel occupancy
Definition at line 911 of file wseries.cc.
|
virtual |
param: WSeries object used for coincidence param: coincidence window in seconds param: threshold on significance return pixel occupancy
Definition at line 1002 of file wseries.cc.
|
virtual |
Reimplemented from wavearray< DataType_t >.
Definition at line 849 of file wseries.cc.
param: n - number of decomposition steps algorithm: 1) do forward wavelet transform with n decomposition steps 2) whiten wavelet layers and calculate noise rms as 1/Sum(1/var) 3) do inverse wavelet transform with n reconstruction steps
Definition at line 1242 of file wseries.cc.
param: wavelet - n is number of steps (-1 means full decomposition)
Definition at line 228 of file wseries.cc.
void WSeries< DataType_t >::Forward | ( | wavearray< DataType_t > & | x, |
int | n = -1 |
||
) |
Definition at line 248 of file wseries.cc.
void WSeries< DataType_t >::Forward | ( | wavearray< DataType_t > & | x, |
Wavelet & | w, | ||
int | n = -1 |
||
) |
Definition at line 261 of file wseries.cc.
|
virtual |
param: t - sub interval duration. If can not divide on integer
param: f - black pixel fraction param: m - mode options: f = 0, m = 0 - returns black pixel occupancy m = 1 - set threshold f m = 2 - random policy m = 0 - random pixel selection
Definition at line 1356 of file wseries.cc.
Definition at line 99 of file wseries.cc.
double WSeries< DataType_t >::Gamma2Gauss | ( | TH1F * | hist = NULL | ) |
Definition at line 558 of file wseries.cc.
|
inline |
Definition at line 99 of file wseries.hh.
|
inline |
Definition at line 118 of file wseries.hh.
int WSeries< DataType_t >::getLayer | ( | wavearray< DataType_t > & | w, |
double | n | ||
) |
param: n - layer number
Definition at line 175 of file wseries.cc.
Definition at line 91 of file wseries.hh.
|
inline |
Definition at line 111 of file wseries.hh.
Definition at line 85 of file wseries.cc.
Definition at line 167 of file wseries.hh.
|
inline |
Definition at line 134 of file wseries.hh.
|
virtual |
param: T - sliding window duration in seconds param: f - black pixel fraction param: t - sliding step in seconds options: f = 0 - returns black pixel occupancy options: t = 0 - sliding step = wavelet time resolution.
Definition at line 1819 of file wseries.cc.
param: n - number of steps (-1 means full reconstruction)
Definition at line 273 of file wseries.cc.
|
inline |
Definition at line 190 of file wseries.hh.
Definition at line 149 of file wseries.cc.
|
virtual |
Reimplemented from wavearray< DataType_t >.
Definition at line 1108 of file wseries.cc.
double WSeries< DataType_t >::maxEnergy | ( | wavearray< DataType_t > & | ts, |
Wavelet & | w, | ||
double | dT = 0 , |
||
int | N = 1 , |
||
int | pattern = 0 , |
||
TH1F * | hist = NULL |
||
) |
Definition at line 486 of file wseries.cc.
|
inline |
Definition at line 131 of file wseries.hh.
Definition at line 121 of file wseries.hh.
|
virtual |
Definition at line 1091 of file wseries.cc.
Definition at line 117 of file wseries.cc.
|
virtual |
Definition at line 753 of file wseries.cc.
|
virtual |
Reimplemented from wavearray< DataType_t >.
Definition at line 826 of file wseries.cc.
|
virtual |
Reimplemented from wavearray< DataType_t >.
Definition at line 728 of file wseries.cc.
|
virtual |
Definition at line 778 of file wseries.cc.
|
virtual |
Reimplemented from wavearray< DataType_t >.
Definition at line 749 of file wseries.cc.
|
virtual |
Reimplemented from wavearray< DataType_t >.
Definition at line 736 of file wseries.cc.
|
virtual |
Definition at line 802 of file wseries.cc.
|
virtual |
Reimplemented from wavearray< DataType_t >.
Definition at line 745 of file wseries.cc.
|
virtual |
Reimplemented from wavearray< DataType_t >.
Definition at line 732 of file wseries.cc.
WSeries< DataType_t > & WSeries< DataType_t >::operator= | ( | const wavearray< DataType_t > & | a | ) |
Definition at line 684 of file wseries.cc.
WSeries< DataType_t > & WSeries< DataType_t >::operator= | ( | const WSeries< DataType_t > & | a | ) |
Definition at line 698 of file wseries.cc.
WSeries< DataType_t > & WSeries< DataType_t >::operator= | ( | const DataType_t | a | ) |
Definition at line 724 of file wseries.cc.
|
virtual |
Reimplemented from wavearray< DataType_t >.
Definition at line 713 of file wseries.cc.
|
virtual |
param: f - black pixel fraction param: m - mode options: f = 0 - returns black pixel occupancy m = 1 - set threshold f, returns percentile amplitudes m =-1 - set threshold f, returns wavelet amplitudes m > 1 - random policy,returns percentile amplitudes m <-1 - random policy,returns wavelet amplitudes m = 0 - random pixel selection if m<0 return wavelet amplitudes instead of the percentile amplitude
Definition at line 2064 of file wseries.cc.
|
virtual |
param: S - threshold on pixel significance return pixel occupancy.
Definition at line 1965 of file wseries.cc.
void WSeries< DataType_t >::print | ( | ) |
param: int n if n<0, zero pixels defined in mask (regression) if n>=0, zero all pixels except ones defined in the mask param: bool - if true, set WSeries data to be positive if pMask.size()=0, mask(0,true) is equivalent to abs(data) return core pixel occupancy
param: true - core pccupancy, false - total occupancy; return wavearray<double> with occupancy.
Definition at line 2334 of file wseries.cc.
void WSeries< DataType_t >::putLayer | ( | wavearray< DataType_t > & | value, |
double | n | ||
) |
param: n - layer number
Definition at line 201 of file wseries.cc.
|
inline |
Definition at line 177 of file wseries.hh.
Reimplemented from wavearray< DataType_t >.
Definition at line 897 of file wseries.cc.
Reimplemented from wavearray< DataType_t >.
Definition at line 883 of file wseries.cc.
Definition at line 137 of file wseries.hh.
|
virtual |
param: n - sub-interval duration in domain units param: f - black pixel fraction options: f = 0 - returns black pixel occupancy
Definition at line 1572 of file wseries.cc.
|
virtual |
param: T - sliding window duration in seconds param: f - black pixel fraction param: t - sliding step in seconds options: f = 0 - returns black pixel occupancy options: t = 0 - sliding step = wavelet time resolution.
Definition at line 1665 of file wseries.cc.
|
inline |
Definition at line 97 of file wseries.hh.
|
inline |
Definition at line 114 of file wseries.hh.
|
inline |
Definition at line 94 of file wseries.hh.
|
inline |
Definition at line 107 of file wseries.hh.
Definition at line 216 of file wseries.cc.
|
virtual |
param: n - sub-interval duration in seconds param: f - black pixel fraction options: f = 0 - returns black pixel occupancy
Definition at line 1469 of file wseries.cc.
|
inline |
Definition at line 126 of file wseries.hh.
|
virtual |
param: first - time window to calculate normalization constants second - low frequency boundary for correction third - high frequency boundary for correction algorithm: 1) sort wavelet amplitudes with the same time stamp 2) calculate left(p) and right(p) amplitudes put (right(p)-left(p))/2 into output array 3) if first parameter >0 - devide WSeries by average variability
Definition at line 1278 of file wseries.cc.
void WSeries< DataType_t >::wavescan | ( | WSeries< DataType_t > ** | pws, |
int | N, | ||
TH1F * | hist = NULL |
||
) |
Definition at line 604 of file wseries.cc.
double WSeries< DataType_t >::wdmPacket | ( | int | pattern, |
char | opt = 'L' , |
||
TH1F * | hist = NULL |
||
) |
Definition at line 358 of file wseries.cc.
|
virtual |
what it does: each wavelet layer is devided into k intervals.
The data for each interval is sorted and the following parameters are calculated: median and the amplitude corresponding to 31% percentile (wp). Wavelet amplitudes (w) are normalized as w' = (w-median(t))/wp(t), where median(t) and wp(t) is a linear interpolation between (median,wp) measurements for each interval.
Definition at line 1128 of file wseries.cc.
|
virtual |
Definition at line 1186 of file wseries.cc.
|
inline |
Definition at line 102 of file wseries.hh.
|
inline |
Definition at line 104 of file wseries.hh.
|
inline |
Definition at line 128 of file wseries.hh.
double WSeries< DataType_t >::bpp |
Definition at line 442 of file wseries.hh.
double WSeries< DataType_t >::f_high |
Definition at line 448 of file wseries.hh.
double WSeries< DataType_t >::f_low |
Definition at line 446 of file wseries.hh.
Definition at line 438 of file wseries.hh.
size_t WSeries< DataType_t >::w_mode |
Definition at line 440 of file wseries.hh.
double WSeries< DataType_t >::wRate |
Definition at line 444 of file wseries.hh.