print('Now starting')
# Not sure it's the correct way to add a new era
options.changeDefaults(era='2016')
print('Changed era')
# 2017 tags from https://twiki.cern.ch/twiki/bin/view/CMS/PdmV2017Analysis
if runOnData:
options.changeDefaults(globalTag='94X_dataRun2_ReReco_EOY17_v6',
process='RECO')
else:
options.changeDefaults(globalTag='94X_mc2017_realistic_v14', process='PAT')
print('Changed defaults')
framework = Framework.Framework(options)
print('Created framework object')
framework.addAnalyzer(
'tth_analyzer',
cms.PSet(
type=cms.string('tth_analyzer'),
prefix=cms.string('tth_'),
enable=cms.bool(True),
categories_parameters=cms.PSet(
# Per-category lepton pt cuts
mumu_leadingLeptonPtCut=cms.untracked.double(20), # muon
mumu_subleadingLeptonPtCut=cms.untracked.double(10), # muon
elel_leadingLeptonPtCut=cms.untracked.double(25), # electron
elel_subleadingLeptonPtCut=cms.untracked.double(15), # electron
muel_leadingLeptonPtCut=cms.untracked.double(25), # muon
import FWCore.ParameterSet.Config as cms
from Configuration.StandardSequences.Eras import eras
from cp3_llbb.Framework import Framework
process = Framework.create(False, eras.Run2_25ns, '74X_mcRun2_asymptotic_v2', redoJEC=True)
process.source.fileNames = cms.untracked.vstring(
'file:///afs/cern.ch/user/s/sbrochet/public/CP3/TTJets_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8_Asympt25ns_MCRUN2_74_V9_reduced.root'
)
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(10))
import FWCore.ParameterSet.Config as cms
from Configuration.StandardSequences.Eras import eras
from cp3_llbb.Framework import Framework
process = Framework.create(True, eras.Run2_50ns, "74X_dataRun2_v2", redoJEC=True, process_name="RECO")
process.framework.analyzers.dilepton = cms.PSet(
type=cms.string("dilepton_analyzer"),
prefix=cms.string("dilepton_"),
enable=cms.bool(True),
categories_parameters=cms.PSet(mll_cut=cms.untracked.double(20)),
parameters=cms.PSet(
standalone=cms.untracked.bool(True),
muons_wp=cms.untracked.string("loose"),
electrons_wp=cms.untracked.string("loose"),
),
)
process.framework.analyzers.test = cms.PSet(
type=cms.string("test_analyzer"), prefix=cms.string("test_"), enable=cms.bool(True)
)
process.framework.analyzers.bTagsLoose = cms.PSet(
type=cms.string("btags_analyzer"),
prefix=cms.string("btags_CSVv2_loose"),
enable=cms.bool(True),
parameters=cms.PSet(
discr_name=cms.untracked.string("pfCombinedInclusiveSecondaryVertexV2BJetTags"),
discr_cut=cms.untracked.double(0.605),
eta_cut=cms.untracked.double(2.4),
process = Framework.create(
True,
eras.Run2_25ns,
'74X_dataRun2_v2',
cms.PSet(tt=cms.PSet(
type=cms.string('tt_analyzer'),
prefix=cms.string('tt_'),
enable=cms.bool(True),
parameters=cms.PSet(
electronPtCut=cms.untracked.double(20),
electronEtaCut=cms.untracked.double(2.5),
electronVetoIDName=cms.untracked.string(
'cutBasedElectronID-Spring15-25ns-V1-standalone-veto'),
electronLooseIDName=cms.untracked.string(
'cutBasedElectronID-Spring15-25ns-V1-standalone-loose'),
electronMediumIDName=cms.untracked.string(
'cutBasedElectronID-Spring15-25ns-V1-standalone-medium'),
electronTightIDName=cms.untracked.string(
'cutBasedElectronID-Spring15-25ns-V1-standalone-tight'),
muonPtCut=cms.untracked.double(20),
muonEtaCut=cms.untracked.double(2.4),
muonLooseIsoCut=cms.untracked.double(
.20), # Loose cut recommended for dilepton analysis
muonTightIsoCut=cms.untracked.double(.12),
jetPtCut=cms.untracked.double(30),
jetEtaCut=cms.untracked.double(2.5),
#jetPUID = cms.untracked.double(-9999999),
jetDRleptonCut=cms.untracked.double(0.3),
jetID=cms.untracked.string(
'loose'
), # not tightLeptonVeto since DeltaR(l,j) cut should be enough
jetCSVv2Name=cms.untracked.string(
'pfCombinedInclusiveSecondaryVertexV2BJetTags'),
jetCSVv2L=cms.untracked.double(0.605),
jetCSVv2M=cms.untracked.double(0.89),
jetCSVv2T=cms.untracked.double(0.97),
hltDRCut=cms.untracked.double(
0.3), # DeltaR cut for trigger matching
hltDPtCut=cms.untracked.double(
0.5), #Delta(Pt)/Pt cut for trigger matching
),
categories_parameters=cms.PSet(
MllCutSF=cms.untracked.double(20),
MllCutDF=cms.untracked.double(20),
MllZVetoCutLow=cms.untracked.double(76),
MllZVetoCutHigh=cms.untracked.double(116),
HLTDoubleMuon=cms.untracked.vstring(
'HLT_Mu17_TrkIsoVVL_(Tk)?Mu8_TrkIsoVVL_DZ_v.*'),
HLTDoubleEG=cms.untracked.vstring(
'HLT_Ele17_Ele12_CaloIdL_TrackIdL_IsoVL_DZ_v.*'),
HLTMuonEG=cms.untracked.vstring(
'HLT_Mu17_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL_v.*',
'HLT_Mu8_TrkIsoVVL_Ele17_CaloIdL_TrackIdL_IsoVL_v.*'),
),
)),
redoJEC=False,
process_name='RECO')
process = Framework.create(False, eras.Run2_25ns, '74X_mcRun2_asymptotic_v2', cms.PSet(
dilepton = cms.PSet(
type = cms.string('dilepton_analyzer'),
prefix = cms.string('dilepton_'),
enable = cms.bool(True),
categories_parameters = cms.PSet(
mll_cut = cms.untracked.double(20)
),
parameters = cms.PSet(
standalone = cms.untracked.bool(True),
muons_wp = cms.untracked.string('loose'),
electrons_wp = cms.untracked.string('loose')
)
),
bTagsLoose = cms.PSet(
type = cms.string('btags_analyzer'),
prefix = cms.string('btags_CSVv2_loose'),
enable = cms.bool(True),
parameters = cms.PSet(
discr_name = cms.untracked.string('pfCombinedInclusiveSecondaryVertexV2BJetTags'),
discr_cut = cms.untracked.double(0.605),
eta_cut = cms.untracked.double(2.4),
pt_cut = cms.untracked.double(30)
)
),
bTagsMedium = cms.PSet(
type = cms.string('btags_analyzer'),
prefix = cms.string('btags_CSVv2_medium'),
enable = cms.bool(True),
parameters = cms.PSet(
discr_name = cms.untracked.string('pfCombinedInclusiveSecondaryVertexV2BJetTags'),
discr_cut = cms.untracked.double(0.89),
eta_cut = cms.untracked.double(2.4),
pt_cut = cms.untracked.double(30)
)
),
bTagsTight = cms.PSet(
type = cms.string('btags_analyzer'),
prefix = cms.string('btags_CSVv2_tight'),
enable = cms.bool(True),
parameters = cms.PSet(
discr_name = cms.untracked.string('pfCombinedInclusiveSecondaryVertexV2BJetTags'),
discr_cut = cms.untracked.double(0.97),
eta_cut = cms.untracked.double(2.4),
pt_cut = cms.untracked.double(30)
)
),
test = cms.PSet(
type = cms.string('test_analyzer'),
prefix = cms.string('test_'),
enable = cms.bool(True)
)
),
redoJEC=True,
)
process = Framework.create(
runOnData,
eras.Run2_25ns,
globalTag,
cms.PSet(
hh_analyzer=cms.PSet(
type=cms.string("hh_analyzer"),
prefix=cms.string("hh_"),
enable=cms.bool(True),
categories_parameters=cms.PSet(mll_cut=cms.untracked.double(10)),
parameters=cms.PSet(
# Here are the default value (just to show what is configurable)
electronIsoCut_EB_Loose=cms.untracked.double(
0.0893
), # https://twiki.cern.ch/twiki/bin/view/CMS/CutBasedElectronIdentificationRun2
electronIsoCut_EE_Loose=cms.untracked.double(0.121),
electronIsoCut_EB_Tight=cms.untracked.double(0.0354),
electronIsoCut_EE_Tight=cms.untracked.double(0.0646),
electronPtCut=cms.untracked.double(20),
electronEtaCut=cms.untracked.double(2.5),
muonIsoCut=cms.untracked.double(0.20), # https://twiki.cern.ch/twiki/bin/view/CMS/TopMUO
muonPtCut=cms.untracked.double(20),
muonEtaCut=cms.untracked.double(2.4),
electrons_loose_wp_name=cms.untracked.string("cutBasedElectronID-Spring15-25ns-V1-standalone-loose"),
electrons_tight_wp_name=cms.untracked.string("cutBasedElectronID-Spring15-25ns-V1-standalone-tight"),
jetEtaCut=cms.untracked.double(2.4),
jetPtCut=cms.untracked.double(20),
discr_name=cms.untracked.string("pfCombinedInclusiveSecondaryVertexV2BJetTags"),
discr_cut_loose=cms.untracked.double(0.605),
discr_cut_medium=cms.untracked.double(0.89),
discr_cut_tight=cms.untracked.double(0.97),
minDR_l_j_Cut=cms.untracked.double(0.3),
hltDRCut=cms.untracked.double(9),
hltDPtCut=cms.untracked.double(1000), # cut will be DPt/Pt < hltDPtCut
),
)
),
redoJEC=True,
process_name=processName,
)
process = Framework.create(False, eras.Run2_25ns, '74X_mcRun2_asymptotic_v2', cms.PSet(
tt = cms.PSet(
type = cms.string('tt_analyzer'),
prefix = cms.string('tt_'),
enable = cms.bool(True),
parameters = cms.PSet(
electronPtCut = cms.untracked.double(20),
electronEtaCut = cms.untracked.double(2.5),
electronVetoIDName = cms.untracked.string('cutBasedElectronID-Spring15-25ns-V1-standalone-veto'),
electronLooseIDName = cms.untracked.string('cutBasedElectronID-Spring15-25ns-V1-standalone-loose'),
electronMediumIDName = cms.untracked.string('cutBasedElectronID-Spring15-25ns-V1-standalone-medium'),
electronTightIDName = cms.untracked.string('cutBasedElectronID-Spring15-25ns-V1-standalone-tight'),
muonPtCut = cms.untracked.double(20),
muonEtaCut = cms.untracked.double(2.4),
muonLooseIsoCut = cms.untracked.double(.20), # Loose cut recommended for dilepton analysis
muonTightIsoCut = cms.untracked.double(.12),
jetPtCut = cms.untracked.double(30),
jetEtaCut = cms.untracked.double(2.5),
#jetPUID = cms.untracked.double(-9999999),
jetDRleptonCut = cms.untracked.double(0.3),
jetID = cms.untracked.string('loose'), # not tightLeptonVeto since DeltaR(l,j) cut should be enough
jetCSVv2Name = cms.untracked.string('pfCombinedInclusiveSecondaryVertexV2BJetTags'),
jetCSVv2L = cms.untracked.double(0.605),
jetCSVv2M = cms.untracked.double(0.89),
jetCSVv2T = cms.untracked.double(0.97),
hltDRCut = cms.untracked.double(0.3), # DeltaR cut for trigger matching
hltDPtCut = cms.untracked.double(0.5), #Delta(Pt)/Pt cut for trigger matching
),
categories_parameters = cms.PSet(
MllCutSF = cms.untracked.double(20),
MllCutDF = cms.untracked.double(20),
MllZVetoCutLow = cms.untracked.double(76),
MllZVetoCutHigh = cms.untracked.double(116),
HLTDoubleMuon = cms.untracked.vstring('HLT_Mu17_TrkIsoVVL_(Tk)?Mu8_TrkIsoVVL_DZ_v.*'),
HLTDoubleEG = cms.untracked.vstring('HLT_Ele17_Ele12_CaloIdL_TrackIdL_IsoVL_DZ_v.*'),
HLTMuonEG = cms.untracked.vstring('HLT_Mu17_TrkIsoVVL_Ele12_CaloIdL_TrackIdL_IsoVL_v.*', 'HLT_Mu8_TrkIsoVVL_Ele17_CaloIdL_TrackIdL_IsoVL_v.*'),
),
)
), redoJEC=False
)
import FWCore.ParameterSet.Config as cms
from Configuration.StandardSequences.Eras import eras
from cp3_llbb.Framework import Framework
process = Framework.create(True, eras.Run2_50ns, '74X_dataRun2_v2', redoJEC=True)
process.framework.analyzers.dilepton = cms.PSet(
type = cms.string('dilepton_analyzer'),
prefix = cms.string('dilepton_'),
enable = cms.bool(True),
categories_parameters = cms.PSet(
mll_cut = cms.untracked.double(20)
),
parameters = cms.PSet(
standalone = cms.untracked.bool(True),
muons_wp = cms.untracked.string('loose'),
electrons_wp = cms.untracked.string('loose')
)
)
process.framework.analyzers.test = cms.PSet(
type = cms.string('test_analyzer'),
prefix = cms.string('test_'),
enable = cms.bool(True)
)
process.framework.analyzers.bTagsLoose = cms.PSet(
type = cms.string('btags_analyzer'),
prefix = cms.string('btags_CSVv2_loose'),
import FWCore.ParameterSet.Config as cms
from Configuration.StandardSequences.Eras import eras
from cp3_llbb.Framework import Framework
process = Framework.create(
False, eras.Run2_25ns, "74X_mcRun2_asymptotic_v2", redoJEC=True, JECDatabase="Summer15_25nsV5_MC.db"
)
process.source.fileNames = cms.untracked.vstring(
"file:///afs/cern.ch/user/s/sbrochet/public/CP3/TTJets_TuneCUETP8M1_13TeV-amcatnloFXFX-pythia8_RunIISpring15MiniAODv2_reduced.root"
)
process.maxEvents = cms.untracked.PSet(input=cms.untracked.int32(10))
import FWCore.ParameterSet.Config as cms
from Configuration.StandardSequences.Eras import eras
from cp3_llbb.Framework import Framework
process = Framework.create(True, eras.Run2_50ns, '74X_dataRun2_v2', redoJEC=True, process_name='RECO')
process.source.fileNames = cms.untracked.vstring(
'file:///afs/cern.ch/user/s/sbrochet/public/CP3/DoubleMuon_Run2015D_PromptReco_reduced.root'
)
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(10))
import FWCore.ParameterSet.Config as cms
from Configuration.StandardSequences.Eras import eras
from cp3_llbb.Framework import Framework
from cp3_llbb.Framework import METProducer
runOnData = False
globalTag_ = '74X_mcRun2_asymptotic_v5'
processName_ = 'PAT'
if runOnData:
globalTag_ = '74X_dataRun2_reMiniAOD_v2'
processName_ = 'RECO'
framework = Framework.Framework(runOnData,
eras.Run2_25ns,
globalTag=globalTag_,
processName=processName_)
framework.addAnalyzer(
'hh_analyzer',
cms.PSet(
type=cms.string('hh_analyzer'),
prefix=cms.string('hh_'),
enable=cms.bool(True),
parameters=cms.PSet(
# Producers
electronsProducer=cms.string('electrons'),
muonsProducer=cms.string('muons'),
jetsProducer=cms.string('jets'),
metProducer=cms.string('met'),
nohfMETProducer=cms.string('nohf_met'),
# Here are the default value (just to show what is configurable)
import FWCore.ParameterSet.Config as cms
from Configuration.StandardSequences.Eras import eras
from cp3_llbb.Framework import Framework
from cp3_llbb.TTWAnalysis.Configuration import addTTWAnalyzer, addTTWCandidatesAnalyzer, customizeProducers
globalTag_ = '76X_dataRun2_16Dec2015_v0'
processName_ = 'RECO'
framework = Framework.Framework(True, eras.Run2_25ns, globalTag=globalTag_, processName=processName_)
## ANALYZERS
addTTWAnalyzer (framework) ## make candidates
addTTWCandidatesAnalyzer(framework) ## fill most branches
## PRODUCERS
customizeProducers(framework)
# framework.redoJEC()
# framework.smearJets()
# framework.doSystematics(['jec', 'jer'])
framework.doSystematics([])
process = framework.create()
process.source.fileNames = cms.untracked.vstring(
'/store/data/Run2015D/DoubleMuon/MINIAOD/16Dec2015-v1/10000/00039A2E-D7A7-E511-98EE-3417EBE64696.root'
)
process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(1000))
