# Make PU corrector from expected data PU distribution # PU corrections .root files from pileupCalc.py pu_distributions = glob.glob( os.path.join( # 'inputs', os.environ['jobid'], 'data_TauPlusX*pu.root')) 'inputs', os.environ['jobid'], 'data_SingleMu*pu.root')) pu_corrector = PileupWeight.PileupWeight('S6' if is7TeV else 'S10', *pu_distributions) muon_pog_PFTight_2011 = MuonPOGCorrections.make_muon_pog_PFTight_2011() muon_pog_PFTight_2012 = MuonPOGCorrections.make_muon_pog_PFTight_2012() muon_pog_PFRelIsoDB02_2011 = MuonPOGCorrections.make_muon_pog_PFRelIsoDB012_2011( ) muon_pog_PFRelIsoDB02_2012 = MuonPOGCorrections.make_muon_pog_PFRelIsoDB012_2012( ) #muon_pog_IsoMu24eta2p1_2011 = MuonPOGCorrections.make_muon_pog_IsoMu24eta2p1_2011() // This does not exist, yet :-) muon_pog_IsoMu24eta2p1_2012 = MuonPOGCorrections.make_muon_pog_IsoMu24eta2p1_2012( ) # Get object ID and trigger corrector functions def mc_corrector_2011(row): if row.run > 2: return 1 pu = pu_corrector(row.nTruePU) #pu = 1 m1id = muon_pog_PFTight_2011(row.mPt, row.mEta)
# Determine MC-DATA corrections is7TeV = bool('7TeV' in os.environ['jobid']) print "Is 7TeV:", is7TeV # Make PU corrector from expected data PU distribution # PU corrections .root files from pileupCalc.py pu_distributions = glob.glob(os.path.join( # 'inputs', os.environ['jobid'], 'data_TauPlusX*pu.root')) 'inputs', os.environ['jobid'], 'data_SingleMu*pu.root')) pu_corrector = PileupWeight.PileupWeight( 'S6' if is7TeV else 'S10', *pu_distributions) muon_pog_PFTight_2011 = MuonPOGCorrections.make_muon_pog_PFTight_2011() muon_pog_PFTight_2012 = MuonPOGCorrections.make_muon_pog_PFTight_2012() muon_pog_PFRelIsoDB02_2011 = MuonPOGCorrections.make_muon_pog_PFRelIsoDB012_2011() muon_pog_PFRelIsoDB02_2012 = MuonPOGCorrections.make_muon_pog_PFRelIsoDB012_2012() #muon_pog_IsoMu24eta2p1_2011 = MuonPOGCorrections.make_muon_pog_IsoMu24eta2p1_2011() // This does not exist, yet :-) muon_pog_IsoMu24eta2p1_2012 = MuonPOGCorrections.make_muon_pog_IsoMu24eta2p1_2012() # Get object ID and trigger corrector functions def mc_corrector_2011(row): if row.run > 2: return 1 pu = pu_corrector(row.nTruePU) #pu = 1 m1id = muon_pog_PFTight_2011(row.mPt, row.mEta) m1iso = muon_pog_PFRelIsoDB02_2011(row.mPt, row.mEta)