array_PH = np.array([ tree.PH1, tree.PH2, tree.PH3, tree.PH4, tree.PH5, tree.PH6, tree.PH7, tree.PH8, tree.PH9, tree.PH10, tree.PH11, tree.PH12, tree.PH13, tree.PH14, tree.PH15, tree.PH16, tree.PH17, tree.PH18, tree.PH19, tree.PH20, tree.PH21, tree.PH22, tree.PH23, tree.PH24, tree.PH25, tree.PH26, tree.PH27, tree.PH28, tree.PH29, tree.PH30, tree.PH31, tree.PH32, tree.PH33, tree.PH34, tree.PH35, tree.PH36, tree.PH37, tree.PH38, tree.PH39, tree.PH40, tree.PH41, tree.PH42, tree.PH43, tree.PH44, tree.PH45, tree.PH46, tree.PH47, tree.PH48, tree.PH49 ]) for i in range(49): if array_PH[i] > 9.279: Xray_PH += array_PH[i] count += 1 print 'count = %s' % count my.Arg = tree.Arg my.Frame = tree.Frame my.Event_th = tree.Event_th my.Split_th = 9.279 my.X = tree.X my.Y = tree.Y my.Xray = Xray_PH my.count = count tr.Fill() array_PH = np.array([]) f.Write() f.Close()