Matlab files can be downloaded from https://github.com/dubajicmilos/LocalOctaTilt. Figure 1: To plot MD data use: just_plot_MD.m (Fig. 1a,b bottom right quadrants) To plot experimental data use: Desy_data_figure_creation.m (Fig. 1a,b top quadrants) For plotting Simulated QEDS (phenomenological model) use Simulation_QEDS.m (Fig. 1a,b bottom left quadrants, Fig. 1c,d) To plot A_cation_orientations use plot_A_cation_orientations.m (Fig. 1g,h) Figure 2: To plot individual nanodomain twins comment out the first line below in and uncomment the 2nd line in local_symmetrized_fun.m: S=1*M+1*Mr+Mr1; % We need to sum all three twin local structure components %S=M; %This line can be used if we want to simulate diffuse pattern if only %1 twin is present Figure 3: Here the data is saved through matlab .fig format. The .fig file has to be loaded to Matlab as a figure, and then the data is extractable through matlab functions shown below. For 1D data: (e.g. if you select MAPbBr3_300K_2.5_2.5_lifetimes.fig or FAPbBr3_300K_2.5_2.5_lifetimes.fig) children = get(gca, 'Children'); X = get(children(1), 'XData'); Y_error_1 = get(children(1), 'XData'); X_data = get(children(2), 'XData'); Y_data = get(children(2), 'XData'); X = get(children(3), 'XData'); Y_error_2 = get(children(3), 'XData'); For 2D data:(e.g. if you select MAPbBr3_300K_2D_2.5_2.5_log.fig or FAPbBr3_300K_2D_2.5_2.5_log.fig) h = findobj(gca, 'Type', 'Surface'); % Get the handle to the surf plot Xdata = h.XData; Ydata = h.YData; Zdata = h.ZData; MA_Elastic_Q_Scans_Taipan.fig (Fig 3.a) FA_Elastic_Q_Scans_Taipan.fig (Fig 3.b) Figure 4: See readme.txt in Figure 4 folder.