add FAR information for the loudest triggers

analyses/STAMPAS/post_processing/background_estimation/background.m

@@ -110,7 +110,8 @@ if config.doZebra
   if sum(strcmp(cut,'Rveto'))
     lcut=[lcut,'Rveto'];
   end
-  get_loudest([searchPath '/figures/loudest.txt'],data,'Cut',cut);
+
+  get_loudest(searchPath,data,'Cut',cut);
 
   %% Rveto DIAGNOSTIC PLOTS
   %------------------------------------------------------------------
@@ -119,9 +120,9 @@ if config.doZebra
     RvetoDiagnostic(data,[searchPath '/figures/'],searchPath);
   end
 
-  %% SNRfrac DIAGNOSTIC PLOTS
-  %------------------------------------------------------------------
-  % only if SNRfrac is used
+%  %% SNRfrac DIAGNOSTIC PLOTS
+%  %------------------------------------------------------------------
+%  % only if SNRfrac is used
   if sum(cell2mat(strfind(cut,'SNRfrac')))
     SNRfracDiagnostic(data,[searchPath '/figures/']);
   end

analyses/STAMPAS/post_processing/background_estimation/src/get_loudest.m

-function [] = get_loudest(File,Data,varargin)  
+function [] = get_loudest(searchPath,Data,varargin)  
+%function [] = get_loudest(File,Data,varargin)  
 % get_loudest :
 % DESCRIPTION :
 % 
@@ -24,6 +25,9 @@ function [] = get_loudest(File,Data,varargin)
   % 
   % parse & check input
   %
+  File = [searchPath  '/figures/loudest.txt'];
+  FAR=load([searchPath  '/figures/FAR_clean.mat']);
+
   p=inputParser;
   addRequired(p,'File');
   addRequired(p,'Data');
@@ -96,7 +100,8 @@ function [] = get_loudest(File,Data,varargin)
   % 
   % get interresting carateristic of the loudNb first trigger
   %
-  loudNb = max(min(sum(S>thr),MaxTriggers),10);
+  loudNb = max(min(sum(S>thr),MaxTriggers),30);
+  F = get_far(S(1:loudNb), FAR);
 
   if doLonetrack
     D=[...
@@ -120,6 +125,7 @@ function [] = get_loudest(File,Data,varargin)
   else
     D=[...
       1:loudNb;...                          % trigger idx
+      F';...                                 % FAR
       S(1:loudNb)';...                      % SNR 
       data(1:loudNb).GPSstart';...          % GPS start ifo1
       data(1:loudNb).GPSstart('ifo2')';...  % GPS start ifo2
@@ -143,7 +149,7 @@ function [] = get_loudest(File,Data,varargin)
 
   % write the loudest triggers property in a text file   
   fid = fopen([folder '/' file '.txt'], 'w');
-  fprintf(fid,['%d %7.2f %12.2f %12.2f %7.2f %7.2f %7.2f %7.2f %7.2f %7.2f ' ...
+  fprintf(fid,['%d %7.2g %7.2f %12.2f %12.2f %7.2f %7.2f %7.2f %7.2f %7.2f %7.2f ' ...
                '%7.4f %3.3f %3.3f %3.3f %d %d\n'],D);
   fclose(fid);
   

analyses/STAMPAS/post_processing/background_estimation/src/html/addLoudestTable.m

@@ -50,6 +50,7 @@ function [] = addLoudestTable(doc,file,params)
   % 
   % <tr>
   %   <th rowspan="2">Index</th>
+  %   <th rowspan="2">FAR</th>
   %   <th rowspan="2">SNR</th>
   %   <th rowspan="2">GPS start</th>
   %   <th rowspan="2">Duration</th>
@@ -78,6 +79,9 @@ function [] = addLoudestTable(doc,file,params)
   doc.createNode('th','Index',...
                  'Attribute',struct('rowspan','2'), ...
                  'Parent',tr);
+  doc.createNode('th','FAR',...
+                 'Attribute',struct('rowspan','2'), ...
+                 'Parent',tr);
   doc.createNode('th','SNR',...
                  'Attribute',struct('rowspan','2'), ...
                  'Parent',tr);