## 3 Results  
In all samples, the curves showed continuous growth consistent in their patterns (Figure 2), while the thermoluminescence curves showed different trends based on the characteristics of each sample, outlining distinct patterns for each typology (Figure 3). Figure 3: Comparison of natural TL curves for each typology. The signals refer to the results of the “Pre-bleach with Blue LED” protocol. Figure 2: TL signal growth as a function of applied beta dose (5.5, 11, 16.5 Gy) in (**a**) sample SE2 and (**b**) sample SE13. The signal refers to the results of the “Pre-bleached with Blue LED” protocol. Examination of the thermoluminescence (TL) curves revealed five distinctive trends, two of which are specifically associated with heated or partially vitrified rocks. It is remarkable that these trends show no correlation with the presence or absence of heating indicators or other intrinsic rock characteristics. In contrast, when examining the natural thermoluminescence (TL) curves of samples from identical blocks, such as SV4V and SV4C, or SV13V and SV13P, no noteworthy distinctions are observed, despite the different initial typology. However, in the case of SE9 and its two associated samples, SV1A and SV1B, a significant increase in the natural TL incandescence curve becomes evident after crossing the 200 \({}^{\circ}\)C threshold. Furthermore, as introduced before, the curves of the samples were evaluated using a classical thermoluminescence measurement, without further steps. A comparison of some of the curves of samples using the \"Pre-bleached with Blue LED\" method and a classical thermoluminescence protocol is shown in Figure 4 to emphasize how the method proposed in this work cuts out all of the signal before 100 \({}^{\circ}\)C.  
Samples showing fading were corrected using the integration method described in Section 2.4 above. It is crucial to emphasize that fading is not uniformly present in all aliquots. Fading was observed in specific aliquots of four samples: SE2, SE13, SV13P, and SV5, except for sample SV5, where all aliquots showed fading. To perform the correction process, an ad hoc solution was created by scripting a code in Matlab [The Math Works, Inc., Natick, MA, USA, MATLAB, version 2020a]. An in-depth analysis of signal loss was carried out, lasting at least one month, and involving at least two aliquots of each sample. Figure 5 shows one of the fading trends for reference. Upon completion of the equivalent dose assessment utilizing the 'pre-bleached with blue LED' protocol and subsequent dose rate calculations, the age of the samples was determined (see Table 5). The uncertainties associated with these age calculations were computed following the principles of error propagation, as illustrated in Figure 6 and Table 6. \begin{table}
\begin{tabular}{l c} \hline \hline
**Sample Code** & **Proposed Age** \\ \hline Heated or partially melted granitoid rocks \\ \hline SE2 & 115 BCE \(\pm\) 230 \\ SE13 & 450 BCE \(\pm\) 150 \\ SV1 3P & 835 BCE \(\pm\) 395 \\ SV 4C & 970 BCE \(\pm\) 205 \\ \hline SV6 & 1010 BCE \(\pm\) 15 \\ \hline Heated or partially melted metamorphic rocks \\ \hline SE7 & 67 CE \(\pm\) 40 \\ SE12 & 470 BCE \(\pm\) 60 \\ \hline SV5 & 935 BCE \(\pm\) 400 \\ \hline Vitreous cements \\ \hline SV 1A & 566 CE \(\pm\) 710 \\ SV 1B & 370 CE \(\pm\) 130 \\ SV 3V & 42 BCE \(\pm\) 150 \\ SV 4V & 543 BCE \(\pm\) 20 \\ \hline \hline \end{tabular}
\end{table} 
Table 6: Age estimation for each sample typology. Figure 5: Example of the fading curve following one month of analysis (**left**). Experimental data have been normalized relative to the imparted dose of 8.25 Gy. The initial data point is derived immediately after the first minute, followed by subsequent measure. (**Right**) In the figure on the right the results of the hyperbolic fading correction curve with the correct age value reading on the x-axis.