def go(): directory = "/afs/cern.ch/user/e/elaird/work/public/d1_utca/" uscFiles = rootFiles(directory, mode="usc") if options.onlysummary: castorFiles = [] else: castorFiles = rootFiles(directory, mode="castor") labels = [] for run in sorted(uscFiles.keys()): if any([options.min and run < int(options.min), options.max and run > int(options.max), options.run and run != int(options.run), #(not options.onlyutca) and (run not in castorFiles), ]): continue label = "Run%d" % run labels.append(label) if options.onlysummary: continue if options.onlyutca: cmsFileName = "" else: #cmsFileName = castorFiles[run] cmsFileName = eosFile(run) analyze.oneRun(utcaFileName="" if options.onlycms else uscFiles[run], cmsFileName=cmsFileName, label=label, useEvn=False, filterEvn=options.filterevn, ornTolerance=1) if options.onlysummary and options.run: labels = ["Run"+options.run] graphs.makeSummaryPdf(labels)
def plot(prefix, roots, feds1s, feds2s): graphs.makeSummaryPdfMulti(inputFiles=roots, feds1s=feds1s, feds2s=feds2s, pdf="output/%s.pdf" % prefix) if not any(feds2s): outputFile = "output/%s.root" % prefix cmd = " ".join(["hadd -f", outputFile] + roots) os.system(cmd) graphs.makeSummaryPdf( inputFile=outputFile, feds1=sorted(sum(feds1s, [])), pdf=outputFile.replace(".root", "_merged.pdf") )
def plot(prefix, roots, feds1s, feds2s): graphs.makeSummaryPdfMulti(inputFiles=roots, feds1s=feds1s, feds2s=feds2s, pdf="output/%s.pdf" % prefix) if not any(feds2s): outputFile = "output/%s.root" % prefix cmd = " ".join(["hadd -f", outputFile] + roots) os.system(cmd) graphs.makeSummaryPdf(inputFile=outputFile, feds1=sorted(sum(feds1s, [])), pdf=outputFile.replace(".root", "_merged.pdf"))
uhtr=True, printEventMap=False, nEvents=nEvents, ) if run <= 200000: # including None if options.file: fileName = options.file else: fileName = baseDir + "/904/B904_Integration_%06d.root" % run analyze.oneRun( utcaFileName=fileName if not options.patterns else "", utcaFedIds=[931], cmsFileName=fileName, cmsFedIds=[702], cmsIsLocal=True, cmsPatternMode=options.patterns, label=label, useEvn=False, filterEvn=False, ornTolerance=1, uhtr=True, printEventMap=False, nEvents=nEvents, ) if not options.patterns: import graphs analyze.printHisto(label) graphs.makeSummaryPdf(labels=[label])