MomentumScaleCalibration/test/template_MuScleFit_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("TEST")
0004
0005 # process.source = cms.Source(
0006 #     "PoolSource",
0007 #     fileNames = cms.untracked.vstring(
0008 #       "file:/home/demattia/MuScleFit/PatMuons/onia2MuMuPAT_Summer10-DESIGN_36_V8-X0MAX-v2.root"
0009 #     )
0010 # )
0011 # process.maxEvents = cms.untracked.PSet(
0012 #     input = cms.untracked.int32(-1)
0013 # )
0014
0015 # Use this when running on a tree
0016 process.source = cms.Source("EmptySource")
0017 process.maxEvents = cms.untracked.PSet(
0018     input = cms.untracked.int32(0)
0019 )
0020
0021 process.looper = cms.Looper(
0022     "MuScleFit",
0023
0024     # Only used when reading events from a root tree
0025     MaxEventsFromRootTree = cms.int32(-1),
0026     # Specify a file if you want to read events from a root tree in a local file.
0027     # In this case the input source should be an empty source with 0 events.
0028
0029     #InputRootTreeFileName = cms.string("/home/castello/7TeV/CMSSW_3_8_5_patch3/src/Tree/Tree_2010AB_INNtk_EOYgeom_BS.root"),
0030     #InputRootTreeFileName = cms.string("/home/castello/7TeV/CMSSW_3_8_5_patch3/src/Tree/Fall10/Tree_MCFall2010_INNtk_CRAFTRealistic_wGEN.root"),
0031     InputRootTreeFileName = cms.string("/home/castello/7TeV/CMSSW_3_8_5_patch3/src/Tree/Tree_2010AB_INNtk_ICHEPgeom_BS.root"),
0032
0033     # Specify the file name where you want to save a root tree with the muon pairs.
0034     # Leave empty if no file should be written.
0035     OutputRootTreeFileName = cms.string(""), #Tree_38X_INNtk_resol42_test7.root
0036
0037     # Choose the kind of muons you want to run on
0038     # -------------------------------------------
0039     MuonLabel = cms.InputTag("patMuons"),
0040     # Defines what type of muons to use:
0041     # -1 = onia guys selection
0042     # -2 = onia guys selection - only GG
0043     # -3 = onia guys selection - only GT
0044     # -4 = onia guys selection - only TT
0045     # Note that the above samples are independent and represent the composition of the inclusive sample
0046     # 1 = global muon
0047     # 2 = standalone muon
0048     # 3 = tracker muon
0049     # 4 = calo muon
0050     # 10 = innerTrack of global muon
0051     MuonType = cms.int32(-1),
0052
0053     # This line allows to switch to PAT muons. Default is false.
0054     # Note that the onia selection works only with onia patTuples.
0055     PATmuons = cms.untracked.bool(True),
0056
0057     # ---------------- #
0058     # Select resonance #
0059     # ---------------- #
0060     # The resonances are to be specified in this order:
0061     # Z0, Y(3S), Y(2S), Y(1S), Psi(2S), J/Psi
0062     # -------------------------------------------------
0063     resfind = cms.vint32(1, 0, 0, 0, 0, 0),
0064
0065     # Likelihood settings
0066     # -------------------
0067     maxLoopNumber = cms.untracked.int32(4),
0068     # Select which fits to do in which loop (0 = do not, 1 = do)
0069     doResolFit =        cms.vint32(1,0,1,0),
0070     doScaleFit =        cms.vint32(0,1,0,0),
0071     doBackgroundFit =   cms.vint32(0,0,0,0),
0072     doCrossSectionFit = cms.vint32(0,0,0,0),
0073
0074
0075
0076
0077     # Use the probability file or not. If not it will perform a simpler selection taking the muon pair with
0078     # invariant mass closer to the pdf value and will crash if some fit is attempted.
0079     UseProbsFile = cms.untracked.bool(True),
0080
0081     # False = use also MC information
0082     speedup = cms.bool(True),
0083     # Set this to false if you do not want to use simTracks.
0084     # (Note that this is skipped anyway if speedup == True).
0085     compareToSimTracks = cms.bool(False),
0086
0087     # Output settings
0088     # ---------------
0089     OutputFileName = cms.untracked.string('MuScleFit.root'),
0090
0091
0092     # Fit parameters and fix flags (1 = use par)
0093     # ==========================================
0094
0095
0096     # BiasType=0 means no bias to muon momenta
0097     # ----------------------------------------
0098     BiasType = cms.int32(0),
0099     parBias = cms.vdouble(),
0100
0101     # SmearType=0 means no smearing applied to muon momenta
0102     # -----------------------------------------------------
0103     SmearType = cms.int32(0),
0104     parSmear = cms.vdouble(),
0105
0106
0107 ### taken from J/Psi #########################
0108     ResolFitType = cms.int32(20),
0109     parResol = cms.vdouble(),
0110
0111     parResolFix = cms.vint32(0, 0, 1, 0,
0112                              0, 0, 0, 0,
0113                              1, 0, 0, 0,
0114                              0, 0, 1),
0115     parResolOrder = cms.vint32(0, 0, 0, 0,
0116                                0, 0, 0, 0,
0117                                0, 0, 0, 0,
0118                                0, 0, 0),
0119
0120     # -------------------- #
0121     # Scale fit parameters #
0122     # -------------------- #
0123
0124     # -----------------------------------------------------------------------------------
0125     ScaleFitType = cms.int32(29),
0126     parScaleOrder = cms.vint32( 0, 0, 0, 0, 0),
0127     parScaleFix =   cms.vint32(0, 1, 0, 0, 0),
0128     parScale = cms.vdouble(0., 0., 0., 0., 0.),
0129
0130     # Scale fit type=11: Linear in pt, sinusoidal in phi with muon sign -->GOOD results in phi
0131     ##  modified for mu+/mu -
0132     # -----------------------------------------------------------------------------------
0133     ##    ScaleFitType = cms.int32(11),
0134     ##     parScaleOrder = cms.vint32(0, 0, 0, 0, 0, 0, 0, 0),
0135     ##     parScaleFix =   cms.vint32(0, 0, 0, 0, 0, 0, 0, 0),
0136     ##     parScale = cms.vdouble(1., 0., 0., 1., 0., 0., 1., 0.),
0137
0138
0139     # ---------------------------- #
0140     # Cross section fit parameters #
0141     # ---------------------------- #
0142     # Note that the cross section fit works differently than the others, it
0143     # fits ratios of parameters. Fix and Order should not be used as is, they
0144     # are there mainly for compatibility.
0145     parCrossSectionOrder = cms.vint32(0, 0, 0, 0, 0, 0),
0146     parCrossSectionFix =   cms.vint32(0, 0, 0, 0, 0, 0),
0147     parCrossSection = cms.vdouble(1.233, 2.07, 6.33, 13.9, 2.169, 127.2),
0148
0149     # ------------------------- #
0150     # Background fit parameters #
0151     # ------------------------- #
0152
0153     # Window factors for: Z, Upsilons and (J/Psi,Psi2S) regions
0154     LeftWindowBorder = cms.vdouble(70., 8., 1.391495),
0155     RightWindowBorder = cms.vdouble(110., 12., 5.391495),
0156
0157     # The two parameters of BgrFitType=2 are respectively:
0158     # bgr fraction, (negative of) bgr exp. slope, bgr constant
0159     # --------------------------------------------------------
0160     # The function types for resonances in a region must be the same
0161     BgrFitType = cms.vint32(2, 2, 2), # regions
0162     # These empty parameters should be used when there is no background
0163     parBgr = cms.vdouble(0., 0.,   0., 0.,   0., 0.,
0164                          0., 0.,   0., 0.,   0., 0.,   0.,0.,   0.,0.,   0.,0.),
0165     parBgrFix = cms.vint32(0, 0,   0, 0,   0, 0,
0166                            # The rest of the parameters is used for the resonance regions. They are automatically fixed in the code
0167                            # because they are never used to fit the background, but only after the rescaling.
0168                            1, 1,   1, 1,   1, 1,   1, 1,   1, 1,   1, 1),
0169     parBgrOrder = cms.vint32(0, 0,   0, 0,   0, 0,
0170                              0, 0,   0, 0,   0, 0,   0, 0,   0, 0,   0, 0),
0171
0172
0173     # ----------------------- #
0174
0175     # Set Minuit fit strategy
0176     FitStrategy = cms.int32(2),
0177
0178     # Fit accuracy and debug parameters
0179     StartWithSimplex = cms.bool(True),
0180     ComputeMinosErrors = cms.bool(False),
0181     MinimumShapePlots = cms.bool(True),
0182
0183     # Set the cuts on muons to be used in the fit
0184
0185     # MinMuonPt = cms.untracked.double(30.),
0186     # MaxMuonPt = cms.untracked.double(50.),
0187     # MinMuonEtaFirstRange = cms.untracked.double(-0.8),
0188     # MaxMuonEtaFirstRange = cms.untracked.double(0.8),
0189     # MinMuonEtaSecondRange = cms.untracked.double(-0.8),
0190     # MaxMuonEtaSecondRange = cms.untracked.double(0.8),
0191
0192     # ProbabilitiesFileInPath = cms.untracked.string("MuonAnalysis/MomentumScaleCalibration/test/Probs_merge.root"),
0193     ProbabilitiesFile = cms.untracked.string("/home/castello/7TeV/CMSSW_3_8_5_patch3/src/MuonAnalysis/MomentumScaleCalibration/test/Probs_merge.root"),
0194
0195     # The following parameters can be used to filter events
0196     TriggerResultsLabel = cms.untracked.string("TriggerResults"),
0197     TriggerResultsProcess = cms.untracked.string("HLT"),
0198     # TriggerPath: "" = No trigger requirements, "All" = No specific path
0199     #TriggerPath = cms.untracked.string("HLT_L1MuOpen"),
0200     TriggerPath = cms.untracked.string("All"),
0201     # Negate the result of the trigger
0202     NegateTrigger = cms.untracked.bool(False),
0203
0204     debug = cms.untracked.int32(0),
0205 )