update from my local mac

plat.txt

@@ -47,5 +47,5 @@ lh	123	140
 ss	123	140
 sh	123	140
 hh	123	155
-t0 	40	160
+t0 	25	170
 #end

src/analysis.jl

@@ -37,8 +37,8 @@ using Revise, juobs, PyPlot, LaTeXStrings, ADerrors, DelimitedFiles
 const path_data = "/Users/ale/Desktop/data"
 const path_plat = "/Users/ale/automation/plat.txt"
 const path_results = "/Users/ale/Desktop/results"
-const ensembles = ["J303", "H400", "N300"]
-const sector = Dict("ll"=>false, "ls"=>false, "lh"=>true, "ss"=>false, "sh"=>true, "hh"=>true )
+const ensembles = ["H400", "J303"]
+const sector = Dict("ll"=>true, "ls"=>true, "lh"=>true, "ss"=>true, "sh"=>true, "hh"=>true )
 const rwf = false
 const compute_t0 = false
 const mass_shift = false #not implemented yet
@@ -64,8 +64,8 @@ end
 #======== ANALYSIS ==========#
 
 println("\n Computing Observables \n")
-
-    
+wpm = Dict{String,Vector{Float64}}()
+wpm["H102"] = [-1.0, 1.0, -1.0, -1.0]
 @time begin
     
     ensobs = Vector{EnsObs}(undef,length(ensembles))
@@ -77,8 +77,8 @@ println("\n Computing Observables \n")
             a1a1 = get_corr(ens, "G1G5", "G1G5", rw=rwf)
             mu_list = getfield.(pp, :mu)
             m_ps = get_meff(pp, ens)
-            m_vec = get_meff(a1a1, ens)
-            f_ps = get_f(pp, m_ps, ens)
+            m_vec = get_meff(a1a1, ens, wpm=wpm)
+            f_ps = get_f(pp, m_ps, ens, wpm=wpm)
             f_vec = get_f(a1a1, m_vec, ens)
             ensobs[i] =  EnsObs(ens, mu_list, m_ps, f_ps, m_vec, f_vec)
             if ens.deg
@@ -100,8 +100,7 @@ println("\n Computing Observables \n")
     end
     println("Total time:")
 end
-    
-
+##  
 #=========== MATCHING =============#
 
 println("\n Matching Procedure \n")