Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:31:41

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