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') self.puJetIDDisc = cfg_ana.puJetIDDisc if hasattr(cfg_ana, 'puJetIDDisc') else 'pileupJetId:fullDiscriminant'
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')
import copy from ROOT import RochCor, RochCor2012, TLorentzVector from CMGTools.RootTools.utils.cmsswRelease import isNewerThan is2012 = isNewerThan('CMSSW_5_2_X') class RochesterCorrections(object): def __init__(self): self.cor = RochCor() self.cor2012 = 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') self.puJetIDDisc = cfg_ana.puJetIDDisc if hasattr( cfg_ana, 'puJetIDDisc') else 'pileupJetId:fullDiscriminant'
import copy from ROOT import RochCor, RochCor2012, TLorentzVector from CMGTools.RootTools.utils.cmsswRelease import isNewerThan is2012 = isNewerThan('CMSSW_5_2_X') class RochesterCorrections(object): def __init__(self): self.cor = RochCor() self.cor2012 = 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: