if (index == -1):
PS3 = 0
else:
PS3 = psActual[i]
if (index == ps5):
if (index == -1):
PS5 = 0
else:
PS5 = psActual[i]
f.close()
'''
SCALER TREE, TSH
'''
s_tree = up.open(rootName)["TSP"]
s_branch = klt.pyBranch(s_tree)
# s_evts = len(s_tree)
s_evts = s_tree.array("P.BCM4A.scaler")
P_BCM4A_scalerCharge = s_tree.array("P.BCM4A.scalerCharge")
P_BCM2_scalerCharge = s_tree.array("P.BCM2.scalerCharge")
P_BCM4B_scalerCharge = s_tree.array("P.BCM4B.scalerCharge")
P_BCM1_scalerCharge = s_tree.array("P.BCM1.scalerCharge")
P_BCM4C_scalerCharge = s_tree.array("P.BCM4C.scalerCharge")
P_BCM4A_scalerCurrent = s_tree.array("P.BCM4A.scalerCurrent")
P_BCM2_scalerCurrent = s_tree.array("P.BCM2.scalerCurrent")
P_BCM4B_scalerCurrent = s_tree.array("P.BCM4B.scalerCurrent")
P_BCM1_scalerCurrent = s_tree.array("P.BCM1.scalerCurrent")
P_BCM4C_scalerCurrent = s_tree.array("P.BCM4C.scalerCurrent")
REPLAYPATH = "/home/%s/work/JLab/hallc_replay_lt" % USER[1]
elif ("trottar" in HOST[1]):
REPLAYPATH = "/home/trottar/Analysis/hallc_replay_lt"
print("Running as %s on %s, hallc_replay_lt path assumed as %s" %
(USER[1], HOST[1], REPLAYPATH))
filename = "%s/UTIL_PION/OUTPUT/Analysis/PID/pid_data.csv" % (REPLAYPATH)
rootName = "%s/UTIL_PION/ROOTfiles/Analysis/PID/pid_coin_offline_%s_%s.root" % (
REPLAYPATH, runNum, MaxEvent)
'''
ANALYSIS TREE, T
'''
tree = up.open(rootName)["T"]
branch = klt.pyBranch(tree)
CTime_ePositronCoinTime_ROC1 = tree.array("CTime.ePositronCoinTime_ROC1")
CTime_eKCoinTime_ROC1 = tree.array("CTime.eKCoinTime_ROC1")
CTime_ePiCoinTime_ROC1 = tree.array("CTime.ePiCoinTime_ROC1")
CTime_epCoinTime_ROC1 = tree.array("CTime.epCoinTime_ROC1")
P_gtr_beta = tree.array("P.gtr.beta")
P_gtr_th = tree.array("P.gtr.th")
P_gtr_ph = tree.array("P.gtr.ph")
H_gtr_beta = tree.array("H.gtr.beta")
H_gtr_th = tree.array("H.gtr.th")
H_gtr_ph = tree.array("H.gtr.ph")
H_cal_etotnorm = tree.array("H.cal.etotnorm")
H_cer_npeSum = tree.array("H.cer.npeSum")
P_cal_etotnorm = tree.array("P.cal.etotnorm")
P_aero_npeSum = tree.array("P.aero.npeSum")