TopKinFitter/python/TtFullHadKinFitProducer_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 # current defaults are resolutions for PF jets derived from Summer11 MC;
0004 # values of udscResolutions and bResolutions (see below) can still be replaced in the user's cff or cfg;
0005 # for the unit tests, we don't mind that CALO jets will be used with PF resolutions
0006 from TopQuarkAnalysis.TopObjectResolutions.stringResolutions_etEtaPhi_Summer11_cff import *
0007
0008 kinFitTtFullHadEvent = cms.EDProducer("TtFullHadKinFitProducer",
0009     jets = cms.InputTag("selectedPatJets"),
0010
0011     # ------------------------------------------------
0012     # maximum number of jets to be considered in the
0013     # jet combinatorics (has to be >= 6, can be set to
0014     # -1 if you want to take all)
0015     # ------------------------------------------------
0016     maxNJets = cms.int32(6),
0017
0018     #-------------------------------------------------
0019     # maximum number of jet combinations finally
0020     # written into the event, starting from the "best"
0021     # (has to be >= 1, can be set to -1 if you want to
0022     # take all)
0023     #-------------------------------------------------
0024     maxNComb = cms.int32(1),
0025
0026     # ------------------------------------------------
0027     # option to take only a given jet combination
0028     # instead of going through the full combinatorics
0029     # ------------------------------------------------
0030     match = cms.InputTag(""),
0031     useOnlyMatch = cms.bool(False),
0032
0033     # ------------------------------------------------
0034     # option to use b-tagging
0035     # ------------------------------------------------
0036     useBTagging         = cms.bool(True),
0037     bTagAlgo            = cms.string("trackCountingHighPurBJetTags"),
0038     minBTagValueBJet    = cms.double(3.0),
0039     maxBTagValueNonBJet = cms.double(3.0),
0040     bTags               = cms.uint32(2), # minimal number of b-tagged
0041                                          # jets, if more are available
0042                                          # they will be used
0043
0044     # ------------------------------------------------
0045     ## specify jet correction level as, Uncorrected,
0046     ## L1Offset, L2Relative, L3Absolute, L4Emf,
0047     ## L5Hadron, L6UE, L7Parton, a flavor specifica-
0048     ## tion will be added automatically, when chosen
0049     # ------------------------------------------------
0050     jetCorrectionLevel = cms.string("L3Absolute"),
0051
0052     # ------------------------------------------------
0053     # settings for the KinFitter
0054     # ------------------------------------------------
0055     maxNrIter = cms.uint32(500),
0056     maxDeltaS = cms.double(5e-05),
0057     maxF      = cms.double(0.0001),
0058
0059     # ------------------------------------------------
0060     # select parametrisation
0061     # 0: EMom, 1: EtEtaPhi, 2: EtThetaPhi
0062     # ------------------------------------------------
0063     jetParametrisation = cms.uint32(1),
0064
0065     # ------------------------------------------------
0066     # set constraints
0067     # 1: W1-mass, 2: W2-mass
0068     # 3: t1-mass, 4: t2-mass
0069     # 5: equal top-masses
0070     # ------------------------------------------------
0071     constraints = cms.vuint32(1, 2, 5),
0072
0073     # ------------------------------------------------
0074     # set mass values used in the constraints
0075     # ------------------------------------------------
0076     mW   = cms.double(80.4),
0077     mTop = cms.double(173.),
0078
0079     # ------------------------------------------------
0080     # resolutions used for the kinematic fit
0081     # (see also comments at the head of this file)
0082     # ------------------------------------------------
0083     udscResolutions = udscResolutionPF.functions,
0084     bResolutions    = bjetResolutionPF.functions,
0085
0086     # ------------------------------------------------
0087     # set correction factor(s) for the jet energy resolution:
0088     # - (optional) eta dependence assumed to be symmetric
0089     #   around eta=0, i.e. parametrized in |eta|
0090     # - any negative value as last bin edge is read as "inf"
0091     # - make sure that number of entries in vector with
0092     #   bin edges = number of scale factors + 1
0093     # ------------------------------------------------
0094     jetEnergyResolutionScaleFactors = cms.vdouble(1.0),
0095     jetEnergyResolutionEtaBinning = cms.vdouble(0.0,-1.0)
0096 )
0097
0098