def fill(the_histos, row, fillNtuple=False): # Get PU weight - fix me weight = 1 if row.run == 1: weight = self.pucorrector(row.nTruePU) * \ mcCorrectors.double_muon_trigger(row,'m1','m2') the_histos['muonPt'].Fill(row.m2Pt, weight) the_histos['muonJetPt'].Fill(max(row.m2JetPt, row.m2Pt), weight) the_histos['muonAbsEta'].Fill(row.m2AbsEta, weight) the_histos['muonPtRatio'].Fill(row.m2Pt/max(row.m2JetPt, row.m2Pt), weight) the_histos['muonPtDiff'].Fill(max(row.m2JetPt, row.m2Pt) - row.m2Pt, weight) the_histos['m1MtToMET'].Fill(row.m1MtToMET, weight) the_histos['m1m2Mass' ].Fill(row.m1_m2_Mass, weight) the_histos['m2MtToMET'].Fill(row.m2MtToMET, weight) the_histos["m1JetBtag"].Fill(row.m1JetBtag, weight) the_histos["m2JetBtag"].Fill(row.m2JetBtag, weight) the_histos['m2JetptD'].Fill(row.m2JetptD, weight) the_histos['m2Jetmult'].Fill(row.m2Jetmult, weight) if fillNtuple: pfidiso02 = float( row.m2PFIDTight and row.m2RelPFIsoDB < 0.2) h2taucuts = float( row.m2PFIDTight and ((row.m2RelPFIsoDB < 0.15 and row.m2AbsEta < 1.479) or row.m2RelPFIsoDB < 0.1 )) h2taucuts020 = float( row.m2PFIDTight and ((row.m2RelPFIsoDB < 0.20 and row.m2AbsEta < 1.479) or row.m2RelPFIsoDB < 0.15)) muon_jet_mass = -1. #inv_mass(row.m1Pt, row.m1Eta, row.m1Phi, row.leadingJetPt, row.leadingJetEta, row.leadingJetPhi) the_histos['muonInfo'].Fill( array("f", [row.m2Pt, max(row.m2JetPt, row.m2Pt), row.m2AbsEta, max(0, row.m2JetCSVBtag), abs(row.m2PVDXY), muon_jet_mass, row.jetVeto20, row.jetVeto40_DR05, weight, pfidiso02, h2taucuts, h2taucuts020] ) )
def event_weight(self, row): if row.run > 2: return 1. return self.pucorrector(row.nTruePU) * \ mcCorrectors.get_muon_corrections(row,'m1','m2') * \ mcCorrectors.get_electron_corrections(row, 'e') * \ mcCorrectors.double_muon_trigger(row,'m1','m2')
def fill(the_histos, row, fillNtuple=False): # Get PU weight - fix me weight = 1 if row.run == 1: weight = self.pucorrector(row.nTruePU) * \ mcCorrectors.double_muon_trigger(row,'m1','m2') the_histos['muonPt'].Fill(row.m2Pt, weight) the_histos['muonJetPt'].Fill(max(row.m2JetPt, row.m2Pt), weight) the_histos['muonAbsEta'].Fill(row.m2AbsEta, weight) the_histos['muonPtRatio'].Fill( row.m2Pt / max(row.m2JetPt, row.m2Pt), weight) the_histos['muonPtDiff'].Fill( max(row.m2JetPt, row.m2Pt) - row.m2Pt, weight) the_histos['m1MtToMET'].Fill(row.m1MtToMET, weight) the_histos['m1m2Mass'].Fill(row.m1_m2_Mass, weight) the_histos['m2MtToMET'].Fill(row.m2MtToMET, weight) the_histos["m1JetBtag"].Fill(row.m1JetBtag, weight) the_histos["m2JetBtag"].Fill(row.m2JetBtag, weight) the_histos['m2JetptD'].Fill(row.m2JetptD, weight) the_histos['m2Jetmult'].Fill(row.m2Jetmult, weight) if fillNtuple: pfidiso02 = float(row.m2PFIDTight and row.m2RelPFIsoDB < 0.2) h2taucuts = float( row.m2PFIDTight and ((row.m2RelPFIsoDB < 0.15 and row.m2AbsEta < 1.479) or row.m2RelPFIsoDB < 0.1)) h2taucuts020 = float( row.m2PFIDTight and ((row.m2RelPFIsoDB < 0.20 and row.m2AbsEta < 1.479) or row.m2RelPFIsoDB < 0.15)) the_histos['muonInfo'].Fill( array("f", [ row.m2Pt, row.m2JetPt, max(0, row.m2JetCSVBtag), abs(row.m2PVDXY), weight, pfidiso02, h2taucuts, h2taucuts020 ]))
def correction(self, row): return self.pucorrector(row.nTruePU) * \ mcCorrectors.get_muon_corrections(row,'m1','m2') * \ mcCorrectors.double_muon_trigger(row,'m1','m2')
def event_weight(self, row): if row.run > 2: return 1. return self.pucorrector(row.nTruePU) * mcCorrectors.double_muon_trigger(row,'m1','m2')\