Adding the benchmark data for extra P-values, for HP, STIM and HP_STIM

        If mask values are provided, then any non-nan values in the mask variable will be set to zero across all variables in the returned output_data object.

    :param input_params:

        A dictionary containing parameters for use by the doublet calculation, If no input_params are provided then the values for the default thermogis scenario are used (the "BaseCase" scenario).

    :param rng_seed:

        A Integer that sets the random seed within the doublet calculation, if not provided then the random number generator picks a seed at random.

    :param p_values:

        A list of p_values for the doublet calculation to perform over; if no p_values are provided then the default value of P50 is used.

    :return output_data:

    doublet.doubletCalc1DData.setReservoirTemp(temperature)

    doublet.doubletCalc1DData.setUseHeatPump(input_params["use_heat_pump"])

    if not ["use_heat_pump"]:

    if not input_params["use_heat_pump"]:

        injectionTemp = np.max([temperature - input_params["max_cooling_temperature_range"], input_params["return_temperature"]])

    elif input_params["use_heat_pump"] and input_params["calculate_cop"] and not input_params["hp_application_mode"]:

        injectionTemp = doublet.calculateInjectionTempWithHeatPump(temperature, input_params["hp_direct_heat_input_temp"])