diff --git a/code/RDGW150914_ptemcee.py b/code/RDGW150914_ptemcee.py
index e27eae7ba8d2f1ab065115f0768194bd1ceca2a8..58cb04a0a222d836e757f0aaae80416b3d41ae93 100755
--- a/code/RDGW150914_ptemcee.py
+++ b/code/RDGW150914_ptemcee.py
@@ -16,6 +16,8 @@ from matplotlib import rc
 plt.rcParams.update({'font.size': 16.5})
 
 import ptemcee
+from pycbc.pool import choose_pool
+#from multiprocessing import Pool
 import math
 import h5py
 import inspect
@@ -32,18 +34,19 @@ from scipy.optimize import minimize
 #tshift: time shift after the strain peak
 #vary_fund: whether you vary the fundamental frequency. Works in the model_dv function.
 
-rootpath="/work/rayne.liu" #"/Users/RayneLiu"
+rootpath="/Users/RayneLiu"# "/work/rayne.liu" 
 nmax=1
 tshift=19
 vary_fund = True
 
 #sampler parameters
 npoints=101 
-nwalkers = 32
+nwalkers = 42
 ntemps=12
 ndim = int(4*(nmax+1))
 burnin = 10 #How many points do you burn before doing the corner plot. You need to watch the convergence of the chain plot a bit.
             #This is trivial but often forgotten: this cannot be more than npoints! Usually 1/5~1/4 npoints is what I observe.
+
 numbins = 21 #corner plot parameter - how many bins you want
 datacolor = '#105670' #'#4fa3a7'
 pkcolor = '#f2c977' #'#ffb45f'
@@ -160,8 +163,8 @@ def log_likelihood(theta):
 # The evidence is just a normalization factor
 def log_probability(theta):
     lp = log_prior(theta)
-    print('lp:')
-    print(lp)
+    #print('lp:')
+    #print(lp)
     if not np.isfinite(lp):
         return -np.inf
     return lp + log_likelihood(theta)
@@ -170,6 +173,9 @@ def log_probability(theta):
 
 
 #Fit with ptemcee
+#Set the number of cores of your processors
+pool = choose_pool(4)
+pool.size = 4
 vary_param = float(vary_fund)
 pos = np.array([[random.uniform(-0.1,0.1), random.uniform(-0.1,0.1), 4.28313743e-01, random.uniform(2.5, 2.6) + (1-vary_param) * np.pi]])
 for i in range (1,nmax+1):
@@ -180,7 +186,8 @@ pos = pos.T.flatten()
 pos = list(pos)
 pos += 1e-5 * np.random.randn(ntemps, nwalkers, ndim)
 #print(pos)
-sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps)
+
+sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps, pool=pool)
 sampler.run_mcmc(pos,npoints)
 
 
@@ -291,13 +298,3 @@ plt.xlabel("t")
 plt.ylabel("h")
 
 figband.savefig(rootpath+'/git/rdstackingproject/plotsmc/vary'+str(vary_fund)+'nmax='+str(nmax)+'_tshift='+str(tshift)+'_'+str(npoints)+'pt_band.pdf', format = 'pdf')
-"""
-
-
-
-
-
-
-
-
-"""
diff --git a/code/RDGW150914_ptemcee1.py b/code/RDGW150914_ptemcee1.py
index bde6934d9ee5018e967755d236e84894aefc73c6..e491c3d0c40e6ee34fec8e552dda7018d6070af4 100755
--- a/code/RDGW150914_ptemcee1.py
+++ b/code/RDGW150914_ptemcee1.py
@@ -16,6 +16,7 @@ from matplotlib import rc
 plt.rcParams.update({'font.size': 16.5})
 
 import ptemcee
+from pycbc.pool import choose_pool
 import math
 import h5py
 import inspect
@@ -170,6 +171,9 @@ def log_probability(theta):
 
 
 #Fit with ptemcee
+#Set the number of cores of your processors
+pool = choose_pool(4)
+pool.size = 4
 vary_param = float(vary_fund)
 pos = np.array([[random.uniform(-0.1,0.1), random.uniform(-0.1,0.1), 4.28313743e-01, random.uniform(2.5, 2.6) + (1-vary_param) * np.pi]])
 for i in range (1,nmax+1):
@@ -180,12 +184,11 @@ pos = pos.T.flatten()
 pos = list(pos)
 pos += 1e-5 * np.random.randn(ntemps, nwalkers, ndim)
 #print(pos)
-sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps)
+sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps, pool=pool)
 sampler.run_mcmc(pos,npoints)
 
 
 
-
 #Define labels and start plotting
 paramlabels_a = [r'$\alpha_'+str(i)+'$' for i in range (nmax+1)]
 paramlabels_b = [r'$\beta_'+str(i)+'$' for i in range (nmax+1)]
diff --git a/code/RDGW150914_ptemcee2.py b/code/RDGW150914_ptemcee2.py
index 3d29a783e399d052c83f276483a5f2ce7112af23..19f8eea566cd166e7118de15e4896eb1b8b0625f 100755
--- a/code/RDGW150914_ptemcee2.py
+++ b/code/RDGW150914_ptemcee2.py
@@ -16,6 +16,7 @@ from matplotlib import rc
 plt.rcParams.update({'font.size': 16.5})
 
 import ptemcee
+from pycbc.pool import choose_pool
 import math
 import h5py
 import inspect
@@ -170,6 +171,9 @@ def log_probability(theta):
 
 
 #Fit with ptemcee
+#Set the number of cores of your processors
+pool = choose_pool(4)
+pool.size = 4
 vary_param = float(vary_fund)
 pos = np.array([[random.uniform(-0.1,0.1), random.uniform(-0.1,0.1), 4.28313743e-01, random.uniform(2.5, 2.6) + (1-vary_param) * np.pi]])
 for i in range (1,nmax+1):
@@ -180,7 +184,7 @@ pos = pos.T.flatten()
 pos = list(pos)
 pos += 1e-5 * np.random.randn(ntemps, nwalkers, ndim)
 #print(pos)
-sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps)
+sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps, pool=pool)
 sampler.run_mcmc(pos,npoints)
 
 
diff --git a/code/RDGW150914_ptemcee3.py b/code/RDGW150914_ptemcee3.py
index 6e2f09697d26434e6b3a4fafddfd886f7a481c30..c49e4c2cbe8b1a37f0618b09546e4145576ba8e2 100755
--- a/code/RDGW150914_ptemcee3.py
+++ b/code/RDGW150914_ptemcee3.py
@@ -16,6 +16,7 @@ from matplotlib import rc
 plt.rcParams.update({'font.size': 16.5})
 
 import ptemcee
+from pycbc.pool import choose_pool
 import math
 import h5py
 import inspect
@@ -170,6 +171,9 @@ def log_probability(theta):
 
 
 #Fit with ptemcee
+#Set the number of cores of your processors
+pool = choose_pool(4)
+pool.size = 4
 vary_param = float(vary_fund)
 pos = np.array([[random.uniform(-0.1,0.1), random.uniform(-0.1,0.1), 4.28313743e-01, random.uniform(2.5, 2.6) + (1-vary_param) * np.pi]])
 for i in range (1,nmax+1):
@@ -180,7 +184,7 @@ pos = pos.T.flatten()
 pos = list(pos)
 pos += 1e-5 * np.random.randn(ntemps, nwalkers, ndim)
 #print(pos)
-sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps)
+sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps, pool=pool)
 sampler.run_mcmc(pos,npoints)
 
 
diff --git a/code/RDGW150914_ptemcee4.py b/code/RDGW150914_ptemcee4.py
index 759575aab9213c5c90a1d69d5fd359056fc0f697..af162ee13be5c64cff10672205638284cdf66546 100755
--- a/code/RDGW150914_ptemcee4.py
+++ b/code/RDGW150914_ptemcee4.py
@@ -16,6 +16,7 @@ from matplotlib import rc
 plt.rcParams.update({'font.size': 16.5})
 
 import ptemcee
+from pycbc.pool import choose_pool
 import math
 import h5py
 import inspect
@@ -170,6 +171,9 @@ def log_probability(theta):
 
 
 #Fit with ptemcee
+#Set the number of cores of your processors
+pool = choose_pool(4)
+pool.size = 4
 vary_param = float(vary_fund)
 pos = np.array([[random.uniform(-0.1,0.1), random.uniform(-0.1,0.1), 4.28313743e-01, random.uniform(2.5, 2.6) + (1-vary_param) * np.pi]])
 for i in range (1,nmax+1):
@@ -180,7 +184,7 @@ pos = pos.T.flatten()
 pos = list(pos)
 pos += 1e-5 * np.random.randn(ntemps, nwalkers, ndim)
 #print(pos)
-sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps)
+sampler = ptemcee.Sampler(nwalkers, ndim, log_likelihood, log_prior, ntemps=ntemps, pool=pool)
 sampler.run_mcmc(pos,npoints)
 
 
diff --git a/code/condor_submit_RdownPtemcee.sub b/code/condor_submit_RdownPtemcee.sub
index ced2df3f21060872efeeb4016a059b0c9ffc257d..478b6f84494b58c0efa2e0e476b34c1b5a68f4c2 100755
--- a/code/condor_submit_RdownPtemcee.sub
+++ b/code/condor_submit_RdownPtemcee.sub
@@ -12,8 +12,8 @@ initialdir = .
 notify_user = rl746@cornell.edu
 notification = Complete
 arguments  = "-processid $(Process)" 
-request_memory = 16GB
-request_cpus = 1
+request_memory = 8GB
+request_cpus = 4
 on_exit_remove = (ExitBySignal == False) || ((ExitBySignal == True) && (ExitSignal != 11))
 accounting_group = aei.dev.test_dynesty
 queue 1
diff --git a/code/condor_submit_RdownPtemcee1.sub b/code/condor_submit_RdownPtemcee1.sub
index 34bcd92279d39c84d9542a061962702e76bbb74d..0358c8491dac5a8a9b425b2264e62d46fd871752 100755
--- a/code/condor_submit_RdownPtemcee1.sub
+++ b/code/condor_submit_RdownPtemcee1.sub
@@ -12,8 +12,8 @@ initialdir = .
 notify_user = rl746@cornell.edu
 notification = Complete
 arguments  = "-processid $(Process)" 
-request_memory = 16GB
-request_cpus = 1
+request_memory = 8GB
+request_cpus = 4
 on_exit_remove = (ExitBySignal == False) || ((ExitBySignal == True) && (ExitSignal != 11))
 accounting_group = aei.dev.test_dynesty
 queue 1
diff --git a/code/condor_submit_RdownPtemcee2.sub b/code/condor_submit_RdownPtemcee2.sub
index 218ef23caa5ac19eb5080fa3390a33bbc16e866f..f780e1f37082973a0be0d329fa06226206523816 100755
--- a/code/condor_submit_RdownPtemcee2.sub
+++ b/code/condor_submit_RdownPtemcee2.sub
@@ -12,8 +12,8 @@ initialdir = .
 notify_user = rl746@cornell.edu
 notification = Complete
 arguments  = "-processid $(Process)" 
-request_memory = 16GB
-request_cpus = 1
+request_memory = 8GB
+request_cpus = 4
 on_exit_remove = (ExitBySignal == False) || ((ExitBySignal == True) && (ExitSignal != 11))
 accounting_group = aei.dev.test_dynesty
 queue 1
diff --git a/code/condor_submit_RdownPtemcee3.sub b/code/condor_submit_RdownPtemcee3.sub
index de8eda0a032e77ef19e3468474c1507def13383d..f9e3d38bf94434f840bc9e10c90ab0c7da1c8488 100755
--- a/code/condor_submit_RdownPtemcee3.sub
+++ b/code/condor_submit_RdownPtemcee3.sub
@@ -12,8 +12,8 @@ initialdir = .
 notify_user = rl746@cornell.edu
 notification = Complete
 arguments  = "-processid $(Process)" 
-request_memory = 16GB
-request_cpus = 1
+request_memory = 8GB
+request_cpus = 4
 on_exit_remove = (ExitBySignal == False) || ((ExitBySignal == True) && (ExitSignal != 11))
 accounting_group = aei.dev.test_dynesty
 queue 1
diff --git a/code/condor_submit_RdownPtemcee4.sub b/code/condor_submit_RdownPtemcee4.sub
index b337b63fcba33f686be5fe6176dc1f1ddd65b580..2780f90f0c714bcd310e17fbeb16414abe6dbf93 100755
--- a/code/condor_submit_RdownPtemcee4.sub
+++ b/code/condor_submit_RdownPtemcee4.sub
@@ -12,8 +12,8 @@ initialdir = .
 notify_user = rl746@cornell.edu
 notification = Complete
 arguments  = "-processid $(Process)" 
-request_memory = 16GB
-request_cpus = 1
+request_memory = 8GB
+request_cpus = 4
 on_exit_remove = (ExitBySignal == False) || ((ExitBySignal == True) && (ExitSignal != 11))
 accounting_group = aei.dev.test_dynesty
 queue 1