float(averages[iring][s]) + float(samples[s]) for s in xrange(len(samples)) ] filledCry[iring] += 1 print "Starting the second loop..." eenewdata = {} partition = '0' total = 0 unchanged = 0 averaged = 0 averagedtoneighbor = 0 for detid, val in detids.geometry.iteritems(): key = (partition, detid) (x, y, z) = detids.xyz(detid) if detids.subdet(detid) == 1: continue # barrel x0 = x - 50 y0 = y - 50 ir = int(math.sqrt(x0 * x0 + y0 * y0)) iring = ringmap.getRing(False, x, y, z) #if iring==-1: # continue # not valid x,y total += 1 if ir > maxIr: if key in eedata: eenewdata[key] = eedata[key] unchanged += 1 else: pass # leave empty crystals that have no measurements in low eta EE else: if iring not in filledCry: