#!/usr/bin/env python # ----------------------------------------------------------------------------------- # # # ROOT macro for reading data for the Applied Statistics exam 2016 problem 4.1. # # Author: Troels C. Petersen (NBI) # Email: petersen@nbi.dk # Date: 11th of January 2017 # # ----------------------------------------------------------------------------------- # from ROOT import * from array import array import math def sqr(a) : return a*a # ----------------------------------------------------------------------------------- # # ----------------------------------------------------------------------------------- # # Setting what to be shown in statistics box: gStyle.SetOptStat("eMR") gStyle.SetOptFit(1111) # Text to write on the plot: text = TLatex() text.SetNDC() pi = TMath.Pi() SavePlots = False # ----------------------------------------------------------------------------------- # # Read data: # ----------------------------------------------------------------------------------- # data = [] with open( 'data_FisherSyndrome.txt', 'r' ) as infile : counter = 0 for line in infile: line = line.strip().split() data.append([float(line[2]), float(line[3]), float(line[4]), int(line[0])]) isIll = data[-1][3] A = data[-1][0] B = data[-1][1] C = data[-1][2] # Print some numbers as a sanity check if (counter < 10) : print" Read data: %d %5.2f %5.2f %5.2f"%(isIll, A, B, C) counter += 1 raw_input('Press Enter to exit') #---------------------------------------------------------------------------------- # # #----------------------------------------------------------------------------------