コード例 #1
0
 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))
コード例 #2
0
ファイル: WHAnalyzeEMT.py プロジェクト: dabelknap/UWHiggs
 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)
コード例 #3
0
 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)
コード例 #4
0
ファイル: WHAnalyzeEET.py プロジェクト: dabelknap/UWHiggs
 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))
コード例 #5
0
ファイル: WHAnalyzeMMT.py プロジェクト: dabelknap/UWHiggs
 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))
コード例 #6
0
 def obj3_weight(self, row, notUsed1=None, notUsed2=None):
     return frfits.tau_fr(row.tPt)
コード例 #7
0
ファイル: WHAnalyzeEMT.py プロジェクト: mverzett/UWHiggs
 def obj3_weight(self, row, notUsed1=None, notUsed2=None):
     return frfits.tau_fr(row.tPt, row.tAbsEta)