from CMGTools.RootTools.RootTools import *
## pathsAndFilters = {
## "HLT_Mu17_Mu8_v16": ("hltDiMuonMu17Mu8DzFiltered0p2","hltDiMuonMu17Mu8DzFiltered0p2"),
## "HLT_Mu17_Mu8_v17": ("hltL3fL1DoubleMu10MuOpenL1f0L2f10L3Filtered17","hltL3fL1DoubleMu10MuOpenL1f0L2f10L3Filtered17"),
## "HLT_Mu17_Mu8_v18": ("hltL3fL1DoubleMu10MuOpenOR3p5L1f0L2f10L3Filtered17","hltL3fL1DoubleMu10MuOpenOR3p5L1f0L2f10L3Filtered17"),
## "HLT_Mu17_Mu8_v19": ("hltL3fL1DoubleMu10MuOpenOR3p5L1f0L2f10L3Filtered17","hltL3fL1DoubleMu10MuOpenOR3p5L1f0L2f10L3Filtered17"),
## "HLT_Mu17_Mu8_v20": ("hltL3fL1DoubleMu10MuOpenOR3p5L1f0L2f10L3Filtered17","hltL3fL1DoubleMu10MuOpenOR3p5L1f0L2f10L3Filtered17"),
## "HLT_Mu17_Mu8_v21": ("hltDiMuonGlb17Glb8DzFiltered0p2","hltDiMuonGlb17Glb8DzFiltered0p2"),
## "HLT_Mu17_Mu8_v22": ("hltDiMuonGlb17Glb8DzFiltered0p2","hltDiMuonGlb17Glb8DzFiltered0p2"),
## }
pathsAndFilters = []
eventSelector = cfg.Analyzer(
'EventSelector',
toSelect = [
# here put the event numbers (actual event numbers from CMSSW)
]
)
jsonAna = cfg.Analyzer(
'JSONAnalyzer',
)
## triggerAna = cfg.Analyzer(
## 'TriggerAnalyzer',
## )
# this analyzer is just there to select a list of good primary vertices.
# in 2012, we don't use it to reweight events for pile-up
vertexAna = cfg.Analyzer(
mc_tauEffWeight = 'effTau2011A'
mc_muEffWeight = 'effMu2011A'
elif period == 'Period_2011B':
mc_vertexWeight = 'vertexWeight2011B'
mc_tauEffWeight = 'effTau2011B'
mc_muEffWeight = 'effMu2011B'
elif period == 'Period_2011AB':
mc_vertexWeight = 'vertexWeight2011AB'
mc_tauEffWeight = 'effTau2011AB'
mc_muEffWeight = 'effMu2011AB'
ZMuMuAna = cfg.Analyzer('ZMuMuAnalyzer',
pt1=20,
pt2=20,
iso1=0.1,
iso2=0.1,
eta1=2,
eta2=2,
m_min=0,
m_max=200)
triggerAna = cfg.Analyzer('TriggerAnalyzer')
jetAna = cfg.Analyzer('SimpleJetAnalyzer', ptCut=0)
effMuAna = cfg.Analyzer(
'EfficiencyAnalyzer',
# recsel = 'cuts_vbtfmuon'
genPdgId=13)
vertexAna = cfg.Analyzer('VertexAnalyzer',
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.RootTools.RootTools import *
# this analyzer finds the initial events before the skim
skimAnalyzer = cfg.Analyzer('skimAnalyzerCount')
jsonAna = cfg.Analyzer('JSONAnalyzer', )
triggerAna = cfg.Analyzer(
#'TriggerAnalyzer',
'triggerBitFilter', )
# this analyzer is just there to select a list of good primary vertices.
ttHVertexAna = cfg.Analyzer('VertexAnalyzer',
goodVertices='offlinePrimaryVertices',
vertexWeight=None,
fixedWeight=1,
verbose=False)
# Gen Info Analyzer
ttHGenAna = cfg.Analyzer(
'ttHGenLevelAnalyzer',
filterHiggsDecays=[0, 15, 23, 24],
verbose=False,
PDFWeights=[],
)
# Lepton Analyzer
useParked = False
MC_trigger_matching = True
puFileMC = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/07-01-13/MC_Summer12_PU_S10-600bins.root'
puFileData = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/07-01-13/Data_Pileup_2012_Moriond-600bins.root'
mc_vertexWeight = None
mc_tauEffWeight_mc = 'effLooseTau15MC'
mc_tauEffWeight = 'eff2012IsoTau19fb'
mc_jetEffWeight = 'eff2012Jet19fb'
if use1prong or useParked:
mc_tauEffWeight='eff2012IsoTau1prong19fb'
triggerAna = cfg.Analyzer(
'TriggerAnalyzer',
verbose = False
)
TauTauAna = cfg.Analyzer(
'TauTauAnalyzer',
pt1 = 0,
pt2 = 0,
iso1 = 999,
iso2 = 999,
eta1 = 999,
eta2 = 999,
m_min = 0,
m_max = 99999,
diLeptonCutString = '',
# No trigger matching for the moment since MC only contains too high threshold trigger
triggerMap = pathsAndFilters,
import glob
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.RootTools.RootTools import *
reader = cfg.Analyzer(
'ObjectReader',
# muonCol = ('cmgMuonSel', 'vector<cmg::Muon>'),
# electronCol = ('cmgElectronSel', 'vector<cmg::Electron>'),
muonCol = None,
electronCol = None,
# only if you have PFJets
jetCol = None,
# jetCol = ('ak5PFJets', 'std::vector<reco::PFJet>'),
simJetCol = ('jets', 'std::vector<reco::LeafCandidate>'),
# genJetCol = ('ak5GenJets', 'std::vector<reco::GenJet>'),
genJetCol = ('genJets', 'std::vector<reco::LeafCandidate>'),
genPartCol = ('genParticles', 'std::vector<reco::GenParticle>'),
simPartCol = ('pfsim', 'std::vector<reco::LeafCandidate>'),
)
jetAna = cfg.Analyzer(
'PFSimJetAnalyzer',
jetPt = 10.,
jetEta = 5.0,
btagSFseed = 123456,
relaxJetId = True,
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
fourJetAna = cfg.Analyzer(
'FourJetAnalyzer',
npfj=(3, 5), # at least 8 PF particle in 3 of the 4 jets
ntkj=(4, 1), # at least 1 PF charged hadron in 4 of the 4 jets
minM=0.5, # No jet with a jet mass smaller than 2.5 GeV
mVis=180., # total visible mass in excess of 180 GeV,
chi2=
1000., # total chi**2 of the beta4 "fit" smaller than 1000. (i.e., all fit energies > 0.)
pair_mass=0., # No jet pair should have a mass smaller than 30 GeV
pair_cijkl=
1000., # Tue smallest sum of cos(ij)+cos(kl) should be larger than -1.1
pair_sumtet=
0., # The sum of the four smallest jet-jet angles must be larger than 350 degrees
h_mass=105., # The Higgs candidate mass should be larger than 90 GeV
z_mass=(
80., 110.
), # The Z candidate jet pair should have a mass closer to the Z mass than 100 GeV
)
def createTreeProducer(ana):
tp = cfg.Analyzer('_'.join(['FourJetTreeProducer', 'FourJetAnalyzer']),
anaName='FourJetAnalyzer')
return tp
# for debugging
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.WMass.triggerMap import triggers_mu
jsonAna = cfg.Analyzer('JSONAnalyzer', )
triggerAna = cfg.Analyzer('triggerBitFilter',
# keepFailingEvents = False
)
vertexAna = cfg.Analyzer(
'VertexAnalyzer',
allVertices='slimmedPrimaryVertices',
goodVertices='slimmedPrimaryVertices',
vertexWeight=None,
fixedWeight=1,
verbose=False,
)
WAna = cfg.Analyzer('WAnalyzer',
recoilcut=1000,
pfmetcut=0,
jetptcut=1000,
pt=30,
eta=2.1,
iso=0.5,
savegenp=False,
verbose=True,
triggerBits={'SingleMu': triggers_mu},
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.H2TauTau.triggerMap import pathsAndFilters
from CMGTools.H2TauTau.proto.weights.weighttable import mu_id_taumu_2012, mu_iso_taumu_2012
from CMGTools.H2TauTau.proto.samples.sampleShift import selectShift
from CMGTools.RootTools.RootTools import *
mc_vertexWeight = 'vertexWeightFall112012AB'
jsonAna = cfg.Analyzer('JSONAnalyzer', )
vertexAna = cfg.Analyzer(
'VertexAnalyzer',
goodVertices='offlinePrimaryVertices',
vertexWeight=mc_vertexWeight,
# fixedWeight = 1,
verbose=False)
ZMuMuAna = cfg.Analyzer('ZMuMuAnalyzer',
pt1=10,
eta1=2.4,
iso1=0.3,
pt2=10,
eta2=2.4,
iso2=0.3,
m_min=50,
m_max=120,
verbose=False)
# minM = 0.5, # No jet with a jet mass smaller than 2.5 GeV
# mVis = 180., # total visible mass in excess of 180 GeV,
# chi2 = 1000., # total chi**2 of the beta4 "fit" smaller than 1000. (i.e., all fit energies > 0.)
# pair_mass = 0., # No jet pair should have a mass smaller than 30 GeV
# pair_cijkl = 1000., # Tue smallest sum of cos(ij)+cos(kl) should be larger than -1.1
# pair_sumtet = 0., # The sum of the four smallest jet-jet angles must be larger than 350 degrees
# h_mass = 105., # The Higgs candidate mass should be larger than 90 GeV
# z_mass = (80.,110.), # The Z candidate jet pair should have a mass closer to the Z mass than 100 GeV
# )
httAna = cfg.Analyzer(
'httanalyzer',
npfj=(2, 3), # at least 8 PF particle in 3 of the 4 jets
ntkj=(2, 1),
minM=0.5,
mVis=100.,
chi2=1000.,
minE=10.,
h_mass=105.,
z_mass=(80., 110.),
)
def createTreeProducer(ana):
tp = cfg.Analyzer('_'.join(['htttreeproducer', 'httanalyzer']),
anaName='httanalyzer')
return tp
# for debugging
stopper = cfg.Analyzer('Stopper')
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.H2TauTau.proto.samples.getFiles import getFiles
from CMGTools.RootTools.RootTools import *
WNJetsAna = cfg.Analyzer('WNJetsAnalyzer', verbose=False)
#########################################################################################
WJets = cfg.MCComponent(
name='WJets',
files=getFiles(
'/WJetsToLNu_TuneZ2_7TeV-madgraph-tauola/Fall11-PU_S6_START42_V14B-v1/AODSIM/V5_B/PAT_CMG_V5_6_0_B',
'cmgtools', 'cmgTuple.*root'),
xSection=31314.,
nGenEvents=1,
triggers=[],
effCorrFactor=1)
#########################################################################################
selectedComponents = [WJets]
sequence = cfg.Sequence([
WNJetsAna,
])
comp = WJets
ttHLepSkim.maxLeptons = 999
ttHLepSkim.ptCuts = [5, 3]
# --- JET-LEPTON CLEANING ---
ttHJetAna.minLepPt = 20
# otherwise with only absIso cut at 10 GeV and no relIso we risk cleaning away good jets
# --- JET-MET SKIMMING ---
ttHJetMETSkim.jetPtCuts = [
100,
]
ttHJetMETSkim.metCut = 100
# Event Analyzer for susy multi-lepton (at the moment, it's the TTH one)
ttHEventAna = cfg.Analyzer(
'ttHLepEventAnalyzer',
minJets25=0,
)
from CMGTools.TTHAnalysis.samples.samples_8TeV_v517 import triggers_1mu, triggers_mumu, triggers_ee, triggers_mue, triggers_MET150, triggers_HT650, triggers_HTMET
# Tree Producer
treeProducer = cfg.Analyzer(
'treeProducerSusySoftlepton',
vectorTree=True,
PDFWeights=PDFWeights,
triggerBits={
'SingleMu': triggers_1mu,
'DoubleMu': triggers_mumu,
'DoubleEl': [t for t in triggers_ee if "Ele15_Ele8_Ele5" not in t],
'TripleEl': [t for t in triggers_ee if "Ele15_Ele8_Ele5" in t],
'MuEG': [t for t in triggers_mue if "Mu" in t and "Ele" in t],
'MET15': triggers_MET150,
import copy
import os
import glob
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.RootTools.RootTools import *
reader = cfg.Analyzer(
'ObjectReader',
# muonCol = ('cmgMuonSel', 'vector<cmg::Muon>'),
# electronCol = ('cmgElectronSel', 'vector<cmg::Electron>'),
muonCol=None,
electronCol=None,
# only if you have PFJets
jetCol=('ak5PFJets', 'std::vector<reco::PFJet>'),
# comment above and uncomment below when doing gensim
# jetCol = None,
simJetCol=('jets', 'std::vector<reco::LeafCandidate>'),
genJetCol=('ak5GenJets', 'std::vector<reco::GenJet>'),
# comment above and uncomment below when doing gensim
# genJetCol = ('genJets', 'std::vector<reco::LeafCandidate>'),
genPartCol=('genParticles', 'std::vector<reco::GenParticle>'),
simPartCol=('pfsim', 'std::vector<reco::LeafCandidate>'),
recPartCol=('particleFlow', 'std::vector<reco::PFCandidate>'))
jetAna = cfg.Analyzer(
'PFSimJetAnalyzer',
jetPt=1.0,
jetEta=2.0,
btagSFseed=123456,
relaxJetId=True,
mc_vertexWeight = 'vertexWeight2011AB'
mc_tauEffWeight = 'effTau2011AB'
mc_muEffWeight = 'effMu2011AB'
jetAna = cfg.Analyzer(
'SimpleJetAnalyzer',
ptCut=0,
genPtCut=0,
## use pat::Jets
jetCollection=('selectedPatJetsAK5NoPUSub', 'std::vector<pat::Jet>'),
## or cmg::Jets
## jetCollection = ('cmgPFJetSel','std::vector<cmg::PFJet>'),
doJetIdHisto=False,
applyPFLooseId=False,
jetIdMva=(0, "%s/src/CMGTools/External/data/mva_JetID.weights.xml" %
os.getenv("CMSSW_BASE"), "JetID"),
## use gen jets in Loic's PATtuples
## genJetsCollection = ('prunedGen','std::vector<reco::GenParticle>'),
## use the gen jets available also in the CMG trees
## genJetsCollection = ('genJetSel','std::vector< cmg::PhysicsObjectWithPtr< edm::Ptr< reco::GenJet > > >'),
## use the gen jets present only in PATtuples
genJetsCollection=(('selectedPatJetsAK5NoPUSub', 'genJets'),
'vector<reco::GenJet>'),
## use level2 gen leptons to make the matching
useGenLeptons=True)
#########################################################################################
DYJetsFall11_PAT = cfg.MCComponent(
def regionNames( self ):
names = []
for dim in self.dimensions:
names.append(dim.regionNames())
# yeah!
allCombs = list(itertools.product(*names))
allCombNames = map( '_'.join, allCombs )
return allCombNames
if __name__ == '__main__':
import CMGTools.RootTools.fwlite.Config as cfg
vbfKwargs = dict( Mjj = 400,
deltaEta = 4.0
)
cuts = cfg.Analyzer(
'Dummy',
# vertexWeight = mc_vertexWeight,
MT_low = 40,
MT_high = 60,
Boosted_JetPt = 150,
**vbfKwargs
)
reg = H2TauTauRegions( cuts )
print reg.regionNames()
'HLT_Ele15_{blah}_LooseIsoPFTau15_v4'.format(blah=blah),
'HLT_Ele15_{blah}_LooseIsoPFTau20_v6'.format(blah=blah),
'HLT_Ele15_{blah}_LooseIsoPFTau20_v8'.format(blah=blah),
'HLT_Ele15_{blah}_LooseIsoPFTau20_v9'.format(blah=blah),
'HLT_Ele15_{blah}_TightIsoPFTau20_v2'.format(blah=blah),
'HLT_Ele18_{blah}_MediumIsoPFTau20_v1'.format(blah=blah)
]
data_triggers_11B = [
'HLT_Ele18_{blah}_MediumIsoPFTau20_v1'.format(blah=blah),
'HLT_Ele20_{blah}_MediumIsoPFTau20_v5'.format(blah=blah),
'HLT_Ele20_{blah}_MediumIsoPFTau20_v6'.format(blah=blah)
]
triggerAna = cfg.Analyzer(
'TriggerAnalyzer'
)
tauEleAna = cfg.Analyzer(
'TauEleAnalyzer',
pt1 = 20,
pt2 = 20,
iso1 = 999,
iso2 = 0.1,
eta1 = 999,
eta2 = 2.1,
# cutString2 = 'cuts_vbtf80ID',
# eta2 = 1.4,
m_min = 10,
m_max = 99999,
diLeptonCutString = 'cuts_baseline'
def createTreeProducer(ana):
tp = cfg.Analyzer('_'.join(['htttreeproducer', 'httanalyzer']),
anaName='httanalyzer')
return tp
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.WMass.triggerMap import pathsAndFilters
#from CMGTools.WMass.analyzers.WAnalyzer import WAnalyzer
#from CMGTools.WMass.analyzers.ZAnalyzer_notrig import ZAnalyzer_notrig
jsonAna = cfg.Analyzer('JSONAnalyzer', )
triggerAna = cfg.Analyzer('TriggerAnalyzer', keepFailingEvents=True)
vertexAna = cfg.Analyzer('VertexAnalyzer',
fixedWeight=1.,
keepFailingEvents=True)
WAna = cfg.Analyzer('WAnalyzer',
recoilcut=1000,
pfmetcut=0,
jetptcut=1000,
pt=30,
eta=2.1,
iso=0.5,
savegenp=True,
verbose=True,
triggerMap=pathsAndFilters)
WtreeProducer = cfg.Analyzer('WTreeProducer')
ZAna = cfg.Analyzer('ZAnalyzer',
recoilcut=1000,
puFileMC = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/12-09-12/MC_Summer12_PU_S7.root'
puFileData = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/12-09-12/Data_Pileup_2012.root'
vertexFileDir = os.environ[
'CMSSW_BASE'] + '/src/CMGTools/RootTools/data/Reweight/2012/Vertices'
vertexFileData = '/'.join(
[vertexFileDir, 'vertices_data_2012A_2012B_start_195947.root'])
mc_vertexWeight = None
mc_tauEffWeight_mc = 'eff2012Tau20MC_TauEle'
mc_eleEffWeight_mc = 'eff2012Ele20MC'
mc_tauEffWeight = 'effTau2012AB_TauEle'
mc_eleEffWeight = 'effEle2012AB'
jsonAna = cfg.Analyzer('JSONAnalyzer', )
triggerAna = cfg.Analyzer('TriggerAnalyzer',
# verbose = True
)
vertexAna = cfg.Analyzer('VertexAnalyzer',
goodVertices='goodPVFilter',
vertexWeight=mc_vertexWeight,
fixedWeight=1,
verbose=False)
embedWeighter = cfg.Analyzer('EmbedWeighter', verbose=False)
pileUpAna = cfg.Analyzer('PileUpAnalyzer', true=True)
def createTreeProducer(ana):
tp = cfg.Analyzer('_'.join(['FourJetTreeProducer', 'FourJetAnalyzer']),
anaName='FourJetAnalyzer')
return tp
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
tautaullAna = cfg.Analyzer('TauTauLLAnalyzer', )
def createTreeProducer(ana):
tp = cfg.Analyzer('_'.join(['TauTauLLTreeProducer', 'TauTauLLAnalyzer']),
anaName='TauTauLLAnalyzer')
return tp
# for debugging
stopper = cfg.Analyzer('Stopper')
sequence = cfg.Sequence([
tautaullAna,
createTreeProducer(tautaullAna),
])
##########################################################
## SUSY CONFIGURATION FOR LEPTON ID STUDIES ##
## makes trees with one entry per lepton, not per event ##
##########################################################
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.RootTools.RootTools import *
PDFWeights = []
#PDFWeights = [ ("CT10",53), ("MSTW2008lo68cl",41), ("NNPDF21_100",101) ]
#PDFWeights = [ ("cteq61",41) ]
# this analyzer finds the initial events before the skim
skimAnalyzer = cfg.Analyzer('skimAnalyzerCount')
eventSelector = cfg.Analyzer(
'EventSelector',
toSelect=[
# here put the event numbers (actual event numbers from CMSSW)
])
jsonAna = cfg.Analyzer('JSONAnalyzer', )
# this analyzer is just there to select a list of good primary vertices.
ttHVertexAna = cfg.Analyzer('VertexAnalyzer',
vertexWeight=None,
fixedWeight=1,
verbose=False)
# this analyzer actually does the pile-up reweighting.
def createTreeProducer(ana):
tp = cfg.Analyzer('_'.join(['TauTauLLTreeProducer', 'TauTauLLAnalyzer']),
anaName='TauTauLLAnalyzer')
return tp
tauScaleShift = 1.0
syncntuple = True
mc_vertexWeight = 'vertexWeightFall112011AB'
mc_tauEffWeight = None
mc_muEffWeight = None
mc_tauEffWeight_mc = 'effMediumIsoTau20MC'
mc_eleEffWeight_mc = 'effEle18MC'
mc_tauEffWeight = 'effTau2011AB_TauEle'
mc_eleEffWeight = 'effEle2011AB'
jsonAna = cfg.Analyzer(
'JSONAnalyzer',
)
triggerAna = cfg.Analyzer(
'TriggerAnalyzer',
verbose = False
)
vertexAna = cfg.Analyzer(
'VertexAnalyzer',
goodVertices = 'goodPVFilter',
vertexWeight = mc_vertexWeight,
# fixedWeight = 1,
verbose = False
)
return lepton.relIso(0.5)<0.1
def idIsoMuon(muon):
return idMuon(muon) and isoLepton(muon)
def passLepton(lepton):
return True
trigMap = { 'HLT_IsoMu15_v5':'hltSingleMuIsoL3IsoFiltered15',
'HLT_IsoMu15_v14':'hltSingleMuIsoL3IsoFiltered15' }
effMuAnaStd = cfg.Analyzer(
'EfficiencyAnalyzer_std',
# recsel = 'cuts_vbtfmuon',
recselFun = isoLepton,
# refselFun = idMuon,
triggerMap = trigMap,
instance = 'cmgMuonSelStdLep',
type = 'std::vector<cmg::Muon>',
genPdgId = 13
)
effMuAnaPF = cfg.Analyzer(
'EfficiencyAnalyzer_pf',
# recsel = 'cuts_vbtfmuon',
recselFun = isoLepton,
# refselFun = idMuon,
triggerMap = trigMap,
instance = 'cmgMuonSel',
type = 'std::vector<cmg::Muon>',
genPdgId = 13
)
# vertexFileDir = os.environ['CMSSW_BASE'] + '/src/CMGTools/RootTools/data/Reweight/2012/Vertices'
# vertexFileData = '/'.join([vertexFileDir, 'vertices_data_2012A_2012B_start_195947.root'])
mc_vertexWeight = None
# mc_tauEffWeight_mc = 'effTau2012MC53X'
# mc_muEffWeight_mc = 'eff_2012_Rebecca_TauMu_IsoMu1753XMC'
# mc_tauEffWeight = 'effTau2012ABC'
# mc_muEffWeight = 'effMu2012_Rebecca_TauMu_ABC'
mc_tauEffWeight_mc = 'effTau_muTau_MC_2012ABCDSummer13'
mc_muEffWeight_mc = 'effMu_muTau_MC_2012ABCD'
mc_tauEffWeight = 'effTau_muTau_Data_2012ABCDSummer13'
mc_muEffWeight = 'effMu_muTau_Data_2012ABCDSummer13'
eventSelector = cfg.Analyzer('EventSelector', toSelect=[105104])
jsonAna = cfg.Analyzer('JSONAnalyzer', )
triggerAna = cfg.Analyzer('TriggerAnalyzer')
vertexAna = cfg.Analyzer(
'VertexAnalyzer',
goodVertices='goodPVFilter',
vertexWeight=mc_vertexWeight,
fixedWeight=1,
verbose=False,
)
embedWeighter = cfg.Analyzer('EmbedWeighter', isRecHit=False, verbose=False)
mc_tauEffWeight = None
mc_tauEffWeight_mc = 'effLooseTau15MC'
if period == 'Period_2011A':
mc_vertexWeight = 'vertexWeightFall112invfb'
mc_tauEffWeight = 'effTau2011A'
elif period == 'Period_2011B':
mc_vertexWeight = 'vertexWeightFall112011B'
mc_tauEffWeight = 'effTau2011B'
elif period == 'Period_2011AB':
mc_vertexWeight = 'vertexWeightFall112011AB'
mc_tauEffWeight = 'effTau2011AB'
if runOnMC:
mc_tauEffWeight = 'effTau1fb'
triggerAna = cfg.Analyzer('TriggerAnalyzer')
TauTauAna = cfg.Analyzer(
'TauTauAnalyzer',
pt1=20,
pt2=20,
iso1=999,
iso2=999,
eta1=999,
eta2=999,
m_min=10,
m_max=99999,
diLeptonCutString='cuts_baseline',
# No trigger matching for the moment since MC only contains too high threshold trigger
#triggerMap = pathsAndFilters
)
#ttHLepAna.loose_muon_absIso5= 10
ttHLepAna.loose_electron_pt = 7
ttHLepAna.loose_electron_relIso = 0.4
#ttHLepAna.loose_electron_absIso = 10
ttHLepAna.ele_isoCorr = "rhoArea"
# --- LEPTON SKIMMING ---
ttHLepSkim.minLeptons = 1
ttHLepSkim.maxLeptons = 999
# --- JET-LEPTON CLEANING ---
ttHJetAna.minLepPt = 10
# Event Analyzer for susy multi-lepton (at the moment, it's the TTH one) // Do we need had W and Top?
ttHEventAna = cfg.Analyzer(
'ttHLepEventAnalyzer',
minJets25=0,
)
ttHIsoTrackAna.setOff = False
susyCoreSequence.insert(susyCoreSequence.index(ttHCoreEventAna), ttHSVAnalyzer)
susyCoreSequence.insert(susyCoreSequence.index(ttHCoreEventAna),
ttHHeavyFlavourHadronAnalyzer)
from CMGTools.TTHAnalysis.samples.samples_8TeV_v517 import triggers_1mu, triggers_1muHT, triggers_1eleHT # need to update the trigger MET pr HTMET?
# Tree Producer
treeProducer = cfg.Analyzer(
'treeProducerSusySingleSoftLepton',
vectorTree=True,
saveTLorentzVectors=
puFileMC = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/28-09-12/MC_Summer12_PU_S10-600bins.root'
puFileData = '/afs/cern.ch/user/a/agilbert/public/HTT_Pileup/28-09-12/Data_Pileup_2012_HCP-600bins.root'
vertexFileDir = os.environ[
'CMSSW_BASE'] + '/src/CMGTools/RootTools/data/Reweight/2012/Vertices'
vertexFileData = '/'.join(
[vertexFileDir, 'vertices_data_2012A_2012B_start_195947.root'])
mc_vertexWeight = None
mc_tauEffWeight_mc = 'effTau2012MC53X'
mc_muEffWeight_mc = 'eff_2012_Rebecca_TauMu_IsoMu1753XMC'
mc_tauEffWeight = 'effTau2012ABC'
mc_muEffWeight = 'effMu2012_Rebecca_TauMu_ABC'
eventSelector = cfg.Analyzer('EventSelector', toSelect=[105104])
jsonAna = cfg.Analyzer('JSONAnalyzer', )
triggerAna = cfg.Analyzer('TriggerAnalyzer')
vertexAna = cfg.Analyzer(
'VertexAnalyzer',
goodVertices='goodPVFilter',
vertexWeight=mc_vertexWeight,
fixedWeight=1,
verbose=False,
)
embedWeighter = cfg.Analyzer('EmbedWeighter', verbose=False)
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.RootTools.RootTools import *
PDFWeights = []
#PDFWeights = [ ("CT10",53), ("MSTW2008lo68cl",41), ("NNPDF21_100",101) ]
#PDFWeights = [ ("cteq61",41) ]
# this analyzer finds the initial events before the skim
skimAnalyzer = cfg.Analyzer(
'skimAnalyzerCount',
useLumiBlocks = False, # since we run on a non-CMG dataset, there's no info on lumi blocks, and this way we use less memory
)
eventSelector = cfg.Analyzer(
'EventSelector',
toSelect = [
# here put the event numbers (actual event numbers from CMSSW)
]
)
jsonAna = cfg.Analyzer(
'JSONAnalyzer',
)
triggerAna = cfg.Analyzer(
#'TriggerAnalyzer',
'triggerBitFilter',
)
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.RootTools.RootTools import *
#Load all common analyzers
from CMGTools.TTHAnalysis.analyzers.susyCore_modules_cff import *
from CMGTools.TTHAnalysis.samples.samples_8TeV_v517 import triggers_mumu, triggers_ee, triggers_mue, triggers_1mu
# Tree Producer
treeProducer = cfg.Analyzer(
'treeProducerSusyCore',
vectorTree=True,
PDFWeights=PDFWeights,
triggerBits={
'SingleMu': triggers_1mu,
'DoubleMu': triggers_mumu,
'DoubleEl': [t for t in triggers_ee if "Ele15_Ele8_Ele5" not in t],
'TripleEl': [t for t in triggers_ee if "Ele15_Ele8_Ele5" in t],
'MuEG': [t for t in triggers_mue if "Mu" in t and "Ele" in t]
})
#-------- SAMPLES AND TRIGGERS -----------
from CMGTools.TTHAnalysis.samples.samples_13TeV_CSA14 import *
selectedComponents = [
SingleMu, DoubleElectron, TTHToWW_PUS14, DYJetsToLL_M50_PU20bx25,
TTJets_PUS14
]
#-------- SEQUENCE