adding \eta_c and J/\psi masses in juobs_const.jl

constants/juobs_const.jl

 using ADerrors
 #PDG
 const hc = 197.3269804 #MeV fm
-const M_values = [1869.65, 2010.26, 1968.34, 2112.2] #MD, MD*, MDs, MDs* (MeV)
-const M_error = [0.05, 0.05, 0.07, 0.4]
+const M_values = [1869.65, 2010.26, 1968.34, 2112.2, 2980.3, 3096.916] #MD, MD*, MDs, MDs*, \eta_c, J/\psi (MeV)
+const M_error = [0.05, 0.05, 0.07, 0.4, 1.2, 0.011]
 #1802.05243 
 const b_values = [3.40, 3.46, 3.55, 3.70, 3.85]
 const ZM_data = [1.9684, 1.9935, 2.0253, 2.0630, 2.0814]
@@ -20,13 +20,15 @@ const t0_ph_error = ones(1,1) .* 5e-3
 const C1 = zeros(5, 5)
 const C2 = zeros(5, 5)
 const C3 = zeros(4, 4)
-const C4 = zeros(4, 4)
-for i = 1:5
+const C4 = zeros(6, 6)
+for i = 1:6
+    C4[i, i] = M_error[i] ^ 2
+    if i<=5
         C1[i, i] = ZM_error[i] ^ 2
         C2[i, i] = ZM_tm_error[i] ^ 2
         if i <= 4
             C3[i, i] = t0_error[i] ^ 2
-        C4[i, i] = M_error[i] ^ 2
+        end
     end
 
 end