# -*- coding: utf-8 -*-
# @Author: Elena Graverini
# @Date:   2017-03-16 15:22:00
# @Last Modified by:   Elena Graverini
# @Last Modified time: 2017-03-17 16:42:38
import os
# import subprocess

chunk_size = 100000
chunks = 100
leptons = ['mu', 'tau']
commands = []
seed = 12345
form_factors = "BBGIOvD2017"
signal_pdf = "Lambda_b->Lambda_c(2625)lnu::dGamma"
model = "SM"
prerun_min = 500000
prerun_update = 2000
kin_limits = {
    'mu': (0.011163612964000001, 8.947815864099999),
    'tau': (3.1570893124, 8.947815864099999),
}

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,
                     ])

for c in commands:
    print(' '.join(c))
    # subprocess.call(c)
    os.system(' '.join(c))