diff --git a/bucoffea/vbfhinv/definitions.py b/bucoffea/vbfhinv/definitions.py index 9726156a..5f84aff9 100644 --- a/bucoffea/vbfhinv/definitions.py +++ b/bucoffea/vbfhinv/definitions.py @@ -376,21 +376,29 @@ def add_lists(orig_list, to_add): ] # Dimuon CR - cr_2m_cuts = ['trig_met','two_muons', 'at_least_one_tight_mu', 'dimuon_mass', 'veto_ele', 'dimuon_charge'] + common_cuts[1:] + ['dpfcalo_cr'] - + cr_2m_cuts = ['trig_met', 'at_least_one_tight_mu', 'two_muons', 'dimuon_charge', 'dimuon_mass'] + common_cuts + ['dpfcalo_cr'] cr_2m_cuts.remove('veto_muo') - regions['cr_2m_vbf'] = cr_2m_cuts # Single muon CR - cr_1m_cuts = ['trig_met','one_muon', 'at_least_one_tight_mu', 'veto_ele'] + common_cuts[1:] + ['dpfcalo_cr'] + cr_1m_cuts = ['trig_met', 'at_least_one_tight_mu', 'one_muon'] + common_cuts + ['dpfcalo_cr'] cr_1m_cuts.remove('veto_muo') regions['cr_1m_vbf'] = cr_1m_cuts # Dielectron CR - cr_2e_cuts = ['trig_ele','two_electrons', 'at_least_one_tight_el', 'dielectron_mass', 'veto_muo', 'dielectron_charge'] + common_cuts[2:] + ['dpfcalo_cr'] - # cr_2e_cuts.remove('veto_ele') + cr_2e_cuts = ['trig_ele', 'at_least_one_tight_el', 'two_electrons', 'dielectron_charge', 'dielectron_mass'] + common_cuts + ['dpfcalo_cr'] + cr_2e_cuts.remove('veto_ele') regions['cr_2e_vbf'] = cr_2e_cuts + + # Single electron CR + cr_1e_cuts = ['trig_ele', 'at_least_one_tight_el', 'one_electron', 'met_el'] + common_cuts + ['dpfcalo_cr', 'no_el_in_hem'] + cr_1e_cuts.remove('veto_ele') + regions['cr_1e_vbf'] = cr_1e_cuts + + # Photon CR + cr_g_cuts = ['trig_photon', 'at_least_one_tight_photon', 'one_photon', 'photon_pt'] + common_cuts + ['dpfcalo_cr'] + cr_g_cuts.remove('veto_photon') + regions['cr_g_vbf'] = cr_g_cuts # Z CRs with CaloMETNoLep cut if cfg.RUN.CALOMET_CHECK: @@ -398,16 +406,6 @@ def add_lists(orig_list, to_add): regions[f'{r}_calocut'] = copy.deepcopy(regions[r]) regions[f'{r}_calocut'].append('calo_metptnolep') - # Single electron CR - cr_1e_cuts = ['trig_ele','one_electron', 'at_least_one_tight_el', 'veto_muo','met_el'] + common_cuts[1:] + ['dpfcalo_cr', 'no_el_in_hem'] - # cr_1e_cuts.remove('veto_ele') - regions['cr_1e_vbf'] = cr_1e_cuts - - # Photon CR - cr_g_cuts = ['trig_photon', 'one_photon', 'at_least_one_tight_photon','photon_pt'] + common_cuts + ['dpfcalo_cr'] - cr_g_cuts.remove('veto_photon') - - regions['cr_g_vbf'] = cr_g_cuts if cfg and cfg.RUN.SYNC: regions['sync_sr_vbf_round1'] = [ @@ -510,6 +508,8 @@ def add_lists(orig_list, to_add): regions[f'tr_g_{trgname}_photon_pt_trig_cut_num'] = tr_g_num_cuts + [trgname, 'photon_pt_trig'] regions[f'tr_g_{trgname}_photon_pt_trig_cut_den'] = tr_g_den_cuts + [trgname, 'photon_pt_trig'] + for reg,cuts in regions.items(): + print(reg, cuts) return regions def ak4_em_frac_weights(weights, diak4, evaluator):