From fbffb0ec424665b761772fbe511b7574f2817ce9 Mon Sep 17 00:00:00 2001
From: Jack Christopher Hutchinson Rolph <jack.rolph@desy.de>
Date: Wed, 10 May 2023 16:06:26 +0200
Subject: [PATCH] Update file Model.py
---
Model.py | 14 +++++++++++---
1 file changed, 11 insertions(+), 3 deletions(-)
diff --git a/Model.py b/Model.py
index e78d026..7651188 100644
--- a/Model.py
+++ b/Model.py
@@ -641,16 +641,24 @@ 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])
dpndQs = []
--
GitLab