Use of the AT Gene f h T G eric Matching in PAT M h P T Roger - - PowerPoint PPT Presentation
Use of the AT Gene f h T G eric Matching in PAT M h P T Roger Wolf Jan V Roger Wolf, Jan V Ve erka Veverka PAT Tutorial Ju une 2010 Object Matching Object Matching What ever Analysis you do you will arrive e at some point where you
Matching for tag & probe methods (STA Muon to Glo Matching of calibration truth and calibration measure
Matching of calibration truth and calibration measure Matching of trigger objects to reco objects Matching of MC truth vs reco object
ement (Photon to Jet) ement (Photon to Jet)
Five-fold template, which can be expande
Input Types Match Config.
ed using several methods/classes d object pre-selection.
Special Selection criteria for MC particles for match collection Allow window in DR and Dpt for matching Match by distance in DR.(closest in DR)
jetPartonMatch = cms.EDFilter("MCMatcher", # cut on src = cms.InputTag("ak5CaloJets"), # RECO o matched = cms.InputTag("genParticles"), # mc-tru mcPdgId = cms.vint32(1, 2, 3, 4, 5, 21), # one or i 32(3) # mcStatus = cms.vint32(3), # PYTHIA checkCharge = cms.bool(False), # False maxDeltaR = cms.double(0.4), # Minimu maxDPtRel = cms.double(3.0), # Minimu resolveAmbiguities = cms.bool(True), # Forbid resolveByMatchQuality = cms.bool(False), # False )
O j C Color Coding: PhysObjectMatcher MatchByDRDPt MCMat
1)to follow the expansion in detail have a look into the links
n deltaR, deltaPt/Pt; pick best by deltaR
uth particle collection r more PDG ID (quarks except top; gluons) d (3 h d i ) A status code (3 = hard scattering) = any value of the charge of MC and RECO is ok um deltaR for the match um deltaPt/Pt for the match d two RECO objects to match to the same GEN object = just match input in order; True = pick lowest dR pair first
S
tchSelector
jetPartonMatch = cms.EDProducer("MCMatcher", # cut src = cms.InputTag("ak5CaloJets"), # RECO matched = cms.InputTag("genParticles"), # mc-t mcPdgId = cms.vint32(1, 2, 3, 4, 5, 21), # one o i 32(3) # mcStatus = cms.vint32(3), # PYTH checkCharge = cms.bool(False), # Fals maxDeltaR = cms.double(0.4), # Mini maxDPtRel = cms.double(3.0), # Mini resolveAmbiguities = cms.bool(True), # Forb resolveByMatchQuality = cms.bool(False), # Fals ) jetGenJetMatch = cms.EDProducer("GenJetMatcher", # cut src = cms.InputTag("ak5CaloJets"), # RECO matched = cms InputTag("ak5GenJets") # GEN matched = cms.InputTag("ak5GenJets"), # GEN mcPdgId = cms.vint32(), # n/a mcStatus = cms.vint32(), # n/a checkCharge = cms.bool(False), # n/a maxDeltaR = cms.double(0.4), # Mini maxDPtRel = cms.double(3.0), # Mini resolveAmbiguities = cms.bool(True), # Forb resolveByMatchQuality = cms.bool(False), # Fals )
O j C Color Coding: PhysObjectMatcher MatchByDRDPt MCMat
1)to follow the expansion in detail have a look into the links
O objects to match truth particle collection
d (3 h d i ) HIA status code (3 = hard scattering) se = any value of the charge of MC and RECO is ok imum deltaR for the match imum deltaPt/Pt for the match bid two RECO objects to match to the same GEN object se = just match input in order; True = pick lowest dR pair first
O jets (any View<Jet> is ok) jets (must be GenJetCollection) jets (must be GenJetCollection) imum deltaR for the match imum deltaPt/Pt for the match bid two RECO objects to match to the same GEN object se = just match input in order; True = pick lowest dR pair first
S tchSelector
Find the muonMatch cfi.py in mcMatchingLayer0
_ py g y Find the muonProducer_cff.py in producersLayer1
Find the jetMatch_cfi.py in mcMatchingLayer0 Find the jetProducer_cff.py in producersLayer1
patMuons = cms.EDProducer("PATMuonProducer", # input muonSource = cms.InputTag("muons"), # use particle flow instead of std reco useParticleFlow = cms.bool( False ), pfMuonSource = cms.InputTag("pfIsolatedMuons"), # add TeV refit tracks ... # add user data … # embedding objects ... # embed IsoDeposits isoDeposits = cms.PSet(), # user defined isolation variables the variables d # via pat::Muon::userIsolation(IsolationKeys key) # DataFormats/PatCandidates/interface/Isolation.h userIsolation = cms.PSet(), # mc matching addGenMatch = cms.bool(True), embedGenMatch = cms.bool(True), genParticleMatch = cms.InputTag("muonMatch"), ## p # efficiencies # efficiencies addEfficiencies = cms.bool(False), efficiencies = cms.PSet(), # resolution configurables addResolutions = cms.bool(False), resolutions = cms.PSet(), # high level selections b d i h l l i b l( ) embedHighLevelSelection = cms.bool(True), beamLineSrc = cms.InputTag("offlineBea )
, defined here will be accessible with the key as defined in particles source to be used for the matching amSpot")
patJets = cms.EDProducer("PATJetProducer", # input jetSource = cms.InputTag("ak5CaloJets"), # add user data ... # embedding of AOD items embedCaloTowers = cms.bool(False), ## ... # jet energy corrections addJetCorrFactors = cms.bool(True), jetCorrFactorsSource = cms.VInputTag(cms.InputTag( # btag information # btag information ... # clone tag infos ATTENTION: these take lots of sp # usually the discriminators from the default algo # are sufficient ... # mc matching addGenPartonMatch = cms.bool(True), embedGenPartonMatch = cms.bool(True), genPartonMatch = cms.InputTag("patJetPartonMa addGenJetMatch = cms.bool(True), embedGenJetMatch = cms.bool(True), genJetMatch = cms.InputTag("patJetGenJetMa ... ... )
("jetCorrFactors") ), pace!
## ... ## ... atch"), ## ... ## ... ## ... atch"), ## ...