#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)
nCores = 8
outputDir = "/raid/raid7/lucien/" + out_path
nEvents = -1
disableProgressBar = False
justEndSequence = True
componentList = sigSamples
for dataset in componentList:
if dataset.isMC:
dataset.lumi = 77.30
for component in dataset.componentList:
component.maxEvents = nEvents
anaSkimmer = AnalysisSkimmer("AnalysisSkimmer")
genWeightCounter = GenWeightCounter("GenWeightCounter")
sequence = Sequence()
sequence.add(anaSkimmer)
sequence.add(genWeightCounter)
efficiencyEndModule = EfficiencyEndModule("/home/lucien/public_html/" +
out_path)
endSequence = EndSequence(skipHadd=justEndSequence)
endSequence.add(efficiencyEndModule)
outputInfo = OutputInfo("OutputInfo")
outputInfo.outputDir = outputDir
outputInfo.TFileName = "EventWeight.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"
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"
fidskimmer = FiducialSkimmer("FiducialSkimmer")
fidCounter = GenWeightCounter("GenWeightCounter", postfix="FidicialCut")
plotter = Plotter("Plotter", validation_plots)
sequence = Sequence()
sequence.add(xsWeighter)
sequence.add(preskimmer)
sequence.add(preskimCounter)
sequence.add(fidskimmer)
sequence.add(fidCounter)
#sequence.add(plotter)
cutflows = [
"Preskim",
"FidicialCut",
]
endModuleOutputDir = "/home/lucien/public_html/" + out_path
cutflowEndModule = CutflowEndModule(endModuleOutputDir,
cutflows=cutflows,
ignoreSumw=True)
endSequence = EndSequence(skipHadd=justEndSequence)
endSequence.add(cutflowEndModule)
#endSequence.add(PlotEndModule(endModuleOutputDir,validation_plots))
outputInfo = OutputInfo("OutputInfo")
outputInfo.outputDir = outputDir
outputInfo.TFileName = "EventWeight.root"
fname = f.replace(".root", "")
tmp = Dataset(
fname,
ComponentList(
[Component(fname, os.path.join(inputDir, f), "tree", False)]),
)
componentList.append(tmp)
rootTree = "tree"
nCores = 10
disableProgressBar = False
nEvents = -1
outputDir = "output/GasGain_2021-03-04/"
justEndSequence = False
sequence = Sequence()
sequence.add(SkimTreeZMuMuSkimmer("ZMuMuSkimmer"))
gasGainPlotter = SkimTreeGasGainPlotter("GasGainPlotter")
skimTreeRunPlotter = SkimTreeRunPlotter("RunPlotter")
sequence.add(gasGainPlotter)
sequence.add(skimTreeRunPlotter)
endSequence = EndSequence(skipHadd=justEndSequence)
endModuleOutputDir = outputDir
skimTreeGasGainEndModule = SkimTreeGasGainEndModule(endModuleOutputDir)
endSequence.add(skimTreeGasGainEndModule)
outputInfo = OutputInfo("OutputInfo")
outputInfo.outputDir = outputDir
outputInfo.TFileName = "test1.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"
elif where == "ihepa":
out_path = "TightLoose/BTagEfficiency/2018-10-16/"
outputDir = "/raid/raid7/lucien/SUSY/RA5/"+out_path
endModuleOutputDir = "/home/lucien/public_html/SUSY/RA5/"+out_path
lepCats = [""]
nCores = 5
nEvents = -1
disableProgressBar = False
justEndSequence = False
verbose = False
componentList = [ c for c in skimComponentDict.values() if c.isMC]
for dataset in componentList:
if dataset.isMC:
dataset.lumi = 35.9
for component in dataset.componentList:
component.maxEvents = nEvents
btagEffAna = BTagEffAnalyzer("BTagEffAnalyzer")
btagEffEnd = BTagEffEndModule(endModuleOutputDir,cutflow="TightLoose",)
sequence = tl_rpv_sequence
sequence.add(btagEffAna)
endSequence = EndSequence(skipHadd=False,haddAllSamples=True)
endSequence.add(btagEffEnd)
outputInfo = OutputInfo("OutputInfo")
outputInfo.outputDir = outputDir
outputInfo.TFileName = "BTagEfficiency.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,
))
