(Temporary?) fix to buf in BayesWavePost.c; printing overlaps to file.

git-svn-id: https://svn.ligo.caltech.edu/svn/bayeswave/trunk@110 c56465c9-8126-4a4f-9d7d-ac845eff4865

src/BayesWavePost.c

@@ -27,7 +27,7 @@ double innerproduct(int imin, int imax, double *a, double *b, double Sn, double
 int main(int argc, char *argv[])
 {
 
-	/*   Variable declaration   */
+  /*   Variable declaration   */
   int i, d, ifo, imin;//, imax;
 
   char filename[100];
@@ -102,7 +102,7 @@ int main(int argc, char *argv[])
   }
   NI=ifo;
 
-	double **freqData = dmatrix(0,NI-1,0,N-1);
+  double **freqData = dmatrix(0,NI-1,0,N-1);
   double **psd      = dmatrix(0,NI-1,0,N/2-1);
 
   struct Data *data = malloc(sizeof(struct Data));
@@ -144,7 +144,6 @@ int main(int argc, char *argv[])
 
   double TwoDeltaToverN = 2.0 * deltaT / ((double) dataPtr->timeData->data->length);
 
-
   //Copy Fourier domain data & PSD from LAL data types into simple arrays
   ifo=0;
   imin = (int)(fmin*Tobs);
@@ -177,7 +176,6 @@ int main(int argc, char *argv[])
   /*
    DATA
    */
-
   initialize_data(data,freqData,N,tsize,psd,Tobs,NI,fmin,fmax);
 
   data->detector  = malloc(NI*sizeof(LALDetector*));
@@ -215,9 +213,15 @@ int main(int argc, char *argv[])
   /******************************************************************************/
 
   //File containing signal-model waveform parameters
-  sprintf(filename,"chains/%s_1_wavechain.dat.0",data->runName);
+  //Was causing segmentation fault. wavechain filename contains the model, not
+  //the run name. Should this file opening command be made more general?
+  //sprintf(filename,"chains/%s_1_wavechain.dat.0",data->runName);
+  sprintf(filename,"chains/signal_1_wavechain.dat.0");
   FILE *signalParams = fopen(filename,"r");
 
+  //Print inj-rec overlaps to file
+  FILE *overlapFile = fopen("injected_reconstructed_overlaps.dat","w");
+
   struct Wavelet *signal = model->signal;
 
   /*
@@ -254,8 +258,6 @@ int main(int argc, char *argv[])
   //TODO:  Should each post-processing quantity have its own chain file?
   double overlap, rr, ii, ri;
 
-
-
   while(!feof(signalParams))
   {
 
@@ -269,22 +271,19 @@ int main(int argc, char *argv[])
      SIGNAL MODEL
     */
 
-    //Zero out signal->templates array (holds the linear combination of wavelets)
-    for(i=0; i<N; i++) signal->templates[i] = 0.0;
-
-    //Get the number of wavelet basis functions for the current sample
-    fscanf(signalParams, "%i", &signal->size);
+   //Zero out signal->templates array (holds the linear combination of wavelets)
+   for(i=0; i<N; i++) signal->templates[i] = 0.0;
 
+   //Get the number of wavelet basis functions for the current sample
+   fscanf(signalParams, "%i", &signal->size);
 
-    //First in the wavechain file are the extrinsic parameters
+   //First in the wavechain file are the extrinsic parameters
     for(i=0; i<NE; i++)
       fscanf(signalParams,"%lg", &model->extParams[i]);
 
-
     //Compute network projection coefficients
     computeProjectionCoeffs(data, model->projection, model->extParams);
 
-
     //Get each wavelets parameters and fill up signal structure
     for(d=1; d<=signal->size; d++)
     {
@@ -303,7 +302,6 @@ int main(int argc, char *argv[])
     //fill hrec with reconstructed waveform
     for(ifo=0; ifo<NI; ifo++) for(i=0; i<N; i++) hrec[ifo][i] = model->response[ifo][i];
 
-
     /* 
      PSD MODEL
     */
@@ -312,8 +310,6 @@ int main(int argc, char *argv[])
     //Fill psd with current model PSD
     for(ifo=0; ifo<NI; ifo++) for(i=0; i<N/2; i++) psd[ifo][i] = data->SnS[ifo][i]; //data->SnS holds the "smooth" part of the PSD
 
-
-
     /******************************************************************************/
     /*                                                                            */
     /*  Compute various moments/overlaps to quantify waveform reconstruction      */
@@ -340,6 +336,7 @@ int main(int argc, char *argv[])
       }
 
       overlap = ri/sqrt(ii*rr);
+      fprintf(overlapFile, "%lg\n", overlap);
     }
 
   }
@@ -354,7 +351,7 @@ int main(int argc, char *argv[])
 
   free(model);
 
-	return 0;
+  return 0;
 
 }