Logo Coherent WaveBurst  
Reference Guide
Logo
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Pages
CreateOverlapCatalog.C
Go to the documentation of this file.
1 {
2  #define OVERLAP_CATALOG "OverlapCatalog_Lev_32_64_128_256_512_1024_iNu_4_Prec_10.bin"
3  #define nRES 6
4 
5  // define resolutions of interest:
6  int layers[nRES] = {32, 64, 128, 256, 512, 1024};
7 
8  // create corresponding WDM transforms:
10  for(int i=0; i<nRES; i++)wdm[i] = new WDM<double>(layers[i], layers[i], 4, 10);
11 
12  // create the catalog:
13  monster x(wdm, nRES);
14 
15  // write the catalog in a file for future use:
17 }
#define OVERLAP_CATALOG
cout<< "skymap size : "<< L<< endl;for(int l=0;l< L;l++) sm.set(l, l);sm > const_cast< char * >("skymap.dat")
WDM< double > * wdm[nRES]
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
int layers
#define nRES
i drho i