def makeDataComponent(self, name, dataset, user, pattern, json):
component = cfg.DataComponent(name=name,
files=self.getFiles(
dataset, user, pattern),
intLumi=1,
triggers=[],
json=json)
return component
def makeMCComponent(self, name, dataset, user, pattern):
component = cfg.MCComponent(
dataset=dataset,
name=name,
files=self.getFiles(dataset, user, pattern),
xSection=1,
nGenEvents=1,
triggers=[],
effCorrFactor=1,
# skimEfficiency = self.getSkimEfficiency(dataset,user)
)
# print 'Skim Efficiency for ',name,'=', component.skimEfficiency
return component
def makePrivateMCComponent(self, name, dataset, files):
if len(files) == 0:
raise RuntimeError, "Trying to make a component %s with no files" % name
# prefix filenames with dataset unless they start with "/"
dprefix = dataset + "/" if files[0][0] != "/" else ""
component = cfg.MCComponent(
dataset=dataset,
name=name,
files=['root://eoscms//eos/cms%s%s' % (dprefix, f) for f in files],
xSection=1,
nGenEvents=1,
triggers=[],
effCorrFactor=1,
# skimEfficiency = self.getSkimEfficiency(dataset,user)
)
# print 'Skim Efficiency for ',name,'=', component.skimEfficiency
return component
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
#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=
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',
# 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)
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
)
from CMGTools.H2TauTau.proto.samples.data_diTau import *
# from CMGTools.H2TauTau.proto.samples.embed import *
from CMGTools.H2TauTau.proto.samples.ewk import *
import copy
import re
import CMGTools.RootTools.fwlite.Config as cfg
#from CMGTools.RootTools.yellowreport.YRParser import yrparser
Higgsgg110 = cfg.MCComponent(name='Higgsgg110',
files=[],
xSection=19.84 * 6.37e-2,
nGenEvents=199988,
triggers=[],
effCorrFactor=1)
Higgsgg115 = cfg.MCComponent(name='Higgsgg115',
files=[],
xSection=18.13 * 6.37e-2,
nGenEvents=196002,
triggers=[],
effCorrFactor=1)
Higgsgg120 = cfg.MCComponent(name='Higgsgg120',
files=[],
xSection=16.63 * 6.37e-2,
nGenEvents=199981,
triggers=[],
effCorrFactor=1)
Higgsgg125 = cfg.MCComponent(name='Higgsgg125',
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)
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
)
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
)
def createTreeProducer(ana):
tp = cfg.Analyzer('_'.join(['TauTauLLTreeProducer', 'TauTauLLAnalyzer']),
anaName='TauTauLLAnalyzer')
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),
])
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)
##########################################################
## 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.
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 = None
)
jetAna = cfg.Analyzer(
'PFSimJetAnalyzer',
jetPt = 10.,
jetEta = 5.0,
btagSFseed = 123456,
relaxJetId = True,
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.H2TauTau.proto.samples.getFiles import getFiles
DYJets = cfg.MCComponent(name='DYJets',
files=[],
xSection=3048.,
nGenEvents=1,
triggers=[],
effCorrFactor=1)
WJets = cfg.MCComponent(name='WJets',
files=[],
xSection=31314.,
nGenEvents=1,
triggers=[],
effCorrFactor=1)
TTJets = cfg.MCComponent(name='TTJets',
files=[],
xSection=165.8,
nGenEvents=1,
triggers=[],
effCorrFactor=1)
mc_basicEWK = [DYJets, WJets, TTJets]
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 copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.RootTools.fwlite.Config import printComps
from CMGTools.RootTools.RootTools import *
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 0.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
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,
'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'
import copy
import os
import CMGTools.RootTools.fwlite.Config as cfg
from CMGTools.H2TauTau.proto.samples.getFiles import getFiles
jsondir = '/afs/cern.ch/cms/CAF/CMSCOMM/COMM_DQM/certification/Collisions11/7TeV'
promptjson = '/'.join([
jsondir, 'Prompt/Cert_160404-180252_7TeV_PromptReco_Collisions11_JSON.txt'
])
embed_Run2011A_May10ReReco_v1 = cfg.EmbedComponent(
name='embed_Run2011A_May10ReReco_v1',
files=[],
triggers=[],
json='/'.join([
jsondir,
'Reprocessing/Cert_160404-163869_7TeV_May10ReReco_Collisions11_JSON_v3.txt'
]))
embed_Run2011A_PromptReco_v4 = cfg.EmbedComponent(
name='embed_Run2011A_PromptReco_v4',
files=[],
triggers=[],
json=promptjson)
embed_Run2011A_05Aug2011_v1 = cfg.EmbedComponent(
name='embed_Run2011A_05Aug2011_v1',
files=[],
triggers=[],
json='/'.join([
jsondir,
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,
import CMGTools.RootTools.fwlite.Config as cfg
import copy
# the luminosities for these components correspond to the
# HLT triggers used in the H->tau tau analysis
# they might differ from the luminosities you want for a
# Z-> mu mu or Z-> ee analysis.
data_Run2012A = cfg.DataComponent(
name='data_Run2012A',
files=[],
intLumi=890., # 806.194 from Colin
triggers=[],
json=None)
data_Run2012B = cfg.DataComponent(name='data_Run2012B',
files=[],
intLumi=4429.,
triggers=[],
json=None)
data_Run2012C = cfg.DataComponent(
name='data_Run2012C',
files=[],
# | HLT_IsoMu17_eta2p1_LooseIsoPFTau20_v3 | 6719 | 495.003(/pb) | 495.003(/pb) |
intLumi=7153.,
triggers=[],
json=None)
data_Run2012D = cfg.DataComponent(
name='data_Run2012D',
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 CMGTools.RootTools.fwlite.Config as cfg
# from CMGTools.H2TauTau.proto.samples.sampleShift import sampleShift
DYJets = cfg.MCComponent(
name = 'DYJets',
files = [],
xSection = 2008.4*3.,
nGenEvents = 1,
triggers = [],
effCorrFactor = 1 )
DY1Jets = cfg.MCComponent(
name = 'DY1Jets',
files = [],
xSection = 561.0,
nGenEvents = 1,
triggers = [],
effCorrFactor = 1 )
DY2Jets = cfg.MCComponent(
name = 'DY2Jets',
files = [],
xSection = 181.,
nGenEvents = 1,
triggers = [],
effCorrFactor = 1 )
DY3Jets = cfg.MCComponent(
name = 'DY3Jets',
files = [],
xSection = 51.1,
def createTreeProducer(ana):
tp = cfg.Analyzer('_'.join(['FourJetTreeProducer', 'FourJetAnalyzer']),
anaName='FourJetAnalyzer')
return tp
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 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},