param.trigger.triggerSrc.setProcessName(dataVersion.getTriggerProcess())
# Set tau selection mode (options: 'tauCandidateSelectionOnly', 'tauCandidateSelectionOnlyReversedRtau')
# other options (use not recommended here): 'standard'
param.setAllTauSelectionOperatingMode('tauCandidateSelectionOnly')
# Set tau sources to trigger matched tau collections
param.setAllTauSelectionSrcSelectedPatTausTriggerMatched()
# Trigger with scale factors (at the moment hard coded)
if applyTriggerScaleFactor and dataVersion.isMC():
param.triggerEfficiencyScaleFactor.mode = "scaleFactor"
# Set the data scenario for vertex/pileup weighting
param.setPileupWeightFor2011(dataVersion, era="EPS", method="intime") # Reweight by true PU distribution
param.setDataTriggerEfficiency(dataVersion, era="EPS")
#Reminder(from HChSignalAnalysisParameters_cff.py):
#def setTriggerPileupFor2011(**kwargs):
# setEfficiencyTriggersFor2011(**kwargs)
# setPileupWeightFor2011()
# Overwrite necessary values here
#param.trigger.hltMetCut = 45.0 # note: 45 is the minimum possible value for which HLT_MET is saved (see histogram hlt_met) attikis
#param.trigger.hltMetCut = 0.0
print "\nhltMetCut:", param.trigger.hltMetCut
param.InvMassVetoOnJets.setTrueToUseModule = False
##############################################################################
process.QCDMeasurement_PASJuly11 = cms.EDProducer("HPlusQCDMeasurement_PASJuly11Producer",
trigger = param.trigger,
#param.setAllTauSelectionSrcSelectedPatTaus()
param.setAllTauSelectionSrcSelectedPatTausTriggerMatched()
# Switch to PF2PAT objects
#param.changeCollectionsToPF2PAT()
param.changeCollectionsToPF2PAT(postfix=PF2PATVersion)
# Trigger with scale factors (at the moment hard coded)
if applyTriggerScaleFactor and dataVersion.isMC():
param.triggerEfficiencyScaleFactor.mode = "scaleFactor"
# Set the data scenario for vertex/pileup weighting
if len(options.puWeightEra) > 0:
puweight = options.puWeightEra
param.setPileupWeight(dataVersion, process=process, commonSequence=process.commonSequence, pset=param.vertexWeight, psetReader=param.pileupWeightReader, era=puweight) # Reweight by true PU distribution
param.setDataTriggerEfficiency(dataVersion, era=puweight)
print "PU weight era =",puweight
#param.trigger.selectionType = "disabled"
import HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.customisations as tauEmbeddingCustomisations
if options.tauEmbeddingInput != 0:
#tauEmbeddingCustomisations.addMuonIsolationEmbeddingForSignalAnalysis(process, process.commonSequence)
tauEmbeddingCustomisations.setCaloMetSum(process, process.commonSequence, options, dataVersion)
tauEmbeddingCustomisations.customiseParamForTauEmbedding(param, options, dataVersion)
if dataVersion.isMC():
process.muonTriggerFixSequence = cms.Sequence()
additionalCounters.extend(tauEmbeddingCustomisations.addMuonTriggerFix(process, dataVersion, process.muonTriggerFixSequence, options))
process.commonSequence.replace(process.patSequence, process.muonTriggerFixSequence*process.patSequence)
if tauEmbeddingFinalizeMuonSelection:
#applyIsolation = not doTauEmbeddingMuonSelectionScan
param.setAllTauSelectionOperatingMode('tauCandidateSelectionOnly')
# Set tau sources to trigger matched tau collections
param.setAllTauSelectionSrcSelectedPatTausTriggerMatched()
# Trigger with scale factors (at the moment hard coded)
if applyTriggerScaleFactor and dataVersion.isMC():
param.triggerEfficiencyScaleFactor.mode = "scaleFactor"
# Set the data scenario for vertex/pileup weighting
puweight = "Run2011A"
if len(options.puWeightEra) > 0:
puweight = options.puWeightEra
param.setPileupWeightFor2011(dataVersion,
era=puweight) # Reweight by true PU distribution
param.setDataTriggerEfficiency(dataVersion, era=puweight)
#Reminder(from HChSignalAnalysisParameters_cff.py):
#def setTriggerPileupFor2011(**kwargs):
# setEfficiencyTriggersFor2011(**kwargs)
# setPileupWeightFor2011()
# Overwrite necessary values here
#param.trigger.hltMetCut = 45.0 # note: 45 is the minimum possible value for which HLT_MET is saved (see histogram hlt_met) attikis
#param.trigger.hltMetCut = 0.0
print "\nhltMetCut:", param.trigger.hltMetCut
param.InvMassVetoOnJets.setTrueToUseModule = False
##############################################################################
process.QCDMeasurement = cms.EDProducer(
"HPlusQCDMeasurementProducer",
