File indexing completed on 2024-04-06 12:22:46
0001 import FWCore.ParameterSet.Config as cms
0002
0003 class SetParameters:
0004 def __init__(self):
0005 self.parFix = cms.vint32()
0006 self.parOrder = cms.vint32()
0007 self.par = cms.vdouble()
0008 self.parStep = cms.untracked.vdouble()
0009 self.parMin = cms.untracked.vdouble()
0010 self.parMax = cms.untracked.vdouble()
0011
0012 def set(self, fix, order, value, step, min, max):
0013 self.parFix.append(fix)
0014 self.parOrder.append(order)
0015 self.par.append(value)
0016 self.parStep.append(step)
0017 self.parMin.append(min)
0018 self.parMax.append(max)
0019
0020
0021 resolSetter = SetParameters()
0022
0023 resolSetter.set( 0 ,0, 0.00354277 , 0.0001, 0.,0. )
0024 resolSetter.set( 0 ,0, 0.00023585 , 0.00001, 0.,0. )
0025 resolSetter.set( 1 ,0, 0 , 0.000001, 0.,0. )
0026 resolSetter.set( 0 ,0, 0.00654324 , 0.0001, 0.,0. )
0027 resolSetter.set( 1 ,0, 0 , 0.00001, 0.,0. )
0028 resolSetter.set( 1 ,0, 0 , 0.0001, 0.,0. )
0029 resolSetter.set( 1 ,0, 0.0418539 , 0.0001, 0.,0. )
0030 resolSetter.set( 1 ,0, -2.4 , 0.0001, 0.,0. )
0031 resolSetter.set( 1 ,0, 0 , 0.00001, 0.,0. )
0032 resolSetter.set( 1 ,0, 0 , 0.0001, 0.,0. )
0033 resolSetter.set( 0 ,0, 0.00526402 , 0.001, 0.,0. )
0034 resolSetter.set( 1 ,0, 0 , 0.001, 0.,0. )
0035 resolSetter.set( 1 ,0, -2.4 , 0.0001, 0.,0. )
0036 resolSetter.set( 1 ,0, 0 , 0.001, 0.,0. )
0037 resolSetter.set( 1 ,0, 0 , 0.0001, 0.,0. )
0038
0039
0040
0041
0042
0043
0044
0045
0046
0047
0048
0049 scaleSetter = SetParameters()
0050
0051 scaleSetter.set( 1 ,0, -6.08453e-05 , 0.0000001, 0.,0. )
0052 scaleSetter.set( 1 ,0, 0. , 0.00001, 0.,0. )
0053 scaleSetter.set( 1 ,0, 0.00174819 , 0.0001, -2.,2. )
0054 scaleSetter.set( 1 ,0, 0.000120695 , 0.00001, 0.,0. )
0055 scaleSetter.set( 1 ,0, 0.254793 , 0.001, -3.1416,3.1416 )
0056
0057 process = cms.Process("TEST")
0058
0059
0060
0061
0062
0063
0064
0065
0066
0067
0068
0069
0070 process.source = cms.Source("EmptySource")
0071 process.maxEvents = cms.untracked.PSet(
0072 input = cms.untracked.int32(0)
0073 )
0074
0075 process.looper = cms.Looper(
0076 "MuScleFit",
0077
0078
0079 MaxEventsFromRootTree = cms.int32(-1),
0080
0081
0082
0083
0084 InputRootTreeFileName = cms.string("/home/castello/7TeV/CMSSW_3_8_5_patch3/src/Tree/Tree_2010AB_INNtk_ICHEPgeom_BS.root"),
0085
0086
0087
0088 OutputRootTreeFileName = cms.string("tree_data_newbinning.root"),
0089
0090
0091
0092 MuonLabel = cms.InputTag("patMuons"),
0093
0094
0095
0096
0097
0098
0099
0100
0101
0102
0103
0104 MuonType = cms.int32(2),
0105
0106
0107
0108 PATmuons = cms.untracked.bool(True),
0109
0110
0111
0112
0113
0114
0115
0116 resfind = cms.vint32(1, 0, 0, 0, 0, 0),
0117
0118
0119
0120 maxLoopNumber = cms.untracked.int32(2),
0121
0122 doResolFit = cms.vint32(0,0),
0123 doScaleFit = cms.vint32(1,0),
0124 doBackgroundFit = cms.vint32(0,0),
0125 doCrossSectionFit = cms.vint32(0,0),
0126
0127
0128
0129 UseProbsFile = cms.untracked.bool(True),
0130
0131
0132 speedup = cms.bool(True),
0133
0134
0135 compareToSimTracks = cms.bool(False),
0136
0137
0138
0139 OutputFileName = cms.untracked.string('MuScleFit_data_newbinning.root'),
0140
0141
0142
0143 BiasType = cms.int32(0),
0144 parBias = cms.vdouble(),
0145
0146
0147
0148 SmearType = cms.int32(0),
0149 parSmear = cms.vdouble(),
0150
0151
0152 ResolFitType = cms.int32(42),
0153 parResolFix = resolSetter.parFix,
0154 parResolOrder = resolSetter.parOrder,
0155 parResol = resolSetter.par,
0156 parResolStep = resolSetter.parStep,
0157 parResolMin = resolSetter.parMin,
0158 parResolMax = resolSetter.parMax,
0159
0160
0161
0162
0163
0164
0165
0166 ScaleFitType = cms.int32(29),
0167 parScaleFix = scaleSetter.parFix,
0168 parScaleOrder = scaleSetter.parOrder,
0169 parScale = scaleSetter.par,
0170 parScaleStep = scaleSetter.parStep,
0171 parScaleMin = scaleSetter.parMin,
0172 parScaleMax = scaleSetter.parMax,
0173
0174
0175
0176
0177
0178
0179
0180
0181 parCrossSectionOrder = cms.vint32(0, 0, 0, 0, 0, 0),
0182 parCrossSectionFix = cms.vint32(0, 0, 0, 0, 0, 0),
0183 parCrossSection = cms.vdouble(1.233, 2.07, 6.33, 13.9, 2.169, 127.2),
0184
0185
0186
0187
0188
0189
0190 LeftWindowBorder = cms.vdouble(70., 8., 1.391495),
0191 RightWindowBorder = cms.vdouble(110., 12., 5.391495),
0192
0193
0194
0195
0196
0197 BgrFitType = cms.vint32(2, 2, 2),
0198
0199 parBgr = cms.vdouble(0., 0., 0., 0., 0., 0.,
0200 0., 0., 0., 0., 0., 0., 0.,0., 0.,0., 0.,0.),
0201 parBgrFix = cms.vint32(0, 0, 0, 0, 0, 0,
0202
0203
0204 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1),
0205 parBgrOrder = cms.vint32(0, 0, 0, 0, 0, 0,
0206 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0),
0207
0208
0209
0210
0211
0212 FitStrategy = cms.int32(2),
0213
0214
0215 StartWithSimplex = cms.bool(True),
0216 ComputeMinosErrors = cms.bool(False),
0217 MinimumShapePlots = cms.bool(True),
0218
0219
0220
0221
0222
0223
0224
0225
0226
0227
0228
0229
0230 RapidityBinsForZ = cms.untracked.bool(False),
0231 ProbabilitiesFile = cms.untracked.string("/home/castello/Development/test/Probs_1001.root"),
0232
0233
0234 TriggerResultsLabel = cms.untracked.string("TriggerResults"),
0235 TriggerResultsProcess = cms.untracked.string("HLT"),
0236
0237
0238 TriggerPath = cms.untracked.string("All"),
0239
0240 NegateTrigger = cms.untracked.bool(False),
0241
0242 debug = cms.untracked.int32(0),
0243 )