OpenREM 1.0.0b2__py3-none-any.whl → 1.0.0b3__py3-none-any.whl

openrem/remapp/views_charts_mg.py

@@ -1,1196 +1,1349 @@
-# pylint: disable=too-many-lines
-import logging
-from datetime import datetime
-from django.conf import settings
-from django.contrib.auth.decorators import login_required
-from django.core.exceptions import ObjectDoesNotExist
-from django.http import JsonResponse
-from remapp.forms import (
-    MGChartOptionsForm,
-    MGChartOptionsFormIncStandard,
-)
-from remapp.interface.mod_filters import (
-    MGSummaryListFilter,
-    MGFilterPlusPid,
-    MGFilterPlusStdNames,
-    MGFilterPlusPidPlusStdNames,
-)
-from remapp.models import (
-    GeneralStudyModuleAttr,
-    create_user_profile,
-    StandardNameSettings,
-)
-from remapp.views_admin import (
-    required_average_choices,
-    initialise_mg_form_data,
-    set_average_chart_options,
-    set_mg_chart_options,
-    set_common_chart_options,
-)
-from .interface.chart_functions import (
-    create_dataframe,
-    create_dataframe_weekdays,
-    plotly_barchart_weekdays,
-    plotly_binned_statistic_barchart,
-    plotly_boxplot,
-    plotly_barchart,
-    plotly_histogram_barchart,
-    plotly_scatter,
-    plotly_frequency_barchart,
-    construct_over_time_charts,
-    plotly_set_default_theme,
-    create_dataframe_aggregates,
-    generate_average_chart_group,
-)
-
-logger = logging.getLogger(__name__)
-
-
-def generate_required_mg_charts_list(profile):
-    # pylint: disable=too-many-branches
-    # pylint: disable=too-many-statements
-    """Obtain a list of dictionaries containing the title string and base
-    variable name for each required chart"""
-
-    # Obtain the system-level enable_standard_names setting
-    try:
-        StandardNameSettings.objects.get()
-    except ObjectDoesNotExist:
-        StandardNameSettings.objects.create()
-    enable_standard_names = StandardNameSettings.objects.values_list(
-        "enable_standard_names", flat=True
-    )[0]
-
-    required_charts = []
-
-    charts_of_interest = [
-        profile.plotMGAcquisitionAGDOverTime,
-    ]
-    if enable_standard_names:
-        charts_of_interest.append(profile.plotMGStandardAcquisitionAGDOverTime)
-
-    if any(charts_of_interest):
-        keys = list(dict(profile.TIME_PERIOD).keys())
-        values = list(dict(profile.TIME_PERIOD).values())
-        time_period = (values[keys.index(profile.plotMGOverTimePeriod)]).lower()
-
-    if profile.plotMGacquisitionFreq:
-        required_charts.append(
-            {
-                "title": "Chart of acquisition protocol frequency",
-                "var_name": "acquisitionFrequency",
-            }
-        )
-
-    if profile.plotMGaverageAGD:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of acquisition protocol mean AGD",
-                    "var_name": "acquisitionMeanAGD",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of acquisition protocol median AGD",
-                    "var_name": "acquisitionMedianAGD",
-                }
-            )
-        if profile.plotBoxplots:
-            required_charts.append(
-                {
-                    "title": "Boxplot of acquisition protocol AGD",
-                    "var_name": "acquisitionBoxplotAGD",
-                }
-            )
-        if profile.plotHistograms:
-            required_charts.append(
-                {
-                    "title": "Histogram of acquisition protocol AGD",
-                    "var_name": "acquisitionHistogramAGD",
-                }
-            )
-
-    if profile.plotMGAGDvsThickness:
-        required_charts.append(
-            {
-                "title": "Chart of acquisition protocol AGD vs compressed breast thickness",
-                "var_name": "acquisitionScatterAGDvsThick",
-            }
-        )
-
-    if profile.plotMGaverageAGDvsThickness:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of acquisition protocol mean AGD vs compressed breast thickness",
-                    "var_name": "acquisitionMeanAGDvsThick",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of acquisition protocol median AGD vs compressed breast thickness",
-                    "var_name": "acquisitionMedianAGDvsThick",
-                }
-            )
-
-    if profile.plotMGAcquisitionAGDOverTime:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of acquisition protocol mean AGD over time ("
-                    + time_period
-                    + ")",
-                    "var_name": "acquisitionMeanAGDOverTime",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of acquisition protocol median AGD over time ("
-                    + time_period
-                    + ")",
-                    "var_name": "acquisitionMedianAGDOverTime",
-                }
-            )
-
-    if profile.plotMGkVpvsThickness:
-        required_charts.append(
-            {
-                "title": "Chart of acquisition protocol kVp vs compressed breast thickness",
-                "var_name": "acquisitionScatterkVpvsThick",
-            }
-        )
-
-    if profile.plotMGmAsvsThickness:
-        required_charts.append(
-            {
-                "title": "Chart of acquisition protocol mAs vs compressed breast thickness",
-                "var_name": "acquisitionScattermAsvsThick",
-            }
-        )
-
-    if profile.plotMGStudyPerDayAndHour:
-        required_charts.append(
-            {
-                "title": "Chart of study description workload",
-                "var_name": "studyWorkload",
-            }
-        )
-
-    if enable_standard_names:
-        if profile.plotMGStandardAcquisitionFreq:
-            required_charts.append(
-                {
-                    "title": "Chart of standard acquisition name frequency",
-                    "var_name": "standardAcquisitionFrequency",
-                }
-            )
-
-        if profile.plotMGStandardAverageAGD:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of standard acquisition name mean AGD",
-                        "var_name": "standardAcquisitionMeanAGD",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of standard acquisition name median AGD",
-                        "var_name": "standardAcquisitionMedianAGD",
-                    }
-                )
-            if profile.plotBoxplots:
-                required_charts.append(
-                    {
-                        "title": "Boxplot of standard acquisition name AGD",
-                        "var_name": "standardAcquisitionBoxplotAGD",
-                    }
-                )
-            if profile.plotHistograms:
-                required_charts.append(
-                    {
-                        "title": "Histogram of standard acquisition name AGD",
-                        "var_name": "standardAcquisitionHistogramAGD",
-                    }
-                )
-
-        if profile.plotMGStandardAGDvsThickness:
-            required_charts.append(
-                {
-                    "title": "Chart of standard acquisition name AGD vs compressed breast thickness",
-                    "var_name": "standardAcquisitionScatterAGDvsThick",
-                }
-            )
-
-        if profile.plotMGStandardAverageAGDvsThickness:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of standard acquisition name mean AGD vs compressed breast thickness",
-                        "var_name": "standardAcquisitionMeanAGDvsThick",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of standard acquisition name median AGD vs compressed breast thickness",
-                        "var_name": "standardAcquisitionMedianAGDvsThick",
-                    }
-                )
-
-        if profile.plotMGStandardAcquisitionAGDOverTime:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of standard acquisition name mean AGD over time ("
-                        + time_period
-                        + ")",
-                        "var_name": "standardAcquisitionMeanAGDOverTime",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of standard acquisition name median AGD over time ("
-                        + time_period
-                        + ")",
-                        "var_name": "standardAcquisitionMedianAGDOverTime",
-                    }
-                )
-
-        if profile.plotMGStandardkVpvsThickness:
-            required_charts.append(
-                {
-                    "title": "Chart of standard acquisition name kVp vs compressed breast thickness",
-                    "var_name": "standardAcquisitionScatterkVpvsThick",
-                }
-            )
-
-        if profile.plotMGStandardmAsvsThickness:
-            required_charts.append(
-                {
-                    "title": "Chart of standard acquisition name mAs vs compressed breast thickness",
-                    "var_name": "standardAcquisitionScattermAsvsThick",
-                }
-            )
-
-        if profile.plotMGStandardStudyPerDayAndHour:
-            required_charts.append(
-                {
-                    "title": "Chart of standard study name workload",
-                    "var_name": "standardStudyWorkload",
-                }
-            )
-
-    return required_charts
-
-
-@login_required
-def mg_summary_chart_data(request):
-    """Obtain data for mammography chart data Ajax view"""
-
-    # Obtain the system-level enable_standard_names setting
-    try:
-        StandardNameSettings.objects.get()
-    except ObjectDoesNotExist:
-        StandardNameSettings.objects.create()
-    enable_standard_names = StandardNameSettings.objects.values_list(
-        "enable_standard_names", flat=True
-    )[0]
-
-    if request.user.groups.filter(name="pidgroup"):
-        if enable_standard_names:
-            f = MGFilterPlusPidPlusStdNames(
-                request.GET,
-                queryset=GeneralStudyModuleAttr.objects.filter(
-                    modality_type__exact="MG"
-                )
-                .order_by()
-                .distinct(),
-            )
-        else:
-            f = MGFilterPlusPid(
-                request.GET,
-                queryset=GeneralStudyModuleAttr.objects.filter(
-                    modality_type__exact="MG"
-                )
-                .order_by()
-                .distinct(),
-            )
-    else:
-        if enable_standard_names:
-            f = MGFilterPlusStdNames(
-                request.GET,
-                queryset=GeneralStudyModuleAttr.objects.filter(
-                    modality_type__exact="MG"
-                )
-                .order_by()
-                .distinct(),
-            )
-        else:
-            f = MGSummaryListFilter(
-                request.GET,
-                queryset=GeneralStudyModuleAttr.objects.filter(
-                    modality_type__exact="MG"
-                )
-                .order_by()
-                .distinct(),
-            )
-    try:
-        # See if the user has plot settings in userprofile
-        user_profile = request.user.userprofile
-    except ObjectDoesNotExist:
-        # Create a default userprofile for the user if one doesn't exist
-        create_user_profile(sender=request.user, instance=request.user, created=True)
-        user_profile = request.user.userprofile
-
-    if settings.DEBUG:
-        start_time = datetime.now()
-
-    return_structure = mg_plot_calculations(f, user_profile)
-
-    if settings.DEBUG:
-        logger.debug(f"Elapsed time is {datetime.now() - start_time}")
-
-    return JsonResponse(return_structure, safe=False)
-
-
-def mg_plot_calculations(f, user_profile, return_as_dict=False):
-    # pylint: disable=too-many-locals
-    # pylint: disable=too-many-branches
-    # pylint: disable=too-many-statements
-    """Calculations for mammography charts"""
-    # Return an empty structure if the queryset is empty
-    if not f.qs.exists():
-        return {}
-
-    # Obtain the system-level enable_standard_names setting
-    try:
-        StandardNameSettings.objects.get()
-    except ObjectDoesNotExist:
-        StandardNameSettings.objects.create()
-    enable_standard_names = StandardNameSettings.objects.values_list(
-        "enable_standard_names", flat=True
-    )[0]
-
-    # Set the Plotly chart theme
-    plotly_set_default_theme(user_profile.plotThemeChoice)
-
-    return_structure = {}
-
-    average_choices = []
-    if user_profile.plotMean:
-        average_choices.append("mean")
-    if user_profile.plotMedian:
-        average_choices.append("median")
-
-    charts_of_interest = [user_profile.plotMGAcquisitionAGDOverTime]
-    if enable_standard_names:
-        charts_of_interest.append(user_profile.plotMGStandardAcquisitionAGDOverTime)
-    if any(charts_of_interest):
-        plot_timeunit_period = user_profile.plotMGOverTimePeriod
-
-    #######################################################################
-    # Prepare acquisition-level Pandas DataFrame to use for charts
-    charts_of_interest = [
-        user_profile.plotMGAGDvsThickness,
-        user_profile.plotMGkVpvsThickness,
-        user_profile.plotMGmAsvsThickness,
-        user_profile.plotMGaverageAGDvsThickness,
-        user_profile.plotMGaverageAGD,
-        user_profile.plotMGacquisitionFreq,
-        user_profile.plotMGAcquisitionAGDOverTime,
-    ]
-    if enable_standard_names:
-        charts_of_interest.append(user_profile.plotMGStandardAGDvsThickness)
-        charts_of_interest.append(user_profile.plotMGStandardkVpvsThickness)
-        charts_of_interest.append(user_profile.plotMGStandardmAsvsThickness)
-        charts_of_interest.append(user_profile.plotMGStandardAverageAGDvsThickness)
-        charts_of_interest.append(user_profile.plotMGStandardAverageAGD)
-        charts_of_interest.append(user_profile.plotMGStandardAcquisitionFreq)
-        charts_of_interest.append(user_profile.plotMGStandardAcquisitionAGDOverTime)
-
-    if any(charts_of_interest):  # pylint: disable=too-many-nested-blocks
-
-        name_fields = []
-        charts_of_interest = [
-            user_profile.plotMGAGDvsThickness,
-            user_profile.plotMGkVpvsThickness,
-            user_profile.plotMGmAsvsThickness,
-            user_profile.plotMGaverageAGDvsThickness,
-            user_profile.plotMGaverageAGD,
-            user_profile.plotMGacquisitionFreq,
-            user_profile.plotMGAcquisitionAGDOverTime,
-        ]
-        if any(charts_of_interest):
-            name_fields.append(
-                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
-            )
-
-        if enable_standard_names:
-            charts_of_interest = [
-                user_profile.plotMGStandardAGDvsThickness,
-                user_profile.plotMGStandardkVpvsThickness,
-                user_profile.plotMGStandardmAsvsThickness,
-                user_profile.plotMGStandardAverageAGDvsThickness,
-                user_profile.plotMGStandardAverageAGD,
-                user_profile.plotMGStandardAcquisitionFreq,
-                user_profile.plotMGStandardAcquisitionAGDOverTime,
-            ]
-            if any(charts_of_interest):
-                name_fields.append(
-                    "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-                )
-
-        value_fields = []
-        value_multipliers = []
-        charts_of_interest = [
-            user_profile.plotMGAGDvsThickness,
-            user_profile.plotMGaverageAGDvsThickness,
-            user_profile.plotMGaverageAGD,
-            user_profile.plotMGAcquisitionAGDOverTime,
-        ]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotMGStandardAGDvsThickness)
-            charts_of_interest.append(user_profile.plotMGStandardAverageAGDvsThickness)
-            charts_of_interest.append(user_profile.plotMGStandardAverageAGD)
-            charts_of_interest.append(user_profile.plotMGStandardAcquisitionAGDOverTime)
-        if any(charts_of_interest):
-            value_fields.append(
-                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose"
-            )
-            value_multipliers.append(1)
-
-        charts_of_interest = [user_profile.plotMGkVpvsThickness]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotMGStandardkVpvsThickness)
-        if any(charts_of_interest):
-            value_fields.append(
-                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp"
-            )
-            value_multipliers.append(1)
-
-        charts_of_interest = [user_profile.plotMGmAsvsThickness]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotMGStandardmAsvsThickness)
-        if any(charts_of_interest):
-            value_fields.append(
-                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure"
-            )
-            value_multipliers.append(0.001)
-
-        charts_of_interest = [
-            user_profile.plotMGAGDvsThickness,
-            user_profile.plotMGkVpvsThickness,
-            user_profile.plotMGmAsvsThickness,
-            user_profile.plotMGaverageAGDvsThickness,
-        ]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotMGStandardAGDvsThickness)
-            charts_of_interest.append(user_profile.plotMGStandardkVpvsThickness)
-            charts_of_interest.append(user_profile.plotMGStandardmAsvsThickness)
-            charts_of_interest.append(user_profile.plotMGStandardAverageAGDvsThickness)
-        if any(charts_of_interest):
-            value_fields.append(
-                "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness"
-            )
-            value_multipliers.append(1)
-
-        date_fields = []
-        charts_of_interest = [user_profile.plotMGAcquisitionAGDOverTime]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotMGStandardAcquisitionAGDOverTime)
-        if any(charts_of_interest):
-            date_fields.append("study_date")
-
-        time_fields = []
-
-        system_field = []
-        if user_profile.plotSeriesPerSystem:
-            system_field.append(
-                "generalequipmentmoduleattr__unique_equipment_name_id__display_name"
-            )
-
-        fields = {
-            "names": name_fields,
-            "values": value_fields,
-            "dates": date_fields,
-            "times": time_fields,
-            "system": system_field,
-        }
-
-        df = create_dataframe(
-            f.qs,
-            fields,
-            data_point_name_lowercase=user_profile.plotCaseInsensitiveCategories,
-            data_point_name_remove_whitespace_padding=user_profile.plotRemoveCategoryWhitespacePadding,
-            data_point_value_multipliers=value_multipliers,
-            char_wrap=user_profile.plotLabelCharWrap,
-            uid="projectionxrayradiationdose__irradeventxraydata__pk",
-        )
-
-        if user_profile.plotMGaverageAGDvsThickness or user_profile.plotMGaverageAGD:
-
-            if user_profile.plotBoxplots and "median" not in average_choices:
-                average_choices = average_choices + ["median"]
-
-            name_field = (
-                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
-            )
-            value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose"  # pylint: disable=line-too-long
-
-            if user_profile.plotMGaverageAGDvsThickness:
-                category_names_col = name_field
-                group_by_col = "x_ray_system_name"
-                legend_title = "Acquisition protocol"
-
-                if user_profile.plotGroupingChoice == "series":
-                    category_names_col = "x_ray_system_name"
-                    group_by_col = name_field
-                    legend_title = "System"
-
-                parameter_dict = {  # pylint: disable=line-too-long
-                    "df_x_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
-                    "df_y_value_col": value_field,  # pylint: disable=line-too-long
-                    "x_axis_title": "Compressed breast thickness (mm)",
-                    "y_axis_title": "AGD (mGy)",
-                    "df_category_col": category_names_col,
-                    "df_facet_col": group_by_col,
-                    "facet_title": legend_title,
-                    "user_bins": [20, 30, 40, 50, 60, 70, 80, 90],
-                    "colourmap": user_profile.plotColourMapChoice,
-                    "filename": "OpenREM CT acquisition protocol average AGD vs thickness",
-                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                    "sorting_choice": [
-                        user_profile.plotInitialSortingDirection,
-                        user_profile.plotMGInitialSortingChoice,
-                    ],
-                    "return_as_dict": return_as_dict,
-                }
-                if user_profile.plotMean:
-                    parameter_dict["stat_name"] = "mean"
-                    return_structure[
-                        "meanAGDvsThickness"
-                    ] = plotly_binned_statistic_barchart(
-                        df,
-                        parameter_dict,
-                    )
-
-                if user_profile.plotMedian:
-                    parameter_dict["stat_name"] = "median"
-                    return_structure[
-                        "medianAGDvsThickness"
-                    ] = plotly_binned_statistic_barchart(
-                        df,
-                        parameter_dict,
-                    )
-
-            if user_profile.plotMGaverageAGD:
-                value_text = "AGD"
-                units_text = "(mGy)"
-                name_text = "Acquisition protocol"
-                variable_name_start = "acquisition"
-                variable_value_name = "AGD"
-                modality_text = "CT"
-                chart_message = ""
-
-                new_charts = generate_average_chart_group(
-                    average_choices,
-                    chart_message,
-                    df,
-                    modality_text,
-                    name_field,
-                    name_text,
-                    return_as_dict,
-                    return_structure,
-                    units_text,
-                    user_profile,
-                    value_field,
-                    value_text,
-                    variable_name_start,
-                    variable_value_name,
-                    user_profile.plotMGInitialSortingChoice,
-                )
-
-                return_structure = {**return_structure, **new_charts}
-
-        if user_profile.plotMGAGDvsThickness:
-            parameter_dict = {  # pylint: disable=line-too-long
-                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
-                "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
-                "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose",  # pylint: disable=line-too-long
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotMGInitialSortingChoice,
-                ],
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "legend_title": "Acquisition protocol",
-                "colourmap": user_profile.plotColourMapChoice,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "x_axis_title": "Compressed breast thickness (mm)",
-                "y_axis_title": "AGD (mGy)",
-                "filename": "OpenREM CT acquisition protocol AGD vs thickness",
-                "return_as_dict": return_as_dict,
-            }
-            return_structure[
-                "AGDvsThickness"
-            ] = plotly_scatter(  # pylint: disable=line-too-long
-                df,
-                parameter_dict,
-            )
-
-        if user_profile.plotMGkVpvsThickness:
-            parameter_dict = {  # pylint: disable=line-too-long
-                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
-                "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
-                "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp",  # pylint: disable=line-too-long
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotMGInitialSortingChoice,
-                ],
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "legend_title": "Acquisition protocol",
-                "colourmap": user_profile.plotColourMapChoice,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "x_axis_title": "Compressed breast thickness (mm)",
-                "y_axis_title": "kVp",
-                "filename": "OpenREM CT acquisition protocol kVp vs thickness",
-                "return_as_dict": return_as_dict,
-            }
-            return_structure[
-                "kVpvsThickness"
-            ] = plotly_scatter(  # pylint: disable=line-too-long
-                df,
-                parameter_dict,
-            )
-
-        if user_profile.plotMGmAsvsThickness:
-            parameter_dict = {  # pylint: disable=line-too-long
-                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
-                "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
-                "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure",  # pylint: disable=line-too-long
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotMGInitialSortingChoice,
-                ],
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "legend_title": "Acquisition protocol",
-                "colourmap": user_profile.plotColourMapChoice,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "x_axis_title": "Compressed breast thickness (mm)",
-                "y_axis_title": "mAs",
-                "filename": "OpenREM CT acquisition protocol mAs vs thickness",
-                "return_as_dict": return_as_dict,
-            }
-            return_structure[
-                "mAsvsThickness"
-            ] = plotly_scatter(  # pylint: disable=line-too-long
-                df,
-                parameter_dict,
-            )
-
-        if user_profile.plotMGacquisitionFreq:
-            parameter_dict = {
-                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotMGInitialSortingChoice,
-                ],
-                "legend_title": "Acquisition",
-                "df_x_axis_col": "x_ray_system_name",
-                "x_axis_title": "System",
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "colourmap": user_profile.plotColourMapChoice,
-                "filename": "OpenREM MG acquisition protocol frequency",
-                "groupby_cols": None,
-                "facet_col": None,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "return_as_dict": return_as_dict,
-            }
-            (
-                return_structure["acquisitionFrequencyData"],
-                return_structure["acquisitionFrequencyDataCSV"],
-            ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
-                df, parameter_dict, csv_name="acquisitionFrequencyData.csv"
-            )
-
-        if user_profile.plotMGAcquisitionAGDOverTime:
-            facet_title = "System"
-
-            if user_profile.plotGroupingChoice == "series":
-                facet_title = "Acquisition protocol"
-
-            parameter_dict = {  # pylint: disable=line-too-long
-                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
-                "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose",  # pylint: disable=line-too-long
-                "df_date_col": "study_date",
-                "name_title": "Acquisition protocol",
-                "value_title": "AGD (mGy)",
-                "date_title": "Study date",
-                "facet_title": facet_title,
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotMGInitialSortingChoice,
-                ],
-                "time_period": plot_timeunit_period,
-                "average_choices": average_choices + ["count"],
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "colourmap": user_profile.plotColourMapChoice,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "filename": "OpenREM MG acquisition protocol AGD over time",
-                "return_as_dict": return_as_dict,
-            }
-            result = construct_over_time_charts(
-                df,
-                parameter_dict,
-            )
-
-            if user_profile.plotMean:
-                return_structure["acquisitionMeanAGDOverTime"] = result["mean"]
-            if user_profile.plotMedian:
-                return_structure["acquisitionMedianAGDOverTime"] = result["median"]
-
-    if enable_standard_names:
-        charts_of_interest = [
-            user_profile.plotMGStandardAGDvsThickness,
-            user_profile.plotMGStandardkVpvsThickness,
-            user_profile.plotMGStandardmAsvsThickness,
-            user_profile.plotMGStandardAverageAGDvsThickness,
-            user_profile.plotMGStandardAverageAGD,
-            user_profile.plotMGStandardAcquisitionFreq,
-            user_profile.plotMGStandardAcquisitionAGDOverTime,
-        ]
-
-        if any(charts_of_interest):
-            # Exclude "Blank" and "blank" standard_acquisition_name data
-            df_without_blanks = df[
-                (
-                    df[
-                        "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-                    ]
-                    != "blank"
-                )
-                & (
-                    df[
-                        "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-                    ]
-                    != "Blank"
-                )
-            ].copy()
-            # Remove any unused categories (this will include "Blank" or "blank")
-            df_without_blanks[
-                "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-            ] = df_without_blanks[
-                "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-            ].cat.remove_unused_categories()
-
-            if (
-                user_profile.plotMGStandardAverageAGDvsThickness
-                or user_profile.plotMGStandardAverageAGD
-            ):
-
-                if user_profile.plotBoxplots and "median" not in average_choices:
-                    average_choices = average_choices + ["median"]
-
-                name_field = "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-                value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose"  # pylint: disable=line-too-long
-
-                if user_profile.plotMGStandardAverageAGDvsThickness:
-                    category_names_col = name_field
-                    group_by_col = "x_ray_system_name"
-                    legend_title = "Standard acquisition name"
-
-                    if user_profile.plotGroupingChoice == "series":
-                        category_names_col = "x_ray_system_name"
-                        group_by_col = name_field
-                        legend_title = "System"
-
-                    parameter_dict = {  # pylint: disable=line-too-long
-                        "df_x_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
-                        "df_y_value_col": value_field,  # pylint: disable=line-too-long
-                        "x_axis_title": "Compressed breast thickness (mm)",
-                        "y_axis_title": "AGD (mGy)",
-                        "df_category_col": category_names_col,
-                        "df_facet_col": group_by_col,
-                        "facet_title": legend_title,
-                        "user_bins": [20, 30, 40, 50, 60, 70, 80, 90],
-                        "colourmap": user_profile.plotColourMapChoice,
-                        "filename": "OpenREM CT standard acquisition name average AGD vs thickness",
-                        "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                        "sorting_choice": [
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotMGInitialSortingChoice,
-                        ],
-                        "return_as_dict": return_as_dict,
-                    }
-                    if user_profile.plotMean:
-                        parameter_dict["stat_name"] = "mean"
-                        return_structure[
-                            "standardMeanAGDvsThickness"
-                        ] = plotly_binned_statistic_barchart(
-                            df_without_blanks,
-                            parameter_dict,
-                        )
-
-                    if user_profile.plotMedian:
-                        parameter_dict["stat_name"] = "median"
-                        return_structure[
-                            "standardMedianAGDvsThickness"
-                        ] = plotly_binned_statistic_barchart(
-                            df_without_blanks,
-                            parameter_dict,
-                        )
-
-                if user_profile.plotMGStandardAverageAGD:
-                    value_text = "AGD"
-                    units_text = "(mGy)"
-                    name_text = "Standard acquisition name"
-                    variable_name_start = "standardAcquisition"
-                    variable_value_name = "AGD"
-                    modality_text = "CT"
-                    chart_message = ""
-
-                    new_charts = generate_average_chart_group(
-                        average_choices,
-                        chart_message,
-                        df_without_blanks,
-                        modality_text,
-                        name_field,
-                        name_text,
-                        return_as_dict,
-                        return_structure,
-                        units_text,
-                        user_profile,
-                        value_field,
-                        value_text,
-                        variable_name_start,
-                        variable_value_name,
-                        user_profile.plotMGInitialSortingChoice,
-                    )
-
-                    return_structure = {**return_structure, **new_charts}
-
-            if user_profile.plotMGStandardAGDvsThickness:
-                parameter_dict = {  # pylint: disable=line-too-long
-                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
-                    "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
-                    "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose",  # pylint: disable=line-too-long
-                    "sorting_choice": [
-                        user_profile.plotInitialSortingDirection,
-                        user_profile.plotMGInitialSortingChoice,
-                    ],
-                    "grouping_choice": user_profile.plotGroupingChoice,
-                    "legend_title": "Standard acquisition name",
-                    "colourmap": user_profile.plotColourMapChoice,
-                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                    "x_axis_title": "Compressed breast thickness (mm)",
-                    "y_axis_title": "AGD (mGy)",
-                    "filename": "OpenREM CT standard acquisition name AGD vs thickness",
-                    "return_as_dict": return_as_dict,
-                }
-                return_structure[
-                    "standardAGDvsThickness"
-                ] = plotly_scatter(  # pylint: disable=line-too-long
-                    df_without_blanks,
-                    parameter_dict,
-                )
-
-            if user_profile.plotMGStandardkVpvsThickness:
-                parameter_dict = {  # pylint: disable=line-too-long
-                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
-                    "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
-                    "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp",  # pylint: disable=line-too-long
-                    "sorting_choice": [
-                        user_profile.plotInitialSortingDirection,
-                        user_profile.plotMGInitialSortingChoice,
-                    ],
-                    "grouping_choice": user_profile.plotGroupingChoice,
-                    "legend_title": "Standard acquisition name",
-                    "colourmap": user_profile.plotColourMapChoice,
-                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                    "x_axis_title": "Compressed breast thickness (mm)",
-                    "y_axis_title": "kVp",
-                    "filename": "OpenREM CT standard acquisition name kVp vs thickness",
-                    "return_as_dict": return_as_dict,
-                }
-                return_structure[
-                    "standardkVpvsThickness"
-                ] = plotly_scatter(  # pylint: disable=line-too-long
-                    df_without_blanks,
-                    parameter_dict,
-                )
-
-            if user_profile.plotMGStandardmAsvsThickness:
-                parameter_dict = {  # pylint: disable=line-too-long
-                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
-                    "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
-                    "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure",  # pylint: disable=line-too-long
-                    "sorting_choice": [
-                        user_profile.plotInitialSortingDirection,
-                        user_profile.plotMGInitialSortingChoice,
-                    ],
-                    "grouping_choice": user_profile.plotGroupingChoice,
-                    "legend_title": "Standard acquisition name",
-                    "colourmap": user_profile.plotColourMapChoice,
-                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                    "x_axis_title": "Compressed breast thickness (mm)",
-                    "y_axis_title": "mAs",
-                    "filename": "OpenREM CT standard acquisition name mAs vs thickness",
-                    "return_as_dict": return_as_dict,
-                }
-                return_structure[
-                    "standardmAsvsThickness"
-                ] = plotly_scatter(  # pylint: disable=line-too-long
-                    df_without_blanks,
-                    parameter_dict,
-                )
-
-            if user_profile.plotMGStandardAcquisitionFreq:
-                parameter_dict = {
-                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
-                    "sorting_choice": [
-                        user_profile.plotInitialSortingDirection,
-                        user_profile.plotMGInitialSortingChoice,
-                    ],
-                    "legend_title": "Acquisition",
-                    "df_x_axis_col": "x_ray_system_name",
-                    "x_axis_title": "System",
-                    "grouping_choice": user_profile.plotGroupingChoice,
-                    "colourmap": user_profile.plotColourMapChoice,
-                    "filename": "OpenREM MG standard acquisition name frequency",
-                    "groupby_cols": None,
-                    "facet_col": None,
-                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                    "return_as_dict": return_as_dict,
-                }
-                (
-                    return_structure["standardAcquisitionFrequencyData"],
-                    return_structure["standardAcquisitionFrequencyDataCSV"],
-                ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
-                    df_without_blanks,
-                    parameter_dict,
-                    csv_name="standardAcquisitionFrequencyData.csv",
-                )
-
-            if user_profile.plotMGStandardAcquisitionAGDOverTime:
-                facet_title = "System"
-
-                if user_profile.plotGroupingChoice == "series":
-                    facet_title = "Standard acquisition name"
-
-                parameter_dict = {  # pylint: disable=line-too-long
-                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
-                    "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose",  # pylint: disable=line-too-long
-                    "df_date_col": "study_date",
-                    "name_title": "Standard acquisition name",
-                    "value_title": "AGD (mGy)",
-                    "date_title": "Study date",
-                    "facet_title": facet_title,
-                    "sorting_choice": [
-                        user_profile.plotInitialSortingDirection,
-                        user_profile.plotMGInitialSortingChoice,
-                    ],
-                    "time_period": plot_timeunit_period,
-                    "average_choices": average_choices + ["count"],
-                    "grouping_choice": user_profile.plotGroupingChoice,
-                    "colourmap": user_profile.plotColourMapChoice,
-                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                    "filename": "OpenREM MG standard acquisition name AGD over time",
-                    "return_as_dict": return_as_dict,
-                }
-                result = construct_over_time_charts(
-                    df_without_blanks,
-                    parameter_dict,
-                )
-
-                if user_profile.plotMean:
-                    return_structure["standardAcquisitionMeanAGDOverTime"] = result[
-                        "mean"
-                    ]
-                if user_profile.plotMedian:
-                    return_structure["standardAcquisitionMedianAGDOverTime"] = result[
-                        "median"
-                    ]
-
-    #######################################################################
-    # Prepare study- and request-level Pandas DataFrame to use for charts
-    charts_of_interest = [user_profile.plotMGStudyPerDayAndHour]
-    if enable_standard_names:
-        charts_of_interest.append(user_profile.plotMGStandardStudyPerDayAndHour)
-
-    if any(charts_of_interest):
-
-        name_fields = []
-        charts_of_interest = [user_profile.plotMGStudyPerDayAndHour]
-        if any(charts_of_interest):
-            name_fields.append("study_description")
-
-        if enable_standard_names:
-            charts_of_interest = [user_profile.plotMGStandardStudyPerDayAndHour]
-            if any(charts_of_interest):
-                name_fields.append("standard_names__standard_name")
-
-        value_fields = []
-
-        date_fields = []
-        time_fields = []
-        charts_of_interest = [user_profile.plotMGStudyPerDayAndHour]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotMGStandardStudyPerDayAndHour)
-        if any(charts_of_interest):
-            date_fields.append("study_date")
-            time_fields.append("study_time")
-
-        system_field = []
-        if user_profile.plotSeriesPerSystem:
-            system_field.append(
-                "generalequipmentmoduleattr__unique_equipment_name_id__display_name"
-            )
-
-        fields = {
-            "names": name_fields,
-            "values": value_fields,
-            "dates": date_fields,
-            "times": time_fields,
-            "system": system_field,
-        }
-
-        # If only standard_names__standard_name is required then exclude all entries where these are None as these are
-        # not required for standard name charts.
-        queryset = f.qs
-        if name_fields == ["standard_names__standard_name"]:
-            queryset = queryset.exclude(standard_names__standard_name__isnull=True)
-
-        df = create_dataframe(
-            queryset,
-            fields,
-            data_point_name_lowercase=user_profile.plotCaseInsensitiveCategories,
-            data_point_name_remove_whitespace_padding=user_profile.plotRemoveCategoryWhitespacePadding,
-            uid="pk",
-        )
-
-        if user_profile.plotMGStudyPerDayAndHour:
-            df_time_series_per_weekday = create_dataframe_weekdays(
-                df, "study_description", df_date_col="study_date"
-            )
-
-            return_structure["studyWorkloadData"] = plotly_barchart_weekdays(
-                df_time_series_per_weekday,
-                "weekday",
-                "study_description",
-                name_axis_title="Weekday",
-                value_axis_title="Frequency",
-                colourmap=user_profile.plotColourMapChoice,
-                filename="OpenREM CT study description workload",
-                facet_col_wrap=user_profile.plotFacetColWrapVal,
-                sorting_choice=[
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotMGInitialSortingChoice,
-                ],
-                return_as_dict=return_as_dict,
-            )
-
-        if enable_standard_names:
-            charts_of_interest = [user_profile.plotMGStandardStudyPerDayAndHour]
-
-            if any(charts_of_interest):
-
-                # Create a standard name data frame - remove any blank standard names
-                standard_name_df = df[
-                    (df["standard_names__standard_name"] != "blank")
-                    & (df["standard_names__standard_name"] != "Blank")
-                ].copy()
-                # Remove any unused categories (this will include "Blank" or "blank")
-                standard_name_df["standard_names__standard_name"] = standard_name_df[
-                    "standard_names__standard_name"
-                ].cat.remove_unused_categories()
-
-                if user_profile.plotMGStandardStudyPerDayAndHour:
-                    df_time_series_per_weekday = create_dataframe_weekdays(
-                        standard_name_df,
-                        "standard_names__standard_name",
-                        df_date_col="study_date",
-                    )
-
-                    return_structure[
-                        "standardStudyWorkloadData"
-                    ] = plotly_barchart_weekdays(
-                        df_time_series_per_weekday,
-                        "weekday",
-                        "standard_names__standard_name",
-                        name_axis_title="Weekday",
-                        value_axis_title="Frequency",
-                        colourmap=user_profile.plotColourMapChoice,
-                        filename="OpenREM CT standard study name workload",
-                        facet_col_wrap=user_profile.plotFacetColWrapVal,
-                        sorting_choice=[
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotMGInitialSortingChoice,
-                        ],
-                        return_as_dict=return_as_dict,
-                    )
-
-    return return_structure
-
-
-def mg_chart_form_processing(request, user_profile):
-    # pylint: disable=too-many-statements
-
-    # Obtain the system-level enable_standard_names setting
-    try:
-        StandardNameSettings.objects.get()
-    except ObjectDoesNotExist:
-        StandardNameSettings.objects.create()
-    enable_standard_names = StandardNameSettings.objects.values_list(
-        "enable_standard_names", flat=True
-    )[0]
-
-    # Obtain the chart options from the request
-    chart_options_form = None
-    if enable_standard_names:
-        chart_options_form = MGChartOptionsFormIncStandard(request.GET)
-    else:
-        chart_options_form = MGChartOptionsForm(request.GET)
-
-    # Check whether the form data is valid
-    if chart_options_form.is_valid():
-        # Use the form data if the user clicked on the submit button
-        if "submit" in request.GET:
-            # process the data in form.cleaned_data as required
-
-            set_common_chart_options(chart_options_form, user_profile)
-
-            set_average_chart_options(chart_options_form, user_profile)
-
-            set_mg_chart_options(chart_options_form, user_profile)
-
-            user_profile.save()
-
-        else:
-            average_choices = required_average_choices(user_profile)
-
-            mg_form_data = initialise_mg_form_data(user_profile)
-
-            form_data = {
-                "plotCharts": user_profile.plotCharts,
-                "plotGrouping": user_profile.plotGroupingChoice,
-                "plotSeriesPerSystem": user_profile.plotSeriesPerSystem,
-                "plotHistograms": user_profile.plotHistograms,
-                "plotInitialSortingDirection": user_profile.plotInitialSortingDirection,
-                "plotAverageChoice": average_choices,
-            }
-
-            form_data = {**form_data, **mg_form_data}
-
-            chart_options_form = None
-            if enable_standard_names:
-                chart_options_form = MGChartOptionsFormIncStandard(form_data)
-            else:
-                chart_options_form = MGChartOptionsForm(form_data)
-
-    return chart_options_form
+# pylint: disable=too-many-lines
+import logging
+from datetime import datetime
+from django.conf import settings
+from django.contrib.auth.decorators import login_required
+from django.core.exceptions import ObjectDoesNotExist
+from django.http import JsonResponse
+from remapp.forms import (
+    MGChartOptionsForm,
+    MGChartOptionsFormIncStandard,
+)
+from remapp.interface.mod_filters import (
+    MGSummaryListFilter,
+    MGFilterPlusPid,
+    MGFilterPlusStdNames,
+    MGFilterPlusPidPlusStdNames,
+)
+from remapp.models import (
+    GeneralStudyModuleAttr,
+    create_user_profile,
+    StandardNameSettings,
+)
+from remapp.views_admin import (
+    required_average_choices,
+    initialise_mg_form_data,
+    set_average_chart_options,
+    set_mg_chart_options,
+    set_common_chart_options,
+)
+from .interface.chart_functions import (
+    create_dataframe,
+    create_dataframe_weekdays,
+    plotly_barchart_weekdays,
+    plotly_binned_statistic_barchart,
+    plotly_boxplot,
+    plotly_barchart,
+    plotly_histogram_barchart,
+    plotly_scatter,
+    plotly_frequency_barchart,
+    construct_over_time_charts,
+    plotly_set_default_theme,
+    create_dataframe_aggregates,
+    generate_average_chart_group,
+)
+
+from .tools.check_standard_name_status import are_standard_names_enabled
+
+logger = logging.getLogger(__name__)
+
+
+def generate_required_mg_charts_list(profile):
+    # pylint: disable=too-many-branches
+    # pylint: disable=too-many-statements
+    """Obtain a list of dictionaries containing the title string and base
+    variable name for each required chart"""
+
+    # Obtain the system-level enable_standard_names setting
+    enable_standard_names = are_standard_names_enabled()
+
+    required_charts = []
+
+    charts_of_interest = [
+        profile.plotMGAcquisitionAGDOverTime,
+    ]
+    if enable_standard_names:
+        charts_of_interest.append(profile.plotMGStandardAcquisitionAGDOverTime)
+
+    if any(charts_of_interest):
+        keys = list(dict(profile.TIME_PERIOD).keys())
+        values = list(dict(profile.TIME_PERIOD).values())
+        time_period = (values[keys.index(profile.plotMGOverTimePeriod)]).lower()
+
+    if profile.plotMGacquisitionFreq:
+        required_charts.append(
+            {
+                "title": "Chart of acquisition protocol frequency",
+                "var_name": "acquisitionFrequency",
+            }
+        )
+
+    if profile.plotMGaverageAGD:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of acquisition protocol mean AGD",
+                    "var_name": "acquisitionMeanAGD",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of acquisition protocol median AGD",
+                    "var_name": "acquisitionMedianAGD",
+                }
+            )
+        if profile.plotBoxplots:
+            required_charts.append(
+                {
+                    "title": "Boxplot of acquisition protocol AGD",
+                    "var_name": "acquisitionBoxplotAGD",
+                }
+            )
+        if profile.plotHistograms:
+            required_charts.append(
+                {
+                    "title": "Histogram of acquisition protocol AGD",
+                    "var_name": "acquisitionHistogramAGD",
+                }
+            )
+
+    if profile.plotMGAGDvsThickness:
+        required_charts.append(
+            {
+                "title": "Chart of acquisition protocol AGD vs compressed breast thickness",
+                "var_name": "acquisitionScatterAGDvsThick",
+            }
+        )
+
+    if profile.plotMGaverageAGDvsThickness:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of acquisition protocol mean AGD vs compressed breast thickness",
+                    "var_name": "acquisitionMeanAGDvsThick",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of acquisition protocol median AGD vs compressed breast thickness",
+                    "var_name": "acquisitionMedianAGDvsThick",
+                }
+            )
+
+    if profile.plotMGAcquisitionAGDOverTime:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of acquisition protocol mean AGD over time ("
+                    + time_period
+                    + ")",
+                    "var_name": "acquisitionMeanAGDOverTime",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of acquisition protocol median AGD over time ("
+                    + time_period
+                    + ")",
+                    "var_name": "acquisitionMedianAGDOverTime",
+                }
+            )
+
+    if profile.plotMGkVpvsThickness:
+        required_charts.append(
+            {
+                "title": "Chart of acquisition protocol kVp vs compressed breast thickness",
+                "var_name": "acquisitionScatterkVpvsThick",
+            }
+        )
+
+    if profile.plotMGmAsvsThickness:
+        required_charts.append(
+            {
+                "title": "Chart of acquisition protocol mAs vs compressed breast thickness",
+                "var_name": "acquisitionScattermAsvsThick",
+            }
+        )
+
+    if profile.plotMGStudyPerDayAndHour:
+        required_charts.append(
+            {
+                "title": "Chart of study description workload",
+                "var_name": "studyWorkload",
+            }
+        )
+
+    if enable_standard_names:
+        if profile.plotMGStandardAcquisitionFreq:
+            required_charts.append(
+                {
+                    "title": "Chart of standard acquisition name frequency",
+                    "var_name": "standardAcquisitionFrequency",
+                }
+            )
+
+        if profile.plotMGStandardAverageAGD:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of standard acquisition name mean AGD",
+                        "var_name": "standardAcquisitionMeanAGD",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of standard acquisition name median AGD",
+                        "var_name": "standardAcquisitionMedianAGD",
+                    }
+                )
+            if profile.plotBoxplots:
+                required_charts.append(
+                    {
+                        "title": "Boxplot of standard acquisition name AGD",
+                        "var_name": "standardAcquisitionBoxplotAGD",
+                    }
+                )
+            if profile.plotHistograms:
+                required_charts.append(
+                    {
+                        "title": "Histogram of standard acquisition name AGD",
+                        "var_name": "standardAcquisitionHistogramAGD",
+                    }
+                )
+
+        if profile.plotMGStandardAGDvsThickness:
+            required_charts.append(
+                {
+                    "title": "Chart of standard acquisition name AGD vs compressed breast thickness",
+                    "var_name": "standardAcquisitionScatterAGDvsThick",
+                }
+            )
+
+        if profile.plotMGStandardAverageAGDvsThickness:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of standard acquisition name mean AGD vs compressed breast thickness",
+                        "var_name": "standardAcquisitionMeanAGDvsThick",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of standard acquisition name median AGD vs compressed breast thickness",
+                        "var_name": "standardAcquisitionMedianAGDvsThick",
+                    }
+                )
+
+        if profile.plotMGStandardAcquisitionAGDOverTime:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of standard acquisition name mean AGD over time ("
+                        + time_period
+                        + ")",
+                        "var_name": "standardAcquisitionMeanAGDOverTime",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of standard acquisition name median AGD over time ("
+                        + time_period
+                        + ")",
+                        "var_name": "standardAcquisitionMedianAGDOverTime",
+                    }
+                )
+
+        if profile.plotMGStandardkVpvsThickness:
+            required_charts.append(
+                {
+                    "title": "Chart of standard acquisition name kVp vs compressed breast thickness",
+                    "var_name": "standardAcquisitionScatterkVpvsThick",
+                }
+            )
+
+        if profile.plotMGStandardmAsvsThickness:
+            required_charts.append(
+                {
+                    "title": "Chart of standard acquisition name mAs vs compressed breast thickness",
+                    "var_name": "standardAcquisitionScattermAsvsThick",
+                }
+            )
+
+        if profile.plotMGStandardStudyPerDayAndHour:
+            required_charts.append(
+                {
+                    "title": "Chart of standard study name workload",
+                    "var_name": "standardStudyWorkload",
+                }
+            )
+
+    return required_charts
+
+
+@login_required
+def mg_summary_chart_data(request):
+    """Obtain data for mammography chart data Ajax view"""
+
+    # Obtain the system-level enable_standard_names setting
+    enable_standard_names = are_standard_names_enabled()
+
+    if request.user.groups.filter(name="pidgroup"):
+        if enable_standard_names:
+            f = MGFilterPlusPidPlusStdNames(
+                request.GET,
+                queryset=GeneralStudyModuleAttr.objects.filter(
+                    modality_type__exact="MG"
+                )
+                .order_by()
+                .distinct(),
+            )
+        else:
+            f = MGFilterPlusPid(
+                request.GET,
+                queryset=GeneralStudyModuleAttr.objects.filter(
+                    modality_type__exact="MG"
+                )
+                .order_by()
+                .distinct(),
+            )
+    else:
+        if enable_standard_names:
+            f = MGFilterPlusStdNames(
+                request.GET,
+                queryset=GeneralStudyModuleAttr.objects.filter(
+                    modality_type__exact="MG"
+                )
+                .order_by()
+                .distinct(),
+            )
+        else:
+            f = MGSummaryListFilter(
+                request.GET,
+                queryset=GeneralStudyModuleAttr.objects.filter(
+                    modality_type__exact="MG"
+                )
+                .order_by()
+                .distinct(),
+            )
+    try:
+        # See if the user has plot settings in userprofile
+        user_profile = request.user.userprofile
+    except ObjectDoesNotExist:
+        # Create a default userprofile for the user if one doesn't exist
+        create_user_profile(sender=request.user, instance=request.user, created=True)
+        user_profile = request.user.userprofile
+
+    if settings.DEBUG:
+        start_time = datetime.now()
+
+    return_structure = mg_plot_calculations(f, user_profile)
+
+    if settings.DEBUG:
+        logger.debug(f"Elapsed time is {datetime.now() - start_time}")
+
+    return JsonResponse(return_structure, safe=False)
+
+
+def mg_plot_calculations(f, user_profile, return_as_dict=False):
+    # pylint: disable=too-many-locals
+    # pylint: disable=too-many-branches
+    # pylint: disable=too-many-statements
+    """Calculations for mammography charts"""
+    # Return an empty structure if the queryset is empty
+    if not f.qs.exists():
+        return {}
+
+    # Obtain the system-level enable_standard_names setting
+    enable_standard_names = are_standard_names_enabled()
+
+    # Set the Plotly chart theme
+    plotly_set_default_theme(user_profile.plotThemeChoice)
+
+    return_structure = {}
+
+    average_choices = []
+    if user_profile.plotMean:
+        average_choices.append("mean")
+    if user_profile.plotMedian:
+        average_choices.append("median")
+
+    charts_of_interest = [user_profile.plotMGAcquisitionAGDOverTime]
+    if enable_standard_names:
+        charts_of_interest.append(user_profile.plotMGStandardAcquisitionAGDOverTime)
+    if any(charts_of_interest):
+        plot_timeunit_period = user_profile.plotMGOverTimePeriod
+
+    #######################################################################
+    # Prepare acquisition-level Pandas DataFrame to use for charts
+    charts_of_interest = [
+        user_profile.plotMGAGDvsThickness,
+        user_profile.plotMGkVpvsThickness,
+        user_profile.plotMGmAsvsThickness,
+        user_profile.plotMGaverageAGDvsThickness,
+        user_profile.plotMGaverageAGD,
+        user_profile.plotMGacquisitionFreq,
+        user_profile.plotMGAcquisitionAGDOverTime,
+    ]
+    if enable_standard_names:
+        charts_of_interest.append(user_profile.plotMGStandardAGDvsThickness)
+        charts_of_interest.append(user_profile.plotMGStandardkVpvsThickness)
+        charts_of_interest.append(user_profile.plotMGStandardmAsvsThickness)
+        charts_of_interest.append(user_profile.plotMGStandardAverageAGDvsThickness)
+        charts_of_interest.append(user_profile.plotMGStandardAverageAGD)
+        charts_of_interest.append(user_profile.plotMGStandardAcquisitionFreq)
+        charts_of_interest.append(user_profile.plotMGStandardAcquisitionAGDOverTime)
+
+    if any(charts_of_interest):  # pylint: disable=too-many-nested-blocks
+
+        name_fields = []
+        charts_of_interest = [
+            user_profile.plotMGAGDvsThickness,
+            user_profile.plotMGkVpvsThickness,
+            user_profile.plotMGmAsvsThickness,
+            user_profile.plotMGaverageAGDvsThickness,
+            user_profile.plotMGaverageAGD,
+            user_profile.plotMGacquisitionFreq,
+            user_profile.plotMGAcquisitionAGDOverTime,
+        ]
+        if any(charts_of_interest):
+            name_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+            )
+            name_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure_control_mode"
+            )
+            name_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__image_view__code_meaning"
+            )
+
+        if enable_standard_names:
+            charts_of_interest = [
+                user_profile.plotMGStandardAGDvsThickness,
+                user_profile.plotMGStandardkVpvsThickness,
+                user_profile.plotMGStandardmAsvsThickness,
+                user_profile.plotMGStandardAverageAGDvsThickness,
+                user_profile.plotMGStandardAverageAGD,
+                user_profile.plotMGStandardAcquisitionFreq,
+                user_profile.plotMGStandardAcquisitionAGDOverTime,
+            ]
+            if any(charts_of_interest):
+                name_fields.append(
+                    "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                )
+
+        value_fields = []
+        value_multipliers = []
+        charts_of_interest = [
+            user_profile.plotMGAGDvsThickness,
+            user_profile.plotMGaverageAGDvsThickness,
+            user_profile.plotMGaverageAGD,
+            user_profile.plotMGAcquisitionAGDOverTime,
+        ]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotMGStandardAGDvsThickness)
+            charts_of_interest.append(user_profile.plotMGStandardAverageAGDvsThickness)
+            charts_of_interest.append(user_profile.plotMGStandardAverageAGD)
+            charts_of_interest.append(user_profile.plotMGStandardAcquisitionAGDOverTime)
+        if any(charts_of_interest):
+            value_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose"
+            )
+            value_multipliers.append(1)
+
+        charts_of_interest = [user_profile.plotMGkVpvsThickness]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotMGStandardkVpvsThickness)
+        if any(charts_of_interest):
+            value_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp"
+            )
+            value_multipliers.append(1)
+
+        charts_of_interest = [user_profile.plotMGmAsvsThickness]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotMGStandardmAsvsThickness)
+        if any(charts_of_interest):
+            value_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure"
+            )
+            value_multipliers.append(0.001)
+
+        charts_of_interest = [
+            user_profile.plotMGAGDvsThickness,
+            user_profile.plotMGkVpvsThickness,
+            user_profile.plotMGmAsvsThickness,
+            user_profile.plotMGaverageAGDvsThickness,
+        ]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotMGStandardAGDvsThickness)
+            charts_of_interest.append(user_profile.plotMGStandardkVpvsThickness)
+            charts_of_interest.append(user_profile.plotMGStandardmAsvsThickness)
+            charts_of_interest.append(user_profile.plotMGStandardAverageAGDvsThickness)
+        if any(charts_of_interest):
+            value_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness"
+            )
+            value_multipliers.append(1)
+
+        date_fields = []
+        charts_of_interest = [user_profile.plotMGAcquisitionAGDOverTime]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotMGStandardAcquisitionAGDOverTime)
+        if any(charts_of_interest):
+            date_fields.append("study_date")
+
+        time_fields = []
+
+        system_field = []
+        if user_profile.plotSeriesPerSystem:
+            system_field.append(
+                "generalequipmentmoduleattr__unique_equipment_name_id__display_name"
+            )
+
+        fields = {
+            "names": name_fields,
+            "values": value_fields,
+            "dates": date_fields,
+            "times": time_fields,
+            "system": system_field,
+        }
+
+        df = create_dataframe(
+            f.qs,
+            fields,
+            data_point_name_lowercase=user_profile.plotCaseInsensitiveCategories,
+            data_point_name_remove_whitespace_padding=user_profile.plotRemoveCategoryWhitespacePadding,
+            data_point_value_multipliers=value_multipliers,
+            char_wrap=user_profile.plotLabelCharWrap,
+            uid="projectionxrayradiationdose__irradeventxraydata__pk",
+        )
+
+        #######################################################################
+        # Addressing issue 1050 (https://bitbucket.org/openrem/openrem/issues/1050)
+        # Filter data frame using acquisition protocol name,  min / max compression thickness,
+        # standard acquisition name, image view code, and exposure control mode filters.
+        if (
+            "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness__range_min"
+            in f.data
+        ):
+            if f.data[
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness__range_min"
+            ]:
+                df = df[
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness"
+                    ]
+                    >= float(
+                        f.data[
+                            "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness__range_min"
+                        ]
+                    )
+                ]
+
+        if (
+            "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness__range_max"
+            in f.data
+        ):
+            if f.data[
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness__range_max"
+            ]:
+                df = df[
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness"
+                    ]
+                    <= float(
+                        f.data[
+                            "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness__range_max"
+                        ]
+                    )
+                ]
+
+        if (
+            "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+            in f.data
+        ):
+            if f.data[
+                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+            ]:
+                df[
+                    "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+                ] = df[
+                    "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+                ].astype(
+                    "object"
+                )
+                df = df[
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+                    ].str.contains(
+                        f.data[
+                            "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+                        ],
+                        case=False,
+                    )
+                ]
+                df[
+                    "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+                ] = df[
+                    "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+                ].astype(
+                    "category"
+                )
+
+        if (
+            "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+            in f.data
+        ):
+            if f.data[
+                "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+            ]:
+                df[
+                    "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                ] = df[
+                    "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                ].astype(
+                    "object"
+                )
+                df = df[
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                    ].str.contains(
+                        f.data[
+                            "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                        ],
+                        case=False,
+                    )
+                ]
+                df[
+                    "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                ] = df[
+                    "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                ].astype(
+                    "category"
+                )
+
+        if (
+            "projectionxrayradiationdose__irradeventxraydata__image_view__code_meaning"
+            in f.data
+        ):
+            if f.data[
+                "projectionxrayradiationdose__irradeventxraydata__image_view__code_meaning"
+            ]:
+                df[
+                    "projectionxrayradiationdose__irradeventxraydata__image_view__code_meaning"
+                ] = df[
+                    "projectionxrayradiationdose__irradeventxraydata__image_view__code_meaning"
+                ].astype(
+                    "object"
+                )
+                df = df[
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__image_view__code_meaning"
+                    ].str.contains(
+                        f.data[
+                            "projectionxrayradiationdose__irradeventxraydata__image_view__code_meaning"
+                        ],
+                        case=False,
+                    )
+                ]
+                df[
+                    "projectionxrayradiationdose__irradeventxraydata__image_view__code_meaning"
+                ] = df[
+                    "projectionxrayradiationdose__irradeventxraydata__image_view__code_meaning"
+                ].astype(
+                    "category"
+                )
+
+        if (
+            "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure_control_mode"
+            in f.data
+        ):
+            if f.data[
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure_control_mode"
+            ]:
+                df[
+                    "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure_control_mode"
+                ] = df[
+                    "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure_control_mode"
+                ].astype(
+                    "object"
+                )
+                df = df[
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure_control_mode"
+                    ].str.contains(
+                        f.data[
+                            "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure_control_mode"
+                        ],
+                        case=False,
+                    )
+                ]
+                df[
+                    "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure_control_mode"
+                ] = df[
+                    "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure_control_mode"
+                ].astype(
+                    "category"
+                )
+        #######################################################################
+
+        if user_profile.plotMGaverageAGDvsThickness or user_profile.plotMGaverageAGD:
+
+            if user_profile.plotBoxplots and "median" not in average_choices:
+                average_choices = average_choices + ["median"]
+
+            name_field = (
+                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+            )
+            value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose"  # pylint: disable=line-too-long
+
+            if user_profile.plotMGaverageAGDvsThickness:
+                category_names_col = name_field
+                group_by_col = "x_ray_system_name"
+                legend_title = "Acquisition protocol"
+
+                if user_profile.plotGroupingChoice == "series":
+                    category_names_col = "x_ray_system_name"
+                    group_by_col = name_field
+                    legend_title = "System"
+
+                parameter_dict = {  # pylint: disable=line-too-long
+                    "df_x_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
+                    "df_y_value_col": value_field,  # pylint: disable=line-too-long
+                    "x_axis_title": "Compressed breast thickness (mm)",
+                    "y_axis_title": "AGD (mGy)",
+                    "df_category_col": category_names_col,
+                    "df_facet_col": group_by_col,
+                    "facet_title": legend_title,
+                    "user_bins": [20, 30, 40, 50, 60, 70, 80, 90],
+                    "colourmap": user_profile.plotColourMapChoice,
+                    "filename": "OpenREM CT acquisition protocol average AGD vs thickness",
+                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                    "sorting_choice": [
+                        user_profile.plotInitialSortingDirection,
+                        user_profile.plotMGInitialSortingChoice,
+                    ],
+                    "return_as_dict": return_as_dict,
+                }
+                if user_profile.plotMean:
+                    parameter_dict["stat_name"] = "mean"
+                    return_structure["meanAGDvsThickness"] = (
+                        plotly_binned_statistic_barchart(
+                            df,
+                            parameter_dict,
+                        )
+                    )
+
+                if user_profile.plotMedian:
+                    parameter_dict["stat_name"] = "median"
+                    return_structure["medianAGDvsThickness"] = (
+                        plotly_binned_statistic_barchart(
+                            df,
+                            parameter_dict,
+                        )
+                    )
+
+            if user_profile.plotMGaverageAGD:
+                value_text = "AGD"
+                units_text = "(mGy)"
+                name_text = "Acquisition protocol"
+                variable_name_start = "acquisition"
+                variable_value_name = "AGD"
+                modality_text = "CT"
+                chart_message = ""
+
+                new_charts = generate_average_chart_group(
+                    average_choices,
+                    chart_message,
+                    df,
+                    modality_text,
+                    name_field,
+                    name_text,
+                    return_as_dict,
+                    return_structure,
+                    units_text,
+                    user_profile,
+                    value_field,
+                    value_text,
+                    variable_name_start,
+                    variable_value_name,
+                    user_profile.plotMGInitialSortingChoice,
+                )
+
+                return_structure = {**return_structure, **new_charts}
+
+        if user_profile.plotMGAGDvsThickness:
+            parameter_dict = {  # pylint: disable=line-too-long
+                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
+                "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
+                "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose",  # pylint: disable=line-too-long
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotMGInitialSortingChoice,
+                ],
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "legend_title": "Acquisition protocol",
+                "colourmap": user_profile.plotColourMapChoice,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "x_axis_title": "Compressed breast thickness (mm)",
+                "y_axis_title": "AGD (mGy)",
+                "filename": "OpenREM CT acquisition protocol AGD vs thickness",
+                "return_as_dict": return_as_dict,
+            }
+            return_structure["AGDvsThickness"] = (
+                plotly_scatter(  # pylint: disable=line-too-long
+                    df,
+                    parameter_dict,
+                )
+            )
+
+        if user_profile.plotMGkVpvsThickness:
+            parameter_dict = {  # pylint: disable=line-too-long
+                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
+                "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
+                "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp",  # pylint: disable=line-too-long
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotMGInitialSortingChoice,
+                ],
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "legend_title": "Acquisition protocol",
+                "colourmap": user_profile.plotColourMapChoice,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "x_axis_title": "Compressed breast thickness (mm)",
+                "y_axis_title": "kVp",
+                "filename": "OpenREM CT acquisition protocol kVp vs thickness",
+                "return_as_dict": return_as_dict,
+            }
+            return_structure["kVpvsThickness"] = (
+                plotly_scatter(  # pylint: disable=line-too-long
+                    df,
+                    parameter_dict,
+                )
+            )
+
+        if user_profile.plotMGmAsvsThickness:
+            parameter_dict = {  # pylint: disable=line-too-long
+                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
+                "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
+                "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure",  # pylint: disable=line-too-long
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotMGInitialSortingChoice,
+                ],
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "legend_title": "Acquisition protocol",
+                "colourmap": user_profile.plotColourMapChoice,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "x_axis_title": "Compressed breast thickness (mm)",
+                "y_axis_title": "mAs",
+                "filename": "OpenREM CT acquisition protocol mAs vs thickness",
+                "return_as_dict": return_as_dict,
+            }
+            return_structure["mAsvsThickness"] = (
+                plotly_scatter(  # pylint: disable=line-too-long
+                    df,
+                    parameter_dict,
+                )
+            )
+
+        if user_profile.plotMGacquisitionFreq:
+            parameter_dict = {
+                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotMGInitialSortingChoice,
+                ],
+                "legend_title": "Acquisition",
+                "df_x_axis_col": "x_ray_system_name",
+                "x_axis_title": "System",
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "colourmap": user_profile.plotColourMapChoice,
+                "filename": "OpenREM MG acquisition protocol frequency",
+                "groupby_cols": None,
+                "facet_col": None,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "return_as_dict": return_as_dict,
+            }
+            (
+                return_structure["acquisitionFrequencyData"],
+                return_structure["acquisitionFrequencyDataCSV"],
+            ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
+                df, parameter_dict, csv_name="acquisitionFrequencyData.csv"
+            )
+
+        if user_profile.plotMGAcquisitionAGDOverTime:
+            facet_title = "System"
+
+            if user_profile.plotGroupingChoice == "series":
+                facet_title = "Acquisition protocol"
+
+            parameter_dict = {  # pylint: disable=line-too-long
+                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
+                "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose",  # pylint: disable=line-too-long
+                "df_date_col": "study_date",
+                "name_title": "Acquisition protocol",
+                "value_title": "AGD (mGy)",
+                "date_title": "Study date",
+                "facet_title": facet_title,
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotMGInitialSortingChoice,
+                ],
+                "time_period": plot_timeunit_period,
+                "average_choices": average_choices + ["count"],
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "colourmap": user_profile.plotColourMapChoice,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "filename": "OpenREM MG acquisition protocol AGD over time",
+                "return_as_dict": return_as_dict,
+            }
+            result = construct_over_time_charts(
+                df,
+                parameter_dict,
+            )
+
+            if user_profile.plotMean:
+                return_structure["acquisitionMeanAGDOverTime"] = result["mean"]
+            if user_profile.plotMedian:
+                return_structure["acquisitionMedianAGDOverTime"] = result["median"]
+
+    if enable_standard_names:
+        charts_of_interest = [
+            user_profile.plotMGStandardAGDvsThickness,
+            user_profile.plotMGStandardkVpvsThickness,
+            user_profile.plotMGStandardmAsvsThickness,
+            user_profile.plotMGStandardAverageAGDvsThickness,
+            user_profile.plotMGStandardAverageAGD,
+            user_profile.plotMGStandardAcquisitionFreq,
+            user_profile.plotMGStandardAcquisitionAGDOverTime,
+        ]
+
+        if any(charts_of_interest):
+            # Exclude "Blank" and "blank" standard_acquisition_name data
+            df_without_blanks = df[
+                (
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                    ]
+                    != "blank"
+                )
+                & (
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                    ]
+                    != "Blank"
+                )
+            ].copy()
+            # Remove any unused categories (this will include "Blank" or "blank")
+            df_without_blanks[
+                "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+            ] = df_without_blanks[
+                "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+            ].cat.remove_unused_categories()
+
+            if (
+                user_profile.plotMGStandardAverageAGDvsThickness
+                or user_profile.plotMGStandardAverageAGD
+            ):
+
+                if user_profile.plotBoxplots and "median" not in average_choices:
+                    average_choices = average_choices + ["median"]
+
+                name_field = "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose"  # pylint: disable=line-too-long
+
+                if user_profile.plotMGStandardAverageAGDvsThickness:
+                    category_names_col = name_field
+                    group_by_col = "x_ray_system_name"
+                    legend_title = "Standard acquisition name"
+
+                    if user_profile.plotGroupingChoice == "series":
+                        category_names_col = "x_ray_system_name"
+                        group_by_col = name_field
+                        legend_title = "System"
+
+                    parameter_dict = {  # pylint: disable=line-too-long
+                        "df_x_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
+                        "df_y_value_col": value_field,  # pylint: disable=line-too-long
+                        "x_axis_title": "Compressed breast thickness (mm)",
+                        "y_axis_title": "AGD (mGy)",
+                        "df_category_col": category_names_col,
+                        "df_facet_col": group_by_col,
+                        "facet_title": legend_title,
+                        "user_bins": [20, 30, 40, 50, 60, 70, 80, 90],
+                        "colourmap": user_profile.plotColourMapChoice,
+                        "filename": "OpenREM CT standard acquisition name average AGD vs thickness",
+                        "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                        "sorting_choice": [
+                            user_profile.plotInitialSortingDirection,
+                            user_profile.plotMGInitialSortingChoice,
+                        ],
+                        "return_as_dict": return_as_dict,
+                    }
+                    if user_profile.plotMean:
+                        parameter_dict["stat_name"] = "mean"
+                        return_structure["standardMeanAGDvsThickness"] = (
+                            plotly_binned_statistic_barchart(
+                                df_without_blanks,
+                                parameter_dict,
+                            )
+                        )
+
+                    if user_profile.plotMedian:
+                        parameter_dict["stat_name"] = "median"
+                        return_structure["standardMedianAGDvsThickness"] = (
+                            plotly_binned_statistic_barchart(
+                                df_without_blanks,
+                                parameter_dict,
+                            )
+                        )
+
+                if user_profile.plotMGStandardAverageAGD:
+                    value_text = "AGD"
+                    units_text = "(mGy)"
+                    name_text = "Standard acquisition name"
+                    variable_name_start = "standardAcquisition"
+                    variable_value_name = "AGD"
+                    modality_text = "CT"
+                    chart_message = ""
+
+                    new_charts = generate_average_chart_group(
+                        average_choices,
+                        chart_message,
+                        df_without_blanks,
+                        modality_text,
+                        name_field,
+                        name_text,
+                        return_as_dict,
+                        return_structure,
+                        units_text,
+                        user_profile,
+                        value_field,
+                        value_text,
+                        variable_name_start,
+                        variable_value_name,
+                        user_profile.plotMGInitialSortingChoice,
+                    )
+
+                    return_structure = {**return_structure, **new_charts}
+
+            if user_profile.plotMGStandardAGDvsThickness:
+                parameter_dict = {  # pylint: disable=line-too-long
+                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
+                    "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
+                    "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose",  # pylint: disable=line-too-long
+                    "sorting_choice": [
+                        user_profile.plotInitialSortingDirection,
+                        user_profile.plotMGInitialSortingChoice,
+                    ],
+                    "grouping_choice": user_profile.plotGroupingChoice,
+                    "legend_title": "Standard acquisition name",
+                    "colourmap": user_profile.plotColourMapChoice,
+                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                    "x_axis_title": "Compressed breast thickness (mm)",
+                    "y_axis_title": "AGD (mGy)",
+                    "filename": "OpenREM CT standard acquisition name AGD vs thickness",
+                    "return_as_dict": return_as_dict,
+                }
+                return_structure["standardAGDvsThickness"] = (
+                    plotly_scatter(  # pylint: disable=line-too-long
+                        df_without_blanks,
+                        parameter_dict,
+                    )
+                )
+
+            if user_profile.plotMGStandardkVpvsThickness:
+                parameter_dict = {  # pylint: disable=line-too-long
+                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
+                    "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
+                    "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp",  # pylint: disable=line-too-long
+                    "sorting_choice": [
+                        user_profile.plotInitialSortingDirection,
+                        user_profile.plotMGInitialSortingChoice,
+                    ],
+                    "grouping_choice": user_profile.plotGroupingChoice,
+                    "legend_title": "Standard acquisition name",
+                    "colourmap": user_profile.plotColourMapChoice,
+                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                    "x_axis_title": "Compressed breast thickness (mm)",
+                    "y_axis_title": "kVp",
+                    "filename": "OpenREM CT standard acquisition name kVp vs thickness",
+                    "return_as_dict": return_as_dict,
+                }
+                return_structure["standardkVpvsThickness"] = (
+                    plotly_scatter(  # pylint: disable=line-too-long
+                        df_without_blanks,
+                        parameter_dict,
+                    )
+                )
+
+            if user_profile.plotMGStandardmAsvsThickness:
+                parameter_dict = {  # pylint: disable=line-too-long
+                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
+                    "df_x_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraymechanicaldata__compression_thickness",  # pylint: disable=line-too-long
+                    "df_y_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure",  # pylint: disable=line-too-long
+                    "sorting_choice": [
+                        user_profile.plotInitialSortingDirection,
+                        user_profile.plotMGInitialSortingChoice,
+                    ],
+                    "grouping_choice": user_profile.plotGroupingChoice,
+                    "legend_title": "Standard acquisition name",
+                    "colourmap": user_profile.plotColourMapChoice,
+                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                    "x_axis_title": "Compressed breast thickness (mm)",
+                    "y_axis_title": "mAs",
+                    "filename": "OpenREM CT standard acquisition name mAs vs thickness",
+                    "return_as_dict": return_as_dict,
+                }
+                return_structure["standardmAsvsThickness"] = (
+                    plotly_scatter(  # pylint: disable=line-too-long
+                        df_without_blanks,
+                        parameter_dict,
+                    )
+                )
+
+            if user_profile.plotMGStandardAcquisitionFreq:
+                parameter_dict = {
+                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
+                    "sorting_choice": [
+                        user_profile.plotInitialSortingDirection,
+                        user_profile.plotMGInitialSortingChoice,
+                    ],
+                    "legend_title": "Acquisition",
+                    "df_x_axis_col": "x_ray_system_name",
+                    "x_axis_title": "System",
+                    "grouping_choice": user_profile.plotGroupingChoice,
+                    "colourmap": user_profile.plotColourMapChoice,
+                    "filename": "OpenREM MG standard acquisition name frequency",
+                    "groupby_cols": None,
+                    "facet_col": None,
+                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                    "return_as_dict": return_as_dict,
+                }
+                (
+                    return_structure["standardAcquisitionFrequencyData"],
+                    return_structure["standardAcquisitionFrequencyDataCSV"],
+                ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
+                    df_without_blanks,
+                    parameter_dict,
+                    csv_name="standardAcquisitionFrequencyData.csv",
+                )
+
+            if user_profile.plotMGStandardAcquisitionAGDOverTime:
+                facet_title = "System"
+
+                if user_profile.plotGroupingChoice == "series":
+                    facet_title = "Standard acquisition name"
+
+                parameter_dict = {  # pylint: disable=line-too-long
+                    "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",
+                    "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__average_glandular_dose",  # pylint: disable=line-too-long
+                    "df_date_col": "study_date",
+                    "name_title": "Standard acquisition name",
+                    "value_title": "AGD (mGy)",
+                    "date_title": "Study date",
+                    "facet_title": facet_title,
+                    "sorting_choice": [
+                        user_profile.plotInitialSortingDirection,
+                        user_profile.plotMGInitialSortingChoice,
+                    ],
+                    "time_period": plot_timeunit_period,
+                    "average_choices": average_choices + ["count"],
+                    "grouping_choice": user_profile.plotGroupingChoice,
+                    "colourmap": user_profile.plotColourMapChoice,
+                    "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                    "filename": "OpenREM MG standard acquisition name AGD over time",
+                    "return_as_dict": return_as_dict,
+                }
+                result = construct_over_time_charts(
+                    df_without_blanks,
+                    parameter_dict,
+                )
+
+                if user_profile.plotMean:
+                    return_structure["standardAcquisitionMeanAGDOverTime"] = result[
+                        "mean"
+                    ]
+                if user_profile.plotMedian:
+                    return_structure["standardAcquisitionMedianAGDOverTime"] = result[
+                        "median"
+                    ]
+
+    #######################################################################
+    # Prepare study- and request-level Pandas DataFrame to use for charts
+    charts_of_interest = [user_profile.plotMGStudyPerDayAndHour]
+    if enable_standard_names:
+        charts_of_interest.append(user_profile.plotMGStandardStudyPerDayAndHour)
+
+    if any(charts_of_interest):
+
+        name_fields = []
+        charts_of_interest = [user_profile.plotMGStudyPerDayAndHour]
+        if any(charts_of_interest):
+            name_fields.append("study_description")
+
+        if enable_standard_names:
+            charts_of_interest = [user_profile.plotMGStandardStudyPerDayAndHour]
+            if any(charts_of_interest):
+                name_fields.append("standard_names__standard_name")
+
+        value_fields = []
+
+        date_fields = []
+        time_fields = []
+        charts_of_interest = [user_profile.plotMGStudyPerDayAndHour]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotMGStandardStudyPerDayAndHour)
+        if any(charts_of_interest):
+            date_fields.append("study_date")
+            time_fields.append("study_time")
+
+        system_field = []
+        if user_profile.plotSeriesPerSystem:
+            system_field.append(
+                "generalequipmentmoduleattr__unique_equipment_name_id__display_name"
+            )
+
+        fields = {
+            "names": name_fields,
+            "values": value_fields,
+            "dates": date_fields,
+            "times": time_fields,
+            "system": system_field,
+        }
+
+        # If only standard_names__standard_name is required then exclude all entries where these are None as these are
+        # not required for standard name charts.
+        queryset = f.qs
+        if name_fields == ["standard_names__standard_name"]:
+            queryset = queryset.exclude(standard_names__standard_name__isnull=True)
+
+        df = create_dataframe(
+            queryset,
+            fields,
+            data_point_name_lowercase=user_profile.plotCaseInsensitiveCategories,
+            data_point_name_remove_whitespace_padding=user_profile.plotRemoveCategoryWhitespacePadding,
+            uid="pk",
+        )
+
+        if user_profile.plotMGStudyPerDayAndHour:
+            df_time_series_per_weekday = create_dataframe_weekdays(
+                df, "study_description", df_date_col="study_date"
+            )
+
+            return_structure["studyWorkloadData"] = plotly_barchart_weekdays(
+                df_time_series_per_weekday,
+                "weekday",
+                "study_description",
+                name_axis_title="Weekday",
+                value_axis_title="Frequency",
+                colourmap=user_profile.plotColourMapChoice,
+                filename="OpenREM CT study description workload",
+                facet_col_wrap=user_profile.plotFacetColWrapVal,
+                sorting_choice=[
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotMGInitialSortingChoice,
+                ],
+                return_as_dict=return_as_dict,
+            )
+
+        if enable_standard_names:
+            charts_of_interest = [user_profile.plotMGStandardStudyPerDayAndHour]
+
+            if any(charts_of_interest):
+
+                # Create a standard name data frame - remove any blank standard names
+                standard_name_df = df[
+                    (df["standard_names__standard_name"] != "blank")
+                    & (df["standard_names__standard_name"] != "Blank")
+                ].copy()
+                # Remove any unused categories (this will include "Blank" or "blank")
+                standard_name_df["standard_names__standard_name"] = standard_name_df[
+                    "standard_names__standard_name"
+                ].cat.remove_unused_categories()
+
+                if user_profile.plotMGStandardStudyPerDayAndHour:
+                    df_time_series_per_weekday = create_dataframe_weekdays(
+                        standard_name_df,
+                        "standard_names__standard_name",
+                        df_date_col="study_date",
+                    )
+
+                    return_structure["standardStudyWorkloadData"] = (
+                        plotly_barchart_weekdays(
+                            df_time_series_per_weekday,
+                            "weekday",
+                            "standard_names__standard_name",
+                            name_axis_title="Weekday",
+                            value_axis_title="Frequency",
+                            colourmap=user_profile.plotColourMapChoice,
+                            filename="OpenREM CT standard study name workload",
+                            facet_col_wrap=user_profile.plotFacetColWrapVal,
+                            sorting_choice=[
+                                user_profile.plotInitialSortingDirection,
+                                user_profile.plotMGInitialSortingChoice,
+                            ],
+                            return_as_dict=return_as_dict,
+                        )
+                    )
+
+    return return_structure
+
+
+def mg_chart_form_processing(request, user_profile):
+    # pylint: disable=too-many-statements
+
+    # Obtain the system-level enable_standard_names setting
+    enable_standard_names = are_standard_names_enabled()
+
+    # Obtain the chart options from the request
+    chart_options_form = None
+    if enable_standard_names:
+        chart_options_form = MGChartOptionsFormIncStandard(request.GET)
+    else:
+        chart_options_form = MGChartOptionsForm(request.GET)
+
+    # Check whether the form data is valid
+    if chart_options_form.is_valid():
+        # Use the form data if the user clicked on the submit button
+        if "submit" in request.GET:
+            # process the data in form.cleaned_data as required
+
+            set_common_chart_options(chart_options_form, user_profile)
+
+            set_average_chart_options(chart_options_form, user_profile)
+
+            set_mg_chart_options(chart_options_form, user_profile)
+
+            user_profile.save()
+
+        else:
+            average_choices = required_average_choices(user_profile)
+
+            mg_form_data = initialise_mg_form_data(user_profile)
+
+            form_data = {
+                "plotCharts": user_profile.plotCharts,
+                "plotGrouping": user_profile.plotGroupingChoice,
+                "plotSeriesPerSystem": user_profile.plotSeriesPerSystem,
+                "plotHistograms": user_profile.plotHistograms,
+                "plotInitialSortingDirection": user_profile.plotInitialSortingDirection,
+                "plotAverageChoice": average_choices,
+            }
+
+            form_data = {**form_data, **mg_form_data}
+
+            chart_options_form = None
+            if enable_standard_names:
+                chart_options_form = MGChartOptionsFormIncStandard(form_data)
+            else:
+                chart_options_form = MGChartOptionsForm(form_data)
+
+    return chart_options_form