Project CMSSW displayed by LXR |
|
|
|
|||
File indexing completed on 2024-04-06 12:22:46
0001 import FWCore.ParameterSet.Config as cms 0002 0003 process = cms.Process("TEST") 0004 #process.load("MuonAnalysis.MomentumScaleCalibration.Summer08_Upsilon1S_cff") 0005 #process.load("MuonAnalysis.MomentumScaleCalibration.Summer08_InclusivePPmuX_old_cff") 0006 # process.load("MuonAnalysis.MomentumScaleCalibration.Summer09_Upsilon1S_prep_cff") 0007 #process.load("MuonAnalysis.MomentumScaleCalibration.Summer08_Z_cff") 0008 0009 process.source = cms.Source( 0010 "PoolSource", 0011 inputCommands = cms.untracked.vstring( 0012 "keep *", 0013 "drop edmGenInfoProduct_*_*_*", 0014 "drop *_TriggerResults_*_*", 0015 ), 0016 fileNames = cms.untracked.vstring( 0017 #Aligned reco 0018 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_1-5000.root', 0019 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_5001-10000.root', 0020 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_10001-15000.root', 0021 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_15001-20000.root', 0022 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_20001-25000.root', 0023 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_25001-30000.root', 0024 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_30001-35000.root', 0025 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_35001-40000.root', 0026 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_40001-45000.root', 0027 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/FullReco0Jet_45001-50000.root' 0028 #Misaligned reco 0029 # 'rfio:/castor/cern.ch/user/m/maborgia/Z_CTEQ6l/Z0jetReco/ZMuMu0jet_ReRecoWMisalignCRAFT.root' 0030 ) 0031 ) 0032 0033 0034 process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi") 0035 process.load("Geometry.CommonTopologies.globalTrackingGeometry_cfi") 0036 process.load("RecoMuon.DetLayers.muonDetLayerGeometry_cfi") 0037 process.load("Geometry.MuonNumbering.muonNumberingInitialization_cfi") 0038 process.load("RecoMuon.TrackingTools.MuonServiceProxy_cff") 0039 process.load("Alignment.CommonAlignmentProducer.GlobalPosition_Fake_cff") 0040 process.load("MagneticField.Engine.uniformMagneticField_cfi") 0041 0042 # process.source = cms.Source("PoolSource", 0043 # fileNames = cms.untracked.vstring() 0044 # ) 0045 0046 #process.poolDBESSource = cms.ESSource("PoolDBESSource", 0047 # BlobStreamerName = cms.untracked.string('TBufferBlobStreamingService'), 0048 # DBParameters = cms.PSet( 0049 # messageLevel = cms.untracked.int32(2), 0050 # authenticationPath = cms.untracked.string('/afs/cern.ch/cms/DB/conddb') 0051 # ), 0052 # timetype = cms.untracked.string('runnumber'), 0053 # connect = cms.string('sqlite_file:/data2/demattia/CMSSW_2_1_12/src/MuonAnalysis/MomentumScaleCalibration/test/dummy2.db'), 0054 # toGet = cms.VPSet(cms.PSet( 0055 # record = cms.string('MuScleFitLikelihoodPdfRcd'), 0056 # tag = cms.string('MuScleFitLikelihoodPdf_2_1_12') 0057 # )) 0058 #) 0059 0060 process.maxEvents = cms.untracked.PSet( 0061 # This are the total background events from InclusivePPmuX (89150) + 0062 # the number of Upsilon1S events. 0063 # input = cms.untracked.int32(89355) 0064 0065 input = cms.untracked.int32(5000) 0066 ) 0067 process.looper = cms.Looper( 0068 "MuScleFit", 0069 process.MuonServiceProxy, 0070 ProbabilitiesFile = cms.untracked.string('/tmp/maborgia/Probs_new_Horace_CTEQ_1000.root'), 0071 # Choose the kind of muons you want to run on 0072 # ------------------------------------------- 0073 0074 # // global muons // 0075 # MuonLabel = cms.InputTag("muons"), 0076 # MuonType = cms.int32(1), 0077 0078 # // inner track // 0079 MuonLabel = cms.InputTag("muons"), 0080 MuonType = cms.int32(10), 0081 0082 MaxMuonPt = cms.untracked.double(50.), 0083 MinMuonPt = cms.untracked.double(20.), 0084 0085 # // standalone muons // 0086 # MuonLabel = cms.InputTag("standAloneMuons:UpdatedAtVtx"), 0087 # muonType = cms.int32(2), 0088 0089 # // tracker tracks // 0090 # MuonLabel = cms.InputTag("generalTracks"), # ctfWithMaterialTracks 0091 # muonType = cms.int32(3), 0092 0093 # Output settings 0094 # --------------- 0095 OutputFileName = cms.untracked.string('MuScleFit.root'), 0096 debug = cms.untracked.int32(0), 0097 0098 # Likelihood settings 0099 # ------------------- 0100 maxLoopNumber = cms.untracked.int32(3), 0101 # Select which fits to do in which loop (0 = do not, 1 = do) 0102 doResolFit = cms.vint32( 1, 0, 0), 0103 doScaleFit = cms.vint32( 0, 1, 0), 0104 doBackgroundFit = cms.vint32( 0, 0, 0), 0105 0106 # Fit parameters and fix flags (1 = use par) 0107 # ========================================== 0108 0109 # BiasType=0 means no bias to muon momenta 0110 # ---------------------------------------- 0111 BiasType = cms.int32(0), 0112 parBias = cms.vdouble(), 0113 0114 # BiasType = 1 means linear bias on the muons Pt 0115 # the two parameters are the constant and the Pt 0116 # coefficient in this order. 0117 # ---------------------------------------------- 0118 # BiasType = cms.int32(1), 0119 # parBias = cms.vdouble(1.015, 0.001), 0120 0121 # Sinusoidal in phi 0122 # ----------------- 0123 # BiasType = 3 means sinusoidal bias on the muons Pt 0124 # the two parameters are defined by: 0125 # pt = (parScale[0] + parScale[1]*sin(phi))*pt; 0126 # BiasType = cms.int32(3), 0127 # parBias = cms.vdouble(1.015, 0.025), 0128 0129 # SmearType=0 means no smearing applied to muon momenta 0130 # ----------------------------------------------------- 0131 # SmearType = cms.int32(0), 0132 # parSmear = cms.vdouble(), 0133 0134 # =======> To use with resolution function type 15 0135 SmearType = cms.int32(6), 0136 parSmear = cms.vdouble(-0.00027357, 0., 0.000319814, 0137 0., 0., 0.00508154, 0138 0., 0., 1.18468, 0.0756988, 0139 0., -0.0349113, 1.13643, 0.0719739, 0140 1.0008, 0.050805), 0141 0142 # ------------------------- # 0143 # Resolution fit parameters # 0144 # ------------------------- # 0145 0146 # ------------------------------------------------- # 0147 # New resolution function derived for low Pt region # 0148 # ------------------------------------------------- # 0149 # The eleven parResol parameters of resolfittype=11 are respectively: 0150 #"offsetEtaCentral", "offsetEtaHigh", "coeffOverPt", "coeffHighPt", "linaerEtaCentral", "parabEtaCentral", "linaerEtaHigh", "parabEtaHigh" }; 0151 #ResolFitType = cms.int32(11), #inner tracks in 31X 0152 #parResol = cms.vdouble(-0.986, -0.986, -0.04, -0.038, -0.0014, 0.006, -0.0025, 0.0185), 0153 #parResolFix = cms.vint32(0, 0, 0, 0, 0, 0, 0, 0), 0154 #parResolOrder = cms.vint32(0, 0, 0, 0, 0, 0, 0, 0), 0155 0156 # ResolFitType = cms.int32(12), #inner tracks in 31X 0157 # parResol = cms.vdouble(-0.986, -0.986, -0.04, -0.038, 0158 # -0.0014, 0.006, -0.0025, 0.0185, 0159 # -0.0014, 0., -0.001), 0160 # parResolFix = cms.vint32(0, 0, 0, 0, 0161 # 0, 0, 0, 0, 0162 # 0, 0, 0), 0163 # parResolOrder = cms.vint32(0, 0, 0, 0, 0164 # 0, 0, 0, 0, 0165 # 0, 0, 0), 0166 0167 # ResolFitType = cms.int32(8), 0168 # parResol = cms.vdouble(0.00, 0.0, 1.38, 0.051, 0169 # 0.00043, 0.0041, 0.0000028, 0.000077, 0170 # 0.00011, 0.0018, -0.00000094, 0.000022), 0171 # parResolFix = cms.vint32(1, 1, 0, 0, 0172 # 1, 1, 1, 1, 0173 # 1, 1, 1, 1), 0174 # parResolOrder = cms.vint32(0, 0, 0, 0, 0175 # 0, 0, 0, 0, 0176 # 0, 0, 0, 0), 0177 0178 #### Resolution function for misaligned data Type 15 0179 ResolFitType = cms.int32(15), 0180 parResol = cms.vdouble(0.02, 0.0, 0.00014, 0181 0.0, 0.0, 0.005759, 0182 0.0, 0.0, 1.38, 0.114, 0183 0.0, 0.0, 1.4856, 0.0954 0184 ), 0185 parResolFix = cms.vint32(0, 1, 0, 0186 1, 1, 0, 0187 1, 1, 0, 0, 0188 1, 0, 0, 0), 0189 parResolOrder = cms.vint32(0, 0, 0, 0190 0, 0, 0, 0191 0, 0, 0, 0, 0192 0, 0, 0, 0), 0193 0194 # -------------------- # 0195 # Scale fit parameters # 0196 # -------------------- # 0197 0198 # Fit a linear Pt scale correction with parameters: 0199 # Pt scale and Pt slope. 0200 # ------------------------------------------------- 0201 # ScaleFitType = cms.int32(1), 0202 # parScaleOrder = cms.vint32(0,0), 0203 # parScaleFix = cms.vint32(0,0), 0204 # parScale = cms.vdouble(1.0, 0.0), 0205 0206 # Scale fit type=14: Pt offset and grade up to three, Eta terms up to the sixth grade 0207 # ----------------------------------------------------------------------------------- 0208 ScaleFitType = cms.int32(14), 0209 parScaleOrder = cms.vint32(0, # scale 0210 0,0,0, # pt up to grade 3 0211 0,0,0,0,0,0), # eta up to grade 6 0212 parScaleFix = cms.vint32(0, 0213 0,1,1, 0214 1,1,1,1,1,1), 0215 parScale = cms.vdouble(1.0, 0216 0., 0., 0., 0217 0., 0., 0., 0., 0., 0.), 0218 0219 # parScale = cms.vdouble(1.0, 0220 # -0.000000315315, 0., 0., 0221 # 0.0000147547, -0.00000836992, 0., 0., 0., 0.), 0222 # parScale = cms.vdouble(1.0, 0223 # -0.00041991, 0., 0., 0224 # 0.000727967, -0.00082597, 0., 0., 0., 0.), 0225 0226 # ------------------------- # 0227 # Background fit parameters # 0228 # ------------------------- # 0229 0230 # Window factors for: Z, Upsilons and (J/Psi,Psi2S) regions 0231 LeftWindowFactor = cms.vdouble(1., 5., 3.), 0232 RightWindowFactor = cms.vdouble(1., 5., 3.), 0233 0234 # The parameter of BgrFitType=1 is the bgr fraction 0235 # ------------------------------------------------- 0236 # BgrFitType = cms.int32(1), 0237 # parBgrFix = cms.vint32(0), 0238 # parBgr = cms.vdouble(0.001), 0239 # parBgrOrder = cms.vint32(0), 0240 0241 # The two parameters of BgrFitType=2 are respectively: 0242 # bgr fraction, (negative of) bgr exp. slope, bgr constant 0243 # -------------------------------------------------------- 0244 # The function types for resonances in a region must be the same 0245 BgrFitType = cms.vint32(2, 2, 2), # resonances 0246 # parBgr = cms.vdouble(0.05, 0.001), 0247 parBgr = cms.vdouble(0., 0., 0., 0., 0., 0., 0248 0., 0., 0., 0., 0., 0., 0.,0., 0.,0., 0.,0.), 0249 # parBgr = cms.vdouble(0., 0., 0.9, 0.0001, 0., 0., 0250 # 0., 0., 0., 0., 0., 0., 0.,0., 0.,0., 0.,0.), 0251 parBgrFix = cms.vint32(0, 0, 0, 0, 0, 0, 0252 # The rest of the parameters is used for the resonance regions. They are automatically fixed in the code 0253 # because they are never used to fit the background, but only after the rescaling. 0254 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1), 0255 parBgrOrder = cms.vint32(0, 0, 0, 0, 0, 0, 0256 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0), 0257 0258 # ---------------- # 0259 # Select resonance # 0260 # ---------------- # 0261 0262 # The resonances are to be specified in this order: 0263 # Z0, Y(3S), Y(2S), Y(1S), Psi(2S), J/Psi 0264 # ------------------------------------------------- 0265 resfind = cms.vint32(1, 0, 0, 0, 0, 0), 0266 FitStrategy = cms.int32(2), 0267 0268 speedup = cms.bool(False), 0269 # speedup = cms.bool(True), 0270 # Set this to false if you do not want to use simTracks. 0271 # (Note that this is skipped anyway if speedup == True). 0272 # compareToSimTracks = cms.bool(True), 0273 compareToSimTracks = cms.bool(False), 0274 Sherpa = cms.untracked.bool(True), 0275 0276 ) 0277 0278 # Timing information 0279 # process.load("FWCore.MessageLogger.MessageLogger_cfi") 0280 # TimingLogFile = cms.untracked.string('timing.log') 0281 # process.Timing = cms.Service("Timing") 0282 0283
| [ Source navigation ] | [ Diff markup ] | [ Identifier search ] | [ general search ] |
|
This page was automatically generated by the 2.2.1 LXR engine. The LXR team |
|