#Plot("nJet", ["TH1D","nJet","",5,-0.5,4.5], LambdaFunc('x: x.njets_pt30_eta2p5[0]'), ),
]
#____________________________________________________________________________________________________
#
plots = muon_plots + general_plots + jet_plots + signif_plots# big list of Plot objects
#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 = darkphoton_signal_unblind_sequence
sequence.add(plotter)
outputInfo = OutputInfo("OutputInfo")
outputInfo.outputDir = outputDir
outputInfo.TFileName = "DataMCDistribution_radiuscut"+radcutstr+"GeV.root"
endSequence = EndSequence(skipHadd=justEndSequence)
endModuleOutputDir = "/home/"+User+"/public_html/Higgs/HToZdZd/"+out_path
endSequence.add(PlotEndModule(endModuleOutputDir,plots,skipSF=True))
## Put an index.php file into Plots dir for easy visualization
copyFile('/home/'+User+'/',phpFile,endModuleOutputDir)
for component in dataset.componentList:
component.maxEvents = nEvents
plots = []
plotter = Plotter("Plotter", plots)
leptonJetProducer = LeptonJetProducer("LeptonJetProducer", "Run2016")
xsWeighter = XSWeighter("XSWeighter")
NJet40Producer = NJet40Producer("NJet40Producer")
CategoryProducer = CategoryProducer("CategoryProducer")
YieldCounter = YieldCounter("YieldCounter")
#baselineSkimmer = BaselineSkimmer("BaselineSkimmer")
#signalRegionSkimmer = SignalRegionSkimmer("SignalRegionSkimmer")
sequence = Sequence()
#sequence.add(leptonJetProducer)
sequence.add(xsWeighter)
#sequence.add(baselineSkimmer)
sequence.add(NJet40Producer)
sequence.add(CategoryProducer)
sequence.add(YieldCounter)
sequence.add(plotter)
endSequence = EndSequence(skipHadd=False)
endModuleOutputDir = "/home/kshi/public_html/RA5/yieldcount/"
endSequence.add(PlotEndModule(endModuleOutputDir, plots))
outputInfo = OutputInfo("OutputInfo")
outputInfo.outputDir = outputDir
outputInfo.TFileName = "MCDistributions.root"
anaProducer = AnalysisProducer("AnaProducer")
eventSkimmer = EventSkimmer("StopToBLepSkim", cutflow="stCR")
hltSkimmer = HLTSkimmer("HLTSkim", cutflow="stCR")
jsonSkimmer = JSONSkimmer("JSONSkim")
sequence = Sequence()
sequence.add(jsonSkimmer)
sequence.add(hltSkimmer)
sequence.add(mediumMuonProducer)
sequence.add(mediumElectronProducer)
sequence.add(jetProducer)
sequence.add(xsWeighter)
sequence.add(puWeighter)
sequence.add(anaProducer)
sequence.add(eventSkimmer)
sequence.add(plotter)
endSequence = EndSequence(skipHadd=False)
if not ratio_switch:
endModuleOutputDir = "/home/kshi/public_html/dataPlot/stCR/"
endSequence.add(PlotEndModule(endModuleOutputDir, data_plots,
ratio_switch))
else:
endModuleOutputDir = "/home/kshi/public_html/ratioPlot/stCR/"
endSequence.add(
PlotEndModule(endModuleOutputDir, ratio_plots, ratio_switch))
outputInfo = OutputInfo("OutputInfo")
outputInfo.outputDir = outputDir
outputInfo.TFileName = "DataMCDistribution.root"
Plot("LepTightPt1" + lepCat,
["TH1D", "LepTightPt1" + lepCat, "", 20, 0., 200.],
LambdaFunc('x: x.firstLep.pt')),
Plot("LepTightEta1" + lepCat,
["TH1D", "LepTightEta1" + lepCat, "", 10, -3., 3.],
LambdaFunc('x: x.firstLep.eta')),
Plot("LepTightPhi1" + lepCat,
["TH1D", "LepTightPhi1" + lepCat, "", 10, -5., 5.],
LambdaFunc('x: x.firstLep.phi')),
#Plot("JetPt1"+lepCat, ["TH1D","JetPt1"+lepCat,"",20,0.,400.], LambdaFunc('x: x.selJets[0].pt') ),
#Plot("JetEta1"+lepCat, ["TH1D","JetEta1"+lepCat,"",10,-3.,3.], LambdaFunc('x: x.selJets[0].eta') ),
#Plot("JetPhi1"+lepCat, ["TH1D","JetPhi1"+lepCat,"",10,-5.,5.], LambdaFunc('x: x.selJets[0].phi') ),
#Plot("JetPt2"+lepCat, ["TH1D","JetPt2"+lepCat,"",20,0.,400.], LambdaFunc('x: x.selJets[1].pt') ),
#Plot("JetEta2"+lepCat, ["TH1D","JetEta2"+lepCat,"",10,-3.,3.], LambdaFunc('x: x.selJets[1].eta') ),
#Plot("JetPhi2"+lepCat, ["TH1D","JetPhi2"+lepCat,"",10,-5.,5.], LambdaFunc('x: x.selJets[1].phi') ),
])
plotter = Plotter("Plotter", plots)
sequence = gamma_cr_sequence
sequence.add(plotter)
endSequence = EndSequence(skipHadd=False, )
endSequence.add(PlotEndModule(endModuleOutputDir, plots, scaleToData=True))
outputInfo = OutputInfo("OutputInfo")
outputInfo.outputDir = outputDir
outputInfo.TFileName = "DataMCDistributions.root"
"mZ2_2mu2e",
"mZ12_4mu",
"mZ12_4e",
"mZ12_2e2mu",
"mZ12_2mu2e",
]:
p.customHistDict["ZPlusX"] = BaseObject(
p.key,
hist=copy.deepcopy(inputShapeFile.Get(p.key + "_shapehist")))
#p.customPdfDict["ZPlusX"] = BaseObject(p.key,hist=copy.deepcopy(inputShapeFile.Get(p.key+"_shapehist")))
#p.customPdfDict = {}
#leptonChannel = p.key.split("_")[-1]
#p.customPdfDict["ZPlusX"] = BaseObject(p.key,hist=inputShapeFile.Get("mZ2"+"_"+leptonChannel+"_shapehist").Clone(p.key))
for sig in sigSamples:
for p in plots:
p.plotSetting.line_style_dict[sig.name] = 10
p.plotSetting.line_width_dict[sig.name] = 4
outputInfo = OutputInfo("OutputInfo")
outputInfo.outputDir = outputDir
outputInfo.TFileName = "DataMCDistribution.root"
endSequence = EndSequence(skipHadd=True)
endModuleOutputDir = system.getPublicHtmlPath() + "/Higgs/HToZdZd/" + out_path
endSequence.add(PlotEndModule(
endModuleOutputDir,
plots,
skipSF=False,
))
