if eMin<100000: E = (like.energies[:-1] + like.energies[1:])/2. # The 'energies' array are the endpoints so we take the midpoint of the bins. plt.figure(figsize=(9,9)) plt.ylim((0.4,1e4)) plt.xlim((200,300000)) sum_model = np.zeros_like(like._srcCnts(like.sourceNames()[0])) for sourceName in like.sourceNames(): sum_model = sum_model + like._srcCnts(sourceName) plt.loglog(E,like._srcCnts(sourceName),label=sourceName[1:]) plt.loglog(E,sum_model,label='Total Model') plt.errorbar(E,like._Nobs(),yerr=np.sqrt(like._Nobs()), fmt='o',label='Counts') plt.legend(bbox_to_anchor=(1.05, 1), loc=2) plt.savefig('results/1.eps',format='eps', bbox_inches='tight') # Save figure! # Plot residuals sum_counts=sum_model # Is this right? Probably not :/ resid = (like._Nobs() - sum_counts)/sum_counts resid_err = (np.sqrt(like._Nobs())/sum_counts) plt.figure(figsize=(9,9)) plt.xscale('log') plt.errorbar(E,resid,yerr=resid_err,fmt='o') plt.axhline(0.0,ls=':') plt.savefig('results/2.eps',format='eps', bbox_inches='tight') # Get indexes and stuff
E = (like.energies[:-1] + like.energies[1:]) / 2. # The 'energies' array are the endpoints so we take the midpoint of the bins. plt.figure(figsize=(9, 9)) plt.ylim((0.4, 1e4)) plt.xlim((200, 300000)) sum_model = np.zeros_like(like._srcCnts(like.sourceNames()[0])) for sourceName in like.sourceNames(): sum_model = sum_model + like._srcCnts(sourceName) plt.loglog(E, like._srcCnts(sourceName), label=sourceName[1:]) plt.loglog(E, sum_model, label='Total Model') plt.errorbar(E, like._Nobs(), yerr=np.sqrt(like._Nobs()), fmt='o', label='Counts') plt.legend(bbox_to_anchor=(1.05, 1), loc=2) plt.savefig('results/1.eps', format='eps', bbox_inches='tight') # Save figure! # Plot residuals sum_counts = sum_model # Is this right? Probably not :/ resid = (like._Nobs() - sum_counts) / sum_counts resid_err = (np.sqrt(like._Nobs()) / sum_counts) plt.figure(figsize=(9, 9)) plt.xscale('log') plt.errorbar(E, resid, yerr=resid_err, fmt='o') plt.axhline(0.0, ls=':') plt.savefig('results/2.eps', format='eps', bbox_inches='tight')