From a32414795585608b0a6a6987d3f2cef60b9444a7 Mon Sep 17 00:00:00 2001 From: Giovanni Marchiori Date: Wed, 25 Oct 2023 10:14:01 +0200 Subject: [PATCH] fix double counting of merged theta cells impact on shield capacitance --- geometry/create_capacitance_file_theta.py | 4 +++- 1 file changed, 3 insertions(+), 1 deletion(-) diff --git a/geometry/create_capacitance_file_theta.py b/geometry/create_capacitance_file_theta.py index c1ed740..2888eeb 100644 --- a/geometry/create_capacitance_file_theta.py +++ b/geometry/create_capacitance_file_theta.py @@ -262,7 +262,9 @@ # analytical formula (nmultShield = 2) #capacitanceShield = nmultShield * nMergedModules[i] * cellLength * tracesPerLayer[i] * 1 / inch2mm * 0.67 * (epsilonR + 1.41) / logMicrostrip # from maxwell (nmultShield = 1) - capacitanceShield = nmultShield * nMergedModules[i] * nMergedThetaCells[i] * cellLength * tracesPerLayer[i] * capa_per_mm[i] + # dont multiply by nMergedThetaCells: the shield/pad capa is reasonably independent of the cell size and the fact that there is some merging + # done for theta cells is already taken into account by the tracesPerLayer + capacitanceShield = nmultShield * nMergedModules[i] * cellLength * tracesPerLayer[i] * capa_per_mm[i] if capacitanceShield > capa_shield_max: capa_shield_max = capacitanceShield hCapShield[i].SetBinContent(index+1, capacitanceShield)