diff --git a/eos-toys/make_toys.py b/eos-toys/make_toys.py
index 1df02ed..d44f987 100644
--- a/eos-toys/make_toys.py
+++ b/eos-toys/make_toys.py
@@ -2,7 +2,7 @@
 # @Author: Elena Graverini
 # @Date:   2017-03-16 15:22:00
 # @Last Modified by:   Elena Graverini
-# @Last Modified time: 2017-03-17 16:42:38
+# @Last Modified time: 2017-04-05 08:24:09
 import os
 # import subprocess
 
@@ -11,7 +11,7 @@
 leptons = ['mu', 'tau']
 commands = []
 seed = 12345
-form_factors = "BBGIOvD2017"
+form_factors = "BBGIOvD2017-linear"
 signal_pdf = "Lambda_b->Lambda_c(2625)lnu::dGamma"
 model = "SM"
 prerun_min = 500000
@@ -23,17 +23,23 @@
 
 commands = []
 for l in leptons:
-    commands.append(['eos-sample-events-mcmc',
-                     # '  --kinematics \"s\" %s %s' % kin_limits[l],
-                     '  --global-option \"form-factors\" \"%s\"' % form_factors,
-                     '  --signal-pdf \"%s,l=%s,model=%s\"' % (signal_pdf, l, model),
-                     '  --seed %s' % seed,
-                     '  --output mcmc_samples-%s.hdf5' % l,
-                     '  --prerun-min %s' % prerun_min,
-                     '  --prerun-update %s' % prerun_update,
-                     '  --chunk-size %s' % chunk_size,
-                     '  --chunks %s' % chunks,
-                     ])
+    rho = {
+        -0.05: 'linear/mcmc_samples-rho-0.05-%s.hdf5' % l,
+        0.00: 'linear/mcmc_samples-rho0.00-%s.hdf5' % l,
+        0.05: 'linear/mcmc_samples-rho0.05-%s.hdf5' % l,
+    }
+    for r in rho.keys():
+        commands.append(['eos-sample-events-mcmc',
+                         '  --fix \"Lambda_b->Lambda_c(2625)::rho@BBGIOvD2017-linear\" %s' % r,
+                         '  --global-option \"form-factors\" \"%s\"' % form_factors,
+                         '  --signal-pdf \"%s,l=%s,model=%s\"' % (signal_pdf, l, model),
+                         '  --seed %s' % seed,
+                         '  --output %s' % rho[r],
+                         '  --prerun-min %s' % prerun_min,
+                         '  --prerun-update %s' % prerun_update,
+                         '  --chunk-size %s' % chunk_size,
+                         '  --chunks %s' % chunks,
+                         ])
 
 for c in commands:
     print(' '.join(c))