Commit b10e83f1 authored by Andrea Giannetti's avatar Andrea Giannetti
Browse files

Fixed typo in reliability computation

parent 9644d03a

etl/mdl/mdl_extract_ch3oh_spectra.py

+1 −1
Changes for etl/mdl/mdl_extract_ch3oh_spectra.py: 1 added line, 1 removed line.
Original line number Diff line number Diff line
@@ -549,7 +549,7 @@ def process_source(source, species, guesses): 
            all_pars = {

                **snr_p,

                'ratio1': np.mean(r87),

                'ratio2': np.mean(r87),

                'ratio2': np.mean(r88),

                'rel_prec_ratio1': 1 - np.std(r87) / np.mean(r87),

                'rel_prec_ratio2': 1 - np.std(r88) / np.mean(r88)

            }

etl/mdl/mdl_extract_clump_spectra_and_densities.py

+9 −4
Changes for etl/mdl/mdl_extract_clump_spectra_and_densities.py: 9 added lines, 4 removed lines.
Original line number Diff line number Diff line
@@ -53,16 +53,21 @@ def process_source(source, species, guesses): 
    densities, fwhms, reliability_score = {}, {}, {}

    region = 'clump'

    coords, components = regions[source][region][:3], regions[source][region][3]

    densities[region] = {'best_fit': [], 'uncertainty': [], 'reliability_score': [], 'fwhm': [], 'vlsr': []}

    fwhms[region], reliability_score[region] = [], {}



    # Initialize storage for this single region

    densities = {region: {'best_fit': [], 'uncertainty': [], 'reliability_score': [], 'fwhm': [], 'vlsr': []}}

    fwhms = {region: []}

    reliability_score = {region: {}}



    apspec = pyspeckit.Spectrum(data=mean_spectrum, xarr=cube.xarr,

                                xarrkwargs={'unit': 'Hz'},

                                unit='K')



    tied, line_guesses, limits, limited = prepare_fit_guesses_and_ties(

        source, region, components, guesses, delta_freq, rms, pars_per_component, freqs_dop,

        channelwidth_hz, coords, max_lw_hz

    )



    spec_rms = run_and_plot_specfit(

        apspec=apspec,

        tied=tied,
@@ -76,6 +81,7 @@ def process_source(source, species, guesses): 
        name_prefix="clump"

    )



    # Process each velocity component

    for vel_comp in range(components):

        idx = vel_comp * 4
@@ -97,7 +103,7 @@ def process_source(source, species, guesses): 
        all_pars = {

            **snr_p,

            'ratio1': np.mean(r87),

            'ratio2': np.mean(r87),

            'ratio2': np.mean(r88),  # Note: this looks like it should be r88?

            'rel_prec_ratio1': 1 - np.std(r87) / np.mean(r87),

            'rel_prec_ratio2': 1 - np.std(r88) / np.mean(r88)

        }
@@ -129,7 +135,6 @@ def process_source(source, species, guesses): 
    with open(os.path.join('mdl', 'inferred_densities', f'clump_densities_{source}.yml'), 'w') as outfile:

        yaml.dump(densities, outfile, default_flow_style=False)





# sources = ['G14', 'G28', 'G332.96']

sources = ['G14', 'G28']

species = 'CH3OH'