Ejemplo n.º 1
0
posTrkD0 = array('d', [0])
posTrkLambda = array('d', [0])
posTrkZ0 = array('d', [0])
posTrkD0Err = array('d', [0])
posTrkLambdaErr = array('d', [0])
posTrkZ0Err = array('d', [0])
eleL1tthetaY = array('d', [0])
eleL2tthetaY = array('d', [0])
eleL1bthetaY = array('d', [0])
eleL2bthetaY = array('d', [0])
posL1tthetaY = array('d', [0])
posL2tthetaY = array('d', [0])
posL1bthetaY = array('d', [0])
posL2bthetaY = array('d', [0])

events.Branch("uncVX", uncVX, "uncVX")
events.Branch("uncVY", uncVY, "uncVY")
events.Branch("uncVZ", uncVZ, "uncVZ")
events.Branch("uncP", uncP, "uncP")
events.Branch("uncChisq", uncChisq, "uncChisq")
events.Branch("uncM", uncM, "uncM")
events.Branch("uncTargProjX", uncTargProjX, "uncTargProjX")
events.Branch("uncTargProjY", uncTargProjY, "uncTargProjY")
events.Branch("uncTargProjXErr", uncTargProjXErr, "uncTargProjXErr")
events.Branch("uncTargProjYErr", uncTargProjYErr, "uncTargProjYErr")
events.Branch("uncCovXX", uncCovXX, "uncCovXX")
events.Branch("uncCovYY", uncCovYY, "uncCovYY")
events.Branch("uncCovZZ", uncCovZZ, "uncCovZZ")
events.Branch("bscChisq", bscChisq, "bscChisq")
events.Branch("tarChisq", tarChisq, "tarChisq")
events.Branch("eleP", eleP, "eleP")
Ejemplo n.º 2
0
tuplename = "ntuple"

options, remainder = getopt.gnu_getopt(sys.argv[1:], 'n:h')

# Parse the command line arguments
for opt, arg in options:
    if opt == '-n':
        tuplename = str(arg)
    if opt == '-h':
        print_usage()
        sys.exit(0)

outfile = remainder[0]
events = TChain("{0}".format(tuplename))
for i in range(1, len(remainder)):
    events.Add(remainder[i])

run = array('d', [0])
event = array('d', [0])

events.Branch("run", run, "run")
events.Branch("event", event, "event")

txtfile = open(outfile + ".txt", "w")

for entry in xrange(events.GetEntries()):
    events.GetEntry(entry)
    txtfile.write(str(events.run) + " " + str(events.event) + "\n")

txtfile.close()