import FWCore.ParameterSet.Config as cms
from HiggsAnalysis.HeavyChHiggsToTauNu.HChOptions import getOptions
from HiggsAnalysis.HeavyChHiggsToTauNu.HChDataVersion import DataVersion
import FWCore.ParameterSet.VarParsing as VarParsing
dataVersion = "36X"
#dataVersion = "36Xspring10"
#dataVersion = "37X"
#dataVersion = "38X"
#dataVersion = "data" # this is for collision data
options = getOptions()
if options.dataVersion != "":
dataVersion = options.dataVersion
options.doPat=1
print "Assuming data is ", dataVersion
dataVersion = DataVersion(dataVersion) # convert string to object
process = cms.Process("MUONNTUPLE")
#process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10000) )
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(10) )
process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
process.GlobalTag.globaltag = cms.string(dataVersion.getGlobalTag())
process.source = cms.Source('PoolSource',
fileNames = cms.untracked.vstring(
#dataVersion.getPatDefaultFileCastor()
#dataVersion.getPatDefaultFileMadhatter()
import FWCore.ParameterSet.Config as cms
import FWCore.ParameterSet.VarParsing as VarParsing
from HiggsAnalysis.HeavyChHiggsToTauNu.HChOptions import getOptions
from HiggsAnalysis.HeavyChHiggsToTauNu.HChDataVersion import DataVersion
options = VarParsing.VarParsing()
options.register('inputFiles', '', options.multiplicity.list,
options.varType.string, "Files to process")
options.register("hltProcess", "", options.multiplicity.singleton,
options.varType.string, "HLT Process name")
options.register("printPaths", 0, options.multiplicity.singleton,
options.varType.int, "Print HLT paths")
options = getOptions(options)
hltProcess = "HLT"
if options.dataVersion != "":
hltProcess = DataVersion(options.dataVersion).getTriggerProcess()
if options.hltProcess != "":
hltProcess = options.hltProcess
print "Using HLT process name", hltProcess
process = cms.Process("HLTTABLEANALYZER")
maxEv = -1
if options.printPaths != 0:
maxEv = 10
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(maxEv))
process.options = cms.untracked.PSet(
Rethrow=cms.untracked.vstring("ProductNotFound"))
options.register("hltProcess",
"",
options.multiplicity.singleton,
options.varType.string,
"HLT Process name")
options.register("printPaths",
0,
options.multiplicity.singleton,
options.varType.int,
"Print HLT paths")
options.register("globalTag",
"",
options.multiplicity.singleton,
options.varType.string,
"Specify GlobalTag")
options = getOptions(options)
hltProcess = "HLT"
if options.dataVersion != "":
hltProcess = DataVersion(options.dataVersion).getTriggerProcess()
if options.hltProcess != "":
hltProcess = options.hltProcess
print "Using HLT process name", hltProcess
process = cms.Process("HLTTABLEANALYZER")
maxEv = -1
if options.printPaths != 0:
maxEv = 10
process.maxEvents = cms.untracked.PSet( input = cms.untracked.int32(maxEv) )
process.options = cms.untracked.PSet(
def customise(process):
options = VarParsing.VarParsing()
options.register(
'overrideBeamSpot',
0, # default value, false
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.int,
"should I override beamspot in globaltag?")
options = getOptions(options)
# Muon isolation
import HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.customisations as customisations
process.muonIsolationSequence = cms.Sequence()
muons = customisations.addMuonIsolationEmbedding(
process,
process.muonIsolationSequence,
muons=process.tightenedMuons.src.value())
process.tightenedMuons.src = muons
process.ProductionFilterSequence.replace(
process.tightenedMuons,
process.muonIsolationSequence * process.tightenedMuons)
# output
outputModule = None
outdict = process.outputModules_()
if len(outdict) == 1:
outputModule = outdict.values()[0]
elif outdict.has_key("RECOSIMoutput"):
outputModule = outdict["RECOSIMoutput"]
#process.TFileService = cms.Service("TFileService", fileName = cms.string("histo_simulation.root") )
print "TAUOLA mdtau =", process.generator.ZTauTau.TauolaOptions.InputCards.mdtau
#process.source.duplicateCheckMode = cms.untracked.string('noDuplicateCheck')
processName = process.name_()
print "Adjusting event content to RAWSIM+misc"
process.load(
"HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.HChEventContent_cff")
#outputModule.outputCommands = cms.untracked.vstring("keep *")
outputModule.outputCommands = cms.untracked.vstring("drop *")
outputModule.outputCommands.extend(
process.RAWSIMEventContent.outputCommands)
#outputModule.outputCommands.extend(process.RAWEventContent.outputCommands)
outputModule.outputCommands.extend(process.AODEventContent.outputCommands)
#outputModule.outputCommands.extend(process.AODSIMEventContent.outputCommands)
outputModule.outputCommands.extend(process.HChEventContent.outputCommands)
outputModule.outputCommands.extend([
"keep *_generalTracks_*_*",
"keep *_adaptedMuonsFromWmunu_*_%s" % processName,
"keep *_dimuonsGlobal_*_%s" % processName,
"keep *_generator_*_%s" % processName
])
# Remove duplicate "drop *"
index = 0
for item in outputModule.outputCommands[:]:
if item == "drop *" and index > 0:
del outputModule.outputCommands[index]
index -= 1
index += 1
# Set the empty event filter source
process.filterEmptyEv.src.setProcessName(processName)
# Do we have to override the beam spot for data?
if options.overrideBeamSpot != 0:
bs = cms.string(
"BeamSpotObjects_2009_LumiBased_SigmaZ_v21_offline") # 42x data gt
#bs = cms.string("BeamSpotObjects_2009_LumiBased_SigmaZ_v18_offline") # 39x data gt
#bs = cms.string("BeamSpotObjects_2009_LumiBased_v17_offline") # 38x data gt
process.GlobalTag.toGet = cms.VPSet(
cms.PSet(record=cms.string("BeamSpotObjectsRcd"),
tag=bs,
connect=cms.untracked.string(
"frontier://FrontierProd/CMS_COND_31X_BEAMSPOT")))
print "BeamSpot in globaltag set to ", bs
else:
print "BeamSpot in globaltag not changed"
if hasattr(process, "hltTrigReport"):
process.hltTrigReport.HLTriggerResults.setProcessName(processName)
if hasattr(process, "DQM_FEDIntegrity_v2"):
process.schedule.remove(process.DQM_FEDIntegrity_v2)
if hasattr(process, "DQM_FEDIntegrity_v3"):
process.schedule.remove(process.DQM_FEDIntegrity_v3)
if hasattr(process, "HLTAnalyzerEndpath"):
process.schedule.remove(process.HLTAnalyzerEndpath)
del process.HLTAnalyzerEndpath
#process.load("HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.printGenParticles_cff")
#process.generation_step *= process.printGenParticles
# Disable gen vertex smearing
process.VtxSmeared = cms.EDProducer("FlatEvtVtxGenerator",
MaxZ=cms.double(0.0),
MaxX=cms.double(0.0),
MaxY=cms.double(0.0),
MinX=cms.double(0.0),
MinY=cms.double(0.0),
MinZ=cms.double(0.0),
TimeOffset=cms.double(0.0),
src=cms.InputTag("generator"))
return process
import FWCore.ParameterSet.Config as cms
from HiggsAnalysis.HeavyChHiggsToTauNu.HChOptions import getOptions
from HiggsAnalysis.HeavyChHiggsToTauNu.HChDataVersion import DataVersion
import FWCore.ParameterSet.VarParsing as VarParsing
#dataVersion = "36X"
#dataVersion = "36Xspring10"
#dataVersion = "37X"
#dataVersion = "38X"
dataVersion = "39Xredigi"
options = getOptions()
if options.dataVersion != "":
dataVersion = options.dataVersion
print "Assuming data is ", dataVersion
dataVersion = DataVersion(dataVersion) # convert string to object
process = cms.Process("PatValidation")
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(10))
process.source = cms.Source(
'PoolSource', fileNames=cms.untracked.vstring("file:pattuple.root"))
process.load("HiggsAnalysis.HeavyChHiggsToTauNu.HChCommon_cfi")
del process.TFileService
process.analyzer = cms.EDAnalyzer(
"HPlusPatTupleValidationAnalyzer",
electrons=cms.untracked.VPSet(
cms.PSet(src=cms.untracked.InputTag("selectedPatElectrons"),
def customise(process):
options = VarParsing.VarParsing()
options.register('overrideBeamSpot',
0, # default value, false
VarParsing.VarParsing.multiplicity.singleton,
VarParsing.VarParsing.varType.int,
"should I override beamspot in globaltag?")
options = getOptions(options)
# Muon isolation
import HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.customisations as customisations
process.muonIsolationSequence = cms.Sequence()
muons = customisations.addMuonIsolationEmbedding(process, process.muonIsolationSequence, muons=process.tightenedMuons.src.value())
process.tightenedMuons.src = muons
process.ProductionFilterSequence.replace(process.tightenedMuons, process.muonIsolationSequence*process.tightenedMuons)
# output
outputModule = None
outdict = process.outputModules_()
if len(outdict) == 1:
outputModule = outdict.values()[0]
elif outdict.has_key("RECOSIMoutput"):
outputModule = outdict["RECOSIMoutput"]
#process.TFileService = cms.Service("TFileService", fileName = cms.string("histo_simulation.root") )
print "TAUOLA mdtau =", process.generator.ZTauTau.TauolaOptions.InputCards.mdtau
#process.source.duplicateCheckMode = cms.untracked.string('noDuplicateCheck')
processName = process.name_()
print "Adjusting event content to RAWSIM+misc"
process.load("HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.HChEventContent_cff")
#outputModule.outputCommands = cms.untracked.vstring("keep *")
outputModule.outputCommands = cms.untracked.vstring("drop *")
outputModule.outputCommands.extend(process.RAWSIMEventContent.outputCommands)
#outputModule.outputCommands.extend(process.RAWEventContent.outputCommands)
outputModule.outputCommands.extend(process.AODEventContent.outputCommands)
#outputModule.outputCommands.extend(process.AODSIMEventContent.outputCommands)
outputModule.outputCommands.extend(process.HChEventContent.outputCommands)
outputModule.outputCommands.extend([
"keep *_generalTracks_*_*",
"keep *_adaptedMuonsFromWmunu_*_%s" % processName,
"keep *_dimuonsGlobal_*_%s" % processName,
"keep *_generator_*_%s" % processName
])
# Remove duplicate "drop *"
index = 0
for item in outputModule.outputCommands[:]:
if item == "drop *" and index > 0:
del outputModule.outputCommands[index]
index -= 1
index += 1
# Set the empty event filter source
process.filterEmptyEv.src.setProcessName(processName)
# Do we have to override the beam spot for data?
if options.overrideBeamSpot != 0:
bs = cms.string("BeamSpotObjects_2009_LumiBased_SigmaZ_v21_offline") # 42x data gt
#bs = cms.string("BeamSpotObjects_2009_LumiBased_SigmaZ_v18_offline") # 39x data gt
#bs = cms.string("BeamSpotObjects_2009_LumiBased_v17_offline") # 38x data gt
process.GlobalTag.toGet = cms.VPSet(
cms.PSet(record = cms.string("BeamSpotObjectsRcd"),
tag = bs,
connect = cms.untracked.string("frontier://FrontierProd/CMS_COND_31X_BEAMSPOT")
)
)
print "BeamSpot in globaltag set to ", bs
else:
print "BeamSpot in globaltag not changed"
if hasattr(process, "hltTrigReport"):
process.hltTrigReport.HLTriggerResults.setProcessName(processName)
if hasattr(process, "DQM_FEDIntegrity_v2"):
process.schedule.remove(process.DQM_FEDIntegrity_v2)
if hasattr(process, "DQM_FEDIntegrity_v3"):
process.schedule.remove(process.DQM_FEDIntegrity_v3)
if hasattr(process, "HLTAnalyzerEndpath"):
process.schedule.remove(process.HLTAnalyzerEndpath)
del process.HLTAnalyzerEndpath
#process.load("HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.printGenParticles_cff")
#process.generation_step *= process.printGenParticles
# Disable gen vertex smearing
process.VtxSmeared = cms.EDProducer("FlatEvtVtxGenerator",
MaxZ = cms.double(0.0),
MaxX = cms.double(0.0),
MaxY = cms.double(0.0),
MinX = cms.double(0.0),
MinY = cms.double(0.0),
MinZ = cms.double(0.0),
TimeOffset = cms.double(0.0),
src = cms.InputTag("generator")
)
return process
