import copy
import PhysicsTools.Heppy.loadlibs
from ROOT import heppy
from PhysicsTools.Heppy.utils.cmsswRelease import isNewerThan
is2012 = isNewerThan('CMSSW_5_2_0')
class RochesterCorrections(object):
def __init__(self):
self.cor = heppy.RochCor()
self.cor2012 = heppy.RochCor2012()
def corrected_p4( self, particle, run ):
'''Returns the corrected p4 for a particle.
The particle remains unchanged.
'''
ptc = particle
p4 = ptc.p4()
tlp4 = TLorentzVector( p4.px(), p4.py(), p4.pz(), p4.energy() )
cortlp4 = copy.copy(tlp4)
if run<100:
if is2012:
self.cor2012.momcor_mc( cortlp4, ptc.charge(), 0.0, 0 )
else:
self.cor.momcor_mc( cortlp4, ptc.charge(), 0.0, 0 )
else: # data
if is2012:
self.cor2012.momcor_data( cortlp4, ptc.charge(), 0.0, 0 )
else:
def __init__(self, cfg_ana, cfg_comp, looperName):
super(JetAnalyzer, self).__init__(cfg_ana, cfg_comp, looperName)
self.btagSF = BTagSF(cfg_ana.btagSFseed)
self.is2012 = isNewerThan('CMSSW_5_2_0')
def __init__(self, cfg_ana, cfg_comp, looperName):
super(JetAnalyzer, self).__init__(cfg_ana, cfg_comp, looperName)
self.btagSF = BTagSF(0)
self.is2012 = isNewerThan('CMSSW_5_2_0')
import copy
from CMGTools.H2TauTau.proto.plotter.categories_common import *
from CMGTools.H2TauTau.proto.plotter.cut import *
from PhysicsTools.Heppy.utils.cmsswRelease import cmsswIs44X, cmsswIs52X
from PhysicsTools.Heppy.utils.cmsswRelease import isNewerThan
pt1 = 30
pt2 = 20 # 2011
if isNewerThan('CMSSW_5_2_0'):
pt2 = 24 #2012
# ELECTRON = lepton 1
# TAU = lepton 2
# this has to be in sync with:
# - H2TauTau/python/proto/analyzers/TauEleAnalyzer.py
# - H2TauTau/Colin/tauEle_2012_cfg.py
# - H2TauTau/python/objects/eleCuts_cff.py
# - H2TauTau/python/objects/tauCuts_cff.py
# - H2TauTau/python/objects/tauEleCuts_cff.py
# inc_sig_tau = Cut('l1_looseMvaIso>0.5 && l1_againstElectronMVA > 0.5 && l1_againstElectronTightMVA2 > 0.5 && l1_againstElectronMedium > 0.5 && l1_againstMuonLoose > 0.5 && l1_dxy<0.045 && l1_dz<0.2 && l1_pt>{pt1}'.format(pt1=pt1))
inc_sig_tau = Cut(
'leptonAccept && thirdLeptonVeto && l1_threeHitIso<1.5 && l1_againstElectronMVA3Medium > 0.5 && l1_againstMuonLoose > 0.5 && l1_dxy<0.045 && l1_dz<0.2 && l1_pt>{pt1}'
.format(pt1=pt1))
inc_sig_ele = Cut(
'l2_relIso05<0.1 && l2_tightId>0.5 && l2_dxy<0.045 && l2_dz<0.2 && l2_pt>{pt2}'
import copy
from ROOT import heppy
from PhysicsTools.Heppy.utils.cmsswRelease import isNewerThan
is2012 = isNewerThan('CMSSW_5_2_0')
class RochesterCorrections(object):
def __init__(self):
self.cor = heppy.RochCor()
self.cor2012 = heppy.RochCor2012()
def corrected_p4(self, particle, run):
'''Returns the corrected p4 for a particle.
The particle remains unchanged.
'''
ptc = particle
p4 = ptc.p4()
tlp4 = TLorentzVector(p4.px(), p4.py(), p4.pz(), p4.energy())
cortlp4 = copy.copy(tlp4)
if run < 100:
if is2012:
self.cor2012.momcor_mc(cortlp4, ptc.charge(), 0.0, 0)
else:
self.cor.momcor_mc(cortlp4, ptc.charge(), 0.0, 0)
else: # data
if is2012:
self.cor2012.momcor_data(cortlp4, ptc.charge(), 0.0, 0)
else:
self.cor.momcor_data(cortlp4, ptc.charge(), 0.0,
import copy
from CMGTools.H2TauTau.proto.plotter.categories_common import *
from CMGTools.H2TauTau.proto.plotter.cut import *
from PhysicsTools.Heppy.utils.cmsswRelease import cmsswIs44X,cmsswIs52X
from PhysicsTools.Heppy.utils.cmsswRelease import isNewerThan
pt1 = 30
pt2 = 20 # 2011
if isNewerThan('CMSSW_5_2_0'):
pt2 = 24 #2012
# ELECTRON = lepton 1
# TAU = lepton 2
# this has to be in sync with:
# - H2TauTau/python/proto/analyzers/TauEleAnalyzer.py
# - H2TauTau/Colin/tauEle_2012_cfg.py
# - H2TauTau/python/objects/eleCuts_cff.py
# - H2TauTau/python/objects/tauCuts_cff.py
# - H2TauTau/python/objects/tauEleCuts_cff.py
# inc_sig_tau = Cut('l1_looseMvaIso>0.5 && l1_againstElectronMVA > 0.5 && l1_againstElectronTightMVA2 > 0.5 && l1_againstElectronMedium > 0.5 && l1_againstMuonLoose > 0.5 && l1_dxy<0.045 && l1_dz<0.2 && l1_pt>{pt1}'.format(pt1=pt1))
inc_sig_tau = Cut('leptonAccept && thirdLeptonVeto && l1_threeHitIso<1.5 && l1_againstElectronMVA3Medium > 0.5 && l1_againstMuonLoose > 0.5 && l1_dxy<0.045 && l1_dz<0.2 && l1_pt>{pt1}'.format(pt1=pt1))
inc_sig_ele = Cut('l2_relIso05<0.1 && l2_tightId>0.5 && l2_dxy<0.045 && l2_dz<0.2 && l2_pt>{pt2}'.format(pt2=pt2))
passleptonvetoes = Cut('leptonAccept > 0.5 && thirdLeptonVeto > 0.5')
#inc_sig = inc_sig_ele & inc_sig_tau
process.recoilCorMETDiTau.enable = True
process.recoilCorMETDiTau.fileCorrectTo = correctFileName
process.recoilCorMETDiTau.leptonLeg = leptonLeg
process.recoilCorMETDiTau.fileZmmData = fileZmmData
process.recoilCorMETDiTau.fileZmmMC = fileZmmMC
process.recoilCorMETDiTau.correctionType = correctionType
else:
print '\tDISABLED'
if runOnMC:
process.metRecoilCorrectionInputSequence.remove(process.genWorZ)
if hasattr(process, 'recoilCorMETTauMu'):
process.recoilCorMETTauMu.enable = False
if hasattr(process, 'recoilCorMETTauEle'):
process.recoilCorMETTauEle.enable = False
if hasattr(process, 'recoilCorMETDiTau'):
process.recoilCorMETDiTau.enable = False
if __name__ == '__main__':
import sys
from PhysicsTools.Heppy.utils.cmsswRelease import isNewerThan
data, mc, type = basicParameters(isNewerThan('CMSSW_5_2_X'))
for line in sys.stdin:
print
line = line.rstrip()
print line
print fileAndLeg(line, True, None, 'di-tau')
process.recoilCorMETDiTau.enable = True
process.recoilCorMETDiTau.fileCorrectTo = correctFileName
process.recoilCorMETDiTau.leptonLeg = leptonLeg
process.recoilCorMETDiTau.fileZmmData = fileZmmData
process.recoilCorMETDiTau.fileZmmMC = fileZmmMC
process.recoilCorMETDiTau.correctionType = correctionType
else:
print '\tDISABLED'
if runOnMC:
process.metRecoilCorrectionInputSequence.remove( process.genWorZ )
if hasattr( process, 'recoilCorMETTauMu'):
process.recoilCorMETTauMu.enable = False
if hasattr( process, 'recoilCorMETTauEle'):
process.recoilCorMETTauEle.enable = False
if hasattr( process, 'recoilCorMETDiTau'):
process.recoilCorMETDiTau.enable = False
if __name__ == '__main__':
import sys
from PhysicsTools.Heppy.utils.cmsswRelease import isNewerThan
data, mc, type = basicParameters(isNewerThan('CMSSW_5_2_X'))
for line in sys.stdin:
print
line = line.rstrip()
print line
print fileAndLeg(line, True, None, 'di-tau')
def __init__(self, cfg_ana, cfg_comp, looperName):
super(VBFAnalyzer, self).__init__(cfg_ana, cfg_comp, looperName)
self.btagSF = BTagSF(cfg_ana.btagSFseed)
# import pdb; pdb.set_trace()
self.is2012 = isNewerThan('CMSSW_5_2_0')