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

openrem/remapp/views_charts_dx.py

@@ -1,1836 +1,1906 @@
-# 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 (
-    DXChartOptionsForm,
-    DXChartOptionsFormIncStandard,
-)
-from remapp.interface.mod_filters import dx_acq_filter
-from remapp.models import (
-    create_user_profile,
-    StandardNameSettings,
-)
-from remapp.views_admin import (
-    set_average_chart_options,
-    required_average_choices,
-    initialise_dx_form_data,
-    set_dx_chart_options,
-    set_common_chart_options,
-)
-from .interface.chart_functions import (
-    create_dataframe,
-    create_dataframe_weekdays,
-    create_dataframe_aggregates,
-    plotly_boxplot,
-    plotly_barchart,
-    plotly_histogram_barchart,
-    plotly_barchart_weekdays,
-    plotly_set_default_theme,
-    plotly_frequency_barchart,
-    plotly_scatter,
-    construct_over_time_charts,
-    generate_average_chart_group,
-)
-
-logger = logging.getLogger(__name__)
-
-
-def generate_required_dx_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.plotDXAcquisitionMeanDAPOverTime,
-        profile.plotDXAcquisitionMeanmAsOverTime,
-        profile.plotDXAcquisitionMeankVpOverTime,
-    ]
-    if enable_standard_names:
-        charts_of_interest.append(profile.plotDXStandardAcquisitionMeanDAPOverTime)
-        charts_of_interest.append(profile.plotDXStandardAcquisitionMeanmAsOverTime)
-        charts_of_interest.append(profile.plotDXStandardAcquisitionMeankVpOverTime)
-
-    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.plotDXAcquisitionMeanDAPOverTimePeriod)]
-        ).lower()
-
-    if profile.plotDXAcquisitionMeanDAP:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of mean DAP for each acquisition protocol",
-                    "var_name": "acquisitionMeanDAP",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of median DAP for each acquisition protocol",
-                    "var_name": "acquisitionMedianDAP",
-                }
-            )
-        if profile.plotBoxplots:
-            required_charts.append(
-                {
-                    "title": "Boxplot of DAP for each acquisition protocol",
-                    "var_name": "acquisitionBoxplotDAP",
-                }
-            )
-        if profile.plotHistograms:
-            required_charts.append(
-                {
-                    "title": "Histogram of DAP for each acquisition protocol",
-                    "var_name": "acquisitionHistogramDAP",
-                }
-            )
-
-    if profile.plotDXAcquisitionFreq:
-        required_charts.append(
-            {
-                "title": "Chart of acquisition protocol frequency",
-                "var_name": "acquisitionFrequency",
-            }
-        )
-
-    if profile.plotDXAcquisitionMeankVp:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of mean kVp for each acquisition protocol",
-                    "var_name": "acquisitionMeankVp",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of median kVp for each acquisition protocol",
-                    "var_name": "acquisitionMediankVp",
-                }
-            )
-        if profile.plotBoxplots:
-            required_charts.append(
-                {
-                    "title": "Boxplot of kVp for each acquisition protocol",
-                    "var_name": "acquisitionBoxplotkVp",
-                }
-            )
-        if profile.plotHistograms:
-            required_charts.append(
-                {
-                    "title": "Histogram of kVp for each acquisition protocol",
-                    "var_name": "acquisitionHistogramkVp",
-                }
-            )
-
-    if profile.plotDXAcquisitionMeanmAs:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of mean mAs for each acquisition protocol",
-                    "var_name": "acquisitionMeanmAs",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of median mAs for each acquisition protocol",
-                    "var_name": "acquisitionMedianmAs",
-                }
-            )
-        if profile.plotBoxplots:
-            required_charts.append(
-                {
-                    "title": "Boxplot of mAs for each acquisition protocol",
-                    "var_name": "acquisitionBoxplotmAs",
-                }
-            )
-        if profile.plotHistograms:
-            required_charts.append(
-                {
-                    "title": "Histogram of mAs for each acquisition protocol",
-                    "var_name": "acquisitionHistogrammAs",
-                }
-            )
-
-    if profile.plotDXAcquisitionMeanDAPOverTime:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of mean DAP per acquisition protocol over time ("
-                    + time_period
-                    + ")",
-                    "var_name": "acquisitionMeanDAPOverTime",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of median DAP per acquisition protocol over time ("
-                    + time_period
-                    + ")",
-                    "var_name": "acquisitionMedianDAPOverTime",
-                }
-            )
-
-    if profile.plotDXAcquisitionMeankVpOverTime:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of mean kVp per acquisition protocol over time ("
-                    + time_period
-                    + ")",
-                    "var_name": "acquisitionMeankVpOverTime",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of median kVp per acquisition protocol over time ("
-                    + time_period
-                    + ")",
-                    "var_name": "acquisitionMediankVpOverTime",
-                }
-            )
-
-    if profile.plotDXAcquisitionMeanmAsOverTime:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of mean mAs per acquisition protocol over time ("
-                    + time_period
-                    + ")",
-                    "var_name": "acquisitionMeanmAsOverTime",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of median mAs per acquisition protocol over time ("
-                    + time_period
-                    + ")",
-                    "var_name": "acquisitionMedianmAsOverTime",
-                }
-            )
-
-    if profile.plotDXAcquisitionDAPvsMass:
-        required_charts.append(
-            {
-                "title": "Chart of acquisition protocol DAP vs patient mass",
-                "var_name": "acquisitionDAPvsMass",
-            }
-        )
-
-    if enable_standard_names:
-        if profile.plotDXStandardAcquisitionMeanDAP:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of mean DAP for each standard acquisition name",
-                        "var_name": "standardAcquisitionMeanDAP",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of median DAP for each standard acquisition name",
-                        "var_name": "standardAcquisitionMedianDAP",
-                    }
-                )
-            if profile.plotBoxplots:
-                required_charts.append(
-                    {
-                        "title": "Boxplot of DAP for each standard acquisition name",
-                        "var_name": "standardAcquisitionBoxplotDAP",
-                    }
-                )
-            if profile.plotHistograms:
-                required_charts.append(
-                    {
-                        "title": "Histogram of DAP for each standard acquisition name",
-                        "var_name": "standardAcquisitionHistogramDAP",
-                    }
-                )
-
-        if profile.plotDXStandardAcquisitionFreq:
-            required_charts.append(
-                {
-                    "title": "Chart of standard acquisition name frequency",
-                    "var_name": "standardAcquisitionFrequency",
-                }
-            )
-
-        if profile.plotDXStandardAcquisitionMeankVp:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of mean kVp for each standard acquisition name",
-                        "var_name": "standardAcquisitionMeankVp",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of median kVp for each standard acquisition name",
-                        "var_name": "standardAcquisitionMediankVp",
-                    }
-                )
-            if profile.plotBoxplots:
-                required_charts.append(
-                    {
-                        "title": "Boxplot of kVp for each standard acquisition name",
-                        "var_name": "standardAcquisitionBoxplotkVp",
-                    }
-                )
-            if profile.plotHistograms:
-                required_charts.append(
-                    {
-                        "title": "Histogram of kVp for each standard acquisition name",
-                        "var_name": "standardAcquisitionHistogramkVp",
-                    }
-                )
-
-        if profile.plotDXStandardAcquisitionMeanmAs:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of mean mAs for each standard acquisition name",
-                        "var_name": "standardAcquisitionMeanmAs",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of median mAs for each standard acquisition name",
-                        "var_name": "standardAcquisitionMedianmAs",
-                    }
-                )
-            if profile.plotBoxplots:
-                required_charts.append(
-                    {
-                        "title": "Boxplot of mAs for each standard acquisition name",
-                        "var_name": "standardAcquisitionBoxplotmAs",
-                    }
-                )
-            if profile.plotHistograms:
-                required_charts.append(
-                    {
-                        "title": "Histogram of mAs for each standard acquisition name",
-                        "var_name": "standardAcquisitionHistogrammAs",
-                    }
-                )
-
-        if profile.plotDXStandardAcquisitionMeanDAPOverTime:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of mean DAP per standard acquisition name over time ("
-                        + time_period
-                        + ")",
-                        "var_name": "standardAcquisitionMeanDAPOverTime",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of median DAP per standard acquisition name over time ("
-                        + time_period
-                        + ")",
-                        "var_name": "standardAcquisitionMedianDAPOverTime",
-                    }
-                )
-
-        if profile.plotDXStandardAcquisitionMeankVpOverTime:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of mean kVp per standard acquisition name over time ("
-                        + time_period
-                        + ")",
-                        "var_name": "standardAcquisitionMeankVpOverTime",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of median kVp per standard acquisition name over time ("
-                        + time_period
-                        + ")",
-                        "var_name": "standardAcquisitionMediankVpOverTime",
-                    }
-                )
-
-        if profile.plotDXStandardAcquisitionMeanmAsOverTime:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of mean mAs per standard acquisition name over time ("
-                        + time_period
-                        + ")",
-                        "var_name": "standardAcquisitionMeanmAsOverTime",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of median mAs per standard acquisition name over time ("
-                        + time_period
-                        + ")",
-                        "var_name": "standardAcquisitionMedianmAsOverTime",
-                    }
-                )
-
-        if profile.plotDXStandardAcquisitionDAPvsMass:
-            required_charts.append(
-                {
-                    "title": "Chart of standard acquisition name DAP vs patient mass",
-                    "var_name": "standardAcquisitionDAPvsMass",
-                }
-            )
-
-        if profile.plotDXStandardStudyMeanDAP:
-            if profile.plotMean:
-                required_charts.append(
-                    {
-                        "title": "Chart of mean DAP for each standard study name",
-                        "var_name": "standardStudyMeanDAP",
-                    }
-                )
-            if profile.plotMedian:
-                required_charts.append(
-                    {
-                        "title": "Chart of median DAP for each standard study name",
-                        "var_name": "standardStudyMedianDAP",
-                    }
-                )
-            if profile.plotBoxplots:
-                required_charts.append(
-                    {
-                        "title": "Boxplot of DAP for each standard study name",
-                        "var_name": "standardStudyBoxplotDAP",
-                    }
-                )
-            if profile.plotHistograms:
-                required_charts.append(
-                    {
-                        "title": "Histogram of DAP for each standard study name",
-                        "var_name": "standardStudyHistogramDAP",
-                    }
-                )
-
-        if profile.plotDXStandardStudyFreq:
-            required_charts.append(
-                {
-                    "title": "Chart of standard study name frequency",
-                    "var_name": "standardStudyFrequency",
-                }
-            )
-
-        if profile.plotDXStandardStudyPerDayAndHour:
-            required_charts.append(
-                {
-                    "title": "Chart of standard study name workload",
-                    "var_name": "standardStudyWorkload",
-                }
-            )
-
-        if profile.plotDXStandardStudyDAPvsMass:
-            required_charts.append(
-                {
-                    "title": "Chart of standard study name DAP vs patient mass",
-                    "var_name": "standardStudyDAPvsMass",
-                }
-            )
-
-    if profile.plotDXStudyMeanDAP:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of mean DAP for each study description",
-                    "var_name": "studyMeanDAP",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of median DAP for each study description",
-                    "var_name": "studyMedianDAP",
-                }
-            )
-        if profile.plotBoxplots:
-            required_charts.append(
-                {
-                    "title": "Boxplot of DAP for each study description",
-                    "var_name": "studyBoxplotDAP",
-                }
-            )
-        if profile.plotHistograms:
-            required_charts.append(
-                {
-                    "title": "Histogram of DAP for each study description",
-                    "var_name": "studyHistogramDAP",
-                }
-            )
-
-    if profile.plotDXStudyFreq:
-        required_charts.append(
-            {
-                "title": "Chart of study description frequency",
-                "var_name": "studyFrequency",
-            }
-        )
-
-    if profile.plotDXStudyPerDayAndHour:
-        required_charts.append(
-            {
-                "title": "Chart of study description workload",
-                "var_name": "studyWorkload",
-            }
-        )
-
-    if profile.plotDXStudyDAPvsMass:
-        required_charts.append(
-            {
-                "title": "Chart of study description DAP vs patient mass",
-                "var_name": "studyDAPvsMass",
-            }
-        )
-
-    if profile.plotDXRequestMeanDAP:
-        if profile.plotMean:
-            required_charts.append(
-                {
-                    "title": "Chart of mean DAP for each requested procedure",
-                    "var_name": "requestMeanDAP",
-                }
-            )
-        if profile.plotMedian:
-            required_charts.append(
-                {
-                    "title": "Chart of median DAP for each requested procedure",
-                    "var_name": "requestMedianDAP",
-                }
-            )
-        if profile.plotBoxplots:
-            required_charts.append(
-                {
-                    "title": "Boxplot of DAP for each requested procedure",
-                    "var_name": "requestBoxplotDAP",
-                }
-            )
-        if profile.plotHistograms:
-            required_charts.append(
-                {
-                    "title": "Histogram of DAP for each requested procedure",
-                    "var_name": "requestHistogramDAP",
-                }
-            )
-
-    if profile.plotDXRequestFreq:
-        required_charts.append(
-            {
-                "title": "Chart of requested procedure frequency",
-                "var_name": "requestFrequency",
-            }
-        )
-
-    if profile.plotDXRequestDAPvsMass:
-        required_charts.append(
-            {
-                "title": "Chart of requested procedure DAP vs patient mass",
-                "var_name": "requestDAPvsMass",
-            }
-        )
-
-    return required_charts
-
-
-@login_required
-def dx_summary_chart_data(request):
-    """Obtain data for Ajax chart call."""
-    pid = bool(request.user.groups.filter(name="pidgroup"))
-    f = dx_acq_filter(request.GET, pid=pid)
-
-    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 = dx_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 dx_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 radiographic 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.plotDXAcquisitionMeanDAPOverTime,
-        user_profile.plotDXAcquisitionMeankVpOverTime,
-        user_profile.plotDXAcquisitionMeanmAsOverTime,
-    ]
-    if enable_standard_names:
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanDAPOverTime)
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeankVpOverTime)
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanmAsOverTime)
-    if any(charts_of_interest):
-        plot_timeunit_period = user_profile.plotDXAcquisitionMeanDAPOverTimePeriod
-
-    #######################################################################
-    # Prepare acquisition-level Pandas DataFrame to use for charts
-    charts_of_interest = [
-        user_profile.plotDXAcquisitionMeanDAP,
-        user_profile.plotDXAcquisitionFreq,
-        user_profile.plotDXAcquisitionMeankVp,
-        user_profile.plotDXAcquisitionMeanmAs,
-        user_profile.plotDXAcquisitionMeankVpOverTime,
-        user_profile.plotDXAcquisitionMeanmAsOverTime,
-        user_profile.plotDXAcquisitionMeanDAPOverTime,
-        user_profile.plotDXAcquisitionDAPvsMass,
-    ]
-    if enable_standard_names:
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanDAP)
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionFreq)
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeankVp)
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanmAs)
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeankVpOverTime)
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanmAsOverTime)
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanDAPOverTime)
-        charts_of_interest.append(user_profile.plotDXStandardAcquisitionDAPvsMass)
-
-    if any(charts_of_interest):
-
-        name_fields = []
-        charts_of_interest = [
-            user_profile.plotDXAcquisitionMeanDAP,
-            user_profile.plotDXAcquisitionFreq,
-            user_profile.plotDXAcquisitionMeankVp,
-            user_profile.plotDXAcquisitionMeanmAs,
-            user_profile.plotDXAcquisitionMeankVpOverTime,
-            user_profile.plotDXAcquisitionMeanmAsOverTime,
-            user_profile.plotDXAcquisitionMeanDAPOverTime,
-            user_profile.plotDXAcquisitionDAPvsMass,
-        ]
-        if any(charts_of_interest):
-            name_fields.append(
-                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
-            )
-
-        if enable_standard_names:
-            charts_of_interest = [
-                user_profile.plotDXStandardAcquisitionMeanDAP,
-                user_profile.plotDXStandardAcquisitionFreq,
-                user_profile.plotDXStandardAcquisitionMeankVp,
-                user_profile.plotDXStandardAcquisitionMeanmAs,
-                user_profile.plotDXStandardAcquisitionMeankVpOverTime,
-                user_profile.plotDXStandardAcquisitionMeanmAsOverTime,
-                user_profile.plotDXStandardAcquisitionMeanDAPOverTime,
-                user_profile.plotDXStandardAcquisitionDAPvsMass,
-            ]
-            if any(charts_of_interest):
-                name_fields.append(
-                    "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-                )
-
-        value_fields = []
-        value_multipliers = []
-        charts_of_interest = [
-            user_profile.plotDXAcquisitionMeanDAP,
-            user_profile.plotDXAcquisitionMeanDAPOverTime,
-            user_profile.plotDXAcquisitionDAPvsMass,
-        ]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanDAP)
-            charts_of_interest.append(
-                user_profile.plotDXStandardAcquisitionMeanDAPOverTime
-            )
-            charts_of_interest.append(user_profile.plotDXStandardAcquisitionDAPvsMass)
-        if any(charts_of_interest):
-            value_fields.append(
-                "projectionxrayradiationdose__irradeventxraydata__dose_area_product"
-            )
-            value_multipliers.append(1000000)
-
-        charts_of_interest = [
-            user_profile.plotDXAcquisitionMeankVp,
-            user_profile.plotDXAcquisitionMeankVpOverTime,
-        ]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeankVp)
-            charts_of_interest.append(
-                user_profile.plotDXStandardAcquisitionMeankVpOverTime
-            )
-        if any(charts_of_interest):
-            value_fields.append(
-                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp"
-            )
-            value_multipliers.append(1)
-
-        charts_of_interest = [
-            user_profile.plotDXAcquisitionMeanmAs,
-            user_profile.plotDXAcquisitionMeanmAsOverTime,
-        ]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanmAs)
-            charts_of_interest.append(
-                user_profile.plotDXStandardAcquisitionMeanmAsOverTime
-            )
-        if any(charts_of_interest):
-            value_fields.append(
-                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure"
-            )
-            value_multipliers.append(0.001)
-
-        charts_of_interest = [user_profile.plotDXAcquisitionDAPvsMass]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotDXStandardAcquisitionDAPvsMass)
-        if any(charts_of_interest):
-            value_fields.append("patientstudymoduleattr__patient_weight")
-            value_multipliers.append(1)
-
-        time_fields = []
-        date_fields = []
-        charts_of_interest = [
-            user_profile.plotDXAcquisitionMeanDAPOverTime,
-            user_profile.plotDXAcquisitionMeankVpOverTime,
-            user_profile.plotDXAcquisitionMeanmAsOverTime,
-        ]
-        if enable_standard_names:
-            charts_of_interest.append(
-                user_profile.plotDXStandardAcquisitionMeanDAPOverTime
-            )
-            charts_of_interest.append(
-                user_profile.plotDXStandardAcquisitionMeankVpOverTime
-            )
-            charts_of_interest.append(
-                user_profile.plotDXStandardAcquisitionMeanmAsOverTime
-            )
-        if any(charts_of_interest):
-            date_fields.append("study_date")
-
-        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.plotDXAcquisitionMeanDAP:
-
-            name_field = (
-                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
-            )
-            value_field = (
-                "projectionxrayradiationdose__irradeventxraydata__dose_area_product"
-            )
-            value_text = "DAP"
-            units_text = "(cGy.cm<sup>2</sup>)"
-            name_text = "Acquisition protocol"
-            variable_name_start = "acquisition"
-            variable_value_name = "DAP"
-            modality_text = "DX"
-            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.plotDXInitialSortingChoice,
-            )
-
-            return_structure = {**return_structure, **new_charts}
-
-        if user_profile.plotDXAcquisitionMeankVp:
-
-            name_field = (
-                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
-            )
-            value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp"
-            value_text = "kVp"
-            units_text = ""
-            name_text = "Acquisition protocol"
-            variable_name_start = "acquisition"
-            variable_value_name = "kVp"
-            modality_text = "DX"
-            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.plotDXInitialSortingChoice,
-            )
-
-            return_structure = {**return_structure, **new_charts}
-
-        if user_profile.plotDXAcquisitionMeanmAs:
-
-            name_field = (
-                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
-            )
-            value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure"  # pylint: disable=line-too-long
-            value_text = "mAs"
-            units_text = ""
-            name_text = "Acquisition protocol"
-            variable_name_start = "acquisition"
-            variable_value_name = "mAs"
-            modality_text = "DX"
-            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.plotDXInitialSortingChoice,
-            )
-
-            return_structure = {**return_structure, **new_charts}
-
-        if user_profile.plotDXAcquisitionFreq:
-            parameter_dict = {
-                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                "legend_title": "Acquisition protocol",
-                "df_x_axis_col": "x_ray_system_name",
-                "x_axis_title": "System",
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "colourmap": user_profile.plotColourMapChoice,
-                "filename": "OpenREM DX 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.plotDXAcquisitionMeanDAPOverTime:
-            facet_title = "System"
-
-            if user_profile.plotGroupingChoice == "series":
-                facet_title = "Acquisition protocol"
-
-            parameter_dict = {
-                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
-                "df_value_col": "projectionxrayradiationdose__irradeventxraydata__dose_area_product",
-                "df_date_col": "study_date",
-                "name_title": "Acquisition protocol",
-                "value_title": "DAP (cGy.cm<sup>2</sup>)",
-                "date_title": "Study date",
-                "facet_title": facet_title,
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                "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 DX acquisition protocol DAP over time",
-                "return_as_dict": return_as_dict,
-            }
-            result = construct_over_time_charts(
-                df,
-                parameter_dict,
-            )
-
-            if user_profile.plotMean:
-                return_structure["acquisitionMeanDAPOverTime"] = result["mean"]
-            if user_profile.plotMedian:
-                return_structure["acquisitionMedianDAPOverTime"] = result["median"]
-
-        if user_profile.plotDXAcquisitionMeankVpOverTime:
-            facet_title = "System"
-
-            if user_profile.plotGroupingChoice == "series":
-                facet_title = "Acquisition protocol"
-
-            parameter_dict = {
-                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
-                "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp",
-                "df_date_col": "study_date",
-                "name_title": "Acquisition protocol",
-                "value_title": "kVp",
-                "date_title": "Study date",
-                "facet_title": facet_title,
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                "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 DX acquisition protocol kVp over time",
-                "return_as_dict": return_as_dict,
-            }
-            result = construct_over_time_charts(
-                df,
-                parameter_dict,
-            )
-
-            if user_profile.plotMean:
-                return_structure["acquisitionMeankVpOverTime"] = result["mean"]
-            if user_profile.plotMedian:
-                return_structure["acquisitionMediankVpOverTime"] = result["median"]
-
-        if user_profile.plotDXAcquisitionMeanmAsOverTime:
-            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__exposure__exposure",  # pylint: disable=line-too-long
-                "df_date_col": "study_date",
-                "name_title": "Acquisition protocol",
-                "value_title": "mAs",
-                "date_title": "Study date",
-                "facet_title": facet_title,
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                "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 DX acquisition protocol mAs over time",
-                "return_as_dict": return_as_dict,
-            }
-            result = construct_over_time_charts(
-                df,
-                parameter_dict,
-            )
-
-            if user_profile.plotMean:
-                return_structure["acquisitionMeanmAsOverTime"] = result["mean"]
-            if user_profile.plotMedian:
-                return_structure["acquisitionMedianmAsOverTime"] = result["median"]
-
-        if user_profile.plotDXAcquisitionDAPvsMass:
-            parameter_dict = {
-                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
-                "df_x_col": "patientstudymoduleattr__patient_weight",
-                "df_y_col": "projectionxrayradiationdose__irradeventxraydata__dose_area_product",
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "legend_title": "Acquisition protocol",
-                "colourmap": user_profile.plotColourMapChoice,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "x_axis_title": "Patient mass (kg)",
-                "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
-                "filename": "OpenREM DX acquisition protocol DAP vs patient mass",
-                "return_as_dict": return_as_dict,
-            }
-            return_structure["acquisitionDAPvsMass"] = plotly_scatter(
-                df,
-                parameter_dict,
-            )
-
-        if enable_standard_names:
-            charts_of_interest = [
-                user_profile.plotDXStandardAcquisitionMeanDAP,
-                user_profile.plotDXStandardAcquisitionFreq,
-                user_profile.plotDXStandardAcquisitionMeankVp,
-                user_profile.plotDXStandardAcquisitionMeanmAs,
-                user_profile.plotDXStandardAcquisitionMeankVpOverTime,
-                user_profile.plotDXStandardAcquisitionMeanmAsOverTime,
-                user_profile.plotDXStandardAcquisitionMeanDAPOverTime,
-                user_profile.plotDXStandardAcquisitionDAPvsMass,
-            ]
-
-            if any(charts_of_interest):
-
-                # Exclude "Blank" and "blank" standard_acqusition_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.plotDXStandardAcquisitionMeanDAP:
-                    name_field = "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-                    value_field = "projectionxrayradiationdose__irradeventxraydata__dose_area_product"
-                    value_text = "DAP"
-                    units_text = "(cGy.cm<sup>2</sup>)"
-                    name_text = "Standard acquisition name"
-                    variable_name_start = "standardAcquisition"
-                    variable_value_name = "DAP"
-                    modality_text = "DX"
-                    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.plotDXInitialSortingChoice,
-                    )
-
-                    return_structure = {**return_structure, **new_charts}
-
-                if user_profile.plotDXStandardAcquisitionMeankVp:
-                    name_field = "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-                    value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp"
-                    value_text = "kVp"
-                    units_text = ""
-                    name_text = "Standard acquisition name"
-                    variable_name_start = "standardAcquisition"
-                    variable_value_name = "kVp"
-                    modality_text = "DX"
-                    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.plotDXInitialSortingChoice,
-                    )
-
-                    return_structure = {**return_structure, **new_charts}
-
-                if user_profile.plotDXStandardAcquisitionMeanmAs:
-                    name_field = "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
-                    value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure"  # pylint: disable=line-too-long
-                    value_text = "mAs"
-                    units_text = ""
-                    name_text = "Standard acquisition name"
-                    variable_name_start = "standardAcquisition"
-                    variable_value_name = "mAs"
-                    modality_text = "DX"
-                    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.plotDXInitialSortingChoice,
-                    )
-
-                    return_structure = {**return_structure, **new_charts}
-
-                if user_profile.plotDXStandardAcquisitionFreq:
-                    parameter_dict = {
-                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
-                        "sorting_choice": [
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotDXInitialSortingChoice,
-                        ],
-                        "legend_title": "Standard acquisition name",
-                        "df_x_axis_col": "x_ray_system_name",
-                        "x_axis_title": "System",
-                        "grouping_choice": user_profile.plotGroupingChoice,
-                        "colourmap": user_profile.plotColourMapChoice,
-                        "filename": "OpenREM DX 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(
-                        df_without_blanks,
-                        parameter_dict,
-                        csv_name="standardAcquisitionFrequencyData.csv",
-                    )
-
-                if user_profile.plotDXStandardAcquisitionMeanDAPOverTime:
-                    facet_title = "System"
-
-                    if user_profile.plotGroupingChoice == "series":
-                        facet_title = "Standard acquisition name"
-
-                    parameter_dict = {
-                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
-                        "df_value_col": "projectionxrayradiationdose__irradeventxraydata__dose_area_product",
-                        "df_date_col": "study_date",
-                        "name_title": "Standard acquisition name",
-                        "value_title": "DAP (cGy.cm<sup>2</sup>)",
-                        "date_title": "Study date",
-                        "facet_title": facet_title,
-                        "sorting_choice": [
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotDXInitialSortingChoice,
-                        ],
-                        "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 DX standard acquisition name DAP over time",
-                        "return_as_dict": return_as_dict,
-                    }
-                    result = construct_over_time_charts(
-                        df_without_blanks,
-                        parameter_dict,
-                    )
-
-                    if user_profile.plotMean:
-                        return_structure["standardAcquisitionMeanDAPOverTime"] = result[
-                            "mean"
-                        ]
-                    if user_profile.plotMedian:
-                        return_structure[
-                            "standardAcquisitionMedianDAPOverTime"
-                        ] = result["median"]
-
-                if user_profile.plotDXStandardAcquisitionMeankVpOverTime:
-                    facet_title = "System"
-
-                    if user_profile.plotGroupingChoice == "series":
-                        facet_title = "Standard acquisition name"
-
-                    parameter_dict = {
-                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
-                        "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp",  # pylint: disable=line-too-long
-                        "df_date_col": "study_date",
-                        "name_title": "Standard acquisition name",
-                        "value_title": "kVp",
-                        "date_title": "Study date",
-                        "facet_title": facet_title,
-                        "sorting_choice": [
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotDXInitialSortingChoice,
-                        ],
-                        "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 DX standard acquisition name kVp over time",
-                        "return_as_dict": return_as_dict,
-                    }
-                    result = construct_over_time_charts(
-                        df_without_blanks,
-                        parameter_dict,
-                    )
-
-                    if user_profile.plotMean:
-                        return_structure["standardAcquisitionMeankVpOverTime"] = result[
-                            "mean"
-                        ]
-                    if user_profile.plotMedian:
-                        return_structure[
-                            "standardAcquisitionMediankVpOverTime"
-                        ] = result["median"]
-
-                if user_profile.plotDXStandardAcquisitionMeanmAsOverTime:
-                    facet_title = "System"
-
-                    if user_profile.plotGroupingChoice == "series":
-                        facet_title = "Standard acquisition name"
-
-                    parameter_dict = {
-                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
-                        "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure",  # pylint: disable=line-too-long
-                        "df_date_col": "study_date",
-                        "name_title": "Standard acquisition name",
-                        "value_title": "mAs",
-                        "date_title": "Study date",
-                        "facet_title": facet_title,
-                        "sorting_choice": [
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotDXInitialSortingChoice,
-                        ],
-                        "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 DX standard acquisition name mAs over time",
-                        "return_as_dict": return_as_dict,
-                    }
-                    result = construct_over_time_charts(
-                        df_without_blanks,
-                        parameter_dict,
-                    )
-
-                    if user_profile.plotMean:
-                        return_structure["standardAcquisitionMeanmAsOverTime"] = result[
-                            "mean"
-                        ]
-                    if user_profile.plotMedian:
-                        return_structure[
-                            "standardAcquisitionMedianmAsOverTime"
-                        ] = result["median"]
-
-                if user_profile.plotDXStandardAcquisitionDAPvsMass:
-                    parameter_dict = {
-                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
-                        "df_x_col": "patientstudymoduleattr__patient_weight",
-                        "df_y_col": "projectionxrayradiationdose__irradeventxraydata__dose_area_product",
-                        "sorting_choice": [
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotDXInitialSortingChoice,
-                        ],
-                        "grouping_choice": user_profile.plotGroupingChoice,
-                        "legend_title": "Standard acquisition name",
-                        "colourmap": user_profile.plotColourMapChoice,
-                        "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                        "x_axis_title": "Patient mass (kg)",
-                        "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
-                        "filename": "OpenREM DX standard acquisition name DAP vs patient mass",
-                        "return_as_dict": return_as_dict,
-                    }
-                    return_structure["standardAcquisitionDAPvsMass"] = plotly_scatter(
-                        df_without_blanks,
-                        parameter_dict,
-                    )
-
-    #######################################################################
-    # Prepare study- and request-level Pandas DataFrame to use for charts
-    charts_of_interest = [
-        user_profile.plotDXStudyMeanDAP,
-        user_profile.plotDXStudyFreq,
-        user_profile.plotDXStudyPerDayAndHour,
-        user_profile.plotDXStudyDAPvsMass,
-        user_profile.plotDXRequestMeanDAP,
-        user_profile.plotDXRequestFreq,
-        user_profile.plotDXRequestDAPvsMass,
-    ]
-    if enable_standard_names:
-        charts_of_interest.append(user_profile.plotDXStandardStudyMeanDAP)
-        charts_of_interest.append(user_profile.plotDXStandardStudyFreq)
-        charts_of_interest.append(user_profile.plotDXStandardStudyPerDayAndHour)
-        charts_of_interest.append(user_profile.plotDXStandardStudyDAPvsMass)
-
-    if any(charts_of_interest):
-
-        name_fields = []
-        charts_of_interest = [
-            user_profile.plotDXStudyMeanDAP,
-            user_profile.plotDXStudyFreq,
-            user_profile.plotDXStudyPerDayAndHour,
-            user_profile.plotDXStudyDAPvsMass,
-        ]
-        if any(charts_of_interest):
-            name_fields.append("study_description")
-
-        charts_of_interest = [
-            user_profile.plotDXRequestMeanDAP,
-            user_profile.plotDXRequestFreq,
-            user_profile.plotDXRequestDAPvsMass,
-        ]
-        if any(charts_of_interest):
-            name_fields.append("requested_procedure_code_meaning")
-
-        if enable_standard_names:
-            charts_of_interest = [
-                user_profile.plotDXStandardStudyMeanDAP,
-                user_profile.plotDXStandardStudyFreq,
-                user_profile.plotDXStandardStudyPerDayAndHour,
-                user_profile.plotDXStandardStudyDAPvsMass,
-            ]
-            if any(charts_of_interest):
-                name_fields.append("standard_names__standard_name")
-
-        value_fields = []
-        value_multipliers = []
-        charts_of_interest = [
-            user_profile.plotDXStudyMeanDAP,
-            user_profile.plotDXRequestMeanDAP,
-            user_profile.plotDXStudyDAPvsMass,
-            user_profile.plotDXRequestDAPvsMass,
-        ]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotDXStandardStudyMeanDAP)
-            charts_of_interest.append(user_profile.plotDXStandardStudyDAPvsMass)
-
-        if any(charts_of_interest):
-            value_fields.append("total_dap")
-            value_multipliers.append(1000000)
-
-        charts_of_interest = [
-            user_profile.plotDXStudyDAPvsMass,
-            user_profile.plotDXRequestDAPvsMass,
-        ]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotDXStandardStudyDAPvsMass)
-
-        if any(charts_of_interest):
-            value_fields.append("patientstudymoduleattr__patient_weight")
-            value_multipliers.append(1)
-
-        date_fields = []
-        time_fields = []
-        charts_of_interest = [user_profile.plotDXStudyPerDayAndHour]
-        if enable_standard_names:
-            charts_of_interest.append(user_profile.plotDXStandardStudyPerDayAndHour)
-
-        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,
-            data_point_value_multipliers=value_multipliers,
-            char_wrap=user_profile.plotLabelCharWrap,
-            uid="pk",
-        )
-        #######################################################################
-
-        if user_profile.plotDXStudyMeanDAP:
-
-            name_field = "study_description"
-            value_field = "total_dap"
-            value_text = "DAP"
-            units_text = "(cGy.cm<sup>2</sup>)"
-            name_text = "Study description"
-            variable_name_start = "study"
-            variable_value_name = "DAP"
-            modality_text = "DX"
-            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.plotDXInitialSortingChoice,
-            )
-
-            return_structure = {**return_structure, **new_charts}
-
-        if user_profile.plotDXStudyFreq:
-            parameter_dict = {
-                "df_name_col": "study_description",
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                "legend_title": "Study description",
-                "df_x_axis_col": "x_ray_system_name",
-                "x_axis_title": "System",
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "colourmap": user_profile.plotColourMapChoice,
-                "filename": "OpenREM DX study description frequency",
-                "groupby_cols": None,
-                "facet_col": None,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "return_as_dict": return_as_dict,
-            }
-            (
-                return_structure["studyFrequencyData"],
-                return_structure["studyFrequencyDataCSV"],
-            ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
-                df,
-                parameter_dict,
-                csv_name="studyFrequencyData.csv",
-            )
-
-        if user_profile.plotDXRequestMeanDAP:
-
-            name_field = "requested_procedure_code_meaning"
-            value_field = "total_dap"
-            value_text = "DAP"
-            units_text = "(cGy.cm<sup>2</sup>)"
-            name_text = "Requested procedure"
-            variable_name_start = "request"
-            variable_value_name = "DAP"
-            modality_text = "DX"
-            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.plotDXInitialSortingChoice,
-            )
-
-            return_structure = {**return_structure, **new_charts}
-
-        if user_profile.plotDXRequestFreq:
-            parameter_dict = {
-                "df_name_col": "requested_procedure_code_meaning",
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                "legend_title": "Requested procedure",
-                "df_x_axis_col": "x_ray_system_name",
-                "x_axis_title": "System",
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "colourmap": user_profile.plotColourMapChoice,
-                "filename": "OpenREM DX requested procedure frequency",
-                "groupby_cols": None,
-                "facet_col": None,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "return_as_dict": return_as_dict,
-            }
-            (
-                return_structure["requestFrequencyData"],
-                return_structure["requestFrequencyDataCSV"],
-            ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
-                df,
-                parameter_dict,
-                csv_name="requestFrequencyData.csv",
-            )
-
-        if user_profile.plotDXStudyPerDayAndHour:
-            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 DX study description workload",
-                facet_col_wrap=user_profile.plotFacetColWrapVal,
-                sorting_choice=[
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                return_as_dict=return_as_dict,
-            )
-
-        if user_profile.plotDXStudyDAPvsMass:
-            parameter_dict = {
-                "df_name_col": "study_description",
-                "df_x_col": "patientstudymoduleattr__patient_weight",
-                "df_y_col": "total_dap",
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "legend_title": "Study description",
-                "colourmap": user_profile.plotColourMapChoice,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "x_axis_title": "Patient mass (kg)",
-                "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
-                "filename": "OpenREM DX study description DAP vs patient mass",
-                "return_as_dict": return_as_dict,
-            }
-            return_structure["studyDAPvsMass"] = plotly_scatter(
-                df,
-                parameter_dict,
-            )
-
-        if user_profile.plotDXRequestDAPvsMass:
-            parameter_dict = {
-                "df_name_col": "requested_procedure_code_meaning",
-                "df_x_col": "patientstudymoduleattr__patient_weight",
-                "df_y_col": "total_dap",
-                "sorting_choice": [
-                    user_profile.plotInitialSortingDirection,
-                    user_profile.plotDXInitialSortingChoice,
-                ],
-                "grouping_choice": user_profile.plotGroupingChoice,
-                "legend_title": "Requested procedure",
-                "colourmap": user_profile.plotColourMapChoice,
-                "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                "x_axis_title": "Patient mass (kg)",
-                "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
-                "filename": "OpenREM DX requested procedure DAP vs patient mass",
-                "return_as_dict": return_as_dict,
-            }
-            return_structure["requestDAPvsMass"] = plotly_scatter(
-                df,
-                parameter_dict,
-            )
-
-        if enable_standard_names:
-            charts_of_interest = [
-                user_profile.plotDXStandardStudyMeanDAP,
-                user_profile.plotDXStandardStudyFreq,
-                user_profile.plotDXStandardStudyPerDayAndHour,
-                user_profile.plotDXStandardStudyDAPvsMass,
-            ]
-
-            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.plotDXStandardStudyMeanDAP:
-
-                    name_field = "standard_names__standard_name"
-                    value_field = "total_dap"
-                    value_text = "DAP"
-                    units_text = "(cGy.cm<sup>2</sup>)"
-                    name_text = "Standard study name"
-                    variable_name_start = "standardStudy"
-                    variable_value_name = "DAP"
-                    modality_text = "DX"
-                    chart_message = ""
-
-                    new_charts = generate_average_chart_group(
-                        average_choices,
-                        chart_message,
-                        standard_name_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.plotDXInitialSortingChoice,
-                    )
-
-                    return_structure = {**return_structure, **new_charts}
-
-                if user_profile.plotDXStandardStudyFreq:
-                    parameter_dict = {
-                        "df_name_col": "standard_names__standard_name",
-                        "sorting_choice": [
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotDXInitialSortingChoice,
-                        ],
-                        "legend_title": "Standard study name",
-                        "df_x_axis_col": "x_ray_system_name",
-                        "x_axis_title": "System",
-                        "grouping_choice": user_profile.plotGroupingChoice,
-                        "colourmap": user_profile.plotColourMapChoice,
-                        "filename": "OpenREM DX standard study name frequency",
-                        "groupby_cols": None,
-                        "facet_col": None,
-                        "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                        "return_as_dict": return_as_dict,
-                    }
-                    (
-                        return_structure["standardStudyFrequencyData"],
-                        return_structure["standardStudyFrequencyDataCSV"],
-                    ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
-                        standard_name_df,
-                        parameter_dict,
-                        csv_name="standardStudyFrequencyData.csv",
-                    )
-
-                if user_profile.plotDXStandardStudyPerDayAndHour:
-                    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 DX standard study name workload",
-                        facet_col_wrap=user_profile.plotFacetColWrapVal,
-                        sorting_choice=[
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotDXInitialSortingChoice,
-                        ],
-                        return_as_dict=return_as_dict,
-                    )
-
-                if user_profile.plotDXStandardStudyDAPvsMass:
-                    parameter_dict = {
-                        "df_name_col": "standard_names__standard_name",
-                        "df_x_col": "patientstudymoduleattr__patient_weight",
-                        "df_y_col": "total_dap",
-                        "sorting_choice": [
-                            user_profile.plotInitialSortingDirection,
-                            user_profile.plotDXInitialSortingChoice,
-                        ],
-                        "grouping_choice": user_profile.plotGroupingChoice,
-                        "legend_title": "Standard study name",
-                        "colourmap": user_profile.plotColourMapChoice,
-                        "facet_col_wrap": user_profile.plotFacetColWrapVal,
-                        "x_axis_title": "Patient mass (kg)",
-                        "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
-                        "filename": "OpenREM DX standard study name DAP vs patient mass",
-                        "return_as_dict": return_as_dict,
-                    }
-                    return_structure["standardStudyDAPvsMass"] = plotly_scatter(
-                        standard_name_df,
-                        parameter_dict,
-                    )
-
-    return return_structure
-
-
-def dx_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 = DXChartOptionsFormIncStandard(request.GET)
-    else:
-        chart_options_form = DXChartOptionsForm(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_dx_chart_options(chart_options_form, user_profile)
-
-            user_profile.save()
-
-        # If submit was not clicked then use the settings already stored in the user's profile
-        else:
-            average_choices = required_average_choices(user_profile)
-
-            dx_form_data = initialise_dx_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, **dx_form_data}
-
-            chart_options_form = None
-            if enable_standard_names:
-                chart_options_form = DXChartOptionsFormIncStandard(form_data)
-            else:
-                chart_options_form = DXChartOptionsForm(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 (
+    DXChartOptionsForm,
+    DXChartOptionsFormIncStandard,
+)
+from remapp.interface.mod_filters import dx_acq_filter
+from remapp.models import (
+    create_user_profile,
+    StandardNameSettings,
+)
+from remapp.views_admin import (
+    set_average_chart_options,
+    required_average_choices,
+    initialise_dx_form_data,
+    set_dx_chart_options,
+    set_common_chart_options,
+)
+from .interface.chart_functions import (
+    create_dataframe,
+    create_dataframe_weekdays,
+    create_dataframe_aggregates,
+    plotly_boxplot,
+    plotly_barchart,
+    plotly_histogram_barchart,
+    plotly_barchart_weekdays,
+    plotly_set_default_theme,
+    plotly_frequency_barchart,
+    plotly_scatter,
+    construct_over_time_charts,
+    generate_average_chart_group,
+)
+from .tools.check_standard_name_status import are_standard_names_enabled
+
+logger = logging.getLogger(__name__)
+
+
+def generate_required_dx_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.plotDXAcquisitionMeanDAPOverTime,
+        profile.plotDXAcquisitionMeanmAsOverTime,
+        profile.plotDXAcquisitionMeankVpOverTime,
+    ]
+    if enable_standard_names:
+        charts_of_interest.append(profile.plotDXStandardAcquisitionMeanDAPOverTime)
+        charts_of_interest.append(profile.plotDXStandardAcquisitionMeanmAsOverTime)
+        charts_of_interest.append(profile.plotDXStandardAcquisitionMeankVpOverTime)
+
+    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.plotDXAcquisitionMeanDAPOverTimePeriod)]
+        ).lower()
+
+    if profile.plotDXAcquisitionMeanDAP:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of mean DAP for each acquisition protocol",
+                    "var_name": "acquisitionMeanDAP",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of median DAP for each acquisition protocol",
+                    "var_name": "acquisitionMedianDAP",
+                }
+            )
+        if profile.plotBoxplots:
+            required_charts.append(
+                {
+                    "title": "Boxplot of DAP for each acquisition protocol",
+                    "var_name": "acquisitionBoxplotDAP",
+                }
+            )
+        if profile.plotHistograms:
+            required_charts.append(
+                {
+                    "title": "Histogram of DAP for each acquisition protocol",
+                    "var_name": "acquisitionHistogramDAP",
+                }
+            )
+
+    if profile.plotDXAcquisitionFreq:
+        required_charts.append(
+            {
+                "title": "Chart of acquisition protocol frequency",
+                "var_name": "acquisitionFrequency",
+            }
+        )
+
+    if profile.plotDXAcquisitionMeankVp:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of mean kVp for each acquisition protocol",
+                    "var_name": "acquisitionMeankVp",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of median kVp for each acquisition protocol",
+                    "var_name": "acquisitionMediankVp",
+                }
+            )
+        if profile.plotBoxplots:
+            required_charts.append(
+                {
+                    "title": "Boxplot of kVp for each acquisition protocol",
+                    "var_name": "acquisitionBoxplotkVp",
+                }
+            )
+        if profile.plotHistograms:
+            required_charts.append(
+                {
+                    "title": "Histogram of kVp for each acquisition protocol",
+                    "var_name": "acquisitionHistogramkVp",
+                }
+            )
+
+    if profile.plotDXAcquisitionMeanmAs:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of mean mAs for each acquisition protocol",
+                    "var_name": "acquisitionMeanmAs",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of median mAs for each acquisition protocol",
+                    "var_name": "acquisitionMedianmAs",
+                }
+            )
+        if profile.plotBoxplots:
+            required_charts.append(
+                {
+                    "title": "Boxplot of mAs for each acquisition protocol",
+                    "var_name": "acquisitionBoxplotmAs",
+                }
+            )
+        if profile.plotHistograms:
+            required_charts.append(
+                {
+                    "title": "Histogram of mAs for each acquisition protocol",
+                    "var_name": "acquisitionHistogrammAs",
+                }
+            )
+
+    if profile.plotDXAcquisitionMeanDAPOverTime:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of mean DAP per acquisition protocol over time ("
+                    + time_period
+                    + ")",
+                    "var_name": "acquisitionMeanDAPOverTime",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of median DAP per acquisition protocol over time ("
+                    + time_period
+                    + ")",
+                    "var_name": "acquisitionMedianDAPOverTime",
+                }
+            )
+
+    if profile.plotDXAcquisitionMeankVpOverTime:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of mean kVp per acquisition protocol over time ("
+                    + time_period
+                    + ")",
+                    "var_name": "acquisitionMeankVpOverTime",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of median kVp per acquisition protocol over time ("
+                    + time_period
+                    + ")",
+                    "var_name": "acquisitionMediankVpOverTime",
+                }
+            )
+
+    if profile.plotDXAcquisitionMeanmAsOverTime:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of mean mAs per acquisition protocol over time ("
+                    + time_period
+                    + ")",
+                    "var_name": "acquisitionMeanmAsOverTime",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of median mAs per acquisition protocol over time ("
+                    + time_period
+                    + ")",
+                    "var_name": "acquisitionMedianmAsOverTime",
+                }
+            )
+
+    if profile.plotDXAcquisitionDAPvsMass:
+        required_charts.append(
+            {
+                "title": "Chart of acquisition protocol DAP vs patient mass",
+                "var_name": "acquisitionDAPvsMass",
+            }
+        )
+
+    if enable_standard_names:
+        if profile.plotDXStandardAcquisitionMeanDAP:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of mean DAP for each standard acquisition name",
+                        "var_name": "standardAcquisitionMeanDAP",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of median DAP for each standard acquisition name",
+                        "var_name": "standardAcquisitionMedianDAP",
+                    }
+                )
+            if profile.plotBoxplots:
+                required_charts.append(
+                    {
+                        "title": "Boxplot of DAP for each standard acquisition name",
+                        "var_name": "standardAcquisitionBoxplotDAP",
+                    }
+                )
+            if profile.plotHistograms:
+                required_charts.append(
+                    {
+                        "title": "Histogram of DAP for each standard acquisition name",
+                        "var_name": "standardAcquisitionHistogramDAP",
+                    }
+                )
+
+        if profile.plotDXStandardAcquisitionFreq:
+            required_charts.append(
+                {
+                    "title": "Chart of standard acquisition name frequency",
+                    "var_name": "standardAcquisitionFrequency",
+                }
+            )
+
+        if profile.plotDXStandardAcquisitionMeankVp:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of mean kVp for each standard acquisition name",
+                        "var_name": "standardAcquisitionMeankVp",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of median kVp for each standard acquisition name",
+                        "var_name": "standardAcquisitionMediankVp",
+                    }
+                )
+            if profile.plotBoxplots:
+                required_charts.append(
+                    {
+                        "title": "Boxplot of kVp for each standard acquisition name",
+                        "var_name": "standardAcquisitionBoxplotkVp",
+                    }
+                )
+            if profile.plotHistograms:
+                required_charts.append(
+                    {
+                        "title": "Histogram of kVp for each standard acquisition name",
+                        "var_name": "standardAcquisitionHistogramkVp",
+                    }
+                )
+
+        if profile.plotDXStandardAcquisitionMeanmAs:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of mean mAs for each standard acquisition name",
+                        "var_name": "standardAcquisitionMeanmAs",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of median mAs for each standard acquisition name",
+                        "var_name": "standardAcquisitionMedianmAs",
+                    }
+                )
+            if profile.plotBoxplots:
+                required_charts.append(
+                    {
+                        "title": "Boxplot of mAs for each standard acquisition name",
+                        "var_name": "standardAcquisitionBoxplotmAs",
+                    }
+                )
+            if profile.plotHistograms:
+                required_charts.append(
+                    {
+                        "title": "Histogram of mAs for each standard acquisition name",
+                        "var_name": "standardAcquisitionHistogrammAs",
+                    }
+                )
+
+        if profile.plotDXStandardAcquisitionMeanDAPOverTime:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of mean DAP per standard acquisition name over time ("
+                        + time_period
+                        + ")",
+                        "var_name": "standardAcquisitionMeanDAPOverTime",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of median DAP per standard acquisition name over time ("
+                        + time_period
+                        + ")",
+                        "var_name": "standardAcquisitionMedianDAPOverTime",
+                    }
+                )
+
+        if profile.plotDXStandardAcquisitionMeankVpOverTime:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of mean kVp per standard acquisition name over time ("
+                        + time_period
+                        + ")",
+                        "var_name": "standardAcquisitionMeankVpOverTime",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of median kVp per standard acquisition name over time ("
+                        + time_period
+                        + ")",
+                        "var_name": "standardAcquisitionMediankVpOverTime",
+                    }
+                )
+
+        if profile.plotDXStandardAcquisitionMeanmAsOverTime:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of mean mAs per standard acquisition name over time ("
+                        + time_period
+                        + ")",
+                        "var_name": "standardAcquisitionMeanmAsOverTime",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of median mAs per standard acquisition name over time ("
+                        + time_period
+                        + ")",
+                        "var_name": "standardAcquisitionMedianmAsOverTime",
+                    }
+                )
+
+        if profile.plotDXStandardAcquisitionDAPvsMass:
+            required_charts.append(
+                {
+                    "title": "Chart of standard acquisition name DAP vs patient mass",
+                    "var_name": "standardAcquisitionDAPvsMass",
+                }
+            )
+
+        if profile.plotDXStandardStudyMeanDAP:
+            if profile.plotMean:
+                required_charts.append(
+                    {
+                        "title": "Chart of mean DAP for each standard study name",
+                        "var_name": "standardStudyMeanDAP",
+                    }
+                )
+            if profile.plotMedian:
+                required_charts.append(
+                    {
+                        "title": "Chart of median DAP for each standard study name",
+                        "var_name": "standardStudyMedianDAP",
+                    }
+                )
+            if profile.plotBoxplots:
+                required_charts.append(
+                    {
+                        "title": "Boxplot of DAP for each standard study name",
+                        "var_name": "standardStudyBoxplotDAP",
+                    }
+                )
+            if profile.plotHistograms:
+                required_charts.append(
+                    {
+                        "title": "Histogram of DAP for each standard study name",
+                        "var_name": "standardStudyHistogramDAP",
+                    }
+                )
+
+        if profile.plotDXStandardStudyFreq:
+            required_charts.append(
+                {
+                    "title": "Chart of standard study name frequency",
+                    "var_name": "standardStudyFrequency",
+                }
+            )
+
+        if profile.plotDXStandardStudyPerDayAndHour:
+            required_charts.append(
+                {
+                    "title": "Chart of standard study name workload",
+                    "var_name": "standardStudyWorkload",
+                }
+            )
+
+        if profile.plotDXStandardStudyDAPvsMass:
+            required_charts.append(
+                {
+                    "title": "Chart of standard study name DAP vs patient mass",
+                    "var_name": "standardStudyDAPvsMass",
+                }
+            )
+
+    if profile.plotDXStudyMeanDAP:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of mean DAP for each study description",
+                    "var_name": "studyMeanDAP",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of median DAP for each study description",
+                    "var_name": "studyMedianDAP",
+                }
+            )
+        if profile.plotBoxplots:
+            required_charts.append(
+                {
+                    "title": "Boxplot of DAP for each study description",
+                    "var_name": "studyBoxplotDAP",
+                }
+            )
+        if profile.plotHistograms:
+            required_charts.append(
+                {
+                    "title": "Histogram of DAP for each study description",
+                    "var_name": "studyHistogramDAP",
+                }
+            )
+
+    if profile.plotDXStudyFreq:
+        required_charts.append(
+            {
+                "title": "Chart of study description frequency",
+                "var_name": "studyFrequency",
+            }
+        )
+
+    if profile.plotDXStudyPerDayAndHour:
+        required_charts.append(
+            {
+                "title": "Chart of study description workload",
+                "var_name": "studyWorkload",
+            }
+        )
+
+    if profile.plotDXStudyDAPvsMass:
+        required_charts.append(
+            {
+                "title": "Chart of study description DAP vs patient mass",
+                "var_name": "studyDAPvsMass",
+            }
+        )
+
+    if profile.plotDXRequestMeanDAP:
+        if profile.plotMean:
+            required_charts.append(
+                {
+                    "title": "Chart of mean DAP for each requested procedure",
+                    "var_name": "requestMeanDAP",
+                }
+            )
+        if profile.plotMedian:
+            required_charts.append(
+                {
+                    "title": "Chart of median DAP for each requested procedure",
+                    "var_name": "requestMedianDAP",
+                }
+            )
+        if profile.plotBoxplots:
+            required_charts.append(
+                {
+                    "title": "Boxplot of DAP for each requested procedure",
+                    "var_name": "requestBoxplotDAP",
+                }
+            )
+        if profile.plotHistograms:
+            required_charts.append(
+                {
+                    "title": "Histogram of DAP for each requested procedure",
+                    "var_name": "requestHistogramDAP",
+                }
+            )
+
+    if profile.plotDXRequestFreq:
+        required_charts.append(
+            {
+                "title": "Chart of requested procedure frequency",
+                "var_name": "requestFrequency",
+            }
+        )
+
+    if profile.plotDXRequestDAPvsMass:
+        required_charts.append(
+            {
+                "title": "Chart of requested procedure DAP vs patient mass",
+                "var_name": "requestDAPvsMass",
+            }
+        )
+
+    return required_charts
+
+
+@login_required
+def dx_summary_chart_data(request):
+    """Obtain data for Ajax chart call."""
+    pid = bool(request.user.groups.filter(name="pidgroup"))
+    f = dx_acq_filter(request.GET, pid=pid)
+
+    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 = dx_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 dx_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 radiographic 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.plotDXAcquisitionMeanDAPOverTime,
+        user_profile.plotDXAcquisitionMeankVpOverTime,
+        user_profile.plotDXAcquisitionMeanmAsOverTime,
+    ]
+    if enable_standard_names:
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanDAPOverTime)
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeankVpOverTime)
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanmAsOverTime)
+    if any(charts_of_interest):
+        plot_timeunit_period = user_profile.plotDXAcquisitionMeanDAPOverTimePeriod
+
+    #######################################################################
+    # Prepare acquisition-level Pandas DataFrame to use for charts
+    charts_of_interest = [
+        user_profile.plotDXAcquisitionMeanDAP,
+        user_profile.plotDXAcquisitionFreq,
+        user_profile.plotDXAcquisitionMeankVp,
+        user_profile.plotDXAcquisitionMeanmAs,
+        user_profile.plotDXAcquisitionMeankVpOverTime,
+        user_profile.plotDXAcquisitionMeanmAsOverTime,
+        user_profile.plotDXAcquisitionMeanDAPOverTime,
+        user_profile.plotDXAcquisitionDAPvsMass,
+    ]
+    if enable_standard_names:
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanDAP)
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionFreq)
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeankVp)
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanmAs)
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeankVpOverTime)
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanmAsOverTime)
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanDAPOverTime)
+        charts_of_interest.append(user_profile.plotDXStandardAcquisitionDAPvsMass)
+
+    if any(charts_of_interest):
+
+        name_fields = []
+        charts_of_interest = [
+            user_profile.plotDXAcquisitionMeanDAP,
+            user_profile.plotDXAcquisitionFreq,
+            user_profile.plotDXAcquisitionMeankVp,
+            user_profile.plotDXAcquisitionMeanmAs,
+            user_profile.plotDXAcquisitionMeankVpOverTime,
+            user_profile.plotDXAcquisitionMeanmAsOverTime,
+            user_profile.plotDXAcquisitionMeanDAPOverTime,
+            user_profile.plotDXAcquisitionDAPvsMass,
+        ]
+        if any(charts_of_interest):
+            name_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+            )
+
+        if enable_standard_names:
+            charts_of_interest = [
+                user_profile.plotDXStandardAcquisitionMeanDAP,
+                user_profile.plotDXStandardAcquisitionFreq,
+                user_profile.plotDXStandardAcquisitionMeankVp,
+                user_profile.plotDXStandardAcquisitionMeanmAs,
+                user_profile.plotDXStandardAcquisitionMeankVpOverTime,
+                user_profile.plotDXStandardAcquisitionMeanmAsOverTime,
+                user_profile.plotDXStandardAcquisitionMeanDAPOverTime,
+                user_profile.plotDXStandardAcquisitionDAPvsMass,
+            ]
+            if any(charts_of_interest):
+                name_fields.append(
+                    "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                )
+
+        value_fields = []
+        value_multipliers = []
+        charts_of_interest = [
+            user_profile.plotDXAcquisitionMeanDAP,
+            user_profile.plotDXAcquisitionMeanDAPOverTime,
+            user_profile.plotDXAcquisitionDAPvsMass,
+        ]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanDAP)
+            charts_of_interest.append(
+                user_profile.plotDXStandardAcquisitionMeanDAPOverTime
+            )
+            charts_of_interest.append(user_profile.plotDXStandardAcquisitionDAPvsMass)
+        if any(charts_of_interest):
+            value_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__dose_area_product"
+            )
+            value_multipliers.append(1000000)
+
+        charts_of_interest = [
+            user_profile.plotDXAcquisitionMeankVp,
+            user_profile.plotDXAcquisitionMeankVpOverTime,
+        ]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeankVp)
+            charts_of_interest.append(
+                user_profile.plotDXStandardAcquisitionMeankVpOverTime
+            )
+        if any(charts_of_interest):
+            value_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp"
+            )
+            value_multipliers.append(1)
+
+        charts_of_interest = [
+            user_profile.plotDXAcquisitionMeanmAs,
+            user_profile.plotDXAcquisitionMeanmAsOverTime,
+        ]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotDXStandardAcquisitionMeanmAs)
+            charts_of_interest.append(
+                user_profile.plotDXStandardAcquisitionMeanmAsOverTime
+            )
+        if any(charts_of_interest):
+            value_fields.append(
+                "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure"
+            )
+            value_multipliers.append(0.001)
+
+        charts_of_interest = [user_profile.plotDXAcquisitionDAPvsMass]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotDXStandardAcquisitionDAPvsMass)
+        if any(charts_of_interest):
+            value_fields.append("patientstudymoduleattr__patient_weight")
+            value_multipliers.append(1)
+
+        time_fields = []
+        date_fields = []
+        charts_of_interest = [
+            user_profile.plotDXAcquisitionMeanDAPOverTime,
+            user_profile.plotDXAcquisitionMeankVpOverTime,
+            user_profile.plotDXAcquisitionMeanmAsOverTime,
+        ]
+        if enable_standard_names:
+            charts_of_interest.append(
+                user_profile.plotDXStandardAcquisitionMeanDAPOverTime
+            )
+            charts_of_interest.append(
+                user_profile.plotDXStandardAcquisitionMeankVpOverTime
+            )
+            charts_of_interest.append(
+                user_profile.plotDXStandardAcquisitionMeanmAsOverTime
+            )
+        if any(charts_of_interest):
+            date_fields.append("study_date")
+
+        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 event DAP,
+        # and standard acquisition name filters.
+        if "event_dap_min" in f.data:
+            if f.data["event_dap_min"]:
+                df = df[
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__dose_area_product"
+                    ]
+                    >= float(f.data["event_dap_min"])
+                ]
+
+        if "event_dap_max" in f.data:
+            if f.data["event_dap_max"]:
+                df = df[
+                    df[
+                        "projectionxrayradiationdose__irradeventxraydata__dose_area_product"
+                    ]
+                    <= float(f.data["event_dap_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 user_profile.plotDXAcquisitionMeanDAP:
+
+            name_field = (
+                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+            )
+            value_field = (
+                "projectionxrayradiationdose__irradeventxraydata__dose_area_product"
+            )
+            value_text = "DAP"
+            units_text = "(cGy.cm<sup>2</sup>)"
+            name_text = "Acquisition protocol"
+            variable_name_start = "acquisition"
+            variable_value_name = "DAP"
+            modality_text = "DX"
+            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.plotDXInitialSortingChoice,
+            )
+
+            return_structure = {**return_structure, **new_charts}
+
+        if user_profile.plotDXAcquisitionMeankVp:
+
+            name_field = (
+                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+            )
+            value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp"
+            value_text = "kVp"
+            units_text = ""
+            name_text = "Acquisition protocol"
+            variable_name_start = "acquisition"
+            variable_value_name = "kVp"
+            modality_text = "DX"
+            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.plotDXInitialSortingChoice,
+            )
+
+            return_structure = {**return_structure, **new_charts}
+
+        if user_profile.plotDXAcquisitionMeanmAs:
+
+            name_field = (
+                "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol"
+            )
+            value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure"  # pylint: disable=line-too-long
+            value_text = "mAs"
+            units_text = ""
+            name_text = "Acquisition protocol"
+            variable_name_start = "acquisition"
+            variable_value_name = "mAs"
+            modality_text = "DX"
+            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.plotDXInitialSortingChoice,
+            )
+
+            return_structure = {**return_structure, **new_charts}
+
+        if user_profile.plotDXAcquisitionFreq:
+            parameter_dict = {
+                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                "legend_title": "Acquisition protocol",
+                "df_x_axis_col": "x_ray_system_name",
+                "x_axis_title": "System",
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "colourmap": user_profile.plotColourMapChoice,
+                "filename": "OpenREM DX 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.plotDXAcquisitionMeanDAPOverTime:
+            facet_title = "System"
+
+            if user_profile.plotGroupingChoice == "series":
+                facet_title = "Acquisition protocol"
+
+            parameter_dict = {
+                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
+                "df_value_col": "projectionxrayradiationdose__irradeventxraydata__dose_area_product",
+                "df_date_col": "study_date",
+                "name_title": "Acquisition protocol",
+                "value_title": "DAP (cGy.cm<sup>2</sup>)",
+                "date_title": "Study date",
+                "facet_title": facet_title,
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                "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 DX acquisition protocol DAP over time",
+                "return_as_dict": return_as_dict,
+            }
+            result = construct_over_time_charts(
+                df,
+                parameter_dict,
+            )
+
+            if user_profile.plotMean:
+                return_structure["acquisitionMeanDAPOverTime"] = result["mean"]
+            if user_profile.plotMedian:
+                return_structure["acquisitionMedianDAPOverTime"] = result["median"]
+
+        if user_profile.plotDXAcquisitionMeankVpOverTime:
+            facet_title = "System"
+
+            if user_profile.plotGroupingChoice == "series":
+                facet_title = "Acquisition protocol"
+
+            parameter_dict = {
+                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
+                "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp",
+                "df_date_col": "study_date",
+                "name_title": "Acquisition protocol",
+                "value_title": "kVp",
+                "date_title": "Study date",
+                "facet_title": facet_title,
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                "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 DX acquisition protocol kVp over time",
+                "return_as_dict": return_as_dict,
+            }
+            result = construct_over_time_charts(
+                df,
+                parameter_dict,
+            )
+
+            if user_profile.plotMean:
+                return_structure["acquisitionMeankVpOverTime"] = result["mean"]
+            if user_profile.plotMedian:
+                return_structure["acquisitionMediankVpOverTime"] = result["median"]
+
+        if user_profile.plotDXAcquisitionMeanmAsOverTime:
+            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__exposure__exposure",  # pylint: disable=line-too-long
+                "df_date_col": "study_date",
+                "name_title": "Acquisition protocol",
+                "value_title": "mAs",
+                "date_title": "Study date",
+                "facet_title": facet_title,
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                "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 DX acquisition protocol mAs over time",
+                "return_as_dict": return_as_dict,
+            }
+            result = construct_over_time_charts(
+                df,
+                parameter_dict,
+            )
+
+            if user_profile.plotMean:
+                return_structure["acquisitionMeanmAsOverTime"] = result["mean"]
+            if user_profile.plotMedian:
+                return_structure["acquisitionMedianmAsOverTime"] = result["median"]
+
+        if user_profile.plotDXAcquisitionDAPvsMass:
+            parameter_dict = {
+                "df_name_col": "projectionxrayradiationdose__irradeventxraydata__acquisition_protocol",
+                "df_x_col": "patientstudymoduleattr__patient_weight",
+                "df_y_col": "projectionxrayradiationdose__irradeventxraydata__dose_area_product",
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "legend_title": "Acquisition protocol",
+                "colourmap": user_profile.plotColourMapChoice,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "x_axis_title": "Patient mass (kg)",
+                "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
+                "filename": "OpenREM DX acquisition protocol DAP vs patient mass",
+                "return_as_dict": return_as_dict,
+            }
+            return_structure["acquisitionDAPvsMass"] = plotly_scatter(
+                df,
+                parameter_dict,
+            )
+
+        if enable_standard_names:
+            charts_of_interest = [
+                user_profile.plotDXStandardAcquisitionMeanDAP,
+                user_profile.plotDXStandardAcquisitionFreq,
+                user_profile.plotDXStandardAcquisitionMeankVp,
+                user_profile.plotDXStandardAcquisitionMeanmAs,
+                user_profile.plotDXStandardAcquisitionMeankVpOverTime,
+                user_profile.plotDXStandardAcquisitionMeanmAsOverTime,
+                user_profile.plotDXStandardAcquisitionMeanDAPOverTime,
+                user_profile.plotDXStandardAcquisitionDAPvsMass,
+            ]
+
+            if any(charts_of_interest):
+
+                # Exclude "Blank" and "blank" standard_acqusition_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.plotDXStandardAcquisitionMeanDAP:
+                    name_field = "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                    value_field = "projectionxrayradiationdose__irradeventxraydata__dose_area_product"
+                    value_text = "DAP"
+                    units_text = "(cGy.cm<sup>2</sup>)"
+                    name_text = "Standard acquisition name"
+                    variable_name_start = "standardAcquisition"
+                    variable_value_name = "DAP"
+                    modality_text = "DX"
+                    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.plotDXInitialSortingChoice,
+                    )
+
+                    return_structure = {**return_structure, **new_charts}
+
+                if user_profile.plotDXStandardAcquisitionMeankVp:
+                    name_field = "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                    value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp"
+                    value_text = "kVp"
+                    units_text = ""
+                    name_text = "Standard acquisition name"
+                    variable_name_start = "standardAcquisition"
+                    variable_value_name = "kVp"
+                    modality_text = "DX"
+                    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.plotDXInitialSortingChoice,
+                    )
+
+                    return_structure = {**return_structure, **new_charts}
+
+                if user_profile.plotDXStandardAcquisitionMeanmAs:
+                    name_field = "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name"
+                    value_field = "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure"  # pylint: disable=line-too-long
+                    value_text = "mAs"
+                    units_text = ""
+                    name_text = "Standard acquisition name"
+                    variable_name_start = "standardAcquisition"
+                    variable_value_name = "mAs"
+                    modality_text = "DX"
+                    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.plotDXInitialSortingChoice,
+                    )
+
+                    return_structure = {**return_structure, **new_charts}
+
+                if user_profile.plotDXStandardAcquisitionFreq:
+                    parameter_dict = {
+                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
+                        "sorting_choice": [
+                            user_profile.plotInitialSortingDirection,
+                            user_profile.plotDXInitialSortingChoice,
+                        ],
+                        "legend_title": "Standard acquisition name",
+                        "df_x_axis_col": "x_ray_system_name",
+                        "x_axis_title": "System",
+                        "grouping_choice": user_profile.plotGroupingChoice,
+                        "colourmap": user_profile.plotColourMapChoice,
+                        "filename": "OpenREM DX 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(
+                        df_without_blanks,
+                        parameter_dict,
+                        csv_name="standardAcquisitionFrequencyData.csv",
+                    )
+
+                if user_profile.plotDXStandardAcquisitionMeanDAPOverTime:
+                    facet_title = "System"
+
+                    if user_profile.plotGroupingChoice == "series":
+                        facet_title = "Standard acquisition name"
+
+                    parameter_dict = {
+                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
+                        "df_value_col": "projectionxrayradiationdose__irradeventxraydata__dose_area_product",
+                        "df_date_col": "study_date",
+                        "name_title": "Standard acquisition name",
+                        "value_title": "DAP (cGy.cm<sup>2</sup>)",
+                        "date_title": "Study date",
+                        "facet_title": facet_title,
+                        "sorting_choice": [
+                            user_profile.plotInitialSortingDirection,
+                            user_profile.plotDXInitialSortingChoice,
+                        ],
+                        "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 DX standard acquisition name DAP over time",
+                        "return_as_dict": return_as_dict,
+                    }
+                    result = construct_over_time_charts(
+                        df_without_blanks,
+                        parameter_dict,
+                    )
+
+                    if user_profile.plotMean:
+                        return_structure["standardAcquisitionMeanDAPOverTime"] = result[
+                            "mean"
+                        ]
+                    if user_profile.plotMedian:
+                        return_structure["standardAcquisitionMedianDAPOverTime"] = (
+                            result["median"]
+                        )
+
+                if user_profile.plotDXStandardAcquisitionMeankVpOverTime:
+                    facet_title = "System"
+
+                    if user_profile.plotGroupingChoice == "series":
+                        facet_title = "Standard acquisition name"
+
+                    parameter_dict = {
+                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
+                        "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__kvp__kvp",  # pylint: disable=line-too-long
+                        "df_date_col": "study_date",
+                        "name_title": "Standard acquisition name",
+                        "value_title": "kVp",
+                        "date_title": "Study date",
+                        "facet_title": facet_title,
+                        "sorting_choice": [
+                            user_profile.plotInitialSortingDirection,
+                            user_profile.plotDXInitialSortingChoice,
+                        ],
+                        "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 DX standard acquisition name kVp over time",
+                        "return_as_dict": return_as_dict,
+                    }
+                    result = construct_over_time_charts(
+                        df_without_blanks,
+                        parameter_dict,
+                    )
+
+                    if user_profile.plotMean:
+                        return_structure["standardAcquisitionMeankVpOverTime"] = result[
+                            "mean"
+                        ]
+                    if user_profile.plotMedian:
+                        return_structure["standardAcquisitionMediankVpOverTime"] = (
+                            result["median"]
+                        )
+
+                if user_profile.plotDXStandardAcquisitionMeanmAsOverTime:
+                    facet_title = "System"
+
+                    if user_profile.plotGroupingChoice == "series":
+                        facet_title = "Standard acquisition name"
+
+                    parameter_dict = {
+                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
+                        "df_value_col": "projectionxrayradiationdose__irradeventxraydata__irradeventxraysourcedata__exposure__exposure",  # pylint: disable=line-too-long
+                        "df_date_col": "study_date",
+                        "name_title": "Standard acquisition name",
+                        "value_title": "mAs",
+                        "date_title": "Study date",
+                        "facet_title": facet_title,
+                        "sorting_choice": [
+                            user_profile.plotInitialSortingDirection,
+                            user_profile.plotDXInitialSortingChoice,
+                        ],
+                        "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 DX standard acquisition name mAs over time",
+                        "return_as_dict": return_as_dict,
+                    }
+                    result = construct_over_time_charts(
+                        df_without_blanks,
+                        parameter_dict,
+                    )
+
+                    if user_profile.plotMean:
+                        return_structure["standardAcquisitionMeanmAsOverTime"] = result[
+                            "mean"
+                        ]
+                    if user_profile.plotMedian:
+                        return_structure["standardAcquisitionMedianmAsOverTime"] = (
+                            result["median"]
+                        )
+
+                if user_profile.plotDXStandardAcquisitionDAPvsMass:
+                    parameter_dict = {
+                        "df_name_col": "projectionxrayradiationdose__irradeventxraydata__standard_protocols__standard_name",  # pylint: disable=line-too-long
+                        "df_x_col": "patientstudymoduleattr__patient_weight",
+                        "df_y_col": "projectionxrayradiationdose__irradeventxraydata__dose_area_product",
+                        "sorting_choice": [
+                            user_profile.plotInitialSortingDirection,
+                            user_profile.plotDXInitialSortingChoice,
+                        ],
+                        "grouping_choice": user_profile.plotGroupingChoice,
+                        "legend_title": "Standard acquisition name",
+                        "colourmap": user_profile.plotColourMapChoice,
+                        "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                        "x_axis_title": "Patient mass (kg)",
+                        "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
+                        "filename": "OpenREM DX standard acquisition name DAP vs patient mass",
+                        "return_as_dict": return_as_dict,
+                    }
+                    return_structure["standardAcquisitionDAPvsMass"] = plotly_scatter(
+                        df_without_blanks,
+                        parameter_dict,
+                    )
+
+    #######################################################################
+    # Prepare study- and request-level Pandas DataFrame to use for charts
+    charts_of_interest = [
+        user_profile.plotDXStudyMeanDAP,
+        user_profile.plotDXStudyFreq,
+        user_profile.plotDXStudyPerDayAndHour,
+        user_profile.plotDXStudyDAPvsMass,
+        user_profile.plotDXRequestMeanDAP,
+        user_profile.plotDXRequestFreq,
+        user_profile.plotDXRequestDAPvsMass,
+    ]
+    if enable_standard_names:
+        charts_of_interest.append(user_profile.plotDXStandardStudyMeanDAP)
+        charts_of_interest.append(user_profile.plotDXStandardStudyFreq)
+        charts_of_interest.append(user_profile.plotDXStandardStudyPerDayAndHour)
+        charts_of_interest.append(user_profile.plotDXStandardStudyDAPvsMass)
+
+    if any(charts_of_interest):
+
+        name_fields = []
+        charts_of_interest = [
+            user_profile.plotDXStudyMeanDAP,
+            user_profile.plotDXStudyFreq,
+            user_profile.plotDXStudyPerDayAndHour,
+            user_profile.plotDXStudyDAPvsMass,
+        ]
+        if any(charts_of_interest):
+            name_fields.append("study_description")
+
+        charts_of_interest = [
+            user_profile.plotDXRequestMeanDAP,
+            user_profile.plotDXRequestFreq,
+            user_profile.plotDXRequestDAPvsMass,
+        ]
+        if any(charts_of_interest):
+            name_fields.append("requested_procedure_code_meaning")
+
+        if enable_standard_names:
+            charts_of_interest = [
+                user_profile.plotDXStandardStudyMeanDAP,
+                user_profile.plotDXStandardStudyFreq,
+                user_profile.plotDXStandardStudyPerDayAndHour,
+                user_profile.plotDXStandardStudyDAPvsMass,
+            ]
+            if any(charts_of_interest):
+                name_fields.append("standard_names__standard_name")
+
+        value_fields = []
+        value_multipliers = []
+        charts_of_interest = [
+            user_profile.plotDXStudyMeanDAP,
+            user_profile.plotDXRequestMeanDAP,
+            user_profile.plotDXStudyDAPvsMass,
+            user_profile.plotDXRequestDAPvsMass,
+        ]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotDXStandardStudyMeanDAP)
+            charts_of_interest.append(user_profile.plotDXStandardStudyDAPvsMass)
+
+        if any(charts_of_interest):
+            value_fields.append("total_dap")
+            value_multipliers.append(1000000)
+
+        charts_of_interest = [
+            user_profile.plotDXStudyDAPvsMass,
+            user_profile.plotDXRequestDAPvsMass,
+        ]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotDXStandardStudyDAPvsMass)
+
+        if any(charts_of_interest):
+            value_fields.append("patientstudymoduleattr__patient_weight")
+            value_multipliers.append(1)
+
+        date_fields = []
+        time_fields = []
+        charts_of_interest = [user_profile.plotDXStudyPerDayAndHour]
+        if enable_standard_names:
+            charts_of_interest.append(user_profile.plotDXStandardStudyPerDayAndHour)
+
+        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,
+            data_point_value_multipliers=value_multipliers,
+            char_wrap=user_profile.plotLabelCharWrap,
+            uid="pk",
+        )
+        #######################################################################
+
+        if user_profile.plotDXStudyMeanDAP:
+
+            name_field = "study_description"
+            value_field = "total_dap"
+            value_text = "DAP"
+            units_text = "(cGy.cm<sup>2</sup>)"
+            name_text = "Study description"
+            variable_name_start = "study"
+            variable_value_name = "DAP"
+            modality_text = "DX"
+            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.plotDXInitialSortingChoice,
+            )
+
+            return_structure = {**return_structure, **new_charts}
+
+        if user_profile.plotDXStudyFreq:
+            parameter_dict = {
+                "df_name_col": "study_description",
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                "legend_title": "Study description",
+                "df_x_axis_col": "x_ray_system_name",
+                "x_axis_title": "System",
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "colourmap": user_profile.plotColourMapChoice,
+                "filename": "OpenREM DX study description frequency",
+                "groupby_cols": None,
+                "facet_col": None,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "return_as_dict": return_as_dict,
+            }
+            (
+                return_structure["studyFrequencyData"],
+                return_structure["studyFrequencyDataCSV"],
+            ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
+                df,
+                parameter_dict,
+                csv_name="studyFrequencyData.csv",
+            )
+
+        if user_profile.plotDXRequestMeanDAP:
+
+            name_field = "requested_procedure_code_meaning"
+            value_field = "total_dap"
+            value_text = "DAP"
+            units_text = "(cGy.cm<sup>2</sup>)"
+            name_text = "Requested procedure"
+            variable_name_start = "request"
+            variable_value_name = "DAP"
+            modality_text = "DX"
+            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.plotDXInitialSortingChoice,
+            )
+
+            return_structure = {**return_structure, **new_charts}
+
+        if user_profile.plotDXRequestFreq:
+            parameter_dict = {
+                "df_name_col": "requested_procedure_code_meaning",
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                "legend_title": "Requested procedure",
+                "df_x_axis_col": "x_ray_system_name",
+                "x_axis_title": "System",
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "colourmap": user_profile.plotColourMapChoice,
+                "filename": "OpenREM DX requested procedure frequency",
+                "groupby_cols": None,
+                "facet_col": None,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "return_as_dict": return_as_dict,
+            }
+            (
+                return_structure["requestFrequencyData"],
+                return_structure["requestFrequencyDataCSV"],
+            ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
+                df,
+                parameter_dict,
+                csv_name="requestFrequencyData.csv",
+            )
+
+        if user_profile.plotDXStudyPerDayAndHour:
+            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 DX study description workload",
+                facet_col_wrap=user_profile.plotFacetColWrapVal,
+                sorting_choice=[
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                return_as_dict=return_as_dict,
+            )
+
+        if user_profile.plotDXStudyDAPvsMass:
+            parameter_dict = {
+                "df_name_col": "study_description",
+                "df_x_col": "patientstudymoduleattr__patient_weight",
+                "df_y_col": "total_dap",
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "legend_title": "Study description",
+                "colourmap": user_profile.plotColourMapChoice,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "x_axis_title": "Patient mass (kg)",
+                "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
+                "filename": "OpenREM DX study description DAP vs patient mass",
+                "return_as_dict": return_as_dict,
+            }
+            return_structure["studyDAPvsMass"] = plotly_scatter(
+                df,
+                parameter_dict,
+            )
+
+        if user_profile.plotDXRequestDAPvsMass:
+            parameter_dict = {
+                "df_name_col": "requested_procedure_code_meaning",
+                "df_x_col": "patientstudymoduleattr__patient_weight",
+                "df_y_col": "total_dap",
+                "sorting_choice": [
+                    user_profile.plotInitialSortingDirection,
+                    user_profile.plotDXInitialSortingChoice,
+                ],
+                "grouping_choice": user_profile.plotGroupingChoice,
+                "legend_title": "Requested procedure",
+                "colourmap": user_profile.plotColourMapChoice,
+                "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                "x_axis_title": "Patient mass (kg)",
+                "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
+                "filename": "OpenREM DX requested procedure DAP vs patient mass",
+                "return_as_dict": return_as_dict,
+            }
+            return_structure["requestDAPvsMass"] = plotly_scatter(
+                df,
+                parameter_dict,
+            )
+
+        if enable_standard_names:
+            charts_of_interest = [
+                user_profile.plotDXStandardStudyMeanDAP,
+                user_profile.plotDXStandardStudyFreq,
+                user_profile.plotDXStandardStudyPerDayAndHour,
+                user_profile.plotDXStandardStudyDAPvsMass,
+            ]
+
+            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.plotDXStandardStudyMeanDAP:
+
+                    name_field = "standard_names__standard_name"
+                    value_field = "total_dap"
+                    value_text = "DAP"
+                    units_text = "(cGy.cm<sup>2</sup>)"
+                    name_text = "Standard study name"
+                    variable_name_start = "standardStudy"
+                    variable_value_name = "DAP"
+                    modality_text = "DX"
+                    chart_message = ""
+
+                    new_charts = generate_average_chart_group(
+                        average_choices,
+                        chart_message,
+                        standard_name_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.plotDXInitialSortingChoice,
+                    )
+
+                    return_structure = {**return_structure, **new_charts}
+
+                if user_profile.plotDXStandardStudyFreq:
+                    parameter_dict = {
+                        "df_name_col": "standard_names__standard_name",
+                        "sorting_choice": [
+                            user_profile.plotInitialSortingDirection,
+                            user_profile.plotDXInitialSortingChoice,
+                        ],
+                        "legend_title": "Standard study name",
+                        "df_x_axis_col": "x_ray_system_name",
+                        "x_axis_title": "System",
+                        "grouping_choice": user_profile.plotGroupingChoice,
+                        "colourmap": user_profile.plotColourMapChoice,
+                        "filename": "OpenREM DX standard study name frequency",
+                        "groupby_cols": None,
+                        "facet_col": None,
+                        "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                        "return_as_dict": return_as_dict,
+                    }
+                    (
+                        return_structure["standardStudyFrequencyData"],
+                        return_structure["standardStudyFrequencyDataCSV"],
+                    ) = plotly_frequency_barchart(  # pylint: disable=line-too-long
+                        standard_name_df,
+                        parameter_dict,
+                        csv_name="standardStudyFrequencyData.csv",
+                    )
+
+                if user_profile.plotDXStandardStudyPerDayAndHour:
+                    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 DX standard study name workload",
+                            facet_col_wrap=user_profile.plotFacetColWrapVal,
+                            sorting_choice=[
+                                user_profile.plotInitialSortingDirection,
+                                user_profile.plotDXInitialSortingChoice,
+                            ],
+                            return_as_dict=return_as_dict,
+                        )
+                    )
+
+                if user_profile.plotDXStandardStudyDAPvsMass:
+                    parameter_dict = {
+                        "df_name_col": "standard_names__standard_name",
+                        "df_x_col": "patientstudymoduleattr__patient_weight",
+                        "df_y_col": "total_dap",
+                        "sorting_choice": [
+                            user_profile.plotInitialSortingDirection,
+                            user_profile.plotDXInitialSortingChoice,
+                        ],
+                        "grouping_choice": user_profile.plotGroupingChoice,
+                        "legend_title": "Standard study name",
+                        "colourmap": user_profile.plotColourMapChoice,
+                        "facet_col_wrap": user_profile.plotFacetColWrapVal,
+                        "x_axis_title": "Patient mass (kg)",
+                        "y_axis_title": "DAP (mGy.cm<sup>2</sub>)",
+                        "filename": "OpenREM DX standard study name DAP vs patient mass",
+                        "return_as_dict": return_as_dict,
+                    }
+                    return_structure["standardStudyDAPvsMass"] = plotly_scatter(
+                        standard_name_df,
+                        parameter_dict,
+                    )
+
+    return return_structure
+
+
+def dx_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 = DXChartOptionsFormIncStandard(request.GET)
+    else:
+        chart_options_form = DXChartOptionsForm(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_dx_chart_options(chart_options_form, user_profile)
+
+            user_profile.save()
+
+        # If submit was not clicked then use the settings already stored in the user's profile
+        else:
+            average_choices = required_average_choices(user_profile)
+
+            dx_form_data = initialise_dx_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, **dx_form_data}
+
+            chart_options_form = None
+            if enable_standard_names:
+                chart_options_form = DXChartOptionsFormIncStandard(form_data)
+            else:
+                chart_options_form = DXChartOptionsForm(form_data)
+
+    return chart_options_form