def calcMETSig(event, sample): if sample.isData: JER = JERData params = paramsData else: JER = JERMC params = paramsMC ev = Event(event, JER, isData=sample.isData) ev.calcMETSig(params) event.MET_significance_rec = ev.MET_sig event.Jet_dpt = [x for x in ev.Jet_dpt]
def calcMETSig( event, sample ): if sample.isData: JER = JERData params = paramsData else: JER = JERMC params = paramsMC ev = Event(event, JER, isData=sample.isData, METPtVar=METPtVar, METPhiVar=METPhiVar, JetCollection=JetCollection, vetoEtaRegion=vetoEtaRegion, jetThreshold=jetThresh, pTdepMetSig=pTdependent) ev.calcMETSig(params) event.MET_significance_rec = ev.MET_sig event.Jet_dpt = [ x*ev.Jet_pt[i] for i,x in enumerate(ev.Jet_dpt) ] if len(ev.Jet_dpt) > 0 else [0]
def calcMETSig( event, sample ): if sample.isData: JER = JERData params = paramsData else: JER = JERMC params = paramsMC if year == 2017: METPtVar = "METFixEE2017_pt" METPhiVar = "METFixEE2017_phi" elif year == 2018: METPtVar = "MET_pt_nom" METPhiVar = "MET_phi" else: METPtVar = "MET_pt" METPhiVar = "MET_phi" if year == 2018: JetCollection = "Jet_pt_nom" else: JetCollection = "Jet_pt" ev = Event(event, JER, isData=sample.isData, METPtVar=METPtVar, METPhiVar=METPhiVar, JetCollection=JetCollection, vetoEtaRegion=vetoEtaRegion, jetThreshold=jetThreshold) ev.calcMETSig(params) event.MET_significance_rec = ev.MET_sig event.Jet_dpt = [ x for x in ev.Jet_dpt ]