process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff") process.GlobalTag.globaltag = cms.string(dataVersion.getGlobalTag()) print "GlobalTag="+dataVersion.getGlobalTag() process.load("HiggsAnalysis.HeavyChHiggsToTauNu.HChCommon_cfi") # Uncomment the following in order to print the counters at the end of # the job (note that if many other modules are being run in the same # job, their INFO messages are printed too) #process.MessageLogger.cerr.threshold = cms.untracked.string("INFO") # Fragment to run PAT on the fly if requested from command line options.doPat = 1 options.trigger="HLT_SingleIsoTau20_Trk15_MET20" from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly process.commonSequence, additionalCounters = addPatOnTheFly(process, options, dataVersion) # Add configuration information to histograms.root #from HiggsAnalysis.HeavyChHiggsToTauNu.HChTools import addConfigInfo #process.infoPath = addConfigInfo(process, options, dataVersion) ################################################################################ # The "golden" version of the signal analysis #import HiggsAnalysis.HeavyChHiggsToTauNu.HChSignalAnalysisParameters_cff as param #param.overrideTriggerFromOptions(options) # Set tau selection mode to 'standard' or 'factorized' #param.setAllTauSelectionOperatingMode('standard') #param.setAllTauSelectionOperatingMode('factorized')
fileNames = cms.untracked.vstring( "/store/group/local/HiggsChToTauNuFullyHadronic/pattuples/CMSSW_4_4_X/TTJets_TuneZ2_Fall11/TTJets_TuneZ2_7TeV-madgraph-tauola/Fall11_PU_S6_START44_V9B_v1_AODSIM_tauembedding_gentauskim_v44_5/9ecb3a23e436fc2ffd8a803eac2a3a15/pattuple_1012_1_LSv.root", ), ) process.load("HiggsAnalysis.HeavyChHiggsToTauNu.HChCommon_cfi") # Fragment to run PAT on the fly if requested from command line from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly patArgs = {"doPatTrigger": False, # "doPatTaus": False, # "doHChTauDiscriminators": False, "doPatElectronID": True, "doTauHLTMatching": False, } process.commonSequence, additionalCounters = addPatOnTheFly(process, options, dataVersion, patArgs=patArgs, doHBHENoiseFilter=False, ) import HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration as AnalysisConfiguration dataEras = [ "Run2011AB", "Run2011A", "Run2011B", ] puWeights = AnalysisConfiguration.addPuWeightProducers(dataVersion, process, process.commonSequence, dataEras) # Add GenTau skim counters import HiggsAnalysis.HeavyChHiggsToTauNu.CustomGenTauSkim as tauSkim additionalCounters = tauSkim.getCounters() + additionalCounters # Add configuration information to histograms.root import HiggsAnalysis.HeavyChHiggsToTauNu.HChTools as HChTools
# the job (note that if many other modules are being run in the same # job, their INFO messages are printed too) #process.MessageLogger.cerr.threshold = cms.untracked.string("INFO") # Fragment to run PAT on the fly if requested from command line from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly from PhysicsTools.PatAlgos.tools.coreTools import removeSpecificPATObjects patArgs = { "doPatTrigger": False, # "doPatTaus": False, # "doHChTauDiscriminators": False, "doPatElectronID": True, "doTauHLTMatching": False, "doPatMuonPFIsolation": True, } process.commonSequence, additionalCounters = addPatOnTheFly( process, options, dataVersion, plainPatArgs=patArgs) #process.commonSequence.remove(process.goodPrimaryVertices10) if options.doPat == 0: process.load("HiggsAnalysis.HeavyChHiggsToTauNu.HChPrimaryVertex_cfi") process.commonSequence *= (process.goodPrimaryVertices * process.goodPrimaryVertices10) from HiggsAnalysis.HeavyChHiggsToTauNu.HChTools import * # Pileup weighting if dataVersion.isMC(): import HiggsAnalysis.HeavyChHiggsToTauNu.HChSignalAnalysisParameters_cff as param # Pileup weighting process.pileupWeight = cms.EDProducer( "HPlusVertexWeightProducer", alias=cms.string("pileupWeight"),
fileNames = cms.untracked.vstring( "/store/group/local/HiggsChToTauNuFullyHadronic/pattuples/CMSSW_4_4_X/TTJets_TuneZ2_Fall11/TTJets_TuneZ2_7TeV-madgraph-tauola/Fall11_PU_S6_START44_V9B_v1_AODSIM_tauembedding_gentauskim_v44_5/9ecb3a23e436fc2ffd8a803eac2a3a15/pattuple_1012_1_LSv.root", ), ) process.load("HiggsAnalysis.HeavyChHiggsToTauNu.HChCommon_cfi") # Fragment to run PAT on the fly if requested from command line from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly patArgs = {"doPatTrigger": False, # "doPatTaus": False, # "doHChTauDiscriminators": False, "doPatElectronID": True, "doTauHLTMatching": False, } process.commonSequence, additionalCounters = addPatOnTheFly(process, options, dataVersion, patArgs=patArgs, doHBHENoiseFilter=False, ) import HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration as AnalysisConfiguration dataEras = [ "Run2011AB", "Run2011A", "Run2011B", ] puWeights = AnalysisConfiguration.addPuWeightProducers(dataVersion, process, process.commonSequence, dataEras) # Add GenTau skim counters import HiggsAnalysis.HeavyChHiggsToTauNu.CustomGenTauSkim as tauSkim additionalCounters = tauSkim.getCounters() + additionalCounters import HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.customisations as customisations process.genTaus = cms.EDFilter("GenParticleSelector", src = cms.InputTag("genParticles"),
fileName=cms.untracked.string('dummy.root'), outputCommands=cms.untracked.vstring()) from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly patArgs = { # "doPatTaus": False, "doPatMET": False, "doPatElectronID": False, "doPatCalo": False, "doBTagging": False, "doPatMuonPFIsolation": True, "doTauHLTMatching": False, } process.commonSequence, counters = addPatOnTheFly(process, options, dataVersion, plainPatArgs=patArgs, doMcPreselection=False) del process.out process.patDefaultSequence.remove(process.countPatTaus) # Triggering if dataVersion.isMC(): process.load("HLTrigger.HLTfilters.triggerResultsFilter_cfi") process.triggerResultsFilter.hltResults = cms.InputTag( "TriggerResults", "", dataVersion.getTriggerProcess()) process.triggerResultsFilter.l1tResults = cms.InputTag("") # dummy process.triggerResultsFilter.throw = cms.bool(True) process.triggerResultsFilter.triggerConditions = cms.vstring(trigger) process.commonSequence *= process.triggerResultsFilter else:
process.out = cms.OutputModule("PoolOutputModule", fileName = cms.untracked.string('dummy.root'), outputCommands = cms.untracked.vstring() ) from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly patArgs = { # "doPatTaus": False, "doPatMET": False, "doPatElectronID": False, "doPatCalo": False, "doBTagging": False, "doPatMuonPFIsolation": True, "doTauHLTMatching": False, } process.commonSequence, counters = addPatOnTheFly(process, options, dataVersion, plainPatArgs=patArgs, doMcPreselection=False) del process.out process.patDefaultSequence.remove(process.countPatTaus) # Triggering if dataVersion.isMC(): process.load("HLTrigger.HLTfilters.triggerResultsFilter_cfi") process.triggerResultsFilter.hltResults = cms.InputTag("TriggerResults", "", dataVersion.getTriggerProcess()) process.triggerResultsFilter.l1tResults = cms.InputTag("") # dummy process.triggerResultsFilter.throw = cms.bool(True) process.triggerResultsFilter.triggerConditions = cms.vstring(trigger) process.commonSequence *= process.triggerResultsFilter else: process.TriggerFilter.triggerConditions = [trigger] process.triggeredCount = cms.EDProducer("EventCountProducer")
from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly patArgs = { # "doPatTaus": False, # "doPatMET": False, "doPatElectronID": False, # "doPatCalo": False, # "doBTagging": False, "doTauHLTMatching": False, } process.out = cms.OutputModule( "PoolOutputModule", fileName=cms.untracked.string('dummy.root'), outputCommands=cms.untracked.vstring(), ) process.commonSequence, counters = addPatOnTheFly(process, options, dataVersion, patArgs=patArgs) del process.out # Pileup weighting from HiggsAnalysis.HeavyChHiggsToTauNu.HChTools import * import HiggsAnalysis.HeavyChHiggsToTauNu.HChSignalAnalysisParameters_cff as param process.pileupWeight = cms.EDProducer( "HPlusVertexWeightProducer", alias=cms.string("pileupWeight"), ) puweight = "Run2011AB" if len(options.puWeightEra) > 0: puweight = options.puWeightEra param.setPileupWeight(dataVersion,
src = cms.InputTag("genMuons"), minNumber = cms.uint32(1), ) process.genMuonsCount = cms.EDProducer("EventCountProducer") process.commonSequence += ( process.allEvents + process.genMuons + process.genMuonsFilter + process.genMuonsCount ) additionalCounters.extend(["allEvents", "genMuonsCount"]) # PAT muons options.doPat = 1 from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly process.fooSequenceSequence, fooCounters = addPatOnTheFly(process, options, dataVersion) process.patMuons.userData.userFloats.src = [] process.selectedPatMuons.src = "patMuons" process.selectedPatMuons.cut = "" process.commonSequence += ( process.pfParticleSelectionSequence + process.muIsoSequence + process.makePatMuons + process.selectedPatMuons ) # MuScleFit correction # https://twiki.cern.ch/twiki/bin/view/CMSPublic/MuScleFitCorrections2012 process.muscleMuons = cms.EDProducer("MuScleFitPATMuonCorrector", src = cms.InputTag("selectedPatMuons"), debug = cms.bool(False),
process.out = cms.OutputModule("PoolOutputModule", SelectEvents = cms.untracked.PSet( SelectEvents = cms.vstring("path") ), fileName = cms.untracked.string('skim.root'), outputCommands = cms.untracked.vstring() ) process.selectionSequence = cms.Sequence() from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly patArgs = {"doTauHLTMatching": False, "doPatTaus": True, "doPatMuonPFIsolation": True, } process.commonSequence, additionalCounters = addPatOnTheFly(process, options, dataVersion, plainPatArgs=patArgs, doPlainPat=True, doMcPreselection=True) # In order to avoid transient references and generalTracks is available anyway if hasattr(process, "patMuons"): process.patMuons.embedTrack = False # Override the outputCommands here, since PAT modifies it process.load("HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.HChEventContent_cff") process.out.outputCommands = cms.untracked.vstring( "keep *", "drop *_MEtoEDMConverter_*_*", # drop DQM histos "drop *_*_*_MUONSKIM", ) import re name_re = re.compile("_\*$") process.out.outputCommands.extend([name_re.sub("_MUONSKIM", x) for x in process.HChEventContent.outputCommands])
# PAT, if requested from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly patArgs = { # "doPatTaus": False, # "doPatMET": False, "doPatElectronID": False, # "doPatCalo": False, # "doBTagging": False, "doTauHLTMatching": False, } process.out = cms.OutputModule("PoolOutputModule", fileName = cms.untracked.string('dummy.root'), outputCommands = cms.untracked.vstring(), ) process.commonSequence, counters = addPatOnTheFly(process, options, dataVersion, patArgs=patArgs) del process.out # Pileup weighting from HiggsAnalysis.HeavyChHiggsToTauNu.HChTools import * import HiggsAnalysis.HeavyChHiggsToTauNu.HChSignalAnalysisParameters_cff as param process.pileupWeight = cms.EDProducer("HPlusVertexWeightProducer", alias = cms.string("pileupWeight"), ) puweight = "Run2011AB" if len(options.puWeightEra) > 0: puweight = options.puWeightEra param.setPileupWeight(dataVersion, process=process, commonSequence=process.commonSequence, era=puweight) insertPSetContentsTo(param.vertexWeight, process.pileupWeight) process.pileupWeight.vertexSrc = "offlinePrimaryVertices"
SelectEvents=cms.untracked.PSet(SelectEvents=cms.vstring("path")), fileName=cms.untracked.string('skim.root'), outputCommands=cms.untracked.vstring()) process.selectionSequence = cms.Sequence() from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly patArgs = { "doTauHLTMatching": False, "doPatTaus": True, "doPatMuonPFIsolation": True, } process.commonSequence, additionalCounters = addPatOnTheFly( process, options, dataVersion, plainPatArgs=patArgs, doPlainPat=True, doMcPreselection=True) # In order to avoid transient references and generalTracks is available anyway if hasattr(process, "patMuons"): process.patMuons.embedTrack = False # Override the outputCommands here, since PAT modifies it process.load( "HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.HChEventContent_cff") process.out.outputCommands = cms.untracked.vstring( "keep *", "drop *_MEtoEDMConverter_*_*", # drop DQM histos "drop *_*_*_MUONSKIM", )
"abs(pdgId()) == 13 && pt() > 40 && abs(eta()) < 2.1 && abs(mother().pdgId()) != 13" )) process.genMuonsFilter = cms.EDFilter( "CandViewCountFilter", src=cms.InputTag("genMuons"), minNumber=cms.uint32(1), ) process.genMuonsCount = cms.EDProducer("EventCountProducer") process.commonSequence += (process.allEvents + process.genMuons + process.genMuonsFilter + process.genMuonsCount) additionalCounters.extend(["allEvents", "genMuonsCount"]) # PAT muons options.doPat = 1 from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly process.fooSequenceSequence, fooCounters = addPatOnTheFly( process, options, dataVersion) process.patMuons.userData.userFloats.src = [] process.selectedPatMuons.src = "patMuons" process.selectedPatMuons.cut = "" process.commonSequence += (process.pfParticleSelectionSequence + process.muIsoSequence + process.makePatMuons + process.selectedPatMuons) # MuScleFit correction # https://twiki.cern.ch/twiki/bin/view/CMSPublic/MuScleFitCorrections2012 process.muscleMuons = cms.EDProducer( "MuScleFitPATMuonCorrector", src=cms.InputTag("selectedPatMuons"), debug=cms.bool(False), identifier=cms.string("Fall11_START44"), applySmearing=cms.bool(False),
outputCommands = cms.untracked.vstring( "keep *", "drop *_MEtoEDMConverter_*_*", # drop DQM histos "drop *_*_*_MUONSKIM", "keep *_tightMuons*_*_MUONSKIM", "keep *_goodJets*_*_MUONSKIM", ) ) process.selectionSequence = cms.Sequence() from HiggsAnalysis.HeavyChHiggsToTauNu.HChPatTuple import addPatOnTheFly patArgs = {"doTauHLTMatching": False, } process.commonSequence, additionalCounters = addPatOnTheFly(process, options, dataVersion, patArgs=patArgs, doHBHENoiseFilter=False, # Only save the HBHE result to event, don't filter calculateEventCleaning=True, # This requires the tags from test/pattuple/checkoutTags.sh ) # In order to avoid transient references and generalTracks is available anyway #if hasattr(process, "patMuons"): # process.patMuons.embedTrack = False # process.patMuonsPFlow.embedTrack = False # process.patMuonsPFlowChs.embedTrack = False # Override the outputCommands here, since PAT modifies it # process.load("HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.HChEventContent_cff") # process.out.outputCommands = cms.untracked.vstring( # ) # import re # name_re = re.compile("_\*$") # process.out.outputCommands.extend([name_re.sub("_MUONSKIM", x) for x in process.HChEventContent.outputCommands])