model[&#39;calibration_veritas&#39;]=veritas_

        # Fitting distributions to lower and veritas thresholds
        dfit = distfit(distr=&#39;lognorm&#39;,verbose=None)
        dfit.fit_transform(lower_)
        df, loc, scale = dfit.model[&#39;params&#39;]
        dist = lognorm(df, loc=loc, scale=scale)
        self.veritas_model[&#39;dist_lower&#39;] = dist        
        self.veritas_model[&#39;LOWER_THRESHOLD&#39;] = dist.ppf(self.threshold_alpha)

        dfitv = distfit(smooth=10, distr=&#39;logn