Reference Manual

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_TriCentralJet30_v*' # 2011B RelVals
00030 triggerSelectionData        = 'HLT_Iso10Mu20_eta2p1_TriCentralPFJet30_v*'
00031 triggerSelectionMC          = 'HLT_IsoMu20_eta2p1_TriCentralPFJet30_v*'
00032 
00033 ### Muon selection
00034 
00035 # Minimal selection for all muons, also basis for signal and veto muons
00036 muonCutBase  =     'pt > 10.'                                                    # transverse momentum
00037 muonCutBase += ' && abs(eta) < 2.5'                                              # pseudo-rapisity range
00038 # veto muon
00039 muonCutVeto  = muonCutBase
00040 muonCutVeto += ' && isGlobalMuon'                                                # general reconstruction property
00041 # standard muon
00042 muonCut  = muonCutVeto
00043 muonCut += ' && (trackIso+caloIso)/pt < 0.2'                                     # relative isolation
00044 # PF muon
00045 muonCutPF  =  muonCutVeto
00046 muonCutPF += ' && (chargedHadronIso+neutralHadronIso+photonIso)/pt < 0.2'        # relative isolation
00047 
00048 # Signal muon selection on top of 'muonCut'
00049 looseMuonCutBase  =     'isTrackerMuon'                                          # general reconstruction property
00050 looseMuonCutBase += ' && pt > 20.'                                               # transverse momentum
00051 looseMuonCutBase += ' && abs(eta) < 2.1'                                         # pseudo-rapisity range
00052 looseMuonCutBase += ' && globalTrack.normalizedChi2 < 10.'                       # muon ID: 'isGlobalMuonPromptTight'
00053 looseMuonCutBase += ' && globalTrack.hitPattern.numberOfValidMuonHits > 0'       # muon ID: 'isGlobalMuonPromptTight'
00054 looseMuonCutBase += ' && abs(dB) < 0.02'                                         # 2-dim impact parameter with respect to beam spot (s. "PAT muon configuration" above)
00055 looseMuonCutBase += ' && innerTrack.numberOfValidHits > 10'                      # tracker reconstruction
00056 looseMuonCutBase += ' && innerTrack.hitPattern.pixelLayersWithMeasurement >= 1'  # tracker reconstruction
00057 looseMuonCutBase += ' && numberOfMatches > 1'                                    # muon chamber reconstruction
00058 #looseMuonCut += ' && ...'                                                        # DeltaZ between muon vertex and PV < 1. not accessible via configuration yet
00059 # standard muon
00060 looseMuonCut  = looseMuonCutBase
00061 looseMuonCut += ' && (trackIso+caloIso)/pt < 0.1'                                # relative isolation
00062 # PF muon
00063 looseMuonCutPF = looseMuonCutBase
00064 looseMuonCutPF += ' && (chargedHadronIso+neutralHadronIso+photonIso)/pt < 0.125' # relative isolation
00065 # Signal muon distance from signal jet
00066 muonJetsDR = 0.3                                                                 # minimal DeltaR of signal muons from any signal jet
00067 
00068 # Tightened signal muon selection on top of 'looseMuonCut'
00069 tightMuonCutBase  =     ''
00070 # standard muon
00071 tightMuonCut  = tightMuonCutBase
00072 tightMuonCut += '(trackIso+caloIso)/pt < 0.05'                                   # relative isolation
00073 # PF muon
00074 tightMuonCutPF  = tightMuonCutBase
00075 tightMuonCutPF += '(chargedHadronIso+neutralHadronIso+photonIso)/pt < 0.125'     # relative isolation
00076 
00077 ### Jet selection
00078 
00079 # Signal jet selection
00080 jetCutBase  =     'pt > 30.'                                             # transverse momentum
00081 jetCutBase += ' && abs(eta) < 2.4'                                       # pseudo-rapisity range
00082 # standard jet
00083 jetCut  = jetCutBase
00084 jetCut += ' && emEnergyFraction > 0.01'                                  # jet ID: electro-magnetic energy fraction
00085 jetCut += ' && jetID.n90Hits > 1'                                        # jet ID: number of RecHits carying 90% of the total energy
00086 jetCut += ' && jetID.fHPD < 0.98'                                        # jet ID: fraction of energy in the hottest readout
00087 # PF jet
00088 jetCutPF  = jetCutBase
00089 jetCutPF += ' && numberOfDaughters > 1'                                  # PF jet ID:
00090 jetCutPF += ' && chargedEmEnergyFraction < 0.99'                         # PF jet ID:
00091 jetCutPF += ' && neutralHadronEnergyFraction < 0.99'                     # PF jet ID:
00092 jetCutPF += ' && neutralEmEnergyFraction < 0.99'                         # PF jet ID:
00093 jetCutPF += ' && (chargedHadronEnergyFraction > 0. || abs(eta) >= 2.4)'  # PF jet ID:
00094 jetCutPF += ' && (chargedMultiplicity > 0 || abs(eta) >= 2.4)'           # PF jet ID:
00095 # Signal jet distance from signal muon
00096 jetMuonsDRPF = 0.1
00097 
00098 ### Electron selection
00099 
00100 # Minimal selection for all electrons, also basis for signal and veto muons
00101 electronCutBase  =     'et > 15.'                                                  # transverse energy
00102 electronCutBase += ' && abs(eta) < 2.5'                                            # pseudo-rapisity range
00103 # veto electron
00104 electronCutVeto = electronCutBase
00105 # standard electron
00106 electronCut  = electronCutVeto
00107 electronCut += ' && (dr03TkSumPt+dr03EcalRecHitSumEt+dr03HcalTowerSumEt)/et < 0.2' # relative isolation
00108 # PF electron
00109 electronCutPF  = electronCutVeto
00110 electronCutPF += ' && (chargedHadronIso+neutralHadronIso+photonIso)/et < 0.2'      # relative isolation
00111 
00112 ### ------------------------------------------------------------------------ ###
00113 
00114 
00115 ### Trigger matching
00116 
00117 # Trigger object selection
00118 triggerObjectSelectionDataRelVals = 'type("TriggerMuon") && ( path("HLT_IsoMu17_eta2p1_TriCentralJet30_v*") )' # 2011B RelVals
00119 triggerObjectSelectionData        = 'type("TriggerMuon") && ( path("HLT_Iso10Mu20_eta2p1_TriCentralPFJet30_v*") )'
00120 triggerObjectSelectionMC          = 'type("TriggerMuon") && ( path("HLT_IsoMu20_eta2p1_TriCentralPFJet30_v*") )'
/data/refman/pasoursint/CMSSW_5_2_9/src/TopQuarkAnalysis/Configuration/python/patRefSel_refMuJets.py
