"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))
