17 sprintf(name,
"%s/wdmXTalk/OverlapCatalog_Lev_8_16_32_64_128_256_iNu_4_Prec_10.bin",
18 gSystem->Getenv(
"HOME_WAT_FILTERS") );
24 for(
int i=0; i<nRES; i++) wdm[i] = new WDM<double>(layers[
i],layers[
i], 4, 8);
31 double Decay2 = 1.5*1.5;
38 double ii =
i/Rate - Duration/2;
39 ts[
i] = 10*rnd0.Gaus()*exp(-ii*ii/2/Decay2);
48 if(r2>r1){
printf(
"Error: r2 must not be greater than r1!\n");
return;}
53 int size1 = pTF[r1].
size()/2;
54 double* map00_r1 = pTF[r1].
data;
59 int size2 = pTF[
r2].
size()/2;
60 double* map00_r2 = subTF.
data;
64 int k = (M1/
M2)*(M2+1);
66 for(
int m=0;
m<=
M1; ++
m){
68 for(
int i=
m;
i<size1;
i+=M1+1){
73 if(
m==0 ||
m==M1)
if(!odd)
continue;
77 if(indx<0 || indx>size2)
continue;
78 map00_r2[indx] -= map00_r1[
i]*tmp.
data[
j].
CC[0];
83 for(
int i=0;
i<size2;
i+=M2+1)map00_r2[
i] = map00_r2[
i+M2] = 0;
virtual void rate(double r)
WDM< double > wdm(nLAYERS, nLAYERS, 6, 10)
void testMonster(int r1, int r2)
virtual size_t size() const
WSeries< double > pTF[nRES]
printf("total live time: non-zero lags = %10.1f \, liveTot)
void read(char *filename)
param: file name
xtalkArray(*** catalog)[2]
wavearray< double > ts(N)
sprintf(tfres,"(1/%g)x(%g) (sec)x(Hz)", 2 *df, df)