From c6fcc0dfb8af339083b13705d50d8b0c6dcf64ae Mon Sep 17 00:00:00 2001
From: Claire Antel <claire.antel@cern.ch>
Date: Thu, 20 Jan 2022 17:06:48 +0100
Subject: [PATCH] adding CP parameter dictionary for ScintDigiAlgs
configuration
---
.../python/ScintDigiAlgsConfig.py | 19 +++++++++++++++----
1 file changed, 15 insertions(+), 4 deletions(-)
diff --git a/Scintillator/ScintDigiAlgs/python/ScintDigiAlgsConfig.py b/Scintillator/ScintDigiAlgs/python/ScintDigiAlgsConfig.py
index 70442751..f6a923ef 100644
--- a/Scintillator/ScintDigiAlgs/python/ScintDigiAlgsConfig.py
+++ b/Scintillator/ScintDigiAlgs/python/ScintDigiAlgsConfig.py
@@ -1,8 +1,19 @@
+# Copyright (C) 2020-2021 CERN for the benefit of the FASER collaboration
+
from AthenaConfiguration.ComponentAccumulator import ComponentAccumulator
from AthenaConfiguration.ComponentFactory import CompFactory
from OutputStreamAthenaPool.OutputStreamConfig import OutputStreamCfg
+# Crystallball parameter dictionary used in simulated digitized wave reconstruction.
+# Crystalball function Parameters estimated from Deion's slides uploaded at
+# https://indico.cern.ch/event/1099652/contributions/4626975/attachments/2352595/4013927/Faser-Physics-run3933-plots.pdf (20/01/2022)
+# Parameters are per scintillator source, but not per channel.
+dict_CB_param = {}
+dict_CB_param["Trigger"]=dict(CB_alpha=-0.38, CB_n=25, CB_mean=815, CB_sigma=7.7)
+dict_CB_param["Timing"]=dict(CB_alpha=-0.32, CB_n=65, CB_mean=846, CB_sigma=5.3) # copy from Preshower; Timing was not in TestBeam
+dict_CB_param["Veto"]=dict(CB_alpha=-0.32, CB_n=65, CB_mean=846, CB_sigma=5.3) # copy from Preshower; Veto was not in TestBeam
+dict_CB_param["Preshower"]=dict(CB_alpha=-0.32, CB_n=65, CB_mean=846, CB_sigma=5.3)
# One stop shopping for normal FASER data
def ScintWaveformDigitizationCfg(flags):
@@ -30,10 +41,10 @@ def ScintWaveformDigiCfg(flags, name="ScintWaveformDigiAlg", source="", **kwargs
kwargs.setdefault("WaveformContainerKey", source+"Waveforms")
digiAlg = CompFactory.ScintWaveformDigiAlg(name, **kwargs)
- digiAlg.CB_alpha = -0.32275
- digiAlg.CB_n = 65
- digiAlg.CB_mean = 846.55
- digiAlg.CB_sigma = 5.372
+ digiAlg.CB_alpha = dict_CB_param[source]["CB_alpha"]
+ digiAlg.CB_n = dict_CB_param[source]["CB_n"]
+ digiAlg.CB_mean = dict_CB_param[source]["CB_mean"]
+ digiAlg.CB_sigma = dict_CB_param[source]["CB_sigma"]
kwargs.setdefault("WaveformDigitisationTool", tool)
--
GitLab