Lv3_analyze_xspec_pars.py
def plot_HID(MJDs):
    """
    Plot the soft color-intensity diagram for a given set of MJDs

    MJDs - list of MJDs used
    """
    mjd_data,soft,soft_err,intensity,intensity_err = np.genfromtxt(Lv0_dirs.NGC300+'soft_color_HID.txt',skip_header=1,usecols=(0,1,2,3,4),unpack=True)

    soft_trunc = [soft[i] for i in range(len(mjd_data)) if str(int(mjd_data[i])) in MJDs]
    soft_err_trunc = [soft_err[i] for i in range(len(mjd_data)) if str(int(mjd_data[i])) in MJDs]
    intensity_trunc = [intensity[i] for i in range(len(mjd_data)) if str(int(mjd_data[i])) in MJDs]
    intensity_err_trunc = [intensity_err[i] for i in range(len(mjd_data)) if str(int(mjd_data[i])) in MJDs]

    pdf_name = Lv0_dirs.NGC300 + 'spectral_fit_' + E1 + '-' + E2 + '/soft_color_HID.pdf'
    with PdfPages(pdf_name) as pdf:
        for i in tqdm(range(len(soft_trunc))):
            plt.figure(figsize=(16,9))
            plt.xlim([-0.1,2.9])
            plt.ylim([-0.9,1.6])
            plt.xlabel('Soft Color: 1.0-2.0 keV/0.4-1.0 keV',fontsize=12)
            plt.ylabel('Intensity (ct/s)',fontsize=12)
            plt.title('MJD: ' + str(MJDs[i]))
            plt.errorbar(x=soft_trunc,y=intensity_trunc,xerr=soft_err_trunc,yerr=intensity_err_trunc,fmt='^')
            plt.errorbar(x=soft_trunc[i],y=intensity_trunc[i],xerr=soft_err_trunc[i],yerr=intensity_err_trunc[i],fmt='+',color='k',markersize=15)
            pdf.savefig()
            plt.close()

    return

    Plot the soft color-intensity diagram for a given set of MJDs

    MJDs - list of MJDs used
    
Plot the soft color-intensity diagram for a given set of MJDs
MJDs - list of MJDs used
Plot the soft color - intensity diagram for a given set of MJDs MJDs - list of MJDs used
def plot_HID(MJDs):
    mjd_data,soft,soft_err,intensity,intensity_err = np.genfromtxt(Lv0_dirs.NGC300+'soft_color_HID.txt',skip_header=1,usecols=(0,1,2,3,4),unpack=True)
    soft_trunc = [soft[i] for i in range(len(mjd_data)) if str(int(mjd_data[i])) in MJDs]
    soft_err_trunc = [soft_err[i] for i in range(len(mjd_data)) if str(int(mjd_data[i])) in MJDs]
    intensity_trunc = [intensity[i] for i in range(len(mjd_data)) if str(int(mjd_data[i])) in MJDs]
    intensity_err_trunc = [intensity_err[i] for i in range(len(mjd_data)) if str(int(mjd_data[i])) in MJDs]
    pdf_name = Lv0_dirs.NGC300 + 'spectral_fit_' + E1 + '-' + E2 + '/soft_color_HID.pdf'
    with PdfPages(pdf_name) as pdf:
        for i in tqdm(range(len(soft_trunc))):
            plt.figure(figsize=(16,9))
            plt.xlim([-0.1,2.9])
            plt.ylim([-0.9,1.6])
            plt.xlabel('Soft Color: 1.0-2.0 keV/0.4-1.0 keV',fontsize=12)
            plt.ylabel('Intensity (ct/s)',fontsize=12)
            plt.title('MJD: ' + str(MJDs[i]))
            plt.errorbar(x=soft_trunc,y=intensity_trunc,xerr=soft_err_trunc,yerr=intensity_err_trunc,fmt='^')
            plt.errorbar(x=soft_trunc[i],y=intensity_trunc[i],xerr=soft_err_trunc[i],yerr=intensity_err_trunc[i],fmt='+',color='k',markersize=15)
            pdf.savefig()
            plt.close()
    return
def plot_HID ( MJDs ) : mjd_data , soft , soft_err , intensity , intensity_err = np . genfromtxt ( Lv0_dirs . NGC300 + 'soft_color_HID.txt' , skip_header = 1 , usecols = ( 0 , 1 , 2 , 3 , 4 ) , unpack = True ) soft_trunc = [ soft [ i ] for i in range ( len ( mjd_data ) ) if str ( int ( mjd_data [ i ] ) ) in MJDs ] soft_err_trunc = [ soft_err [ i ] for i in range ( len ( mjd_data ) ) if str ( int ( mjd_data [ i ] ) ) in MJDs ] intensity_trunc = [ intensity [ i ] for i in range ( len ( mjd_data ) ) if str ( int ( mjd_data [ i ] ) ) in MJDs ] intensity_err_trunc = [ intensity_err [ i ] for i in range ( len ( mjd_data ) ) if str ( int ( mjd_data [ i ] ) ) in MJDs ] pdf_name = Lv0_dirs . NGC300 + 'spectral_fit_' + E1 + '-' + E2 + '/soft_color_HID.pdf' with PdfPages ( pdf_name ) as pdf : for i in tqdm ( range ( len ( soft_trunc ) ) ) : plt . figure ( figsize = ( 16 , 9 ) ) plt . xlim ( [ - 0.1 , 2.9 ] ) plt . ylim ( [ - 0.9 , 1.6 ] ) plt . xlabel ( 'Soft Color: 1.0-2.0 keV/0.4-1.0 keV' , fontsize = 12 ) plt . ylabel ( 'Intensity (ct/s)' , fontsize = 12 ) plt . title ( 'MJD: ' + str ( MJDs [ i ] ) ) plt . errorbar ( x = soft_trunc , y = intensity_trunc , xerr = soft_err_trunc , yerr = intensity_err_trunc , fmt = '^' ) plt . errorbar ( x = soft_trunc [ i ] , y = intensity_trunc [ i ] , xerr = soft_err_trunc [ i ] , yerr = intensity_err_trunc [ i ] , fmt = '+' , color = 'k' , markersize = 15 ) pdf . savefig ( ) plt . close ( ) return
Plot the soft color-intensity diagram for a given set of MJDs
MJDs - list of MJDs used
Plot the soft color - intensity diagram for a given set of MJDs MJDs - list of MJDs used
"""
    Plot the soft color-intensity diagram for a given set of MJDs

    MJDs - list of MJDs used
    """
MJDs