vbtfElectronIdFilter = Electron_IDFilter(vbtfelectronidWP95ps.ps()) ra3PhotonIdFilter = Photon_IDFilter(ra3photonidps.ps()) def addCutFlowMC(b): # b.AddWeightFilter("Weight", vertex_reweight) b.AddMuonFilter("PreCC", vbtfMuonId_cff) b.AddPhotonFilter("PreCC", ra3PhotonIdFilter) b.AddElectronFilter("PreCC", vbtfElectronIdFilter) b += cutTreeMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC = Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) outDir = "../results_" + strftime("%d_%b_%H") + "//NoSmear/" ensure_dir(outDir) from montecarlo.Summer11.QCD_Summer11_madgraph_All import * from montecarlo.Summer11.TTJets_TuneZ2_7TeV_madgraph_tauola_Summer11_PU_S4_START42_V11_v1_V15_03_19_jetCorrections_L1FastJet_L2Relative_L3Absolute_jetCollections_ak5calo_ak5pf_hbheNoiseFilterDefaultIsoReq_1 import * anal_ak5_caloMC.Run(outDir, conf_ak5_caloMC, [ TTJets_TuneZ2_7TeV_madgraph_tauola_Summer11_PU_S4_START42_V11_v1_V15_03_19_jetCorrections_L1FastJet_L2Relative_L3Absolute_jetCollections_ak5calo_ak5pf_hbheNoiseFilterDefaultIsoReq_1 ])
elif PUS == "PUS6" and DataSet == "MC": b.AddWeightFilter("Weight", vertex_reweight_PUS6) elif PUS != "PUS4" or PUS != "PUS6" or DataSet == "Data": print "Vertex reweighting is applied" if MuonForm == None: b.AddMuonFilter("PreCC", vbtfMuonId_cff) else: b.AddMuonFilter("PreCC", CustomMuID) b.AddPhotonFilter("PreCC", ra3PhotonIdFilter) b.AddElectronFilter("PreCC", vbtfElectronIdFilter) b += cutTreeMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC = Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) #outDir = "../results_"+strftime("%d_%b")+"//MuonNoSmearTT/" outDir = "../results_" + strftime("%d_%b") + ( "//%s%s_%s%s/" % (RunScriptName, SampleIndicated, PUS, ExtraInstruction)) ensure_dir(outDir) anal_ak5_caloMC.Run(outDir, conf_ak5_caloMC, sample)
cut_flow.TAttach(numComPhotons, oddPhoton) cut_flow.TAttach(oddPhoton, numComElectrons) cut_flow.TAttach(numComElectrons, oddElectron) cut_flow.TAttach(oddElectron, numComMuons) cut_flow.TAttach(numComMuons, oddMuon) cut_flow.TAttach(oddMuon, test) # ----------------------------------------------------------------------------- # Dataset tmp = PSet( Name="LM1", Events=100000, Format=("ICF", 2), File=[ "/Users/bainbrid/Desktop/work/data/LM1.root", ], CrossSection=4.888, #FirstEntry = 1, LastEntry=1000, ) # ----------------------------------------------------------------------------- # Analysis anal = Analysis("Ratio") anal.AddPhotonFilter("PreCC", phot_filter) anal.AddElectronFilter("PreCC", ele_filter) anal += cut_flow anal.Run(".", conf_ak5_calo, [tmp])
def addCutFlowBtagData(a) : a.AddPhotonFilter("PreCC",ra3PhotonIdFilter) a.AddElectronFilter("PreCC",vbtfElectronIdFilter) a.AddMuonFilter("PreCC",CustomMuID) a+=cutBtagTree # AK5 Calo conf_ak5_caloData = deepcopy(defaultConfig) conf_ak5_caloData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloData.XCleaning = deepcopy(default_cc) conf_ak5_caloData.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloData=Analysis("AK5Calo") addCutFlowData(anal_ak5_caloData) conf_ak5_caloBtagData = deepcopy(defaultConfig) conf_ak5_caloBtagData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloBtagData.XCleaning = deepcopy(default_cc) conf_ak5_caloBtagData.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloBtagData=Analysis("AK5Calo") addCutFlowBtagData(anal_ak5_caloBtagData) from data.Run2011.HT_Run2011AB import * from data.Run2011.MuHad_Run2011AB import * from HT_L1OffSet import *
a.AddWeightFilter("Weight", PreScaleWeights) # a.AddJetFilter("PreCC",JetCorrections) a.AddPhotonFilter("PreCC", ra3PhotonIdFilter) a.AddElectronFilter("PreCC", vbtfElectronIdFilter) a.AddMuonFilter("PreCC", vbtfMuonId_cff) a += cutTree # AK5 Calo conf_ak5_caloData = deepcopy(defaultConfig) conf_ak5_caloData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloData.XCleaning = deepcopy(default_cc) conf_ak5_caloData.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloData = Analysis("AK5Calo") addCutFlowData(anal_ak5_caloData) # AK5 PF conf_ak5_pfData = deepcopy(defaultConfig) conf_ak5_pfData.Ntuple = deepcopy(ak5_pf) conf_ak5_pfData.XCleaning = deepcopy(default_cc) conf_ak5_pfData.Common = deepcopy(default_common) anal_ak5_pfData = Analysis("AK5PF") addCutFlowData(anal_ak5_pfData) tedSkim = PSet(Name="TedSkim", Format=("ICF", 3), Weight=1.0, File="../resultsWithSingleTop//Data/AK5Calo_Jets.root")
cutTreeMC.TAttach(htCut350, TriggerPlots) from ra1objectid.vbtfElectronId_cff import * from ra1objectid.vbtfMuonId_cff import * from ra1objectid.ra3PhotonId_cff import * #JetSmear = JetSmear(0.1,30) vbtfElectronIdFilter = Electron_IDFilter(vbtfelectronidWP95ps.ps()) ra3PhotonIdFilter = Photon_IDFilter(ra3photonidps.ps()) def addCutFlowMC(b): b.AddPhotonFilter("PreCC", ra3PhotonIdFilter) b.AddElectronFilter("PreCC", vbtfElectronIdFilter) b += cutTreeMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC = Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) Triggersamples = [LM0, LM1, LM2, LM3, LM4] ensure_dir("../results/Triggers/") anal_ak5_caloMC.Run("../results/Triggers/", conf_ak5_caloMC, Triggersamples)
data = PSet( Name="Data", Format=("ICF", 3), File=[ "dcap://gfe02.grid.hep.ph.ic.ac.uk:22128//pnfs/hep.ph.ic.ac.uk/data/cms/store/user/georgia//ICF/automated/2011_04_09_18_07_10///SusyCAF_Tree_100_1_OjA.root", #"/vols/cms02/bainbrid/qcd/data/HT_Run2011A_PromptReco_v1.root", ], Weight=1., #FirstEntry = 1, #LastEntry = 10, ) from montecarlo.LMx import LM0 from montecarlo.QCD_Pythia6_384patch3_V14_00_02_ALL import * from data.Jet_35pb_WithTP_json051110 import Jet_35pb_WithTP_json051110 # ----------------------------------------------------------------------------- # Analysis anal = Analysis("Validation") anal.AddPhotonFilter("PreCC", phot_filter) anal.AddElectronFilter("PreCC", ele_filter) #anal.AddJetFilter("PreCC",JetCorrections) anal += cut_flow anal.Run(".", conf_ak5_calo, [data]) #anal.Run(".",conf_ak5_calo,[lm1]) #anal.Run(".",conf_ak5_calo,[LM0]) #anal.Run(".",conf_ak5_calo,[Jet_35pb_WithTP_json051110]) #anal.Run(".",conf_ak5_calo,QCD_Pythia6_384patch3_V14_00_02_ALL)
#JetSmear = JetSmear(0.1,30) vbtfElectronIdFilter = Electron_IDFilter( vbtfelectronidWP95ps.ps() ) ra3PhotonIdFilter = Photon_IDFilter( ra3photonidps.ps() ) def addCutFlowMC(b) : b.AddPhotonFilter("PreCC",ra3PhotonIdFilter) b.AddElectronFilter("PreCC",vbtfElectronIdFilter) b+=cutTreeMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC=Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) # AK5 PF conf_ak5_pfMC = deepcopy(defaultConfig) conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf) conf_ak5_pfMC.XCleaning = deepcopy(default_cc) conf_ak5_pfMC.Common = deepcopy(default_common) anal_ak5_pfMC=Analysis("AK5PF") addCutFlowMC(anal_ak5_pfMC) # AK5 JPT conf_ak5_jptMC = deepcopy(defaultConfig) conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)
conf = defaultConfig.copy() conf.Ntuple.Jets.Prefix = "ak5JetPF" conf.Common.Jets.EtCut = 20.0 conf.Common.Jets.EtaCut = 5.0 conf.Common.Muons.EtaCut = 2.4 conf.Common.Muons.PtCut = 15.0 conf.Common.Muons.TrkIsoCut = -1. conf.Common.Muons.CombIsoCut = 0.1 conf.Common.Electrons.PtCut = 15.0 conf.Common.Electrons.EtaCut = 2.4 conf.Common.Photons.EtCut = 30. conf.Common.Electrons.ApplyID = False # Create the analysis a = Analysis("KinCorrection") tree = Tree("Main") Trigger = PSet(McAnal=False, MSugraScan=0., TriggerBits=PSet(bit1="")) ssdlTrigg = SSDLTrigger("Trigger", Trigger.ps()) ZeroMuons = OP_NumComMuons("==", 0) KinCorrParPlus = PSet( zMassMin=71., zMassMax=111., BarrelEtaMax=1.4442, EndCapEtaMin=1.56, EndCapEtaMax=2.5, MinElecRescaledPt=0., c_ErNu_pt=0.,
b+=cutTreeMC def addCutFlowData(a): a.AddPhotonFilter("PreCC",ra3PhotonIdFilter) a.AddElectronFilter("PreCC",vbtfElectronIdFilter) a+=cutTreeData #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) # conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC=Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) #AK5 Calo conf_ak5_caloData = deepcopy(defaultConfig) conf_ak5_caloData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloData.XCleaning = deepcopy(default_cc) conf_ak5_caloData.Common = deepcopy(default_common) # conf_ak5_caloMC.Common.print_out() anal_ak5_caloData=Analysis("AK5Calo") addCutFlowData(anal_ak5_caloData) testFile = PSet( Name = "TestFile",
# JESUncert(-0.02,true) means -2% per unit eta JESUncert = JESUncert(-0.02, True) def addCutFlowMC(b): b.AddJetFilter("PreCC", JESUncert) b += cutTreeMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC = Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) # AK5 PF conf_ak5_pfMC = deepcopy(defaultConfig) conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf) conf_ak5_pfMC.XCleaning = deepcopy(default_cc) conf_ak5_pfMC.Common = deepcopy(default_common) anal_ak5_pfMC = Analysis("AK5PF") addCutFlowMC(anal_ak5_pfMC) # AK5 JPT conf_ak5_jptMC = deepcopy(defaultConfig) conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)
def addCutFlowMC(b): b.AddPhotonFilter("PreCC", ra3PhotonIdFilter) b.AddElectronFilter("PreCC", vbtfElectronIdFilter) b.AddJetFilter("PreCC", JESUncert) b += cutTreeMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC = Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) # AK5 PF conf_ak5_pfMC = deepcopy(defaultConfig) conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf) conf_ak5_pfMC.XCleaning = deepcopy(default_cc) conf_ak5_pfMC.Common = deepcopy(default_common) anal_ak5_pfMC = Analysis("AK5PF") addCutFlowMC(anal_ak5_pfMC) # AK5 JPT conf_ak5_jptMC = deepcopy(defaultConfig) conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)
a += kinsuitecomplot a += myplots a += HadAlphaT # ----------------------------------------------------------------------------- # Analyses # AK5 Calo conf_ak5_calo = deepcopy(defaultConfig) conf_ak5_calo.Ntuple = deepcopy(ak5_calo) conf_ak5_calo.XCleaning = deepcopy(default_cc) conf_ak5_calo.Common = deepcopy(default_common) anal_ak5_calo = Analysis("AK5Calo") addCutFlow(anal_ak5_calo) # AK5 JPT conf_ak5_jpt = deepcopy(defaultConfig) conf_ak5_jpt.Ntuple = deepcopy(ak5_jpt) conf_ak5_jpt.XCleaning = deepcopy(default_cc) conf_ak5_jpt.Common = deepcopy(default_common) anal_ak5_jpt = Analysis("AK5JPT") addCutFlow(anal_ak5_jpt) # Run analyses #anal_ak5_calo.Run("results",conf_ak5_calo,[lm0])
lm1 = PSet( Name="lm1", Events=100000, Format=("ICF", 2), File=[ #"dcap://gfe02.grid.hep.ph.ic.ac.uk:22128//pnfs/hep.ph.ic.ac.uk/data/cms/store/user/bainbrid/ICF/automated/2010_07_12_17_52_54/LM1.Spring10-START3X_V26_S09-v1.GEN-SIM-RECO/SusyCAF_Tree_1_1.root", "/Users/bainbrid/Desktop/work/data/LM1.root", ], CrossSection=4.888, #FirstEntry = 1, #LastEntry = 10, ) from data.Run2011.HT_Run2011A_PromptReco_v1 import HT_Run2011A_PromptReco_v1 from data.Run2011.HT_Run2011_promptReco_DCS import HT_Run2011_promptReco_DCS # ----------------------------------------------------------------------------- # Analysis anal = Analysis("TurnOn") anal.AddMuonFilter("PreCC", vbtf_muon) anal.AddElectronFilter("PreCC", vbtf_electron) anal.AddPhotonFilter("PreCC", ra3_photon) anal += cut_flow #anal.Run(".",conf,[data]) #anal.Run(".",conf,[tmp]) #anal.Run(".",conf,[lm1]) #anal.Run(".",conf_ak5_calo,[HT_Run2011A_PromptReco_v1]) anal.Run(".", conf_ak5_calo, [HT_Run2011_promptReco_DCS])
from montecarlo.Spring11.TToBLNu_TuneZ2_tW_channel_7TeV_madgraph_Spring11_PU_S1_START311_V1G1_v1 import * from montecarlo.Spring11.TToBLNu_TuneZ2_s_channel_7TeV_madgraph_Spring11_PU_S1_START311_V1G1_v1 import * from montecarlo.Spring11.LMx_SUSY_sftsht_7TeV_pythia6_Spring11_PU_S1_START311_V1G1_v1 import * from data.Jet_35pb_WithTP_json051110 import Jet_35pb_WithTP_json051110 from data.Run2011.HT42_incomplete import HT42_incomplete from data.Run2011.HT_Run2011_promptReco_DCS import HT_Run2011_promptReco_DCS from data.Run2011.HT_Run2011A import HT_Run2011A from data.Run2011.HT_Run2011A_PromptReco_v1 import HT_Run2011A_PromptReco_v1 from data.Run2011.HT_Run2011A_AllReco_17June import HT_Run2011A_AllReco_17June from data.Run2011.HT_Run2011A_L1OffSet import HT_Run2011A_L1OffSet # ----------------------------------------------------------------------------- # Analysis anal=Analysis("Ratio") if ( not DataOption ) : anal.AddWeightFilter("Weight", vertex_reweight) anal.AddMuonFilter("PreCC",vbtf_muon) anal.AddElectronFilter("PreCC",vbtf_electron) anal.AddPhotonFilter("PreCC",ra3_photon) #anal.AddJetFilter("PreCC",JetCorrections) anal+=cut_flow if ( dataset == 0 ) : # Data 2010 anal.Run(output_dir,conf_ak5_calo,[Jet_35pb_WithTP_json051110]) # 160 files elif ( dataset == 1 ) : # Data 2011 #anal.Run(output_dir,conf_ak5_calo,[HT_Run2011_promptReco_DCS]) # 1543 files (and growing) #anal.Run(output_dir,conf_ak5_calo,[HT42_incomplete]) # 2579 files (and growing) anal.Run(output_dir,conf_ak5_calo,[HT_Run2011A]) # 2579 files (and growing)
a+=plots_hadronic a+=plots_comjet a+=plots_objkin a+=plots_kinsuite # ----------------------------------------------------------------------------- # Definition of analyses # Analyses # IC5 Calo conf_ic5_calo = deepcopy(defaultConfig) conf_ic5_calo.Ntuple = deepcopy(ic5_calo) conf_ic5_calo.XCleaning = deepcopy(default_cc) conf_ic5_calo.Common = deepcopy(default_common) # conf_ic5_calo.Common.print_out() anal_ic5_calo=Analysis("IC5Calo") addCutFlow(anal_ic5_calo) # AK5 Calo conf_ak5_calo = deepcopy(defaultConfig) conf_ak5_calo.Ntuple = deepcopy(ak5_calo) conf_ak5_calo.XCleaning = deepcopy(default_cc) conf_ak5_calo.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_calo=Analysis("AK5Calo") addCutFlow(anal_ak5_calo) # AK5 JPT conf_ak5_jpt = deepcopy(defaultConfig)
# JESUncert(0.1,false) means +10%, independent of eta # JESUncert(-0.1,false) means -10%, independent of eta # JESUncert(0.02,true) means +2% per unit eta # JESUncert(-0.02,true) means -2% per unit eta JESUncert = JESUncert(-0.1,False) def addCutFlowMC(b) : b.AddJetFilter("PreCC",JESUncert) b+=cutTreeMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC=Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) # AK5 PF conf_ak5_pfMC = deepcopy(defaultConfig) conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf) conf_ak5_pfMC.XCleaning = deepcopy(default_cc) conf_ak5_pfMC.Common = deepcopy(default_common) anal_ak5_pfMC=Analysis("AK5PF") addCutFlowMC(anal_ak5_pfMC) # AK5 JPT conf_ak5_jptMC = deepcopy(defaultConfig) conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)
#"/Users/bainbrid/Desktop/work/data/HT_Run2011A_PromptReco_v1.root", ], Weight=1., LastEntry=10, ) from data.Run2011.HT_Run2011A_PromptReco_v1 import HT_Run2011A_PromptReco_v1 from data.Run2011.HT_Run2011_promptReco_DCS import HT_Run2011_promptReco_DCS from data.Run2011.HT42_incomplete import HT42_incomplete from data.Run2011.HT_Run2011A_L1OffSet import HT_Run2011A_L1OffSet from data.Run2011.HT_Run2011A import HT_Run2011A from data.Run2011.HT_Run2011A_AllReco_17June import HT_Run2011A_AllReco_17June # ----------------------------------------------------------------------------- # Analysis anal = Analysis("Trigger") #anal.AddPhotonFilter("PreCC",phot_filter) #anal.AddElectronFilter("PreCC",ele_filter) #anal.AddJetFilter("PreCC",JetCorrections) anal += cut_flow #anal.Run(".",conf,[data]) #anal.Run(".",conf,[Jet_35pb_WithTP_json051110]) #anal.Run(".",conf_ak5_calo,[HT_Run2011A_PromptReco_v1]) #anal.Run(".",conf_ak5_calo,[HT_Run2011_promptReco_DCS]) #anal.Run(output_dir,conf_ak5_calo,[HT42_incomplete]) #anal.Run(output_dir,conf_ak5_calo,[HT_Run2011A_L1OffSet]) anal.Run(output_dir, conf_ak5_calo, [HT_Run2011A]) #anal.Run(output_dir,conf_ak5_calo,[HT_Run2011A_AllReco_17June])
def addCutFlowMC(b): b.AddJetFilter("PreCC", JetSmear) b += cutTreeMC # ----------------------------------------------------------------------------- # Definition of analyses # Analyses #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC = Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) # AK5 Calo conf_ak5_caloData = deepcopy(defaultConfig) conf_ak5_caloData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloData.XCleaning = deepcopy(default_cc) conf_ak5_caloData.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloData = Analysis("AK5Calo") addCutFlowData(anal_ak5_caloData) # AK5 PF conf_ak5_pfData = deepcopy(defaultConfig)
a += cutTreeData # ----------------------------------------------------------------------------- # Definition of analyses # Analyses # IC5 Calo conf_ic5_calo = deepcopy(defaultConfig) conf_ic5_calo.Ntuple = deepcopy(ic5_calo) conf_ic5_calo.XCleaning = deepcopy(default_cc) conf_ic5_calo.Common = deepcopy(default_common) # conf_ic5_calo.Common.print_out() anal_ic5_calo = Analysis("IC5Calo") addCutDataFlow(anal_ic5_calo) # AK5 Calo # conf_ak5_caloMC = deepcopy(defaultCon# fig) # conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) # conf_ak5_caloMC.XCleaning = deepcopy(default_cc) # conf_ak5_caloMC.Common = deepcopy(default_common) # # conf_ak5_calo.Common.print_out() # anal_ak5_caloMC=Analysis("AK5Calo") # addCutMCFlow(anal_ak5_caloMC) # # # AK5 Calo # # conf_ak5_caloData = deepcopy(defaultConfig)
cutTreeMC, junkVar, l = MakeMCTree(100. * (325. / 375.), Muon=True) vbtfElectronIdFilter = Electron_IDFilter(vbtfelectronidWP95ps.ps()) ra3PhotonIdFilter = Photon_IDFilter(ra3photonidps.ps()) CustomMuID = OL_CustomVBTFMuID(mu_id_lower.ps()) def addCutFlowMC(b): b.AddWeightFilter("Weight", vertex_reweight) b.AddMuonFilter("PreCC", CustomMuID) b.AddPhotonFilter("PreCC", ra3PhotonIdFilter) b.AddElectronFilter("PreCC", vbtfElectronIdFilter) b += cutTreeMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC = Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) outDir = "../results_" + strftime("%d_%b") + "//MuonNoSmear43/" ensure_dir(outDir) anal_ak5_caloMC.Run(outDir, conf_ak5_caloMC, L1OffSet_MC_Lower_Bins) #anal_ak5_caloMC.Run(outDir,conf_ak5_caloMC,Summer11_MC_Lower_Bins)
# a.AddWeightFilter("Weight",PreScaleWeights) # a.AddJetFilter("PreCC",JetCorrections) a.AddPhotonFilter("PreCC", ra3PhotonIdFilter) a.AddElectronFilter("PreCC", vbtfElectronIdFilter) a.AddMuonFilter("PreCC", vbtfMuonId_cff) a += cutTreeData # AK5 Calo conf_ak5_caloData = deepcopy(defaultConfig) conf_ak5_caloData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloData.XCleaning = deepcopy(default_cc) conf_ak5_caloData.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloData = Analysis("AK5Calo") addCutFlowData(anal_ak5_caloData) # AK5 PF conf_ak5_pfData = deepcopy(defaultConfig) conf_ak5_pfData.Ntuple = deepcopy(ak5_pf) conf_ak5_pfData.XCleaning = deepcopy(default_cc) conf_ak5_pfData.Common = deepcopy(default_common) anal_ak5_pfData = Analysis("AK5PF") addCutFlowData(anal_ak5_pfData) from data.Run2011.HT_Run2011_promptReco_DCS import * outDir = "../results/TriggerTurnOns/" ensure_dir(outDir) anal_ak5_caloData.Run(outDir, conf_ak5_caloData, [HT_Run2011_promptReco_DCS])
def addCutFlowData(a): a.AddPhotonFilter("PreCC", ra3PhotonIdFilter) a.AddElectronFilter("PreCC", vbtfElectronIdFilter) a.AddMuonFilter("PreCC", vbtfMuonId_cff) a += cutTree # AK5 Calo conf_ak5_caloData = deepcopy(defaultConfig) conf_ak5_caloData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloData.XCleaning = deepcopy(default_cc) conf_ak5_caloData.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloData = Analysis("AK5Calo") addCutFlowData(anal_ak5_caloData) from HT_L1OffSet import * from data.Run2011.HT_Run2011AB import * #from RA2b_Skims import * #outDir = "../Split_Jsons_"+strftime("%d_%b")+"/Data43/" #outDir = "../RA2B_"+strftime("%d_%b")+"/Data43/" outDir = "../results_" + strftime("%d_%b") + "//Data43" ensure_dir(outDir) #anal_ak5_caloData.Run(outDir,conf_ak5_caloData,[RA2B_Skim])#HT_Run2011A]) #anal_ak5_caloData.Run(outDir,conf_ak5_caloData,[HT_Run2011AB]) anal_ak5_caloData.Run(outDir, conf_ak5_caloData, [HT_L1OffSet]) #HT_Run2011A])
def addCutFlowMC(b) : b+=cutTreeMC # ----------------------------------------------------------------------------- # Definition of analyses # Analyses # AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloMC=Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) # AK5 Calo conf_ak5_caloData = deepcopy(defaultConfig) conf_ak5_caloData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloData.XCleaning = deepcopy(default_cc) conf_ak5_caloData.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloData=Analysis("AK5Calo") addCutFlowData(anal_ak5_caloData) # ----------------------------------------------------------------------------- # Run analyses #IMPERIAL COLLEGE
a.AddWeightFilter("Weight", PreScaleWeights) # a.AddJetFilter("PreCC",JetCorrections) a.AddPhotonFilter("PreCC", ra3PhotonIdFilter) a.AddElectronFilter("PreCC", vbtfElectronIdFilter) a.AddMuonFilter("PreCC", vbtfMuonId_cff) a += cutTree # AK5 Calo conf_ak5_caloData = deepcopy(defaultConfig) conf_ak5_caloData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloData.XCleaning = deepcopy(default_cc) conf_ak5_caloData.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloData = Analysis("AK5Calo") addCutFlowData(anal_ak5_caloData) # AK5 PF conf_ak5_pfData = deepcopy(defaultConfig) conf_ak5_pfData.Ntuple = deepcopy(ak5_pf) conf_ak5_pfData.XCleaning = deepcopy(default_cc) conf_ak5_pfData.Common = deepcopy(default_common) anal_ak5_pfData = Analysis("AK5PF") addCutFlowData(anal_ak5_pfData) tedSkim = PSet(Name="TedSkim", Format=("ICF", 3), Weight=1.0, File="../results_" + strftime("%d_%b_%H") +
b+=cutTreeMC def addCutFlowBtagMC(a) : a.AddWeightFilter("Weight",vertex_reweight) a.AddPhotonFilter("PreCC",ra3PhotonIdFilter) a.AddElectronFilter("PreCC",vbtfElectronIdFilter) a.AddMuonFilter("PreCC",CustomMuID) a+=cutTreeBtagMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) #conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC=Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) conf_ak5_caloBtagData = deepcopy(defaultConfig) conf_ak5_caloBtagData.Ntuple = deepcopy(ak5_calo) conf_ak5_caloBtagData.XCleaning = deepcopy(default_cc) conf_ak5_caloBtagData.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_caloBtagData=Analysis("AK5Calo") addCutFlowBtagMC(anal_ak5_caloBtagData) outDir = "../results_"+strftime("%d_%b")+"//MuonNoSmear_Btag/" #outDir = "../results_"+strftime("%d_%b")+"//MuonNoSmear/" ensure_dir(outDir)
def addCutFlowMC(b) : b.AddJetFilter("PreCC",JetAdd) # b.AddWeightFilter("Weight", vertex_reweight) b.AddMuonFilter("PreCC",vbtfMuonId_cff) b.AddPhotonFilter("PreCC",ra3PhotonIdFilter) b.AddElectronFilter("PreCC",vbtfElectronIdFilter) b.AddJetFilter("PreCC",JESUncert) b+=cutTreeMC #AK5 Calo conf_ak5_caloMC = deepcopy(defaultConfig) conf_ak5_caloMC.Ntuple = deepcopy(ak5_calo) conf_ak5_caloMC.XCleaning = deepcopy(default_cc) conf_ak5_caloMC.Common = deepcopy(default_common) conf_ak5_caloMC.Common.print_out() anal_ak5_caloMC=Analysis("AK5Calo") addCutFlowMC(anal_ak5_caloMC) # AK5 PF conf_ak5_pfMC = deepcopy(defaultConfig) conf_ak5_pfMC.Ntuple = deepcopy(ak5_pf) conf_ak5_pfMC.XCleaning = deepcopy(default_cc) conf_ak5_pfMC.Common = deepcopy(default_common) anal_ak5_pfMC=Analysis("AK5PF") addCutFlowMC(anal_ak5_pfMC) # AK5 JPT conf_ak5_jptMC = deepcopy(defaultConfig) conf_ak5_jptMC.Ntuple = deepcopy(ak5_jpt)
a+=greaterThreeJet a+=PlotsMetG3 a+=plots_comjet a+=plots_kinsuite # ----------------------------------------------------------------------------- # Definition of analyses # Analyses # IC5 Calo conf_ic5_calo = deepcopy(defaultConfig) conf_ic5_calo.Ntuple = deepcopy(ic5_calo) conf_ic5_calo.XCleaning = deepcopy(default_cc) conf_ic5_calo.Common = deepcopy(default_common) # conf_ic5_calo.Common.print_out() anal_ic5_calo=Analysis("IC5Calo") addCutFlow(anal_ic5_calo) # AK5 Calo conf_ak5_calo = deepcopy(defaultConfig) conf_ak5_calo.Ntuple = deepcopy(ak5_calo) conf_ak5_calo.XCleaning = deepcopy(default_cc) conf_ak5_calo.Common = deepcopy(default_common) # conf_ak5_calo.Common.print_out() anal_ak5_calo=Analysis("AK5Calo") addCutFlow(anal_ak5_calo) # AK5 JPT conf_ak5_jpt = deepcopy(defaultConfig)