def Pat(self) :
from PhysicsTools.PatAlgos.tools.coreTools import removeMCMatching,removeCleaning,restrictInputToAOD,removeAllPATObjectsBut
from PhysicsTools.PatAlgos.tools.jetTools import switchJetCollection
from PhysicsTools.PatAlgos.tools.jetTools import addJetCollection
if self.options.isData: removeMCMatching(self.process, ['All'])
restrictInputToAOD(self.process)
removeCleaning(self.process)
removeAllPATObjectsBut(self.process, ['Jets','METs','Muons'])
addJetCollection(self.process,
cms.InputTag('ak5CaloJets'),
'AK5',
'Calo',
jetCorrLabel = ('AK5Calo',self.options.jetCorrections),
genJetCollection = cms.InputTag('ak5GenJets'),
doType1MET = False,
)
switchJetCollection(self.process,
cms.InputTag('ak5PFJets'),
doJTA = True,
doBTagging = True,
doJetID = True,
jetCorrLabel = ('AK5PF',self.options.jetCorrections),
genJetCollection = cms.InputTag('ak5GenJets'),
doType1MET = False,
jetIdLabel = 'ak5pf'
)
self.process.selectedPatJets.cut = cms.string("pt > 30 && \
abs(eta) < 3.0 && \
neutralHadronEnergyFraction < 0.9 && \
neutralEmEnergyFraction < 0.90 && \
nConstituents > 1 && \
(? abs(eta)<2.4 ? chargedHadronEnergyFraction : 1) > 0 && \
(? abs(eta)<2.4 ? chargedHadronMultiplicity : 1) > 0 && \
(? abs(eta)<2.4 ? chargedEmEnergyFraction : 0) < 0.99")
del self.process.out
del self.process.outpath
return self.process.patDefaultSequence
import PhysicsTools.PatAlgos.tools.coreTools as patcore
import HtoZg.CommonAnalysis.process_cfi
process = copy.deepcopy(HtoZg.CommonAnalysis.process_cfi.process)
## Input files
process.source.fileNames.append(
'/store/mc/Summer12/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball/AODSIM/PU_S7_START52_V9-v2/0003/EAF43999-8D9B-E111-A418-003048D4610E.root'
)
## Maximal Number of Events
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(-1) )
## Standard PAT Configuration File
# process.load("PhysicsTools.PatAlgos.patSequences_cff")
## The HtoZg sequence
process.load('HtoZg.MuonAnalysis.skimSequence_cff')
patcore.removeAllPATObjectsBut(process, ['Muons', 'Photons'])
## TFileService for the ntuple output
process.load('HtoZg.CommonAnalysis.TFileService_cfi')
process.TFileService.fileName = 'mc.root'
process.p = cms.Path(process.skimSequence)
if __name__ == '__main__':
## Adds tab-completion and history for interactive testing.
import user
## Input files
process.source.fileNames.append(
"/store/mc/Summer12/DYJetsToLL_M-50_TuneZ2Star_8TeV-madgraph-tarball/AODSIM/PU_S7_START52_V9-v2/0003/EAF43999-8D9B-E111-A418-003048D4610E.root"
)
process.source.eventsToProcess = cms.untracked.VEventRange("1:10847674", "1:14658698")
## Maximal Number of Events
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(-1))
## The full skim sequence
process.load("HtoZg.MuonAnalysis.skimSequence_cff")
## PATification
# process.load('HtoZg.MuonAnalysis.patSequence_cff')
patcore.removeAllPATObjectsBut(process, ["Muons"])
## The muon tree maker
process.load("HtoZg.MuonAnalysis.muonTree_cfi")
## TFileService for the ntuple output
process.load("HtoZg.CommonAnalysis.TFileService_cfi")
process.TFileService.fileName = "zg_pick.root"
process.p = cms.Path(
# process.filterSequence +
process.patSequence
+ process.muonTree
)
if __name__ == "__main__":