7 int cwb_xtalk_low_level=0;
8 if(gSystem->Getenv(
"CWB_XTALK_LOW_LEVEL")==NULL) {
9 cout <<
"cwb_xtalk.C - Error : environment CWB_XTALK_LOW_LEVEL is not defined!!!" << endl;
exit(1);
11 if(
TString(gSystem->Getenv(
"CWB_XTALK_LOW_LEVEL")).IsDigit()) {
12 cwb_xtalk_low_level=
TString(gSystem->Getenv(
"CWB_XTALK_LOW_LEVEL")).Atoi();
14 cout <<
"cwb_xtalk.C - Error : environment CWB_XTALK_LOW_LEVEL is not defined!!!" << endl;
exit(1);
19 int cwb_xtalk_high_level=0;
20 if(gSystem->Getenv(
"CWB_XTALK_HIGH_LEVEL")==NULL) {
21 cout <<
"cwb_xtalk.C - Error : environment CWB_XTALK_HIGH_LEVEL is not defined!!!" << endl;
exit(1);
23 if(
TString(gSystem->Getenv(
"CWB_XTALK_HIGH_LEVEL")).IsDigit()) {
24 cwb_xtalk_high_level=
TString(gSystem->Getenv(
"CWB_XTALK_HIGH_LEVEL")).Atoi();
26 cout <<
"cwb_xtalk.C - Error : environment CWB_XTALK_HIGH_LEVEL is not defined!!!" << endl;
exit(1);
30 int nRes = cwb_xtalk_high_level-cwb_xtalk_low_level+1;
32 cout <<
"cwb_xtalk.C - Error : low res Level must be < high res level" << endl;
exit(1);
35 cout <<
"cwb_xtalk.C - Error : number of max resolutions is : " <<
NMAX_RES << endl;
exit(1);
40 if(gSystem->Getenv(
"CWB_XTALK_INU")!=NULL) {
41 if(
TString(gSystem->Getenv(
"CWB_XTALK_INU")).IsDigit()) {
42 cwb_xtalk_iNu=
TString(gSystem->Getenv(
"CWB_XTALK_INU")).Atoi();
44 cout <<
"cwb_xtalk.C - Error : environment CWB_XTALK_INU is not defined!!!" << endl;
exit(1);
50 if(gSystem->Getenv(
"CWB_XTALK_PRECISION")!=NULL) {
51 if(
TString(gSystem->Getenv(
"CWB_XTALK_HIGH_LEVEL")).IsDigit()) {
52 cwb_xtalk_precision=
TString(gSystem->Getenv(
"CWB_XTALK_PRECISION")).Atoi();
54 cout <<
"cwb_xtalk.C - Error : environment CWB_XTALK_PRECISION is not defined!!!" << endl;
exit(1);
59 cout <<
"cwb_xtalk.C - define wdm ..." << endl;
61 for(
int level=cwb_xtalk_low_level; level<=cwb_xtalk_high_level; level++) {
62 int layers = level>0 ? 1<<level : 0;
67 cout <<
"cwb_xtalk.C - generate xtalk catalog : be patient, it takes a while ..." << endl;
71 cout <<
"cwb_xtalk.C - write xtalk catalog ..." << endl;
73 sprintf(fName,
"OverlapCatalog-ilLev%d-hLev%d-iNu%d-P%d.xbin",
74 cwb_xtalk_low_level,cwb_xtalk_high_level,cwb_xtalk_iNu,cwb_xtalk_precision);
76 cout <<
"OverlapCatalog Name : " << fName << endl;
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
cout<< "cwb_xtalk.C - define wdm ..."<< endl;WDM< double > * wdm[NMAX_RES]
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)