Merge pull request #29 from giovannimarchiori/gmarchio-main-20240704

update sampling fractions as obtained from ddsim

geometry/create_noise_file_chargePreAmp_theta.py

@@ -54,9 +54,8 @@ def get_ref_charge(SF, E_dep = 1 * pow(10, 6)): #E_dep en eV, choose 1 MeV
                                                      # Assumption: shaping time is similar or bigger to drift time
 
 # Sampling fraction in each layer
-# numbers updated for v03 model with 1536 modules and 11 layers
-# SFfcc = [0.3864252122990472] * 1 + [0.13597644835735828] * 1 + [0.14520427829645913] * 1 + [0.1510076084632846] * 1 + [0.1552347580991012] * 1 + [0.159694330729184] * 1 + [0.1632954482794191] * 1 + [0.16720711037339814] * 1 + [0.17047749048884808] * 1 + [0.17461698117974286] * 1 + [0.1798984163980135] * 1 + [0.17920355117405806] * 1
-SFfcc = [0.3775596654349802] * 1 + [0.13400227700041234] * 1 + [0.14390509963164044] * 1 + [0.14998482026270935] * 1 + [0.15457673722531148] * 1 + [0.15928098152159675] * 1 + [0.1635367867767212] * 1 + [0.16801070646031507] * 1 + [0.1713409944779989] * 1 + [0.17580195406064622] * 1 + [0.17966699467772812] * 1
+# numbers updated for v03 model with 1536 modules and 11 layers, calculated with ddsim
+SFfcc = [0.3800493723322256] * 1 + [0.13494147915064658] * 1 + [0.142866851721152] * 1 + [0.14839315921940666] * 1 + [0.15298362570665006] * 1 + [0.15709704561942747] * 1 + [0.16063717490147533] * 1 + [0.1641723795419055] * 1 + [0.16845490287689746] * 1 + [0.17111520115997653] * 1 + [0.1730605163148862] * 1
 nLayers = len(SFfcc)
 
 print("##########: ", get_ref_charge(0.16))