#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, ))