if max(row) > 0: psf_matias_single = temp_psf del indices_single, new_data, temp_psf, row break del best_cells ################### #for chi2 expected = sample_size*phi_tilde_pdf*.01 ################ z_dist = df.max_finder(psf_base_single,phi_tilde_point) base_point_dist.append(z_dist) ################ z_dist, c_val = df.integral_calc(initial_guess,psf_base_t,phi_tilde_pdf) z_temp = z_dist for j in xrange(max_iters): #expected = sample_size*phi_tilde_pdf*.01 observed = sample_size*c_val*.01 chisq_stat = chisquare(observed,expected)[0] if chisq_stat < 233: base_pdf_dist.append(z_temp)