# signalAnalysis.metTriggerEfficiencyScaleFactor.mode = "mcEfficiency"
# signalAnalysis.metTriggerEfficiencyScaleFactor.dataSelect = ["runs_190456_208686"]
print "Customisations done"
from HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration import ConfigBuilder
builder = ConfigBuilder(dataVersion, dataEras,
maxEvents=-1, # default is -1
customizeLightAnalysis=customize,
doQCDTailKillerScenarios=True,
applyTauTriggerScaleFactor=True,
#applyTauTriggerLowPurityScaleFactor=True,
applyMETTriggerScaleFactor=False,
applyL1ETMScaleFactor=True,
#doAgainstElectronScan=True,
doTauIDandMisIDSystematicsAsShapes=True,
doSystematics=True,
doAsymmetricTriggerUncertainties=True,
histogramAmbientLevel = "Vital",
allowTooManyAnalyzers=True,
#doOptimisation=True, optimisationScheme="metScenarios",
#doOptimisation=True, optimisationScheme="jetScenarios",
#doOptimisation=True, optimisationScheme="btagSymmetricScenarios",
#doOptimisation=True, optimisationScheme="myOptimisation"
)
process = builder.buildSignalAnalysis()
# An example how to use a non-default file(s)
#process.source.fileNames = [
# "Run2011B",
]
def customize(analysis):
analysis.trigger.selectionType = "disabled"
# analysis.embeddingMuonSrc = "tauEmbeddingMuons"
# analysis.tauPtCut = 41.0
# analysis.tauEtaCut = 2.1
from HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration import ConfigBuilder
builder = ConfigBuilder(dataVersion, dataEras,
pickEvents=False,
# maxEvents=1000, # default is -1
customizeLightAnalysis=customize,
#doAgainstElectronScan=True,
#doSystematics=True,
histogramAmbientLevel = "Systematics",
#doOptimisation=True, optimisationScheme=myOptimisation
)
process = builder.buildEwkBackgroundCoverageAnalysis()
if builder.options.doPat != 0:
process.source.fileNames = ["file:/mnt/flustre/wendland/AODSIM_PU_S6_START44_V9B_7TeV/Fall11_TTJets_TuneZ2_7TeV-madgraph-tauola_AODSIM_PU_S6_START44_V9B-v1_testfile.root"]
# We want some kind of preselection before running PAT
import HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.customisations as customisations
import HiggsAnalysis.HeavyChHiggsToTauNu.ewkBackgroundCoverageAnalysis as ewkBackgroundCoverageAnalysis
s = cms.Sequence()
process.preGenTausSequence = s
def customize(treeAnalysis):
# Disable filling of jet energy fractions as it crashes.
treeAnalysis.Tree.fillJetEnergyFractions=False
treeAnalysis.Tree.fillNonIsoLeptonVars=True
# Apply beta cut for jets to reject PU jets
#treeAnalysis.jetSelection.betaCut = 0.2 # Disable by setting to 0.0; if you want to enable, set to 0.2
# treeAnalysis.tauSelection.ptCut = 80.0 #
# treeAnalysis.MET.METCut = 100.0
from HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration import ConfigBuilder
builder = ConfigBuilder(dataVersion, dataEras,
maxEvents=-1, # default is -1
customizeLightAnalysis=customize,
doFillTree = True, #False (crashes if True)
tauSelectionOperatingMode="tauCandidateSelectionOnly",
#doAgainstElectronScan=True,
#doSystematics=True,
#histogramAmbientLevel = "Vital",
#doOptimisation=True, optimisationScheme="myOptimisation"
)
process = builder.buildSignalAnalysis()
if builder.options.tauEmbeddingInput != 0:
if builder.dataVersion.isMC():
process.source.fileNames = [
]
else:
# HLT_Mu40_eta2p1_v1
signalAnalysis.tauSelection.rtauCut = 0.0 # this disables the R_{tau} cut for efficiency measurement (gives better statistics)
print "Customisations done"
from HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration import ConfigBuilder
builder = ConfigBuilder(
dataVersion,
dataEras,
maxEvents=-1, # default is -1
customizeLightAnalysis=customize,
#doHeavyAnalysis=True,
#customizeHeavyAnalysis=customize,
#applyTauTriggerScaleFactor=False,
#applyTauTriggerLowPurityScaleFactor=True,
#applyMETTriggerScaleFactor=True,
#doTriggerMatching=False,
#useCHSJets=True,
#doQCDTailKillerScenarios=True,
#doAgainstElectronScan=True,
#doTauIsolationAndJetPUScan=True,
doBTagScan=True,
#doSystematics=True,
histogramAmbientLevel=
"Informative", # at least "Informative" required for eff. measurement
#doOptimisation=True, optimisationScheme="myOptimisation"
)
process = builder.buildSignalAnalysis()
# An example how to use a non-default file(s)
#process.source.fileNames = [
#signalAnalysis.bTagging.subleadingDiscriminatorCut = 0.244
#signalAnalysis.MET.METCut = 50.0
#signalAnalysis.MET.preMETCut = 30.0
print "Customisation applied"
from HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration import ConfigBuilder
builder = ConfigBuilder(dataVersion, dataEras,
maxEvents=-1, # default is -1
customizeLightAnalysis=customize,
doQCDTailKillerScenarios=True,
applyTauTriggerScaleFactor=True,
#applyTauTriggerLowPurityScaleFactor=True,
applyMETTriggerScaleFactor=False,
applyL1ETMScaleFactor=True,
#doAgainstElectronScan=True,
doTauIDandMisIDSystematicsAsShapes=True,
doSystematics=True,
doAsymmetricTriggerUncertainties=True,
histogramAmbientLevel = "Vital",
allowTooManyAnalyzers=True,
#doOptimisation=True, optimisationScheme="metScenarios",
#doOptimisation=True, optimisationScheme="jetScenarios",
#doOptimisation=True, optimisationScheme="btagSymmetricScenarios",
#doOptimisation=True, optimisationScheme="myOptimisation"
)
process = builder.buildSignalAnalysisInvertedTau()
# An example how to use a non-default file(s)
#process.source.fileNames = [
# "store/group/local/HiggsChToTauNuFullyHadronic/pattuples/CMSSW_5_3_X/TTJets_TuneZ2star_Summer12/TTJets_MassiveBinDECAY_TuneZ2star_8TeV-madgraph-tauola/Summer12_DR53X_PU_S10_START53_V7A_v1_AODSIM_pattuple_v53_1/cad8d1056ca20d363262a3efa1d97a74/pattuple_570_1_k4M.root"
import FWCore.ParameterSet.Config as cms
### FIXME: This file has not been tested after migration to AnalysisConfiguration
# Select the version of the data (needed only for interactice running
#dataVersion="44XmcS6" # Fall11 MC
dataVersion="44Xdata" # Run2011 08Nov and 19Nov ReRecos
dataEras = ["Run2011AB"]
from HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration import ConfigBuilder
builder = ConfigBuilder(dataVersion, dataEras)
# Always run PAT on the fly when running on PickEvents AOD file
builder.options.doPat=1
import HiggsAnalysis.HeavyChHiggsToTauNu.HChSignalAnalysisParameters_cff as param
builder.options.trigger = param.singleTauMetTriggerPaths
process = builder.buildSignalAnalysis()
# Set input file(s)
process.source.fileNames = [
"file:pickevents_Tau_Run2011A-May10ReReco-v1_AOD.root",
# "file:pickevents_Tau_Run2011A-PromptReco-v4_AOD.root",
]
# Add output module
process.out = cms.OutputModule("PoolOutputModule",
SelectEvents = cms.untracked.PSet(
SelectEvents = cms.vstring("signalAnalysisPath")
),
print "- Nprongs cut included in tau candidate selections:",signalAnalysis.applyNprongsCutForTauCandidate.value()
print "- Rtau cut included in tau candidate selections:",signalAnalysis.applyRtauCutForTauCandidate.value()
print "- Analysis method used:",signalAnalysis.analysisMode.value()
print "- MET cut:",signalAnalysis.MET.METCut.value()
print "- pre-MET cut:",signalAnalysis.MET.preMETCut.value()
print "- Tail-Killer:", signalAnalysis.QCDTailKiller
from HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration import ConfigBuilder
builder = ConfigBuilder(dataVersion, dataEras,
maxEvents=1000, # default is -1
customizeLightAnalysis=customize,
#doHeavyAnalysis=True,
#customizeHeavyAnalysis=customize,
#useCHSJets=True,
applyTauTriggerScaleFactor=True,
#applyTauTriggerLowPurityScaleFactor=True,
#applyMETTriggerScaleFactor=True,
tauSelectionOperatingMode="tauCandidateSelectionOnly",
#doSystematics=True,
#doQCDTailKillerScenarios=True,
doFillTree=False,
#histogramAmbientLevel = "Vital", # Informative by default
#doOptimisation=True, optimisationScheme="myOptimisation"
)
process = builder.buildQCDMeasurementFactorised()
if builder.options.tauEmbeddingInput != 0:
process.source.fileNames = [
#"file:/mnt/flustre/wendland/embedded_latest.root"
"file:/home/wendland/v25_embed/CMSSW_4_4_4/src/HiggsAnalysis/HeavyChHiggsToTauNu/test/tauEmbedding/embedded.root"
# For testing data
def customize(analysis):
analysis.trigger.selectionType = "disabled"
# analysis.embeddingMuonSrc = "tauEmbeddingMuons"
# analysis.tauPtCut = 41.0
# analysis.tauEtaCut = 2.1
from HiggsAnalysis.HeavyChHiggsToTauNu.AnalysisConfiguration import ConfigBuilder
builder = ConfigBuilder(
dataVersion,
dataEras,
pickEvents=False,
# maxEvents=1000, # default is -1
customizeLightAnalysis=customize,
#doAgainstElectronScan=True,
#doSystematics=True,
histogramAmbientLevel="Systematics",
#doOptimisation=True, optimisationScheme=myOptimisation
)
process = builder.buildEwkBackgroundCoverageAnalysis()
if builder.options.doPat != 0:
process.source.fileNames = [
"file:/mnt/flustre/wendland/AODSIM_PU_S6_START44_V9B_7TeV/Fall11_TTJets_TuneZ2_7TeV-madgraph-tauola_AODSIM_PU_S6_START44_V9B-v1_testfile.root"
]
# We want some kind of preselection before running PAT
import HiggsAnalysis.HeavyChHiggsToTauNu.tauEmbedding.customisations as customisations