Skip to content
Snippets Groups Projects

Check in old MC conditions updates

Merged Check in old MC conditions updates
torrence/calypso:torrence-digi into master
Merged Eric Torrence requested to merge torrence/calypso:torrence-digi into master
Compare
and
2 files
+ 156
− 55
Compare changes
  • Side-by-side
  • Inline
Files
2
Waveform/WaveformConditions/WaveCondUtils/scripts/makeDigiDB.py
+ 126
− 46
@@ -57,7 +57,77 @@ if not folder:
# 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)
# Fill values into record
# Updated after the geometry updates November 2023
# This also requires different testbeam and TI12 values
testbeam_iov = (1000 << 32)
# Start with the testbeam
base_mean = 15350.
base_rms = 3.
# Calorimeter default (low gain/filters)
calo_record = cool.Record(spec)
calo_record['base_mean'] = base_mean
calo_record['base_rms'] = base_rms
calo_record['cb_norm'] = 0.69
calo_record['cb_mean'] = 815 # Time in ns
calo_record['cb_sigma'] = 4.14
calo_record['cb_alpha'] = -0.32
calo_record['cb_n'] = 1000.
vetonu_record = cool.Record(spec)
vetonu_record['base_mean'] = base_mean
vetonu_record['base_rms'] = base_rms
vetonu_record['cb_norm'] = 19360
vetonu_record['cb_mean'] = 810 # Time in ns
vetonu_record['cb_sigma'] = 7.68
vetonu_record['cb_alpha'] = -0.38
vetonu_record['cb_n'] = 1000.
veto_record = vetonu_record # not sure if testbeam reco uses vetoNu or veto digi parameters
preshower_record = cool.Record(spec)
preshower_record['base_mean'] = base_mean
preshower_record['base_rms'] = base_rms
preshower_record['cb_norm'] = 500
preshower_record['cb_mean'] = 845 # Time in ns
preshower_record['cb_sigma'] = 5.2
preshower_record['cb_alpha'] = -0.30
preshower_record['cb_n'] = 1000.
# Low gain, filters in default
tag = "WAVE-Digitization-02-LG"
calo_record['cb_norm'] = 0.69
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, calo_record, 0, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, vetonu_record, 1, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, veto_record, 2, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, preshower_record, 4, tag)
# Low gain, filters out
tag = "WAVE-Digitization-02-LG-nofilt"
calo_record['cb_norm'] = 6.9
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, calo_record, 0, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, vetonu_record, 1, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, veto_record, 2, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, preshower_record, 4, tag)
# High gain
tag = "WAVE-Digitization-02-HG"
calo_record['cb_norm'] = 20.7
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, calo_record, 0, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, vetonu_record, 1, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, veto_record, 2, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, preshower_record, 4, tag)
# High gain, filters out
tag = "WAVE-Digitization-02-HG-nofilt"
calo_record['cb_norm'] = 207
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, calo_record, 0, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, vetonu_record, 1, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, veto_record, 2, tag)
folder.storeObject( cool.ValidityKeyMin, testbeam_iov, preshower_record, 4, tag)
# Now fill TI12 values
# Values in ADC counts
base_mean = 15650.
base_rms = 3.
@@ -66,7 +136,7 @@ base_rms = 3.
calo_record = cool.Record(spec)
calo_record['base_mean'] = base_mean
calo_record['base_rms'] = base_rms
calo_record['cb_norm'] = 0.56
calo_record['cb_norm'] = 0.55
calo_record['cb_mean'] = 820 # Time in ns
calo_record['cb_sigma'] = 3.67
calo_record['cb_alpha'] = -0.32
@@ -75,8 +145,8 @@ calo_record['cb_n'] = 1000.
vetonu_record = cool.Record(spec)
vetonu_record['base_mean'] = base_mean
vetonu_record['base_rms'] = base_rms
vetonu_record['cb_norm'] = 7040
vetonu_record['cb_mean'] = 815 # Time in ns
vetonu_record['cb_norm'] = 10420
vetonu_record['cb_mean'] = 800 # Time in ns
vetonu_record['cb_sigma'] = 5.0
vetonu_record['cb_alpha'] = -0.28
vetonu_record['cb_n'] = 1000.
@@ -84,8 +154,8 @@ vetonu_record['cb_n'] = 1000.
veto_record = cool.Record(spec)
veto_record['base_mean'] = base_mean
veto_record['base_rms'] = base_rms
veto_record['cb_norm'] = 6840
veto_record['cb_mean'] = 815 # Time in ns
veto_record['cb_norm'] = 8890
veto_record['cb_mean'] = 825 # Time in ns
veto_record['cb_sigma'] = 3.35
veto_record['cb_alpha'] = -0.32
veto_record['cb_n'] = 9.
@@ -93,8 +163,8 @@ veto_record['cb_n'] = 9.
timing_record = cool.Record(spec)
timing_record['base_mean'] = base_mean
timing_record['base_rms'] = base_rms
timing_record['cb_norm'] = 4240
timing_record['cb_mean'] = 815 # Time in ns
timing_record['cb_norm'] = 6360
timing_record['cb_mean'] = 830 # Time in ns
timing_record['cb_sigma'] = 3.21
timing_record['cb_alpha'] = -0.424
timing_record['cb_n'] = 6.14
@@ -102,8 +172,8 @@ timing_record['cb_n'] = 6.14
preshower_record = cool.Record(spec)
preshower_record['base_mean'] = base_mean
preshower_record['base_rms'] = base_rms
preshower_record['cb_norm'] = 400
preshower_record['cb_mean'] = 846 # Time in ns
preshower_record['cb_norm'] = 530
preshower_record['cb_mean'] = 850 # Time in ns
preshower_record['cb_sigma'] = 4.0
preshower_record['cb_alpha'] = -0.32
preshower_record['cb_n'] = 1000.
@@ -112,50 +182,60 @@ preshower_record['cb_n'] = 1000.
# Values: 0.56, 5.6, 16.8, 168
# Dec. 19, 2022
# Lowered by 10% to 0.51, 5.1, 15.3, 153
# Nov. 21, 2023
# Adjusted to 0.55, 5.5, 16.5, 165
# Low gain, filters in default
calo_record['cb_norm'] = 0.51
tag = "WAVE-Digitization-01-LG"
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, calo_record, 0, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, vetonu_record, 1, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, veto_record, 2, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, timing_record, 3, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, preshower_record, 4, tag)
folder.createTagRelation("WAVE-01", tag)
# Also associate this default with all of the global tags
# We will use a folder override to set other conditions
wave_folder.createTagRelation("OFLCOND-FASER-TB00", "WAVE-01")
wave_folder.createTagRelation("OFLCOND-FASER-01", "WAVE-01")
wave_folder.createTagRelation("OFLCOND-FASER-02", "WAVE-01")
wave_folder.createTagRelation("OFLCOND-FASER-03", "WAVE-01")
calo_record['cb_norm'] = 0.55
tag = "WAVE-Digitization-02-LG"
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, calo_record, 0, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, vetonu_record, 1, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, veto_record, 2, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, timing_record, 3, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, preshower_record, 4, tag)
# Low gain, filters out
tag = "WAVE-Digitization-01-LG-nofilt"
calo_record['cb_norm'] = 5.1
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, calo_record, 0, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, vetonu_record, 1, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, veto_record, 2, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, timing_record, 3, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, preshower_record, 4, tag)
tag = "WAVE-Digitization-02-LG-nofilt"
calo_record['cb_norm'] = 5.5
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, calo_record, 0, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, vetonu_record, 1, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, veto_record, 2, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, timing_record, 3, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, preshower_record, 4, tag)
# High gain
tag = "WAVE-Digitization-01-HG"
calo_record['cb_norm'] = 15.3
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, calo_record, 0, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, vetonu_record, 1, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, veto_record, 2, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, timing_record, 3, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, preshower_record, 4, tag)
tag = "WAVE-Digitization-02-HG"
calo_record['cb_norm'] = 16.5
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, calo_record, 0, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, vetonu_record, 1, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, veto_record, 2, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, timing_record, 3, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, preshower_record, 4, tag)
# High gain, filters out
tag = "WAVE-Digitization-01-HG-nofilt"
calo_record['cb_norm'] = 153
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, calo_record, 0, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, vetonu_record, 1, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, veto_record, 2, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, timing_record, 3, tag)
folder.storeObject( cool.ValidityKeyMin, cool.ValidityKeyMax, preshower_record, 4, tag)
tag = "WAVE-Digitization-02-HG-nofilt"
calo_record['cb_norm'] = 165
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, calo_record, 0, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, vetonu_record, 1, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, veto_record, 2, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, timing_record, 3, tag)
folder.storeObject( testbeam_iov, cool.ValidityKeyMax, preshower_record, 4, tag)
# Make the tag releations
# This is the default folder tag
folder.createTagRelation("WAVE-02", "WAVE-Digitization-02-LG")
# Also associate this default with all of the global tags
# We will use a folder override to set other conditions
# Assign to global tags used with new geometry
wave_folder.createTagRelation("OFLCOND-FASER-TB01", "WAVE-02")
wave_folder.createTagRelation("OFLCOND-FASER-04", "WAVE-02")
# Previous
#wave_folder.createTagRelation("OFLCOND-FASER-TB00", "WAVE-01")
#wave_folder.createTagRelation("OFLCOND-FASER-01", "WAVE-01")
#wave_folder.createTagRelation("OFLCOND-FASER-02", "WAVE-01")
#wave_folder.createTagRelation("OFLCOND-FASER-03", "WAVE-01")
# All done
db.closeDatabase()