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