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__":