fix double counting of merged theta cells impact on shield capacitance

BrieucF · Oct 25, 2023 · a324147 · a324147
1 parent dc99d4a
commit a324147
Showing 1 changed file with 3 additions and 1 deletion.
diff --git 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)