|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
|
import FWCore.ParameterSet.Config as cms
process = cms.Process("ProdTPGParam")
# Calo geometry service model
process.load("Geometry.CaloEventSetup.CaloGeometry_cfi")
process.load("Geometry.CaloEventSetup.EcalTrigTowerConstituents_cfi")
process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
# ecal mapping
process.load("Geometry.EcalMapping.EcalMapping_cfi")
process.eegeom = cms.ESSource("EmptyESSource",
recordName = cms.string('EcalMappingRcd'),
iovIsRunNotTime = cms.bool(True),
firstValid = cms.vuint32(1)
)
# Get hardcoded conditions the same used for standard digitization
process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")
#########################
process.source = cms.Source("EmptySource")
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(1)
)
process.TPGParamProducer = cms.EDAnalyzer("EcalTPGParamBuilder",
#### inputs/ouputs control ####
writeToDB = cms.bool(False),
allowDBEE = cms.bool(False),
DBsid = cms.string('ecalh4db'),
DBuser = cms.string('test09'),
DBpass = cms.string('oratest09'),
DBport = cms.uint32(1521),
writeToFiles = cms.bool(True),
outFile = cms.string('TPG_cosmics.txt'),
#### TPG config tag and version (if not given it will be automatically given ) ####
TPGtag = cms.string('CRAFT'),
TPGversion = cms.uint32(1),
#### TPG calculation parameters ####
#######################################################
# Warning: for Cosmics, APD Gain = 200 (i.e. 50x4) #
# ==> All numbers in GeV below must be divided by 4 #
# for the real energy correspondance #
#######################################################
useTransverseEnergy = cms.bool(False), ## true when TPG computes transverse energy, false for energy
Et_sat_EB = cms.double(35.84), ## Saturation value (in GeV) of the TPG before the compressed-LUT (rem: with 35.84 the TPG_LSB = crystal_LSB)
Et_sat_EE = cms.double(61.44), ## Saturation value (in GeV) of the TPG before the compressed-LUT (rem: with 35.84 the TPG_LSB = crystal_LSB)
sliding = cms.uint32(0), ## Parameter used for the FE data format, should'nt be changed
weight_sampleMax = cms.uint32(3), ## position of the maximum among the 5 samples used by the TPG amplitude filter
forcedPedestalValue = cms.int32(160), ## use this value instead of getting it from DB or MC (-1 means use DB or MC)
forceEtaSlice = cms.bool(False), ## when true, same linearization coeff for all crystals belonging to a given eta slice (tower). Implemented only for EB
LUT_option = cms.string('Identity'), ## compressed LUT option can be: "Identity", "Linear", "EcalResolution"
LUT_threshold_EB = cms.double(0.875), ## All Trigger Primitives <= threshold (in GeV) will be set to 0
LUT_threshold_EE = cms.double(0.875), ## All Trigger Primitives <= threshold (in GeV) will be set to 0
LUT_stochastic_EB = cms.double(0.03), ## Stochastic term of the ECAL-EB ET resolution (used only if LUT_option="EcalResolution")
LUT_noise_EB = cms.double(0.2), ## noise term (GeV) of the ECAL-EB ET resolution (used only if LUT_option="EcalResolution")
LUT_constant_EB = cms.double(0.005), ## constant term of the ECAL-EB ET resolution (used only if LUT_option="EcalResolution")
LUT_stochastic_EE = cms.double(0.03), ## Stochastic term of the ECAL-EE ET resolution (used only if LUT_option="EcalResolution")
LUT_noise_EE = cms.double(0.2), ## noise term (GeV) of the ECAL-EE ET resolution (used only if LUT_option="EcalResolution")
LUT_constant_EE = cms.double(0.005), ## constant term of the ECAL-EE ET resolution (used only if LUT_option="EcalResolution")
TTF_lowThreshold_EB = cms.double(0.385), ## EB Trigger Tower Flag low threshold in GeV
TTF_highThreshold_EB = cms.double(0.385), ## EB Trigger Tower Flag high threshold in GeV
TTF_lowThreshold_EE = cms.double(0.385), ## EE Trigger Tower Flag low threshold in GeV
TTF_highThreshold_EE = cms.double(0.385), ## EE Trigger Tower Flag high threshold in GeV
FG_lowThreshold_EB = cms.double(0.7), ## EB Fine Grain Et low threshold in GeV
FG_highThreshold_EB = cms.double(0.7), ## EB Fine Grain Et high threshold in GeV
FG_lowRatio_EB = cms.double(0.8), ## EB Fine Grain low-ratio
FG_highRatio_EB = cms.double(0.8), ## EB Fine Grain high-ratio
FG_lut_EB = cms.uint32(0), ## EB Fine Grain Look-up table. Put something != 0 if you really know what you do!
FG_Threshold_EE = cms.double(0.0), ## EE Fine threshold in GeV
FG_lut_strip_EE = cms.uint32(0), ## EE Fine Grain strip Look-up table
FG_lut_tower_EE = cms.uint32(0) ## EE Fine Grain tower Look-up table
)
process.p = cms.Path(process.TPGParamProducer)
|