diff --git a/Model.py b/Model.py
index 6fdc8b28c042b8c791c644b84762d35a6a196e3c..c4b4a12b1613e7717a4d4764dab11a736654dac2 100644
--- a/Model.py
+++ b/Model.py
@@ -588,11 +588,6 @@ def dpDCR1dQ(Q, tau, t_gate, t_0, dx, h_order=1):
 
 
 
-
-
-
-
-
 ########################################################################################################
 ## DRM(Q, Q_0, G, mu, lbda, sigma_0, sigma_1, tau_Ap, p_Ap, DCR, tau, tau_R, t_gate, t_0, k_max, len_pad) 
 ########################################################################################################
@@ -605,7 +600,7 @@ def dpDCR1dQ(Q, tau, t_gate, t_0, dx, h_order=1):
 ##
 ########################################################################################################
 
-def DRM(Q, Q_0, G, mu, lbda, sigma_0, sigma_1, tau_Ap, p_Ap, DCR, tau, tau_R, t_gate, t_0, k_max, k_max_dcr, len_pad):
+def DRM(Q, Q_0, G, mu, lbda, sigma_0, sigma_1, tau_Ap, p_Ap, DCR, tau, tau_R, t_gate, t_0, k_max, k_max_dcr, len_pad, A):
     
 
     t_0 = abs(t_0)
@@ -808,10 +803,11 @@ def DRM(Q, Q_0, G, mu, lbda, sigma_0, sigma_1, tau_Ap, p_Ap, DCR, tau, tau_R, t_
         
     
     pdf_drm = pdf_drm[Q_idx_orig]/G
+    count_drm = A*pdf_drm
     
-    
+
         
-    return pdf_drm
+    return count_drm