subStep = 2 * deltaSpot / nSteps
regionBins[iHisto] = []
for ix in range(nSteps):
iPhi = phiSpot - deltaSpot + ix * subStep
tmpBin = histo[iHisto].FindBin(iPhi)
if (tmpBin not in regionBins[iHisto]):
regionBins[iHisto].append(tmpBin)
elif (histoType == "1d_integralAbove"):
regionBins[iHisto] = []
for iBin in range(histo[iHisto].FindBin(integralAbove),
histo[iHisto].GetNbinsX()):
regionBins[iHisto].append(iBin)
# Extract all the unmerged files available with the LB range
lbFilePathList = pathExtract.returnEosHistPathLB(runNumber, lowerLumiBlock,
upperLumiBlock, stream,
amiTag, tag)
nbHitInHot = []
nLB = 2500
nbHitInHot = {}
for iHisto in histoKeys:
nbHitInHot[iHisto] = [0.] * nLB
lowerLB = 2500
upperLB = 0
lbCanvas = []
histoLBNoisy = []
fLB = {}
print "I have found the merged HIST file %s" % (runFilePath)
print "I have found %d unmerged HIST files" % (len(lbFilePathList))
}
run_info = s.get_run_information(run_spec)
if '%d' % runNumber not in run_info.keys() or len(
run_info['%d' % runNumber]) < 2:
print "Unable to retrieve the data project tag via atlasdqm... Please double check your atlasdqmpass.txt or define it by hand with -t option"
sys.exit()
tag = run_info['%d' % runNumber][1]
amiTag = args.amiTag
path = []
if args.lumiblock == 0:
path.append("root://eosatlas.cern.ch/%s" % (pathExtract.returnEosHistPath(
args.runNumber, args.stream, args.amiTag, tag)).rstrip())
else:
allUnmerged = pathExtract.returnEosHistPathLB(runNumber, LB, LB, stream,
amiTag, tag)
for iFile in allUnmerged:
path.append("root://eosatlas.cern.ch/%s" % (iFile).rstrip())
# path.add("root://eosatlas.cern.ch/%s"%(.rstrip())
file = []
for iPath in path:
if ("NO FILE" not in iPath):
print "I am opening %s" % (iPath)
file.append(TFile.Open(iPath))
gStyle.SetPalette(1)
tb = TBrowser()