def filter_electron(el_pt=' > 25 ', do_cutflow=False, do_hists=False, apply_corrections=False, evalPID='medium', invertIso=False): filt = Filter('FilterElectron') filt.cut_pt = el_pt filt.cut_eta = ' < 2.5' filt.cut_abssceta = ' <2.5 ' #filt.cut_tight = ' == True ' #filt.cut_medium = ' == True ' #filt.cut_vid_tight = ' == True ' filt.cut_vid_medium = ' == True ' filt.cut_muon_dr = ' > 0.4 ' filt.add_var('triggerMatchBits', '58,109') filt.cut_d0_barrel = ' < 0.05 ' filt.cut_d0_endcap = ' < 0.10 ' filt.cut_dz_barrel = ' < 0.10 ' filt.cut_dz_endcap = ' < 0.20 ' filt.cut_hovere_94x = ' == True' filt.cut_isorho_94x = ' == True' ### 94X V2 PID cuts ### filt.cut_sigmaIEIE_barrel_tight = ' < 0.0104 ' filt.cut_absdEtaIn_barrel_tight = ' < 0.00255 ' filt.cut_absdPhiIn_barrel_tight = ' < 0.022 ' filt.cut_hovere_barrel_tight = ' < 0.026' filt.cut_isoRho_barrel_tight = ' < 0.0287 ' filt.cut_ooEmooP_barrel_tight = ' < 0.0159 ' filt.cut_misshits_barrel_tight = ' < 2 ' filt.cut_passConvVeto_barrel_tight = ' == 1 ' filt.cut_sigmaIEIE_barrel_medium = ' < 0.0106 ' filt.cut_absdEtaIn_barrel_medium = ' < 0.0032 ' filt.cut_absdPhiIn_barrel_medium = ' < 0.0547 ' filt.cut_hovere_barrel_medium = ' < 0.046 ' filt.cut_isoRho_barrel_medium = ' < 0.0478 ' filt.cut_ooEmooP_barrel_medium = ' < 0.184 ' filt.cut_misshits_barrel_medium = ' < 2 ' filt.cut_passConvVeto_barrel_medium = ' == 1 ' filt.cut_sigmaIEIE_barrel_loose = ' < 0.0112 ' filt.cut_absdEtaIn_barrel_loose = ' < 0.00377 ' filt.cut_absdPhiIn_barrel_loose = ' < 0.0884 ' filt.cut_hovere_barrel_loose = ' < 0.05 ' filt.cut_isoRho_barrel_loose = ' < 0.112 ' filt.cut_ooEmooP_barrel_loose = ' < 0.193 ' filt.cut_misshits_barrel_loose = ' < 2 ' filt.cut_passConvVeto_barrel_loose = ' == 1 ' filt.cut_sigmaIEIE_barrel_veryloose = ' < 0.0126 ' filt.cut_absdEtaIn_barrel_veryloose = ' < 0.00463 ' filt.cut_absdPhiIn_barrel_veryloose = ' < 0.148 ' filt.cut_hovere_barrel_veryloose = ' < 0.05 ' filt.cut_isoRho_barrel_veryloose = ' < 0.198 ' filt.cut_ooEmooP_barrel_veryloose = ' < 0.209 ' filt.cut_misshits_barrel_veryloose = ' < 3 ' filt.cut_passConvVeto_barrel_veryloose = ' == 1 ' filt.cut_sigmaIEIE_endcap_tight = ' < 0.0353 ' filt.cut_absdEtaIn_endcap_tight = ' < 0.00501 ' filt.cut_absdPhiIn_endcap_tight = ' < 0.0236 ' filt.cut_hovere_endcap_tight = ' < 0.0188' filt.cut_isoRho_endcap_tight = ' < 0.0445 ' filt.cut_ooEmooP_endcap_tight = ' < 0.0197 ' filt.cut_misshits_endcap_tight = ' < 2 ' filt.cut_passConvVeto_endcap_tight = ' == 1 ' filt.cut_sigmaIEIE_endcap_medium = ' < 0.0387 ' filt.cut_absdEtaIn_endcap_medium = ' < 0.00632 ' filt.cut_absdPhiIn_endcap_medium = ' < 0.0394 ' filt.cut_hovere_endcap_medium = ' < 0.0275 ' filt.cut_isoRho_endcap_medium = ' < 0.0658 ' filt.cut_ooEmooP_endcap_medium = ' < 0.0721 ' filt.cut_misshits_endcap_medium = ' < 2 ' filt.cut_passConvVeto_endcap_medium = ' == 1 ' filt.cut_sigmaIEIE_endcap_loose = ' < 0.0425 ' filt.cut_absdEtaIn_endcap_loose = ' < 0.00674 ' filt.cut_absdPhiIn_endcap_loose = ' < 0.169 ' filt.cut_hovere_endcap_loose = ' < 0.0441 ' filt.cut_isoRho_endcap_loose = ' < 0.108 ' filt.cut_ooEmooP_endcap_loose = ' < 0.111 ' filt.cut_misshits_endcap_loose = ' < 2 ' filt.cut_passConvVeto_endcap_loose = ' == 1 ' filt.cut_sigmaIEIE_endcap_veryloose = ' < 0.0457 ' filt.cut_absdEtaIn_endcap_veryloose = ' < 0.00814 ' filt.cut_absdPhiIn_endcap_veryloose = ' < 0.19 ' filt.cut_hovere_endcap_veryloose = ' < 0.05 ' filt.cut_isoRho_endcap_veryloose = ' < 0.203 ' filt.cut_ooEmooP_endcap_veryloose = ' < 0.132 ' filt.cut_misshits_endcap_veryloose = ' < 4 ' filt.cut_passConvVeto_endcap_veryloose = ' == 1 ' if invertIso: for key, var in vars(filt).items(): if 'iso' in key: print('Inverting iso', key, var) filt.invert(key) # ### 80X VID cuts ### # filt.cut_sigmaIEIE_barrel_tight = ' < 0.00998 ' # filt.cut_absdEtaIn_barrel_tight = ' < 0.00308 ' # filt.cut_absdPhiIn_barrel_tight = ' < 0.0816 ' # filt.cut_hovere_barrel_tight = ' < 0.0414 ' # filt.cut_isoRho_barrel_tight = ' < 0.0588 ' # filt.cut_ooEmooP_barrel_tight = ' < 0.0129 ' # #filt.cut_d0_barrel_tight = ' < 0.0111 ' # #filt.cut_z0_barrel_tight = ' < 0.0466 ' # filt.cut_misshits_barrel_tight = ' < 2 ' # filt.cut_passConvVeto_barrel_tight = ' == 1 ' # # filt.cut_sigmaIEIE_barrel_medium = ' < 0.00998 ' # filt.cut_absdEtaIn_barrel_medium = ' < 0.00311 ' # filt.cut_absdPhiIn_barrel_medium = ' < 0.103 ' # filt.cut_hovere_barrel_medium = ' < 0.253 ' # filt.cut_isoRho_barrel_medium = ' < 0.0695 ' # filt.cut_ooEmooP_barrel_medium = ' < 0.134 ' # #filt.cut_d0_barrel_medium = ' < 0.0118 ' # #filt.cut_z0_barrel_medium = ' < 0.373 ' # filt.cut_misshits_barrel_medium = ' < 2 ' # filt.cut_passConvVeto_barrel_medium = ' == 1 ' # # filt.cut_sigmaIEIE_barrel_loose = ' < 0.011 ' # filt.cut_absdEtaIn_barrel_loose = ' < 0.00477 ' # filt.cut_absdPhiIn_barrel_loose = ' < 0.222 ' # filt.cut_hovere_barrel_loose = ' < 0.298 ' # filt.cut_isoRho_barrel_loose = ' < 0.0994 ' # filt.cut_ooEmooP_barrel_loose = ' < 0.241 ' # #filt.cut_d0_barrel_loose = ' < 0.0261 ' # #filt.cut_z0_barrel_loose = ' < 0.41 ' # filt.cut_misshits_barrel_loose = ' < 2 ' # filt.cut_passConvVeto_barrel_loose = ' == 1 ' # # filt.cut_sigmaIEIE_barrel_veryloose = ' < 0.0115 ' # filt.cut_absdEtaIn_barrel_veryloose = ' < 0.00749 ' # filt.cut_absdPhiIn_barrel_veryloose = ' < 0.228 ' # filt.cut_hovere_barrel_veryloose = ' < 0.356 ' # filt.cut_isoRho_barrel_veryloose = ' < 0.175 ' # filt.cut_ooEmooP_barrel_veryloose = ' < 0.299 ' # #filt.cut_d0_barrel_veryloose = ' < 0.0564 ' # #filt.cut_z0_barrel_veryloose = ' < 0.472 ' # filt.cut_misshits_barrel_veryloose = ' < 3 ' # filt.cut_passConvVeto_barrel_veryloose = ' == 1 ' # # filt.cut_sigmaIEIE_endcap_tight = ' < 0.0292 ' # filt.cut_absdEtaIn_endcap_tight = ' < 0.00605 ' # filt.cut_absdPhiIn_endcap_tight = ' < 0.0394 ' # filt.cut_hovere_endcap_tight = ' < 0.0641 ' # filt.cut_isoRho_endcap_tight = ' < 0.0571 ' # filt.cut_ooEmooP_endcap_tight = ' < 0.0129 ' # #filt.cut_d0_endcap_tight = ' < 0.0351 ' # #filt.cut_z0_endcap_tight = ' < 0.417 ' # filt.cut_misshits_endcap_tight = ' < 2 ' # filt.cut_passConvVeto_endcap_tight = ' == 1 ' # # filt.cut_sigmaIEIE_endcap_medium = ' < 0.0298 ' # filt.cut_absdEtaIn_endcap_medium = ' < 0.00609 ' # filt.cut_absdPhiIn_endcap_medium = ' < 0.045 ' # filt.cut_hovere_endcap_medium = ' < 0.0878 ' # filt.cut_isoRho_endcap_medium = ' < 0.0821 ' # filt.cut_ooEmooP_endcap_medium = ' < 0.13 ' # #filt.cut_d0_endcap_medium = ' < 0.0739 ' # #filt.cut_z0_endcap_medium = ' < 0.602 ' # filt.cut_misshits_endcap_medium = ' < 2 ' # filt.cut_passConvVeto_endcap_medium = ' == 1 ' # # filt.cut_sigmaIEIE_endcap_loose = ' < 0.0314 ' # filt.cut_absdEtaIn_endcap_loose = ' < 0.00868 ' # filt.cut_absdPhiIn_endcap_loose = ' < 0.213 ' # filt.cut_hovere_endcap_loose = ' < 0.101 ' # filt.cut_isoRho_endcap_loose = ' < 0.107 ' # filt.cut_ooEmooP_endcap_loose = ' < 0.14 ' # #filt.cut_d0_endcap_loose = ' < 0.118 ' # #filt.cut_z0_endcap_loose = ' < 0.822 ' # filt.cut_misshits_endcap_loose = ' < 2 ' # filt.cut_passConvVeto_endcap_loose = ' == 1 ' # # filt.cut_sigmaIEIE_endcap_veryloose = ' < 0.037 ' # filt.cut_absdEtaIn_endcap_veryloose = ' < 0.00895 ' # filt.cut_absdPhiIn_endcap_veryloose = ' < 0.213 ' # filt.cut_hovere_endcap_veryloose = ' < 0.211 ' # filt.cut_isoRho_endcap_veryloose = ' < 0.159 ' # filt.cut_ooEmooP_endcap_veryloose = ' < 0.15 ' # #filt.cut_d0_endcap_veryloose = ' < 0.222 ' # #filt.cut_z0_endcap_veryloose = ' < 0.921 ' # filt.cut_misshits_endcap_veryloose = ' < 4 ' # filt.cut_passConvVeto_endcap_veryloose = ' == 1 ' if do_cutflow: filt.do_cutflow = True filt.add_var('evalPID', evalPID) if apply_corrections: workarea = os.getenv('WorkArea') filt.add_var('applyCorrections', 'true') filt.add_var( 'correctionFile', '%s/TreeFilter/RecoWgg/data/step2-invMass_SC-loose-Et_20-trigger-noPF-HggRunEtaR9.dat' % workarea) filt.add_var( 'smearingFile', '%s/TreeFilter/RecoWgg/data/outFile-step4-invMass_SC-loose-Et_20-trigger-noPF-HggRunEtaR9-smearEle.dat' % workarea) if do_hists: filt.add_hist('cut_pt', 100, 0, 200) filt.add_hist('cut_eta', 50, 0, 5) filt.add_hist('cut_abssceta', 50, 0, 5) filt.add_hist('cut_abseta_crack', 50, 0, 5) filt.add_hist('cut_absdEtaIn_barrel_medium', 100, -0.1, 0.1) filt.add_hist('cut_absdPhiIn_barrel_medium', 100, -0.1, 0.1) filt.add_hist('cut_sigmaIEIE_barrel_medium', 100, 0, 0.05) filt.add_hist('cut_hovere_barrel_medium', 100, -1, 1) filt.add_hist('cut_d0_barrel_medium', 100, -1, 1) filt.add_hist('cut_z0_barrel_medium', 100, -1, 1) filt.add_hist('cut_ooEmooP_barrel_medium', 100, 0, 1) filt.add_hist('cut_pfIso30_barrel_medium', 100, 0, 10) filt.add_hist('cut_passConvVeto_barrel_medium', 2, 0, 2) filt.add_hist('cut_misshits_barrel_medium', 10, 0, 10) return filt
def build_electron( do_cutflow=False, do_hists=False, filtPID=None, evalPID=None, applyCorrections=False ) : filt = Filter('BuildElectron') filt.do_cutflow = do_cutflow filt.cut_pt = ' > 5' filt.cut_abssceta = ' <2.5 ' # no crack for now #filt.cut_abssceta_crack = ' > 1.44 & < 1.57 ' #filt.invert('cut_abssceta_crack') filt.cut_absdEtaIn_barrel_tight = ' < 0.004 ' filt.cut_absdPhiIn_barrel_tight = ' < 0.03 ' filt.cut_sigmaIEIE_barrel_tight = ' < 0.01 ' filt.cut_hovere_barrel_tight = ' < 0.12 ' filt.cut_d0_barrel_tight = ' < 0.02 ' filt.cut_z0_barrel_tight = ' < 0.1 ' filt.cut_eoverp_barrel_tight = ' < 0.05 ' filt.cut_pfIso30_barrel_tight = ' < 0.1 ' filt.cut_convfit_barrel_tight = ' == 0 ' filt.cut_misshits_barrel_tight = ' == 0 ' filt.cut_absdEtaIn_barrel_medium = ' < 0.004 ' filt.cut_absdPhiIn_barrel_medium = ' < 0.06 ' filt.cut_sigmaIEIE_barrel_medium = ' < 0.01 ' filt.cut_hovere_barrel_medium = ' < 0.12 ' filt.cut_d0_barrel_medium = ' < 0.02 ' filt.cut_z0_barrel_medium = ' < 0.1 ' filt.cut_eoverp_barrel_medium = ' < 0.05 ' filt.cut_pfIso30_barrel_medium = ' < 0.15 ' filt.cut_convfit_barrel_medium = ' == 0 ' filt.cut_misshits_barrel_medium = ' <= 1 ' filt.cut_absdEtaIn_barrel_loose = ' < 0.007 ' filt.cut_absdPhiIn_barrel_loose = ' < 0.15 ' filt.cut_sigmaIEIE_barrel_loose = ' < 0.01 ' filt.cut_hovere_barrel_loose = ' < 0.12 ' filt.cut_d0_barrel_loose = ' < 0.02 ' filt.cut_z0_barrel_loose = ' < 0.2 ' filt.cut_eoverp_barrel_loose = ' < 0.05 ' filt.cut_pfIso30_barrel_loose = ' < 0.15 ' filt.cut_convfit_barrel_loose = ' == 0 ' filt.cut_misshits_barrel_loose = ' <= 1 ' filt.cut_absdEtaIn_barrel_veryloose = ' < 0.007 ' filt.cut_absdPhiIn_barrel_veryloose = ' < 0.8 ' filt.cut_sigmaIEIE_barrel_veryloose = ' < 0.01 ' filt.cut_hovere_barrel_veryloose = ' < 0.15 ' filt.cut_d0_barrel_veryloose = ' < 0.04 ' filt.cut_z0_barrel_veryloose = ' < 0.2 ' #filt.cut_eoverp_barrel_veryloose = ' < 0.05 ' #no cut filt.cut_pfIso30_barrel_veryloose = ' < 0.1 ' #filt.cut_convfit_barrel_veryloose = ' == 0 ' #no cut #filt.cut_misshits_barrel_veryloose = ' == 0 ' #no cut filt.cut_absdEtaIn_barrel_tightTrig = ' < 0.007 ' filt.cut_absdPhiIn_barrel_tightTrig = ' < 0.15 ' filt.cut_sigmaIEIE_barrel_tightTrig = ' < 0.01 ' filt.cut_hovere_barrel_tightTrig = ' < 0.12 ' filt.cut_ecalIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_hcalIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_trkIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_absdEtaIn_endcap_tight = ' < 0.005 ' filt.cut_absdPhiIn_endcap_tight = ' < 0.02 ' filt.cut_sigmaIEIE_endcap_tight = ' < 0.03 ' filt.cut_hovere_endcap_tight = ' < 0.1 ' filt.cut_d0_endcap_tight = ' < 0.02 ' filt.cut_z0_endcap_tight = ' < 0.1 ' filt.cut_eoverp_endcap_tight = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_tight = ' < 0.1 ' filt.cut_pfIso30_endcap_highPt_tight = ' < 0.15 ' filt.cut_convfit_endcap_tight = ' == 0 ' filt.cut_misshits_endcap_tight = ' == 0 ' filt.cut_absdEtaIn_endcap_medium = ' < 0.007 ' filt.cut_absdPhiIn_endcap_medium = ' < 0.03 ' filt.cut_sigmaIEIE_endcap_medium = ' < 0.03 ' filt.cut_hovere_endcap_medium = ' < 0.1 ' filt.cut_d0_endcap_medium = ' < 0.02 ' filt.cut_z0_endcap_medium = ' < 0.1 ' filt.cut_eoverp_endcap_medium = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_medium = ' < 0.1 ' filt.cut_pfIso30_endcap_highPt_medium = ' < 0.15 ' filt.cut_convfit_endcap_medium = ' == 0 ' filt.cut_misshits_endcap_medium = ' <= 1 ' filt.cut_absdEtaIn_endcap_loose = ' < 0.009 ' filt.cut_absdPhiIn_endcap_loose = ' < 0.1 ' filt.cut_sigmaIEIE_endcap_loose = ' < 0.03 ' filt.cut_hovere_endcap_loose = ' < 0.1 ' filt.cut_d0_endcap_loose = ' < 0.02 ' filt.cut_z0_endcap_loose = ' < 0.2 ' filt.cut_eoverp_endcap_loose = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_loose = ' < 0.10 ' filt.cut_pfIso30_endcap_highPt_loose = ' < 0.15 ' filt.cut_convfit_endcap_loose = ' == 0 ' filt.cut_misshits_endcap_loose = ' <= 1 ' filt.cut_absdEtaIn_endcap_veryloose = ' < 0.01 ' filt.cut_absdPhiIn_endcap_veryloose = ' < 0.7 ' filt.cut_sigmaIEIE_endcap_veryloose = ' < 0.03 ' #filt.cut_hovere_endcap_veryloose = ' < 0.15 ' #no cut filt.cut_d0_endcap_veryloose = ' < 0.04 ' filt.cut_z0_endcap_veryloose = ' < 0.2 ' #filt.cut_eoverp_endcap_veryloose = ' < 0.05 ' #no cut filt.cut_pfIso30_endcap_veryloose = ' < 0.15 ' #filt.cut_convfit_endcap_veryloose = ' == 0 ' #no cut #filt.cut_misshits_endcap_veryloose = ' == 0 ' #no cut filt.cut_absdEtaIn_endcap_tightTrig = ' < 0.009 ' filt.cut_absdPhiIn_endcap_tightTrig = ' < 0.10 ' filt.cut_sigmaIEIE_endcap_tightTrig = ' < 0.03 ' filt.cut_hovere_endcap_tightTrig = ' < 0.1 ' filt.cut_ecalIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_hcalIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_trkIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_mva_central_lowpt_mvatrig = ' > 0.0 ' filt.cut_mva_crack_lowpt_mvatrig = ' > 0.1 ' filt.cut_mva_endcap_lowpt_mvatrig = ' > 0.62 ' filt.cut_mva_central_highpt_mvatrig = ' > 0.94 ' filt.cut_mva_crack_highpt_mvatrig = ' > 0.85 ' filt.cut_mva_endcap_highpt_mvatrig = ' > 0.92 ' filt.cut_relpfiso_mvatrig = ' < 0.15 ' filt.cut_misshits_mvatrig = ' == 0 ' filt.cut_convfit_mvatrig = ' == 0 ' filt.cut_mva_central_lowpt_mvanontrig = ' > 0.47 ' filt.cut_mva_crack_lowpt_mvanontrig = ' > 0.004 ' filt.cut_mva_endcap_lowpt_mvanontrig = ' > 0.295 ' filt.cut_mva_central_highpt_mvanontrig = ' > -0.34 ' filt.cut_mva_crack_highpt_mvanontrig = ' > -0.65 ' filt.cut_mva_endcap_highpt_mvanontrig = ' > 0.6 ' filt.cut_relpfiso_mvanontrig = ' < 0.4 ' filt.cut_misshits_mvanontrig = ' < 2 ' filt.cut_sip_mvanontrig = ' < 4' if filtPID is not None : setattr(filt, 'cut_pid_%s' %filtPID, ' == True' ) if evalPID is not None : filt.add_var( 'evalPID', evalPID ) if applyCorrections : workarea = os.getenv('WorkArea') filt.add_var('applyCorrections', 'true' ) filt.add_var('correctionFile', '%s/TreeFilter/RecoWgg/data/step2-invMass_SC-loose-Et_20-trigger-noPF-HggRunEtaR9.dat' %workarea ) filt.add_var('smearingFile', '%s/TreeFilter/RecoWgg/data/outFile-step4-invMass_SC-loose-Et_20-trigger-noPF-HggRunEtaR9-smearEle.dat' %workarea ) if do_hists : filt.add_hist( 'cut_pt', 100, 0, 500 ) filt.add_hist( 'cut_abseta', 50, 0, 5 ) filt.add_hist( 'cut_abseta_crack', 50, 0, 5 ) filt.add_hist( 'cut_absdEtaIn_barrel_tight', 100, -0.1, 0.1 ) filt.add_hist( 'cut_absdPhiIn_barrel_tight', 100, -0.1, 0.1 ) filt.add_hist( 'cut_sigmaIEIE_barrel_tight', 100, 0, 0.05 ) filt.add_hist( 'cut_hovere_barrel_tight', 100, -1, 1 ) filt.add_hist( 'cut_d0_barrel_tight', 100, -1, 1 ) filt.add_hist( 'cut_z0_barrel_tight', 100, -1, 1 ) filt.add_hist( 'cut_eoverp_barrel_tight', 100, 0, 1 ) filt.add_hist( 'cut_pfIso30_barrel_tight', 100, 0, 10 ) filt.add_hist( 'cut_convfit_barrel_tight', 2, 0, 2 ) filt.add_hist( 'cut_misshits_barrel_tight', 10, 0, 10 ) filt.add_hist( 'cut_mva_central_lowpt_mvatrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_crack_lowpt_mvatrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_endcap_lowpt_mvatrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_central_highpt_mvatrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_crack_highpt_mvatrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_endcap_highpt_mvatrig', 50, -1, 1 ) filt.add_hist( 'cut_relpfiso_mvatrig', 50, 0, 5 ) filt.add_hist( 'cut_misshits_mvatrig', 10, 0, 10 ) filt.add_hist( 'cut_convfit_mvatrig', 2, 0, 1 ) filt.add_hist( 'cut_mva_central_lowpt_mvanontrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_crack_lowpt_mvanontrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_endcap_lowpt_mvanontrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_central_highpt_mvanontrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_crack_highpt_mvanontrig', 50, -1, 1 ) filt.add_hist( 'cut_mva_endcap_highpt_mvanontrig', 50, -1, 1 ) filt.add_hist( 'cut_relpfiso_mvanontrig', 50, 0, 5 ) filt.add_hist( 'cut_misshits_mvanontrig', 10, 0, 10 ) filt.add_hist( 'cut_sip_mvanontrig', 50, 0, 10 ) return filt
def build_electron( do_cutflow=False, do_hists=False, filtPID=None, evalPID=None ) : filt = Filter('BuildElectron') filt.do_cutflow = do_cutflow filt.cut_pt = ' > 25' filt.cut_abssceta = ' <2.5 ' # no crack for now #filt.cut_abssceta_crack = ' > 1.44 & < 1.57 ' #filt.invert('cut_abssceta_crack') filt.cut_absdEtaIn_barrel_tight = ' < 0.004 ' filt.cut_absdPhiIn_barrel_tight = ' < 0.03 ' filt.cut_sigmaIEIE_barrel_tight = ' < 0.01 ' filt.cut_hovere_barrel_tight = ' < 0.12 ' filt.cut_d0_barrel_tight = ' < 0.02 ' filt.cut_z0_barrel_tight = ' < 0.1 ' filt.cut_eoverp_barrel_tight = ' < 0.05 ' filt.cut_pfIso30_barrel_tight = ' < 0.1 ' filt.cut_convfit_barrel_tight = ' < 0.000001 ' filt.cut_misshits_barrel_tight = ' == 0 ' filt.cut_absdEtaIn_barrel_medium = ' < 0.004 ' filt.cut_absdPhiIn_barrel_medium = ' < 0.06 ' filt.cut_sigmaIEIE_barrel_medium = ' < 0.01 ' filt.cut_hovere_barrel_medium = ' < 0.12 ' filt.cut_d0_barrel_medium = ' < 0.02 ' filt.cut_z0_barrel_medium = ' < 0.1 ' filt.cut_eoverp_barrel_medium = ' < 0.05 ' filt.cut_pfIso30_barrel_medium = ' < 0.15 ' filt.cut_convfit_barrel_medium = ' < 0.000001 ' filt.cut_misshits_barrel_medium = ' <= 1 ' filt.cut_absdEtaIn_barrel_loose = ' < 0.007 ' filt.cut_absdPhiIn_barrel_loose = ' < 0.15 ' filt.cut_sigmaIEIE_barrel_loose = ' < 0.01 ' filt.cut_hovere_barrel_loose = ' < 0.12 ' filt.cut_d0_barrel_loose = ' < 0.02 ' filt.cut_z0_barrel_loose = ' < 0.2 ' filt.cut_eoverp_barrel_loose = ' < 0.05 ' filt.cut_pfIso30_barrel_loose = ' < 0.15 ' filt.cut_convfit_barrel_loose = ' < 0.000001 ' filt.cut_misshits_barrel_loose = ' <= 1 ' filt.cut_absdEtaIn_barrel_veryloose = ' < 0.007 ' filt.cut_absdPhiIn_barrel_veryloose = ' < 0.8 ' filt.cut_sigmaIEIE_barrel_veryloose = ' < 0.01 ' filt.cut_hovere_barrel_veryloose = ' < 0.15 ' filt.cut_d0_barrel_veryloose = ' < 0.04 ' filt.cut_z0_barrel_veryloose = ' < 0.2 ' #filt.cut_eoverp_barrel_veryloose = ' < 0.05 ' #no cut filt.cut_pfIso30_barrel_veryloose = ' < 0.1 ' #filt.cut_convfit_barrel_veryloose = ' < 0.000001 ' #no cut #filt.cut_misshits_barrel_veryloose = ' == 0 ' #no cut filt.cut_absdEtaIn_barrel_tightTrig = ' < 0.007 ' filt.cut_absdPhiIn_barrel_tightTrig = ' < 0.15 ' filt.cut_sigmaIEIE_barrel_tightTrig = ' < 0.01 ' filt.cut_hovere_barrel_tightTrig = ' < 0.12 ' filt.cut_ecalIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_hcalIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_trkIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_absdEtaIn_endcap_tight = ' < 0.005 ' filt.cut_absdPhiIn_endcap_tight = ' < 0.02 ' filt.cut_sigmaIEIE_endcap_tight = ' < 0.03 ' filt.cut_hovere_endcap_tight = ' < 0.1 ' filt.cut_d0_endcap_tight = ' < 0.02 ' filt.cut_z0_endcap_tight = ' < 0.1 ' filt.cut_eoverp_endcap_tight = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_tight = ' < 0.1 ' filt.cut_pfIso30_endcap_highPt_tight = ' < 0.15 ' filt.cut_convfit_endcap_tight = ' < 0.000001 ' filt.cut_misshits_endcap_tight = ' == 0 ' filt.cut_absdEtaIn_endcap_medium = ' < 0.007 ' filt.cut_absdPhiIn_endcap_medium = ' < 0.03 ' filt.cut_sigmaIEIE_endcap_medium = ' < 0.03 ' filt.cut_hovere_endcap_medium = ' < 0.1 ' filt.cut_d0_endcap_medium = ' < 0.02 ' filt.cut_z0_endcap_medium = ' < 0.1 ' filt.cut_eoverp_endcap_medium = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_medium = ' < 0.1 ' filt.cut_pfIso30_endcap_highPt_medium = ' < 0.15 ' filt.cut_convfit_endcap_medium = ' < 0.000001 ' filt.cut_misshits_endcap_medium = ' <= 1 ' filt.cut_absdEtaIn_endcap_loose = ' < 0.009 ' filt.cut_absdPhiIn_endcap_loose = ' < 0.1 ' filt.cut_sigmaIEIE_endcap_loose = ' < 0.03 ' filt.cut_hovere_endcap_loose = ' < 0.1 ' filt.cut_d0_endcap_loose = ' < 0.02 ' filt.cut_z0_endcap_loose = ' < 0.2 ' filt.cut_eoverp_endcap_loose = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_loose = ' < 0.10 ' filt.cut_pfIso30_endcap_highPt_loose = ' < 0.15 ' filt.cut_convfit_endcap_loose = ' < 0.000001 ' filt.cut_misshits_endcap_loose = ' <= 1 ' filt.cut_absdEtaIn_endcap_veryloose = ' < 0.01 ' filt.cut_absdPhiIn_endcap_veryloose = ' < 0.7 ' filt.cut_sigmaIEIE_endcap_veryloose = ' < 0.03 ' #filt.cut_hovere_endcap_veryloose = ' < 0.15 ' #no cut filt.cut_d0_endcap_veryloose = ' < 0.04 ' filt.cut_z0_endcap_veryloose = ' < 0.2 ' #filt.cut_eoverp_endcap_veryloose = ' < 0.05 ' #no cut filt.cut_pfIso30_endcap_veryloose = ' < 0.15 ' #filt.cut_convfit_endcap_veryloose = ' < 0.000001 ' #no cut #filt.cut_misshits_endcap_veryloose = ' == 0 ' #no cut filt.cut_absdEtaIn_endcap_tightTrig = ' < 0.009 ' filt.cut_absdPhiIn_endcap_tightTrig = ' < 0.10 ' filt.cut_sigmaIEIE_endcap_tightTrig = ' < 0.03 ' filt.cut_hovere_endcap_tightTrig = ' < 0.1 ' filt.cut_ecalIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_hcalIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_trkIso30_endcap_tightTrig = ' < 0.2 ' if filtPID is not None : setattr(filt, 'cut_pid_%s' %filtPID, ' == True' ) if evalPID is not None : filt.add_var( 'evalPID', evalPID ) if do_hists : filt.add_hist( 'cut_pt', 100, 0, 500 ) filt.add_hist( 'cut_abseta', 50, 0, 5 ) filt.add_hist( 'cut_abseta_crack', 50, 0, 5 ) filt.add_hist( 'cut_absdEtaIn_barrel_tight', 100, -0.1, 0.1 ) filt.add_hist( 'cut_absdPhiIn_barrel_tight', 100, -0.1, 0.1 ) filt.add_hist( 'cut_sigmaIEIE_barrel_tight', 100, 0, 0.05 ) filt.add_hist( 'cut_hovere_barrel_tight', 100, -1, 1 ) filt.add_hist( 'cut_d0_barrel_tight', 100, -1, 1 ) filt.add_hist( 'cut_z0_barrel_tight', 100, -1, 1 ) filt.add_hist( 'cut_eoverp_barrel_tight', 100, 0, 1 ) filt.add_hist( 'cut_pfIso30_barrel_tight', 100, 0, 10 ) filt.add_hist( 'cut_convfit_barrel_tight', 2, 0, 2 ) filt.add_hist( 'cut_misshits_barrel_tight', 10, 0, 10 ) return filt
def build_electron(do_cutflow=False, do_hists=False, filtPID=None, evalPID=None, applyCorrections=False): filt = Filter('BuildElectron') filt.do_cutflow = do_cutflow filt.cut_pt = ' > 10' filt.cut_abssceta = ' <2.5 ' # no crack for now #filt.cut_abssceta_crack = ' > 1.44 & < 1.57 ' #filt.invert('cut_abssceta_crack') filt.cut_absdEtaIn_barrel_tight = ' < 0.004 ' filt.cut_absdPhiIn_barrel_tight = ' < 0.03 ' filt.cut_sigmaIEIE_barrel_tight = ' < 0.01 ' filt.cut_hovere_barrel_tight = ' < 0.12 ' filt.cut_d0_barrel_tight = ' < 0.02 ' filt.cut_z0_barrel_tight = ' < 0.1 ' filt.cut_eoverp_barrel_tight = ' < 0.05 ' filt.cut_pfIso30_barrel_tight = ' < 0.1 ' filt.cut_convfit_barrel_tight = ' == 0 ' filt.cut_misshits_barrel_tight = ' == 0 ' filt.cut_absdEtaIn_barrel_medium = ' < 0.004 ' filt.cut_absdPhiIn_barrel_medium = ' < 0.06 ' filt.cut_sigmaIEIE_barrel_medium = ' < 0.01 ' filt.cut_hovere_barrel_medium = ' < 0.12 ' filt.cut_d0_barrel_medium = ' < 0.02 ' filt.cut_z0_barrel_medium = ' < 0.1 ' filt.cut_eoverp_barrel_medium = ' < 0.05 ' filt.cut_pfIso30_barrel_medium = ' < 0.15 ' filt.cut_convfit_barrel_medium = ' == 0 ' filt.cut_misshits_barrel_medium = ' <= 1 ' filt.cut_absdEtaIn_barrel_loose = ' < 0.007 ' filt.cut_absdPhiIn_barrel_loose = ' < 0.15 ' filt.cut_sigmaIEIE_barrel_loose = ' < 0.01 ' filt.cut_hovere_barrel_loose = ' < 0.12 ' filt.cut_d0_barrel_loose = ' < 0.02 ' filt.cut_z0_barrel_loose = ' < 0.2 ' filt.cut_eoverp_barrel_loose = ' < 0.05 ' filt.cut_pfIso30_barrel_loose = ' < 0.15 ' filt.cut_convfit_barrel_loose = ' == 0 ' filt.cut_misshits_barrel_loose = ' <= 1 ' filt.cut_absdEtaIn_barrel_veryloose = ' < 0.007 ' filt.cut_absdPhiIn_barrel_veryloose = ' < 0.8 ' filt.cut_sigmaIEIE_barrel_veryloose = ' < 0.01 ' filt.cut_hovere_barrel_veryloose = ' < 0.15 ' filt.cut_d0_barrel_veryloose = ' < 0.04 ' filt.cut_z0_barrel_veryloose = ' < 0.2 ' #filt.cut_eoverp_barrel_veryloose = ' < 0.05 ' #no cut filt.cut_pfIso30_barrel_veryloose = ' < 0.1 ' #filt.cut_convfit_barrel_veryloose = ' == 0 ' #no cut #filt.cut_misshits_barrel_veryloose = ' == 0 ' #no cut filt.cut_absdEtaIn_barrel_tightTrig = ' < 0.007 ' filt.cut_absdPhiIn_barrel_tightTrig = ' < 0.15 ' filt.cut_sigmaIEIE_barrel_tightTrig = ' < 0.01 ' filt.cut_hovere_barrel_tightTrig = ' < 0.12 ' filt.cut_ecalIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_hcalIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_trkIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_absdEtaIn_endcap_tight = ' < 0.005 ' filt.cut_absdPhiIn_endcap_tight = ' < 0.02 ' filt.cut_sigmaIEIE_endcap_tight = ' < 0.03 ' filt.cut_hovere_endcap_tight = ' < 0.1 ' filt.cut_d0_endcap_tight = ' < 0.02 ' filt.cut_z0_endcap_tight = ' < 0.1 ' filt.cut_eoverp_endcap_tight = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_tight = ' < 0.1 ' filt.cut_pfIso30_endcap_highPt_tight = ' < 0.15 ' filt.cut_convfit_endcap_tight = ' == 0 ' filt.cut_misshits_endcap_tight = ' == 0 ' filt.cut_absdEtaIn_endcap_medium = ' < 0.007 ' filt.cut_absdPhiIn_endcap_medium = ' < 0.03 ' filt.cut_sigmaIEIE_endcap_medium = ' < 0.03 ' filt.cut_hovere_endcap_medium = ' < 0.1 ' filt.cut_d0_endcap_medium = ' < 0.02 ' filt.cut_z0_endcap_medium = ' < 0.1 ' filt.cut_eoverp_endcap_medium = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_medium = ' < 0.1 ' filt.cut_pfIso30_endcap_highPt_medium = ' < 0.15 ' filt.cut_convfit_endcap_medium = ' == 0 ' filt.cut_misshits_endcap_medium = ' <= 1 ' filt.cut_absdEtaIn_endcap_loose = ' < 0.009 ' filt.cut_absdPhiIn_endcap_loose = ' < 0.1 ' filt.cut_sigmaIEIE_endcap_loose = ' < 0.03 ' filt.cut_hovere_endcap_loose = ' < 0.1 ' filt.cut_d0_endcap_loose = ' < 0.02 ' filt.cut_z0_endcap_loose = ' < 0.2 ' filt.cut_eoverp_endcap_loose = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_loose = ' < 0.10 ' filt.cut_pfIso30_endcap_highPt_loose = ' < 0.15 ' filt.cut_convfit_endcap_loose = ' == 0 ' filt.cut_misshits_endcap_loose = ' <= 1 ' filt.cut_absdEtaIn_endcap_veryloose = ' < 0.01 ' filt.cut_absdPhiIn_endcap_veryloose = ' < 0.7 ' filt.cut_sigmaIEIE_endcap_veryloose = ' < 0.03 ' #filt.cut_hovere_endcap_veryloose = ' < 0.15 ' #no cut filt.cut_d0_endcap_veryloose = ' < 0.04 ' filt.cut_z0_endcap_veryloose = ' < 0.2 ' #filt.cut_eoverp_endcap_veryloose = ' < 0.05 ' #no cut filt.cut_pfIso30_endcap_veryloose = ' < 0.15 ' #filt.cut_convfit_endcap_veryloose = ' == 0 ' #no cut #filt.cut_misshits_endcap_veryloose = ' == 0 ' #no cut filt.cut_absdEtaIn_endcap_tightTrig = ' < 0.009 ' filt.cut_absdPhiIn_endcap_tightTrig = ' < 0.10 ' filt.cut_sigmaIEIE_endcap_tightTrig = ' < 0.03 ' filt.cut_hovere_endcap_tightTrig = ' < 0.1 ' filt.cut_ecalIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_hcalIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_trkIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_mva_central_lowpt_mvatrig = ' > 0.0 ' filt.cut_mva_crack_lowpt_mvatrig = ' > 0.1 ' filt.cut_mva_endcap_lowpt_mvatrig = ' > 0.62 ' filt.cut_mva_central_highpt_mvatrig = ' > 0.94 ' filt.cut_mva_crack_highpt_mvatrig = ' > 0.85 ' filt.cut_mva_endcap_highpt_mvatrig = ' > 0.92 ' filt.cut_relpfiso_mvatrig = ' < 0.15 ' filt.cut_misshits_mvatrig = ' == 0 ' filt.cut_convfit_mvatrig = ' == 0 ' filt.cut_mva_central_lowpt_mvanontrig = ' > 0.47 ' filt.cut_mva_crack_lowpt_mvanontrig = ' > 0.004 ' filt.cut_mva_endcap_lowpt_mvanontrig = ' > 0.295 ' filt.cut_mva_central_highpt_mvanontrig = ' > -0.34 ' filt.cut_mva_crack_highpt_mvanontrig = ' > -0.65 ' filt.cut_mva_endcap_highpt_mvanontrig = ' > 0.6 ' filt.cut_relpfiso_mvanontrig = ' < 0.4 ' filt.cut_misshits_mvanontrig = ' < 2 ' filt.cut_sip_mvanontrig = ' < 4' if filtPID is not None: setattr(filt, 'cut_pid_%s' % filtPID, ' == True') if evalPID is not None: filt.add_var('evalPID', evalPID) if applyCorrections: workarea = os.getenv('WorkArea') filt.add_var('applyCorrections', 'true') filt.add_var( 'correctionFile', '%s/TreeFilter/RecoWgg/data/scalesNewReg-May2013.csv' % workarea) filt.add_var( 'linCorrectionFile', '%s/TreeFilter/RecoWgg/data/linearityNewReg-May2013.csv' % workarea) if do_hists: filt.add_hist('cut_pt', 100, 0, 500) filt.add_hist('cut_abseta', 50, 0, 5) filt.add_hist('cut_abseta_crack', 50, 0, 5) filt.add_hist('cut_absdEtaIn_barrel_tight', 100, -0.1, 0.1) filt.add_hist('cut_absdPhiIn_barrel_tight', 100, -0.1, 0.1) filt.add_hist('cut_sigmaIEIE_barrel_tight', 100, 0, 0.05) filt.add_hist('cut_hovere_barrel_tight', 100, -1, 1) filt.add_hist('cut_d0_barrel_tight', 100, -1, 1) filt.add_hist('cut_z0_barrel_tight', 100, -1, 1) filt.add_hist('cut_eoverp_barrel_tight', 100, 0, 1) filt.add_hist('cut_pfIso30_barrel_tight', 100, 0, 10) filt.add_hist('cut_convfit_barrel_tight', 2, 0, 2) filt.add_hist('cut_misshits_barrel_tight', 10, 0, 10) filt.add_hist('cut_mva_central_lowpt_mvatrig', 50, -1, 1) filt.add_hist('cut_mva_crack_lowpt_mvatrig', 50, -1, 1) filt.add_hist('cut_mva_endcap_lowpt_mvatrig', 50, -1, 1) filt.add_hist('cut_mva_central_highpt_mvatrig', 50, -1, 1) filt.add_hist('cut_mva_crack_highpt_mvatrig', 50, -1, 1) filt.add_hist('cut_mva_endcap_highpt_mvatrig', 50, -1, 1) filt.add_hist('cut_relpfiso_mvatrig', 50, 0, 5) filt.add_hist('cut_misshits_mvatrig', 10, 0, 10) filt.add_hist('cut_convfit_mvatrig', 2, 0, 1) filt.add_hist('cut_mva_central_lowpt_mvanontrig', 50, -1, 1) filt.add_hist('cut_mva_crack_lowpt_mvanontrig', 50, -1, 1) filt.add_hist('cut_mva_endcap_lowpt_mvanontrig', 50, -1, 1) filt.add_hist('cut_mva_central_highpt_mvanontrig', 50, -1, 1) filt.add_hist('cut_mva_crack_highpt_mvanontrig', 50, -1, 1) filt.add_hist('cut_mva_endcap_highpt_mvanontrig', 50, -1, 1) filt.add_hist('cut_relpfiso_mvanontrig', 50, 0, 5) filt.add_hist('cut_misshits_mvanontrig', 10, 0, 10) filt.add_hist('cut_sip_mvanontrig', 50, 0, 10) return filt
def build_electron(do_cutflow=False, do_hists=False, filtPID=None, evalPID=None): filt = Filter('BuildElectron') filt.do_cutflow = do_cutflow filt.cut_pt = ' > 25' filt.cut_abssceta = ' <2.5 ' # no crack for now #filt.cut_abssceta_crack = ' > 1.44 & < 1.57 ' #filt.invert('cut_abssceta_crack') filt.cut_absdEtaIn_barrel_tight = ' < 0.004 ' filt.cut_absdPhiIn_barrel_tight = ' < 0.03 ' filt.cut_sigmaIEIE_barrel_tight = ' < 0.01 ' filt.cut_hovere_barrel_tight = ' < 0.12 ' filt.cut_d0_barrel_tight = ' < 0.02 ' filt.cut_z0_barrel_tight = ' < 0.1 ' filt.cut_eoverp_barrel_tight = ' < 0.05 ' filt.cut_pfIso30_barrel_tight = ' < 0.1 ' filt.cut_convfit_barrel_tight = ' < 0.000001 ' filt.cut_misshits_barrel_tight = ' == 0 ' filt.cut_absdEtaIn_barrel_medium = ' < 0.004 ' filt.cut_absdPhiIn_barrel_medium = ' < 0.06 ' filt.cut_sigmaIEIE_barrel_medium = ' < 0.01 ' filt.cut_hovere_barrel_medium = ' < 0.12 ' filt.cut_d0_barrel_medium = ' < 0.02 ' filt.cut_z0_barrel_medium = ' < 0.1 ' filt.cut_eoverp_barrel_medium = ' < 0.05 ' filt.cut_pfIso30_barrel_medium = ' < 0.15 ' filt.cut_convfit_barrel_medium = ' < 0.000001 ' filt.cut_misshits_barrel_medium = ' <= 1 ' filt.cut_absdEtaIn_barrel_loose = ' < 0.007 ' filt.cut_absdPhiIn_barrel_loose = ' < 0.15 ' filt.cut_sigmaIEIE_barrel_loose = ' < 0.01 ' filt.cut_hovere_barrel_loose = ' < 0.12 ' filt.cut_d0_barrel_loose = ' < 0.02 ' filt.cut_z0_barrel_loose = ' < 0.2 ' filt.cut_eoverp_barrel_loose = ' < 0.05 ' filt.cut_pfIso30_barrel_loose = ' < 0.15 ' filt.cut_convfit_barrel_loose = ' < 0.000001 ' filt.cut_misshits_barrel_loose = ' <= 1 ' filt.cut_absdEtaIn_barrel_veryloose = ' < 0.007 ' filt.cut_absdPhiIn_barrel_veryloose = ' < 0.8 ' filt.cut_sigmaIEIE_barrel_veryloose = ' < 0.01 ' filt.cut_hovere_barrel_veryloose = ' < 0.15 ' filt.cut_d0_barrel_veryloose = ' < 0.04 ' filt.cut_z0_barrel_veryloose = ' < 0.2 ' #filt.cut_eoverp_barrel_veryloose = ' < 0.05 ' #no cut filt.cut_pfIso30_barrel_veryloose = ' < 0.1 ' #filt.cut_convfit_barrel_veryloose = ' < 0.000001 ' #no cut #filt.cut_misshits_barrel_veryloose = ' == 0 ' #no cut filt.cut_absdEtaIn_barrel_tightTrig = ' < 0.007 ' filt.cut_absdPhiIn_barrel_tightTrig = ' < 0.15 ' filt.cut_sigmaIEIE_barrel_tightTrig = ' < 0.01 ' filt.cut_hovere_barrel_tightTrig = ' < 0.12 ' filt.cut_ecalIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_hcalIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_trkIso30_barrel_tightTrig = ' < 0.2 ' filt.cut_absdEtaIn_endcap_tight = ' < 0.005 ' filt.cut_absdPhiIn_endcap_tight = ' < 0.02 ' filt.cut_sigmaIEIE_endcap_tight = ' < 0.03 ' filt.cut_hovere_endcap_tight = ' < 0.1 ' filt.cut_d0_endcap_tight = ' < 0.02 ' filt.cut_z0_endcap_tight = ' < 0.1 ' filt.cut_eoverp_endcap_tight = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_tight = ' < 0.1 ' filt.cut_pfIso30_endcap_highPt_tight = ' < 0.15 ' filt.cut_convfit_endcap_tight = ' < 0.000001 ' filt.cut_misshits_endcap_tight = ' == 0 ' filt.cut_absdEtaIn_endcap_medium = ' < 0.007 ' filt.cut_absdPhiIn_endcap_medium = ' < 0.03 ' filt.cut_sigmaIEIE_endcap_medium = ' < 0.03 ' filt.cut_hovere_endcap_medium = ' < 0.1 ' filt.cut_d0_endcap_medium = ' < 0.02 ' filt.cut_z0_endcap_medium = ' < 0.1 ' filt.cut_eoverp_endcap_medium = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_medium = ' < 0.1 ' filt.cut_pfIso30_endcap_highPt_medium = ' < 0.15 ' filt.cut_convfit_endcap_medium = ' < 0.000001 ' filt.cut_misshits_endcap_medium = ' <= 1 ' filt.cut_absdEtaIn_endcap_loose = ' < 0.009 ' filt.cut_absdPhiIn_endcap_loose = ' < 0.1 ' filt.cut_sigmaIEIE_endcap_loose = ' < 0.03 ' filt.cut_hovere_endcap_loose = ' < 0.1 ' filt.cut_d0_endcap_loose = ' < 0.02 ' filt.cut_z0_endcap_loose = ' < 0.2 ' filt.cut_eoverp_endcap_loose = ' < 0.05 ' filt.cut_pfIso30_endcap_lowPt_loose = ' < 0.10 ' filt.cut_pfIso30_endcap_highPt_loose = ' < 0.15 ' filt.cut_convfit_endcap_loose = ' < 0.000001 ' filt.cut_misshits_endcap_loose = ' <= 1 ' filt.cut_absdEtaIn_endcap_veryloose = ' < 0.01 ' filt.cut_absdPhiIn_endcap_veryloose = ' < 0.7 ' filt.cut_sigmaIEIE_endcap_veryloose = ' < 0.03 ' #filt.cut_hovere_endcap_veryloose = ' < 0.15 ' #no cut filt.cut_d0_endcap_veryloose = ' < 0.04 ' filt.cut_z0_endcap_veryloose = ' < 0.2 ' #filt.cut_eoverp_endcap_veryloose = ' < 0.05 ' #no cut filt.cut_pfIso30_endcap_veryloose = ' < 0.15 ' #filt.cut_convfit_endcap_veryloose = ' < 0.000001 ' #no cut #filt.cut_misshits_endcap_veryloose = ' == 0 ' #no cut filt.cut_absdEtaIn_endcap_tightTrig = ' < 0.009 ' filt.cut_absdPhiIn_endcap_tightTrig = ' < 0.10 ' filt.cut_sigmaIEIE_endcap_tightTrig = ' < 0.03 ' filt.cut_hovere_endcap_tightTrig = ' < 0.1 ' filt.cut_ecalIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_hcalIso30_endcap_tightTrig = ' < 0.2 ' filt.cut_trkIso30_endcap_tightTrig = ' < 0.2 ' if filtPID is not None: setattr(filt, 'cut_pid_%s' % filtPID, ' == True') if evalPID is not None: filt.add_var('evalPID', evalPID) if do_hists: filt.add_hist('cut_pt', 100, 0, 500) filt.add_hist('cut_abseta', 50, 0, 5) filt.add_hist('cut_abseta_crack', 50, 0, 5) filt.add_hist('cut_absdEtaIn_barrel_tight', 100, -0.1, 0.1) filt.add_hist('cut_absdPhiIn_barrel_tight', 100, -0.1, 0.1) filt.add_hist('cut_sigmaIEIE_barrel_tight', 100, 0, 0.05) filt.add_hist('cut_hovere_barrel_tight', 100, -1, 1) filt.add_hist('cut_d0_barrel_tight', 100, -1, 1) filt.add_hist('cut_z0_barrel_tight', 100, -1, 1) filt.add_hist('cut_eoverp_barrel_tight', 100, 0, 1) filt.add_hist('cut_pfIso30_barrel_tight', 100, 0, 10) filt.add_hist('cut_convfit_barrel_tight', 2, 0, 2) filt.add_hist('cut_misshits_barrel_tight', 10, 0, 10) return filt
def build_electron( do_cutflow=False, do_hists=False, filtPID=None, evalPID=None, applyCorrections=False ) : filt = Filter('BuildElectron') filt.do_cutflow = do_cutflow # using values here #https://twiki.cern.ch/twiki/bin/view/CMS/CutBasedElectronIdentificationRun2 # d0 and z0 cuts not recommeded # as of Dec 2016 filt.cut_pt = ' > 10' filt.cut_abssceta = ' <2.5 ' #filt.cut_abssceta = ' <1.479 ' #filt.cut_abssceta = ' >= 1.479 & < 2.5' # no crack for now #filt.cut_abssceta_crack = ' > 1.4442 & < 1.566 ' #filt.invert('cut_abssceta_crack') #filt.invert('cut_abssceta') filt.cut_sigmaIEIE_barrel_tight = ' < 0.00998 ' filt.cut_absdEtaIn_barrel_tight = ' < 0.00308 ' filt.cut_absdPhiIn_barrel_tight = ' < 0.0816 ' filt.cut_hovere_barrel_tight = ' < 0.0414 ' filt.cut_isoRho_barrel_tight = ' < 0.0588 ' filt.cut_ooEmooP_barrel_tight = ' < 0.0129 ' #filt.cut_d0_barrel_tight = ' < 0.0111 ' #filt.cut_z0_barrel_tight = ' < 0.0466 ' filt.cut_misshits_barrel_tight = ' < 1 ' filt.cut_passConvVeto_barrel_tight = ' == 1 ' filt.cut_sigmaIEIE_barrel_medium = ' < 0.00998 ' filt.cut_absdEtaIn_barrel_medium = ' < 0.00311 ' filt.cut_absdPhiIn_barrel_medium = ' < 0.103 ' filt.cut_hovere_barrel_medium = ' < 0.253 ' filt.cut_isoRho_barrel_medium = ' < 0.0695 ' filt.cut_ooEmooP_barrel_medium = ' < 0.134 ' #filt.cut_d0_barrel_medium = ' < 0.0118 ' #filt.cut_z0_barrel_medium = ' < 0.373 ' filt.cut_misshits_barrel_medium = ' < 1 ' filt.cut_passConvVeto_barrel_medium = ' == 1 ' filt.cut_sigmaIEIE_barrel_loose = ' < 0.011 ' filt.cut_absdEtaIn_barrel_loose = ' < 0.00477 ' filt.cut_absdPhiIn_barrel_loose = ' < 0.222 ' filt.cut_hovere_barrel_loose = ' < 0.298 ' filt.cut_isoRho_barrel_loose = ' < 0.0994 ' filt.cut_ooEmooP_barrel_loose = ' < 0.241 ' #filt.cut_d0_barrel_loose = ' < 0.0261 ' #filt.cut_z0_barrel_loose = ' < 0.41 ' filt.cut_misshits_barrel_loose = ' < 1 ' filt.cut_passConvVeto_barrel_loose = ' == 1 ' filt.cut_sigmaIEIE_barrel_veryloose = ' < 0.0115 ' filt.cut_absdEtaIn_barrel_veryloose = ' < 0.00749 ' filt.cut_absdPhiIn_barrel_veryloose = ' < 0.228 ' filt.cut_hovere_barrel_veryloose = ' < 0.356 ' filt.cut_isoRho_barrel_veryloose = ' < 0.175 ' filt.cut_ooEmooP_barrel_veryloose = ' < 0.299 ' #filt.cut_d0_barrel_veryloose = ' < 0.0564 ' #filt.cut_z0_barrel_veryloose = ' < 0.472 ' filt.cut_misshits_barrel_veryloose = ' < 2 ' filt.cut_passConvVeto_barrel_veryloose = ' == 1 ' filt.cut_sigmaIEIE_endcap_tight = ' < 0.0292 ' filt.cut_absdEtaIn_endcap_tight = ' < 0.00605 ' filt.cut_absdPhiIn_endcap_tight = ' < 0.0394 ' filt.cut_hovere_endcap_tight = ' < 0.0641 ' filt.cut_isoRho_endcap_tight = ' < 0.0571 ' filt.cut_ooEmooP_endcap_tight = ' < 0.0129 ' #filt.cut_d0_endcap_tight = ' < 0.0351 ' #filt.cut_z0_endcap_tight = ' < 0.417 ' filt.cut_misshits_endcap_tight = ' < 1 ' filt.cut_passConvVeto_endcap_tight = ' == 1 ' filt.cut_sigmaIEIE_endcap_medium = ' < 0.0298 ' filt.cut_absdEtaIn_endcap_medium = ' < 0.00609 ' filt.cut_absdPhiIn_endcap_medium = ' < 0.045 ' filt.cut_hovere_endcap_medium = ' < 0.0878 ' filt.cut_isoRho_endcap_medium = ' < 0.0821 ' filt.cut_ooEmooP_endcap_medium = ' < 0.13 ' #filt.cut_d0_endcap_medium = ' < 0.0739 ' #filt.cut_z0_endcap_medium = ' < 0.602 ' filt.cut_misshits_endcap_medium = ' < 1 ' filt.cut_passConvVeto_endcap_medium = ' == 1 ' filt.cut_sigmaIEIE_endcap_loose = ' < 0.0314 ' filt.cut_absdEtaIn_endcap_loose = ' < 0.00868 ' filt.cut_absdPhiIn_endcap_loose = ' < 0.213 ' filt.cut_hovere_endcap_loose = ' < 0.101 ' filt.cut_isoRho_endcap_loose = ' < 0.107 ' filt.cut_ooEmooP_endcap_loose = ' < 0.14 ' #filt.cut_d0_endcap_loose = ' < 0.118 ' #filt.cut_z0_endcap_loose = ' < 0.822 ' filt.cut_misshits_endcap_loose = ' < 1 ' filt.cut_passConvVeto_endcap_loose = ' == 1 ' filt.cut_sigmaIEIE_endcap_veryloose = ' < 0.037 ' filt.cut_absdEtaIn_endcap_veryloose = ' < 0.00895 ' filt.cut_absdPhiIn_endcap_veryloose = ' < 0.213 ' filt.cut_hovere_endcap_veryloose = ' < 0.211 ' filt.cut_isoRho_endcap_veryloose = ' < 0.159 ' filt.cut_ooEmooP_endcap_veryloose = ' < 0.15 ' #filt.cut_d0_endcap_veryloose = ' < 0.222 ' #filt.cut_z0_endcap_veryloose = ' < 0.921 ' filt.cut_misshits_endcap_veryloose = ' < 3 ' filt.cut_passConvVeto_endcap_veryloose = ' == 1 ' if filtPID is not None : setattr(filt, 'cut_pid_%s' %filtPID, ' == True' ) if evalPID is not None : filt.add_var( 'evalPID', evalPID ) if applyCorrections : workarea = os.getenv('WorkArea') filt.add_var('applyCorrections', 'true' ) filt.add_var('correctionFile', '%s/TreeFilter/RecoWgg/data/step2-invMass_SC-loose-Et_20-trigger-noPF-HggRunEtaR9.dat' %workarea ) filt.add_var('smearingFile', '%s/TreeFilter/RecoWgg/data/outFile-step4-invMass_SC-loose-Et_20-trigger-noPF-HggRunEtaR9-smearEle.dat' %workarea ) if do_hists : filt.add_hist( 'cut_pt', 100, 0, 500 ) filt.add_hist( 'cut_abseta', 50, 0, 5 ) filt.add_hist( 'cut_abseta_crack', 50, 0, 5 ) filt.add_hist( 'cut_absdEtaIn_barrel_tight', 100, -0.1, 0.1 ) filt.add_hist( 'cut_absdPhiIn_barrel_tight', 100, -0.1, 0.1 ) filt.add_hist( 'cut_sigmaIEIE_barrel_tight', 100, 0, 0.05 ) filt.add_hist( 'cut_hovere_barrel_tight', 100, -1, 1 ) filt.add_hist( 'cut_d0_barrel_tight', 100, -1, 1 ) filt.add_hist( 'cut_z0_barrel_tight', 100, -1, 1 ) filt.add_hist( 'cut_ooEmooP_barrel_tight', 100, 0, 1 ) filt.add_hist( 'cut_pfIso30_barrel_tight', 100, 0, 10 ) filt.add_hist( 'cut_passConvVeto_barrel_tight', 2, 0, 2 ) filt.add_hist( 'cut_misshits_barrel_tight', 10, 0, 10 ) return filt