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