diff --git a/README.md b/README.md
new file mode 100644
index 0000000..fb12afa
--- /dev/null
+++ b/README.md
@@ -0,0 +1,38 @@
+Quick tutorial:
+1. Browse the .mat file you want to open. (put it into /data folder)
+2. Choose either Load (fast) or Load (Presets)
+    - Load(fast):
+        - will load ONLY the name of all saved data in your .mat file
+        - select data you want to plot in custom subplots -> select data (left side) 
+        - all features will work for this option
+        - for customizing look under the next chapter
+
+    - Load (Presets)
+        - will load a pre configured set of data in your .mat file (takes a bit of time)
+        - choose a preset and a season, then press "Plot Presets"
+            -> will plot pre configured graphs
+        - all features in "Custom Subplots" do NOT work with this
+        - the idea is that only sortdata.m and the plot_settings need to be configured each season
+        - for customizing look under the next chapter (only features within Set Timestamps panel work)
+
+Features:
+1. Set Timestamps:
+    - works for both
+    - entire stint (all data recorded) is the default option that will be plotted
+    - otherwise you can select either the detected laps (will plot, but laps might be inaccurately recorded) or custom timestamps (remember to tick the Use box)
+2. Reset Panel and To Workspace Buttons:
+    - work only for Presets
+    - should something not work try pressing reset and plot preset again
+    - to workspace throws the entire preset dataset into the matlab workspace (for further manipulation)
+        - currently limited in capabilites
+3. Custom Subplots:
+    - only works for the selected data (if you press load preset you can also use the displayed data, but somethings like cell voltages might not be avaible)
+    - columns and rows to set the number of plots you want in one figure
+    - how it works:
+![selected subplot is 1](tutorial/1.png)
+![this time a lap, and a different dataset was selected](tutorial/2.png)
+![only if you want a second -yaxis select the option "plot with 2nd yAxis](tutorial/3.png)
+    - press reset subplot to delete the selected graph
+
+not working:
+- x Axis might not be avaible for everything
diff --git a/app.mlapp b/app.mlapp
index 00be9e9..35bd298 100644
Binary files a/app.mlapp and b/app.mlapp differ
diff --git a/assets/loadvector.m b/assets/loadvector.m
deleted file mode 100644
index 6a5e5d0..0000000
--- a/assets/loadvector.m
+++ /dev/null
@@ -1,79 +0,0 @@
-function [pltData,laps_strArray] = loadvector(filepath)
-
-    % Load dataset:
-    fprintf("Loading: %s\n",filepath)
-    data = vector2timetable(filepath,"ABX*","Shunt*","AMS_Status*","INV*","XSens*");
-    fprintf("Finished loading\n------------------------\n")
-    
-    % get date / time
-    date = dateshift(data.Time(1),'start','day');
-    fprintf("Date: %s\n",date)
-    fprintf("Start: %d:%d:%02.0f\n",hour(data.Time(1)),minute(data.Time(1)),second(data.Time(1)))
-    fprintf("End: %d:%d:%02.0f\n",hour(data.Time(end)),minute(data.Time(end)),second(data.Time(end)))
-    fprintf("------------------------\n")
-    
-    % get lapcounter / laptrigger
-    laps_counter = data.ABX_Driver_ABX_Lapcounter;
-    timestamp_hms(1) = data.Time(1) - date;
-    timestamp_index(1) = 1;
-    j = 2;
-    for i = 2:length(laps_counter)
-        if laps_counter(i) ~= laps_counter(i-1)
-            timestamp_hms(j) = data.Time(i) - date;
-            timestamp_index(j) = i;
-            j = j + 1;
-        end
-    end
-
-    % calculate laptimes
-    timestamp_s = seconds(timestamp_hms);
-    for i = 1:length(timestamp_s)-1
-        laps_s(i) = timestamp_s(i+1) - timestamp_s(i);
-    end
-    
-    % if laps could be detected find fastest lap and create list of
-    % laptimes, otherwise only "Entire stint" will be selectable
-    laps_strArray(1) = sprintf("Entire stint");
-    if (exist("laps_s","var") == true)
-        % create list of laptimes
-        for i = 1:length(laps_s)
-            if laps_s(i) == min(laps_s)
-                laps_strArray(i+1) = sprintf("Lap %d: %.3fs (best)",i,laps_s(i));
-            else
-                laps_strArray(i+1) = sprintf("Lap %d: %.3fs",i,laps_s(i));
-            end
-        end
-        % display fastest lap
-        [laps_best_time, laps_best] = min(laps_s);
-        fprintf("Fastest lap: %.2fs (%d)\n",laps_best_time,laps_best)
-    end
-    
-    % create time arrays in hh:mm:ss and just seconds
-    time_hms = data.Time-date;
-    time_s = seconds(time_hms);
-    time_s = time_s - time_s(1);    % set beginning of session to 0
-
-    % calcualte energy consumption and total distance
-    power_kw = data.Shunt_Voltage1_Shunt_Voltage1.*data.Shunt_Current_Shunt_Current/1000;
-    power_regen_kw = power_kw;
-    power_regen_kw(power_regen_kw > 0) = 0;
-    power_used_kw = power_kw;
-    power_used_kw(power_used_kw < 0) = 0;
-    energy_kwh = trapz(time_s, power_kw)/3600;
-    energy_regen_kwh = trapz(time_s, power_regen_kw)/3600;
-    energy_used_kwh = trapz(time_s, power_used_kw)/3600;
-    fprintf("Energy used: %.2fkWh\n",energy_used_kwh)
-    fprintf("Energy regen: %.2fkWh\n",energy_regen_kwh)
-    fprintf("Energy total: %.2fkWh\n",energy_kwh)
-    distance_km = trapz(time_s, movmean(data.ABX_Driver_ABX_Speed,10))/1000;
-    fprintf("Distance driven: %.2fkm\n",distance_km)
-    fprintf("------------------------\n")
-
-    pltData(1) = sortdata(data,1,length(data.Time));
-    % sort data for plotting:
-    for i = 1:length(timestamp_index)-1
-        pltData(i+1) = sortdata(data,timestamp_index(i),timestamp_index(i+1));
-    end
-    
-end
-
diff --git a/assets/sortdata.m b/assets/sortdata.m
deleted file mode 100644
index dedf2af..0000000
--- a/assets/sortdata.m
+++ /dev/null
@@ -1,123 +0,0 @@
-function [pltData] = sortdata(data,start,stop)
-    date = dateshift(data.Time(1),'start','day');
-    time_hms = data.Time-date;
-    time_s = seconds(time_hms);
-    time_s = time_s - time_s(1);    % set beginning of session to 0
-
-    % x-Axis:
-    pltData.Time = data.Time(start:stop);
-    pltData.time_s = time_s(start:stop)-time_s(start);
-    pltData.time_hms = time_hms(start:stop);
-
-    % distance calculation
-    speed_mps = movmean(data.ABX_Driver_ABX_Speed(start:stop),50); % TODO: test different speed filters for distance calculation
-    pltData.distance(1) = 0;
-    for i = 2:length(speed_mps)
-        pltData.distance(i) = speed_mps(i)*(pltData.time_s(i)-pltData.time_s(i-1)) + pltData.distance(i-1);
-    end
-    pltData.xAxis = pltData.time_s;     % set time as default x-Axis for plotting
-
-    %% Misc
-    pltData.app_percent = data.ABX_Driver_ABX_APPS_percent(start:stop);
-    pltData.speed_kph = 3.6*movmean(data.ABX_Driver_ABX_Speed(start:stop),50);  % same filter as for distance calculation ???
-    pltData.steering_deg = data.ABX_Driver_ABX_Steering_Angle(start:stop);
-
-    %% AMS:
-    pltData.ams_soc = data.AMS_Status_SOC(start:stop);
-    pltData.ams_tmax = data.AMS_Status_Max_cell_temp(start:stop);
-    pltData.ams_utot = data.Shunt_Voltage1_Shunt_Voltage1(start:stop);
-    pltData.ams_itot = data.Shunt_Current_Shunt_Current(start:stop);
-    % calculations:
-    pltData.ams_ptot = pltData.ams_utot.*pltData.ams_itot/1000;
-
-    %% Brakes
-    % brake pressure
-    pltData.brakePFront_bar = data.ABX_Driver_ABX_BrakeP_F(start:stop);
-    pltData.brakePRear_bar = data.ABX_Driver_ABX_BrakeP_R(start:stop);
-    % brake disc temperatures
-    pltData.brakeTFrontLeft_degC = data.ABX_BrakeT_ABX_BrakeT_FL(start:stop);
-    pltData.brakeTFrontRight_degC = data.ABX_BrakeT_ABX_BrakeT_FR(start:stop);
-    pltData.brakeTRearLeft_degC = data.ABX_BrakeT_ABX_BrakeT_RL(start:stop);
-    pltData.brakeTRearRight_degC = data.ABX_BrakeT_ABX_BrakeT_RR(start:stop);
-    % calculate brake bias [%]
-    minBrakeP = 5;  % minimum brake pressure to avoid artifacts due to sensor noise near 0 bar
-    brakePFront = pltData.brakePFront_bar;
-    brakePFront(brakePFront < minBrakeP) = minBrakeP;
-    brakePRear = pltData.brakePRear_bar;
-    brakePRear(brakePRear < minBrakeP) = minBrakeP;
-    pltData.brakeBias_perc = 100*brakePFront./brakePRear;  % check calculation!
-
-    %% Cooling system
-
-    %% Dampers
-    % damper positions
-    pltData.damper_FL_mm = data.ABX_Dampers_ABX_Damper_FL(start:stop); %Heave_F
-    pltData.damper_FR_mm = data.ABX_Dampers_ABX_Damper_FR(start:stop); %Roll_F
-    pltData.damper_RL_mm = data.ABX_Dampers_ABX_Damper_RL(start:stop); %Heave_R
-    pltData.damper_RR_mm = data.ABX_Dampers_ABX_Damper_RR(start:stop); %Roll_R
-    % calculate damper velocities
-    pltData.velocity_FL_mmps(1) = 0;
-    pltData.velocity_FR_mmps(1) = 0;
-    pltData.velocity_RL_mmps(1) = 0;
-    pltData.velocity_RR_mmps(1) = 0;
-    for i = 2:length(pltData.time_s)
-        timestep = pltData.time_s(i)-pltData.time_s(i-1);
-        pltData.velocity_FL_mmps(i) = pltData.damper_FL_mm(i)-pltData.damper_FL_mm(i-1)/timestep;
-        pltData.velocity_FR_mmps(i) = pltData.damper_FR_mm(i)-pltData.damper_FR_mm(i-1)/timestep;
-        pltData.velocity_RL_mmps(i) = pltData.damper_RL_mm(i)-pltData.damper_RL_mm(i-1)/timestep;
-        pltData.velocity_RR_mmps(i) = pltData.damper_RR_mm(i)-pltData.damper_RR_mm(i-1)/timestep;
-    end
-    % filter damper velocities ??? bessere Berechnung über mittelwert aus mehreren werten? Vorfilterung?
-    pltData.velocity_FL_mmps = movmean(pltData.velocity_FL_mmps,100);
-    pltData.velocity_FR_mmps = movmean(pltData.velocity_FR_mmps,100);
-    pltData.velocity_RL_mmps = movmean(pltData.velocity_RL_mmps,100);
-    pltData.velocity_RR_mmps = movmean(pltData.velocity_RR_mmps,100);
-
-    %% IMU
-    % Acceleration
-    % pltData.acc_long_g = movmean(data.XSens_Acceleration_XSens_accX(start:stop),100)/9.81;
-    % pltData.acc_lat_g = movmean(data.XSens_Acceleration_XSens_accY(start:stop),100)/9.81;
-    % Rate of turn
-    % pltData.rot_roll_degps = movmean(data.XSens_RateOfTurn_XSens_gyrX(start:stop),100);
-    % pltData.rot_pitch_degps = movmean(data.XSens_RateOfTurn_XSens_gyrY(start:stop),100);
-    % pltData.rot_yaw_degps = movmean(data.XSens_RateOfTurn_XSens_gyrZ(start:stop),100);
-
-    %% Inverters
-    % inverter temperatures
-    pltData.invL_temp = data.INV_L_TxPDO_1_T_Inv_L(start:stop);
-    pltData.invR_temp = data.INV_R_TxPDO_1_T_Inv_R(start:stop);
-    % motor temperatures
-    pltData.motL_temp = data.INV_L_TxPDO_1_T_Mot_L(start:stop);
-    pltData.motR_temp = data.INV_R_TxPDO_1_T_Mot_R(start:stop);
-    % motor velocities
-    pltData.motL_vel_rpm = 60*data.INV_L_TxPDO_4_Velocity_L(start:stop);
-    pltData.motR_vel_rpm = 60*data.INV_R_TxPDO_4_Velocity_R(start:stop);
-    % inverter torque demand
-    pltData.invL_torqueDemand = data.INV_L_TxPDO_3_DemandedTorque_L(start:stop)/10; % /10 to match autobox torque
-    pltData.invR_torqueDemand = data.INV_R_TxPDO_3_DemandedTorque_R(start:stop)/10;
-    % inverter actual torque
-    pltData.invL_torqueActual = data.INV_L_TxPDO_3_ActualTorque_L(start:stop)/10;
-    pltData.invR_torqueActual = data.INV_R_TxPDO_3_ActualTorque_R(start:stop)/10;
-
-    %% Wheelspeed
-
-
-    
-
-    %% Statistics:
-    power_regen_kw = pltData.ams_ptot;
-    power_regen_kw(power_regen_kw > 0) = 0;
-    power_used_kw = pltData.ams_ptot;
-    power_used_kw(power_used_kw < 0) = 0;
-    pltData.energy_kwh = trapz(pltData.time_s, pltData.ams_ptot)/3600;
-    pltData.energy_regen_kwh = trapz(pltData.time_s, power_regen_kw)/3600;
-    pltData.energy_used_kwh = trapz(pltData.time_s, power_used_kw)/3600;
-    pltData.distanceTotal_km = trapz(pltData.time_s, pltData.speed_kph./3.6)/1000; % 1/3.6 for km/h to m/s and /1000 for km output
-    pltData.peakPower_kw = max(pltData.ams_ptot);
-    pltData.peakPowerMean_kw = max(movmean(pltData.ams_ptot,500));
-
-    pltData.maxSpeed_kph = max(pltData.speed_kph);
-    pltData.startTime = time_hms(start);
-    pltData.stopTime = time_hms(stop);
-end
-
diff --git a/cut_data.m b/cut_data.m
deleted file mode 100644
index b52d186..0000000
--- a/cut_data.m
+++ /dev/null
@@ -1,39 +0,0 @@
-close all
-clear all
-
-% user input:
-input_path = "data\TriggerF.mat";
-output_path = "data\TriggerF_cut.mat";
-startTime = [10, 15, 00];
-stopTime = [17, 34, 00];
-
-% load data
-data_in = load(input_path);
-% get time
-time = days(data_in.ABX_Driver_ABX_Speed(:,1)) + datetime(0,1,0);
-% plot speed for visualization
-plot(time,data_in.ABX_Driver_ABX_Speed(:,2))
-grid on
-% set start/stop timestamps
-date = dateshift(time(1),'start','day');
-date_d = datenum(date)+1;
-startTime_d = date_d + convertHMStoDays(startTime);
-stopTime_d = date_d + convertHMStoDays(stopTime);
-
-fn = fieldnames(data_in);
-for i=1:numel(fn)
-    clear var_out
-    var = data_in.(fn{i});
-    index1 = find(var(:,1)>startTime_d,1);
-    index2 = find(var(:,1)>stopTime_d,1);
-    var_out(:,1) = var(index1:index2,1);
-    var_out(:,2) = var(index1:index2,2);
-    data_out.(fn{i}) = var_out;
-end
-
-save(output_path,"-struct","data_out")
-
-%% Functions
-function time_d = convertHMStoDays(hmsArray)
-    time_d = hmsArray(1)/24 + hmsArray(2)/(24*60) + hmsArray(3)/(24*60*60);
-end
\ No newline at end of file
diff --git a/lap_analysis.m b/lap_analysis.m
deleted file mode 100644
index d24760b..0000000
--- a/lap_analysis.m
+++ /dev/null
@@ -1,229 +0,0 @@
-% Lap analysis script based on vector CAN Data
-% ToDo:
-% - add Driver - Steering -> Oversteer calculation (steering/acc_lat)
-% - fix Driver - Statistics -> Braking distribution to start at >1 bar
-% - calculation damper velocities!!
-% - add size selection for output graphics
-% - add legend for multiple laps selection (instead of statistics?)
-% - automatic plot_settings detection??
-
-close all
-clear all
-format shortEng
-
-%% Add required paths
-addpath("assets\");
-addpath("plot_settings\");
-save_path = "plots\";
-
-% intital variables:
-global pltData;
-pltData = 0;
-global filepath_last;
-filepath_last = "none";
-
-%% create uifigure
-fig = uifigure;
-fig.Position = [100 100 800 600];
-fig.Name = "FT23-Ultra Data Analysis Tool";
-fig.Icon = "assets/FT_icon.png";
-
-% create 5x3 grid layout (row,column)
-g = uigridlayout(fig);
-g.RowHeight = {22, 22,'1x',22,22};
-g.ColumnWidth = {150,'1x',200};
-
-% (1,2) Load path
-filepathField = uieditfield(g);
-filepathField.Value = "data/";
-filepathField.Layout.Row = 1;
-filepathField.Layout.Column = 2;
-% (2,1) label "Plot type:"
-labelPlot = uilabel(g);
-labelPlot.Text = "Plot type:";
-labelPlot.Layout.Row = 2;
-labelPlot.Layout.Column = 1;
-labelPlot.HorizontalAlignment = "right";
-% (2,2) create dropdown selection for plot type
-dd = uidropdown(g);
-dd.Items = ["Driver - General", "Driver - Steering", "Driver - Braking", "Driver - Statistics" ...
-    "Powertrain", "Brakes", "Accumulator", "Inverter", ...
-    "Suspension - Positions", "Suspension - Velocities", "Suspension - Histogram", ...
-    "Tires - Temperatures", "Tires - Friction Circle"];
-dd.Layout.Row = 2;
-dd.Layout.Column = 2;
-% ([3 4],1) create list box with laptimes and multiselection
-lb = uilistbox(g);
-lb.Multiselect = "on";
-lb.Items = "No data loaded!";
-lb.ItemsData = 1;
-lb.Layout.Row = [3 4];
-lb.Layout.Column = 1;
-% (3,2) create plot window
-pltpanel = uipanel(g);
-pltpanel.Layout.Row = 3;
-pltpanel.Layout.Column = 2;
-% (3,3) lap/stint statistics (peak power, max speed, engergy used)
-labelStats = uilabel(g);
-labelStats.Text = "Statistics:";
-labelStats.Layout.Row = 3;
-labelStats.Layout.Column = 3;
-labelStats.VerticalAlignment = "top";
-% (4,2) create dropdown selection for x-Axis (time/distance)
-ddAxis = uidropdown(g);
-ddAxis.Items = ["Time [s]", "Distance [m]", "Time [h:m:s]"];
-ddAxis.ItemsData = linspace(1,3,3);
-ddAxis.Layout.Row = 4;
-ddAxis.Layout.Column = 2;
-% (4,3) ???
-
-% (5,1) dropdown selection for plot output format
-ddFormat = uidropdown(g);
-ddFormat.Items = ["png", "pdf", "fig", "m"];
-ddFormat.ItemsData = ddFormat.Items;
-ddFormat.Layout.Row = 5;
-ddFormat.Layout.Column = 1;
-% (5,2) edit field for output plot name
-editPlotName = uieditfield(g);
-editPlotName.Value = "plot";
-editPlotName.Layout.Row = 5;
-editPlotName.Layout.Column = 2;
-
-% (1,1) browse button
-btnbrowse = uibutton(g,"ButtonPushedFcn", @(src,event) browseButtonPushed(filepathField));
-btnbrowse.Text = "Browse";
-btnbrowse.Layout.Row = 1;
-btnbrowse.Layout.Column = 1;
-% (1,3) Load button
-btnLoad = uibutton(g,"ButtonPushedFcn", @(src,event) loadButtonPushed(lb, filepathField.Value));
-btnLoad.Text = "Load";
-btnLoad.Layout.Row = 1;
-btnLoad.Layout.Column = 3;
-% (2,3) create button for plotting new selection
-btnPlot = uibutton(g,"ButtonPushedFcn", @(src,event) plotButtonPushed(pltpanel, labelStats, lb.Value, dd.Value, ddAxis.Value));
-btnPlot.Text = "Plot";
-btnPlot.Layout.Row = 2;
-btnPlot.Layout.Column = 3;
-% (5,3) create button for saving displayed plot
-btnSave = uibutton(g,"ButtonPushedFcn", @(src,event) saveButtonPushed(pltpanel, append(save_path, editPlotName.Value), ddFormat.Value));
-btnSave.Text = "Save";
-btnSave.Layout.Row = 5;
-btnSave.Layout.Column = 3;
-
-% Plot window FT Logo for style points (and placeholder)
-tl = tiledlayout(pltpanel,"vertical");
-pltpanelDefault = nexttile(tl);
-[imageData,~,imageAlpha] = imread("assets\FT_logo.png");
-image(imageData, "Parent",pltpanelDefault, "alphadata", imageAlpha);
-pltpanelDefault.Visible = "off";
-pltpanelDefault.DataAspectRatio = [1,1,1];
-pltpanelDefault.Toolbar.Visible = "off";
-
-
-%% Functions
-
-function plotButtonPushed(panel, labelStats, selected_laps, plottype, xplot)
-
-    global pltData;
-
-    % If complete stint selected -> disable plotting multiple laps
-    if ismember(1, selected_laps)
-        selected_laps = 1;
-        % xplot = 3;    % select hh:mm:ss format automatically??
-    end
-
-    if size(selected_laps) == 1
-        % update statistics
-        sidebarLabel = sidebar_stats(pltData(selected_laps));
-    else
-        % clear stats
-        sidebarLabel = "not available in multiselection (yet)";
-    end
-    labelStats.Text = sidebarLabel;
-
-    % change x-Axis based on drop-down selection:
-    switch xplot
-        case 2
-            for i = 1:length(pltData)
-                pltData(i).xAxis = pltData(i).distance;
-            end
-        case 3
-            for i = 1:length(pltData)
-                pltData(i).xAxis = pltData(i).time_hms;
-            end
-        otherwise
-            for i = 1:length(pltData)
-                pltData(i).xAxis = pltData(i).time_s;
-            end
-    end
-
-    switch plottype
-        case "Accumulator"
-            plot_accumulator(panel, selected_laps, pltData);
-        case "Brakes"
-            plot_brakes(panel, selected_laps, pltData);
-        case "Driver - Braking"
-            plot_driver_braking(panel, selected_laps, pltData);
-        case "Driver - General"
-            plot_driver_general(panel, selected_laps, pltData);
-        case "Driver - Statistics"
-            plot_driver_statistics(panel, selected_laps, pltData);
-        case "Driver - Steering"
-            plot_driver_steering(panel, selected_laps, pltData);
-        case "Inverter"
-            plot_inverter(panel, selected_laps, pltData);
-        case "Powertrain"
-            plot_powertrain(panel, selected_laps, pltData);
-        case "Suspension - Histogram"
-            plot_suspension_histogram(panel, selected_laps, pltData);
-        case "Suspension - Positions"
-            plot_suspension_positions(panel, selected_laps, pltData);
-        case "Suspension - Velocities"
-            plot_suspension_velocities(panel, selected_laps, pltData);
-        case "Tires - Friction Circle"
-            plot_tires_firctionCircle(panel, selected_laps, pltData);
-        case "Tires - Temperatures"
-            % plot_tires_temperatures(panel, selected_laps, pltData);
-        otherwise
-            error("Error plotting")
-    end
-
-    % legend(ax1,string(selected_laps),'Location','northeastoutside'); 
-end
-
-function loadButtonPushed(listbox, filepath)
-    global pltData;
-    global filepath_last;
-    
-    % check if data is already loaded
-    if strcmp(filepath,filepath_last) == false
-        [pltData,laps_strArray] = loadvector(filepath);
-        % update UI elements
-        listbox.Items = laps_strArray;
-        listbox.ItemsData = linspace(1,length(laps_strArray),length(laps_strArray));
-    else
-        msgbox("Data already loaded!");
-    end
-
-    filepath_last = filepath;
-end
-
-function browseButtonPushed(filepathField)
-    [name, location] = uigetfile(...
-            {'*.m;*.mat',...
-            'MATLAB Files (*.m,*.mat)';
-            '*.*',  'All Files (*.*)'}, ...
-        "Select File","data\");
-    filepath = append(location,name);
-    filepathField.Value = filepath;
-end
-
-function saveButtonPushed(panel, savepath, format)
-    try
-        figure = panel.Children;
-        savePlot(figure,savepath,format);
-        fprintf("Plot saved to: %s.%s\n",savepath,format);
-    catch
-        msgbox("No plot available!")
-    end
-end
\ No newline at end of file
diff --git a/lap_analysis_guide.pdf b/lap_analysis_guide.pdf
deleted file mode 100644
index 5a759e1..0000000
Binary files a/lap_analysis_guide.pdf and /dev/null differ
diff --git a/plot_settings/plot_accumulator.m b/plot_settings/plot_accumulator.m
deleted file mode 100644
index 3fc63e8..0000000
--- a/plot_settings/plot_accumulator.m
+++ /dev/null
@@ -1,37 +0,0 @@
-function [outputArg] = plot_accumulator(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    
-    % plot 1: speed [km/h]
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Speed [km/h]")
-    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
-
-    % plot 2: power [kW]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "Power [kW]")
-    plot(ax2,pltData.xAxis(start:stop),pltData.ams_ptot(start:stop))
-
-    % plot 3: Max Cell Temp [°C]
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Max Cell Temp [°C]")
-    plot(ax3,pltData.xAxis(start:stop),pltData.ams_tmax(start:stop))
-    % plot 4: State of charge [%]
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "SOC")
-    plot(ax4,pltData.xAxis(start:stop),pltData.ams_soc(start:stop))
-    % link all x axes
-    linkaxes([ax1, ax2, ax3, ax4],"x")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_brakes.m b/plot_settings/plot_brakes.m
deleted file mode 100644
index a5f8587..0000000
--- a/plot_settings/plot_brakes.m
+++ /dev/null
@@ -1,53 +0,0 @@
-function [outputArg] = plot_brakes(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    
-    % plot 1: speed
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Speed [km/h]")
-    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
-
-    % plot 2: brake pressure front/rear [bar]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "Brake Pressure F/R [bar]")
-    plot(ax2,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
-    plot(ax2,pltData.xAxis(start:stop),pltData.brakePRear_bar(start:stop))
-    legend(ax2, "Front", "Rear")
-
-    % plot 3: longitudinal acceleration [g]
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Long Acc [g]")
-    plot(ax3,pltData.xAxis(start:stop),pltData.acc_long_g(start:stop))
-
-    % plot 4: Brake Temp [°C]
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "Brake Temp [°C]")
-    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTFrontLeft_degC(start:stop))
-    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTFrontRight_degC(start:stop))
-    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTRearLeft_degC(start:stop))
-    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTRearRight_degC(start:stop))
-    legend(ax4, "FL", "FR", "RL", "RR")
-    
-    % plot 5: brake bias [%]
-    ax5 = nexttile(tl);
-    hold(ax5, "on")
-    grid(ax5, "on")
-    title(ax5, "Brake Bias [%]")
-    plot(ax5,pltData.xAxis(start:stop),pltData.brakeBias_perc(start:stop))
-
-    % link all x axes
-    linkaxes([ax1, ax2, ax3, ax4, ax5],"x")
-
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_driver_braking.m b/plot_settings/plot_driver_braking.m
deleted file mode 100644
index 0648479..0000000
--- a/plot_settings/plot_driver_braking.m
+++ /dev/null
@@ -1,35 +0,0 @@
-function [outputArg] = plot_driver_braking(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    
-    % plot 1: speed
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Speed [km/h]")
-    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
-    % plot 2: accelerator pedal position [%]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "APP [%]")
-    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
-    % plot 3: brake pressure front [bar]
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Brake Pressure Front [bar]")
-    plot(ax3,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
-    % plot 4: longitudinal acceleration [g]
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "Long Acc [g]")
-    plot(ax4,pltData.xAxis(start:stop),pltData.acc_long_g(start:stop))
-    % link all x axes
-    linkaxes([ax1, ax2, ax3, ax4],"x")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_driver_general.m b/plot_settings/plot_driver_general.m
deleted file mode 100644
index e33b68c..0000000
--- a/plot_settings/plot_driver_general.m
+++ /dev/null
@@ -1,35 +0,0 @@
-function [outputArg] = plot_driver_general(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    
-    % plot 1: speed
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Speed [km/h]")
-    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
-    % plot 2: accelerator pedal position [%]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "APP [%]")
-    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
-    % plot 3: brake pressure front [bar]
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Brake Pressure Front [bar]")
-    plot(ax3,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
-    % plot 4: steering angle [°]
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "Steering Angle [°]")
-    plot(ax4,pltData.xAxis(start:stop),pltData.steering_deg(start:stop))
-    % link all x axes
-    linkaxes([ax1, ax2, ax3, ax4],"x")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_driver_statistics.m b/plot_settings/plot_driver_statistics.m
deleted file mode 100644
index a72408e..0000000
--- a/plot_settings/plot_driver_statistics.m
+++ /dev/null
@@ -1,38 +0,0 @@
-function [outputArg] = plot_driver_statistics(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-
-    % plot 1: speed
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "speed [km/h]")
-    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
-    % plot 2: accelerator pedal position [%]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "APP [%]")
-    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
-    linkaxes([ax1, ax2],"x")
-    % plot 3: accelerator pedal position histogram
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "APP distribution")
-    ylabel(ax3, "\sigma [%]")
-    xlabel(ax3, "APP [%]")
-    histogram(ax3,pltData.app_percent(start:stop),"Normalization","percentage")
-    % plot 4: brake pressure histogram
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "BrakeP distribution")
-    ylabel(ax4, "\sigma [%]")
-    xlabel(ax4, "Brake Pressure Front [bar]")
-    histogram(ax4,pltData.brakePFront_bar(start:stop),"Normalization","percentage")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_driver_steering.m b/plot_settings/plot_driver_steering.m
deleted file mode 100644
index dc1eb6f..0000000
--- a/plot_settings/plot_driver_steering.m
+++ /dev/null
@@ -1,36 +0,0 @@
-function [outputArg] = plot_driver_steering(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    
-    % plot 1: speed
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Speed [km/h]")
-    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
-    % plot 2: steering angle [°]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "Steering Angle [°]")
-    plot(ax2,pltData.xAxis(start:stop),pltData.steering_deg(start:stop))
-    % plot 3: oversteer
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Oversteer")
-    % TODO!!
-    
-    % plot 4: lateral acceleration [g]
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "Lateral Acc [g]")
-    plot(ax4,pltData.xAxis(start:stop),pltData.acc_lat_g)(start:stop)
-    % link all x axes
-    linkaxes([ax1, ax2, ax3, ax4],"x")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_inverter.m b/plot_settings/plot_inverter.m
deleted file mode 100644
index 6d57e10..0000000
--- a/plot_settings/plot_inverter.m
+++ /dev/null
@@ -1,50 +0,0 @@
-function [outputArg] = plot_inverter(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    
-    % plot 1: speed [km/h]
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Speed [km/h]")
-    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
-    % plot 2: Inverter Temps [°C]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "Inverter Temp [°C]")
-    plot(ax2,pltData.xAxis(start:stop),pltData.invL_temp(start:stop))
-    plot(ax2,pltData.xAxis(start:stop),pltData.invR_temp(start:stop))
-    legend(ax2, "Left", "Right")
-    % plot 3: Torque request / actual inverter left
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Left Torque")
-    plot(ax3,pltData.xAxis(start:stop),pltData.invL_torqueDemand(start:stop))
-    plot(ax3,pltData.xAxis(start:stop),pltData.invL_torqueActual(start:stop))
-    legend(ax3, "Demand", "Actual")
-    % plot 4: Torque request / actual inverter right
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "Right Torque")
-    plot(ax4,pltData.xAxis(start:stop),pltData.invR_torqueDemand(start:stop))
-    plot(ax4,pltData.xAxis(start:stop),pltData.invR_torqueActual(start:stop))
-    legend(ax4, "Demand", "Actual")
-    % plot 5: Motor velocities
-    ax5 = nexttile(tl);
-    hold(ax5, "on")
-    grid(ax5, "on")
-    title(ax5, "Motor Velocities [1/min]")
-    plot(ax5,pltData.xAxis(start:stop),pltData.motL_vel_rpm(start:stop))
-    plot(ax5,pltData.xAxis(start:stop),pltData.motR_vel_rpm(start:stop))
-    legend(ax5, "Left", "Right")
-
-    % link all x axes
-    linkaxes([ax1, ax2, ax3, ax4, ax5],"x")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_powertrain.m b/plot_settings/plot_powertrain.m
deleted file mode 100644
index d5199f6..0000000
--- a/plot_settings/plot_powertrain.m
+++ /dev/null
@@ -1,48 +0,0 @@
-function [outputArg] = plot_powertrain(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    
-    % plot 1: speed [km/h]
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Speed [km/h]")
-    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
-    % plot 2: power [kW]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "Power [kW]")
-    plot(ax2,pltData.xAxis(start:stop),pltData.ams_ptot(start:stop))
-    % plot 3: Inverter Temps [°C]
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Inverter Temp [°C]")
-    plot(ax3,pltData.xAxis(start:stop),pltData.invL_temp(start:stop))
-    plot(ax3,pltData.xAxis(start:stop),pltData.invR_temp(start:stop))
-    legend(ax3, "Left", "Right")
-    % plot 4: Motor Temps
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "Motor Temp [°C]")
-    plot(ax4,pltData.xAxis(start:stop),pltData.motL_temp(start:stop))
-    plot(ax4,pltData.xAxis(start:stop),pltData.motR_temp(start:stop))
-    legend(ax4, "Left", "Right")
-    % plot 5: Motor velocities
-    ax5 = nexttile(tl);
-    hold(ax5, "on")
-    grid(ax5, "on")
-    title(ax5, "Motor Velocities [1/min]")
-    plot(ax5,pltData.xAxis(start:stop),pltData.motL_vel_rpm(start:stop))
-    plot(ax5,pltData.xAxis(start:stop),pltData.motR_vel_rpm(start:stop))
-    legend(ax5, "Left", "Right")
-
-    % link all x axes
-    linkaxes([ax1, ax2, ax3, ax4, ax5],"x")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_suspension_histogram.m b/plot_settings/plot_suspension_histogram.m
deleted file mode 100644
index ed6363f..0000000
--- a/plot_settings/plot_suspension_histogram.m
+++ /dev/null
@@ -1,33 +0,0 @@
-function [outputArg] = plot_suspension_histogram(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    % plot 1: position heave front [mm/s]
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Heave Front [mm/s]")
-    histogram(ax1, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
-    % plot 2: position roll front [mm/s]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "Roll Front [mm/s]")
-    histogram(ax2, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
-    % plot 3: position heave front [mm/s]
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Heave Rear [mm/s]")
-    histogram(ax3, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
-    % plot 4: position roll front [mm/s]
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "Roll Rear [mm/s]")
-    histogram(ax4, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
-    linkaxes([ax1, ax2, ax3, ax4],"x")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_suspension_positions.m b/plot_settings/plot_suspension_positions.m
deleted file mode 100644
index 49c1ea6..0000000
--- a/plot_settings/plot_suspension_positions.m
+++ /dev/null
@@ -1,33 +0,0 @@
-function [outputArg] = plot_suspension_positions(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    % plot 1: position heave front [mm]
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Heave Front [mm]")
-    plot(ax1, pltData.xAxis(start:stop),pltData.damperHeaveFront_mm(start:stop))
-    % plot 2: position roll front [mm]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "Roll Front [mm]")
-    plot(ax2, pltData.xAxis(start:stop),pltData.damperRollFront_mm(start:stop))
-    % plot 3: position heave front [mm]
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Heave Rear [mm]")
-    plot(ax3, pltData.xAxis(start:stop),pltData.damperHeaveRear_mm(start:stop))
-    % plot 4: position roll front [mm]
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "Roll Rear [mm]")
-    plot(ax4, pltData.xAxis(start:stop),pltData.damperRollRear_mm(start:stop))
-    linkaxes([ax1, ax2, ax3, ax4],"x")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_suspension_velocities.m b/plot_settings/plot_suspension_velocities.m
deleted file mode 100644
index a64a10e..0000000
--- a/plot_settings/plot_suspension_velocities.m
+++ /dev/null
@@ -1,33 +0,0 @@
-function [outputArg] = plot_suspension_velocities(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    tl = tiledlayout(panel,"vertical");
-    % plot 1: velocity heave front [mm/s]
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Heave Front [mm/s]")
-    plot(ax1, pltData.xAxis(start:stop), pltData.velocityHeaveFront_mmps(start:stop))
-    % plot 2: velocity roll front [mm/s]
-    ax2 = nexttile(tl);
-    hold(ax2, "on")
-    grid(ax2, "on")
-    title(ax2, "Roll Front [mm/s]")
-    plot(ax2, pltData.xAxis(start:stop), pltData.velocityRollFront_mmps(start:stop))
-    % plot 3: velocity heave front [mm/s]
-    ax3 = nexttile(tl);
-    hold(ax3, "on")
-    grid(ax3, "on")
-    title(ax3, "Heave Rear [mm/s]")
-    plot(ax3, pltData.xAxis(start:stop), pltData.velocityHeaveRear_mmps(start:stop))
-    % plot 4: velocity roll front [mm/s]
-    ax4 = nexttile(tl);
-    hold(ax4, "on")
-    grid(ax4, "on")
-    title(ax4, "Roll Rear [mm/s]")
-    plot(ax4, pltData.xAxis(start:stop), pltData.velocityRollRear_mmps(start:stop))
-    linkaxes([ax1, ax2, ax3, ax4],"x")
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/plot_tires_firctionCircle.m b/plot_settings/plot_tires_firctionCircle.m
deleted file mode 100644
index ef2830c..0000000
--- a/plot_settings/plot_tires_firctionCircle.m
+++ /dev/null
@@ -1,20 +0,0 @@
-function [outputArg] = plot_tires_firctionCircle(panel, start, stop, pltData)
-    % create tiledlayout (R2023a and newer)
-    
-    tl = tiledlayout(panel,"flow");
-
-    ax1 = nexttile(tl);
-    hold(ax1, "on")
-    grid(ax1, "on")
-    title(ax1, "Friction Circle")
-    ylabel(ax1, "Longitudinal Acc [g]")
-    xlabel(ax1, "Lateral Acc [g]")
-    colors = colororder(ax1);
-    plot(ax1,pltData.acc_lat_g(start:stop), pltData.acc_long_g(start:stop),"Color",colors(i,:))
-    K = convhull(pltData.acc_lat_g(start:stop), pltData.acc_long_g(start:stop));
-    plot(ax1,pltData.acc_lat_g(K), pltData.acc_long_g(K),"Color",colors(i,:),"LineStyle","--","LineWidth",2)
-
-    % return null (not relevant for plots!)
-    outputArg = [];
-end
-
diff --git a/plot_settings/sidebar_stats.m b/plot_settings/sidebar_stats.m
deleted file mode 100644
index 018b9b8..0000000
--- a/plot_settings/sidebar_stats.m
+++ /dev/null
@@ -1,21 +0,0 @@
-function [sidebarLabel] = sidebar_stats(pltData)
-%SIDEBAR_STATS Summary of this function goes here
-%   Detailed explanation goes here
-    sidebarLabel = sprintf(...
-        "Statistics:\nEnergy\n  Used: %.2f kWh\n" + ...
-        "  Regen: %.2f kWh\n" + ...
-        "  Total: %.2f kWh\n" + ...
-        "\nPower\n  Peak: %.1f kW\n" + ...
-        "  Peak (mean): %1.f kW\n" + ...
-        "\nDistance: %.2f km\n" + ...
-        "\nTopspeed: %.1f km/h\n", ...
-        pltData.energy_used_kwh, ...
-        pltData.energy_regen_kwh, ...
-        pltData.energy_kwh, ...
-        pltData.peakPower_kw, ...
-        pltData.peakPowerMean_kw, ...
-        pltData.distanceTotal_km, ...
-        pltData.maxSpeed_kph ...
-        );
-end
-
diff --git a/season/FT23/plot_settings/plot_accumulator.m b/season/FT23/plot_settings/plot_accumulator.m
index e3a8d9a..3fc63e8 100644
--- a/season/FT23/plot_settings/plot_accumulator.m
+++ b/season/FT23/plot_settings/plot_accumulator.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_accumulator(panel, selected_laps, pltData)
+function [outputArg] = plot_accumulator(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,33 +7,27 @@ function [outputArg] = plot_accumulator(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
+
     % plot 2: power [kW]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Power [kW]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_ptot)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.ams_ptot(start:stop))
+
     % plot 3: Max Cell Temp [°C]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Max Cell Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_tmax)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.ams_tmax(start:stop))
     % plot 4: State of charge [%]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "SOC")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_soc)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.ams_soc(start:stop))
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT23/plot_settings/plot_brakes.m b/season/FT23/plot_settings/plot_brakes.m
index 7f6b922..a5f8587 100644
--- a/season/FT23/plot_settings/plot_brakes.m
+++ b/season/FT23/plot_settings/plot_brakes.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_brakes(panel, selected_laps, pltData)
+function [outputArg] = plot_brakes(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,47 +7,41 @@ function [outputArg] = plot_brakes(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
+
     % plot 2: brake pressure front/rear [bar]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Brake Pressure F/R [bar]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePFront_bar)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePRear_bar)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
+    plot(ax2,pltData.xAxis(start:stop),pltData.brakePRear_bar(start:stop))
     legend(ax2, "Front", "Rear")
+
     % plot 3: longitudinal acceleration [g]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Long Acc [g]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).acc_long_g)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.acc_long_g(start:stop))
+
     % plot 4: Brake Temp [°C]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Brake Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTFrontLeft_degC)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTFrontRight_degC)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTRearLeft_degC)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTRearRight_degC)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTFrontLeft_degC(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTFrontRight_degC(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTRearLeft_degC(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTRearRight_degC(start:stop))
     legend(ax4, "FL", "FR", "RL", "RR")
+    
     % plot 5: brake bias [%]
     ax5 = nexttile(tl);
     hold(ax5, "on")
     grid(ax5, "on")
     title(ax5, "Brake Bias [%]")
-    for i = 1:length(selected_laps)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeBias_perc)
-    end
+    plot(ax5,pltData.xAxis(start:stop),pltData.brakeBias_perc(start:stop))
 
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4, ax5],"x")
diff --git a/plot_settings/plot_custom.m b/season/FT23/plot_settings/plot_custom.m
similarity index 100%
rename from plot_settings/plot_custom.m
rename to season/FT23/plot_settings/plot_custom.m
diff --git a/season/FT23/plot_settings/plot_driver_braking.m b/season/FT23/plot_settings/plot_driver_braking.m
index c5643d5..0648479 100644
--- a/season/FT23/plot_settings/plot_driver_braking.m
+++ b/season/FT23/plot_settings/plot_driver_braking.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_braking(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_braking(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,33 +7,25 @@ function [outputArg] = plot_driver_braking(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: accelerator pedal position [%]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "APP [%]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).app_percent)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
     % plot 3: brake pressure front [bar]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Brake Pressure Front [bar]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePFront_bar)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
     % plot 4: longitudinal acceleration [g]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Long Acc [g]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).acc_long_g)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.acc_long_g(start:stop))
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT23/plot_settings/plot_driver_general.m b/season/FT23/plot_settings/plot_driver_general.m
index fbe881a..e33b68c 100644
--- a/season/FT23/plot_settings/plot_driver_general.m
+++ b/season/FT23/plot_settings/plot_driver_general.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_general(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_general(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,33 +7,25 @@ function [outputArg] = plot_driver_general(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: accelerator pedal position [%]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "APP [%]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).app_percent)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
     % plot 3: brake pressure front [bar]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Brake Pressure Front [bar]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePFront_bar)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
     % plot 4: steering angle [°]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Steering Angle [°]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).steering_deg)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.steering_deg(start:stop))
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT23/plot_settings/plot_driver_statistics.m b/season/FT23/plot_settings/plot_driver_statistics.m
index 04cf635..a72408e 100644
--- a/season/FT23/plot_settings/plot_driver_statistics.m
+++ b/season/FT23/plot_settings/plot_driver_statistics.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_statistics(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
 
@@ -7,17 +7,13 @@ function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis,pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: accelerator pedal position [%]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "APP [%]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis,pltData(selected_laps(i)).app_percent)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
     linkaxes([ax1, ax2],"x")
     % plot 3: accelerator pedal position histogram
     ax3 = nexttile(tl);
@@ -26,9 +22,7 @@ function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
     title(ax3, "APP distribution")
     ylabel(ax3, "\sigma [%]")
     xlabel(ax3, "APP [%]")
-    for i = 1:length(selected_laps)
-        histogram(ax3,pltData(selected_laps(i)).app_percent,"Normalization","percentage")
-    end
+    histogram(ax3,pltData.app_percent(start:stop),"Normalization","percentage")
     % plot 4: brake pressure histogram
     ax4 = nexttile(tl);
     hold(ax4, "on")
@@ -36,9 +30,7 @@ function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
     title(ax4, "BrakeP distribution")
     ylabel(ax4, "\sigma [%]")
     xlabel(ax4, "Brake Pressure Front [bar]")
-    for i = 1:length(selected_laps)
-        histogram(ax4,pltData(selected_laps(i)).brakePFront_bar,"Normalization","percentage")
-    end
+    histogram(ax4,pltData.brakePFront_bar(start:stop),"Normalization","percentage")
 
     % return null (not relevant for plots!)
     outputArg = [];
diff --git a/season/FT23/plot_settings/plot_driver_steering.m b/season/FT23/plot_settings/plot_driver_steering.m
index d8066ea..dc1eb6f 100644
--- a/season/FT23/plot_settings/plot_driver_steering.m
+++ b/season/FT23/plot_settings/plot_driver_steering.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_steering(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_steering(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,17 +7,13 @@ function [outputArg] = plot_driver_steering(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: steering angle [°]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Steering Angle [°]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).steering_deg)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.steering_deg(start:stop))
     % plot 3: oversteer
     ax3 = nexttile(tl);
     hold(ax3, "on")
@@ -30,9 +26,7 @@ function [outputArg] = plot_driver_steering(panel, selected_laps, pltData)
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Lateral Acc [g]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).acc_lat_g)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.acc_lat_g)(start:stop)
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT23/plot_settings/plot_inverter.m b/season/FT23/plot_settings/plot_inverter.m
index 878f72f..6d57e10 100644
--- a/season/FT23/plot_settings/plot_inverter.m
+++ b/season/FT23/plot_settings/plot_inverter.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_inverter(panel, selected_laps, pltData)
+function [outputArg] = plot_inverter(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,48 +7,38 @@ function [outputArg] = plot_inverter(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: Inverter Temps [°C]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Inverter Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_temp)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_temp)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.invL_temp(start:stop))
+    plot(ax2,pltData.xAxis(start:stop),pltData.invR_temp(start:stop))
     legend(ax2, "Left", "Right")
     % plot 3: Torque request / actual inverter left
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Left Torque")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_torqueDemand)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_torqueActual)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.invL_torqueDemand(start:stop))
+    plot(ax3,pltData.xAxis(start:stop),pltData.invL_torqueActual(start:stop))
     legend(ax3, "Demand", "Actual")
     % plot 4: Torque request / actual inverter right
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Right Torque")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_torqueDemand)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_torqueActual)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.invR_torqueDemand(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.invR_torqueActual(start:stop))
     legend(ax4, "Demand", "Actual")
     % plot 5: Motor velocities
     ax5 = nexttile(tl);
     hold(ax5, "on")
     grid(ax5, "on")
     title(ax5, "Motor Velocities [1/min]")
-    for i = 1:length(selected_laps)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motL_vel_rpm)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motR_vel_rpm)
-    end
+    plot(ax5,pltData.xAxis(start:stop),pltData.motL_vel_rpm(start:stop))
+    plot(ax5,pltData.xAxis(start:stop),pltData.motR_vel_rpm(start:stop))
     legend(ax5, "Left", "Right")
 
     % link all x axes
diff --git a/season/FT23/plot_settings/plot_powertrain.m b/season/FT23/plot_settings/plot_powertrain.m
index dfb51d8..d5199f6 100644
--- a/season/FT23/plot_settings/plot_powertrain.m
+++ b/season/FT23/plot_settings/plot_powertrain.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_powertrain(panel, selected_laps, pltData)
+function [outputArg] = plot_powertrain(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,46 +7,36 @@ function [outputArg] = plot_powertrain(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: power [kW]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Power [kW]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_ptot)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.ams_ptot(start:stop))
     % plot 3: Inverter Temps [°C]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Inverter Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_temp)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_temp)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.invL_temp(start:stop))
+    plot(ax3,pltData.xAxis(start:stop),pltData.invR_temp(start:stop))
     legend(ax3, "Left", "Right")
     % plot 4: Motor Temps
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Motor Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motL_temp)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motR_temp)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.motL_temp(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.motR_temp(start:stop))
     legend(ax4, "Left", "Right")
     % plot 5: Motor velocities
     ax5 = nexttile(tl);
     hold(ax5, "on")
     grid(ax5, "on")
     title(ax5, "Motor Velocities [1/min]")
-    for i = 1:length(selected_laps)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motL_vel_rpm)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motR_vel_rpm)
-    end
+    plot(ax5,pltData.xAxis(start:stop),pltData.motL_vel_rpm(start:stop))
+    plot(ax5,pltData.xAxis(start:stop),pltData.motR_vel_rpm(start:stop))
     legend(ax5, "Left", "Right")
 
     % link all x axes
diff --git a/season/FT23/plot_settings/plot_suspension_histogram.m b/season/FT23/plot_settings/plot_suspension_histogram.m
index fa1e6b4..ed6363f 100644
--- a/season/FT23/plot_settings/plot_suspension_histogram.m
+++ b/season/FT23/plot_settings/plot_suspension_histogram.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_suspension_histogram(panel, selected_laps, pltData)
+function [outputArg] = plot_suspension_histogram(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     % plot 1: position heave front [mm/s]
@@ -6,33 +6,25 @@ function [outputArg] = plot_suspension_histogram(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Heave Front [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax1, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax1, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     % plot 2: position roll front [mm/s]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Roll Front [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax2, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax2, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     % plot 3: position heave front [mm/s]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Heave Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax3, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax3, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     % plot 4: position roll front [mm/s]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Roll Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax4, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax4, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
     % return null (not relevant for plots!)
diff --git a/season/FT23/plot_settings/plot_suspension_positions.m b/season/FT23/plot_settings/plot_suspension_positions.m
index 844242f..49c1ea6 100644
--- a/season/FT23/plot_settings/plot_suspension_positions.m
+++ b/season/FT23/plot_settings/plot_suspension_positions.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_suspension_positions(panel, selected_laps, pltData)
+function [outputArg] = plot_suspension_positions(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     % plot 1: position heave front [mm]
@@ -6,33 +6,25 @@ function [outputArg] = plot_suspension_positions(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Heave Front [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax1, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperHeaveFront_mm)
-    end
+    plot(ax1, pltData.xAxis(start:stop),pltData.damperHeaveFront_mm(start:stop))
     % plot 2: position roll front [mm]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Roll Front [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax2, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperRollFront_mm)
-    end
+    plot(ax2, pltData.xAxis(start:stop),pltData.damperRollFront_mm(start:stop))
     % plot 3: position heave front [mm]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Heave Rear [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax3, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperHeaveRear_mm)
-    end
+    plot(ax3, pltData.xAxis(start:stop),pltData.damperHeaveRear_mm(start:stop))
     % plot 4: position roll front [mm]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Roll Rear [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax4, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperRollRear_mm)
-    end
+    plot(ax4, pltData.xAxis(start:stop),pltData.damperRollRear_mm(start:stop))
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
     % return null (not relevant for plots!)
diff --git a/season/FT23/plot_settings/plot_suspension_velocities.m b/season/FT23/plot_settings/plot_suspension_velocities.m
index aed80e5..a64a10e 100644
--- a/season/FT23/plot_settings/plot_suspension_velocities.m
+++ b/season/FT23/plot_settings/plot_suspension_velocities.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_suspension_velocities(panel, selected_laps, pltData)
+function [outputArg] = plot_suspension_velocities(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     % plot 1: velocity heave front [mm/s]
@@ -6,33 +6,25 @@ function [outputArg] = plot_suspension_velocities(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Heave Front [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax1, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityHeaveFront_mmps)
-    end
+    plot(ax1, pltData.xAxis(start:stop), pltData.velocityHeaveFront_mmps(start:stop))
     % plot 2: velocity roll front [mm/s]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Roll Front [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax2, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityRollFront_mmps)
-    end
+    plot(ax2, pltData.xAxis(start:stop), pltData.velocityRollFront_mmps(start:stop))
     % plot 3: velocity heave front [mm/s]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Heave Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax3, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityHeaveRear_mmps)
-    end
+    plot(ax3, pltData.xAxis(start:stop), pltData.velocityHeaveRear_mmps(start:stop))
     % plot 4: velocity roll front [mm/s]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Roll Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax4, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityRollRear_mmps)
-    end
+    plot(ax4, pltData.xAxis(start:stop), pltData.velocityRollRear_mmps(start:stop))
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
     % return null (not relevant for plots!)
diff --git a/season/FT23/plot_settings/plot_tires_firctionCircle.m b/season/FT23/plot_settings/plot_tires_firctionCircle.m
index baa117e..ef2830c 100644
--- a/season/FT23/plot_settings/plot_tires_firctionCircle.m
+++ b/season/FT23/plot_settings/plot_tires_firctionCircle.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_tires_firctionCircle(panel, selected_laps, pltData)
+function [outputArg] = plot_tires_firctionCircle(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     
     tl = tiledlayout(panel,"flow");
@@ -10,13 +10,9 @@ function [outputArg] = plot_tires_firctionCircle(panel, selected_laps, pltData)
     ylabel(ax1, "Longitudinal Acc [g]")
     xlabel(ax1, "Lateral Acc [g]")
     colors = colororder(ax1);
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).acc_lat_g, ...
-            pltData(selected_laps(i)).acc_long_g,"Color",colors(i,:))
-        K = convhull(pltData(selected_laps(i)).acc_lat_g, pltData(selected_laps(i)).acc_long_g);
-        plot(ax1,pltData(selected_laps(i)).acc_lat_g(K), ...
-            pltData(selected_laps(i)).acc_long_g(K),"Color",colors(i,:),"LineStyle","--","LineWidth",2)
-    end
+    plot(ax1,pltData.acc_lat_g(start:stop), pltData.acc_long_g(start:stop),"Color",colors(i,:))
+    K = convhull(pltData.acc_lat_g(start:stop), pltData.acc_long_g(start:stop));
+    plot(ax1,pltData.acc_lat_g(K), pltData.acc_long_g(K),"Color",colors(i,:),"LineStyle","--","LineWidth",2)
 
     % return null (not relevant for plots!)
     outputArg = [];
diff --git a/season/FT23/sortdata.m b/season/FT23/sortdata.m
index 9832753..dedf2af 100644
--- a/season/FT23/sortdata.m
+++ b/season/FT23/sortdata.m
@@ -5,6 +5,7 @@ function [pltData] = sortdata(data,start,stop)
     time_s = time_s - time_s(1);    % set beginning of session to 0
 
     % x-Axis:
+    pltData.Time = data.Time(start:stop);
     pltData.time_s = time_s(start:stop)-time_s(start);
     pltData.time_hms = time_hms(start:stop);
 
diff --git a/assets/sortdata_XSens.m b/season/FT23/sortdata_XSens.m
similarity index 100%
rename from assets/sortdata_XSens.m
rename to season/FT23/sortdata_XSens.m
diff --git a/season/FT24/plot_settings/plot_accumulator.m b/season/FT24/plot_settings/plot_accumulator.m
index e3a8d9a..3fc63e8 100644
--- a/season/FT24/plot_settings/plot_accumulator.m
+++ b/season/FT24/plot_settings/plot_accumulator.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_accumulator(panel, selected_laps, pltData)
+function [outputArg] = plot_accumulator(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,33 +7,27 @@ function [outputArg] = plot_accumulator(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
+
     % plot 2: power [kW]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Power [kW]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_ptot)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.ams_ptot(start:stop))
+
     % plot 3: Max Cell Temp [°C]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Max Cell Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_tmax)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.ams_tmax(start:stop))
     % plot 4: State of charge [%]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "SOC")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_soc)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.ams_soc(start:stop))
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT24/plot_settings/plot_brakes.m b/season/FT24/plot_settings/plot_brakes.m
index 7f6b922..a5f8587 100644
--- a/season/FT24/plot_settings/plot_brakes.m
+++ b/season/FT24/plot_settings/plot_brakes.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_brakes(panel, selected_laps, pltData)
+function [outputArg] = plot_brakes(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,47 +7,41 @@ function [outputArg] = plot_brakes(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
+
     % plot 2: brake pressure front/rear [bar]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Brake Pressure F/R [bar]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePFront_bar)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePRear_bar)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
+    plot(ax2,pltData.xAxis(start:stop),pltData.brakePRear_bar(start:stop))
     legend(ax2, "Front", "Rear")
+
     % plot 3: longitudinal acceleration [g]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Long Acc [g]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).acc_long_g)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.acc_long_g(start:stop))
+
     % plot 4: Brake Temp [°C]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Brake Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTFrontLeft_degC)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTFrontRight_degC)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTRearLeft_degC)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTRearRight_degC)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTFrontLeft_degC(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTFrontRight_degC(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTRearLeft_degC(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTRearRight_degC(start:stop))
     legend(ax4, "FL", "FR", "RL", "RR")
+    
     % plot 5: brake bias [%]
     ax5 = nexttile(tl);
     hold(ax5, "on")
     grid(ax5, "on")
     title(ax5, "Brake Bias [%]")
-    for i = 1:length(selected_laps)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeBias_perc)
-    end
+    plot(ax5,pltData.xAxis(start:stop),pltData.brakeBias_perc(start:stop))
 
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4, ax5],"x")
diff --git a/season/FT24/plot_settings/plot_custom.m b/season/FT24/plot_settings/plot_custom.m
new file mode 100644
index 0000000..ee8707f
--- /dev/null
+++ b/season/FT24/plot_settings/plot_custom.m
@@ -0,0 +1,44 @@
+function [outputArg] = plot_custom(panel, start, stop, pltData)
+    % create tiledlayout (R2023a and newer)
+    
+    tl = tiledlayout(panel,"flow");
+
+    ax1 = nexttile(tl);
+    hold(ax1, "on")
+    grid(ax1, "on")
+    title(ax1, "Title")
+    ylabel(ax1, "y-axis label")
+    xlabel(ax1, "x-axis label")
+    %"XXX = enter name of data you want to plot"
+    plot(ax1,pltData.xAxis(start:stop),pltData.XXX(start:stop)) 
+    
+    % more plots for customizing
+
+    % plot 2: 
+    % ax2 = nexttile(tl);
+    % hold(ax2, "on")
+    % grid(ax2, "on")
+    % title(ax2, "Power [kW]")
+    % plot(ax2,pltData.xAxis(start:stop),pltData. XXXX (start:stop))
+
+    % plot 3:
+    % ax3 = nexttile(tl);
+    % hold(ax3, "on")
+    % grid(ax3, "on")
+    % title(ax3, "Max Cell Temp [°C]")
+    % plot(ax3,pltData.xAxis(start:stop),pltData. XXXX (start:stop))
+
+    % plot 4:
+    % ax4 = nexttile(tl);
+    % hold(ax4, "on")
+    % grid(ax4, "on")
+    % title(ax4, "SOC")
+    % plot(ax4,pltData.xAxis(start:stop),pltData. XXXX (start:stop))
+
+
+    % link all x axes
+    linkaxes([ax1, ax2, ax3, ax4],"x")
+    % return null (not relevant for plots!)
+    outputArg = [];
+end
+
diff --git a/season/FT24/plot_settings/plot_driver_braking.m b/season/FT24/plot_settings/plot_driver_braking.m
index c5643d5..0648479 100644
--- a/season/FT24/plot_settings/plot_driver_braking.m
+++ b/season/FT24/plot_settings/plot_driver_braking.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_braking(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_braking(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,33 +7,25 @@ function [outputArg] = plot_driver_braking(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: accelerator pedal position [%]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "APP [%]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).app_percent)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
     % plot 3: brake pressure front [bar]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Brake Pressure Front [bar]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePFront_bar)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
     % plot 4: longitudinal acceleration [g]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Long Acc [g]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).acc_long_g)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.acc_long_g(start:stop))
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT24/plot_settings/plot_driver_general.m b/season/FT24/plot_settings/plot_driver_general.m
index fbe881a..e33b68c 100644
--- a/season/FT24/plot_settings/plot_driver_general.m
+++ b/season/FT24/plot_settings/plot_driver_general.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_general(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_general(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,33 +7,25 @@ function [outputArg] = plot_driver_general(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: accelerator pedal position [%]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "APP [%]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).app_percent)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
     % plot 3: brake pressure front [bar]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Brake Pressure Front [bar]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePFront_bar)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
     % plot 4: steering angle [°]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Steering Angle [°]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).steering_deg)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.steering_deg(start:stop))
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT24/plot_settings/plot_driver_statistics.m b/season/FT24/plot_settings/plot_driver_statistics.m
index 04cf635..a72408e 100644
--- a/season/FT24/plot_settings/plot_driver_statistics.m
+++ b/season/FT24/plot_settings/plot_driver_statistics.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_statistics(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
 
@@ -7,17 +7,13 @@ function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis,pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: accelerator pedal position [%]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "APP [%]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis,pltData(selected_laps(i)).app_percent)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
     linkaxes([ax1, ax2],"x")
     % plot 3: accelerator pedal position histogram
     ax3 = nexttile(tl);
@@ -26,9 +22,7 @@ function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
     title(ax3, "APP distribution")
     ylabel(ax3, "\sigma [%]")
     xlabel(ax3, "APP [%]")
-    for i = 1:length(selected_laps)
-        histogram(ax3,pltData(selected_laps(i)).app_percent,"Normalization","percentage")
-    end
+    histogram(ax3,pltData.app_percent(start:stop),"Normalization","percentage")
     % plot 4: brake pressure histogram
     ax4 = nexttile(tl);
     hold(ax4, "on")
@@ -36,9 +30,7 @@ function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
     title(ax4, "BrakeP distribution")
     ylabel(ax4, "\sigma [%]")
     xlabel(ax4, "Brake Pressure Front [bar]")
-    for i = 1:length(selected_laps)
-        histogram(ax4,pltData(selected_laps(i)).brakePFront_bar,"Normalization","percentage")
-    end
+    histogram(ax4,pltData.brakePFront_bar(start:stop),"Normalization","percentage")
 
     % return null (not relevant for plots!)
     outputArg = [];
diff --git a/season/FT24/plot_settings/plot_driver_steering.m b/season/FT24/plot_settings/plot_driver_steering.m
index d8066ea..dc1eb6f 100644
--- a/season/FT24/plot_settings/plot_driver_steering.m
+++ b/season/FT24/plot_settings/plot_driver_steering.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_steering(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_steering(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,17 +7,13 @@ function [outputArg] = plot_driver_steering(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: steering angle [°]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Steering Angle [°]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).steering_deg)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.steering_deg(start:stop))
     % plot 3: oversteer
     ax3 = nexttile(tl);
     hold(ax3, "on")
@@ -30,9 +26,7 @@ function [outputArg] = plot_driver_steering(panel, selected_laps, pltData)
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Lateral Acc [g]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).acc_lat_g)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.acc_lat_g)(start:stop)
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT24/plot_settings/plot_inverter.m b/season/FT24/plot_settings/plot_inverter.m
index 878f72f..6d57e10 100644
--- a/season/FT24/plot_settings/plot_inverter.m
+++ b/season/FT24/plot_settings/plot_inverter.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_inverter(panel, selected_laps, pltData)
+function [outputArg] = plot_inverter(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,48 +7,38 @@ function [outputArg] = plot_inverter(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: Inverter Temps [°C]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Inverter Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_temp)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_temp)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.invL_temp(start:stop))
+    plot(ax2,pltData.xAxis(start:stop),pltData.invR_temp(start:stop))
     legend(ax2, "Left", "Right")
     % plot 3: Torque request / actual inverter left
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Left Torque")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_torqueDemand)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_torqueActual)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.invL_torqueDemand(start:stop))
+    plot(ax3,pltData.xAxis(start:stop),pltData.invL_torqueActual(start:stop))
     legend(ax3, "Demand", "Actual")
     % plot 4: Torque request / actual inverter right
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Right Torque")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_torqueDemand)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_torqueActual)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.invR_torqueDemand(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.invR_torqueActual(start:stop))
     legend(ax4, "Demand", "Actual")
     % plot 5: Motor velocities
     ax5 = nexttile(tl);
     hold(ax5, "on")
     grid(ax5, "on")
     title(ax5, "Motor Velocities [1/min]")
-    for i = 1:length(selected_laps)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motL_vel_rpm)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motR_vel_rpm)
-    end
+    plot(ax5,pltData.xAxis(start:stop),pltData.motL_vel_rpm(start:stop))
+    plot(ax5,pltData.xAxis(start:stop),pltData.motR_vel_rpm(start:stop))
     legend(ax5, "Left", "Right")
 
     % link all x axes
diff --git a/season/FT24/plot_settings/plot_powertrain.m b/season/FT24/plot_settings/plot_powertrain.m
index dfb51d8..d5199f6 100644
--- a/season/FT24/plot_settings/plot_powertrain.m
+++ b/season/FT24/plot_settings/plot_powertrain.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_powertrain(panel, selected_laps, pltData)
+function [outputArg] = plot_powertrain(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,46 +7,36 @@ function [outputArg] = plot_powertrain(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: power [kW]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Power [kW]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_ptot)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.ams_ptot(start:stop))
     % plot 3: Inverter Temps [°C]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Inverter Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_temp)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_temp)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.invL_temp(start:stop))
+    plot(ax3,pltData.xAxis(start:stop),pltData.invR_temp(start:stop))
     legend(ax3, "Left", "Right")
     % plot 4: Motor Temps
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Motor Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motL_temp)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motR_temp)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.motL_temp(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.motR_temp(start:stop))
     legend(ax4, "Left", "Right")
     % plot 5: Motor velocities
     ax5 = nexttile(tl);
     hold(ax5, "on")
     grid(ax5, "on")
     title(ax5, "Motor Velocities [1/min]")
-    for i = 1:length(selected_laps)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motL_vel_rpm)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motR_vel_rpm)
-    end
+    plot(ax5,pltData.xAxis(start:stop),pltData.motL_vel_rpm(start:stop))
+    plot(ax5,pltData.xAxis(start:stop),pltData.motR_vel_rpm(start:stop))
     legend(ax5, "Left", "Right")
 
     % link all x axes
diff --git a/season/FT24/plot_settings/plot_suspension_histogram.m b/season/FT24/plot_settings/plot_suspension_histogram.m
index fa1e6b4..ed6363f 100644
--- a/season/FT24/plot_settings/plot_suspension_histogram.m
+++ b/season/FT24/plot_settings/plot_suspension_histogram.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_suspension_histogram(panel, selected_laps, pltData)
+function [outputArg] = plot_suspension_histogram(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     % plot 1: position heave front [mm/s]
@@ -6,33 +6,25 @@ function [outputArg] = plot_suspension_histogram(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Heave Front [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax1, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax1, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     % plot 2: position roll front [mm/s]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Roll Front [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax2, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax2, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     % plot 3: position heave front [mm/s]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Heave Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax3, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax3, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     % plot 4: position roll front [mm/s]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Roll Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax4, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax4, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
     % return null (not relevant for plots!)
diff --git a/season/FT24/plot_settings/plot_suspension_positions.m b/season/FT24/plot_settings/plot_suspension_positions.m
index 844242f..49c1ea6 100644
--- a/season/FT24/plot_settings/plot_suspension_positions.m
+++ b/season/FT24/plot_settings/plot_suspension_positions.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_suspension_positions(panel, selected_laps, pltData)
+function [outputArg] = plot_suspension_positions(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     % plot 1: position heave front [mm]
@@ -6,33 +6,25 @@ function [outputArg] = plot_suspension_positions(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Heave Front [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax1, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperHeaveFront_mm)
-    end
+    plot(ax1, pltData.xAxis(start:stop),pltData.damperHeaveFront_mm(start:stop))
     % plot 2: position roll front [mm]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Roll Front [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax2, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperRollFront_mm)
-    end
+    plot(ax2, pltData.xAxis(start:stop),pltData.damperRollFront_mm(start:stop))
     % plot 3: position heave front [mm]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Heave Rear [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax3, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperHeaveRear_mm)
-    end
+    plot(ax3, pltData.xAxis(start:stop),pltData.damperHeaveRear_mm(start:stop))
     % plot 4: position roll front [mm]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Roll Rear [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax4, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperRollRear_mm)
-    end
+    plot(ax4, pltData.xAxis(start:stop),pltData.damperRollRear_mm(start:stop))
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
     % return null (not relevant for plots!)
diff --git a/season/FT24/plot_settings/plot_suspension_velocities.m b/season/FT24/plot_settings/plot_suspension_velocities.m
index aed80e5..a64a10e 100644
--- a/season/FT24/plot_settings/plot_suspension_velocities.m
+++ b/season/FT24/plot_settings/plot_suspension_velocities.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_suspension_velocities(panel, selected_laps, pltData)
+function [outputArg] = plot_suspension_velocities(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     % plot 1: velocity heave front [mm/s]
@@ -6,33 +6,25 @@ function [outputArg] = plot_suspension_velocities(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Heave Front [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax1, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityHeaveFront_mmps)
-    end
+    plot(ax1, pltData.xAxis(start:stop), pltData.velocityHeaveFront_mmps(start:stop))
     % plot 2: velocity roll front [mm/s]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Roll Front [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax2, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityRollFront_mmps)
-    end
+    plot(ax2, pltData.xAxis(start:stop), pltData.velocityRollFront_mmps(start:stop))
     % plot 3: velocity heave front [mm/s]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Heave Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax3, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityHeaveRear_mmps)
-    end
+    plot(ax3, pltData.xAxis(start:stop), pltData.velocityHeaveRear_mmps(start:stop))
     % plot 4: velocity roll front [mm/s]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Roll Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax4, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityRollRear_mmps)
-    end
+    plot(ax4, pltData.xAxis(start:stop), pltData.velocityRollRear_mmps(start:stop))
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
     % return null (not relevant for plots!)
diff --git a/season/FT24/plot_settings/plot_tires_firctionCircle.m b/season/FT24/plot_settings/plot_tires_firctionCircle.m
index baa117e..ef2830c 100644
--- a/season/FT24/plot_settings/plot_tires_firctionCircle.m
+++ b/season/FT24/plot_settings/plot_tires_firctionCircle.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_tires_firctionCircle(panel, selected_laps, pltData)
+function [outputArg] = plot_tires_firctionCircle(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     
     tl = tiledlayout(panel,"flow");
@@ -10,13 +10,9 @@ function [outputArg] = plot_tires_firctionCircle(panel, selected_laps, pltData)
     ylabel(ax1, "Longitudinal Acc [g]")
     xlabel(ax1, "Lateral Acc [g]")
     colors = colororder(ax1);
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).acc_lat_g, ...
-            pltData(selected_laps(i)).acc_long_g,"Color",colors(i,:))
-        K = convhull(pltData(selected_laps(i)).acc_lat_g, pltData(selected_laps(i)).acc_long_g);
-        plot(ax1,pltData(selected_laps(i)).acc_lat_g(K), ...
-            pltData(selected_laps(i)).acc_long_g(K),"Color",colors(i,:),"LineStyle","--","LineWidth",2)
-    end
+    plot(ax1,pltData.acc_lat_g(start:stop), pltData.acc_long_g(start:stop),"Color",colors(i,:))
+    K = convhull(pltData.acc_lat_g(start:stop), pltData.acc_long_g(start:stop));
+    plot(ax1,pltData.acc_lat_g(K), pltData.acc_long_g(K),"Color",colors(i,:),"LineStyle","--","LineWidth",2)
 
     % return null (not relevant for plots!)
     outputArg = [];
diff --git a/season/FT24/sortdata.m b/season/FT24/sortdata.m
index 9832753..dedf2af 100644
--- a/season/FT24/sortdata.m
+++ b/season/FT24/sortdata.m
@@ -5,6 +5,7 @@ function [pltData] = sortdata(data,start,stop)
     time_s = time_s - time_s(1);    % set beginning of session to 0
 
     % x-Axis:
+    pltData.Time = data.Time(start:stop);
     pltData.time_s = time_s(start:stop)-time_s(start);
     pltData.time_hms = time_hms(start:stop);
 
diff --git a/season/FT25/plot_settings/plot_accumulator.m b/season/FT25/plot_settings/plot_accumulator.m
index e3a8d9a..3fc63e8 100644
--- a/season/FT25/plot_settings/plot_accumulator.m
+++ b/season/FT25/plot_settings/plot_accumulator.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_accumulator(panel, selected_laps, pltData)
+function [outputArg] = plot_accumulator(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,33 +7,27 @@ function [outputArg] = plot_accumulator(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
+
     % plot 2: power [kW]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Power [kW]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_ptot)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.ams_ptot(start:stop))
+
     % plot 3: Max Cell Temp [°C]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Max Cell Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_tmax)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.ams_tmax(start:stop))
     % plot 4: State of charge [%]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "SOC")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_soc)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.ams_soc(start:stop))
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT25/plot_settings/plot_brakes.m b/season/FT25/plot_settings/plot_brakes.m
index 7f6b922..a5f8587 100644
--- a/season/FT25/plot_settings/plot_brakes.m
+++ b/season/FT25/plot_settings/plot_brakes.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_brakes(panel, selected_laps, pltData)
+function [outputArg] = plot_brakes(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,47 +7,41 @@ function [outputArg] = plot_brakes(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
+
     % plot 2: brake pressure front/rear [bar]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Brake Pressure F/R [bar]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePFront_bar)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePRear_bar)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
+    plot(ax2,pltData.xAxis(start:stop),pltData.brakePRear_bar(start:stop))
     legend(ax2, "Front", "Rear")
+
     % plot 3: longitudinal acceleration [g]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Long Acc [g]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).acc_long_g)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.acc_long_g(start:stop))
+
     % plot 4: Brake Temp [°C]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Brake Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTFrontLeft_degC)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTFrontRight_degC)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTRearLeft_degC)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeTRearRight_degC)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTFrontLeft_degC(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTFrontRight_degC(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTRearLeft_degC(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.brakeTRearRight_degC(start:stop))
     legend(ax4, "FL", "FR", "RL", "RR")
+    
     % plot 5: brake bias [%]
     ax5 = nexttile(tl);
     hold(ax5, "on")
     grid(ax5, "on")
     title(ax5, "Brake Bias [%]")
-    for i = 1:length(selected_laps)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakeBias_perc)
-    end
+    plot(ax5,pltData.xAxis(start:stop),pltData.brakeBias_perc(start:stop))
 
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4, ax5],"x")
diff --git a/season/FT25/plot_settings/plot_custom.m b/season/FT25/plot_settings/plot_custom.m
new file mode 100644
index 0000000..ee8707f
--- /dev/null
+++ b/season/FT25/plot_settings/plot_custom.m
@@ -0,0 +1,44 @@
+function [outputArg] = plot_custom(panel, start, stop, pltData)
+    % create tiledlayout (R2023a and newer)
+    
+    tl = tiledlayout(panel,"flow");
+
+    ax1 = nexttile(tl);
+    hold(ax1, "on")
+    grid(ax1, "on")
+    title(ax1, "Title")
+    ylabel(ax1, "y-axis label")
+    xlabel(ax1, "x-axis label")
+    %"XXX = enter name of data you want to plot"
+    plot(ax1,pltData.xAxis(start:stop),pltData.XXX(start:stop)) 
+    
+    % more plots for customizing
+
+    % plot 2: 
+    % ax2 = nexttile(tl);
+    % hold(ax2, "on")
+    % grid(ax2, "on")
+    % title(ax2, "Power [kW]")
+    % plot(ax2,pltData.xAxis(start:stop),pltData. XXXX (start:stop))
+
+    % plot 3:
+    % ax3 = nexttile(tl);
+    % hold(ax3, "on")
+    % grid(ax3, "on")
+    % title(ax3, "Max Cell Temp [°C]")
+    % plot(ax3,pltData.xAxis(start:stop),pltData. XXXX (start:stop))
+
+    % plot 4:
+    % ax4 = nexttile(tl);
+    % hold(ax4, "on")
+    % grid(ax4, "on")
+    % title(ax4, "SOC")
+    % plot(ax4,pltData.xAxis(start:stop),pltData. XXXX (start:stop))
+
+
+    % link all x axes
+    linkaxes([ax1, ax2, ax3, ax4],"x")
+    % return null (not relevant for plots!)
+    outputArg = [];
+end
+
diff --git a/season/FT25/plot_settings/plot_driver_braking.m b/season/FT25/plot_settings/plot_driver_braking.m
index c5643d5..0648479 100644
--- a/season/FT25/plot_settings/plot_driver_braking.m
+++ b/season/FT25/plot_settings/plot_driver_braking.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_braking(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_braking(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,33 +7,25 @@ function [outputArg] = plot_driver_braking(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: accelerator pedal position [%]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "APP [%]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).app_percent)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
     % plot 3: brake pressure front [bar]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Brake Pressure Front [bar]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePFront_bar)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
     % plot 4: longitudinal acceleration [g]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Long Acc [g]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).acc_long_g)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.acc_long_g(start:stop))
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT25/plot_settings/plot_driver_general.m b/season/FT25/plot_settings/plot_driver_general.m
index fbe881a..e33b68c 100644
--- a/season/FT25/plot_settings/plot_driver_general.m
+++ b/season/FT25/plot_settings/plot_driver_general.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_general(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_general(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,33 +7,25 @@ function [outputArg] = plot_driver_general(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: accelerator pedal position [%]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "APP [%]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).app_percent)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
     % plot 3: brake pressure front [bar]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Brake Pressure Front [bar]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).brakePFront_bar)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.brakePFront_bar(start:stop))
     % plot 4: steering angle [°]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Steering Angle [°]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).steering_deg)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.steering_deg(start:stop))
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT25/plot_settings/plot_driver_statistics.m b/season/FT25/plot_settings/plot_driver_statistics.m
index 04cf635..a72408e 100644
--- a/season/FT25/plot_settings/plot_driver_statistics.m
+++ b/season/FT25/plot_settings/plot_driver_statistics.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_statistics(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
 
@@ -7,17 +7,13 @@ function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis,pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: accelerator pedal position [%]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "APP [%]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis,pltData(selected_laps(i)).app_percent)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.app_percent(start:stop))
     linkaxes([ax1, ax2],"x")
     % plot 3: accelerator pedal position histogram
     ax3 = nexttile(tl);
@@ -26,9 +22,7 @@ function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
     title(ax3, "APP distribution")
     ylabel(ax3, "\sigma [%]")
     xlabel(ax3, "APP [%]")
-    for i = 1:length(selected_laps)
-        histogram(ax3,pltData(selected_laps(i)).app_percent,"Normalization","percentage")
-    end
+    histogram(ax3,pltData.app_percent(start:stop),"Normalization","percentage")
     % plot 4: brake pressure histogram
     ax4 = nexttile(tl);
     hold(ax4, "on")
@@ -36,9 +30,7 @@ function [outputArg] = plot_driver_statistics(panel, selected_laps, pltData)
     title(ax4, "BrakeP distribution")
     ylabel(ax4, "\sigma [%]")
     xlabel(ax4, "Brake Pressure Front [bar]")
-    for i = 1:length(selected_laps)
-        histogram(ax4,pltData(selected_laps(i)).brakePFront_bar,"Normalization","percentage")
-    end
+    histogram(ax4,pltData.brakePFront_bar(start:stop),"Normalization","percentage")
 
     % return null (not relevant for plots!)
     outputArg = [];
diff --git a/season/FT25/plot_settings/plot_driver_steering.m b/season/FT25/plot_settings/plot_driver_steering.m
index d8066ea..dc1eb6f 100644
--- a/season/FT25/plot_settings/plot_driver_steering.m
+++ b/season/FT25/plot_settings/plot_driver_steering.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_driver_steering(panel, selected_laps, pltData)
+function [outputArg] = plot_driver_steering(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,17 +7,13 @@ function [outputArg] = plot_driver_steering(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: steering angle [°]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Steering Angle [°]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).steering_deg)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.steering_deg(start:stop))
     % plot 3: oversteer
     ax3 = nexttile(tl);
     hold(ax3, "on")
@@ -30,9 +26,7 @@ function [outputArg] = plot_driver_steering(panel, selected_laps, pltData)
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Lateral Acc [g]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).acc_lat_g)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.acc_lat_g)(start:stop)
     % link all x axes
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
diff --git a/season/FT25/plot_settings/plot_inverter.m b/season/FT25/plot_settings/plot_inverter.m
index 878f72f..6d57e10 100644
--- a/season/FT25/plot_settings/plot_inverter.m
+++ b/season/FT25/plot_settings/plot_inverter.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_inverter(panel, selected_laps, pltData)
+function [outputArg] = plot_inverter(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,48 +7,38 @@ function [outputArg] = plot_inverter(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: Inverter Temps [°C]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Inverter Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_temp)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_temp)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.invL_temp(start:stop))
+    plot(ax2,pltData.xAxis(start:stop),pltData.invR_temp(start:stop))
     legend(ax2, "Left", "Right")
     % plot 3: Torque request / actual inverter left
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Left Torque")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_torqueDemand)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_torqueActual)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.invL_torqueDemand(start:stop))
+    plot(ax3,pltData.xAxis(start:stop),pltData.invL_torqueActual(start:stop))
     legend(ax3, "Demand", "Actual")
     % plot 4: Torque request / actual inverter right
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Right Torque")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_torqueDemand)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_torqueActual)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.invR_torqueDemand(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.invR_torqueActual(start:stop))
     legend(ax4, "Demand", "Actual")
     % plot 5: Motor velocities
     ax5 = nexttile(tl);
     hold(ax5, "on")
     grid(ax5, "on")
     title(ax5, "Motor Velocities [1/min]")
-    for i = 1:length(selected_laps)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motL_vel_rpm)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motR_vel_rpm)
-    end
+    plot(ax5,pltData.xAxis(start:stop),pltData.motL_vel_rpm(start:stop))
+    plot(ax5,pltData.xAxis(start:stop),pltData.motR_vel_rpm(start:stop))
     legend(ax5, "Left", "Right")
 
     % link all x axes
diff --git a/season/FT25/plot_settings/plot_powertrain.m b/season/FT25/plot_settings/plot_powertrain.m
index dfb51d8..d5199f6 100644
--- a/season/FT25/plot_settings/plot_powertrain.m
+++ b/season/FT25/plot_settings/plot_powertrain.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_powertrain(panel, selected_laps, pltData)
+function [outputArg] = plot_powertrain(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     
@@ -7,46 +7,36 @@ function [outputArg] = plot_powertrain(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Speed [km/h]")
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).speed_kph)
-    end
+    plot(ax1,pltData.xAxis(start:stop),pltData.speed_kph(start:stop))
     % plot 2: power [kW]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Power [kW]")
-    for i = 1:length(selected_laps)
-        plot(ax2,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).ams_ptot)
-    end
+    plot(ax2,pltData.xAxis(start:stop),pltData.ams_ptot(start:stop))
     % plot 3: Inverter Temps [°C]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Inverter Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invL_temp)
-        plot(ax3,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).invR_temp)
-    end
+    plot(ax3,pltData.xAxis(start:stop),pltData.invL_temp(start:stop))
+    plot(ax3,pltData.xAxis(start:stop),pltData.invR_temp(start:stop))
     legend(ax3, "Left", "Right")
     % plot 4: Motor Temps
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Motor Temp [°C]")
-    for i = 1:length(selected_laps)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motL_temp)
-        plot(ax4,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motR_temp)
-    end
+    plot(ax4,pltData.xAxis(start:stop),pltData.motL_temp(start:stop))
+    plot(ax4,pltData.xAxis(start:stop),pltData.motR_temp(start:stop))
     legend(ax4, "Left", "Right")
     % plot 5: Motor velocities
     ax5 = nexttile(tl);
     hold(ax5, "on")
     grid(ax5, "on")
     title(ax5, "Motor Velocities [1/min]")
-    for i = 1:length(selected_laps)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motL_vel_rpm)
-        plot(ax5,pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).motR_vel_rpm)
-    end
+    plot(ax5,pltData.xAxis(start:stop),pltData.motL_vel_rpm(start:stop))
+    plot(ax5,pltData.xAxis(start:stop),pltData.motR_vel_rpm(start:stop))
     legend(ax5, "Left", "Right")
 
     % link all x axes
diff --git a/season/FT25/plot_settings/plot_suspension_histogram.m b/season/FT25/plot_settings/plot_suspension_histogram.m
index fa1e6b4..ed6363f 100644
--- a/season/FT25/plot_settings/plot_suspension_histogram.m
+++ b/season/FT25/plot_settings/plot_suspension_histogram.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_suspension_histogram(panel, selected_laps, pltData)
+function [outputArg] = plot_suspension_histogram(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     % plot 1: position heave front [mm/s]
@@ -6,33 +6,25 @@ function [outputArg] = plot_suspension_histogram(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Heave Front [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax1, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax1, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     % plot 2: position roll front [mm/s]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Roll Front [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax2, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax2, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     % plot 3: position heave front [mm/s]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Heave Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax3, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax3, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     % plot 4: position roll front [mm/s]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Roll Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        histogram(ax4, pltData(selected_laps(i)).velocityHeaveFront_mmps,"Normalization","percentage")
-    end
+    histogram(ax4, pltData.velocityHeaveFront_mmps(start:stop),"Normalization","percentage")
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
     % return null (not relevant for plots!)
diff --git a/season/FT25/plot_settings/plot_suspension_positions.m b/season/FT25/plot_settings/plot_suspension_positions.m
index 844242f..49c1ea6 100644
--- a/season/FT25/plot_settings/plot_suspension_positions.m
+++ b/season/FT25/plot_settings/plot_suspension_positions.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_suspension_positions(panel, selected_laps, pltData)
+function [outputArg] = plot_suspension_positions(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     % plot 1: position heave front [mm]
@@ -6,33 +6,25 @@ function [outputArg] = plot_suspension_positions(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Heave Front [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax1, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperHeaveFront_mm)
-    end
+    plot(ax1, pltData.xAxis(start:stop),pltData.damperHeaveFront_mm(start:stop))
     % plot 2: position roll front [mm]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Roll Front [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax2, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperRollFront_mm)
-    end
+    plot(ax2, pltData.xAxis(start:stop),pltData.damperRollFront_mm(start:stop))
     % plot 3: position heave front [mm]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Heave Rear [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax3, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperHeaveRear_mm)
-    end
+    plot(ax3, pltData.xAxis(start:stop),pltData.damperHeaveRear_mm(start:stop))
     % plot 4: position roll front [mm]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Roll Rear [mm]")
-    for i = 1:length(selected_laps)
-        plot(ax4, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).damperRollRear_mm)
-    end
+    plot(ax4, pltData.xAxis(start:stop),pltData.damperRollRear_mm(start:stop))
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
     % return null (not relevant for plots!)
diff --git a/season/FT25/plot_settings/plot_suspension_velocities.m b/season/FT25/plot_settings/plot_suspension_velocities.m
index aed80e5..a64a10e 100644
--- a/season/FT25/plot_settings/plot_suspension_velocities.m
+++ b/season/FT25/plot_settings/plot_suspension_velocities.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_suspension_velocities(panel, selected_laps, pltData)
+function [outputArg] = plot_suspension_velocities(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     tl = tiledlayout(panel,"vertical");
     % plot 1: velocity heave front [mm/s]
@@ -6,33 +6,25 @@ function [outputArg] = plot_suspension_velocities(panel, selected_laps, pltData)
     hold(ax1, "on")
     grid(ax1, "on")
     title(ax1, "Heave Front [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax1, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityHeaveFront_mmps)
-    end
+    plot(ax1, pltData.xAxis(start:stop), pltData.velocityHeaveFront_mmps(start:stop))
     % plot 2: velocity roll front [mm/s]
     ax2 = nexttile(tl);
     hold(ax2, "on")
     grid(ax2, "on")
     title(ax2, "Roll Front [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax2, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityRollFront_mmps)
-    end
+    plot(ax2, pltData.xAxis(start:stop), pltData.velocityRollFront_mmps(start:stop))
     % plot 3: velocity heave front [mm/s]
     ax3 = nexttile(tl);
     hold(ax3, "on")
     grid(ax3, "on")
     title(ax3, "Heave Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax3, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityHeaveRear_mmps)
-    end
+    plot(ax3, pltData.xAxis(start:stop), pltData.velocityHeaveRear_mmps(start:stop))
     % plot 4: velocity roll front [mm/s]
     ax4 = nexttile(tl);
     hold(ax4, "on")
     grid(ax4, "on")
     title(ax4, "Roll Rear [mm/s]")
-    for i = 1:length(selected_laps)
-        plot(ax4, pltData(selected_laps(i)).xAxis, pltData(selected_laps(i)).velocityRollRear_mmps)
-    end
+    plot(ax4, pltData.xAxis(start:stop), pltData.velocityRollRear_mmps(start:stop))
     linkaxes([ax1, ax2, ax3, ax4],"x")
 
     % return null (not relevant for plots!)
diff --git a/season/FT25/plot_settings/plot_tires_firctionCircle.m b/season/FT25/plot_settings/plot_tires_firctionCircle.m
index baa117e..ef2830c 100644
--- a/season/FT25/plot_settings/plot_tires_firctionCircle.m
+++ b/season/FT25/plot_settings/plot_tires_firctionCircle.m
@@ -1,4 +1,4 @@
-function [outputArg] = plot_tires_firctionCircle(panel, selected_laps, pltData)
+function [outputArg] = plot_tires_firctionCircle(panel, start, stop, pltData)
     % create tiledlayout (R2023a and newer)
     
     tl = tiledlayout(panel,"flow");
@@ -10,13 +10,9 @@ function [outputArg] = plot_tires_firctionCircle(panel, selected_laps, pltData)
     ylabel(ax1, "Longitudinal Acc [g]")
     xlabel(ax1, "Lateral Acc [g]")
     colors = colororder(ax1);
-    for i = 1:length(selected_laps)
-        plot(ax1,pltData(selected_laps(i)).acc_lat_g, ...
-            pltData(selected_laps(i)).acc_long_g,"Color",colors(i,:))
-        K = convhull(pltData(selected_laps(i)).acc_lat_g, pltData(selected_laps(i)).acc_long_g);
-        plot(ax1,pltData(selected_laps(i)).acc_lat_g(K), ...
-            pltData(selected_laps(i)).acc_long_g(K),"Color",colors(i,:),"LineStyle","--","LineWidth",2)
-    end
+    plot(ax1,pltData.acc_lat_g(start:stop), pltData.acc_long_g(start:stop),"Color",colors(i,:))
+    K = convhull(pltData.acc_lat_g(start:stop), pltData.acc_long_g(start:stop));
+    plot(ax1,pltData.acc_lat_g(K), pltData.acc_long_g(K),"Color",colors(i,:),"LineStyle","--","LineWidth",2)
 
     % return null (not relevant for plots!)
     outputArg = [];
diff --git a/season/FT25/sortdata.m b/season/FT25/sortdata.m
index 9832753..dedf2af 100644
--- a/season/FT25/sortdata.m
+++ b/season/FT25/sortdata.m
@@ -5,6 +5,7 @@ function [pltData] = sortdata(data,start,stop)
     time_s = time_s - time_s(1);    % set beginning of session to 0
 
     % x-Axis:
+    pltData.Time = data.Time(start:stop);
     pltData.time_s = time_s(start:stop)-time_s(start);
     pltData.time_hms = time_hms(start:stop);
 
diff --git a/set_timestamps.m b/set_timestamps.m
deleted file mode 100644
index e69de29..0000000
diff --git a/tutorial/1.png b/tutorial/1.png
new file mode 100644
index 0000000..73b0c1a
Binary files /dev/null and b/tutorial/1.png differ
diff --git a/tutorial/2.png b/tutorial/2.png
new file mode 100644
index 0000000..d3c3626
Binary files /dev/null and b/tutorial/2.png differ
diff --git a/tutorial/3.png b/tutorial/3.png
new file mode 100644
index 0000000..7a227a5
Binary files /dev/null and b/tutorial/3.png differ