forked from cms-sw/cmssw
/
ecalMultiFitUncalibRecHit_cfi.py
82 lines (71 loc) · 3.98 KB
/
ecalMultiFitUncalibRecHit_cfi.py
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
import FWCore.ParameterSet.Config as cms
from RecoLocalCalo.EcalRecProducers.ecalPulseShapeParameters_cff import *
ecalMultiFitUncalibRecHit = cms.EDProducer("EcalUncalibRecHitProducer",
EBdigiCollection = cms.InputTag("ecalDigis","ebDigis"),
EEdigiCollection = cms.InputTag("ecalDigis","eeDigis"),
EBhitCollection = cms.string("EcalUncalibRecHitsEB"),
EEhitCollection = cms.string('EcalUncalibRecHitsEE'),
algo = cms.string("EcalUncalibRecHitWorkerMultiFit"),
algoPSet = cms.PSet(
# for multifit method
EcalPulseShapeParameters = cms.PSet( ecal_pulse_shape_parameters ),
activeBXs = cms.vint32(-5,-4,-3,-2,-1,0,1,2,3,4),
ampErrorCalculation = cms.bool(True),
useLumiInfoRunHeader = cms.bool(True),
doPrefitEB = cms.bool(False),
doPrefitEE = cms.bool(False),
prefitMaxChiSqEB = cms.double(25.),
prefitMaxChiSqEE = cms.double(10.),
dynamicPedestalsEB = cms.bool(False),
dynamicPedestalsEE = cms.bool(False),
mitigateBadSamplesEB = cms.bool(False),
mitigateBadSamplesEE = cms.bool(False),
gainSwitchUseMaxSampleEB = cms.bool(True),
gainSwitchUseMaxSampleEE = cms.bool(False),
selectiveBadSampleCriteriaEB = cms.bool(False),
selectiveBadSampleCriteriaEE = cms.bool(False),
simplifiedNoiseModelForGainSwitch = cms.bool(True),
addPedestalUncertaintyEB = cms.double(0.),
addPedestalUncertaintyEE = cms.double(0.),
# decide which algorithm to be use to calculate the jitter
timealgo = cms.string("RatioMethod"),
# for ratio method
EBtimeFitParameters = cms.vdouble(-2.015452e+00, 3.130702e+00, -1.234730e+01, 4.188921e+01, -8.283944e+01, 9.101147e+01, -5.035761e+01, 1.105621e+01),
EEtimeFitParameters = cms.vdouble(-2.390548e+00, 3.553628e+00, -1.762341e+01, 6.767538e+01, -1.332130e+02, 1.407432e+02, -7.541106e+01, 1.620277e+01),
EBamplitudeFitParameters = cms.vdouble(1.138,1.652),
EEamplitudeFitParameters = cms.vdouble(1.890,1.400),
EBtimeFitLimits_Lower = cms.double(0.2),
EBtimeFitLimits_Upper = cms.double(1.4),
EEtimeFitLimits_Lower = cms.double(0.2),
EEtimeFitLimits_Upper = cms.double(1.4),
# for time error
EBtimeConstantTerm= cms.double(.6),
EEtimeConstantTerm= cms.double(1.0),
# for kOutOfTime flag
EBtimeNconst = cms.double(28.5),
EEtimeNconst = cms.double(31.8),
outOfTimeThresholdGain12pEB = cms.double(5), # times estimated precision
outOfTimeThresholdGain12mEB = cms.double(5), # times estimated precision
outOfTimeThresholdGain61pEB = cms.double(5), # times estimated precision
outOfTimeThresholdGain61mEB = cms.double(5), # times estimated precision
outOfTimeThresholdGain12pEE = cms.double(1000), # times estimated precision
outOfTimeThresholdGain12mEE = cms.double(1000), # times estimated precision
outOfTimeThresholdGain61pEE = cms.double(1000), # times estimated precision
outOfTimeThresholdGain61mEE = cms.double(1000), # times estimated precision
amplitudeThresholdEB = cms.double(10),
amplitudeThresholdEE = cms.double(10),
ebSpikeThreshold = cms.double(1.042),
# these are now taken from DB. Here the MC parameters for backward compatibility
ebPulseShape = cms.vdouble( 5.2e-05,-5.26e-05 , 6.66e-05, 0.1168, 0.7575, 1., 0.8876, 0.6732, 0.4741, 0.3194 ),
eePulseShape = cms.vdouble( 5.2e-05,-5.26e-05 , 6.66e-05, 0.1168, 0.7575, 1., 0.8876, 0.6732, 0.4741, 0.3194 ),
# for kPoorReco flag
kPoorRecoFlagEB = cms.bool(True),
kPoorRecoFlagEE = cms.bool(False),
chi2ThreshEB_ = cms.double(65.0),
chi2ThreshEE_ = cms.double(50.0),
# for crossCorrelationMethod
crossCorrelationStartTime = cms.double(-25),
crossCorrelationStopTime = cms.double(25),
crossCorrelationTargetTimePrecision = cms.double(0.01),
)
)