"ggZZ": 0.4, "Data": 1.0, } p.plotSetting.scatter_density = "1.0" p.selectedSamples = [ "Higgs", "ZZ", "Data", ] #p.plotSetting.cms_lumi = True #p.plotSetting.tdr_style = True #p.plotSetting.SetNColumns = 4 else: p.plotSetting.divideByBinWidth = True plotter = Plotter("Plotter", plots) sequence = darkphoton_signal_unblind_sequence sequence.add(plotter) outputInfo = OutputInfo("OutputInfo") outputInfo.outputDir = outputDir outputInfo.TFileName = "DataMCDistribution.root" endSequence = EndSequence( skipHadd=skipHadd, haddDataSamples=True, ) endModuleOutputDir = system.getPublicHtmlPath() + "/Higgs/HToZdZd/" + out_path endSequence.add(PlotEndModule(endModuleOutputDir, plots, skipSF=False))
makeTFile = True makePlots = True fitlist = [ "pol5", "pol6", "pol7", "pol8", "pol9", "landau", ] in_path = "DarkPhotonSR/StatInput/2020-02-29_SignalInterpolation2D_Run2017/" inputDir = system.getStoragePath() + "/" + User + "/Higgs/HToZdZd/" + in_path TFileName = "StatInput.root" outputDir = system.getPublicHtmlPath( ) + "/Higgs/HToZdZd/Interpolation/" + os.path.basename( os.path.normpath(in_path)) y_range = [0., 0.5] draw_option = "AP" # ________________________________________________________________________________________________ || signals = [SignalModel("HToZdZd_M" + str(m), m) for m in mass_points] bins = [ BaseObject("MuMu", histName="MuMu", x_width=0.02, y_width=0.02), BaseObject("ElMu", histName="ElMu", x_width=0.05, y_width=0.02), BaseObject("ElEl", histName="ElEl", x_width=0.05, y_width=0.05), BaseObject("MuEl", histName="MuEl", x_width=0.02, y_width=0.05), ]
for sig in ppZZdSamples: for p in plots: p.plotSetting.line_style_dict[sig.name] = 9 p.plotSetting.line_width_dict[sig.name] = 4 p.plotSetting.line_color_dict[sig.name] = ROOT.kOrange #for plot in plots: # plot.plotSetting.divideByBinWidth = True for dataset in componentList: if dataset.isMC: dataset.lumi = lumi for component in dataset.componentList: component.maxEvents = nEvents plotter = Plotter("Plotter", plots) variableProducer = VariableProducer("VariableProducer") sequence = darkphoton_signal_sequence #sequence = higgs_m4lNarrowWindow_sequence sequence.add(variableProducer) sequence.add(plotter) outputInfo = OutputInfo("OutputInfo") outputInfo.outputDir = outputDir outputInfo.TFileName = "DataMCDistribution.root" endSequence = EndSequence(skipHadd=justEndSequence) endModuleOutputDir = system.getPublicHtmlPath() + "/Higgs/DarkZ/" + out_path endSequence.add(PlotEndModule(endModuleOutputDir, plots, skipSF=False))
# )), #Plot("SummaryCategory", ["TH1D","SummaryCategory","",10,0.5,10.5], LambdaFunc('x: x.cateogryNumber'), selFunc=LambdaFunc('x: x.region == \"CR\"')), ] for dataset in componentList: if dataset.isMC: dataset.lumi = lumi for component in dataset.componentList: component.maxEvents = nEvents plotter = Plotter("Plotter", plots) xsWeighter = XSWeighter("XSWeighter") dataMCWeighter = DataMCWeighter("DataMCWeighter") anaSkimmer = AnalysisSkimmer("AnalysisSkimmer") catProducer = CategoryProducer("CategoryProducer") sequence = Sequence() sequence.add(anaSkimmer) sequence.add(xsWeighter) sequence.add(dataMCWeighter) sequence.add(catProducer) sequence.add(plotter) outputInfo = OutputInfo("OutputInfo") outputInfo.outputDir = outputDir outputInfo.TFileName = "DataMCDistribution.root" endSequence = EndSequence(skipHadd=False, ) endModuleOutputDir = system.getPublicHtmlPath() + "/LJMet/B2G/" + out_path endSequence.add(PlotEndModule(endModuleOutputDir, plots))
] justEndSequence = False plots = lepton_plots + photon_plots + general_plots for p in plots: p.plotSetting.line_color_dict[ALP_HToZa_M15.name] = ROOT.kRed p.plotSetting.line_color_dict[ALP_HToZa_M30.name] = ROOT.kBlue #for plot in plots: # plot.plotSetting.divideByBinWidth = True for dataset in componentList: if dataset.isMC: dataset.lumi = lumi for component in dataset.componentList: component.maxEvents = nEvents plotter = Plotter("Plotter", plots) sequence = sr_resolved_sequence sequence.add(plotter) outputInfo = OutputInfo("OutputInfo") outputInfo.outputDir = outputDir outputInfo.TFileName = "DataMCDistribution.root" endSequence = EndSequence(skipHadd=justEndSequence) endModuleOutputDir = system.getPublicHtmlPath() + "/HToZaToLLGG/" + out_path endSequence.add(PlotEndModule(endModuleOutputDir, plots, skipSF=False))
for sig in sigSamples: for p in plots: p.plotSetting.line_style_dict[sig.name] = 10 p.plotSetting.line_width_dict[sig.name] = 4 p.plotSetting.line_color_dict[sig.name] = ROOT.kRed for sig in ppZZdSamples: for p in plots: p.plotSetting.line_style_dict[sig.name] = 9 p.plotSetting.line_width_dict[sig.name] = 4 p.plotSetting.line_color_dict[sig.name] = ROOT.kOrange for s in sample2016.bkgSamples: s.input = BaseObject("2016",inputDir=os.path.join(baseDir,"2019-11-21_Run2016/"),postfix="_Run2016",) for s in sample2017.bkgSamples: s.input = BaseObject("2017",inputDir=os.path.join(baseDir,"2019-11-21_Run2017/"),postfix="_Run2017",) for s in sample2018.bkgSamples: s.input = BaseObject("2018",inputDir=os.path.join(baseDir,"2019-11-21_Run2018/"),postfix="_Run2018",) for s in sample2016.dataSamples: s.input = BaseObject("2016",inputDir=os.path.join(baseDir,"2019-11-21_Run2016/"),postfix="",) for s in sample2017.dataSamples: s.input = BaseObject("2017",inputDir=os.path.join(baseDir,"2019-11-21_Run2017/"),postfix="",) for s in sample2018.dataSamples: s.input = BaseObject("2018",inputDir=os.path.join(baseDir,"2019-11-21_Run2018/"),postfix="",) for cmp in componentList: targetDir = os.path.join(outputDir,cmp.name) sourceDir = os.path.join(cmp.input.inputDir,cmp.name.replace(cmp.input.postfix,"")) shutil.copytree(sourceDir,targetDir) outputInfo = OutputInfo("OutputInfo") outputInfo.outputDir = outputDir outputInfo.TFileName = "DataMCDistribution.root" endSequence = EndSequence(skipHadd=False) endModuleOutputDir = system.getPublicHtmlPath()+"/Higgs/DarkZ/DarkPhotonSR/ShapeTemplate/"+outputBaseDir endSequence.add(PlotEndModule(endModuleOutputDir,plots,skipSF=False))