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: