def f_par_prob(e1Pt, e1JetPt, e2Pt, e2JetPt, tPt): p_e1 = (( frfits.highpt_ee_fr[lead_iso](electronJetPt=max(e1JetPt, e1Pt), electronPt=e1Pt) +\ frfits.highpt_ee_qcd_fr[lead_iso](electronJetPt=max(e1JetPt, e1Pt), electronPt=e1Pt) )/2) p_e2 = (( frfits.lowpt_ee_fr[sublead_iso](electronJetPt=max(e2JetPt, e2Pt), electronPt=e2Pt) + \ frfits.lowpt_ee_qcd_fr[sublead_iso](electronJetPt=max(e2JetPt, e2Pt), electronPt=e2Pt))/2) p_t = frfits.tau_fr(tPt) return (p_e1 + p_e2 * (1 - p_e1) + p_t * (1 - p_e1) * (1 - p_e2))
def f_prob(row, weight): p_m = ( self.obj1_weight(row, lead_iso, sublead_iso) + \ self.obj1_qcd_weight(row, lead_iso, sublead_iso))/2 p_e = ( self.obj2_weight(row, lead_iso, sublead_iso) + \ self.obj2_qcd_weight(row, lead_iso, sublead_iso))/2 p_t = frfits.tau_fr(row.tPt) return ((p_m + p_e*(1 - p_m) + p_t*(1 - p_m)*(1 - p_e)), weight)
def f_prob(row, weight): p_m = ( self.obj1_weight(row, lead_iso, sublead_iso) + \ self.obj1_qcd_weight(row, lead_iso, sublead_iso))/2 p_e = ( self.obj2_weight(row, lead_iso, sublead_iso) + \ self.obj2_qcd_weight(row, lead_iso, sublead_iso))/2 p_t = frfits.tau_fr(row.tPt) return ((p_m + p_e * (1 - p_m) + p_t * (1 - p_m) * (1 - p_e)), weight)
def f_par_prob(e1Pt, e1JetPt, e2Pt, e2JetPt, tPt): p_e1 = (( frfits.highpt_ee_fr[lead_iso](electronJetPt=max(e1JetPt, e1Pt), electronPt=e1Pt) +\ frfits.highpt_ee_qcd_fr[lead_iso](electronJetPt=max(e1JetPt, e1Pt), electronPt=e1Pt) )/2) p_e2 = (( frfits.lowpt_ee_fr[sublead_iso](electronJetPt=max(e2JetPt, e2Pt), electronPt=e2Pt) + \ frfits.lowpt_ee_qcd_fr[sublead_iso](electronJetPt=max(e2JetPt, e2Pt), electronPt=e2Pt))/2) p_t = frfits.tau_fr(tPt) return (p_e1 + p_e2*(1 - p_e1) + p_t*(1 - p_e1)*(1 - p_e2))
def f_par_prob(m1Pt, m1JetPt, m2Pt, m2JetPt, tPt): p_m1 = (( frfits.highpt_mu_fr[lead_iso](muonJetPt=max(m1JetPt, m1Pt), muonPt=m1Pt) +\ frfits.highpt_mu_qcd_fr[lead_iso](muonJetPt=max(m1JetPt, m1Pt), muonPt=m1Pt) )/2) p_m2 = (( frfits.lowpt_mu_fr[sublead_iso](muonJetPt=max(m2JetPt, m2Pt), muonPt=m2Pt) + \ frfits.lowpt_mu_qcd_fr[sublead_iso](muonJetPt=max(m2JetPt, m2Pt), muonPt=m2Pt))/2) p_t = frfits.tau_fr(tPt) return (p_m1 + p_m2*(1 - p_m1) + p_t*(1 - p_m1)*(1 - p_m2))
def obj3_weight(self, row, notUsed1=None, notUsed2=None): return frfits.tau_fr(row.tPt)
def obj3_weight(self, row, notUsed1=None, notUsed2=None): return frfits.tau_fr(row.tPt, row.tAbsEta)