|
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
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
|
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 before CMSSW_3_1_x
## process.load("CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetriever_cfi")
# or Get DB parameters
# process.load('Configuration/StandardSequences/FrontierConditions_GlobalTag_cff')
# process.GlobalTag.globaltag = "GR09_31X_V2H::All"
process.load("CondCore.DBCommon.CondDBCommon_cfi")
process.CondDBCommon.connect = 'oracle://cms_orcon_prod/CMS_COND_31X_ECAL'
process.CondDBCommon.DBParameters.authenticationPath = '/nfshome0/popcondev/conddb'
process.PoolDBESSource = cms.ESSource("PoolDBESSource",
process.CondDBCommon,
toGet = cms.VPSet(
cms.PSet(
record = cms.string('EcalPedestalsRcd'),
#tag = cms.string('EcalPedestals_v5_online')
tag = cms.string('EcalPedestals_2009runs_hlt') ### obviously diff w.r.t previous
),
cms.PSet(
record = cms.string('EcalADCToGeVConstantRcd'),
#tag = cms.string('EcalADCToGeVConstant_EBg50_EEnoB_new')
tag = cms.string('EcalADCToGeVConstant_2009runs_express') ### the 2 ADCtoGEV in EB and EE are diff w.r.t previous
),
cms.PSet(
record = cms.string('EcalIntercalibConstantsRcd'),
#tag = cms.string('EcalIntercalibConstants_EBg50_EEnoB_new')
tag = cms.string('EcalIntercalibConstants_2009runs_express') ### differs from previous
),
cms.PSet(
record = cms.string('EcalGainRatiosRcd'),
#tag = cms.string('EcalGainRatios_TestPulse_online')
tag = cms.string('EcalGainRatios_TestPulse_express') ### no diff w.r.t previous
),
cms.PSet(
record = cms.string('EcalMappingElectronicsRcd'),
tag = cms.string('EcalMappingElectronics_EEMap')
)
)
)
#########################
process.source = cms.Source("EmptySource",
firstRun = cms.untracked.uint32(100000000) ### need to use latest run to pick-up update values from DB
)
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(1)
)
process.TPGParamProducer = cms.EDFilter("EcalTPGParamBuilder",
#### inputs/ouputs control ####
writeToDB = cms.bool(False),
allowDBEE = cms.bool(True),
DBsid = cms.string('cms_omds_lb'),
## P5 online DB
DBuser = cms.string('cms_ecal_conf'),
DBpass = cms.string('0r4cms_3c4lc0nf'),
## test DB
##DBuser = cms.string('cms_ecal_conf_test'),
##DBpass = cms.string('0r4cms_3c4l'),
DBport = cms.uint32(10121),
TPGWritePed = cms.uint32(1),
TPGWriteLin = cms.uint32(1),
TPGWriteSli = cms.uint32(1),
TPGWriteWei = cms.uint32(1),
TPGWriteLut = cms.uint32(1),
TPGWriteFgr = cms.uint32(1),
TPGWriteBxt = cms.uint32(0),
TPGWriteBtt = cms.uint32(0), #do not change
writeToFiles = cms.bool(True),
outFile = cms.string('TPG_beamv3_FG.txt'), #modif
#### TPG config tag and version (if not given it will be automatically given ) ####
TPGtag = cms.string('BEAMV3FORFG'),
TPGversion = cms.uint32(1),
#### TPG calculation parameters ####
useTransverseEnergy = cms.bool(True), ## true when TPG computes transverse energy, false for energy
Et_sat_EB = cms.double(64.0), ## 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(64.0), ## 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_timeShift = cms.double(0.), ## weights are computed shifting the timing of the shape by this amount in ns: val>0 => shape shifted to the right
weight_sampleMax = cms.uint32(3), ## position of the maximum among the 5 samples used by the TPG amplitude filter
weight_unbias_recovery = cms.bool(True), ## true if weights after int conversion are forced to have sum=0. Pb, in that case it can't have sum f*w = 1
forcedPedestalValue = cms.int32(-3), ## use this value instead of getting it from DB or MC
## -1: means use value from DB or MC.
## -2: ped12 = 0 used to cope with FENIX bug
## -3: used with sFGVB: baseline subtracted is pedestal-offset*sin(theta)/G with G=mult*2^-(shift+2)
pedestal_offset = cms.uint32(300), ## pedestal offset used with option forcedPedestalValue = -3
useInterCalibration = cms.bool(True), ## use or not values from DB. If not, 1 is assumed
SFGVB_Threshold = cms.uint32(15), ## used with option forcedPedestalValue = -3
SFGVB_lut = cms.uint32(0xfffefee8), ## used with option forcedPedestalValue = -3
forceEtaSlice = cms.bool(False), ## when true, same linearization coeff for all crystals belonging to a given eta slice (tower)
LUT_option = cms.string('Linear'), ## compressed LUT option can be: "Identity", "Linear", "EcalResolution"
LUT_threshold_EB = cms.double(0.250), ## All Trigger Primitives <= threshold (in GeV) will be set to 0
LUT_threshold_EE = cms.double(0.250), ## 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(1.0), ## EB Trigger Tower Flag low threshold in GeV
TTF_highThreshold_EB = cms.double(2.0), ## EB Trigger Tower Flag high threshold in GeV
TTF_lowThreshold_EE = cms.double(1.0), ## EE Trigger Tower Flag low threshold in GeV
TTF_highThreshold_EE = cms.double(2.0), ## EE Trigger Tower Flag high threshold in GeV
FG_lowThreshold_EB = cms.double(3.9), ## EB Fine Grain Et low threshold in GeV
FG_highThreshold_EB = cms.double(3.9), ## EB Fine Grain Et high threshold in GeV
FG_lowRatio_EB = cms.double(0.95), ## EB Fine Grain low-ratio
FG_highRatio_EB = cms.double(0.95), ## EB Fine Grain high-ratio
FG_lut_EB = cms.uint32(0x08), ## EB Fine Grain Look-up table. Put something != 0 if you really know what you do!
FG_Threshold_EE = cms.double(18.75), ## EE Fine threshold in GeV
FG_lut_strip_EE = cms.uint32(0xfffefee8), ## 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)
|