import theory_errors

for i in [110, 120, 130, 140, 150, 160]:
    print 'mass:', i, 'pdf:', theory_errors.get_pdf_err_str(i, 4), \
            'scale:', theory_errors.get_scale_err_str(i, 4)
    mass = int(options.mass.replace('FF', ''))
elif 'SM4' in options.mass:
    is_sm4 = True
    mass = int(options.mass.replace('SM4', ''))
elif 'WW' in options.mass:
    is_ww_only = True
    mass = int(options.mass.replace('WW', ''))
elif 'tt' in options.mass:
    is_tt_only = True
    mass = int(options.mass.replace('tt', ''))
else:
    mass = int(options.mass)

# Get the theory errors
pdf_err = theory_errors.get_pdf_err_str(mass)
scale_err = theory_errors.get_scale_err_str(mass)

# Figure out what the scale systematics are
scale_systematics = {}
with open('scale_systematics.json') as scale_sys_file:
    scale_systematics = json.load(scale_sys_file)

# Define which histograms are signal histograms
signal_datasets = ['VH%i' % mass, 'VH%iWW' % mass]

# If we are fermiophobic, we only use HWW
if is_fermiophobic:
    signal_datasets = ['VH%iWWFF' % mass]
elif is_sm4:
    signal_datasets = ['VH%iSM4' % mass, 'VH%iWWSM4' % mass]
elif is_tt_only: