Extreme TeV BL Lacs, a unique class of blazars with distinct spectral and temporal features, cannot be explained by single-zone models that rely on single shock acceleration. To address this, we developed a dual acceleration model where a recollimation shock and downstream turbulence energize non-thermal electrons. To assess the model's robustness, we compared it with data from some Extreme TeV BL Lacs. We chose to automate the parameter space exploration, focusing on understanding the distributions and cross-correlations of these parameters. This automation of parameter space exploration was facilitated by a Markov Chain Monte Carlo (MCMC) sampler, specifically emcee. The results showed that the parameter distributions closely align with theoretical predictions. However, the cross-correlation among parameters is complex, leading us to conclude that employing an MCMC sampler for parameter space exploration is crucial in capturing the complexity of time-dependent phenomenological models.