CMS 3D CMS Logo

/data/refman/pasoursint/CMSSW_5_3_10_patch2/src/TopQuarkAnalysis/Configuration/python/patRefSel_refMuJets.py

Go to the documentation of this file.
00001 #
00002 # This file contains the Top PAG reference selection for mu + jets analysis.
00003 #
00004 
00005 ### Muon configuration
00006 
00007 # PAT muons
00008 muonsUsePV     = False # use beam spot rather than PV, which is necessary for 'dB' cut
00009 muonEmbedTrack = True  # embedded track needed for muon ID cuts
00010 
00011 ### Jet configuration
00012 
00013 # Jet algorithm
00014 jetAlgo = 'AK5'
00015 
00016 # JEC sets
00017 jecSetBase = jetAlgo
00018 
00019 
00020 ### ------------------------- Reference selection -------------------------- ###
00021 
00022 
00023 # PF2PAT settings
00024 from TopQuarkAnalysis.Configuration.patRefSel_PF2PAT import *
00025 
00026 ### Trigger selection
00027 
00028 # HLT selection
00029 triggerSelectionDataRelVals = 'HLT_IsoMu17_eta2p1_TriCentralPFJet30_v*' # 2012A RelVals
00030 triggerSelectionData        = 'HLT_IsoMu17_eta2p1_TriCentralPFNoPUJet30_30_20_v*'
00031 triggerSelectionMC          = 'HLT_*' # not recommended
00032 
00033 ### Muon selection
00034 
00035 # Minimal selection for all muons, also basis for signal and veto muons
00036 muonCut  =     'isPFMuon'                                                                      # general reconstruction property
00037 muonCut += ' && (isGlobalMuon || isTrackerMuon)'                                               # general reconstruction property
00038 muonCut += ' && pt > 10.'                                                                      # transverse momentum
00039 muonCut += ' && abs(eta) < 2.5'                                                                # pseudo-rapisity range
00040 muonCut += ' && (chargedHadronIso+neutralHadronIso+photonIso-0.5*puChargedHadronIso)/pt < 0.2' # relative isolation w/ Delta beta corrections (factor 0.5)
00041 
00042 # Signal muon selection on top of 'muonCut'
00043 signalMuonCut  =     'isPFMuon'                                                                       # general reconstruction property
00044 signalMuonCut += ' && isGlobalMuon'                                                                   # general reconstruction property
00045 signalMuonCut += ' && pt > 26.'                                                                       # transverse momentum
00046 signalMuonCut += ' && abs(eta) < 2.1'                                                                 # pseudo-rapisity range
00047 signalMuonCut += ' && globalTrack.normalizedChi2 < 10.'                                               # muon ID: 'isGlobalMuonPromptTight'
00048 signalMuonCut += ' && track.hitPattern.trackerLayersWithMeasurement > 5'                              # muon ID: 'isGlobalMuonPromptTight'
00049 signalMuonCut += ' && globalTrack.hitPattern.numberOfValidMuonHits > 0'                               # muon ID: 'isGlobalMuonPromptTight'
00050 signalMuonCut += ' && abs(dB) < 0.2'                                                                  # 2-dim impact parameter with respect to beam spot (s. "PAT muon configuration" above)
00051 signalMuonCut += ' && innerTrack.hitPattern.numberOfValidPixelHits > 0'                               # tracker reconstruction
00052 signalMuonCut += ' && numberOfMatchedStations > 1'                                                    # muon chamber reconstruction
00053 signalMuonCut += ' && (chargedHadronIso+neutralHadronIso+photonIso-0.5*puChargedHadronIso)/pt < 0.12' # relative isolation w/ Delta beta corrections (factor 0.5)
00054 
00055 muonVertexMaxDZ = 0.5 # DeltaZ between muon vertex and PV
00056 
00057 ### Jet selection
00058 
00059 # Signal jet selection
00060 jetCut  =     'abs(eta) < 2.5'                                        # pseudo-rapisity range
00061 jetCut += ' && numberOfDaughters > 1'                                 # PF jet ID:
00062 jetCut += ' && neutralHadronEnergyFraction < 0.99'                    # PF jet ID:
00063 jetCut += ' && neutralEmEnergyFraction < 0.99'                        # PF jet ID:
00064 jetCut += ' && (chargedEmEnergyFraction < 0.99 || abs(eta) >= 2.4)'   # PF jet ID:
00065 jetCut += ' && (chargedHadronEnergyFraction > 0. || abs(eta) >= 2.4)' # PF jet ID:
00066 jetCut += ' && (chargedMultiplicity > 0 || abs(eta) >= 2.4)'          # PF jet ID:
00067 # varying jet pt thresholds
00068 veryLooseJetCut = 'pt > 20.' # transverse momentum (all jets)
00069 looseJetCut     = 'pt > 35.' # transverse momentum (3rd jet, optional)
00070 tightJetCut     = 'pt > 45.' # transverse momentum (leading jets)
00071 
00072 ### Electron selection
00073 
00074 # Minimal selection for all electrons, also basis for signal and veto muons
00075 electronCut  =     'pt > 20.'                                                                              # transverse energy
00076 electronCut += ' && abs(eta) < 2.5'                                                                        # pseudo-rapisity range
00077 electronCut += ' && electronID("mvaTrigV0") > 0.'                                                          # MVA electrons ID
00078 electronCut += ' && (chargedHadronIso+max(0.,neutralHadronIso)+photonIso-0.5*puChargedHadronIso)/et < 0.2' # relative isolation with Delta beta corrections
00079 
00080 ### ------------------------------------------------------------------------ ###
00081 
00082 
00083 ### Trigger matching
00084 
00085 # Trigger object selection
00086 triggerObjectSelectionDataRelVals = 'type("TriggerMuon") && ( path("HLT_IsoMu17_eta2p1_TriCentralPFJet30_v*") )' # 2012A RelVals
00087 triggerObjectSelectionData        = 'type("TriggerMuon") && ( path("HLT_IsoMu17_eta2p1_TriCentralPFNoPUJet30_30_20_v*") )'
00088 triggerObjectSelectionMC          = 'type("TriggerMuon") && ( path("HLT_*") )' # not recommended