MomentumScaleCalibration/test/ResolutionAnalyzer_cfg.py

0001 # -*- coding: utf-8 -*-
0002 import FWCore.ParameterSet.Config as cms
0003
0004 process = cms.Process("RESOLUTIONANALYZER")
0005 process.load("MuonAnalysis.MomentumScaleCalibration.local_CSA08_JPsi_cff")
0006
0007 process.load("Configuration.StandardSequences.MagneticField_cff")
0008
0009 process.load("Geometry.CMSCommonData.cmsIdealGeometryXML_cfi")
0010
0011 process.load("Geometry.CommonTopologies.globalTrackingGeometry_cfi")
0012
0013 process.load("RecoMuon.DetLayers.muonDetLayerGeometry_cfi")
0014
0015 process.load("Geometry.MuonNumbering.muonNumberingInitialization_cfi")
0016
0017 process.load("RecoMuon.TrackingTools.MuonServiceProxy_cff")
0018
0019 # process.source = cms.Source("PoolSource",
0020 #     fileNames = cms.untracked.vstring()
0021 # )
0022
0023 process.maxEvents = cms.untracked.PSet(
0024     input = cms.untracked.int32(-1)
0025 )
0026 process.ResolutionAnalyzerModule = cms.EDAnalyzer(
0027     "ResolutionAnalyzer",
0028     process.MuonServiceProxy,
0029     # // standalone muons
0030     # InputTag MuonLabel = standAloneMuons:UpdatedAtVtx
0031     # int32 muonType = 2
0032     # // tracker tracks
0033     # InputTag MuonLabel = generalTracks //ctfWithMaterialTracks
0034     # int32 muonType = 3
0035
0036     # No fit is done, but the resolution function needs parameters
0037     # ResolFitType = cms.int32(6),
0038     # parResol = cms.vdouble(0.002, -0.0015, 0.000056, -0.00000085, 0.0000000046, -0.000027, 0.0037,
0039     #                        0.005, 0.00027, 0.0000027, 0.000094,
0040     #                        0.002, 0.00016, -0.00000051, 0.000022),
0041
0042     # The eleven parResol parameters of resolfittype=8 are respectively:
0043     # constant of sigmaPt, Pt dep. of sigmaPt,
0044     # scale of the eta dep. made by points with values derived from MuonGun.
0045     # constant of sigmaCotgTheta, 1/Pt dep. of sigmaCotgTheta, Eta dep. of
0046     # sigmaCotgTheta, Eta^2 dep of sigmaCotgTheta;
0047     # constant of sigmaPhi, 1/Pt dep. of sigmaPhi, Eta dep. of sigmaPhi,
0048     # Eta^2 dep. of sigmaPhi.
0049     # This parameters are taken directly from the MuonGun (5<Pt<100, |eta|<3)
0050     ResolFitType = cms.int32(8),
0051     parResol = cms.vdouble(-0.007, 0.0001, 1.0,
0052                            0.00043, 0.0041, 0.0000028, 0.000077,
0053                            0.00011, 0.0018, -0.00000094, 0.000022),
0054     parResolFix = cms.vint32(0, 0, 0,
0055                              0, 0, 0, 0,
0056                              0, 0, 0, 0),
0057     parResolOrder = cms.vint32(0, 0, 0,
0058                                1, 1, 1, 1,
0059                                2, 2, 2, 2),
0060
0061
0062     # The resonances are to be specified in this order:
0063     # Z0, Y(3S), Y(2S), Y(1S), Psi(2S), J/Psi
0064     # -------------------------------------------------
0065     ResFind = cms.vint32(0, 0, 0, 0, 0, 1),
0066
0067     # Tree settings
0068     MaxEvents = cms.uint32(-1),
0069     InputTreeName = cms.string("tree.root"),
0070
0071     # Output settings
0072     # ---------------
0073     OutputFileName = cms.untracked.string('ResolutionAnalyzer_JPsi.root'),
0074     Debug = cms.untracked.bool(False),
0075     # Choose the kind of muons you want to run on
0076     # -------------------------------------------
0077     # global muons
0078     # MuonType = cms.int32(1),
0079     # MuonLabel = cms.InputTag("muons"),
0080     # // standalone muons
0081     # MuonType = cms.int32(2),
0082     # MuonLabel = cms.InputTag("standAloneMuons:UpdatedAtVtx"),
0083     # tracker tracks
0084     MuonType = cms.int32(3),
0085     MuonLabel = cms.InputTag("generalTracks"),
0086     Resonance = cms.untracked.bool(True),
0087     ReadCovariances = cms.untracked.bool(False),
0088     # This is used only when the ReadCovariances bool == True
0089     InputFileName = cms.untracked.string('ResolutionAnalyzer_JPsi.root')
0090 )
0091
0092 process.p1 = cms.Path(process.ResolutionAnalyzerModule)
0093