Project CMSSW displayed by LXR CMSSW_14_2_X_2024-09-06-2300/​PhysicsTools/​PatExamples/​test/​analyzePatMCMatchingExtended_cfg.py	Source navigation Diff markup Identifier search General search
	Version: CMSSW_14_2_X_2024-09-06-2300 CMSSW_14_2_X_2024-09-05-2300 CMSSW_14_2_X_2024-09-04-2300 CMSSW_14_2_X_2024-09-03-2300 CMSSW_14_2_X_2024-09-01-2300 Revert   Change

File indexing completed on 2024-04-06 12:24:02

 import FWCore.ParameterSet.Config as cms
 
 process = cms.Process("Test")
 
 #------------------------------------------------------------------------------------------------------
 # To configure the Matching, we have to configure the PAT-Workflow starting from the patDefaultSequence:
 #------------------------------------------------------------------------------------------------------
 
 from PhysicsTools.PatAlgos.patTemplate_cfg import *
 
 ## prepare several clones of match associations for status 1, 3 and in flight muons (status -1)
 ## add inFlightMuons
 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
 process.inFlightMuons = cms.EDProducer("PATGenCandsFromSimTracksProducer",
         src           = cms.InputTag("g4SimHits"),   ## use "fastSimProducer" for FastSim
         setStatus     = cms.int32(-1),
         particleTypes = cms.vstring("mu+"),          ## picks also mu-, of course
         filter        = cms.vstring("pt > 0.5"),     ## just for testing
         makeMotherLink = cms.bool(True),
         writeAncestors = cms.bool(True),             ## save also the intermediate GEANT ancestors of the muons
         genParticles   = cms.InputTag("genParticles"),
 )
 
 process.muMatch3 = process.muonMatch.clone(mcStatus = [3]) # hard scattering
 process.muMatch1 = process.muonMatch.clone(mcStatus = [1]) # stable
 process.muMatchF = process.muonMatch.clone(mcStatus = [-1], matched = cms.InputTag("inFlightMuons"))
 
 process.muMatch1.checkCharge = False
 process.muMatch3.checkCharge = False
 
 #process.muMatch3.resolveByMatchQuality = True
 #process.muMatch1.resolveByMatchQuality = True
 
 process.muMatch3.maxDeltaR = 0.05
 process.muMatch3.maxDPtRel = 0.1
 
 process.muMatch1.maxDeltaR = 0.05
 process.muMatch1.maxDPtRel = 0.1
 
 process.muMatchF.maxDeltaR = 0.3
 process.muMatchF.maxDPtRel = 0.2
 
 
 process.muonMatchByPt = cms.EDProducer("MCMatcherByPt", # cut on deltaR, deltaPt/Pt; pick best by deltaPt
     src     = cms.InputTag("muons"), # RECO objects to match
     matched = cms.InputTag("genParticles"),   # mc-truth particle collection
     mcPdgId     = cms.vint32(13), # one or more PDG ID (13 = muon); absolute values (see below)
     checkCharge = cms.bool(True), # True = require RECO and MC objects to have the same charge
     mcStatus = cms.vint32(1),     # PYTHIA status code (1 = stable, 2 = shower, 3 = hard scattering)
     maxDeltaR = cms.double(0.5),  # Minimum deltaR for the match
     maxDPtRel = cms.double(0.5),  # Minimum deltaPt/Pt for the match
     resolveAmbiguities = cms.bool(True),     # Forbid two RECO objects to match to the same GEN object
     resolveByMatchQuality = cms.bool(False), # False = just match input in order; True = pick lowest deltaR pair first
 )
 
 
 ## add the new matches to the default sequence
 process.patDefaultSequence.replace(process.muonMatch,
                                    process.muMatch1 +
                                    process.muMatch3 +
                                    process.muMatchF
                                   #+ process.muonMatchByPt
 )
 
 process.patMuons.genParticleMatch = cms.VInputTag(
     cms.InputTag("muMatch3"),
     cms.InputTag("muMatch1"),
     cms.InputTag("muMatchF")
     #, cms.InputTag("muonMatchByPt")
 )
 
 
 #-----------------------------------------
 # As usual add those two usefull things:
 #----------------------------------------
 
 process.TFileService = cms.Service("TFileService",
   fileName = cms.string('analyzePatMCMatchingExtended.root')
 )
 
 process.MessageLogger = cms.Service("MessageLogger")
 
 
 #----------------------------------------------------------------------
 # Finally let's analyze the matching and run all that in correct order:
 #----------------------------------------------------------------------
 
 process.analyzePatMCMatching = cms.EDAnalyzer("PatMCMatchingExtended",
   muonSrc     = cms.untracked.InputTag("cleanPatMuons")
 )
 
 process.out.outputCommands = cms.untracked.vstring('keep *')
 process.outpath.remove(process.out)
 
 
 process.p = cms.Path(process.inFlightMuons + process.patDefaultSequence + process.analyzePatMCMatching)
 
 
 #----------------------------------------------------------------------
 # Change the input file to compare
 #----------------------------------------------------------------------
 process.maxEvents.input = -1
 
 ## Source
 from PhysicsTools.PatAlgos.tools.cmsswVersionTools import pickRelValInputFiles
 process.source = cms.Source("PoolSource",
     fileNames = cms.untracked.vstring(
      pickRelValInputFiles( cmsswVersion  = 'CMSSW_4_2_8'
                         #, relVal        =  'RelValZMM'
                         , relVal        =  'RelValJpsiMM'
                         , globalTag     = 'START42_V12'
                         , numberOfFiles = 1
                         )
     )
 )
 
 del(process.out)
 del(process.outpath)

This page was automatically generated by the 2.2.1 LXR engine. The LXR team