Update file Model.py

Model.py

@@ -641,15 +641,23 @@ def DRM(Q, Q_0, G, mu, lbda, sigma_0, sigma_1, tau_Ap, p_Ap, DCR, tau, tau_R, t_
     
     Pns = np.zeros((int(k_max), int(k_max)))
     
+    
+    #############################################
+    # COMMENT
+    #############################################
     # alpha is either one of two things:
     # the total probability of an afterpulse occuring within the gate  
     
     alpha = p_Ap
 
     # Or the probability that an electron-hole pair is produced in the silicon bulk.  
-    # alpha = P_Ap(tau_R, tau_Ap, t_gate, p_Ap)
+    # alpha = P_Ap(tau_R, tau_Ap, t_gate, p_Ap) -> here, p_Ap is as is in the LightSimtastic simulation program.
     # 
-    # To replicate the results in the paper, please keep this variable as it is. 
+    # To replicate the results in the paper, please keep this variable as it is, and compare it to the value of 
+    #  P_Ap(tau_R, tau_Ap, t_gate, p_Ap) for the input values of the simulation. 
+    #############################################
+
+
 
     dp0dQs = np.vstack([norm.pdf(k_pad, loc=_k, scale=sigmak(_k, sigma_0/G, sigma_1/G)) for _k in k_pg])