From 6f368ff1c0a923f664af125bc78e48de81920b5e Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Fri, 28 Jun 2013 08:35:49 -0500 Subject: [PATCH 01/34] first update of ewk17 --- samples/__ewk17__.py | 33 ++++++++++++--------------------- 1 file changed, 12 insertions(+), 21 deletions(-) diff --git a/samples/__ewk17__.py b/samples/__ewk17__.py index 976787c..d008bd0 100644 --- a/samples/__ewk17__.py +++ b/samples/__ewk17__.py @@ -3,37 +3,28 @@ pnfs = pnfs() ewk17 = SampleHolder() -four = "WJetsToLNu_HT-400ToInf_8TeV-madgraph.Summer12-PU_S7_START52_V9-v1.AODSIM/" -incl = "WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball.Summer12-PU_S7_START52_V9-v1.AODSIM/" -ewk17.add("wj_lv_mg_ht_0_250_other_reqs", '%s/dburton/WJet_Skim_v2/")'%pnfs, xs = 36140.2*4135747.0/18273090.0) #wrong XS! -ewk17.add("wj_lv_mg_ht_250_300", '%s/clucas/ICF/automated/2012_06_20_16_29_08/")'%pnfs, xs = 48.01) -ewk17.add("wj_lv_mg_ht_300_400", '%s/clucas/ICF/automated/2012_06_21_02_04_54/")'%pnfs, xs = 38.30) -ewk17.add("wj_lv_mg_ht_400_inf", '%s/clucas/ICF/automated/2012_06_05_09_30_46/%s")'%(pnfs, four), xs = 25.22) -ewk17.add("wj_lv_mg_ht_incl", '%s/clucas/ICF/automated/2012_06_05_09_30_46/%s")'%(pnfs, incl), xs = {"LO":30400.0, "NLO":36257.2}["NLO"]) - -ewk17.add("zinv_mg_ht_50_100.job214", '%s/clucas//ICF/automated/2012_06_10_22_27_16/")'%pnfs, xs = {"NNLO":452.75, "LO":381.2}["NNLO"]) -ewk17.add("zinv_mg_ht_100_200.job234", '%s/clucas//ICF/automated/2012_06_17_22_29_52/")'%pnfs, xs = {"NNLO":190.39, "LO":160.3}["NNLO"]) -ewk17.add("zinv_mg_ht_200_400.job233", '%s/clucas//ICF/automated/2012_06_15_00_32_54/")'%pnfs, xs = {"NNLO":49.2776, "LO":41.49}["NNLO"]) -ewk17.add("zinv_mg_ht_400_inf.job213", '%s/clucas//ICF/automated/2012_06_10_22_08_24/")'%pnfs, xs = {"NNLO":6.2639, "LO":5.274}["NNLO"]) +### Wjets ### -vv = '%s/clucas//ICF/automated/2012_06_05_09_30_46/'%pnfs -ewk17.add("ww_py.job188", '%s/WW_TuneZ2star_8TeV_pythia6_tauola.Summer12-PU_S7_START52_V9-v1.AODSIM/")'%vv, xs = {"NNLO":57.1097, "LO":33.61}["NNLO"]) -ewk17.add("wz_py.job188", '%s/WZ_TuneZ2star_8TeV_pythia6_tauola.Summer12-PU_S7_START52_V9-v1.AODSIM/")'%vv, xs = {"NNLO":32.3161, "LO":12.63}["NNLO"]) -ewk17.add("zz_py.job188", '%s/ZZ_TuneZ2star_8TeV_pythia6_tauola.Summer12-PU_S7_START52_V9-v1.AODSIM/")'%vv, xs = {"NNLO":8.25561, "LO":5.196}["NNLO"]) - -#https://twiki.cern.ch/twiki/bin/view/CMS/HiggsMCProductionSummer12#VBF_H_bb_POWHEG -ewk17.add("zinv_hbb_125_powheg.job342", '%s/yeshaq//ICF/automated/2012_08_31_15_26_30/")'%pnfs, xs = 0.0361) +ewk17.add("wj_lv_mg_ht_0_250_other_reqs", '%s/dburton/WJet_Skim_v2/")'%pnfs, xs = 36140.2*4135747.0/18273090.0) #wrong XS! S10_incl = "WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" S10_excl = "/WJetsToLNu_HT-%s_8TeV-madgraph_v2.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" - ewk17.add("wj_lv_mg_ht_250_300.job498", '%s/karage//ICF/automated/2012_11_29_20_11_51/%s")'%(pnfs, S10_excl)%"250To300", xs = {"LO":48.01, "NLO":57.26}["NLO"]) ewk17.add("wj_lv_mg_ht_300_400.job498", '%s/karage//ICF/automated/2012_11_29_20_11_51/%s")'%(pnfs, S10_excl)%"300To400", xs = {"LO":38.30, "NLO":45.68}["NLO"]) ewk17.add("wj_lv_mg_ht_400_inf.job498", '%s/karage//ICF/automated/2012_11_29_20_11_51/%s")'%(pnfs, S10_excl)%"400ToInf", xs = {"LO":25.22, "NLO":30.08}["NLO"]) ewk17.add("wj_lv_mg_ht_incl.job363", '%s/clucas/ICF//automated/2012_09_21_09_36_56/%s")'%(pnfs, S10_incl), xs = {"LO":30400.0, "NLO":36257.2}["NLO"]) +ewk17.add("wj_lv_mg_ht_5_250_skim", '%s/yeshaq/ICF/supy-output/wj_skim/")'%pnfs, xs = {"LO":30400.0, "NLO":36257.2}["NLO"]) + +for part in [1,2,3,4,5] : + ewk17.add("wj_lv_mg_ht_incl_v2.job673_part%i"%part, '%s/yeshaq//ICF/automated/2013_05_23_14_13_07/part%i")'%(pnfs,part), xs = {"LO":30400.0, "NLO":36257.2}["NLO"]) -Zinv_S10_excl = "ZJetsToNuNu_%s_HT_%s_TuneZ2Star_8TeV_madgraph.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" +### Zinv ### + +#https://twiki.cern.ch/twiki/bin/view/CMS/HiggsMCProductionSummer12#VBF_H_bb_POWHEG +ewk17.add("zinv_hbb_125_powheg.job342", '%s/yeshaq//ICF/automated/2012_08_31_15_26_30/")'%pnfs, xs = 0.0361) + +Zinv_S10_excl = "ZJetsToNuNu_%s_HT_%s_TuneZ2Star_8TeV_madgraph.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" ewk17.add("zinv_mg_ht_50_100.job407", '%s/clucas//ICF/automated/2012_09_23_19_53_51/")'%pnfs, xs = {"NNLO":452.75, "LO":381.2}["NNLO"]) ewk17.add("zinv_mg_ht_100_200.job365", '%s/clucas//ICF/automated/2012_09_21_10_32_59/%s")'%(pnfs, Zinv_S10_excl)%(100,200), xs = {"NNLO":190.39, "LO":160.3}["NNLO"]) ewk17.add("zinv_mg_ht_200_400.job365", '%s/clucas//ICF/automated/2012_09_21_10_32_59/%s")'%(pnfs, Zinv_S10_excl)%(200,400), xs = {"NNLO":49.2776, "LO":41.49}["NNLO"]) From 3cfccc3501581f302c13bb717de0df85ad459150 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Fri, 28 Jun 2013 09:40:49 -0500 Subject: [PATCH 02/34] update of ewk17 --- samples/__ewk17__.py | 25 +++++++++++++++++++++++-- 1 file changed, 23 insertions(+), 2 deletions(-) diff --git a/samples/__ewk17__.py b/samples/__ewk17__.py index d008bd0..e4cd929 100644 --- a/samples/__ewk17__.py +++ b/samples/__ewk17__.py @@ -4,20 +4,37 @@ ewk17 = SampleHolder() ### Wjets ### - +ewk17.add("wj_lv_mg_ht_5_250_skim", '%s/yeshaq/ICF/supy-output/wj_skim/")'%pnfs, xs = {"LO":30400.0, "NLO":36257.2}["NLO"]) ewk17.add("wj_lv_mg_ht_0_250_other_reqs", '%s/dburton/WJet_Skim_v2/")'%pnfs, xs = 36140.2*4135747.0/18273090.0) #wrong XS! S10_incl = "WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" S10_excl = "/WJetsToLNu_HT-%s_8TeV-madgraph_v2.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" + +ewk17.add("wj_lv_mg_ht_10_150_skim", '%s/yeshaq/ICF/supy-output/wj_skim_v2/")'%pnfs, xs = {"NLO":9090.43}["NLO"]) +ewk17.add("wj_lv_mg_ht_150_200.job663", '%s/agapitos//ICF/automated/2013_05_15_18_11_50/")'%pnfs, xs = 290.69475) # see comment 1 +ewk17.add("wj_lv_mg_ht_200_250.job672", '%s/agapitos//ICF/automated/2013_05_15_17_18_45/")'%pnfs, xs = 123.3417) # see comment 2 ewk17.add("wj_lv_mg_ht_250_300.job498", '%s/karage//ICF/automated/2012_11_29_20_11_51/%s")'%(pnfs, S10_excl)%"250To300", xs = {"LO":48.01, "NLO":57.26}["NLO"]) ewk17.add("wj_lv_mg_ht_300_400.job498", '%s/karage//ICF/automated/2012_11_29_20_11_51/%s")'%(pnfs, S10_excl)%"300To400", xs = {"LO":38.30, "NLO":45.68}["NLO"]) ewk17.add("wj_lv_mg_ht_400_inf.job498", '%s/karage//ICF/automated/2012_11_29_20_11_51/%s")'%(pnfs, S10_excl)%"400ToInf", xs = {"LO":25.22, "NLO":30.08}["NLO"]) ewk17.add("wj_lv_mg_ht_incl.job363", '%s/clucas/ICF//automated/2012_09_21_09_36_56/%s")'%(pnfs, S10_incl), xs = {"LO":30400.0, "NLO":36257.2}["NLO"]) -ewk17.add("wj_lv_mg_ht_5_250_skim", '%s/yeshaq/ICF/supy-output/wj_skim/")'%pnfs, xs = {"LO":30400.0, "NLO":36257.2}["NLO"]) + for part in [1,2,3,4,5] : ewk17.add("wj_lv_mg_ht_incl_v2.job673_part%i"%part, '%s/yeshaq//ICF/automated/2013_05_23_14_13_07/part%i")'%(pnfs,part), xs = {"LO":30400.0, "NLO":36257.2}["NLO"]) +#comment 1 +# CrossSection = ( NNLO_W_inclusive_from[1] / LO_W_inclusive_from_PREP[2] ) * LO_W_exclusive_from_PREP: [3] + +#[1] https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat8TeV +#[2] http://cms.cern.ch/iCMS/prep/requestmanagement?dsn=WJetsToLNu_TuneZ2Star_8TeV-madgraph*&campid=Summer12 +#[3] http://cms.cern.ch/iCMS/prep/requestmanagement?dsn=WJetsToLNu_HT*200*madgraph*&campid=Summer12 + +#comment 2 +#CrossSection = 123.3417 #NNLO correct with the disconnection correction factor see: https://twiki.cern.ch/twiki/pub/CMS/SusyRA1Material/FixSmapleProblems_24May2013.pdf, correction factor is 1.1074 +#[1] https://twiki.cern.ch/twiki/bin/viewauth/CMS/StandardModelCrossSectionsat8TeV +#[2] http://cms.cern.ch/iCMS/prep/requestmanagement?dsn=WJetsToLNu_TuneZ2Star_8TeV-madgraph*&campid=Summer12 +#[3] http://cms.cern.ch/iCMS/prep/requestmanagement?dsn=WJetsToLNu_HT*200*madgraph*&campid=Summer12 + ### Zinv ### @@ -30,3 +47,7 @@ ewk17.add("zinv_mg_ht_200_400.job365", '%s/clucas//ICF/automated/2012_09_21_10_32_59/%s")'%(pnfs, Zinv_S10_excl)%(200,400), xs = {"NNLO":49.2776, "LO":41.49}["NNLO"]) ewk17.add("zinv_mg_ht_400_inf.job365", '%s/clucas//ICF/automated/2012_09_21_10_32_59/%s")'%(pnfs, Zinv_S10_excl)%(400,"inf"), xs = {"NNLO":6.2639, "LO":5.274}["NNLO"]) +ZJets = [ ZJetsToNuNu_50_HT_100_TuneZ2Star_8TeV_madgraph_Summer12_DR53X_PU_S10_START53_V7A_v1_V17_5_taus_0_clucasJob407CombinkarageJob500, ZJetsToNuNu_100_HT_200_TuneZ2Star_8TeV_madgraph_Summer12_DR53X_PU_S10_START53_V7A_v1_V17\ + _5_taus_0_clucasJob365, ZJetsToNuNu_200_HT_400_TuneZ2Star_8TeV_madgraph_Summer12_DR53X_PU_S10_START53_V7A_v1_V17_5_taus_0_clucasJob365CombinkarageJob500, ZJetsToNuNu_400_HT_inf_TuneZ2Star_8TeV_madgraph_Summer12_DR53X_PU_S10_STA\ + RT53_V7A_v1_V17_5_taus_0_clucasJob365CombinkarageJob500 ] + From f6bb6082390d83fb4cdb70463166c03d54b7b20d Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Fri, 28 Jun 2013 09:51:29 -0500 Subject: [PATCH 03/34] update of ewk17 --- samples/__ewk17__.py | 7 +------ 1 file changed, 1 insertion(+), 6 deletions(-) diff --git a/samples/__ewk17__.py b/samples/__ewk17__.py index e4cd929..935c3ab 100644 --- a/samples/__ewk17__.py +++ b/samples/__ewk17__.py @@ -9,8 +9,7 @@ S10_incl = "WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" S10_excl = "/WJetsToLNu_HT-%s_8TeV-madgraph_v2.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" - -ewk17.add("wj_lv_mg_ht_10_150_skim", '%s/yeshaq/ICF/supy-output/wj_skim_v2/")'%pnfs, xs = {"NLO":9090.43}["NLO"]) +ewk17.add("wj_lv_mg_ht_10_150", '%s/yeshaq/ICF/supy-output/wj_skim_v2/")'%pnfs, xs = {"NLO":9090.43}["NLO"]) ewk17.add("wj_lv_mg_ht_150_200.job663", '%s/agapitos//ICF/automated/2013_05_15_18_11_50/")'%pnfs, xs = 290.69475) # see comment 1 ewk17.add("wj_lv_mg_ht_200_250.job672", '%s/agapitos//ICF/automated/2013_05_15_17_18_45/")'%pnfs, xs = 123.3417) # see comment 2 ewk17.add("wj_lv_mg_ht_250_300.job498", '%s/karage//ICF/automated/2012_11_29_20_11_51/%s")'%(pnfs, S10_excl)%"250To300", xs = {"LO":48.01, "NLO":57.26}["NLO"]) @@ -47,7 +46,3 @@ ewk17.add("zinv_mg_ht_200_400.job365", '%s/clucas//ICF/automated/2012_09_21_10_32_59/%s")'%(pnfs, Zinv_S10_excl)%(200,400), xs = {"NNLO":49.2776, "LO":41.49}["NNLO"]) ewk17.add("zinv_mg_ht_400_inf.job365", '%s/clucas//ICF/automated/2012_09_21_10_32_59/%s")'%(pnfs, Zinv_S10_excl)%(400,"inf"), xs = {"NNLO":6.2639, "LO":5.274}["NNLO"]) -ZJets = [ ZJetsToNuNu_50_HT_100_TuneZ2Star_8TeV_madgraph_Summer12_DR53X_PU_S10_START53_V7A_v1_V17_5_taus_0_clucasJob407CombinkarageJob500, ZJetsToNuNu_100_HT_200_TuneZ2Star_8TeV_madgraph_Summer12_DR53X_PU_S10_START53_V7A_v1_V17\ - _5_taus_0_clucasJob365, ZJetsToNuNu_200_HT_400_TuneZ2Star_8TeV_madgraph_Summer12_DR53X_PU_S10_START53_V7A_v1_V17_5_taus_0_clucasJob365CombinkarageJob500, ZJetsToNuNu_400_HT_inf_TuneZ2Star_8TeV_madgraph_Summer12_DR53X_PU_S10_STA\ - RT53_V7A_v1_V17_5_taus_0_clucasJob365CombinkarageJob500 ] - From fe2133dd8d27655c61437fe2ec324dd0bce001ba Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Fri, 28 Jun 2013 10:44:57 -0500 Subject: [PATCH 04/34] adding binned wjets samples --- analyses/hadronicLook.py | 45 +++++++++++++++++++++------------------- 1 file changed, 24 insertions(+), 21 deletions(-) diff --git a/analyses/hadronicLook.py b/analyses/hadronicLook.py index d7cf283..0f19872 100644 --- a/analyses/hadronicLook.py +++ b/analyses/hadronicLook.py @@ -48,19 +48,20 @@ def parameters(self): "rechitComp": "PF", }, }), - "nJetsMinMax": self.vary({#"ge2j": (2, None), + "nJetsMinMax": self.vary({"ge2j": (2, None), #"le3j": (2, 3), - "ge4j": (4, None), + #"ge4j": (4, None), }), - "nBTagJets": self.vary({#"eq0b": (0, 0), + #"nBTagJets": self.vary({#"eq0b": (0, 0), #"eq1b": (1, 1), - "eq2b": (2, 2), + #"eq2b": (2, 2), #"eq3b": (3, 3), #"ge4b": (4, None), - }), - "thresholds": self.vary({#"275s": (275.0, 325.0, 73.3, 36.7), + #}), + "thresholds": self.vary({"200s": (200.0, None, 73.3, 36.7), + #"275s": (275.0, 325.0, 73.3, 36.7), #"325s": (325.0, 375.0, 86.7, 43.3), - "375": (375.0, None, 100.0, 50.0), + #"375": (375.0, None, 100.0, 50.0), #"675": (675.0, None, 100.0, 50.0), #"875": (875.0, None, 100.0, 50.0), }), @@ -225,8 +226,8 @@ def stepsTrigger(self, params): ]+([] if params["thresholds"][1]!=325.0 else [steps.trigger.lowestUnPrescaledTriggerFilter().onlyData()]) #apply trigger in lowest HT bin def stepsGenValidation(self) : - return [supy.steps.histos.histogrammer("genpthat",200,0,2000,title=";#hat{p_{T}} (GeV);events / bin").onlySim(), - #supy.steps.histos.histogrammer("genPartonHT",200,0,1000,title=";parton H_{T} (GeV);events / bin").onlySim(), + return [#supy.steps.histos.histogrammer("genpthat",200,0,2000,title=";#hat{p_{T}} (GeV);events / bin").onlySim(), + supy.steps.histos.histogrammer("genPartonHT",200,0,1000,title=";parton H_{T} (GeV);events / bin").onlySim(), ] def stepsEvent(self) : @@ -412,18 +413,18 @@ def listOfSteps(self, params) : return ([supy.steps.printer.progressPrinter()] + #self.stepsEventCount(params, label = "scanBefore") + - #self.stepsGenValidation() + + self.stepsGenValidation() + self.stepsEvent() + - self.stepsTrigger(params) + + #self.stepsTrigger(params) + self.stepsHtLeadingJets(params) + - self.stepsXclean(params) + - self.stepsBtagJets(params) + + #self.stepsXclean(params) + + #self.stepsBtagJets(params) + #self.stepsPlotsOne(params) + self.stepsQcdRejection(params) + self.stepsPlotsTwo(params) + #self.stepsMbb(params) + #self.stepsDisplayer(params) + - self.stepsOptional(params) + + #self.stepsOptional(params) + #self.stepsHtBins(params) + #self.stepsEventCount(params, label = "scanAfter") + []) @@ -481,10 +482,12 @@ def g_jets_mg(eL) : def w_binned() : out = [] - #out += specify(names = "wj_lv_mg_ht_0_250_other_reqs", nFilesMax = 1, nEventsMax = 20000, color = r.kRed) - out += specify(names = "wj_lv_mg_ht_250_300", color = r.kBlue) - out += specify(names = "wj_lv_mg_ht_300_400", color = r.kGreen) - out += specify(names = "wj_lv_mg_ht_400_inf", color = r.kCyan) + #out += specify(names = "wj_lv_mg_ht_10_150", color = r.kBlue, nFilesMax = 1, nEventsMax = 20000) + out += specify(names = "wj_lv_mg_ht_150_200.job663", color = r.kGreen, nFilesMax = 1, nEventsMax = 20000) + out += specify(names = "wj_lv_mg_ht_200_250.job672", color = r.kCyan, nFilesMax = 1, nEventsMax = 20000) + out += specify(names = "wj_lv_mg_ht_250_300.job498", color = r.kOrange, nFilesMax = 1, nEventsMax = 20000) + out += specify(names = "wj_lv_mg_ht_300_400.job498", color = r.kViolet, nFilesMax = 1, nEventsMax = 20000) + out += specify(names = "wj_lv_mg_ht_400_inf.job498", color = r.kAzure, nFilesMax = 1, nEventsMax = 20000) return out def z_binned() : @@ -535,10 +538,10 @@ def sms() : return out return ( - data_53X() + + #data_53X() + #data_52X() + #data_52X_2b_skim() + - #w_binned() + + w_binned() + #z_binned() + #top() + #vv() + @@ -576,7 +579,7 @@ def md(x, y) : "tt_8_mg.job315", "ttz_8_mg.job269", "t_s_powheg.job200", "t_t_powheg.job187", "t_tw_powheg.job187", "tbar_t_powheg.job194", "tbar_tw_powheg.job187"]) org.mergeSamples(targetSpec = md({"name":"Z + jets", "color": r.kBlue}, mcOps), allWithPrefix = "zinv_mg_ht") - org.mergeSamples(targetSpec = md({"name":"W + jets", "color": r.kOrange-3}, mcOps), allWithPrefix = "wj_lv_mg_ht_") + #org.mergeSamples(targetSpec = md({"name":"W + jets", "color": r.kOrange-3}, mcOps), allWithPrefix = "wj_lv_mg_ht_") org.mergeSamples(targetSpec = md({"name":"VV", "color": r.kOrange+3}, mcOps), sources = ["ww_py.job188", "wz_py.job188", "zz_py.job188"]) org.mergeSamples(targetSpec = md({"name":"ZH", "color":r.kMagenta}, mcOps), sources = ["zinv_hbb_125_powheg.job342"]) org.mergeSamples(targetSpec = md({"name":"LM6", "color":r.kMagenta}, mcOps), allWithPrefix = "lm6") From c6149067f2244ec3c170390df46a3c2f9742da76 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Fri, 28 Jun 2013 10:45:24 -0500 Subject: [PATCH 05/34] update of ewk17 --- samples/__ewk17__.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/samples/__ewk17__.py b/samples/__ewk17__.py index 935c3ab..53eb12d 100644 --- a/samples/__ewk17__.py +++ b/samples/__ewk17__.py @@ -9,7 +9,7 @@ S10_incl = "WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" S10_excl = "/WJetsToLNu_HT-%s_8TeV-madgraph_v2.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" -ewk17.add("wj_lv_mg_ht_10_150", '%s/yeshaq/ICF/supy-output/wj_skim_v2/")'%pnfs, xs = {"NLO":9090.43}["NLO"]) +ewk17.add("wj_lv_mg_ht_10_150", '%s/yeshaq/ICF/supy-output/wj_skim_v2")'%pnfs, xs = {"NLO":9090.43}["NLO"]) ewk17.add("wj_lv_mg_ht_150_200.job663", '%s/agapitos//ICF/automated/2013_05_15_18_11_50/")'%pnfs, xs = 290.69475) # see comment 1 ewk17.add("wj_lv_mg_ht_200_250.job672", '%s/agapitos//ICF/automated/2013_05_15_17_18_45/")'%pnfs, xs = 123.3417) # see comment 2 ewk17.add("wj_lv_mg_ht_250_300.job498", '%s/karage//ICF/automated/2012_11_29_20_11_51/%s")'%(pnfs, S10_excl)%"250To300", xs = {"LO":48.01, "NLO":57.26}["NLO"]) From fcf04e1c4018d6faac2186ffede5d7d4ce9a649b Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Sat, 29 Jun 2013 11:11:06 -0500 Subject: [PATCH 06/34] don't prune skimmed sample --- samples/__ewk17__.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/samples/__ewk17__.py b/samples/__ewk17__.py index 53eb12d..7d40ce2 100644 --- a/samples/__ewk17__.py +++ b/samples/__ewk17__.py @@ -9,7 +9,7 @@ S10_incl = "WJetsToLNu_TuneZ2Star_8TeV-madgraph-tarball.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" S10_excl = "/WJetsToLNu_HT-%s_8TeV-madgraph_v2.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM" -ewk17.add("wj_lv_mg_ht_10_150", '%s/yeshaq/ICF/supy-output/wj_skim_v2")'%pnfs, xs = {"NLO":9090.43}["NLO"]) +ewk17.add("wj_lv_mg_ht_10_150", '%s/yeshaq/ICF/supy-output/wj_skim_v2", pruneList = False)'%pnfs, xs = {"NLO":9090.43}["NLO"]) ewk17.add("wj_lv_mg_ht_150_200.job663", '%s/agapitos//ICF/automated/2013_05_15_18_11_50/")'%pnfs, xs = 290.69475) # see comment 1 ewk17.add("wj_lv_mg_ht_200_250.job672", '%s/agapitos//ICF/automated/2013_05_15_17_18_45/")'%pnfs, xs = 123.3417) # see comment 2 ewk17.add("wj_lv_mg_ht_250_300.job498", '%s/karage//ICF/automated/2012_11_29_20_11_51/%s")'%(pnfs, S10_excl)%"250To300", xs = {"LO":48.01, "NLO":57.26}["NLO"]) From f96bf1fa1823b06c222da240495465210dca675b Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Sat, 29 Jun 2013 11:37:12 -0500 Subject: [PATCH 07/34] adding HTMHTParked samples --- analyses/hadronicLook.py | 40 ++++++++++++++-------------------------- samples/__ht17__.py | 15 +++------------ 2 files changed, 17 insertions(+), 38 deletions(-) diff --git a/analyses/hadronicLook.py b/analyses/hadronicLook.py index 0f19872..d76fa8c 100644 --- a/analyses/hadronicLook.py +++ b/analyses/hadronicLook.py @@ -439,24 +439,12 @@ def listOfSamples(self,params) : from supy.samples import specify def data_53X() : - jw2012 = calculables.other.jsonWeight("cert/Cert_190456-203002_8TeV_PromptReco_Collisions12_JSON.txt") + jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt") out = [] - #out += specify(names = "HTMHT.Run2012B-13Jul2012-v1.AOD.job358", weights = jw2012)#, nFilesMax = 1, nEventsMax = 1000) - #out += specify(names = "HTMHT.Run2012C-24Aug2012-v1.AOD.job361", weights = jw2012)#, nFilesMax = 1, nEventsMax = 1000) - #out += specify(names = "HTMHT.Run2012C-PromptReco-v2.AOD.job360", weights = jw2012)#, nFilesMax = 1, nEventsMax = 1000) - #out += specify(names = "HT.Run2012A-13Jul2012-v1.AOD.job358", weights = jw2012)#, nFilesMax = 1, nEventsMax = 1000) - #out += specify(names = "HT.Run2012A-recover-06Aug2012-v1.AOD.job359", weights = jw2012)#, nFilesMax = 1, nEventsMax = 1000) - out += specify(names = "375_ge2b") - return out - - def data_52X() : - jw2012 = calculables.other.jsonWeight("cert/Cert_190456-196531_8TeV_PromptReco_Collisions12_JSON.txt") - - out = [] - out += specify(names = "HT.Run2012A-PromptReco-v1.AOD.job229", weights = jw2012, overrideLumi = 707.3810) - out += specify(names = "HTMHT.Run2012B-PromptReco-v1.AOD.job228", weights = jw2012, overrideLumi = 3354.0000) - out += specify(names = "HTMHT.Run2012B-PromptReco-v1.AOD.job238", weights = jw2012, overrideLumi = 923.7680) + out += specify(names = "HTMHTParked.Run2012B-22Jan2013-v1.job649", weights = jw2012, nFilesMax = 1, nEventsMax = 20000) + out += specify(names = "HTMHTParked.Run2012C-22Jan2013-v1.job649", weights = jw2012, nFilesMax = 1, nEventsMax = 20000) + out += specify(names = "HTMHTParked.Run2012D-22Jan2013-v1.job649", weights = jw2012, nFilesMax = 1, nEventsMax = 20000) return out def data_52X_2b_skim() : @@ -482,7 +470,7 @@ def g_jets_mg(eL) : def w_binned() : out = [] - #out += specify(names = "wj_lv_mg_ht_10_150", color = r.kBlue, nFilesMax = 1, nEventsMax = 20000) + out += specify(names = "wj_lv_mg_ht_10_150", color = r.kBlue, nFilesMax = 1, nEventsMax = 20000) out += specify(names = "wj_lv_mg_ht_150_200.job663", color = r.kGreen, nFilesMax = 1, nEventsMax = 20000) out += specify(names = "wj_lv_mg_ht_200_250.job672", color = r.kCyan, nFilesMax = 1, nEventsMax = 20000) out += specify(names = "wj_lv_mg_ht_250_300.job498", color = r.kOrange, nFilesMax = 1, nEventsMax = 20000) @@ -499,8 +487,10 @@ def z_binned() : return out def w_inclusive() : - return specify(names = "wj_lv_mg_ht_incl", color = r.kOrange) - + out = [] + for part in range(1,6): + out += specify(names = "wj_lv_mg_ht_incl_v2.job673_part%i"%part, nFilesMax = 1, nEventsMax = 1000) + return out def vv() : out = [] out += specify("ww_py.job188") @@ -538,16 +528,15 @@ def sms() : return out return ( - #data_53X() + - #data_52X() + + data_53X() + #data_52X_2b_skim() + - w_binned() + + #w_binned() + #z_binned() + #top() + #vv() + #qcd_py6(30.0e3) + #qcd_b_py6(30.0e3) + - ##w_inclusive() + + #w_inclusive() + #sms() + [] ) @@ -557,7 +546,6 @@ def md(x, y) : x.update(y) return x - #org.mergeSamples(targetSpec = {"name":"2012 Data", "color":r.kBlack, "markerStyle":20}, allWithPrefix = "HT") org.mergeSamples(targetSpec = {"name":"2012 Data", "color":r.kBlack, "markerStyle":20}, allWithPrefix = "HT") mcOps = {"markerStyle":1, "lineWidth":3, "goptions":"hist"} @@ -579,7 +567,7 @@ def md(x, y) : "tt_8_mg.job315", "ttz_8_mg.job269", "t_s_powheg.job200", "t_t_powheg.job187", "t_tw_powheg.job187", "tbar_t_powheg.job194", "tbar_tw_powheg.job187"]) org.mergeSamples(targetSpec = md({"name":"Z + jets", "color": r.kBlue}, mcOps), allWithPrefix = "zinv_mg_ht") - #org.mergeSamples(targetSpec = md({"name":"W + jets", "color": r.kOrange-3}, mcOps), allWithPrefix = "wj_lv_mg_ht_") + org.mergeSamples(targetSpec = md({"name":"W + jets", "color": r.kOrange-3}, mcOps), allWithPrefix = "wj_lv_mg_ht_") org.mergeSamples(targetSpec = md({"name":"VV", "color": r.kOrange+3}, mcOps), sources = ["ww_py.job188", "wz_py.job188", "zz_py.job188"]) org.mergeSamples(targetSpec = md({"name":"ZH", "color":r.kMagenta}, mcOps), sources = ["zinv_hbb_125_powheg.job342"]) org.mergeSamples(targetSpec = md({"name":"LM6", "color":r.kMagenta}, mcOps), allWithPrefix = "lm6") @@ -611,7 +599,7 @@ def makeStandardPlots(self, org) : showStatBox = True, rowColors = [r.kBlack, r.kViolet+4], #whiteList = ["lowestUnPrescaledTrigger"], - doLog = False, + doLog = True, #pegMinimum = 0.1, linYAfter = ("variableGreaterFilter", "xcak5JetAlphaTEtPat>=0.550 "), blackList = ["lumiHisto","xsHisto","nJobsHisto"], diff --git a/samples/__ht17__.py b/samples/__ht17__.py index 3420b4f..f50e736 100644 --- a/samples/__ht17__.py +++ b/samples/__ht17__.py @@ -4,11 +4,6 @@ ht17 = SampleHolder() a = ", alwaysUseLastAttempt = True" -#L1FJL2L3Residual -ht17.add("HT.Run2012A-PromptReco-v1.AOD.job229", '%s/zmeng//ICF/automated/2012_06_14_14_26_12/HT.Run2012A-PromptReco-v1.AOD/")'%pnfs, lumi = 1.0) #job 217, 226/226 comp. - -ht17.add("HTMHT.Run2012B-PromptReco-v1.AOD.job228", '%s/zmeng//ICF/automated/2012_06_14_11_22_04/HTMHT.Run2012B-PromptReco-v1.AOD/")'%pnfs, lumi = 1.0) -ht17.add("HTMHT.Run2012B-PromptReco-v1.AOD.job238", '%s/zmeng//ICF/automated/2012_06_22_14_25_16/HTMHT.Run2012B-PromptReco-v1.AOD/")'%pnfs, lumi = 1.0) #53X ht17.add("HTMHT.Run2012B-13Jul2012-v1.AOD.job358", '%s/yeshaq//ICF/automated/2012_09_19_01_03_03/HTMHT.Run2012B-13Jul2012-v1.AOD/")'%pnfs, lumi = 1.0) @@ -18,10 +13,6 @@ ht17.add("HT.Run2012A-13Jul2012-v1.AOD.job358", '%s/yeshaq//ICF/automated/2012_09_19_01_03_03/HT.Run2012A-13Jul2012-v1.AOD/")'%pnfs, lumi = 1.0) ht17.add("HT.Run2012A-recover-06Aug2012-v1.AOD.job359", '%s/yeshaq//ICF/automated/2012_09_19_01_25_46/HT.Run2012A-recover-06Aug2012-v1.AOD/")'%pnfs, lumi = 1.0) - - - - - - -#ht.add("HT.Run2011B-PromptReco-v1.AOD.job570", '%s/bm409//ICF/automated/2011_10_17_12_55_58/HT.Run2011B-PromptReco-v1.AOD")'%srm, lumi = 1.0) #job 570, 82/432 comp. +for era in ["B","C","D"] : + ht17.add("HTMHTParked.Run2012%s-22Jan2013-v1.job649"%era, '%s/yeshaq//ICF/automated/2013_04_19_17_30_49/HTMHTParked.Run2012%s-22Jan2013-v1.AOD/")'%(pnfs,era), lumi = 1.0) + From f66e862ee99683880eb5de9fc5f0512341c082da Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Tue, 2 Jul 2013 10:30:48 -0500 Subject: [PATCH 08/34] adding zinv samples, enabling METFilters, moving towards pep8 standard --- analyses/hadronicLook.py | 662 ++++++++++++++++++++------------------- 1 file changed, 333 insertions(+), 329 deletions(-) diff --git a/analyses/hadronicLook.py b/analyses/hadronicLook.py index d76fa8c..04f9660 100644 --- a/analyses/hadronicLook.py +++ b/analyses/hadronicLook.py @@ -17,7 +17,8 @@ def triggerTuple(l=[], keys=[]): out.append("%s_v%d" % (stem, version)) return tuple(out) -triggers_2012 = [{"HT": 250, "AlphaT": 0.55, "v": range(1, 10)}, +triggers_2012 = [{"HT": 200, "AlphaT": 0.57, "v": [5, 6, 8]}, + {"HT": 250, "AlphaT": 0.55, "v": range(1, 10)}, {"HT": 300, "AlphaT": 0.53, "v": range(1, 4)}, {"HT": 300, "AlphaT": 0.54, "v": range(1, 10)}, {"HT": 350, "AlphaT": 0.52, "v": range(1, 4)}, @@ -34,13 +35,13 @@ def parameters(self): {"muon": ("muon", "Pat"), "electron": ("electron", "Pat"), "photon": ("photon", "Pat"), - # + "jet": ("xcak5Jet", "Pat"), "jetId": "JetIDloose", "muonsInJets": False, "met": "metP4TypeIPF", "rechit": "Calo", - # + "jetComp": ("xcak5JetPF", "Pat"), "jetIdComp": "JetIDtight", "muonsInJetsComp": True, @@ -52,16 +53,16 @@ def parameters(self): #"le3j": (2, 3), #"ge4j": (4, None), }), - #"nBTagJets": self.vary({#"eq0b": (0, 0), + "nBTagJets": self.vary({"eq0b": (0, 0), #"eq1b": (1, 1), #"eq2b": (2, 2), #"eq3b": (3, 3), #"ge4b": (4, None), - #}), - "thresholds": self.vary({"200s": (200.0, None, 73.3, 36.7), + }), + "thresholds": self.vary({#"200s": (200.0, 275.0, 73.3, 36.7), #"275s": (275.0, 325.0, 73.3, 36.7), #"325s": (325.0, 375.0, 86.7, 43.3), - #"375": (375.0, None, 100.0, 50.0), + "375": (375.0, None, 100.0, 50.0), #"675": (675.0, None, 100.0, 50.0), #"875": (875.0, None, 100.0, 50.0), }), @@ -70,7 +71,7 @@ def parameters(self): "lowPtName": "lowPt", "highPtThreshold": 50.0, "highPtName": "highPt", - "signalScan": True, + "signalScan": False, #required to be sorted "triggerList": triggerTuple(l=triggers_2012, keys=("HT", "AlphaT"), @@ -223,27 +224,28 @@ def stepsTrigger(self, params): #steps.trigger.hltPrescaleHistogrammer(params["triggerList"]).onlyData(), steps.trigger.l1Filter("L1Tech_BPTX_plus_AND_minus.v0").onlyData(), steps.trigger.physicsDeclaredFilter().onlyData(), - ]+([] if params["thresholds"][1]!=325.0 else [steps.trigger.lowestUnPrescaledTriggerFilter().onlyData()]) #apply trigger in lowest HT bin + ]+([] if params["thresholds"][1] != 275.0 else [steps.trigger.lowestUnPrescaledTriggerFilter().onlyData()]) # apply trigger in lowest HT bin - def stepsGenValidation(self) : - return [#supy.steps.histos.histogrammer("genpthat",200,0,2000,title=";#hat{p_{T}} (GeV);events / bin").onlySim(), - supy.steps.histos.histogrammer("genPartonHT",200,0,1000,title=";parton H_{T} (GeV);events / bin").onlySim(), - ] + def stepsGenValidation(self): + return [ + #supy.steps.histos.histogrammer("genpthat",200,0,2000,title=";#hat{p_{T}} (GeV);events / bin").onlySim(), + supy.steps.histos.histogrammer("genPartonHT", 200, 0, 1000, title=";parton H_{T} (GeV);events / bin").onlySim(), + ] - def stepsEvent(self) : + def stepsEvent(self, params): return [steps.filters.monster(), steps.filters.hbheNoise().onlyData(), - #supy.steps.filters.value("beamHaloCSCTightHaloId").invert().onlyData(), - #supy.steps.filters.value("trackingFailureFilterFlag").onlyData(), - #supy.steps.filters.value("hcalLaserEventFilterFlag").onlyData(), - #supy.steps.filters.value("ecalDeadCellTPFilterFlag").onlyData(), + supy.steps.filters.value("beamHaloCSCTightHaloId", max=0).onlyData(), + supy.steps.filters.value("trackingFailureFilterFlag", min=1).onlyData(), + supy.steps.filters.value("hcalLaserEventFilterFlag", min=1).onlyData(), + supy.steps.filters.value("ecalDeadCellTPFilterFlag", min=1).onlyData(), #supy.steps.histos.histogrammer("logErrorTooManySeeds", 2, 0.0, 1.0, title = ";logErrorTooManySeeds;events / bin"), #supy.steps.histos.histogrammer("logErrorTooManyClusters", 2, 0.0, 1.0, title = ";logErrorTooManyClusters;events / bin"), - supy.steps.histos.multiplicity("vertexIndices", max = 30), - supy.steps.filters.multiplicity("vertexIndices", min = 1), + supy.steps.histos.multiplicity("vertexIndices", max=30), + supy.steps.filters.multiplicity("vertexIndices", min=1), ] - def stepsHtLeadingJets(self, params) : + def stepsHtLeadingJets(self, params): jet = params["objects"]["jet"] et = "Et" if params["etRatherThanPt"] else "Pt" @@ -256,107 +258,108 @@ def stepsHtLeadingJets(self, params) : out = [steps.jet.jetPtSelector(jet, params["thresholds"][2], 0), steps.jet.jetPtSelector(jet, params["thresholds"][2], 1), steps.jet.jetEtaSelector(jet, 2.5, 0), - supy.steps.filters.value("%sSum%s%s"%(jet[0], et, jet[1]), min = params["thresholds"][0]), + supy.steps.filters.value("%sSum%s%s" % (jet[0], et, jet[1]), min=params["thresholds"][0]), ] - if params["thresholds"][1]!=None : - out.append(supy.steps.filters.value("%sSum%s%s"%(jet[0], et, jet[1]), max = params["thresholds"][1])) + if params["thresholds"][1] is not None: + out.append(supy.steps.filters.value("%sSum%s%s" % (jet[0], et, jet[1]), max=params["thresholds"][1])) return out - def stepsBtagJets(self, params) : + def stepsBtagJets(self, params): _jet = params["objects"]["jet"] - out = [supy.steps.filters.multiplicity("%sIndicesBtagged2%s"%_jet, min = params["nBTagJets"][0], max = params["nBTagJets"][1]), - supy.steps.histos.multiplicity("%sIndicesBtagged2%s"%_jet), + out = [supy.steps.filters.multiplicity("%sIndicesBtagged2%s" % _jet, min=params["nBTagJets"][0], max=params["nBTagJets"][1]), + supy.steps.histos.multiplicity("%sIndicesBtagged2%s" % _jet), ] return out - - def stepsXclean(self, params) : + + def stepsXclean(self, params): return [ #steps.other.iterHistogrammer("ecalDeadTowerTrigPrimP4", 256, 0.0, 128.0, title=";E_{T} of ECAL TP in each dead region (GeV);TPs / bin", funcString="lambda x:x.Et()"), - supy.steps.filters.multiplicity("%sIndices%s"%params["objects"]["muon"], max = 0), - supy.steps.filters.multiplicity("%sIndices%s"%params["objects"]["electron"], max = 0), - supy.steps.filters.multiplicity("%sIndices%s"%params["objects"]["photon"], max = 0), - supy.steps.filters.multiplicity("%sIndicesOther%s"%params["objects"]["jet"], max = 0), - supy.steps.filters.multiplicity("%sIndices%s"%params["objects"]["jet"], min = params["nJetsMinMax"][0], max = params["nJetsMinMax"][1]), + supy.steps.filters.multiplicity("%sIndices%s" % params["objects"]["muon"], max=0), + supy.steps.filters.multiplicity("%sIndices%s" % params["objects"]["electron"], max=0), + supy.steps.filters.multiplicity("%sIndices%s" % params["objects"]["photon"], max=0), + supy.steps.filters.multiplicity("%sIndicesOther%s" % params["objects"]["jet"], max=0), + supy.steps.filters.multiplicity("%sIndices%s" % params["objects"]["jet"], min=params["nJetsMinMax"][0], max=params["nJetsMinMax"][1]), #steps.jet.uniquelyMatchedNonisoMuons(_jet), ] - def stepsPlotsOne(self, params) : + def stepsPlotsOne(self, params): _jet = params["objects"]["jet"] _met = params["objects"]["met"] _et = "Et" if params["etRatherThanPt"] else "Pt" return [ - supy.steps.histos.multiplicity("%sIndicesBtagged2%s"%_jet), - supy.steps.histos.histogrammer("%sSum%s%s"%(_jet[0], _et, _jet[1]), 50, 0, 2500, - title = ";H_{T} (GeV) from %s%s %ss;events / bin"%(_jet[0], _jet[1], _et)), - supy.steps.histos.histogrammer("%sSum%s%s"%(_jet[0], _et, _jet[1]), 60, 675, 1275, - title = ";H_{T} (GeV) from %s%s %ss;events / bin"%(_jet[0], _jet[1], _et)), - supy.steps.histos.histogrammer("%sSum%s%s"%(_jet[0], _et, _jet[1]), 100, 0, 1000, - title = ";H_{T} (GeV) from %s%s %ss;events / bin"%(_jet[0], _jet[1], _et)), - - supy.steps.histos.histogrammer("%sSumP4%s"%_jet, 50, 0, 500, title = ";MHT from %s%s (GeV);events / bin"%_jet, funcString = "lambda x:x.pt()"), + + supy.steps.histos.multiplicity("%sIndicesBtagged2%s" % _jet), + supy.steps.histos.histogrammer("%sSum%s%s" % (_jet[0], _et, _jet[1]), 50, 0, 2500, + title=";H_{T} (GeV) from %s%s %ss;events / bin" % (_jet[0], _jet[1], _et)), + supy.steps.histos.histogrammer("%sSum%s%s" % (_jet[0], _et, _jet[1]), 60, 675, 1275, + title=";H_{T} (GeV) from %s%s %ss;events / bin" % (_jet[0], _jet[1], _et)), + supy.steps.histos.histogrammer("%sSum%s%s" % (_jet[0], _et, _jet[1]), 100, 0, 1000, + title=";H_{T} (GeV) from %s%s %ss;events / bin" % (_jet[0], _jet[1], _et)), + + supy.steps.histos.histogrammer("%sSumP4%s" % _jet, 50, 0, 500, title=";MHT from %s%s (GeV);events / bin" % _jet, funcString="lambda x:x.pt()"), supy.steps.histos.histogrammer("vertexIndices", 20, -0.5, 19.5, title=";N vertices;events / bin", funcString="lambda x:len(x)"), - supy.steps.histos.histogrammer("vertexSumPt", 100, 0.0, 1.0e3, title = ";SumPt of 2nd vertex (GeV);events / bin", funcString = "lambda x:([0.0,0.0]+sorted(x))[-2]"), - + supy.steps.histos.histogrammer("vertexSumPt", 100, 0.0, 1.0e3, title=";SumPt of 2nd vertex (GeV);events / bin", funcString="lambda x:([0.0,0.0]+sorted(x))[-2]"), + #steps.trigger.lowestUnPrescaledTriggerHistogrammer(), #supy.steps.filters.label("singleJetPlots1"), - #steps.jet.singleJetHistogrammer(_jet), - supy.steps.filters.label("jetSumPlots1"), - steps.jet.cleanJetHtMhtHistogrammer(_jet,params["etRatherThanPt"]), - supy.steps.histos.histogrammer("%sDeltaPhiStar%s%s"%(_jet[0], params["lowPtName"], _jet[1]), 20, 0.0, r.TMath.Pi(), title = ";#Delta#phi*;events / bin", funcString = 'lambda x:x[0][0]'), - supy.steps.histos.histogrammer("%sDeltaPhiStar%s"%(_jet[0], _jet[1]), 20, 0.0, r.TMath.Pi(), title = ";#Delta#phi*;events / bin", funcString = 'lambda x:x[0][0]'), - supy.steps.histos.histogrammer("%sMaxEmEnergyFraction%s"%(_jet[0], _jet[1]), 20, 0.0, 1.0, title = ";MaxEmEnergyFraction;events / bin"), - supy.steps.histos.histogrammer(_met,100,0.0,500.0,title=";"+_met+" (GeV);events / bin", funcString = "lambda x: x.pt()"), + steps.jet.singleJetHistogrammer(_jet), + supy.steps.filters.label("jetSumPlots1"), + steps.jet.cleanJetHtMhtHistogrammer(_jet, params["etRatherThanPt"]), + supy.steps.histos.histogrammer("%sDeltaPhiStar%s%s" % (_jet[0], params["lowPtName"], _jet[1]), 20, 0.0, r.TMath.Pi(), title=";#Delta#phi*;events / bin", funcString='lambda x:x[0][0]'), + supy.steps.histos.histogrammer("%sDeltaPhiStar%s" % (_jet[0], _jet[1]), 20, 0.0, r.TMath.Pi(), title=";#Delta#phi*;events / bin", funcString='lambda x:x[0][0]'), + supy.steps.histos.histogrammer("%sMaxEmEnergyFraction%s" % (_jet[0], _jet[1]), 20, 0.0, 1.0, title=";MaxEmEnergyFraction;events / bin"), + supy.steps.histos.histogrammer(_met, 100, 0.0, 500.0, title=";"+_met+" (GeV);events / bin", funcString="lambda x: x.pt()"), supy.steps.filters.label("kinematicPlots1"), - steps.jet.alphaHistogrammer(cs = _jet, deltaPhiStarExtraName = params["lowPtName"], etRatherThanPt = params["etRatherThanPt"]), + steps.jet.alphaHistogrammer(cs=_jet, deltaPhiStarExtraName=params["lowPtName"], etRatherThanPt=params["etRatherThanPt"]), ] - def stepsQcdRejection(self, params) : + def stepsQcdRejection(self, params): _jet = params["objects"]["jet"] _met = params["objects"]["met"] _et = "Et" if params["etRatherThanPt"] else "Pt" return [ - supy.steps.histos.histogrammer("%sMht%sOver%s"%(_jet[0],params["highPtName"]+_jet[1],_met), 100, 0.0, 3.0, - title = ";MHT %s%s / %s;events / bin"%(_jet[0],params["highPtName"]+_jet[1],_met)), - supy.steps.filters.value("%sMht%sOver%s"%(_jet[0],params["highPtName"]+_jet[1],_met), max = 1.25), - steps.other.deadEcalFilter(jets = _jet, extraName = params["lowPtName"], dR = 0.3, dPhiStarCut = 0.5), - + supy.steps.histos.histogrammer("%sMht%sOver%s" % (_jet[0], params["highPtName"]+_jet[1], _met), 100, 0.0, 3.0, + title=";MHT %s%s / %s;events / bin" % (_jet[0], params["highPtName"]+_jet[1], _met)), + supy.steps.filters.value("%sMht%sOver%s" % (_jet[0], params["highPtName"]+_jet[1], _met), max=1.25), + steps.other.deadEcalFilter(jets=_jet, extraName=params["lowPtName"], dR=0.3, dPhiStarCut=0.5), + steps.jet.cleanJetHtMhtHistogrammer(_jet, params["etRatherThanPt"]), - steps.jet.alphaHistogrammer(cs = _jet, deltaPhiStarExtraName = params["lowPtName"], etRatherThanPt = params["etRatherThanPt"]), + steps.jet.alphaHistogrammer(cs=_jet, deltaPhiStarExtraName=params["lowPtName"], etRatherThanPt=params["etRatherThanPt"]), #steps.jet.alphaMetHistogrammer(cs = _jet, deltaPhiStarExtraName = params["lowPtName"], etRatherThanPt = _etRatherThanPt, metName = _met), - supy.steps.histos.histogrammer("%sRecHitSumPt"%params["objects"]["rechit"], 30, 0, 300, title = ";Sum of HBHE (sev.#geq10), EB,EE (sev.#geq2) RecHit p_{T} (GeV);events / bin"), - supy.steps.filters.value("%sRecHitSumPt"%params["objects"]["rechit"], max = 30.0), + supy.steps.histos.histogrammer("%sRecHitSumPt" % params["objects"]["rechit"], 30, 0, 300, title=";Sum of HBHE (sev.#geq10), EB,EE (sev.#geq2) RecHit p_{T} (GeV);events / bin"), + supy.steps.filters.value("%sRecHitSumPt" % params["objects"]["rechit"], max=30.0), - supy.steps.filters.value("%sAlphaT%s%s"%(_jet[0], _et, _jet[1]), min = 0.55), + supy.steps.filters.value("%sAlphaT%s%s" % (_jet[0], _et, _jet[1]), min=0.55), - supy.steps.histos.histogrammer("%sMaxEmEnergyFraction%s"%(_jet[0], _jet[1]), 20, 0.0, 1.0, title = ";MaxEmEnergyFraction;events / bin"), + supy.steps.histos.histogrammer("%sMaxEmEnergyFraction%s" % (_jet[0], _jet[1]), 20, 0.0, 1.0, title=";MaxEmEnergyFraction;events / bin"), #supy.steps.filters.value("%sMaxEmEnergyFraction%s"%(_jet[0],_jet[1]), max = .1), ] - def stepsPlotsTwo(self, params) : + def stepsPlotsTwo(self, params): _jet = params["objects"]["jet"] _met = params["objects"]["met"] return [ supy.steps.histos.histogrammer("vertexIndices", 20, -0.5, 19.5, title=";N vertices;events / bin", funcString="lambda x:len(x)"), - supy.steps.histos.histogrammer("%sIndices%s"%_jet, 20, -0.5, 19.5, title=";number of %s%s passing ID#semicolon p_{T}#semicolon #eta cuts;events / bin"%_jet, funcString="lambda x:len(x)"), - steps.jet.cleanJetHtMhtHistogrammer(_jet,params["etRatherThanPt"]), - supy.steps.histos.histogrammer("%sDeltaPhiStar%s%s"%(_jet[0], params["lowPtName"], _jet[1]), 20, 0.0, r.TMath.Pi(), title = ";#Delta#phi*;events / bin", funcString = 'lambda x:x[0][0]'), - supy.steps.histos.histogrammer("%sDeltaPhiStar%s"%(_jet[0], _jet[1]), 20, 0.0, r.TMath.Pi(), title = ";#Delta#phi*;events / bin", funcString = 'lambda x:x[0][0]'), - supy.steps.histos.histogrammer("%sMht%sOver%s"%(_jet[0],params["highPtName"]+_jet[1],_met), 100, 0.0, 3.0, - title = ";MHT %s%s / %s;events / bin"%(_jet[0],params["highPtName"]+_jet[1],_met)), + supy.steps.histos.histogrammer("%sIndices%s" % _jet, 20, -0.5, 19.5, title=";number of %s%s passing ID#semicolon p_{T}#semicolon #eta cuts;events / bin" % _jet, funcString="lambda x:len(x)"), + steps.jet.cleanJetHtMhtHistogrammer(_jet, params["etRatherThanPt"]), + supy.steps.histos.histogrammer("%sDeltaPhiStar%s%s" % (_jet[0], params["lowPtName"], _jet[1]), 20, 0.0, r.TMath.Pi(), title=";#Delta#phi*;events / bin", funcString='lambda x:x[0][0]'), + supy.steps.histos.histogrammer("%sDeltaPhiStar%s" % (_jet[0], _jet[1]), 20, 0.0, r.TMath.Pi(), title=";#Delta#phi*;events / bin", funcString='lambda x:x[0][0]'), + supy.steps.histos.histogrammer("%sMht%sOver%s" % (_jet[0], params["highPtName"]+_jet[1], _met), 100, 0.0, 3.0, + title=";MHT %s%s / %s;events / bin" % (_jet[0], params["highPtName"]+_jet[1], _met)), ] - def stepsHtBins(self, params) : + def stepsHtBins(self, params): _jet = params["objects"]["jet"] - return [supy.steps.filters.value("%sSumEt%s"%_jet, min = bin) for bin in [475, 575, 675, 775, 875]] + return [supy.steps.filters.value("%sSumEt%s" % _jet, min=bin) for bin in [475, 575, 675, 775, 875]] - def stepsOptional(self, params) : + def stepsOptional(self, params): return [ - supy.steps.other.skimmer(), - #steps.other.duplicateEventCheck(), + #supy.steps.other.skimmer(), + steps.other.duplicateEventCheck(), #steps.other.pickEventSpecMaker(), #steps.other.cutBitHistogrammer(self.togglePfJet(_jet), self.togglePfMet(_met)), #steps.Print.jetPrinter(_jet), @@ -365,133 +368,136 @@ def stepsOptional(self, params) : #steps.Gen.genParticlePrinter(minPt = 10.0, minStatus = 3), ] - def stepsDisplayer(self, params) : + def stepsDisplayer(self, params): jet = params["objects"]["jet"] return [ - steps.displayer.displayer(jets = jet, - muons = params["objects"]["muon"], - met = params["objects"]["met"], - electrons = params["objects"]["electron"], - photons = params["objects"]["photon"], - recHits = params["objects"]["rechit"], recHitPtThreshold = 1.0,#GeV - scale = 400.0,#GeV - etRatherThanPt = params["etRatherThanPt"], - deltaPhiStarExtraName = params["lowPtName"], - deltaPhiStarCut = 0.5, - deltaPhiStarDR = 0.3, - j2Factor = params["thresholds"][2]/params["thresholds"][0], - mhtOverMetName = "%sMht%sOver%s"%(jet[0],params["highPtName"]+jet[1],params["objects"]["met"]), - metOtherAlgo = params["objects"]["metComp"], - jetsOtherAlgo = params["objects"]["jetComp"], - #doGenJets = True, + steps.displayer.displayer(jets=jet, + muons=params["objects"]["muon"], + met=params["objects"]["met"], + electrons=params["objects"]["electron"], + photons=params["objects"]["photon"], + recHits=params["objects"]["rechit"], recHitPtThreshold=1.0, # GeV + scale=400.0, # GeV + etRatherThanPt=params["etRatherThanPt"], + deltaPhiStarExtraName=params["lowPtName"], + deltaPhiStarCut=0.5, + deltaPhiStarDR=0.3, + j2Factor=params["thresholds"][2]/params["thresholds"][0], + mhtOverMetName="%sMht%sOver%s" % (jet[0], params["highPtName"]+jet[1], params["objects"]["met"]), + metOtherAlgo=params["objects"]["metComp"], + jetsOtherAlgo=params["objects"]["jetComp"], + #doGenJets=True, ), ] - def stepsMbb(self, params) : + def stepsMbb(self, params): _jet = params["objects"]["jet"] return [ - #supy.steps.filters.multiplicity("genIndicesStatus3b", min = 4, max = 4), + #supy.steps.filters.multiplicity("genIndicesStatus3b", min=4, max=4), #steps.printer.eventPrinter(), #steps.printer.jetPrinter(_jet), #steps.gen.particlePrinter(), - supy.steps.filters.multiplicity("%sIndicesBtagged2%s"%_jet, min = 2, max = 2), - supy.steps.histos.multiplicity("%sIndicesBtagged2%s"%_jet), - supy.steps.histos.eta("%sCorrectedP4%s"%_jet, 24, -3.0, 3.0, indices = "%sIndicesBtagged2%s"%_jet, index = 0, xtitle = "b jet 0"), - supy.steps.histos.eta("%sCorrectedP4%s"%_jet, 24, -3.0, 3.0, indices = "%sIndicesBtagged2%s"%_jet, index = 1, xtitle = "b jet 1"), - supy.steps.histos.mass("%sSumP4Btagged2%s"%_jet, 24, 0.0, 1200.0, xtitle = "sum P4 {b jets}"), - supy.steps.histos.pt("%sSumP4Btagged2%s"%_jet, 24, 0.0, 1200.0, xtitle = "sum P4 {b jets}"), - supy.steps.filters.mass("%sSumP4Btagged2%s"%_jet, min = 450.0), - supy.steps.histos.pt("%sSumP4Btagged2%s"%_jet, 24, 0.0, 1200.0, xtitle = "sum P4 {b jets}"), - #steps.jet.mbbHistogrammer(_jet, drMatch = 0.2, bZDaughters = "genIndicesStatus3bZDaughters"), + supy.steps.filters.multiplicity("%sIndicesBtagged2%s" % _jet, min=2, max=2), + supy.steps.histos.multiplicity("%sIndicesBtagged2%s" % _jet), + supy.steps.histos.eta("%sCorrectedP4%s" % _jet, 24, -3.0, 3.0, indices="%sIndicesBtagged2%s" % _jet, index=0, xtitle="b jet 0"), + supy.steps.histos.eta("%sCorrectedP4%s" % _jet, 24, -3.0, 3.0, indices="%sIndicesBtagged2%s" % _jet, index=1, xtitle="b jet 1"), + supy.steps.histos.mass("%sSumP4Btagged2%s" % _jet, 24, 0.0, 1200.0, xtitle="sum P4 {b jets}"), + supy.steps.histos.pt("%sSumP4Btagged2%s" % _jet, 24, 0.0, 1200.0, xtitle="sum P4 {b jets}"), + supy.steps.filters.mass("%sSumP4Btagged2%s" % _jet, min=450.0), + supy.steps.histos.pt("%sSumP4Btagged2%s" % _jet, 24, 0.0, 1200.0, xtitle="sum P4 {b jets}"), + #steps.jet.mbbHistogrammer(_jet, drMatch=0.2, bZDaughters="genIndicesStatus3bZDaughters"), ] - def listOfSteps(self, params) : - _jet = params["objects"]["jet"] - _met = params["objects"]["met"] + def listOfSteps(self, params): + _jet = params["objects"]["jet"] + _met = params["objects"]["met"] _etRatherThanPt = params["etRatherThanPt"] _et = "Et" if _etRatherThanPt else "Pt" return ([supy.steps.printer.progressPrinter()] + - #self.stepsEventCount(params, label = "scanBefore") + + #self.stepsEventCount(params, label="scanBefore") + self.stepsGenValidation() + - self.stepsEvent() + - #self.stepsTrigger(params) + + self.stepsEvent(params) + + self.stepsTrigger(params) + self.stepsHtLeadingJets(params) + - #self.stepsXclean(params) + - #self.stepsBtagJets(params) + - #self.stepsPlotsOne(params) + + self.stepsXclean(params) + + self.stepsOptional(params) + + self.stepsBtagJets(params) + + self.stepsPlotsOne(params) + self.stepsQcdRejection(params) + self.stepsPlotsTwo(params) + #self.stepsMbb(params) + #self.stepsDisplayer(params) + - #self.stepsOptional(params) + + #self.stepsHtBins(params) + - #self.stepsEventCount(params, label = "scanAfter") + + #self.stepsEventCount(params, label="scanAfter") + []) - def listOfSampleDictionaries(self) : + def listOfSampleDictionaries(self): sh = supy.samples.SampleHolder() - sh.add("275_ge2b", '["/uscms/home/elaird/08_mbb/02_skim/2012_5fb_275_ge2b.root"]', lumi = 5.0e3) - sh.add("375_ge2b", '["/uscms/home/yeshaq/nobackup/supy-output/hadronicLook/375_calo_ge2/HadronicRegion.root"]', lumi = 5.0e3) + sh.add("275_ge2b", '["/uscms/home/elaird/08_mbb/02_skim/2012_5fb_275_ge2b.root"]', lumi=5.0e3) + sh.add("375_ge2b", '["/uscms/home/yeshaq/nobackup/supy-output/hadronicLook/375_calo_ge2/HadronicRegion.root"]', lumi=5.0e3) return [samples.ht17, samples.top17, samples.ewk17, samples.qcd17, sh, samples.susy17] - - def listOfSamples(self,params) : - from supy.samples import specify - def data_53X() : - jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt") + def listOfSamples(self, params): + from supy.samples import specify + def data_53X(): + jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptRecon_Collisions12_JSON.txt") out = [] - out += specify(names = "HTMHTParked.Run2012B-22Jan2013-v1.job649", weights = jw2012, nFilesMax = 1, nEventsMax = 20000) - out += specify(names = "HTMHTParked.Run2012C-22Jan2013-v1.job649", weights = jw2012, nFilesMax = 1, nEventsMax = 20000) - out += specify(names = "HTMHTParked.Run2012D-22Jan2013-v1.job649", weights = jw2012, nFilesMax = 1, nEventsMax = 20000) + out += specify(names="HTMHTParked.Run2012B-22Jan2013-v1.job649", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="HTMHTParked.Run2012C-22Jan2013-v1.job649", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="HTMHTParked.Run2012D-22Jan2013-v1.job649", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) + #out += specify(names="HTMHTParked_ICF_sync_test", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) return out - def data_52X_2b_skim() : - return specify(names = "275_ge2b") - - def qcd_py6(eL) : - low = map(lambda x:x[0],samples.__qcd17__.binsXs)[:-1] + def qcd_py6(eL): + low = map(lambda x: x[0], samples.__qcd17__.binsXs)[:-1] out = [] - for pt in low[low.index(80):] : - out += specify("qcd_py6_pt_%d"%pt, effectiveLumi = eL) + for pt in low[low.index(80):]: + out += specify("qcd_py6_pt_%d" % pt, effectiveLumi=eL) return out - def qcd_b_py6(eL) : + def qcd_b_py6(eL): out = [] - out += specify("qcd_b_py6_pt_50", effectiveLumi = eL) - out += specify("qcd_b_py6_pt_150", effectiveLumi = eL) + out += specify("qcd_b_py6_pt_50", effectiveLumi=eL) + out += specify("qcd_b_py6_pt_150", effectiveLumi=eL) return out - def g_jets_mg(eL) : - gM = [40,100,200] - return specify( effectiveLumi = eL, color = r.kGreen, - names = [("g_jets_mg_ht_%d_%d")[:None if t[1] else -2] for t in zip(gM,gM[1:]+["inf"])] ) + def g_jets_mg(eL): + gM = [40, 100, 200] + return specify(effectiveLumi=eL, color=r.kGreen, + names=[("g_jets_mg_ht_%d_%d")[:None if t[1] else -2] for t in zip(gM, gM[1:]+["inf"])]) - def w_binned() : + def w_binned(): out = [] - out += specify(names = "wj_lv_mg_ht_10_150", color = r.kBlue, nFilesMax = 1, nEventsMax = 20000) - out += specify(names = "wj_lv_mg_ht_150_200.job663", color = r.kGreen, nFilesMax = 1, nEventsMax = 20000) - out += specify(names = "wj_lv_mg_ht_200_250.job672", color = r.kCyan, nFilesMax = 1, nEventsMax = 20000) - out += specify(names = "wj_lv_mg_ht_250_300.job498", color = r.kOrange, nFilesMax = 1, nEventsMax = 20000) - out += specify(names = "wj_lv_mg_ht_300_400.job498", color = r.kViolet, nFilesMax = 1, nEventsMax = 20000) - out += specify(names = "wj_lv_mg_ht_400_inf.job498", color = r.kAzure, nFilesMax = 1, nEventsMax = 20000) + out += specify(names="wj_lv_mg_ht_10_150", color=r.kBlue) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="wj_lv_mg_ht_150_200.job663", color=r.kGreen) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="wj_lv_mg_ht_200_250.job672", color=r.kCyan) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="wj_lv_mg_ht_250_300.job498", color=r.kOrange) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="wj_lv_mg_ht_300_400.job498", color=r.kViolet) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="wj_lv_mg_ht_400_inf.job498", color=r.kAzure) # , nFilesMax=1, nEventsMax=20000) return out - def z_binned() : + def z_binned(): out = [] - out += specify("zinv_mg_ht_50_100.job214") - out += specify("zinv_mg_ht_100_200.job234") - out += specify("zinv_mg_ht_200_400.job233") - out += specify("zinv_mg_ht_400_inf.job213") + out += specify("zinv_mg_ht_50_100.job407", color=r.kBlue) # , nFilesMax=1, nEventsMax=20000) + out += specify("zinv_mg_ht_100_200.job365", color=r.kGreen) # , nFilesMax=1, nEventsMax=20000) + out += specify("zinv_mg_ht_200_400.job365", color=r.kOrange) # , nFilesMax=1, nEventsMax=20000) + out += specify("zinv_mg_ht_400_inf.job365", color=r.kViolet) # , nFilesMax=1, nEventsMax=20000) + out += specify("zinv_mg_ht_50_100_ext.job500", color=r.kBlue) # , nFilesMax=1, nEventsMax=20000) + out += specify("zinv_mg_ht_100_200_ext.job680", color=r.kGreen) # , nFilesMax=1, nEventsMax=20000) + out += specify("zinv_mg_ht_200_400_ext.job500", color=r.kOrange) # , nFilesMax=1, nEventsMax=20000) + out += specify("zinv_mg_ht_400_inf_ext.job500", color=r.kViolet) # , nFilesMax=1, nEventsMax=20000) return out - def w_inclusive() : + def w_inclusive(): out = [] - for part in range(1,6): - out += specify(names = "wj_lv_mg_ht_incl_v2.job673_part%i"%part, nFilesMax = 1, nEventsMax = 1000) + for part in range(1, 6): + out += specify(names="wj_lv_mg_ht_incl_v2.job673_part%i" % part, nFilesMax=1, nEventsMax=1000) return out - def vv() : + + def vv(): out = [] out += specify("ww_py.job188") out += specify("wz_py.job188") @@ -499,39 +505,38 @@ def vv() : out += specify("zinv_hbb_125_powheg.job342") return out - def top() : + def top(): out = [] out += specify(names="ttbar_powheg_v1.job410") - #out += specify(names = "tt_8_mg.job188") - #out += specify(names = "tt_8_mg.job315") - out += specify(names = "ttz_8_mg.job269", nFilesMax = 1) - - out += specify("t_s_powheg.job200" ) - #out += specify("t_t_powheg.job187" ) #low MC stats - out += specify("t_tw_powheg.job187" ) - out += specify("tbar_t_powheg.job194" ) + #out += specify(names="tt_8_mg.job188") + #out += specify(names="tt_8_mg.job315") + out += specify(names="ttz_8_mg.job269", nFilesMax=1) + + out += specify("t_s_powheg.job200") + #oAut += specify("t_t_powheg.job187" ) #low MC stats + out += specify("t_tw_powheg.job187") + out += specify("tbar_t_powheg.job194") out += specify("tbar_tw_powheg.job187") - + return out - def susy(eL) : - return specify(names = "lm6", effectiveLumi = eL, color = r.kRed) + def susy(eL): + return specify(names="lm6", effectiveLumi=eL, color=r.kRed) - def sms() : + def sms(): out = [] - #out += specify(names = "t2bb.job418")#, nFilesMax = 1, nEventsMax = 500) - out += specify(names = "t1tttt.job442")#, nFilesMax = 1, nEventsMax = 500) - #out += specify(names = "t1bbbb.job443", nFilesMax = 1, nEventsMax = 500) - #out += specify(names = "t1.job444", nFilesMax = 1, nEventsMax = 500) - #out += specify(names = "t2tt.job445")#, nFilesMax = 1, nEventsMax = 500) - #out += specify(names = "t2.job446", nFilesMax = 1, nEventsMax = 500) + #out += specify(names="t2bb.job418")#, nFilesMax=1, nEventsMax=500) + out += specify(names="t1tttt.job442") # , nFilesMax=1, nEventsMax=500) + #out += specify(names="t1bbbb.job443", nFilesMax=1, nEventsMax=500) + #out += specify(names="t1.job444", nFilesMax=1, nEventsMax=500) + #out += specify(names="t2tt.job445")#, nFilesMax=1, nEventsMax=500) + #out += specify(names="t2.job446", nFilesMax=1, nEventsMax=500) return out return ( data_53X() + - #data_52X_2b_skim() + - #w_binned() + - #z_binned() + +# w_binned() + +# z_binned() + #top() + #vv() + #qcd_py6(30.0e3) + @@ -541,212 +546,211 @@ def sms() : [] ) - def mergeSamples(self, org) : - def md(x, y) : + def mergeSamples(self, org): + def md(x, y): x.update(y) return x - - org.mergeSamples(targetSpec = {"name":"2012 Data", "color":r.kBlack, "markerStyle":20}, allWithPrefix = "HT") - mcOps = {"markerStyle":1, "lineWidth":3, "goptions":"hist"} + org.mergeSamples(targetSpec={"name": "2012 Data", "color": r.kBlack, "markerStyle": 20}, allWithPrefix="HT") + + mcOps = {"markerStyle": 1, "lineWidth": 3, "goptions": "hist"} qcdSources = [] - org.mergeSamples(targetSpec = md({"name":"QCD Multijet", "color":r.kGreen+3}, mcOps), allWithPrefix = "qcd_py6") + org.mergeSamples(targetSpec=md({"name": "QCD Multijet", "color": r.kGreen+3}, mcOps), allWithPrefix="qcd_py6") qcdSources = ["QCD Multijet"] - #org.mergeSamples(targetSpec = md({"name":"QCD Multijet (b-en.)", "color":r.kGreen+3}, mcOps), allWithPrefix = "qcd_b_py6") - #qcdSources = ["QCD Multijet (b-en.)"] + #org.mergeSamples(targetSpec=md({"name":"QCD Multijet (b-en.)", "color":r.kGreen+3}, mcOps), allWithPrefix="qcd_b_py6") + #qcdSources=["QCD Multijet (b-en.)"] + + #org.mergeSamples(targetSpec=md({"name":"ttz", "color": r.kYellow}, mcOps), allWithPrefix="ttz") + #org.mergeSamples(targetSpec=md({"name":"tt", "color": r.kRed+1}, mcOps), allWithPrefix="tt_") + #org.mergeSamples(targetSpec=md({"name":"t", "color": r.kGreen}, mcOps), + # sources=["t_s_powheg.job200", "t_t_powheg.job187", "t_tw_powheg.job187", "tbar_t_powheg.job194", "tbar_tw_powheg.job187"]) - #org.mergeSamples(targetSpec = md({"name":"ttz", "color": r.kYellow}, mcOps), allWithPrefix = "ttz") - #org.mergeSamples(targetSpec = md({"name":"tt", "color": r.kRed+1}, mcOps), allWithPrefix = "tt_") - #org.mergeSamples(targetSpec = md({"name":"t", "color": r.kGreen}, mcOps), - # sources = ["t_s_powheg.job200", "t_t_powheg.job187", "t_tw_powheg.job187", "tbar_t_powheg.job194", "tbar_tw_powheg.job187"]) - - org.mergeSamples(targetSpec = md({"name":"tt/t/ttz", "color":r.kRed+1}, mcOps), sources = [ + org.mergeSamples(targetSpec=md({"name": "tt/t/ttz", "color": r.kRed+1}, mcOps), sources=[ "tt_8_mg.job315", "ttz_8_mg.job269", "t_s_powheg.job200", "t_t_powheg.job187", "t_tw_powheg.job187", "tbar_t_powheg.job194", "tbar_tw_powheg.job187"]) - org.mergeSamples(targetSpec = md({"name":"Z + jets", "color": r.kBlue}, mcOps), allWithPrefix = "zinv_mg_ht") - org.mergeSamples(targetSpec = md({"name":"W + jets", "color": r.kOrange-3}, mcOps), allWithPrefix = "wj_lv_mg_ht_") - org.mergeSamples(targetSpec = md({"name":"VV", "color": r.kOrange+3}, mcOps), sources = ["ww_py.job188", "wz_py.job188", "zz_py.job188"]) - org.mergeSamples(targetSpec = md({"name":"ZH", "color":r.kMagenta}, mcOps), sources = ["zinv_hbb_125_powheg.job342"]) - org.mergeSamples(targetSpec = md({"name":"LM6", "color":r.kMagenta}, mcOps), allWithPrefix = "lm6") + org.mergeSamples(targetSpec=md({"name": "Z + jets", "color": r.kBlue}, mcOps), allWithPrefix="zinv_mg_ht") + org.mergeSamples(targetSpec=md({"name": "W + jets", "color": r.kOrange-3}, mcOps), allWithPrefix="wj_lv_mg_ht_") + org.mergeSamples(targetSpec=md({"name": "VV", "color": r.kOrange+3}, mcOps), sources=["ww_py.job188", "wz_py.job188", "zz_py.job188"]) + org.mergeSamples(targetSpec=md({"name": "ZH", "color": r.kMagenta}, mcOps), sources=["zinv_hbb_125_powheg.job342"]) + org.mergeSamples(targetSpec=md({"name": "LM6", "color": r.kMagenta}, mcOps), allWithPrefix="lm6") ewkSources = ["tt/t/ttz", "Z + jets", "W + jets", "VV"] - org.mergeSamples(targetSpec = md({"name":"Standard Model ", "color":r.kAzure+6}, mcOps), sources = ewkSources + qcdSources, keepSources = True) - - def conclude(self, conf) : - org = self.organizer(conf) - ##for skimming only - #utils.printSkimResults(org) + org.mergeSamples(targetSpec=md({"name": "Standard Model ", "color": r.kAzure+6}, mcOps), sources=ewkSources + qcdSources, keepSources=True) + def conclude(self, conf): + org = self.organizer(conf, verbose=True) self.mergeSamples(org) + ##for skimming only + #utils.printSkimResults(org) org.scale() if not self.parameters()["signalScan"] else org.scale(100.0) - self.makeStandardPlots(org) #self.makeIndividualPlots(org) - self.makeEfficiencyPlots(org) + #self.makeEfficiencyPlots(org) - def makeStandardPlots(self, org) : + def makeStandardPlots(self, org): #plot pl = supy.plotter(org, - pdfFileName = self.pdfFileName(org.tag), - samplesForRatios = ("2012 Data","Standard Model "), - sampleLabelsForRatios = ("data","s.m."), - #samplesForRatios = ("2012 Data","tt"), - #sampleLabelsForRatios = ("data","tt"), - printRatios = True, - showStatBox = True, - rowColors = [r.kBlack, r.kViolet+4], - #whiteList = ["lowestUnPrescaledTrigger"], - doLog = True, - #pegMinimum = 0.1, - linYAfter = ("variableGreaterFilter", "xcak5JetAlphaTEtPat>=0.550 "), - blackList = ["lumiHisto","xsHisto","nJobsHisto"], + pdfFileName=self.pdfFileName(org.tag), + samplesForRatios=("2012 Data", "Standard Model "), + sampleLabelsForRatios=("data", "s.m."), + #samplesForRatios=("2012 Data","tt"), + #sampleLabelsForRatios=("data","tt"), + printRatios=True, + showStatBox=True, + rowColors=[r.kBlack, r.kViolet+4], + #whiteList=["lowestUnPrescaledTrigger"], + doLog=True, + #pegMinimum=0.1, + linYAfter=("variableGreaterFilter", "xcak5JetAlphaTEtPat>=0.550 "), + blackList=["lumiHisto", "xsHisto", "nJobsHisto"], ) pl.plotAll() - def makeEfficiencyPlots(self, org) : - for dct in org.samples : - self.makeEffPlots(org, sampleName = dct["name"]) + def makeEfficiencyPlots(self, org): + for dct in org.samples: + self.makeEffPlots(org, sampleName=dct["name"]) - def makeIndividualPlots(self, org) : + def makeIndividualPlots(self, org): ht = "H_{T}^{#color[0]{T}}" - if "275" in org.tag : - htLabel = "275 < %s < 325"%ht - if "325" in org.tag : - htLabel = "325 < %s < 375"%ht - if "375" in org.tag : - htLabel = "375 < %s"%ht + if "275" in org.tag: + htLabel = "275 < %s < 325" % ht + if "325" in org.tag: + htLabel = "325 < %s < 375" % ht + if "375" in org.tag: + htLabel = "375 < %s" % ht #plot all pl = supy.plotter(org, - pdfFileName = self.pdfFileName(org.tag), - showStatBox = False, - doLog = False, - #pegMinimum = 0.1, - anMode = False, + pdfFileName=self.pdfFileName(org.tag), + showStatBox=False, + doLog=False, + #pegMinimum=0.1, + anMode=False, ) - pl.individualPlots(plotSpecs = [{"plotName":"xcak5JetAlphaTRoughPat", - "stepName" :"alphaHistogrammer", - "stepDesc" :"xcak5JetPat", - "newTitle":";#alpha_{T};events / bin", - "legendCoords": (0.55, 0.60, 0.85, 0.90), - "stampCoords": (0.75, 0.55) - }, - {"plotName":"jetMultiplicity", - "stepName":"singleJetHistogrammer", - "stepDesc":"xcak5JetPat through index 2", - "newTitle":";N_{jets};events / bin", - "legendCoords": (0.7, 0.7, 0.92, 0.92), - "stampCoords": (0.5, 0.28), - }, - {"plotName":"xcak5JetHtPat", - "stepName":"cleanJetHtMhtHistogrammer", - "stepDesc":"xcak5JetPat", - "newTitle":";H_{T} (GeV);events / bin", - "legendCoords": (0.6, 0.60, 0.92, 0.92), - "stampCoords": (0.45, 0.88) - }, - ##after alphaT - {"plotName":"xcak5JetDeltaPhiStarPat", - "stepName":"histogrammer", - "stepDesc" :"(lambda x:x[0][0])(xcak5JetDeltaPhiStarPat)", - "index" : -1, - "newTitle":";#Delta#phi*;events / bin", - "legendCoords": (0.6, 0.6, 0.92, 0.92), - "stampCoords": (0.33, 0.88), - }, - {"plotName":"xcak5JetHtPlusMhtRoughPat", - "stepName":"cleanJetHtMhtHistogrammer", - "stepDesc":"xcak5JetPat", - "index":-1, - "newTitle":";M_{eff} (GeV);events / bin", - "legendCoords": (0.7, 0.7, 0.92, 0.92), - "stampCoords": (0.75, 0.4), - }, - {"plotName":"xcak5JetIndicesPat", - "stepName":"histogrammer", - "stepDesc" :"(lambda x:len(x))(xcak5JetIndicesPat)", - "index" : -1, - "newTitle":";N_{jets};events / bin", - "legendCoords": (0.6, 0.6, 0.92, 0.92), - "stampCoords": (0.6, 0.38), - }, - {"plotName":"xcak5JetMbbListPat_2b", - "stepName":"mbbHistogrammer", - "stepDesc" :"mbbHistogrammer", - "index" : -1, - "stamp": False, - "newTitle":"%s;m_{bb} (GeV) [2 b-jets];events / bin / 5.0 fb^{-1}"%htLabel, - "legendCoords": (0.55, 0.6, 0.8, 0.75), - "reBinFactor": 5, - }, - ][-1:], - newSampleNames = {"tt": "Madgraph t#bar{t}", - "2012 Data": "Data"}, - #newSampleNames = {"qcd_mg_nVtx": "Madgraph QCD", + pl.individualPlots(plotSpecs=[{"plotName": "xcak5JetAlphaTRoughPat", + "stepName": "alphaHistogrammer", + "stepDesc": "xcak5JetPat", + "newTitle": ";#alpha_{T};events / bin", + "legendCoords": (0.55, 0.60, 0.85, 0.90), + "stampCoords": (0.75, 0.55) + }, + {"plotName": "jetMultiplicity", + "stepName": "singleJetHistogrammer", + "stepDesc": "xcak5JetPat through index 2", + "newTitle": ";N_{jets};events / bin", + "legendCoords": (0.7, 0.7, 0.92, 0.92), + "stampCoords": (0.5, 0.28), + }, + {"plotName": "xcak5JetHtPat", + "stepName": "cleanJetHtMhtHistogrammer", + "stepDesc": "xcak5JetPat", + "newTitle": ";H_{T} (GeV);events / bin", + "legendCoords": (0.6, 0.60, 0.92, 0.92), + "stampCoords": (0.45, 0.88) + }, + ##after alphaT + {"plotName": "xcak5JetDeltaPhiStarPat", + "stepName": "histogrammer", + "stepDesc": "(lambda x:x[0][0])(xcak5JetDeltaPhiStarPat)", + "index": -1, + "newTitle": ";#Delta#phi*;events / bin", + "legendCoords": (0.6, 0.6, 0.92, 0.92), + "stampCoords": (0.33, 0.88), + }, + {"plotName": "xcak5JetHtPlusMhtRoughPat", + "stepName": "cleanJetHtMhtHistogrammer", + "stepDesc": "xcak5JetPat", + "index": -1, + "newTitle": ";M_{eff} (GeV);events / bin", + "legendCoords": (0.7, 0.7, 0.92, 0.92), + "stampCoords": (0.75, 0.4), + }, + {"plotName": "xcak5JetIndicesPat", + "stepName": "histogrammer", + "stepDesc": "(lambda x:len(x))(xcak5JetIndicesPat)", + "index": -1, + "newTitle": ";N_{jets};events / bin", + "legendCoords": (0.6, 0.6, 0.92, 0.92), + "stampCoords": (0.6, 0.38), + }, + {"plotName": "xcak5JetMbbListPat_2b", + "stepName": "mbbHistogrammer", + "stepDesc": "mbbHistogrammer", + "index": -1, + "stamp": False, + "newTitle": "%s;m_{bb} (GeV) [2 b-jets];events / bin / 5.0 fb^{-1}" % htLabel, + "legendCoords": (0.55, 0.6, 0.8, 0.75), + "reBinFactor": 5, + }, + ][-1:], + newSampleNames={"tt": "Madgraph t#bar{t}", + "2012 Data": "Data"}, + #newSampleNames={"qcd_mg_nVtx": "Madgraph QCD", # "g_jets_mg_nVtx": "Madgraph #gamma + jets", # "2011 Data": "Data", # "standard_model_nVtx": "Standard Model", # }, - preliminary = True, - tdrStyle = False, + preliminary=True, + tdrStyle=False, ) + def makeEffPlots(self, org, sampleName): - def makeEffPlots(self, org, sampleName) : - - def sampleIndex(org, name) : - for iSample,sample in enumerate(org.samples) : - if sample["name"]==name : + def sampleIndex(org, name): + for iSample, sample in enumerate(org.samples): + if sample["name"] == name: return iSample - assert False, "could not find sample %s"%name - - def numerAndDenom(org, var) : + assert False, "could not find sample %s" % name + + def numerAndDenom(org, var): d = {} - - for selection in org.steps : - if "scanBefore" in selection.title : label = "before" - elif "scanAfter" in selection.title : label = "after" - if selection.name!= "scanHistogrammer" : continue + for selection in org.steps: + if "scanBefore" in selection.title: + label = "before" + elif "scanAfter" in selection.title: + label = "after" + if selection.name != "scanHistogrammer": + continue d[label] = selection[var][sampleIndex(org, sampleName)].Clone(label) return d - + keep = [] - file = r.TFile("%s_%s.root"%(sampleName, org.tag), "RECREATE") + file = r.TFile("%s_%s.root" % (sampleName, org.tag), "RECREATE") canvas = r.TCanvas() canvas.SetRightMargin(0.2) canvas.SetTickx() canvas.SetTicky() - psFileName = "%s_%s.ps"%(sampleName, org.tag) - canvas.Print(psFileName+"[","Lanscape") - - assert len(self.parameters()["objects"])==1 - for key,value in self.parameters()["objects"].iteritems() : + psFileName = "%s_%s.ps" % (sampleName, org.tag) + canvas.Print(psFileName+"[", "Lanscape") + + assert len(self.parameters()["objects"]) == 1 + for key, value in self.parameters()["objects"].iteritems(): jet = value["jet"] - - for variable in ["nEvents"] : + + for variable in ["nEvents"]: histos = numerAndDenom(org, variable) - if "before" not in histos or "after" not in histos : continue + if "before" not in histos or "after" not in histos: + continue result = histos["after"].Clone(variable) result.Divide(histos["before"]) result.SetMarkerStyle(20) result.SetStats(False) - if result.ClassName()[2]=="1" : - result.GetYaxis().SetRangeUser(0.0,0.35) + if result.ClassName()[2] == "1": + result.GetYaxis().SetRangeUser(0.0, 0.35) result.GetYaxis().SetTitle("efficiency") result.Draw() - else : - result.GetZaxis().SetRangeUser(0.0,0.35) + else: + result.GetZaxis().SetRangeUser(0.0, 0.35) result.GetZaxis().SetTitle("efficiency") result.Draw("colz") - canvas.Print(psFileName,"Lanscape") + canvas.Print(psFileName, "Lanscape") result.Write() - canvas.Print(psFileName+"]","Lanscape") - temp = psFileName.replace(".chr","_chr").replace(".ps",".pdf") - os.system("ps2pdf "+ psFileName +" "+ temp) + canvas.Print(psFileName+"]", "Lanscape") + temp = psFileName.replace(".chr", "_chr").replace(".ps", ".pdf") + cmd = "ps2pdfwr " + psFileName + " " + temp + os.system(cmd) os.remove(psFileName) file.Close() - - From 16b5cf55f1ca30cb0a495254246527aa65aa239e Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Tue, 2 Jul 2013 10:31:19 -0500 Subject: [PATCH 09/34] HTMHTParked --- analyses/jsonMaker.py | 59 +++++++++++++++++-------------------------- 1 file changed, 23 insertions(+), 36 deletions(-) diff --git a/analyses/jsonMaker.py b/analyses/jsonMaker.py index 8107e14..71e62ab 100644 --- a/analyses/jsonMaker.py +++ b/analyses/jsonMaker.py @@ -5,29 +5,27 @@ def parameters(self) : jwPrompt = calculables.other.jsonWeight("cert/Cert_160404-178677_7TeV_PromptReco_Collisions11_JSON.sub.txt") jwMay = calculables.other.jsonWeight("cert/Cert_160404-163869_7TeV_May10ReReco_Collisions11_JSON_v3.txt") jwAug = calculables.other.jsonWeight("cert/Cert_170249-172619_7TeV_ReReco5Aug_Collisions11_JSON_v3.txt") - - jw2012 = calculables.other.jsonWeight("cert/Cert_190456-196531_8TeV_PromptReco_Collisions12_JSON.txt") - jw2012_v2 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt") + jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt") group = self.vary() - group['SingleMu'] = [(["SingleMu.Run2012A-13Jul2012-v1.AOD.job358", - "SingleMu.Run2012A-13Jul2012-v1.AOD.job467", - - "SingleMu.Run2012A-recover-06Aug2012-v1.AOD.job359", - "SingleMu.Run2012A-recover-06Aug2012-v1.AOD.job477", - - "SingleMu.Run2012B-13Jul2012-v1.AOD.job358", - "SingleMu.Run2012B-13Jul2012-v1.AOD.job461", - "SingleMu.Run2012C-24Aug2012-v1.AOD.job361", - "SingleMu.Run2012C-24Aug2012-v1.AOD.job470", - - "SingleMu.Run2012C-PromptReco-v2.AOD.job360", - "SingleMu.Run2012C-PromptReco-v2.AOD.job747", - "SingleMu.Run2012D-PromptReco-v1.AOD.job508", - "SingleMu.Run2012D-PromptReco-v1.AOD.job525" - ], jw2012_v2)] - +# group['SingleMu'] = [(["SingleMu.Run2012A-13Jul2012-v1.AOD.job358", +# "SingleMu.Run2012A-13Jul2012-v1.AOD.job467", +# +# "SingleMu.Run2012A-recover-06Aug2012-v1.AOD.job359", +# "SingleMu.Run2012A-recover-06Aug2012-v1.AOD.job477", +# +# "SingleMu.Run2012B-13Jul2012-v1.AOD.job358", +# "SingleMu.Run2012B-13Jul2012-v1.AOD.job461", +# "SingleMu.Run2012C-24Aug2012-v1.AOD.job361", +# "SingleMu.Run2012C-24Aug2012-v1.AOD.job470", +# +# "SingleMu.Run2012C-PromptReco-v2.AOD.job360", +# "SingleMu.Run2012C-PromptReco-v2.AOD.job747", +# "SingleMu.Run2012D-PromptReco-v1.AOD.job508", +# "SingleMu.Run2012D-PromptReco-v1.AOD.job525" +# ], jw2012_v2)] +# # group['SingleEl'] = [(['SingleElectron.Run2012A-PromptReco-v1.AOD.job229', # 'SingleElectron.Run2012B-PromptReco-v1.AOD.job228', # 'SingleElectron.Run2012B-PromptReco-v1.AOD.job238', @@ -43,22 +41,11 @@ def parameters(self) : # 'DoubleMu.Run2012B-PromptReco-v1.AOD.job239', # ], jw2012)] # -# group['HT1'] = [(['HT.Run2012A-PromptReco-v1.AOD.job229', -# 'HTMHT.Run2012B-PromptReco-v1.AOD.job228', -# 'HTMHT.Run2012B-PromptReco-v1.AOD.job238', -# 'JetHT.Run2012B-PromptReco-v1.AOD.job228', -# 'JetHT.Run2012B-PromptReco-v1.AOD.job238', -# ], jw2012)] -# -# group['HT2'] = [('HT.Run2011A-May10ReReco-v1.AOD.Darren1',jwMay), -# ('HT.Run2011A-05Aug2011-v1.AOD.Bryn1',jwAug), -# (['HT.Run2011A-PromptReco-v4.AOD.Bryn1', -# 'HT.Run2011A-PromptReco-v6.AOD.Bryn1', -# 'HT.Run2011B-PromptReco-v1.AOD.Bryn1', -# 'HT.Run2011B-PromptReco-v1.AOD.Bryn2', -# 'HT.Run2011B-PromptReco-v1.AOD.Bryn3', -# ], jwPrompt)] -# + group['HT1'] = [(["HTMHTParked.Run2012B-22Jan2013-v1.job649", + "HTMHTParked.Run2012C-22Jan2013-v1.job649", + "HTMHTParked.Run2012D-22Jan2013-v1.job649"], jw2012)] + + return {'group':group} def listOfSteps(self,pars) : From ae3f9f1d75595418aeed7e30310827ef31330b7d Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Tue, 2 Jul 2013 10:32:19 -0500 Subject: [PATCH 10/34] adding zinv samples --- samples/__ewk17__.py | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/samples/__ewk17__.py b/samples/__ewk17__.py index 7d40ce2..9dc25c7 100644 --- a/samples/__ewk17__.py +++ b/samples/__ewk17__.py @@ -46,3 +46,10 @@ ewk17.add("zinv_mg_ht_200_400.job365", '%s/clucas//ICF/automated/2012_09_21_10_32_59/%s")'%(pnfs, Zinv_S10_excl)%(200,400), xs = {"NNLO":49.2776, "LO":41.49}["NNLO"]) ewk17.add("zinv_mg_ht_400_inf.job365", '%s/clucas//ICF/automated/2012_09_21_10_32_59/%s")'%(pnfs, Zinv_S10_excl)%(400,"inf"), xs = {"NNLO":6.2639, "LO":5.274}["NNLO"]) +ext_suffix = "ZJetsToNuNu_%s_HT_%s_TuneZ2Star_8TeV_madgraph_ext.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM/" +ewk17.add("zinv_mg_ht_50_100_ext.job500", '%s/karage//ICF/automated/2012_11_29_23_55_17/%s")' % (pnfs, ext_suffix % (50, 100)), xs={"NNLO":452.75, "LO":381.2}["NNLO"]) +ewk17.add("zinv_mg_ht_100_200_ext.job680", '%s/clucas//ICF/automated/2013_06_17_17_08_07/%s")' % (pnfs, + ext_suffix.replace("V7A","V7C") % (100, 200)),xs={"NNLO":190.39, + "LO":160.3}["NNLO"]) +ewk17.add("zinv_mg_ht_200_400_ext.job500", '%s/karage//ICF/automated/2012_11_29_23_55_17/%s")' % (pnfs, ext_suffix % (200, 400)), xs={"NNLO":49.2776, "LO":41.49}["NNLO"]) +ewk17.add("zinv_mg_ht_400_inf_ext.job500", '%s/karage//ICF/automated/2012_11_29_23_55_17/%s")' % (pnfs, ext_suffix % (400, "inf")), xs={"NNLO":6.2639, "LO":5.274}["NNLO"]) From 17778b3c1b21b1a4b0e8ee636bd11947e03d24fb Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Tue, 2 Jul 2013 10:33:06 -0500 Subject: [PATCH 11/34] adding lumi to parked data --- samples/__ht17__.py | 7 ++++--- 1 file changed, 4 insertions(+), 3 deletions(-) diff --git a/samples/__ht17__.py b/samples/__ht17__.py index f50e736..60608c0 100644 --- a/samples/__ht17__.py +++ b/samples/__ht17__.py @@ -13,6 +13,7 @@ ht17.add("HT.Run2012A-13Jul2012-v1.AOD.job358", '%s/yeshaq//ICF/automated/2012_09_19_01_03_03/HT.Run2012A-13Jul2012-v1.AOD/")'%pnfs, lumi = 1.0) ht17.add("HT.Run2012A-recover-06Aug2012-v1.AOD.job359", '%s/yeshaq//ICF/automated/2012_09_19_01_25_46/HT.Run2012A-recover-06Aug2012-v1.AOD/")'%pnfs, lumi = 1.0) -for era in ["B","C","D"] : - ht17.add("HTMHTParked.Run2012%s-22Jan2013-v1.job649"%era, '%s/yeshaq//ICF/automated/2013_04_19_17_30_49/HTMHTParked.Run2012%s-22Jan2013-v1.AOD/")'%(pnfs,era), lumi = 1.0) - +for era,lum in zip(["B","C","D"],[4427.0,6893.0,7263.0]) : + ht17.add("HTMHTParked.Run2012%s-22Jan2013-v1.job649"%era, '%s/yeshaq//ICF/automated/2013_04_19_17_30_49/HTMHTParked.Run2012%s-22Jan2013-v1.AOD/")'%(pnfs,era), lumi=lum) + +ht17.add("HTMHTParked_ICF_sync_test", '%s/yeshaq//ICF/supy-output/ICF_sync_test/")'%pnfs, lumi = 1.0) From 36d307246aee830e72c4a6dd5a7462c958a6b862 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Tue, 2 Jul 2013 10:33:49 -0500 Subject: [PATCH 12/34] exending single jet pt histogram range --- steps/jet.py | 2 +- 1 file changed, 1 insertion(+), 1 deletion(-) diff --git a/steps/jet.py b/steps/jet.py index d3a2d0c..853319f 100644 --- a/steps/jet.py +++ b/steps/jet.py @@ -305,7 +305,7 @@ def uponAcceptance (self,eventVars) : self.book.fill(eta2, "%s%s%sEta2mom" %(self.cs+(jetLabel,)), 50, 0.0, mom2Max, title=";jet%s #sigma_{#eta}^{2};events / bin"%jetLabel) self.book.fill(phi2, "%s%s%sPhi2mom" %(self.cs+(jetLabel,)), 50, 0.0, mom2Max, title=";jet%s #sigma_{#phi}^{2};events / bin"%jetLabel) - self.book.fill(pt, "%s%s%sPt" %(self.cs+(jetLabel,)), 50, 0.0, 500.0, title=";jet%s p_{T} (GeV);events / bin"%jetLabel) + self.book.fill(pt, "%s%s%sPt" %(self.cs+(jetLabel,)), 100, 0.0, 1000.0, title=";jet%s p_{T} (GeV);events / bin"%jetLabel) self.book.fill(eta, "%s%s%seta"%(self.cs+(jetLabel,)), 50, -5.0, 5.0, title=";jet%s #eta;events / bin"%jetLabel) if i>self.maxIndex: continue for j,jJet in list(enumerate(cleanJetIndices))[i+1:self.maxIndex+1] : From 3026240f329566198c7fd6be60aa4fc14fe801f6 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Tue, 2 Jul 2013 10:37:15 -0500 Subject: [PATCH 13/34] removed extraneous char --- analyses/hadronicLook.py | 6 +++--- 1 file changed, 3 insertions(+), 3 deletions(-) diff --git a/analyses/hadronicLook.py b/analyses/hadronicLook.py index 04f9660..cd538b0 100644 --- a/analyses/hadronicLook.py +++ b/analyses/hadronicLook.py @@ -443,7 +443,7 @@ def listOfSamples(self, params): from supy.samples import specify def data_53X(): - jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptRecon_Collisions12_JSON.txt") + jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt") out = [] out += specify(names="HTMHTParked.Run2012B-22Jan2013-v1.job649", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) out += specify(names="HTMHTParked.Run2012C-22Jan2013-v1.job649", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) @@ -535,8 +535,8 @@ def sms(): return ( data_53X() + -# w_binned() + -# z_binned() + + w_binned() + + z_binned() + #top() + #vv() + #qcd_py6(30.0e3) + From 352b5f2aa9ae2671c55df83edf340ce5c61602a9 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Mon, 8 Jul 2013 03:06:45 -0500 Subject: [PATCH 14/34] updating rest of SM MC --- analyses/hadronicLook.py | 84 ++++++++++++++++++---------------------- 1 file changed, 37 insertions(+), 47 deletions(-) diff --git a/analyses/hadronicLook.py b/analyses/hadronicLook.py index cd538b0..2781651 100644 --- a/analyses/hadronicLook.py +++ b/analyses/hadronicLook.py @@ -35,13 +35,13 @@ def parameters(self): {"muon": ("muon", "Pat"), "electron": ("electron", "Pat"), "photon": ("photon", "Pat"), - + "jet": ("xcak5Jet", "Pat"), "jetId": "JetIDloose", "muonsInJets": False, "met": "metP4TypeIPF", "rechit": "Calo", - + "jetComp": ("xcak5JetPF", "Pat"), "jetIdComp": "JetIDtight", "muonsInJetsComp": True, @@ -308,7 +308,6 @@ def stepsPlotsOne(self, params): steps.jet.cleanJetHtMhtHistogrammer(_jet, params["etRatherThanPt"]), supy.steps.histos.histogrammer("%sDeltaPhiStar%s%s" % (_jet[0], params["lowPtName"], _jet[1]), 20, 0.0, r.TMath.Pi(), title=";#Delta#phi*;events / bin", funcString='lambda x:x[0][0]'), supy.steps.histos.histogrammer("%sDeltaPhiStar%s" % (_jet[0], _jet[1]), 20, 0.0, r.TMath.Pi(), title=";#Delta#phi*;events / bin", funcString='lambda x:x[0][0]'), - supy.steps.histos.histogrammer("%sMaxEmEnergyFraction%s" % (_jet[0], _jet[1]), 20, 0.0, 1.0, title=";MaxEmEnergyFraction;events / bin"), supy.steps.histos.histogrammer(_met, 100, 0.0, 500.0, title=";"+_met+" (GeV);events / bin", funcString="lambda x: x.pt()"), supy.steps.filters.label("kinematicPlots1"), @@ -358,8 +357,8 @@ def stepsHtBins(self, params): def stepsOptional(self, params): return [ - #supy.steps.other.skimmer(), - steps.other.duplicateEventCheck(), + supy.steps.other.skimmer(), + #steps.other.duplicateEventCheck(), #steps.other.pickEventSpecMaker(), #steps.other.cutBitHistogrammer(self.togglePfJet(_jet), self.togglePfMet(_met)), #steps.Print.jetPrinter(_jet), @@ -423,12 +422,11 @@ def listOfSteps(self, params): self.stepsXclean(params) + self.stepsOptional(params) + self.stepsBtagJets(params) + - self.stepsPlotsOne(params) + self.stepsQcdRejection(params) + + self.stepsPlotsOne(params) + self.stepsPlotsTwo(params) + #self.stepsMbb(params) + #self.stepsDisplayer(params) + - #self.stepsHtBins(params) + #self.stepsEventCount(params, label="scanAfter") + []) @@ -437,7 +435,7 @@ def listOfSampleDictionaries(self): sh = supy.samples.SampleHolder() sh.add("275_ge2b", '["/uscms/home/elaird/08_mbb/02_skim/2012_5fb_275_ge2b.root"]', lumi=5.0e3) sh.add("375_ge2b", '["/uscms/home/yeshaq/nobackup/supy-output/hadronicLook/375_calo_ge2/HadronicRegion.root"]', lumi=5.0e3) - return [samples.ht17, samples.top17, samples.ewk17, samples.qcd17, sh, samples.susy17] + return [samples.ht17, samples.top17, samples.ewk17, samples.qcd17, sh, samples.susy17, samples.photon17] def listOfSamples(self, params): from supy.samples import specify @@ -445,9 +443,9 @@ def listOfSamples(self, params): def data_53X(): jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt") out = [] - out += specify(names="HTMHTParked.Run2012B-22Jan2013-v1.job649", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) - out += specify(names="HTMHTParked.Run2012C-22Jan2013-v1.job649", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) - out += specify(names="HTMHTParked.Run2012D-22Jan2013-v1.job649", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="HTMHTParked.Run2012B-22Jan2013-v1.job649", weights=jw2012) + out += specify(names="HTMHTParked.Run2012C-22Jan2013-v1.job649", weights=jw2012) + out += specify(names="HTMHTParked.Run2012D-22Jan2013-v1.job649", weights=jw2012) #out += specify(names="HTMHTParked_ICF_sync_test", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) return out @@ -455,7 +453,7 @@ def qcd_py6(eL): low = map(lambda x: x[0], samples.__qcd17__.binsXs)[:-1] out = [] for pt in low[low.index(80):]: - out += specify("qcd_py6_pt_%d" % pt, effectiveLumi=eL) + out += specify("qcd_py6_pt_v2_%d" % pt, effectiveLumi=eL) return out def qcd_b_py6(eL): @@ -465,59 +463,50 @@ def qcd_b_py6(eL): return out def g_jets_mg(eL): - gM = [40, 100, 200] - return specify(effectiveLumi=eL, color=r.kGreen, - names=[("g_jets_mg_ht_%d_%d")[:None if t[1] else -2] for t in zip(gM, gM[1:]+["inf"])]) + out = [] + out += specify("g_jets_mg_ht_200_400.job501", effectiveLumi=eL, color=r.kGreen) + out += specify("g_jets_mg_ht_400_inf.job501", effectiveLumi=eL, color=r.kGreen) + return out def w_binned(): out = [] - out += specify(names="wj_lv_mg_ht_10_150", color=r.kBlue) # , nFilesMax=1, nEventsMax=20000) - out += specify(names="wj_lv_mg_ht_150_200.job663", color=r.kGreen) # , nFilesMax=1, nEventsMax=20000) - out += specify(names="wj_lv_mg_ht_200_250.job672", color=r.kCyan) # , nFilesMax=1, nEventsMax=20000) - out += specify(names="wj_lv_mg_ht_250_300.job498", color=r.kOrange) # , nFilesMax=1, nEventsMax=20000) - out += specify(names="wj_lv_mg_ht_300_400.job498", color=r.kViolet) # , nFilesMax=1, nEventsMax=20000) - out += specify(names="wj_lv_mg_ht_400_inf.job498", color=r.kAzure) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="wj_lv_mg_ht_10_150", color=r.kBlue) + out += specify(names="wj_lv_mg_ht_150_200.job663", color=r.kGreen) + out += specify(names="wj_lv_mg_ht_200_250.job672", color=r.kCyan) + out += specify(names="wj_lv_mg_ht_250_300.job498", color=r.kOrange) + out += specify(names="wj_lv_mg_ht_300_400.job498", color=r.kViolet) + out += specify(names="wj_lv_mg_ht_400_inf.job498", color=r.kAzure) return out def z_binned(): out = [] - out += specify("zinv_mg_ht_50_100.job407", color=r.kBlue) # , nFilesMax=1, nEventsMax=20000) - out += specify("zinv_mg_ht_100_200.job365", color=r.kGreen) # , nFilesMax=1, nEventsMax=20000) - out += specify("zinv_mg_ht_200_400.job365", color=r.kOrange) # , nFilesMax=1, nEventsMax=20000) - out += specify("zinv_mg_ht_400_inf.job365", color=r.kViolet) # , nFilesMax=1, nEventsMax=20000) - out += specify("zinv_mg_ht_50_100_ext.job500", color=r.kBlue) # , nFilesMax=1, nEventsMax=20000) - out += specify("zinv_mg_ht_100_200_ext.job680", color=r.kGreen) # , nFilesMax=1, nEventsMax=20000) - out += specify("zinv_mg_ht_200_400_ext.job500", color=r.kOrange) # , nFilesMax=1, nEventsMax=20000) - out += specify("zinv_mg_ht_400_inf_ext.job500", color=r.kViolet) # , nFilesMax=1, nEventsMax=20000) + out += specify("zinv_mg_ht_50_100.job407", color=r.kBlue) + out += specify("zinv_mg_ht_100_200.job365", color=r.kGreen) + out += specify("zinv_mg_ht_200_400.job365", color=r.kOrange) + out += specify("zinv_mg_ht_400_inf.job365", color=r.kViolet) + out += specify("zinv_mg_ht_50_100_ext.job500", color=r.kBlue) + out += specify("zinv_mg_ht_100_200_ext.job680", color=r.kGreen) + out += specify("zinv_mg_ht_200_400_ext.job500", color=r.kOrange) + out += specify("zinv_mg_ht_400_inf_ext.job500", color=r.kViolet) return out def w_inclusive(): out = [] for part in range(1, 6): - out += specify(names="wj_lv_mg_ht_incl_v2.job673_part%i" % part, nFilesMax=1, nEventsMax=1000) + out += specify(names="wj_lv_mg_ht_incl_v2.job673_part%i" % part) return out def vv(): out = [] - out += specify("ww_py.job188") - out += specify("wz_py.job188") - out += specify("zz_py.job188") - out += specify("zinv_hbb_125_powheg.job342") + #out += specify("zinv_hbb_125_powheg.job342") + for diBos in ["ZZ", "WZ", "WW"]: out += specify("%s_pythia6.job370" % diBos) return out def top(): out = [] out += specify(names="ttbar_powheg_v1.job410") - #out += specify(names="tt_8_mg.job188") - #out += specify(names="tt_8_mg.job315") - out += specify(names="ttz_8_mg.job269", nFilesMax=1) - - out += specify("t_s_powheg.job200") - #oAut += specify("t_t_powheg.job187" ) #low MC stats - out += specify("t_tw_powheg.job187") - out += specify("tbar_t_powheg.job194") - out += specify("tbar_tw_powheg.job187") - + for sTop in ["T_s","T_t", "T_tW", "Tbar_s", "Tbar_t", "Tbar_tW"]: + out += specify("%s_powheg.job368" % sTop) return out def susy(eL): @@ -535,12 +524,13 @@ def sms(): return ( data_53X() + - w_binned() + - z_binned() + + #w_binned() + + #z_binned() + #top() + #vv() + #qcd_py6(30.0e3) + #qcd_b_py6(30.0e3) + + #g_jets_mg(30e3) + #w_inclusive() + #sms() + [] @@ -581,7 +571,7 @@ def md(x, y): org.mergeSamples(targetSpec=md({"name": "Standard Model ", "color": r.kAzure+6}, mcOps), sources=ewkSources + qcdSources, keepSources=True) def conclude(self, conf): - org = self.organizer(conf, verbose=True) + org = self.organizer(conf) self.mergeSamples(org) ##for skimming only #utils.printSkimResults(org) From d37711ee78095ae48ff7d1ea1b0cb3e0925fda68 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Mon, 8 Jul 2013 03:08:43 -0500 Subject: [PATCH 15/34] adding VV samples --- samples/__ewk17__.py | 5 +++++ 1 file changed, 5 insertions(+) diff --git a/samples/__ewk17__.py b/samples/__ewk17__.py index 9dc25c7..b04c4df 100644 --- a/samples/__ewk17__.py +++ b/samples/__ewk17__.py @@ -53,3 +53,8 @@ "LO":160.3}["NNLO"]) ewk17.add("zinv_mg_ht_200_400_ext.job500", '%s/karage//ICF/automated/2012_11_29_23_55_17/%s")' % (pnfs, ext_suffix % (200, 400)), xs={"NNLO":49.2776, "LO":41.49}["NNLO"]) ewk17.add("zinv_mg_ht_400_inf_ext.job500", '%s/karage//ICF/automated/2012_11_29_23_55_17/%s")' % (pnfs, ext_suffix % (400, "inf")), xs={"NNLO":6.2639, "LO":5.274}["NNLO"]) + +diBos = "%s_TuneZ2star_8TeV_pythia6_tauola.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM/" +diBosons ={"ZZ":8.25561, "WZ":32.3161, "WW":57.1097} +for key in diBosons: + ewk17.add("%s_pythia6.job370" % key, '%s/zmeng/ICF/automated/2012_09_21_17_43_41/%s")' % (pnfs, diBos % key), xs = diBosons[key]) From abc5dea605277f11cbbbb6a51b16087b02f9ad9b Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Mon, 8 Jul 2013 03:11:09 -0500 Subject: [PATCH 16/34] updating gJets MC, cleaning up file --- samples/__photon17__.py | 17 +++++++---------- 1 file changed, 7 insertions(+), 10 deletions(-) diff --git a/samples/__photon17__.py b/samples/__photon17__.py index d90ab30..064a1cf 100644 --- a/samples/__photon17__.py +++ b/samples/__photon17__.py @@ -3,20 +3,17 @@ pnfs = pnfs() photon17 = SampleHolder() -#photon17.add("Photon.Run2012A-PromptReco-v1.AOD.job29", '%s/zmeng//ICF/automated/2012_04_13_19_13_29/")'%pnfs, lumi = 1.0) -#photon17.add("Photon.Run2012A-PromptReco-v1.AOD.job44", '%s/zmeng//ICF/automated/2012_04_16_09_55_55/")'%pnfs, lumi = 1.0) -#photon17.add("Photon.Run2012A-PromptReco-v1.AOD.job57", '%s/zmeng//ICF/automated/2012_04_22_16_58_56/")'%pnfs, lumi = 1.0) -#photon17.add("Photon.Run2012A-PromptReco-v1.AOD.job69", '%s/zmeng//ICF/automated/2012_04_26_16_22_08/")'%pnfs, lumi = 1.0) -#photon17.add("Photon.Run2012A-PromptReco-v1.AOD.job74", '%s/zmeng//ICF/automated/2012_04_27_00_07_16/")'%pnfs, lumi = 1.0) -#photon17.add("Photon.Run2012A-PromptReco-v1.AOD.job81", '%s/clucas//ICF/automated/2012_05_06_13_23_30/")'%pnfs,lumi = 1.0) - a = "%s/yeshaq/ICF/automated/2012_05_25_16_29_39"%pnfs photon17.add("Photon.Run2012A-PromptReco-v1.AOD.job171", '%s/Photon.Run2012A-PromptReco-v1.AOD")'%a, lumi = 1.0) photon17.add("SinglePhoton.Run2012B-PromptReco-v1.AOD.job171", '%s/SinglePhoton.Run2012B-PromptReco-v1.AOD")'%a, lumi = 1.0) -photon17.add("GJets_HT400.job92", '%s/clucas//ICF/automated/2012_05_08_11_07_51/")'%pnfs, xs = {"LO":107.5}["LO"]) -photon17.add("GJets_HT400.job174", '%s/yeshaq//ICF/automated/2012_05_30_22_41_51/")'%pnfs, xs = {"LO":107.5}["LO"]) - photon17.add("Photon.Run2012A-PromptReco-v1.AOD.job228", '%s/zmeng//ICF/automated/2012_06_14_11_22_04/Photon.Run2012A-PromptReco-v1.AOD")'%pnfs, lumi = 1.0) photon17.add("SinglePhoton.Run2012B-PromptReco-v1.AOD.job229", '%s/zmeng//ICF/automated/2012_06_14_14_26_12/SinglePhoton.Run2012B-PromptReco-v1.AOD")'%pnfs, lumi = 1.0) photon17.add("SinglePhoton.Run2012B-PromptReco-v1.AOD.job238", '%s/zmeng//ICF/automated/2012_06_22_14_25_16/SinglePhoton.Run2012B-PromptReco-v1.AOD")'%pnfs, lumi = 1.0) + +gJet = "GJets_HT-%s_8TeV-madgraph_v2.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM/" + +photon17.add("g_jets_mg_ht_200_400.job501", '%s/agapitos//ICF/automated/2012_11_29_18_15_22/%s")' % (pnfs, gJet % "200To400"), + xs={"LO": 960.5, "NLO":1140.78}["NLO"]) +photon17.add("g_jets_mg_ht_400_inf.job501", '%s/agapitos//ICF/automated/2012_11_29_18_15_22/%s")' % (pnfs, gJet % "400ToInf"), + xs={"LO": 107.5, "NLO":124.68}["NLO"]) From e00d2f0995e55e4737fd1fe9e3d750171600417d Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Mon, 8 Jul 2013 03:11:48 -0500 Subject: [PATCH 17/34] updating QCD sample --- samples/__qcd17__.py | 7 +++++++ 1 file changed, 7 insertions(+) diff --git a/samples/__qcd17__.py b/samples/__qcd17__.py index 815e1ef..a16d203 100644 --- a/samples/__qcd17__.py +++ b/samples/__qcd17__.py @@ -33,3 +33,10 @@ qcd17.add("qcd_b_py6_pt_30", '%s/QCD_Pt-30To50_bEnriched_TuneZ2star_8TeV-pythia6-evtgen.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM")'%bEn, xs = 6.677E7) qcd17.add("qcd_b_py6_pt_50", '%s/QCD_Pt-50To150_bEnriched_TuneZ2star_8TeV-pythia6-evtgen.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM")'%bEn, xs = 9355000.0) qcd17.add("qcd_b_py6_pt_150", '%s/QCD_Pt-150_bEnriched_TuneZ2star_8TeV-pythia6-evtgen.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM")'%bEn, xs = 67340.0) + +for i,(low,xs) in enumerate(binsXs[:-1]) : + hi = binsXs[i+1][0] + dir = "_.QCD_Pt-%d%s_TuneZ2star_8TeV_pythia6.Summer12-PU_S7_START52_V9-v1.AODSIM_/"%(low, "to%d"%hi if hi else "") + if i < 1 : dir = dir.replace("_.QCD", "QCD").replace("AODSIM_","AODSIM") + if low == 1800.: dir = dir.replace("AODSIM_","AODSIM") + qcd17.add("qcd_py6_pt_v2_%d"%low, '%s/agapitos//ICF/automated/2012_09_22_14_28_53/%s")'%(pnfs, dir), xs = xs) From 51ffae3168606055340a579c4de5cf67f0a50f5d Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Mon, 8 Jul 2013 03:13:25 -0500 Subject: [PATCH 18/34] updating single top samples, removing old samples --- samples/__top17__.py | 29 +++++++++++++---------------- 1 file changed, 13 insertions(+), 16 deletions(-) diff --git a/samples/__top17__.py b/samples/__top17__.py index 2430216..91e5f4d 100644 --- a/samples/__top17__.py +++ b/samples/__top17__.py @@ -1,31 +1,28 @@ from supy.samples import SampleHolder from supy.sites import pnfs +import itertools pnfs = pnfs() top17 = SampleHolder() -top17.add("tt_8_mg.job188", '%s/clucas//ICF/automated/2012_06_05_09_30_46/TTJets_TuneZ2star_8TeV-madgraph-tauola.Summer12-PU_S7_START52_V9-v1.AODSIM")'%pnfs, xs = 234.) - -top17.add("tt_8_mg.job315", '%s/clucas/ICF/automated/2012_08_17_14_42_02/")'%pnfs, xs = 234.) top17.add("tt_8_mg.job315_1", '["/uscms/home/elaird/141_ntuples/job315/SusyCAF_Tree_100_1_h3D.root"]', xs = 234.) top17.add("tt_8_mg.job315_zeroNu_50Met", '["/uscms/home/elaird/141_ntuples/job315/tt_8_mg.job315_zeroNu_50Met.root"]', xs = 234.*(44005./6923750)) top17.add("tt_8_mg.job315_1_oneNu_100Met", '["/uscms/home/elaird/141_ntuples/job315/tt_8_mg.job315_oneNu_100Met.root"]', xs = 234.*(1465./20000.)) top17.add("tt_8_mg.job315_1_twoNu_100Met", '["/uscms/home/elaird/141_ntuples/job315/tt_8_mg.job315_twoNu_100Met.root"]', xs = 234.*(510./20000.)) loc = "/uscms/home/elaird/141_ntuples/job315" -top17.add("tt_8_mg.job315_oneNu_100Met", 'utils.fileListFromDisk(location = "%s/oneNu_100Met/tt_8_mg.job315_*_skim.root", isDirectory = False)'%loc, xs = 7.037891e-02 * 2.340000e+02) -top17.add("tt_8_mg.job315_twoNu_100Met", 'utils.fileListFromDisk(location = "%s/twoNu_100Met/tt_8_mg.job315_*_skim.root", isDirectory = False)'%loc, xs = 2.542509e-02 * 2.340000e+02) - -top17.add("ttz_8_mg.job269", '%s/zmeng//ICF/automated/2012_07_20_15_24_50")'%pnfs, xs = 0.172) #xs from PREP +top17.add("tt_8_mg.job315_oneNu_100Met", 'utils.fileListFromDisk(location = "%s/oneNu_100Met/tt_8_mg.job315_*_skim.root", isDirectory = False)'%loc, + xs = 7.037891e-02 * 2.340000e+02) +top17.add("tt_8_mg.job315_twoNu_100Met", 'utils.fileListFromDisk(location = "%s/twoNu_100Met/tt_8_mg.job315_*_skim.root", isDirectory = False)'%loc, + xs = 2.542509e-02 * 2.340000e+02) top17.add("ttz_8_mg.job269_1", '["/uscms/home/elaird/141_ntuples/job269/SusyCAF_Tree_11_1_7LX.root"]', xs = 0.172) #xs from PREP - -y = '%s/yeshaq//ICF/automated/2012_06_05_02_07_12/'%pnfs -top17.add("t_s_powheg.job200", '%s/clucas//ICF/automated/2012_06_06_18_49_32/")'%pnfs, xs = {"NNLO":5.55, "LO":2.82}["NNLO"]) -top17.add("t_t_powheg.job187", '%s/T_t-channel_TuneZ2star_8TeV-powheg-tauola.Summer12-PU_S7_START52_V9-v1.AODSIM")'%y, xs = {"NNLO":56.4, "LO":47.0}["NNLO"]) -top17.add("t_tw_powheg.job187", '%s/T_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola.Summer12-PU_S7_START52_V9-v1.AODSIM")'%y, xs = {"NNLO":11.1, "LO":10.7}["NNLO"]) - -top17.add("tbar_t_powheg.job194", '%s/clucas//ICF/automated/2012_06_06_11_53_47/")'%pnfs, xs = {"NNLO":30.7, "LO":25.0}["NNLO"]) -top17.add("tbar_tw_powheg.job187", '%s/Tbar_tW-channel-DR_TuneZ2star_8TeV-powheg-tauola.Summer12-PU_S7_START52_V9-v1.AODSIM")'%y, xs = {"NNLO":11.1, "LO":10.7}["NNLO"]) - top17.add("ttbar_powheg_v1.job410", '%s/clucas//ICF/automated/2012_09_26_14_02_42/")'%pnfs, xs = 211.) top17.add("ttbar_powheg_v2.job404", '%s/zmeng//ICF/automated/2012_09_21_17_07_53/")'%pnfs, xs = 211.) + +sTop = "-channel_TuneZ2star_8TeV-powheg-tauola.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM/" +combs = {"T_s":{"NLO":3.79, "LO":2.82},"T_t":{"NLO":56.4, "LO":47.0},"T_tW":{"NLO":11.1, "LO":10.7}, + "Tbar_s":{"NLO":1.76, "LO":1.57},"Tbar_t": {"NLO":30.7, "LO":25.0},"Tbar_tW":{"NLO":11.1, "LO":10.7}} + +for key in combs: + top17.add("%s_powheg.job368"%key, + '%s/zmeng//ICF/automated/2012_09_21_17_19_15/%s")'%(pnfs, (key + sTop).replace("tW-channel","tW-channel-DR")),combs[key]["NLO"]) From 66e5fca4ae8ef924ef0db651d166370bad9004a5 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Mon, 8 Jul 2013 03:18:24 -0500 Subject: [PATCH 19/34] comforming to pep8 as best as possible for now --- analyses/hadronicLook.py | 11 ++++++----- 1 file changed, 6 insertions(+), 5 deletions(-) diff --git a/analyses/hadronicLook.py b/analyses/hadronicLook.py index 2781651..b1f9e27 100644 --- a/analyses/hadronicLook.py +++ b/analyses/hadronicLook.py @@ -35,13 +35,13 @@ def parameters(self): {"muon": ("muon", "Pat"), "electron": ("electron", "Pat"), "photon": ("photon", "Pat"), - + "jet": ("xcak5Jet", "Pat"), "jetId": "JetIDloose", "muonsInJets": False, "met": "metP4TypeIPF", "rechit": "Calo", - + "jetComp": ("xcak5JetPF", "Pat"), "jetIdComp": "JetIDtight", "muonsInJetsComp": True, @@ -59,7 +59,7 @@ def parameters(self): #"eq3b": (3, 3), #"ge4b": (4, None), }), - "thresholds": self.vary({#"200s": (200.0, 275.0, 73.3, 36.7), + "thresholds": self.vary({"200s": (200.0, 275.0, 73.3, 36.7), #"275s": (275.0, 325.0, 73.3, 36.7), #"325s": (325.0, 375.0, 86.7, 43.3), "375": (375.0, None, 100.0, 50.0), @@ -499,13 +499,14 @@ def w_inclusive(): def vv(): out = [] #out += specify("zinv_hbb_125_powheg.job342") - for diBos in ["ZZ", "WZ", "WW"]: out += specify("%s_pythia6.job370" % diBos) + for diBos in ["ZZ", "WZ", "WW"]: + out += specify("%s_pythia6.job370" % diBos) return out def top(): out = [] out += specify(names="ttbar_powheg_v1.job410") - for sTop in ["T_s","T_t", "T_tW", "Tbar_s", "Tbar_t", "Tbar_tW"]: + for sTop in ["T_s", "T_t", "T_tW", "Tbar_s", "Tbar_t", "Tbar_tW"]: out += specify("%s_powheg.job368" % sTop) return out From 2359e5e2875764efb188799f4fcab41e2731048c Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Wed, 4 Dec 2013 10:09:09 -0600 Subject: [PATCH 20/34] adding jan22 SingleMu --- samples/__muon17__.py | 9 ++++++++- 1 file changed, 8 insertions(+), 1 deletion(-) diff --git a/samples/__muon17__.py b/samples/__muon17__.py index 735b4ab..96a42c3 100644 --- a/samples/__muon17__.py +++ b/samples/__muon17__.py @@ -1,8 +1,10 @@ from supy.samples import SampleHolder -from supy.sites import pnfs +from supy.sites import pnfs, eos pnfs = pnfs() +eos = eos() muon17 = SampleHolder() + muon17.add("SingleMu.Run2012A-PromptReco-v1.AOD.job228", '%s/zmeng//ICF/automated/2012_06_14_11_22_04/SingleMu.Run2012A-PromptReco-v1.AOD/", alwaysUseLastAttempt = True)'%pnfs, lumi = 1 ) #/pb muon17.add("SingleMu.Run2012B-PromptReco-v1.AOD.job229", '%s/zmeng//ICF/automated/2012_06_14_14_26_12/SingleMu.Run2012B-PromptReco-v1.AOD/", alwaysUseLastAttempt = True)'%pnfs, lumi = 1 ) #/pb muon17.add("SingleMu.Run2012B-PromptReco-v1.AOD.job238", '%s/zmeng//ICF/automated/2012_06_22_14_25_16/SingleMu.Run2012B-PromptReco-v1.AOD/", alwaysUseLastAttempt = True)'%pnfs, lumi = 1 ) #/pb @@ -24,3 +26,8 @@ muon17.add("SingleMu.Run2012D-PromptReco-v1.AOD.job508", '%s/yeshaq/ICF/automated/2012_12_04_16_51_31/SingleMu.Run2012D-PromptReco-v1.AOD/")'%pnfs, lumi = 1 ) #/pb muon17.add("SingleMu.Run2012D-PromptReco-v1.AOD.job525", '%s/yeshaq/ICF/automated/2013_01_18_14_00_11/")'%pnfs, lumi = 1 ) #/pb + +muon17.add("SingleMu.Run2012A-22Jan2013", '%s/clucas/Parked13/SingleMu_Run2012A_22Jan2013/")'%eos, lumi = 1) #/pb +muon17.add("SingleMu.Run2012B-22Jan2013", '%s/clucas/Parked13/SingleMu_Run2012B_22Jan2013/")'%eos, lumi = 1) #/pb +muon17.add("SingleMu.Run2012C-22Jan2013", '%s/clucas/Parked13/SingleMu_Run2012C_22Jan2013/")'%eos, lumi = 1) #/pb +muon17.add("SingleMu.Run2012D-22Jan2013", '%s/clucas/Parked13/SingleMu_Run2012D_22Jan2013/")'%eos, lumi = 1) #/pb From 87f6c1fe4a337b8797297fa3cfa3ad46be13cbc3 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Wed, 4 Dec 2013 10:30:17 -0600 Subject: [PATCH 21/34] adding HT,HTMHT eos Jan22 --- samples/__ht17__.py | 7 ++++++- 1 file changed, 6 insertions(+), 1 deletion(-) diff --git a/samples/__ht17__.py b/samples/__ht17__.py index 60608c0..26125bb 100644 --- a/samples/__ht17__.py +++ b/samples/__ht17__.py @@ -1,6 +1,7 @@ from supy.samples import SampleHolder -from supy.sites import pnfs +from supy.sites import pnfs, eos pnfs = pnfs() +eos = eos() ht17 = SampleHolder() a = ", alwaysUseLastAttempt = True" @@ -17,3 +18,7 @@ ht17.add("HTMHTParked.Run2012%s-22Jan2013-v1.job649"%era, '%s/yeshaq//ICF/automated/2013_04_19_17_30_49/HTMHTParked.Run2012%s-22Jan2013-v1.AOD/")'%(pnfs,era), lumi=lum) ht17.add("HTMHTParked_ICF_sync_test", '%s/yeshaq//ICF/supy-output/ICF_sync_test/")'%pnfs, lumi = 1.0) + +for era,lum in zip(["B","C","D"],[4427.0,6893.0,7263.0]) : + ht17.add("HTMHTParked.Run2012%s-22Jan2013"%era, '%s/clucas/Parked13/HTMHTParked_Run2012%s_22Jan2013/")'%(eos,era), lumi=lum) +ht17.add("HT.Run2012A-22Jan2013", '%s/clucas/Parked13/HT_Run2012A_22Jan2013/")'%eos, lumi=lum) From 95290a68b13c251a275d9450b0baf50743718435 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Wed, 4 Dec 2013 11:05:43 -0600 Subject: [PATCH 22/34] adding eos photon samples --- samples/__photon17__.py | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/samples/__photon17__.py b/samples/__photon17__.py index 064a1cf..9db6e66 100644 --- a/samples/__photon17__.py +++ b/samples/__photon17__.py @@ -1,6 +1,7 @@ from supy.samples import SampleHolder -from supy.sites import pnfs +from supy.sites import pnfs, eos pnfs = pnfs() +eos = eos() photon17 = SampleHolder() a = "%s/yeshaq/ICF/automated/2012_05_25_16_29_39"%pnfs @@ -17,3 +18,8 @@ xs={"LO": 960.5, "NLO":1140.78}["NLO"]) photon17.add("g_jets_mg_ht_400_inf.job501", '%s/agapitos//ICF/automated/2012_11_29_18_15_22/%s")' % (pnfs, gJet % "400ToInf"), xs={"LO": 107.5, "NLO":124.68}["NLO"]) + +for era in ["B","C","D"]: + photon17.add("SinglePhotonParked.Run2012%s-22Jan2013"%era, '%s/clucas/Parked13/SinglePhoton_Run2012%s_22Jan2013/")' % (eos, era), lumi = 1.0) + +photon17.add("Photon.Run2012A-22Jan2013", '%s/clucas/Parked13/Photon_Run2012A_22Jan2013/")' % eos, lumi = 1.0) From 0e4f0683b4491ee116cf4af1eafb503bb1e6a4f7 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Wed, 4 Dec 2013 11:18:49 -0600 Subject: [PATCH 23/34] adding GJets eos samples --- samples/__photon17__.py | 2 ++ 1 file changed, 2 insertions(+) diff --git a/samples/__photon17__.py b/samples/__photon17__.py index 9db6e66..87cc1b4 100644 --- a/samples/__photon17__.py +++ b/samples/__photon17__.py @@ -23,3 +23,5 @@ photon17.add("SinglePhotonParked.Run2012%s-22Jan2013"%era, '%s/clucas/Parked13/SinglePhoton_Run2012%s_22Jan2013/")' % (eos, era), lumi = 1.0) photon17.add("Photon.Run2012A-22Jan2013", '%s/clucas/Parked13/Photon_Run2012A_22Jan2013/")' % eos, lumi = 1.0) +photon17.add("GJets_HT200to400", '%s/clucas/Parked13/GJets_200to400/")' % eos, xs={"LO": 960.5, "NLO":1140.78}["NLO"]) +photon17.add("GJets_HT400toinf", '%s/clucas/Parked13/GJets_400toinf/")' % eos, xs={"LO": 960.5, "NLO":1140.78}["NLO"]) From 97180ae93e0820ad3aa88a2b2f7d220e567adc1b Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Mon, 9 Dec 2013 08:46:09 -0600 Subject: [PATCH 24/34] adding eos dy samples --- samples/__dyll17__.py | 8 ++++++-- 1 file changed, 6 insertions(+), 2 deletions(-) diff --git a/samples/__dyll17__.py b/samples/__dyll17__.py index d8322f8..3cd12f6 100644 --- a/samples/__dyll17__.py +++ b/samples/__dyll17__.py @@ -1,6 +1,7 @@ from supy.samples import SampleHolder -from supy.sites import pnfs +from supy.sites import pnfs, eos pnfs = pnfs() +eos = eos() dyll17 = SampleHolder() dyChris = "clucas//ICF/automated/2012_09_21_09_51_06" @@ -9,4 +10,7 @@ dyll17.add("dyll_HT-400ToInf_mg", '%s/%s/DYJetsToLL_HT-200To400_TuneZ2Star_8TeV-madgraph.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM/")'%(pnfs,dyChris), xs = 19.73 ) dyll17.add("dyll_M-50_mg", '%s/%s/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM/")'%(pnfs,dyChris), xs = 2950.0) - +dyll17.add("dyll_HT_10To200_M-10To50", '%s/clucas/Parked13/DYJets_10to200_M_10to50/")' % eos, xs = 13124.07) +dyll17.add("dyll_HT_10To200_M-50", '%s/clucas/Parked13/DYJets_10to200/")'% eos, xs = 3503.71) +dyll17.add("dyll_HT_200To400_M-50", '%s/clucas/Parked13/DYJets_200to400_Combined/")'% eos, xs = 24.29) +dyll17.add("dyll_HT_400ToInf_M-50", '%s/clucas/Parked13/DYJets_400toinf_Combined/")'% eos, xs = 3.3564) From d6a7ba7b0289fd05a20f7749389ff003859692b8 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Mon, 9 Dec 2013 08:47:48 -0600 Subject: [PATCH 25/34] adding eos qcd samples --- samples/__qcd17__.py | 12 ++++++++++-- 1 file changed, 10 insertions(+), 2 deletions(-) diff --git a/samples/__qcd17__.py b/samples/__qcd17__.py index a16d203..13d41ec 100644 --- a/samples/__qcd17__.py +++ b/samples/__qcd17__.py @@ -1,7 +1,8 @@ from supy.samples import SampleHolder -from supy.sites import pnfs +from supy.sites import pnfs, eos qcd17 = SampleHolder() pnfs = pnfs() +eos = eos() binsXs = [ ( 0, 4.9586188E10 ), ( 5, 4.2639499E10 ), @@ -24,7 +25,7 @@ for i,(low,xs) in enumerate(binsXs[:-1]) : hi = binsXs[i+1][0] dir = "QCD_Pt-%d%s_TuneZ2star_8TeV_pythia6.Summer12-PU_S7_START52_V9-v1.AODSIM"%(low, "to%d"%hi if hi else "") - qcd17.add("qcd_py6_pt_%d"%low, '%s/clucas//ICF/automated/2012_07_03_11_15_20/%s")'%(pnfs, dir), xs = xs) + qcd17.add("qcd_py6_pt_%d_dCache"%low, '%s/clucas//ICF/automated/2012_07_03_11_15_20/%s")'%(pnfs, dir), xs = xs) bEn = '%s/yeshaq//ICF/automated/2012_09_03_15_56_26/'%pnfs qcd17.add("qcd_b_mg_ht_250", '%s/BJets_HT-250To500_8TeV-madgraph.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM")'%bEn, xs = 5828.0) @@ -40,3 +41,10 @@ if i < 1 : dir = dir.replace("_.QCD", "QCD").replace("AODSIM_","AODSIM") if low == 1800.: dir = dir.replace("AODSIM_","AODSIM") qcd17.add("qcd_py6_pt_v2_%d"%low, '%s/agapitos//ICF/automated/2012_09_22_14_28_53/%s")'%(pnfs, dir), xs = xs) + +for i,(low,xs) in enumerate(binsXs[4:-1]) : + hi = binsXs[i+1][0] + dir = "QCD_%d%s/"%(low, "to%d"%hi if hi else "") + if low in (300,170) : dir = dir.replace("/","_merged/") + qcd17.add("qcd_py6_pt_%d"% low, '%s/clucas/Parked13/%s")'%(eos, dir), xs = xs) + From bb8090e05b9ef39521d6401804e23791e20bd432 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Mon, 9 Dec 2013 10:07:17 -0600 Subject: [PATCH 26/34] adding eos wjet samples --- samples/__ewk17__.py | 14 +++++++++++++- 1 file changed, 13 insertions(+), 1 deletion(-) diff --git a/samples/__ewk17__.py b/samples/__ewk17__.py index b04c4df..5a702d6 100644 --- a/samples/__ewk17__.py +++ b/samples/__ewk17__.py @@ -1,6 +1,7 @@ from supy.samples import SampleHolder -from supy.sites import pnfs +from supy.sites import pnfs, eos pnfs = pnfs() +eos = eos() ewk17 = SampleHolder() ### Wjets ### @@ -21,6 +22,17 @@ for part in [1,2,3,4,5] : ewk17.add("wj_lv_mg_ht_incl_v2.job673_part%i"%part, '%s/yeshaq//ICF/automated/2013_05_23_14_13_07/part%i")'%(pnfs,part), xs = {"LO":30400.0, "NLO":36257.2}["NLO"]) +for part in [1,2,3,4] : + ewk17.add("wj_lv_mg_ht_incl_v2.job50_part%i"%part, '%s/zmeng/ICF/automated/2013_09_02_15_42_12/part%i")'%(pnfs,part), xs = {"LO":30400.0, "NLO":37509}["NLO"]) + +ewk17.add("wj_lv_mg_ht_incl" , '%s/clucas/Parked13/WJets_inc/")' % eos, xs = {"NLO":37509}["NLO"]) +ewk17.add("wj_lv_mg_ht_10To150" , '%s/clucas/Parked13/WJets_10to150/")' % eos, xs = {"NLO":37509.0}["NLO"]) +ewk17.add("wj_lv_mg_ht_150To200", '%s/clucas/Parked13/WJets_150to200/")' % eos, xs = {"NLO":253.84}["NLO"]) +ewk17.add("wj_lv_mg_ht_200To250", '%s/clucas/Parked13/WJets_200to250/")' % eos, xs = {"NLO":116.46}["NLO"]) +ewk17.add("wj_lv_mg_ht_250To300", '%s/clucas/Parked13/WJets_250to300/")' % eos, xs = {"NLO":57.6358}["NLO"]) +ewk17.add("wj_lv_mg_ht_300To400", '%s/clucas/Parked13/WJets_300to400/")' % eos, xs = {"NLO":48.3664}["NLO"]) +ewk17.add("wj_lv_mg_ht_400ToInf", '%s/clucas/Parked13/WJets_400toinf/")' % eos, xs = {"NLO":30.7777}["NLO"]) + #comment 1 # CrossSection = ( NNLO_W_inclusive_from[1] / LO_W_inclusive_from_PREP[2] ) * LO_W_exclusive_from_PREP: [3] From c3458a19fb2092848905f6847581d3b6ef461322 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Tue, 10 Dec 2013 09:00:41 -0600 Subject: [PATCH 27/34] added zjet, VV, eos samples --- samples/__ewk17__.py | 9 +++++++++ 1 file changed, 9 insertions(+) diff --git a/samples/__ewk17__.py b/samples/__ewk17__.py index 5a702d6..23b0d3f 100644 --- a/samples/__ewk17__.py +++ b/samples/__ewk17__.py @@ -66,7 +66,16 @@ ewk17.add("zinv_mg_ht_200_400_ext.job500", '%s/karage//ICF/automated/2012_11_29_23_55_17/%s")' % (pnfs, ext_suffix % (200, 400)), xs={"NNLO":49.2776, "LO":41.49}["NNLO"]) ewk17.add("zinv_mg_ht_400_inf_ext.job500", '%s/karage//ICF/automated/2012_11_29_23_55_17/%s")' % (pnfs, ext_suffix % (400, "inf")), xs={"NNLO":6.2639, "LO":5.274}["NNLO"]) +ht = [(50,100),(100,200),(200,400),(400,"inf")] +xss = [{"NNLO":452.75, "LO":381.2}["NNLO"], + {"NNLO":190.39, "LO":160.3}["NNLO"], + {"NNLO":49.2776, "LO":41.49}["NNLO"], + {"NNLO":6.2639, "LO":5.274}["NNLO"]] +for ht,xs in zip(ht,xss): + ewk17.add("zinv_mg_ht_%s_%s" % ht, '%s/clucas/Parked13/ZJets_%sto%s_Combined/")' % (eos, ht[0], ht[1]), xs=xs) + diBos = "%s_TuneZ2star_8TeV_pythia6_tauola.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM/" diBosons ={"ZZ":8.25561, "WZ":32.3161, "WW":57.1097} for key in diBosons: ewk17.add("%s_pythia6.job370" % key, '%s/zmeng/ICF/automated/2012_09_21_17_43_41/%s")' % (pnfs, diBos % key), xs = diBosons[key]) + ewk17.add("%s_py6" % key, '%s/clucas/Parked13/%s")' % (eos, key), xs = diBosons[key]) From 5af194a6f5335c5c863bac14b7399450cc6af114 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Tue, 10 Dec 2013 09:24:53 -0600 Subject: [PATCH 28/34] adding top eos samples --- samples/__top17__.py | 8 +++++++- 1 file changed, 7 insertions(+), 1 deletion(-) diff --git a/samples/__top17__.py b/samples/__top17__.py index 91e5f4d..d16671d 100644 --- a/samples/__top17__.py +++ b/samples/__top17__.py @@ -1,7 +1,8 @@ from supy.samples import SampleHolder -from supy.sites import pnfs +from supy.sites import pnfs, eos import itertools pnfs = pnfs() +eos = eos() top17 = SampleHolder() top17.add("tt_8_mg.job315_1", '["/uscms/home/elaird/141_ntuples/job315/SusyCAF_Tree_100_1_h3D.root"]', xs = 234.) @@ -19,6 +20,10 @@ top17.add("ttbar_powheg_v1.job410", '%s/clucas//ICF/automated/2012_09_26_14_02_42/")'%pnfs, xs = 211.) top17.add("ttbar_powheg_v2.job404", '%s/zmeng//ICF/automated/2012_09_21_17_07_53/")'%pnfs, xs = 211.) +top17.add("ttbar_CT10_powheg", '%s/clucas/Parked13/TT_CT10_Combined/")' % eos, xs=234) +top17.add("ttz_mg", '%s/clucas/Parked13/TTZJets/")' % eos, xs=0.172) + + sTop = "-channel_TuneZ2star_8TeV-powheg-tauola.Summer12_DR53X-PU_S10_START53_V7A-v1.AODSIM/" combs = {"T_s":{"NLO":3.79, "LO":2.82},"T_t":{"NLO":56.4, "LO":47.0},"T_tW":{"NLO":11.1, "LO":10.7}, "Tbar_s":{"NLO":1.76, "LO":1.57},"Tbar_t": {"NLO":30.7, "LO":25.0},"Tbar_tW":{"NLO":11.1, "LO":10.7}} @@ -26,3 +31,4 @@ for key in combs: top17.add("%s_powheg.job368"%key, '%s/zmeng//ICF/automated/2012_09_21_17_19_15/%s")'%(pnfs, (key + sTop).replace("tW-channel","tW-channel-DR")),combs[key]["NLO"]) + top17.add("%s_powheg"%key, '%s/clucas/Parked13/%s_channel/")'%(eos, key),combs[key]["NLO"]) From 96d393d72eca0d31c014564a6dc9a227375dd60c Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Wed, 15 Jan 2014 05:28:50 -0600 Subject: [PATCH 29/34] fix --- samples/__photon17__.py | 3 ++- 1 file changed, 2 insertions(+), 1 deletion(-) diff --git a/samples/__photon17__.py b/samples/__photon17__.py index 87cc1b4..879e2be 100644 --- a/samples/__photon17__.py +++ b/samples/__photon17__.py @@ -19,8 +19,9 @@ photon17.add("g_jets_mg_ht_400_inf.job501", '%s/agapitos//ICF/automated/2012_11_29_18_15_22/%s")' % (pnfs, gJet % "400ToInf"), xs={"LO": 107.5, "NLO":124.68}["NLO"]) -for era in ["B","C","D"]: +for era in ["B","C"]: photon17.add("SinglePhotonParked.Run2012%s-22Jan2013"%era, '%s/clucas/Parked13/SinglePhoton_Run2012%s_22Jan2013/")' % (eos, era), lumi = 1.0) +photon17.add("SinglePhotonParked.Run2012D-22Jan2013", '%s/clucas/Parked13/SinglePhotonParked_Run2012D_22Jan2013/")' % eos, lumi = 1.0) photon17.add("Photon.Run2012A-22Jan2013", '%s/clucas/Parked13/Photon_Run2012A_22Jan2013/")' % eos, lumi = 1.0) photon17.add("GJets_HT200to400", '%s/clucas/Parked13/GJets_200to400/")' % eos, xs={"LO": 960.5, "NLO":1140.78}["NLO"]) From 85df3c94d219b2170d4f947c873329732d12727e Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Wed, 15 Jan 2014 06:54:11 -0600 Subject: [PATCH 30/34] newest JSON --- ...Cert_190456-208686_8TeV_22Jan2013ReReco_Collisions12_JSON.txt | 1 + 1 file changed, 1 insertion(+) create mode 100644 cert/Cert_190456-208686_8TeV_22Jan2013ReReco_Collisions12_JSON.txt diff --git a/cert/Cert_190456-208686_8TeV_22Jan2013ReReco_Collisions12_JSON.txt b/cert/Cert_190456-208686_8TeV_22Jan2013ReReco_Collisions12_JSON.txt new file mode 100644 index 0000000..913cbcd --- /dev/null +++ b/cert/Cert_190456-208686_8TeV_22Jan2013ReReco_Collisions12_JSON.txt @@ -0,0 +1 @@ +{"190645": [[10, 110]], "190646": [[1, 111]], "190659": [[33, 167]], "190679": [[1, 55]], "190688": [[69, 249]], "190702": [[51, 53], [55, 122], [124, 169]], "190703": [[1, 252]], "190704": [[1, 3]], "190705": [[1, 5], [7, 65], [81, 336], [338, 350], [353, 383]], "190706": [[1, 126]], "190707": [[1, 237], [239, 257]], "190708": [[1, 189]], "190733": [[71, 96], [99, 389], [392, 460]], "190736": [[1, 80], [83, 185]], "190738": [[1, 130], [133, 226], [229, 349]], "190782": [[55, 181], [184, 233], [236, 399], [401, 409]], "190895": [[64, 202], [210, 302], [305, 584], [587, 948]], "190906": [[73, 256], [259, 354], [356, 496]], "190945": [[124, 207]], "190949": [[1, 81]], "191043": [[45, 46]], "191046": [[1, 21], [24, 82], [84, 88], [92, 116], [119, 180], [183, 183], [185, 239]], "191056": [[1, 1], [4, 9], [16, 17], [19, 19]], "191057": [[1, 1], [4, 40]], "191062": [[1, 1], [3, 3], [5, 214], [216, 541]], "191090": [[1, 55]], "191201": [[38, 49], [52, 79]], "191202": [[1, 64], [66, 68], [87, 105], [108, 118]], "191226": [[77, 78], [81, 831], [833, 1454], [1456, 1466], [1469, 1507], [1510, 1686]], "191247": [[1, 153], [156, 280], [283, 606], [608, 620], [622, 818], [821, 834], [837, 1031], [1034, 1046], [1049, 1140], [1143, 1187], [1190, 1214], [1217, 1224]], "191248": [[1, 103]], "191264": [[59, 79], [82, 152], [155, 189]], "191271": [[56, 223], [225, 363]], "191276": [[1, 16]], "191277": [[1, 28], [30, 164], [167, 253], [255, 457], [460, 535], [537, 576], [579, 775], [778, 811], [813, 849]], "191367": [[1, 2]], "191411": [[1, 23]], "191695": [[1, 1]], "191718": [[43, 95], [98, 207]], "191720": [[1, 1], [3, 15], [17, 181]], "191721": [[1, 1], [3, 34], [36, 183], [186, 189]], "191726": [[1, 13]], "191810": [[15, 15], [22, 49], [52, 92]], "191830": [[54, 242], [245, 301], [304, 393]], "191833": [[1, 1], [3, 103]], "191834": [[1, 30], [33, 74], [77, 299], [302, 352]], "191837": [[1, 44], [47, 53], [56, 65]], "191856": [[1, 133]], "191859": [[1, 28], [31, 126]], "193093": [[1, 33]], "193123": [[1, 27]], "193124": [[1, 52]], "193192": [[58, 86]], "193193": [[1, 6], [8, 8], [11, 83], [86, 120], [122, 160], [162, 274], [276, 495], [497, 506]], "193207": [[54, 182]], "193334": [[29, 172]], "193336": [[1, 264], [267, 492], [495, 684], [687, 729], [732, 951]], "193541": [[77, 101], [103, 413], [416, 575], [578, 619]], "193556": [[41, 83]], "193557": [[1, 84]], "193575": [[48, 173], [176, 349], [351, 394], [397, 415], [417, 658], [660, 752]], "193621": [[60, 570], [573, 769], [772, 976], [979, 1053], [1056, 1137], [1139, 1193], [1195, 1371], [1373, 1654]], "193834": [[1, 35]], "193835": [[1, 20], [22, 26]], "193836": [[1, 2]], "193998": [[66, 113], [115, 278]], "193999": [[1, 45]], "194027": [[57, 113]], "194050": [[53, 113], [116, 273], [275, 355], [357, 369], [372, 391], [394, 490], [492, 814], [816, 1435], [1437, 1735], [1760, 1888]], "194051": [[1, 12]], "194052": [[1, 99], [102, 166]], "194075": [[48, 101], [103, 103], [105, 107], [109, 109], [111, 111]], "194076": [[1, 9], [11, 55], [58, 163], [165, 228], [230, 264], [267, 507], [509, 527], [530, 538], [541, 562], [565, 748]], "194108": [[81, 161], [164, 264], [266, 373], [376, 396], [398, 433], [436, 452], [454, 577], [579, 590], [593, 668], [671, 872]], "194115": [[66, 184], [186, 338], [340, 346], [348, 493], [496, 731], [819, 857]], "194117": [[1, 38]], "194119": [[1, 229], [232, 261]], "194120": [[1, 162], [165, 406]], "194150": [[42, 127], [129, 261], [264, 311]], "194151": [[47, 72], [75, 191], [193, 238], [240, 617], [619, 619], [621, 621], [623, 623]], "194153": [[1, 115]], "194199": [[96, 227], [229, 336], [339, 402]], "194210": [[3, 195], [198, 217], [220, 359], [361, 555]], "194223": [[61, 112]], "194224": [[1, 126], [129, 206], [208, 250], [253, 309], [312, 386], [389, 412]], "194225": [[1, 23], [26, 47], [49, 85], [88, 149]], "194270": [[56, 68]], "194303": [[56, 66], [69, 102]], "194304": [[1, 43], [46, 46]], "194305": [[1, 84]], "194314": [[52, 130], [133, 300]], "194315": [[1, 10], [13, 314], [317, 428], [431, 452], [455, 467]], "194317": [[1, 20]], "194424": [[63, 141], [144, 195], [198, 266], [268, 421], [424, 478], [481, 531], [534, 553], [556, 706], [708, 708]], "194428": [[1, 85], [87, 122], [125, 294], [296, 465]], "194429": [[1, 4], [7, 54], [57, 147], [150, 411], [413, 742], [745, 986], [988, 1019]], "194439": [[46, 77], [79, 106]], "194455": [[45, 64], [67, 140], [142, 255], [293, 303]], "194464": [[1, 127], [130, 142], [145, 210]], "194479": [[1, 44], [165, 232], [235, 262], [265, 374], [377, 431], [434, 489], [492, 529], [531, 566]], "194480": [[1, 32], [34, 205], [207, 375], [377, 387], [389, 759], [762, 956], [959, 1402]], "194533": [[46, 379], [382, 415], [417, 618], [620, 872]], "194619": [[31, 110]], "194631": [[1, 42], [44, 100], [102, 169], [171, 222]], "194643": [[1, 287]], "194644": [[1, 168], [171, 181], [184, 185], [187, 319], [321, 421]], "194691": [[61, 104], [107, 155], [158, 251], [254, 268], [271, 272], [275, 289], [292, 313]], "194699": [[1, 30], [32, 52], [55, 64], [67, 71], [73, 154], [157, 215], [218, 238], [241, 259]], "194702": [[1, 138], [141, 191]], "194704": [[1, 41], [44, 545], [548, 592]], "194711": [[1, 7], [9, 619]], "194712": [[1, 56], [61, 418], [420, 625], [627, 759]], "194735": [[44, 71], [74, 101], [104, 130]], "194778": [[60, 118], [120, 219]], "194789": [[1, 18], [21, 32], [34, 80], [82, 166], [168, 269], [272, 405], [409, 414], [417, 427], [430, 566]], "194790": [[1, 45]], "194825": [[72, 117], [120, 221]], "194896": [[34, 55], [58, 79], [82, 103]], "194897": [[1, 6], [8, 78], [80, 96], [98, 102]], "194912": [[53, 70], [72, 96], [98, 444], [446, 450], [453, 467], [470, 561], [564, 660], [663, 813], [815, 840], [843, 864], [866, 1004], [1007, 1025], [1027, 1067], [1069, 1137], [1140, 1166], [1168, 1249], [1251, 1304], [1307, 1444], [1447, 1487], [1489, 1503], [1506, 1662]], "194914": [[1, 38]], "194915": [[1, 74]], "195013": [[94, 144], [146, 185], [187, 206], [208, 299], [302, 324], [326, 366], [369, 447], [450, 526], [528, 541]], "195014": [[1, 6], [9, 119], [121, 148]], "195015": [[1, 13]], "195016": [[1, 21], [23, 55], [58, 63], [65, 174], [177, 184], [186, 241], [243, 246], [248, 251], [254, 367], [370, 422], [425, 560], [563, 569]], "195099": [[70, 144], [147, 186], [189, 208], [211, 224], [227, 248]], "195109": [[98, 241]], "195112": [[1, 12], [15, 26]], "195113": [[1, 209], [212, 388], [391, 403], [406, 419], [422, 492], [495, 579]], "195114": [[1, 69], [72, 103]], "195115": [[1, 7], [10, 22]], "195147": [[132, 282], [285, 294], [297, 331], [334, 363], [366, 442], [445, 536], [539, 559]], "195163": [[72, 138], [140, 224], [227, 240], [243, 243], [246, 347]], "195164": [[1, 64]], "195165": [[1, 4], [7, 41], [44, 54], [56, 153], [156, 260], [263, 266]], "195251": [[1, 131], [134, 137], [140, 152], [154, 165], [167, 242]], "195303": [[109, 191], [194, 277], [280, 310], [312, 316], [318, 409]], "195304": [[1, 3], [6, 22], [27, 80], [83, 100], [103, 154], [157, 341], [344, 588], [590, 727], [729, 1003], [1006, 1079], [1083, 1140], [1143, 1229]], "195378": [[90, 117], [120, 127], [130, 185], [187, 204], [206, 302], [305, 542], [544, 565], [567, 645], [647, 701], [703, 734], [737, 1120], [1122, 1133]], "195390": [[1, 1], [4, 27], [30, 145], [147, 183], [186, 187], [190, 208], [210, 213], [215, 400]], "195396": [[49, 55], [58, 63], [66, 131]], "195397": [[1, 10], [12, 89], [92, 120], [123, 141], [143, 251], [253, 253], [256, 475], [478, 525], [527, 608], [611, 776], [779, 970], [972, 1121], [1123, 1181], [1184, 1198], [1200, 1209]], "195398": [[3, 137], [139, 494], [497, 585], [587, 817], [820, 824], [827, 1225], [1228, 1307], [1309, 1712], [1721, 1736], [1741, 1752], [1767, 1795]], "195399": [[1, 192], [194, 382]], "195530": [[1, 80], [82, 104], [107, 156], [159, 300], [302, 405]], "195540": [[68, 123], [126, 137], [140, 283], [286, 319]], "195551": [[91, 106]], "195552": [[1, 21], [23, 27], [30, 147], [149, 155], [158, 182], [185, 287], [290, 349], [352, 469], [472, 815], [818, 823], [825, 883], [885, 1152], [1154, 1300], [1303, 1789]], "195633": [[40, 42]], "195647": [[1, 41]], "195649": [[1, 69], [72, 151], [154, 181], [183, 247]], "195655": [[1, 129], [131, 184], [186, 260], [263, 350], [353, 446], [448, 483], [485, 498]], "195656": [[1, 362]], "195658": [[1, 37], [40, 362], [364, 382], [384, 386]], "195749": [[1, 8], [10, 33], [36, 131]], "195757": [[1, 82], [85, 115], [118, 161], [163, 206]], "195758": [[1, 18]], "195774": [[1, 13], [16, 137], [139, 151], [154, 162], [164, 256], [258, 276], [279, 362], [365, 466], [469, 618], [620, 649], [651, 830]], "195775": [[1, 57], [60, 100], [103, 170]], "195776": [[1, 63], [66, 283], [286, 337], [340, 399], [401, 409], [411, 477]], "195841": [[74, 85]], "195868": [[1, 88], [90, 107], [110, 205]], "195915": [[1, 109], [111, 275], [278, 390], [393, 417], [419, 429], [432, 505], [507, 747], [749, 785], [787, 828], [830, 850]], "195916": [[1, 16], [19, 68], [71, 212]], "195917": [[1, 4]], "195918": [[1, 44], [46, 46], [49, 64]], "195919": [[1, 15]], "195923": [[1, 14]], "195925": [[1, 12]], "195926": [[1, 1], [3, 19], [21, 34]], "195929": [[1, 29]], "195930": [[1, 77], [80, 176], [179, 526], [529, 596]], "195937": [[1, 28], [31, 186], [188, 396]], "195947": [[23, 62], [64, 88]], "195948": [[51, 116], [119, 144], [147, 147], [150, 352], [355, 369], [372, 402], [404, 500], [503, 540], [543, 565], [567, 602], [605, 615]], "195950": [[1, 71], [73, 138], [141, 169], [172, 332], [335, 350], [353, 382], [385, 421], [424, 450], [453, 483], [485, 616], [619, 715], [718, 787], [789, 800], [803, 829], [831, 831], [833, 1587]], "195963": [[54, 58]], "195970": [[44, 49], [51, 85]], "196019": [[54, 68]], "196027": [[1, 55], [58, 119], [121, 155], [158, 186]], "196046": [[12, 40]], "196047": [[1, 64], [70, 75]], "196048": [[1, 44], [46, 48]], "196197": [[58, 122], [125, 179], [181, 311], [313, 516], [519, 562]], "196199": [[1, 33], [36, 83], [86, 118], [121, 147], [150, 237], [239, 285], [287, 534]], "196200": [[1, 68]], "196202": [[3, 61], [64, 108]], "196203": [[1, 102], [107, 117]], "196218": [[55, 199], [201, 224], [226, 393], [396, 494], [496, 741], [744, 752], [754, 757], [759, 820]], "196239": [[1, 59], [62, 154], [157, 272], [274, 373], [375, 432], [435, 465], [468, 647], [650, 706], [709, 1025]], "196249": [[63, 77], [80, 99]], "196250": [[1, 2], [5, 265], [267, 426]], "196252": [[1, 35]], "196334": [[59, 111], [113, 123], [126, 132], [135, 167], [170, 193], [196, 257], [259, 267], [270, 289], [292, 342]], "196349": [[65, 84], [86, 154], [157, 244], [246, 258]], "196357": [[1, 4]], "196359": [[1, 2]], "196362": [[1, 88]], "196363": [[1, 8], [11, 34]], "196364": [[1, 93], [96, 136], [139, 365], [368, 380], [382, 601], [603, 795], [798, 884], [887, 1196], [1199, 1200], [1203, 1299]], "196437": [[1, 1], [3, 74], [77, 169]], "196438": [[1, 181], [184, 699], [701, 1269]], "196452": [[82, 112], [114, 490], [493, 586], [589, 618], [622, 668], [671, 716], [718, 726], [728, 956], [958, 1004], [1007, 1091]], "196453": [[1, 74], [77, 145], [147, 669], [673, 714], [717, 799], [802, 988], [991, 1178], [1180, 1180], [1182, 1248], [1250, 1528], [1531, 1647]], "196495": [[114, 180], [182, 272]], "196509": [[1, 68]], "196531": [[62, 150], [152, 253], [256, 285], [288, 302], [305, 422], [425, 440]], "198049": [[1, 11], [14, 57]], "198050": [[2, 155]], "198063": [[1, 37], [40, 72], [74, 124], [127, 294]], "198116": [[36, 52], [54, 55], [58, 96], [98, 112]], "198207": [[1, 97]], "198208": [[1, 92], [94, 134], [137, 147], [150, 209]], "198210": [[1, 221]], "198212": [[1, 574]], "198213": [[1, 107]], "198215": [[1, 12]], "198230": [[1, 33], [36, 57], [60, 235], [237, 324], [326, 388], [390, 459], [462, 625], [627, 651], [653, 805], [808, 811], [814, 948], [950, 1090], [1093, 1103], [1106, 1332], [1335, 1380]], "198249": [[1, 7]], "198269": [[3, 198]], "198271": [[1, 91], [93, 170], [173, 299], [301, 450], [453, 513], [516, 616], [619, 628], [631, 791], [793, 797]], "198272": [[1, 185], [188, 245], [248, 314], [317, 433], [436, 444], [454, 620]], "198346": [[44, 47]], "198372": [[57, 110]], "198485": [[68, 109], [112, 134], [136, 181], [184, 239]], "198487": [[1, 145], [147, 514], [517, 668], [671, 733], [736, 757], [760, 852], [854, 994], [997, 1434], [1437, 1610]], "198522": [[65, 144], [147, 208]], "198941": [[102, 189], [191, 220], [222, 241], [243, 249], [252, 284]], "198954": [[108, 156], [159, 277]], "198955": [[1, 45], [47, 50], [53, 220], [223, 269], [271, 284], [286, 338], [340, 580], [583, 742], [744, 910], [913, 946], [949, 1162], [1165, 1169], [1172, 1182], [1185, 1188], [1190, 1246], [1249, 1304], [1306, 1467], [1470, 1485], [1487, 1552]], "198969": [[58, 81], [84, 247], [249, 323], [325, 365], [367, 413], [416, 466], [468, 643], [646, 918], [920, 1011], [1013, 1175], [1178, 1236], [1239, 1253]], "199008": [[75, 93], [95, 121], [124, 208], [211, 331], [333, 373], [376, 482], [485, 605], [608, 644]], "199011": [[1, 11], [13, 24]], "199021": [[59, 88], [91, 128], [130, 133], [136, 309], [311, 333], [335, 410], [414, 469], [471, 533], [535, 563], [565, 1223], [1226, 1479], [1481, 1494]], "199318": [[65, 138]], "199319": [[1, 7], [9, 223], [226, 277], [280, 348], [351, 358], [360, 422], [424, 490], [492, 493], [496, 612], [615, 642], [645, 720], [723, 728], [730, 731], [734, 741], [744, 752], [754, 943], [945, 997]], "199336": [[1, 33], [36, 122], [125, 231], [234, 614], [617, 789], [791, 977]], "199356": [[95, 121], [123, 168], [171, 205], [208, 231]], "199409": [[25, 54], [56, 89], [91, 204], [206, 290], [293, 583], [586, 602], [604, 1014], [1016, 1300]], "199428": [[61, 197], [200, 210], [212, 382], [387, 414], [417, 436], [439, 530], [533, 648]], "199429": [[1, 28], [30, 36], [39, 55], [58, 101], [103, 148], [151, 154]], "199435": [[63, 106], [109, 261], [263, 579], [582, 654], [656, 696], [699, 1034], [1037, 1144], [1147, 1327], [1330, 1411], [1414, 1431], [1434, 1441], [1444, 1487], [1489, 1610]], "199436": [[1, 113], [116, 254], [257, 675], [678, 748]], "199564": [[1, 3]], "199569": [[1, 2], [5, 136], [139, 367]], "199570": [[1, 17]], "199571": [[1, 184], [186, 360], [363, 561]], "199572": [[1, 317]], "199573": [[1, 22]], "199574": [[1, 53], [56, 153], [156, 246]], "199608": [[60, 157], [159, 209], [211, 341], [344, 390], [392, 461], [464, 800], [802, 1064], [1067, 1392], [1395, 1630], [1633, 1904], [1907, 1962], [1965, 2252], [2255, 2422]], "199698": [[72, 94], [96, 127]], "199699": [[1, 154], [157, 169], [172, 410], [412, 756]], "199703": [[1, 94], [97, 482], [485, 529]], "199739": [[66, 133]], "199751": [[103, 119], [121, 127]], "199752": [[1, 141], [144, 180], [182, 186], [188, 211], [214, 322]], "199753": [[1, 59]], "199754": [[1, 203], [205, 325], [328, 457], [459, 607], [610, 613], [615, 806], [808, 998]], "199804": [[78, 88], [90, 181], [183, 235], [238, 278], [281, 290], [292, 519], [522, 575], [577, 628], [631, 632]], "199812": [[70, 141], [144, 163], [182, 211], [214, 471], [474, 505], [508, 557], [560, 571], [574, 623], [626, 751], [754, 796]], "199832": [[58, 62], [65, 118], [121, 139], [142, 286]], "199833": [[1, 13], [16, 103], [105, 250], [253, 493], [496, 794], [797, 1032], [1034, 1185], [1188, 1239]], "199834": [[1, 9], [11, 11], [14, 18], [21, 54], [56, 57], [62, 65], [69, 284], [286, 503], [505, 942]], "199862": [[59, 141]], "199864": [[1, 87], [89, 89], [92, 103], [106, 372], [374, 385], [388, 486]], "199867": [[1, 134], [136, 172], [174, 218], [221, 320]], "199868": [[1, 21]], "199875": [[70, 150], [152, 334]], "199876": [[1, 19], [22, 95], [97, 249], [252, 272], [274, 340], [343, 362], [365, 376]], "199877": [[1, 173], [175, 605], [607, 701], [703, 871]], "199960": [[72, 139], [141, 197], [204, 232], [235, 363], [365, 367], [370, 380], [383, 459], [461, 466], [469, 485]], "199961": [[1, 211], [213, 287]], "199967": [[60, 120], [122, 170], [172, 198]], "199973": [[73, 89]], "200041": [[62, 83], [85, 157], [162, 274], [277, 318], [321, 335], [337, 386], [388, 389], [392, 400], [402, 568], [571, 593], [595, 646], [649, 728], [731, 860], [862, 930], [932, 1096]], "200042": [[1, 110], [112, 536]], "200049": [[1, 177]], "200075": [[76, 139], [142, 232], [256, 326], [329, 422], [425, 431], [434, 500], [502, 605]], "200091": [[67, 67], [70, 151], [154, 172], [174, 187], [190, 196], [199, 201], [204, 425], [428, 535], [537, 607], [610, 879], [881, 943], [946, 999], [1001, 1025], [1027, 1132], [1135, 1339], [1341, 1433], [1435, 1450], [1453, 1523], [1526, 1664], [1667, 1680], [1683, 1710]], "200152": [[74, 116]], "200160": [[52, 68]], "200161": [[1, 97], [100, 112]], "200174": [[81, 84]], "200177": [[1, 56]], "200178": [[1, 38]], "200180": [[1, 18]], "200186": [[1, 3], [6, 24]], "200188": [[1, 24], [27, 28], [31, 76], [79, 271], [274, 352]], "200190": [[1, 4], [6, 76], [79, 143], [146, 159], [162, 256], [258, 321], [324, 401], [403, 453], [456, 457], [460, 565], [567, 588], [591, 591], [593, 595], [597, 646], [649, 878]], "200229": [[1, 33], [41, 219], [222, 244], [247, 290], [293, 624], [627, 629]], "200243": [[69, 103], [106, 139]], "200244": [[3, 304], [307, 442], [445, 507], [510, 619]], "200245": [[1, 103], [105, 128], [131, 248], [251, 357]], "200368": [[72, 180]], "200369": [[1, 5], [8, 61], [64, 360], [363, 439], [441, 578], [580, 603], [606, 684], [686, 686]], "200381": [[8, 15], [18, 36], [38, 89], [91, 195]], "200466": [[134, 274]], "200473": [[96, 157], [159, 224], [226, 304], [306, 469], [472, 524], [527, 542], [545, 619], [622, 688], [691, 730], [733, 738], [740, 1324]], "200491": [[87, 107], [110, 149], [152, 157], [160, 197], [199, 237], [240, 270], [273, 273], [276, 334], [336, 360], [363, 419]], "200515": [[97, 183]], "200519": [[1, 111], [114, 126], [129, 136], [138, 224], [227, 258], [261, 350], [353, 611], [613, 747]], "200525": [[77, 149], [151, 164], [166, 190], [193, 276], [278, 311], [314, 464], [467, 488], [491, 674], [676, 704], [707, 755], [757, 895], [898, 937], [939, 990]], "200532": [[1, 37]], "200599": [[75, 129], [132, 137]], "200600": [[1, 183], [186, 299], [302, 313], [316, 324], [327, 334], [336, 397], [399, 417], [420, 526], [529, 591], [594, 596], [598, 609], [611, 660], [663, 823], [826, 900], [902, 943], [945, 1139]], "200961": [[1, 115]], "200976": [[94, 164]], "200990": [[75, 143]], "200991": [[1, 42], [44, 44], [47, 80], [83, 175], [178, 181], [184, 252], [255, 632], [635, 916], [918, 1017], [1019, 1048]], "200992": [[1, 405], [408, 434], [436, 581]], "201062": [[78, 268]], "201097": [[83, 136], [138, 245], [248, 300], [303, 370], [372, 429], [432, 497]], "201114": [[1, 14]], "201115": [[1, 73]], "201159": [[70, 211]], "201164": [[1, 8], [10, 94], [96, 125], [128, 178], [180, 198], [200, 271], [274, 416], [418, 418]], "201168": [[1, 37], [39, 275], [278, 481], [483, 558], [560, 730]], "201173": [[1, 194], [197, 586]], "201174": [[1, 214], [216, 263], [265, 339], [342, 451]], "201191": [[75, 98], [100, 216], [218, 389], [392, 492], [494, 506], [509, 585], [587, 594], [597, 607], [609, 794], [796, 838], [841, 974], [977, 1105], [1108, 1117], [1120, 1382], [1385, 1386]], "201193": [[1, 19]], "201196": [[1, 238], [241, 278], [286, 299], [302, 338], [341, 515], [518, 720], [723, 789], [803, 841]], "201197": [[1, 23]], "201202": [[1, 437]], "201229": [[1, 5], [8, 26], [29, 73]], "201278": [[62, 163], [166, 229], [232, 256], [259, 316], [318, 595], [598, 938], [942, 974], [976, 1160], [1163, 1304], [1306, 1793], [1796, 1802], [1805, 1906], [1909, 1929], [1932, 2174]], "201554": [[70, 86], [88, 114], [116, 126]], "201602": [[76, 81], [83, 194], [196, 494], [496, 614], [617, 635]], "201611": [[87, 145], [149, 182], [184, 186]], "201613": [[1, 42], [44, 49], [53, 210], [213, 215], [218, 225], [228, 646]], "201624": [[83, 92], [95, 240], [270, 270]], "201625": [[211, 312], [315, 348], [351, 416], [418, 588], [591, 671], [673, 758], [760, 791], [793, 944]], "201657": [[77, 93], [95, 108], [110, 118]], "201658": [[1, 19], [21, 118], [121, 136], [139, 288]], "201668": [[78, 157]], "201669": [[1, 9], [12, 136], [139, 141], [143, 165]], "201671": [[1, 120], [122, 174], [177, 462], [464, 482], [485, 499], [501, 545], [547, 571], [574, 614], [617, 766], [768, 896], [899, 911], [914, 1007]], "201678": [[1, 120]], "201679": [[1, 110], [112, 241], [244, 298], [302, 321], [324, 461], [463, 483]], "201692": [[78, 81], [83, 179]], "201705": [[65, 73], [75, 109], [111, 187]], "201706": [[1, 62]], "201707": [[1, 23], [26, 42], [45, 115], [118, 130], [133, 160], [163, 276], [279, 471], [473, 511], [514, 545], [547, 570], [572, 622], [625, 735], [738, 806], [809, 876], [879, 964]], "201708": [[1, 79]], "201718": [[58, 108]], "201727": [[67, 185]], "201729": [[6, 20], [22, 75], [77, 126], [129, 154], [156, 216], [219, 244]], "201794": [[58, 94]], "201802": [[68, 209], [211, 214], [216, 220], [223, 288], [290, 296]], "201816": [[1, 72], [74, 105], [107, 157]], "201817": [[1, 274]], "201818": [[1, 1]], "201819": [[1, 94], [96, 241]], "201824": [[1, 139], [141, 176], [179, 286], [289, 492]], "202012": [[98, 121], [126, 131]], "202013": [[1, 2], [5, 35], [38, 57]], "202014": [[1, 5], [8, 14], [16, 18], [20, 77], [79, 102], [104, 174], [177, 190], [192, 196]], "202016": [[1, 48], [51, 134], [137, 177], [179, 743], [745, 831], [834, 890], [893, 896], [898, 932], [934, 1010]], "202044": [[84, 101], [104, 266], [268, 461], [463, 466]], "202045": [[1, 30], [33, 72], [75, 528], [531, 601], [603, 785], [788, 809], [822, 823]], "202054": [[6, 266], [268, 489], [492, 605], [608, 631]], "202060": [[76, 142], [144, 154], [156, 244], [246, 497], [499, 642], [644, 682], [684, 743], [746, 936]], "202074": [[66, 174]], "202075": [[1, 18], [21, 187], [189, 214], [217, 247], [250, 342], [345, 406], [409, 497], [500, 537], [539, 539], [542, 560], [562, 615], [618, 628]], "202084": [[83, 156], [159, 177], [179, 180], [182, 239]], "202087": [[1, 25], [28, 208], [210, 357], [359, 652], [655, 853], [856, 1093]], "202088": [[1, 286]], "202093": [[1, 104], [107, 320], [322, 360]], "202116": [[59, 60]], "202178": [[67, 78], [80, 88], [91, 177], [180, 186], [188, 337], [340, 377], [379, 425], [428, 475], [478, 548], [551, 717], [720, 965], [967, 1444], [1447, 1505], [1508, 1519], [1522, 1555]], "202205": [[94, 114]], "202209": [[1, 48], [51, 142]], "202237": [[39, 128], [131, 131], [134, 219], [222, 235], [238, 275], [277, 289], [291, 316], [319, 419], [422, 538], [540, 936], [939, 950], [952, 976], [979, 1079]], "202272": [[76, 112], [115, 141], [144, 185], [188, 205], [208, 305], [307, 313], [315, 371], [436, 480], [483, 555], [558, 577], [579, 683], [686, 705], [707, 740], [742, 890], [937, 1295], [1299, 1481]], "202299": [[68, 84], [87, 141], [143, 193], [196, 358], [361, 379], [382, 414], [416, 452], [455, 555]], "202305": [[1, 89], [92, 130], [133, 323]], "202314": [[67, 104], [107, 265], [268, 278]], "202328": [[46, 89], [92, 156], [158, 276], [278, 291], [294, 434], [437, 460], [463, 586], [588, 610], [612, 614]], "202333": [[1, 235]], "202389": [[81, 182], [185, 190], [192, 199]], "202469": [[87, 158], [160, 174], [177, 352]], "202472": [[1, 96], [99, 112]], "202477": [[1, 129], [131, 150]], "202478": [[1, 177], [180, 183], [186, 219], [222, 360], [362, 506], [509, 531], [534, 718], [720, 927], [929, 973], [975, 1029], [1031, 1186], [1189, 1212], [1215, 1248]], "202504": [[77, 96], [99, 133], [135, 182], [184, 211], [213, 241], [243, 392], [395, 527], [529, 617], [620, 715], [718, 763], [766, 1172], [1174, 1247], [1250, 1471], [1474, 1679], [1682, 1704]], "202972": [[1, 30], [33, 184], [186, 290], [292, 295], [298, 371], [374, 429], [431, 544]], "202973": [[1, 234], [237, 305], [308, 437], [439, 530], [532, 541], [544, 552], [555, 851], [853, 1408]], "203002": [[77, 128], [130, 141], [144, 207], [209, 267], [270, 360], [362, 501], [504, 641], [643, 669], [671, 671], [674, 717], [720, 1034], [1037, 1070], [1073, 1370], [1372, 1392], [1395, 1410], [1413, 1596]], "203709": [[1, 121]], "203742": [[1, 29]], "203777": [[103, 113]], "203830": [[82, 182]], "203832": [[1, 11]], "203833": [[1, 70], [73, 128]], "203834": [[1, 40]], "203835": [[1, 70], [73, 358]], "203853": [[122, 222]], "203894": [[82, 272], [275, 477], [480, 902], [905, 1319]], "203909": [[79, 113], [116, 117], [120, 140], [143, 382]], "203912": [[1, 306], [308, 566], [569, 609], [611, 698], [701, 820], [823, 865], [867, 1033], [1035, 1321]], "203987": [[1, 9], [12, 241], [243, 339], [342, 781], [784, 1014]], "203992": [[1, 15]], "203994": [[1, 56], [59, 136], [139, 304], [306, 342], [344, 425]], "204100": [[117, 139]], "204101": [[1, 74]], "204113": [[82, 96], [98, 102], [105, 127], [129, 191], [194, 258], [261, 327], [329, 388], [390, 400], [402, 583], [585, 690]], "204114": [[1, 358]], "204238": [[23, 52], [55, 55]], "204250": [[92, 118], [121, 177], [179, 285], [287, 336], [339, 400], [403, 521], [524, 543], [546, 682], [684, 801]], "204511": [[1, 56]], "204541": [[5, 39], [42, 42], [44, 139], [142, 149], [151, 204]], "204544": [[1, 11], [13, 93], [96, 195], [197, 224], [226, 334], [337, 426]], "204552": [[1, 9]], "204553": [[1, 51], [53, 60], [63, 101]], "204554": [[1, 5], [7, 221], [224, 455], [458, 470], [472, 481], [483, 514]], "204555": [[1, 329], [331, 334]], "204563": [[91, 99], [102, 178], [180, 219], [222, 229], [231, 364], [366, 366], [369, 470], [473, 524], [527, 571]], "204564": [[1, 84], [87, 89], [92, 159], [161, 187], [190, 191], [193, 293], [296, 315], [317, 340], [343, 427], [429, 434], [437, 735], [737, 855], [858, 1206], [1209, 1248], [1251, 1284]], "204565": [[1, 48]], "204566": [[1, 12]], "204567": [[1, 38]], "204576": [[49, 192], [195, 301]], "204577": [[1, 46], [49, 64], [67, 105], [107, 170], [173, 181], [183, 193], [196, 653], [656, 669], [671, 740], [742, 913], [915, 1057], [1059, 1115], [1117, 1282]], "204599": [[73, 83], [85, 94], [97, 121], [124, 125], [128, 173], [175, 240], [243, 245], [248, 264], [266, 292], [294, 334]], "204601": [[1, 25], [28, 62], [65, 80], [83, 89], [92, 290], [292, 563], [565, 591], [593, 652], [655, 780], [783, 812], [814, 892], [894, 984], [986, 1003], [1006, 1038], [1040, 1088], [1091, 1102], [1105, 1161], [1164, 1250]], "205086": [[95, 149]], "205111": [[88, 390], [392, 441], [444, 446]], "205158": [[81, 289], [292, 313], [315, 473], [476, 591], [594, 595], [597, 612], [615, 663], [665, 667], [672, 685], [687, 733]], "205193": [[80, 109], [111, 349], [352, 486], [488, 650], [652, 712], [714, 902]], "205217": [[1, 12], [16, 111], [113, 171], [174, 250], [253, 318]], "205233": [[94, 153]], "205236": [[1, 190], [193, 207], [209, 260], [263, 331], [334, 352]], "205238": [[1, 6], [9, 199], [202, 254], [256, 304], [306, 355], [358, 381], [384, 596], [598, 617]], "205303": [[35, 54], [90, 132], [135, 144]], "205310": [[76, 306], [309, 313], [316, 316], [319, 321], [324, 457], [460, 559]], "205311": [[1, 85], [88, 92], [95, 183], [186, 395], [397, 592], [595, 910], [913, 1260]], "205339": [[71, 175], [178, 213], [216, 230], [233, 262], [265, 404]], "205344": [[1, 83], [86, 104], [106, 359], [362, 431], [433, 949], [951, 967], [969, 1127], [1129, 1346], [1348, 1586]], "205515": [[82, 201], [203, 216]], "205519": [[1, 47], [50, 172], [175, 367], [370, 386], [389, 472]], "205526": [[1, 269], [272, 277], [280, 332]], "205614": [[1, 4], [7, 40]], "205617": [[1, 29], [32, 102], [105, 123], [125, 140], [143, 264], [266, 448], [451, 532], [534, 547]], "205618": [[1, 12]], "205620": [[1, 175]], "205666": [[60, 119], [122, 165], [168, 259], [261, 322], [325, 578], [580, 594], [597, 721], [724, 739]], "205667": [[1, 165], [168, 282], [285, 318], [321, 412], [415, 689], [692, 751], [754, 774], [777, 1109]], "205683": [[76, 82], [85, 178], [181, 198], [201, 305]], "205690": [[1, 40]], "205694": [[1, 205], [208, 230], [233, 347], [350, 452], [455, 593], [595, 890]], "205718": [[49, 75], [78, 97], [100, 103], [105, 176], [178, 338], [341, 361], [363, 524], [527, 531], [534, 589], [591, 694]], "205774": [[1, 80]], "205777": [[1, 8]], "205781": [[1, 89], [91, 197], [200, 502]], "205826": [[80, 232], [235, 303], [306, 468]], "205833": [[84, 86], [89, 121], [123, 155], [157, 165], [167, 173], [176, 219], [221, 267], [270, 312], [315, 346], [350, 355], [360, 366]], "205834": [[1, 12], [14, 195]], "205908": [[68, 200], [202, 209]], "205921": [[22, 73], [76, 268], [271, 394], [397, 401], [410, 428], [431, 498], [500, 571], [574, 779], [782, 853]], "206066": [[89, 146]], "206088": [[86, 159], [161, 178], [181, 199], [202, 286]], "206102": [[83, 116], [120, 130], [133, 208], [211, 235], [238, 246], [249, 278], [281, 349]], "206187": [[107, 169], [172, 242], [245, 288], [290, 340], [343, 427], [429, 435], [437, 486], [489, 569], [571, 647], [649, 662], [664, 708]], "206188": [[1, 40], [42, 55]], "206199": [[1, 75], [77, 82], [85, 114]], "206207": [[82, 130], [132, 176], [179, 194], [196, 388], [390, 419], [422, 447], [450, 569], [572, 690]], "206208": [[1, 470], [472, 518]], "206210": [[11, 25], [28, 275], [277, 298], [300, 383], [386, 466]], "206243": [[62, 169], [172, 196], [199, 354], [357, 433], [435, 448], [451, 533], [536, 554], [557, 723], [726, 905]], "206245": [[1, 62]], "206246": [[1, 14], [16, 237], [240, 285], [288, 407], [412, 676], [678, 704], [706, 785], [787, 962], [965, 997], [1000, 1198], [1201, 1290]], "206257": [[1, 29]], "206258": [[1, 36], [39, 223], [226, 249]], "206302": [[1, 8], [11, 33], [36, 44], [47, 82], [84, 108], [110, 149], [151, 186], [189, 229], [231, 232], [234, 241], [243, 276]], "206303": [[1, 19], [23, 286]], "206304": [[1, 4], [6, 62]], "206331": [[91, 222], [225, 312]], "206389": [[88, 185], [187, 249], [252, 272], [275, 392]], "206391": [[1, 55], [57, 91]], "206401": [[69, 90], [92, 194], [197, 210], [212, 249], [251, 265], [267, 409]], "206446": [[92, 141], [143, 159], [162, 205], [208, 301], [304, 442], [445, 445], [448, 474], [476, 616], [619, 872], [874, 910], [912, 948], [950, 989], [992, 1030], [1033, 1075], [1109, 1149]], "206448": [[1, 143], [145, 559], [561, 1170], [1173, 1231], [1235, 1237]], "206466": [[24, 137], [140, 277], [280, 296], [299, 303], [306, 405], [407, 419], [422, 477], [480, 511], [514, 676]], "206476": [[73, 129], [133, 137], [140, 141], [143, 219]], "206477": [[1, 14], [16, 31], [33, 41], [44, 51], [53, 70], [73, 75], [77, 89], [91, 94], [97, 115], [118, 184]], "206478": [[1, 27], [29, 136], [139, 144]], "206484": [[73, 95], [98, 133], [136, 163], [166, 186], [189, 384], [387, 463], [465, 551], [554, 554], [556, 669]], "206512": [[91, 123], [125, 133], [136, 161], [163, 190], [193, 201], [203, 212], [214, 332], [334, 584], [587, 604], [607, 1005], [1008, 1123], [1126, 1163], [1165, 1211]], "206513": [[3, 39], [42, 188], [191, 234], [237, 238], [241, 323]], "206542": [[1, 115], [117, 165], [168, 511], [514, 547], [550, 603], [606, 668], [671, 727], [730, 739], [741, 833]], "206550": [[77, 132], [135, 144]], "206572": [[37, 47]], "206573": [[2, 14]], "206574": [[1, 87]], "206575": [[1, 7], [10, 10], [12, 69]], "206594": [[72, 107], [110, 246], [249, 281]], "206595": [[1, 34], [37, 42], [45, 193]], "206596": [[1, 13], [15, 220], [222, 228], [231, 236], [239, 292], [295, 695], [697, 728], [730, 810]], "206598": [[1, 81], [83, 103], [105, 588], [591, 657], [659, 719]], "206605": [[1, 36], [39, 78]], "206744": [[49, 157], [160, 192], [195, 395], [398, 452]], "206745": [[1, 81], [84, 199], [202, 224], [227, 237], [240, 304], [306, 318], [321, 720], [723, 796], [799, 894], [897, 944], [946, 1106], [1108, 1524], [1527, 1862], [1988, 1996]], "206859": [[79, 210], [212, 258], [260, 323], [325, 356], [359, 609], [612, 681], [684, 732], [734, 768], [771, 808], [811, 827], [830, 848]], "206866": [[1, 30], [33, 113], [115, 274]], "206868": [[1, 3], [10, 16]], "206869": [[1, 251], [253, 271], [274, 502], [507, 520], [522, 566], [568, 752]], "206897": [[1, 34], [38, 61], [63, 102], [109, 109], [111, 112], [114, 131], [133, 137]], "206901": [[1, 98]], "206906": [[1, 31], [38, 94], [96, 136], [138, 139], [142, 149], [151, 175], [177, 206]], "206940": [[1, 151], [153, 153], [155, 298], [301, 382], [384, 712], [715, 803], [805, 960], [963, 1027]], "207099": [[83, 134], [137, 172], [175, 213], [216, 314], [316, 320], [323, 330], [333, 367], [370, 481], [484, 602], [605, 755], [757, 1046], [1048, 1171]], "207100": [[1, 91], [94, 94]], "207214": [[57, 112], [114, 177], [179, 181], [184, 196], [199, 220], [223, 262], [265, 405], [408, 482], [485, 640], [643, 708], [718, 757], [759, 808], [811, 829]], "207217": [[1, 32]], "207219": [[1, 112]], "207220": [[1, 160]], "207221": [[1, 102]], "207222": [[1, 17], [20, 289]], "207231": [[70, 84], [86, 121], [123, 184], [187, 189], [192, 303], [306, 354], [357, 481], [484, 504], [508, 549], [552, 626], [628, 690], [693, 875], [878, 1000], [1003, 1170], [1173, 1187], [1189, 1227], [1229, 1415], [1418, 1445], [1447, 1505]], "207233": [[1, 119], [121, 148]], "207269": [[80, 394], [397, 436], [439, 463], [466, 551], [568, 577]], "207273": [[3, 877]], "207279": [[68, 138], [141, 149], [151, 237], [240, 266], [269, 307], [309, 416], [498, 551], [554, 640], [643, 961], [963, 1095], [1098, 1160]], "207320": [[1, 110], [112, 350]], "207371": [[72, 117], [120, 124]], "207372": [[1, 27], [30, 113], [116, 154], [156, 174], [176, 478], [480, 496]], "207397": [[32, 77], [80, 140], [143, 179]], "207398": [[1, 14], [16, 33]], "207454": [[79, 95], [98, 123], [126, 259], [261, 363], [365, 458], [461, 498], [501, 609], [612, 632], [635, 781], [784, 866], [869, 974], [977, 1064], [1067, 1079], [1081, 1321], [1323, 1464], [1467, 1569], [1571, 1604], [1607, 1712], [1714, 1988]], "207469": [[1, 31], [34, 45]], "207477": [[76, 104], [107, 111], [114, 147], [150, 295], [298, 483], [486, 494], [497, 527], [530, 563], [565, 570]], "207487": [[50, 98], [101, 311], [313, 359], [363, 468], [471, 472]], "207488": [[1, 63], [66, 92], [95, 113], [116, 198], [200, 250], [252, 288], [291, 365], [368, 377], [379, 440]], "207490": [[1, 48], [51, 111]], "207491": [[1, 176], [179, 458]], "207492": [[1, 20], [23, 298]], "207515": [[79, 109], [112, 132], [134, 208], [211, 225], [228, 320], [322, 381], [383, 498], [500, 730], [733, 849], [851, 954], [957, 994], [997, 1052], [1055, 1143], [1145, 1211]], "207517": [[1, 12], [15, 57]], "207518": [[1, 59], [61, 83]], "207882": [[22, 45]], "207883": [[1, 1], [3, 4], [7, 75]], "207884": [[1, 106], [108, 183]], "207885": [[1, 90]], "207886": [[1, 30], [32, 90], [92, 156], [158, 166], [168, 171]], "207889": [[1, 43], [47, 57], [60, 303], [306, 442], [445, 445], [447, 551], [553, 731], [733, 907], [910, 945]], "207898": [[1, 33], [36, 57], [60, 235], [239, 257], [260, 277]], "207905": [[75, 196], [198, 281], [284, 329], [331, 402], [404, 565], [568, 672], [675, 805], [807, 850], [852, 861], [864, 884], [886, 1180], [1183, 1283], [1285, 1331], [1333, 1515], [1518, 1734], [1737, 1796]], "207920": [[84, 146], [149, 241], [243, 261], [264, 291], [294, 486], [489, 518], [520, 598], [600, 708], [710, 826]], "207921": [[1, 37], [40, 58]], "207922": [[1, 69], [71, 100], [103, 126], [129, 242], [274, 291]], "207924": [[1, 52], [54, 171], [173, 178], [181, 339]], "208307": [[2, 42], [45, 45], [47, 70], [72, 147], [150, 252], [256, 259], [262, 275], [278, 342], [345, 450], [453, 527], [530, 583], [586, 605], [608, 616], [618, 667], [670, 761], [763, 798], [800, 889], [891, 893], [896, 1055], [1057, 1205], [1208, 1294], [1297, 1328]], "208339": [[77, 89], [91, 122], [125, 208], [211, 346], [349, 363]], "208341": [[1, 84], [87, 117], [120, 513], [515, 685], [688, 693], [695, 775], [777, 824]], "208351": [[83, 97], [100, 356], [359, 367], [369, 369]], "208352": [[1, 15], [17, 17], [19, 19]], "208353": [[1, 76], [78, 269], [271, 348]], "208357": [[1, 70], [73, 507]], "208390": [[72, 128], [130, 169]], "208391": [[52, 82], [84, 162], [164, 216], [219, 493], [495, 498], [500, 523], [526, 533], [535, 588], [591, 660], [663, 869]], "208427": [[49, 89], [92, 161], [164, 164], [166, 173], [175, 268], [271, 312], [315, 315], [317, 335], [337, 361], [364, 402], [404, 422], [425, 577], [580, 647]], "208428": [[1, 58], [61, 68], [70, 156], [159, 227]], "208429": [[1, 56], [59, 139], [141, 159], [162, 237], [240, 440], [442, 452], [455, 589], [592, 712], [715, 922]], "208487": [[2, 26], [29, 159], [161, 307], [309, 459], [462, 476], [479, 621]], "208509": [[71, 232]], "208538": [[2, 43]], "208540": [[1, 26], [29, 98]], "208541": [[1, 57], [59, 173], [175, 376], [378, 413]], "208551": [[119, 193], [195, 212], [215, 300], [303, 354], [356, 554], [557, 580]], "208686": [[73, 79], [82, 181], [183, 224], [227, 243], [246, 311], [313, 459]]} \ No newline at end of file From 3b1df3da35b8fd4bbd2a91c66c416bb650690cda Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Thu, 16 Jan 2014 06:40:14 -0600 Subject: [PATCH 31/34] postponing pu weight changes --- analyses/hadronicLook.py | 103 +++++++++++++++++---------------------- 1 file changed, 45 insertions(+), 58 deletions(-) diff --git a/analyses/hadronicLook.py b/analyses/hadronicLook.py index b1f9e27..e4ba5ce 100644 --- a/analyses/hadronicLook.py +++ b/analyses/hadronicLook.py @@ -18,14 +18,10 @@ def triggerTuple(l=[], keys=[]): return tuple(out) triggers_2012 = [{"HT": 200, "AlphaT": 0.57, "v": [5, 6, 8]}, - {"HT": 250, "AlphaT": 0.55, "v": range(1, 10)}, + {"HT": 250, "AlphaT": 0.55, "v": range(1, 9)}, {"HT": 300, "AlphaT": 0.53, "v": range(1, 4)}, - {"HT": 300, "AlphaT": 0.54, "v": range(1, 10)}, - {"HT": 350, "AlphaT": 0.52, "v": range(1, 4)}, - {"HT": 350, "AlphaT": 0.53, "v": range(1, 15)}, - {"HT": 400, "AlphaT": 0.51, "v": range(1, 15)}, - {"HT": 400, "AlphaT": 0.52, "v": range(1, 10)}, - {"HT": 450, "AlphaT": 0.51, "v": range(1, 10)}, + {"HT": 350, "AlphaT": 0.52, "v": range(1, 9)}, + {"HT": 400, "AlphaT": 0.51, "v": range(13,20)}, ] @@ -59,13 +55,13 @@ def parameters(self): #"eq3b": (3, 3), #"ge4b": (4, None), }), - "thresholds": self.vary({"200s": (200.0, 275.0, 73.3, 36.7), - #"275s": (275.0, 325.0, 73.3, 36.7), - #"325s": (325.0, 375.0, 86.7, 43.3), - "375": (375.0, None, 100.0, 50.0), - #"675": (675.0, None, 100.0, 50.0), - #"875": (875.0, None, 100.0, 50.0), - }), + "thresholds": self.vary(dict([("200s", (200.0, 275.0, 73.3, 36.7)), + ("275s", (275.0, 325.0, 73.3, 36.7)), + ("325s", (325.0, 375.0, 86.7, 43.3)), + ("375", (375.0, None, 100.0, 50.0)), + ("675", (675.0, None, 100.0, 50.0)), + ("875", (875.0, None, 100.0, 50.0)), + ][3:4])), "etRatherThanPt": True, "lowPtThreshold": 30.0, "lowPtName": "lowPt", @@ -220,8 +216,8 @@ def stepsEventCount(self, params, label=""): def stepsTrigger(self, params): return [ - #steps.trigger.lowestUnPrescaledTriggerFilter().onlyData(), - #steps.trigger.hltPrescaleHistogrammer(params["triggerList"]).onlyData(), + steps.trigger.lowestUnPrescaledTriggerFilter().onlyData(), + steps.trigger.hltPrescaleHistogrammer(params["triggerList"]).onlyData(), steps.trigger.l1Filter("L1Tech_BPTX_plus_AND_minus.v0").onlyData(), steps.trigger.physicsDeclaredFilter().onlyData(), ]+([] if params["thresholds"][1] != 275.0 else [steps.trigger.lowestUnPrescaledTriggerFilter().onlyData()]) # apply trigger in lowest HT bin @@ -342,7 +338,7 @@ def stepsPlotsTwo(self, params): _jet = params["objects"]["jet"] _met = params["objects"]["met"] return [ - supy.steps.histos.histogrammer("vertexIndices", 20, -0.5, 19.5, title=";N vertices;events / bin", funcString="lambda x:len(x)"), + supy.steps.histos.histogrammer("vertexIndices", 60, -0.5, 59.5, title=";N vertices;events / bin", funcString="lambda x:len(x)"), supy.steps.histos.histogrammer("%sIndices%s" % _jet, 20, -0.5, 19.5, title=";number of %s%s passing ID#semicolon p_{T}#semicolon #eta cuts;events / bin" % _jet, funcString="lambda x:len(x)"), steps.jet.cleanJetHtMhtHistogrammer(_jet, params["etRatherThanPt"]), supy.steps.histos.histogrammer("%sDeltaPhiStar%s%s" % (_jet[0], params["lowPtName"], _jet[1]), 20, 0.0, r.TMath.Pi(), title=";#Delta#phi*;events / bin", funcString='lambda x:x[0][0]'), @@ -353,12 +349,12 @@ def stepsPlotsTwo(self, params): def stepsHtBins(self, params): _jet = params["objects"]["jet"] - return [supy.steps.filters.value("%sSumEt%s" % _jet, min=bin) for bin in [475, 575, 675, 775, 875]] + return [supy.steps.filters.value("%sSumEt%s" % _jet, min=bin) for bin in [475, 575, 675, 775, 875, 975, 1075]] def stepsOptional(self, params): return [ - supy.steps.other.skimmer(), - #steps.other.duplicateEventCheck(), + #supy.steps.other.skimmer(), + steps.other.duplicateEventCheck(), #steps.other.pickEventSpecMaker(), #steps.other.cutBitHistogrammer(self.togglePfJet(_jet), self.togglePfMet(_met)), #steps.Print.jetPrinter(_jet), @@ -420,42 +416,39 @@ def listOfSteps(self, params): self.stepsTrigger(params) + self.stepsHtLeadingJets(params) + self.stepsXclean(params) + - self.stepsOptional(params) + + #self.stepsOptional(params) + self.stepsBtagJets(params) + self.stepsQcdRejection(params) + self.stepsPlotsOne(params) + self.stepsPlotsTwo(params) + #self.stepsMbb(params) + #self.stepsDisplayer(params) + - #self.stepsHtBins(params) + + self.stepsHtBins(params) + #self.stepsEventCount(params, label="scanAfter") + []) def listOfSampleDictionaries(self): - sh = supy.samples.SampleHolder() - sh.add("275_ge2b", '["/uscms/home/elaird/08_mbb/02_skim/2012_5fb_275_ge2b.root"]', lumi=5.0e3) - sh.add("375_ge2b", '["/uscms/home/yeshaq/nobackup/supy-output/hadronicLook/375_calo_ge2/HadronicRegion.root"]', lumi=5.0e3) - return [samples.ht17, samples.top17, samples.ewk17, samples.qcd17, sh, samples.susy17, samples.photon17] + return [samples.ht17, samples.top17, samples.ewk17, samples.qcd17, samples.susy17, samples.photon17, samples.dyll17] def listOfSamples(self, params): from supy.samples import specify def data_53X(): - jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt") + jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_22Jan2013ReReco_Collisions12_JSON.txt") out = [] - out += specify(names="HTMHTParked.Run2012B-22Jan2013-v1.job649", weights=jw2012) - out += specify(names="HTMHTParked.Run2012C-22Jan2013-v1.job649", weights=jw2012) - out += specify(names="HTMHTParked.Run2012D-22Jan2013-v1.job649", weights=jw2012) - #out += specify(names="HTMHTParked_ICF_sync_test", weights=jw2012) # , nFilesMax=1, nEventsMax=20000) + out += specify(names="HT.Run2012A-22Jan2013", weights=jw2012, nFilesMax=1, nEventsMax=10) + for era in ["B","C","D"]: + out += specify(names="HTMHTParked.Run2012%s-22Jan2013"%era, weights=jw2012, nFilesMax=1, nEventsMax=10) return out def qcd_py6(eL): - low = map(lambda x: x[0], samples.__qcd17__.binsXs)[:-1] + low = map(lambda x: x[0], samples.__qcd17__.binsXs)[4:-1] out = [] - for pt in low[low.index(80):]: - out += specify("qcd_py6_pt_v2_%d" % pt, effectiveLumi=eL) + for pt in low: + out += specify("qcd_py6_pt_%d" % pt, nFilesMax=1, nEventsMax=10) return out + def qcd_b_py6(eL): out = [] out += specify("qcd_b_py6_pt_50", effectiveLumi=eL) @@ -470,44 +463,38 @@ def g_jets_mg(eL): def w_binned(): out = [] - out += specify(names="wj_lv_mg_ht_10_150", color=r.kBlue) - out += specify(names="wj_lv_mg_ht_150_200.job663", color=r.kGreen) - out += specify(names="wj_lv_mg_ht_200_250.job672", color=r.kCyan) - out += specify(names="wj_lv_mg_ht_250_300.job498", color=r.kOrange) - out += specify(names="wj_lv_mg_ht_300_400.job498", color=r.kViolet) - out += specify(names="wj_lv_mg_ht_400_inf.job498", color=r.kAzure) + out += specify(names="wj_lv_mg_ht_10To150", color=r.kBlue, nFilesMax=1, nEventsMax=10) + out += specify(names="wj_lv_mg_ht_150To200", color=r.kBlue, nFilesMax=1, nEventsMax=10) + out += specify(names="wj_lv_mg_ht_200To250", color=r.kBlue, nFilesMax=1, nEventsMax=10) + out += specify(names="wj_lv_mg_ht_250To300", color=r.kBlue, nFilesMax=1, nEventsMax=10) + out += specify(names="wj_lv_mg_ht_300To400", color=r.kBlue, nFilesMax=1, nEventsMax=10) + out += specify(names="wj_lv_mg_ht_400ToInf", color=r.kBlue, nFilesMax=1, nEventsMax=10) return out def z_binned(): out = [] - out += specify("zinv_mg_ht_50_100.job407", color=r.kBlue) - out += specify("zinv_mg_ht_100_200.job365", color=r.kGreen) - out += specify("zinv_mg_ht_200_400.job365", color=r.kOrange) - out += specify("zinv_mg_ht_400_inf.job365", color=r.kViolet) - out += specify("zinv_mg_ht_50_100_ext.job500", color=r.kBlue) - out += specify("zinv_mg_ht_100_200_ext.job680", color=r.kGreen) - out += specify("zinv_mg_ht_200_400_ext.job500", color=r.kOrange) - out += specify("zinv_mg_ht_400_inf_ext.job500", color=r.kViolet) + out += specify("zinv_mg_ht_50_100", color=r.kBlue, nFilesMax=1, nEventsMax=10) + out += specify("zinv_mg_ht_100_200", color=r.kGreen, nFilesMax=1, nEventsMax=10) + out += specify("zinv_mg_ht_200_400", color=r.kOrange, nFilesMax=1, nEventsMax=10) + out += specify("zinv_mg_ht_400_inf", color=r.kViolet, nFilesMax=1, nEventsMax=10) return out def w_inclusive(): out = [] - for part in range(1, 6): - out += specify(names="wj_lv_mg_ht_incl_v2.job673_part%i" % part) + out += specify(names="wj_lv_mg_ht_incl", color=r.kBlue, nFilesMax=1, nEventsMax=10) return out def vv(): out = [] - #out += specify("zinv_hbb_125_powheg.job342") for diBos in ["ZZ", "WZ", "WW"]: - out += specify("%s_pythia6.job370" % diBos) + out += specify("%s_py6" % diBos, nFilesMax=1, nEventsMax=10) return out def top(): out = [] - out += specify(names="ttbar_powheg_v1.job410") + out += specify(names="ttbar_CT10_powheg", nFilesMax=1, nEventsMax=10) for sTop in ["T_s", "T_t", "T_tW", "Tbar_s", "Tbar_t", "Tbar_tW"]: - out += specify("%s_powheg.job368" % sTop) + out += specify("%s_powheg" % sTop, nFilesMax=1, nEventsMax=10) return out def susy(eL): @@ -524,8 +511,8 @@ def sms(): return out return ( - data_53X() + - #w_binned() + + #data_53X() + + w_binned() + #z_binned() + #top() + #vv() + @@ -573,10 +560,10 @@ def md(x, y): def conclude(self, conf): org = self.organizer(conf) - self.mergeSamples(org) + #self.mergeSamples(org) ##for skimming only #utils.printSkimResults(org) - org.scale() if not self.parameters()["signalScan"] else org.scale(100.0) + org.scale(toPdf=True) if not self.parameters()["signalScan"] else org.scale(100.0) self.makeStandardPlots(org) #self.makeIndividualPlots(org) #self.makeEfficiencyPlots(org) From cf4d8198fca99376d05f125e70e9a207e495f35e Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Fri, 24 Jan 2014 06:32:13 -0600 Subject: [PATCH 32/34] updating lumis --- samples/__muon17__.py | 8 ++++---- 1 file changed, 4 insertions(+), 4 deletions(-) diff --git a/samples/__muon17__.py b/samples/__muon17__.py index 96a42c3..5b20d80 100644 --- a/samples/__muon17__.py +++ b/samples/__muon17__.py @@ -27,7 +27,7 @@ muon17.add("SingleMu.Run2012D-PromptReco-v1.AOD.job508", '%s/yeshaq/ICF/automated/2012_12_04_16_51_31/SingleMu.Run2012D-PromptReco-v1.AOD/")'%pnfs, lumi = 1 ) #/pb muon17.add("SingleMu.Run2012D-PromptReco-v1.AOD.job525", '%s/yeshaq/ICF/automated/2013_01_18_14_00_11/")'%pnfs, lumi = 1 ) #/pb -muon17.add("SingleMu.Run2012A-22Jan2013", '%s/clucas/Parked13/SingleMu_Run2012A_22Jan2013/")'%eos, lumi = 1) #/pb -muon17.add("SingleMu.Run2012B-22Jan2013", '%s/clucas/Parked13/SingleMu_Run2012B_22Jan2013/")'%eos, lumi = 1) #/pb -muon17.add("SingleMu.Run2012C-22Jan2013", '%s/clucas/Parked13/SingleMu_Run2012C_22Jan2013/")'%eos, lumi = 1) #/pb -muon17.add("SingleMu.Run2012D-22Jan2013", '%s/clucas/Parked13/SingleMu_Run2012D_22Jan2013/")'%eos, lumi = 1) #/pb +eras = {"A":887.6540,"B":4419.0,"C":7119.0,"D":7295.0} +for era in eras: + muon17.add("SingleMu.Run2012%s-22Jan2013" % era, '%s/clucas/Parked13/SingleMu_Run2012%s_22Jan2013/")'% (eos, era), lumi = eras[era]) #/pb + From ba2e8258f75d094bc77b23674f2bd222de336986 Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Fri, 24 Jan 2014 06:33:24 -0600 Subject: [PATCH 33/34] updating lumis --- samples/__photon17__.py | 13 +++++++------ 1 file changed, 7 insertions(+), 6 deletions(-) diff --git a/samples/__photon17__.py b/samples/__photon17__.py index 879e2be..a2f2648 100644 --- a/samples/__photon17__.py +++ b/samples/__photon17__.py @@ -19,10 +19,11 @@ photon17.add("g_jets_mg_ht_400_inf.job501", '%s/agapitos//ICF/automated/2012_11_29_18_15_22/%s")' % (pnfs, gJet % "400ToInf"), xs={"LO": 107.5, "NLO":124.68}["NLO"]) -for era in ["B","C"]: - photon17.add("SinglePhotonParked.Run2012%s-22Jan2013"%era, '%s/clucas/Parked13/SinglePhoton_Run2012%s_22Jan2013/")' % (eos, era), lumi = 1.0) -photon17.add("SinglePhotonParked.Run2012D-22Jan2013", '%s/clucas/Parked13/SinglePhotonParked_Run2012D_22Jan2013/")' % eos, lumi = 1.0) +eras = {"B":4419.0,"C":7114.0} +for era in eras: + photon17.add("SinglePhotonParked.Run2012%s-22Jan2013"%era, '%s/clucas/Parked13/SinglePhoton_Run2012%s_22Jan2013/")' % (eos, era), lumi = eras[era]) +photon17.add("SinglePhotonParked.Run2012D-22Jan2013", '%s/clucas/Parked13/SinglePhotonParked_Run2012D_22Jan2013/")' % eos, lumi = 7289.0) -photon17.add("Photon.Run2012A-22Jan2013", '%s/clucas/Parked13/Photon_Run2012A_22Jan2013/")' % eos, lumi = 1.0) -photon17.add("GJets_HT200to400", '%s/clucas/Parked13/GJets_200to400/")' % eos, xs={"LO": 960.5, "NLO":1140.78}["NLO"]) -photon17.add("GJets_HT400toinf", '%s/clucas/Parked13/GJets_400toinf/")' % eos, xs={"LO": 960.5, "NLO":1140.78}["NLO"]) +photon17.add("Photon.Run2012A-22Jan2013", '%s/clucas/Parked13/Photon_Run2012A_22Jan2013/")' % eos, lumi = 886.5740) +photon17.add("GJets_HT200to400", '%s/clucas/Parked13/GJets_200to400/")' % eos, xs={"LO": 960.5, "NLO":1140.78}["LO"]) +photon17.add("GJets_HT400toinf", '%s/clucas/Parked13/GJets_400toinf/")' % eos, xs={"LO": 107.5, "NLO":124.68}["NLO"]) From bb15479f6f786a03c01752191ac31b5466bd17de Mon Sep 17 00:00:00 2001 From: Yossof Eshaq Date: Fri, 24 Jan 2014 06:46:39 -0600 Subject: [PATCH 34/34] updating jsonMaker --- analyses/jsonMaker.py | 58 ++++++++++++------------------------------- 1 file changed, 16 insertions(+), 42 deletions(-) diff --git a/analyses/jsonMaker.py b/analyses/jsonMaker.py index 71e62ab..9fb3362 100644 --- a/analyses/jsonMaker.py +++ b/analyses/jsonMaker.py @@ -2,52 +2,26 @@ class jsonMaker(supy.analysis) : def parameters(self) : - jwPrompt = calculables.other.jsonWeight("cert/Cert_160404-178677_7TeV_PromptReco_Collisions11_JSON.sub.txt") - jwMay = calculables.other.jsonWeight("cert/Cert_160404-163869_7TeV_May10ReReco_Collisions11_JSON_v3.txt") - jwAug = calculables.other.jsonWeight("cert/Cert_170249-172619_7TeV_ReReco5Aug_Collisions11_JSON_v3.txt") - jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_PromptReco_Collisions12_JSON.txt") - + jw2012 = calculables.other.jsonWeight("cert/Cert_190456-208686_8TeV_22Jan2013ReReco_Collisions12_JSON.txt") group = self.vary() -# group['SingleMu'] = [(["SingleMu.Run2012A-13Jul2012-v1.AOD.job358", -# "SingleMu.Run2012A-13Jul2012-v1.AOD.job467", -# -# "SingleMu.Run2012A-recover-06Aug2012-v1.AOD.job359", -# "SingleMu.Run2012A-recover-06Aug2012-v1.AOD.job477", -# -# "SingleMu.Run2012B-13Jul2012-v1.AOD.job358", -# "SingleMu.Run2012B-13Jul2012-v1.AOD.job461", -# "SingleMu.Run2012C-24Aug2012-v1.AOD.job361", -# "SingleMu.Run2012C-24Aug2012-v1.AOD.job470", -# -# "SingleMu.Run2012C-PromptReco-v2.AOD.job360", -# "SingleMu.Run2012C-PromptReco-v2.AOD.job747", -# "SingleMu.Run2012D-PromptReco-v1.AOD.job508", -# "SingleMu.Run2012D-PromptReco-v1.AOD.job525" -# ], jw2012_v2)] -# -# group['SingleEl'] = [(['SingleElectron.Run2012A-PromptReco-v1.AOD.job229', -# 'SingleElectron.Run2012B-PromptReco-v1.AOD.job228', -# 'SingleElectron.Run2012B-PromptReco-v1.AOD.job238', -# ], jw2012)] -# -# group['Photon'] = [(["Photon.Run2012A-PromptReco-v1.AOD.job228", -# "SinglePhoton.Run2012B-PromptReco-v1.AOD.job229", -# "SinglePhoton.Run2012B-PromptReco-v1.AOD.job238", -# ], jw2012)] -# -# group['Mumu'] = [(['DoubleMu.Run2012A-PromptReco-v1.AOD.job229', -# 'DoubleMu.Run2012B-PromptReco-v1.AOD.job228', -# 'DoubleMu.Run2012B-PromptReco-v1.AOD.job239', -# ], jw2012)] -# - group['HT1'] = [(["HTMHTParked.Run2012B-22Jan2013-v1.job649", - "HTMHTParked.Run2012C-22Jan2013-v1.job649", - "HTMHTParked.Run2012D-22Jan2013-v1.job649"], jw2012)] - + group['SingleMu'] = [(["SingleMu.Run2012%s-22Jan2013" + % era for era in ["A","B","C","D"]], + jw2012)] - return {'group':group} + group['Photon']=[(["Photon.Run2012A-22Jan2013", + "SinglePhotonParked.Run2012B-22Jan2013", + "SinglePhotonParked.Run2012C-22Jan2013", + "SinglePhotonParked.Run2012D-22Jan2013"], + jw2012)] + + group['HT'] = [(["HT.Run2012A-22Jan2013"] + + ["HTMHTParked.Run2012%s-22Jan2013" + % era for era in ["B","C","D"]], + jw2012)] + return {'group':group} + def listOfSteps(self,pars) : return [ supy.steps.printer.progressPrinter(2,300), steps.other.jsonMaker(pixelLumi = False),