fout = ROOT.TFile("%s/fits/rootFiles.root"%rootDirec,"recreate") ipwCan = ROOT.TCanvas() pinCan = ROOT.TCanvas() ipwCan.SetCanvasSize(1000,400) pinCan.SetCanvasSize(1000,400) ipwCan.Divide(2,1) pinCan.Divide(2,1) # Loop over all data using standarised directory structure resultsList = [] boxes = return_boxes("%s/broad_sweep/"%rootDirec) for box in boxes: broadFiles = return_files("%s/broad_sweep"%rootDirec, box) lowFiles = return_files("%s/low_intensity"%rootDirec, box) for j in range(len(lowFiles)): broadVals = check_data(plot_ipw.read_scope_scan(broadFiles[j])) lowVals = check_data(plot_ipw.read_scope_scan(lowFiles[j])) # Creat plots photonVsPIN_broad = ROOT.TGraphErrors() photonVsIPW_broad = ROOT.TGraphErrors() photonVsPIN_low = ROOT.TGraphErrors() photonVsIPW_low = ROOT.TGraphErrors() print len(broadVals), len(lowVals), broadFiles[j], lowFiles[j] for i in range(len(broadVals)): photonBroad = plot_ipw.get_photons(broadVals[i]["area"], 0.5) photonErrBroad = plot_ipw.get_photons(broadVals[i]["area_err"], 0.5) # Fill plots with data # Note: Data points are returned as mean and stdev(rms) # for fitting, uncertainties should be given as standard error. photonVsPIN_broad.SetPoint(i,broadVals[i]["pin"],photonBroad) photonVsPIN_broad.SetPointError(i,broadVals[i]["pin_err"]/np.sqrt(100),photonErrBroad/np.sqrt(100))
fout = ROOT.TFile("%s/fits/rootFiles.root" % rootDirec, "recreate") ipwCan = ROOT.TCanvas() pinCan = ROOT.TCanvas() ipwCan.SetCanvasSize(1000, 400) pinCan.SetCanvasSize(1000, 400) ipwCan.Divide(2, 1) pinCan.Divide(2, 1) # Loop over all data using standarised directory structure resultsList = [] boxes = return_boxes("%s/broad_sweep/" % rootDirec) for box in boxes: broadFiles = return_files("%s/broad_sweep" % rootDirec, box) lowFiles = return_files("%s/low_intensity" % rootDirec, box) for j in range(len(lowFiles)): broadVals = check_data(plot_ipw.read_scope_scan(broadFiles[j])) lowVals = check_data(plot_ipw.read_scope_scan(lowFiles[j])) # Creat plots photonVsPIN_broad = ROOT.TGraphErrors() photonVsIPW_broad = ROOT.TGraphErrors() photonVsPIN_low = ROOT.TGraphErrors() photonVsIPW_low = ROOT.TGraphErrors() print len(broadVals), len(lowVals), broadFiles[j], lowFiles[j] for i in range(len(broadVals)): photonBroad = plot_ipw.get_photons(broadVals[i]["area"], 0.5) photonErrBroad = plot_ipw.get_photons(broadVals[i]["area_err"], 0.5) # Fill plots with data # Note: Data points are returned as mean and stdev(rms) # for fitting, uncertainties should be given as standard error. photonVsPIN_broad.SetPoint(i, broadVals[i]["pin"], photonBroad)
# Date : 11/05/15 ####################################### import calc_utils as calc import plot_ipw import optparse import os if __name__ == "__main__": parser = optparse.OptionParser("Usage: python checkMeasurements.py [option]") parser.add_option("-f", "--file", dest="file", default=False, help="Pass path to data file to be checked") (options,args) = parser.parse_args() res_list = plot_ipw.read_scope_scan(options.file) p = options.file.split('/') print p box = int(p[2][4:6]) chan = int(p[3][4:6]) print chan overallChan = ((box-1)*8)+chan dataPath = os.path.join(p[0],p[1],"Box_%02d/raw_data/channel_%02d/" % (box,overallChan)) files = [ f for f in os.listdir(dataPath) if os.path.isfile(os.path.join(dataPath,f)) ] output_filename = "%s_check.dat" % options.file[0:-4] output_file = file(output_filename,'w') output_file.write("#PWIDTH\tPWIDTH Error\tPIN\tPIN Error\tWIDTH\tWIDTH Error\tRISE\tRISE Error\tFALL\tFALL Error\tAREA\tAREA Error\t\ Minimum\tMinimum Error\n")