loudest list now excludes triggers vetoed by SNRfrac and FAR computed with background FAR

analyses/STAMPAS/post_processing/background_estimation/openbox.m

@@ -42,7 +42,6 @@ end
 %% INIT
 %------------------------------------------------------------------
 % 
-
 if ~exist([searchPath '/figures'])
   system(['mkdir -p ' searchPath '/figures']);
 end
@@ -116,7 +115,7 @@ end
 %% LOUDEST
 %------------------------------------------------------------------
 % comput the loudest cut : cut apply in the list of loudest
-% only TFveto & Rveto if exist
+% only TFveto & Rveto & SNRfrac if exist
 %
 
 lcut={'noCut'};
@@ -126,6 +125,9 @@ end
 if sum(strcmp(cut,'Rveto'))
   lcut=[lcut,'Rveto'];
 end
+if ~isempty(SNRfracCut)
+  lcut=[lcut, ['SNRfrac:' num2str(SNRfracCut)]];
+end
 
 if config.doLonetrack
   jobs=load([config.searchPath '/tmp/joblistTrue.txt']);
@@ -136,7 +138,8 @@ if config.doLonetrack
                            load([bkgPath '/figures/pproc_FAR_clean.mat']),...
                            ZLTotTime);
 else
-  get_loudest([config.searchPath '/figures/loudest.txt'],data,'Cut',lcut);
+  far=load([bkgPath '/figures/FAR_clean.mat']);
+  get_loudest_ZL(config.searchPath,data,far,'Cut',lcut);
 end
 
 if config.doZebra
@@ -161,9 +164,9 @@ if config.doZebra
   %------------------------------------------------------------------
   % only if Rveto is used
   if sum(cell2mat(strfind(cut,'Rveto')))
-  RvetoDiagnostic(data,[config.searchPath '/figures/'],searchPath);
+    RvetoDiagnostic(data,[config.searchPath '/figures/'],searchPath);
   end
-  
+ 
   %% SNRfrac DIAGNOSTIC PLOTS
   %------------------------------------------------------------------
   % only if SNRfrac is used