# Mjj = 700,
# deltaEta = 4.0,
cjvPtCut=30.)
treeProducer = cfg.Analyzer('GenTreeProducer')
multiProducer = cfg.Analyzer('MultiTreeProducer')
###############################################################################
testComp = cfg.MCComponent(
name='Test',
# files = ['gen.root'],
files=['/data1/ytakahas/data/gen_atlas.root'],
# files = ['/data1/ytakahas/data/gen_cms.root'],
xSection=1,
nGenEvents=1,
triggers=[],
effCorrFactor=1)
testComp.splitFactor = 1
#testComp.splitFactor = 6
selectedComponents = [testComp]
sequence = cfg.Sequence([reader, jetAna, vbfAna, treeProducer, multiProducer])
config = cfg.Config(components=selectedComponents, sequence=sequence)
printComps(config.components, True)
data_2011A = [
data_Run2011A_May10ReReco_v1,
data_Run2011A_PromptReco_v4,
data_Run2011A_05Aug2011_v1,
data_Run2011A_03Oct2011_v1,
]
data_2011B = [data_Run2011B_PromptReco_v1]
selectedComponents = MC
if period == 'Period_2011A':
selectedComponents.extend(data_2011A)
elif period == 'Period_2011B':
selectedComponents.extend(data_2011B)
elif period == 'Period_2011AB':
selectedComponents.extend(data_2011A)
selectedComponents.extend(data_2011B)
# selectedComponents = data_2011A
# selectedComponents = [embed_Run2011A_PromptReco_v4]
selectedComponents = [DYJets]
sequence = cfg.Sequence([
# triggerAna,
# vertexAna,
# ZMuMuAna,
effMuAna
])
config = cfg.Config(components=selectedComponents, sequence=sequence)
JetMonABCD.splitFactor = 500
selectedComponents = mcSamples + dataSamples1Mu + dataSamplesMu + dataSamplesE + dataSamplesMuE + [
JetMonABCD
]
SingleMuD.splitFactor = 1200
SingleMuC.splitFactor = 1200
QCDMuPt15.splitFactor = 500
#-------- SEQUENCE
sequence = cfg.Sequence([
skimAnalyzer,
jsonAna,
triggerAna,
ttHGenAna,
ttHVertexAna,
ttHLepAna,
ttHLepMCAna,
ttHJetAna,
ttHJetMCAna,
treeProducer,
])
#-------- HOW TO RUN
# set test = 0 to run all jobs, in case you are using pybatch.py
selectedComponents.remove(DoubleMuD)
selectedComponents.remove(DYJetsM50)
selectedComponents.remove(WJets)
#selectedComponents.remove(TTJets)
test = 1
selectedComponents = allsamples
#selectedComponents = [HiggsGGH125]
#HiggsGGH125.splitFactor=10
#HiggsGGH125.files = getFiles('/GluGluToHToTauTau_M-125_8TeV-powheg-pythia6/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM/V5_B/PAT_CMG_V5_8_0/HTT_Nom_Mar6_hinzmann', 'hinzmann', 'diTau.*fullsel.*root', True)
#HiggsGGH125.files = getFiles('/GluGluToHToTauTau_M-125_8TeV-powheg-pythia6/Summer12_DR53X-PU_S10_START53_V7A-v1/AODSIM/V5_B/PAT_CMG_V5_8_0/HTT_Mar7_manzoni_Nom', 'manzoni', 'diTau.*fullsel.*root', True)
print [c.name for c in selectedComponents]
sequence = cfg.Sequence( [
triggerAna,
vertexAna,
TauTauAna,
vbfAna,
pileUpAna,
embedWeighter,
tau1Weighter,
tau2Weighter,
jetWeighter,
treeProducer,
treeProducerXcheck,
] )
if runOnEmbedded:
sequence.remove(triggerAna)
if use1prong or useParked:
sequence.remove(jetWeighter)
if use1prong or useParked:
for data in data_2012A:
triggers = [],
effCorrFactor = 1 )
WWH = cfg.MCComponent(
name = 'WWH',
files = ['wwh_gensim.root'],
xSection = 0.0267,
nGenEvents = 5000,
triggers = [],
effCorrFactor = 1 )
selectedComponents = [ZH]
sequence = cfg.Sequence( [
# effAndSmear,
reader,
jetAna,
missMassAna,
treeProducer
] )
# creation of the processing configuration.
# we define here on which components to run, and
# what is the sequence of analyzers to run on each event.
config = cfg.Config( components = selectedComponents,
sequence = sequence )
printComps(config.components, True)
selectedComponents.extend( data_2011A )
selectedComponents.extend( embed_2011A )
elif period == 'Period_2011B':
selectedComponents.extend( data_2011B )
selectedComponents.extend( embed_2011B )
elif period == 'Period_2011AB':
selectedComponents.extend( data_2011A )
selectedComponents.extend( data_2011B )
selectedComponents.extend( embed_2011A )
selectedComponents.extend( embed_2011B )
sequence = cfg.Sequence( [
triggerAna,
vertexAna,
tauEleAna,
tauWeighter,
eleWeighter,
vbfAna,
eventSorter
] )
DYJets.fakes = True
TTJets.splitFactor = 2
DYJets.splitFactor = 8
WJets.splitFactor = 3
data_Run2011B_PromptReco_v1.splitFactor = 8
data_Run2011A_PromptReco_v4.splitFactor = 2
embed_Run2011A_May10ReReco_v1.splitFactor = 2
embed_Run2011A_PromptReco_v4.splitFactor = 4
embed_Run2011A_05Aug2011_v1.splitFactor = 2
pfmetcut=0,
jetptcut=1000,
pt=30,
eta=2.1,
iso=0.5,
savegenp=False,
verbose=True,
triggerBits={'SingleMu': triggers_mu},
keepFailingEvents=False)
ZtreeProducer = cfg.Analyzer('ZTreeProducer')
genAna = cfg.Analyzer('GenParticleAnalyzerFSR', src='genParticlesPruned')
sequence = cfg.Sequence([
genAna, jsonAna, triggerAna, vertexAna, WAna, WtreeProducer, ZAna,
ZtreeProducer
])
from CMGTools.H2TauTau.proto.samples.ewk import TTJets
from CMGTools.H2TauTau.proto.samples.getFiles import getFiles
TTJets.files = getFiles(
'/TTJets_TuneZ2_7TeV-madgraph-tauola/Summer11LegDR-PU_S13_START53_LV6-v1/AODSIM/V5_B/PAT_CMG_V5_18_0',
'cmgtools', '.*root')
TTJets.triggers = triggers_mu
# TTJets.triggers = ["HLT_IsoMu24_v1","HLT_IsoMu24_v2","HLT_IsoMu24_v3","HLT_IsoMu24_v4","HLT_IsoMu24_v5","HLT_IsoMu24_v6","HLT_IsoMu24_v7",\
# "HLT_IsoMu24_v8","HLT_IsoMu24_v9","HLT_IsoMu24_v10","HLT_IsoMu24_v11","HLT_IsoMu24_v12","HLT_IsoMu24_v13","HLT_IsoMu24_v14",\
# "HLT_IsoMu24_eta2p1_v1","HLT_IsoMu24_eta2p1_v2","HLT_IsoMu24_eta2p1_v3","HLT_IsoMu24_eta2p1_v4","HLT_IsoMu24_eta2p1_v5",\
# "HLT_IsoMu24_eta2p1_v6","HLT_IsoMu24_eta2p1_v7","HLT_IsoMu24_eta2p1_v8"
# ]
selectedComponents = [sam for sam in allsamples if len(sam.files)>0]
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'
for mc in MC_list:
mc.puFileMC = puFileMC
mc.puFileData = puFileData
sequence = cfg.Sequence( [
# eventSelector,
jsonAna,
# triggerAna,
vertexAna,
pileUpAna,
ZMuMuAna,
# ZMuMuRochCor,
jetAna,
# vbfAna,
treeProducer
] )
# set test = 0 to run all jobs, in case you are using pybatch.py
test = 1
if test==1:
# test a single component, using a single thread.
# necessary to debug the code, until it doesn't crash anymore
comp = DYJets
comp.files = comp.files[:10]
selectedComponents = [comp]
triggers=[],
effCorrFactor=1)
DYJetsChamonix_PAT = cfg.MCComponent(
name='DYJetsChamonix',
files=getFiles(
'/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola/Fall11-PU_Chamonix12_START44_V10-v2/AODSIM/PAT_CMG_V3_0_0',
'cmgtools', 'pat.*root')[:50],
# files = createDataset('LOCAL','/DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola/Fall11-PU_Chamonix12_START44_V10-v2/AODSIM/PAT_CMG_V3_0_0', '.*root', True).listOfGoodFiles(),
xSection=3048.,
nGenEvents=34915945,
triggers=[],
effCorrFactor=1)
#########################################################################################
MC = [DYJetsFall11_PAT, DYJetsChamonix_PAT]
for mc in MC:
# could handle the weights in the same way
mc.jetScale = mc_jet_scale
mc.jetSmear = mc_jet_smear
DYJetsChamonix_PAT.splitFactor = 50
DYJetsFall11_PAT.splitFactor = 50
selectedComponents = [DYJetsFall11_PAT]
sequence = cfg.Sequence([jetAna])
config = cfg.Config(components=selectedComponents, sequence=sequence)
'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
sequence = cfg.Sequence(susyCoreSequence + [
treeProducer,
])
#-------- HOW TO RUN
test = 1
if test == 1:
# test a single component, using a single thread.
comp = TTJets_PU20bx25
comp.files = comp.files[:1]
selectedComponents = [comp]
comp.splitFactor = 1
elif test == 2:
# test all components (1 thread per component).
for comp in selectedComponents:
comp.splitFactor = 1
comp.files = comp.files[:1]
selectedComponents = [WJetsSoup, TTJets, DYJets]
selectedComponents = [WJets, W1Jets, W2Jets, W3Jets, W4Jets, TTJets, DYJets]
higgs = mc_higgs
selectedComponents.extend(higgs)
selectedComponents.extend(diboson_list)
selectedComponents.extend(data_list_2011)
selectedComponents.extend(embed_list_2011)
sequence = cfg.Sequence([
# eventSelector,
jsonAna,
triggerAna,
vertexAna,
TauMuAna,
dyJetsFakeAna,
WNJetsAna,
higgsWeighter,
vbfAna,
embedWeighter,
tauWeighter,
muonWeighter,
treeProducer
])
if syncntuple:
sequence.append(treeProducerXCheck)
test = 1
if test == 1:
comp = HiggsVBF125
# comp.files = comp.files[:2]
dataSequence=[
jsonFilter,
puAna,
triggerAna,
vertexAna,
fakeAna,
fakeTree
]
#for full production use data sequence
sequence = cfg.Sequence(dataSequence)
test = 0
if test==1:
dataset = DoubleMu
selectedComponents = [dataset]
dataset.splitFactor = 1
if test ==2:
#########################################################################################
for mc in MC_list:
mc.puFileMC = puFileMC
mc.puFileData = puFileData
for emb in embed_list:
emb.puFileData = None
emb.puFileMC = None
selectedComponents = allsamples
sequence = cfg.Sequence([
eventSelector, jsonAna, triggerAna, vertexAna, TauEleAna, dyJetsFakeAna,
WNJetsAna, higgsWeighter, vbfAna, pileUpAna, embedWeighter, tauWeighter,
electronWeighter, treeProducer
])
if syncntuple:
sequence.append(treeProducerXCheck)
test = 1
if test == 1:
# comp = DYJets
# comp = data_Run2012A
# comp = data_Run2012C_v2
comp = HiggsVBF125
# comp = HiggsVH125
# comp = Tbar_tW
# HiggsVBF125.triggers = []
# mc.triggers = triggers_HT650+triggers_MET150+triggers_HTMET
#for data in dataSamples:
# data.triggers = triggers_HT
mcSamples_mainFullHad = [
T2qq_2J, T2bb_2J, TTJetsSem, TTJetsLep, ZNuNu50HT100, ZNuNu100HT200,
ZNuNu200HT400, ZNuNu400
]
selectedComponents = mcSamples_mainFullHad
#-------- SEQUENCE
sequence = cfg.Sequence(susyCoreSequence + [
ttHIsoTrackAna,
ttHTopoJetAna,
treeProducer,
])
#-------- HOW TO RUN
test = 1
if test == 1:
# test a single component, using a single thread.
comp = ZNuNu50HT100
## comp = T2qq_2J_1
## comp = TTJetsSem
## comp = TTJetslep
comp.files = comp.files[:1]
selectedComponents = [comp]
comp.splitFactor = 1
if runOnData:
selectedComponents = data_2011
if runOnMC:
selectedComponents = [DYJets, WJets, TTJets, WW, WZ, ZZ]
selectedComponents += [
Higgsgg110, Higgsgg115, Higgsgg120, Higgsgg125, Higgsgg130, Higgsgg135,
Higgsgg140, Higgsgg145, HiggsVBF110, HiggsVBF115, HiggsVBF120,
HiggsVBF125, HiggsVBF130, HiggsVBF135, HiggsVBF140, HiggsVBF145
]
print[c.name for c in selectedComponents]
if runOnMC:
sequence = cfg.Sequence([
TauTauAna, vbfAna, vertexAna, tau1Weighter, tau2Weighter, treeProducer
])
else:
sequence = cfg.Sequence([
triggerAna, TauTauAna, vbfAna, vertexAna, tau1Weighter, tau2Weighter,
treeProducer
])
DYJets.fakes = True
DYJets.splitFactor = 50
WJets.fakes = True
WJets.splitFactor = 50
TTJets.splitFactor = 50
WW.splitFactor = 10
WZ.splitFactor = 10
ZZ.splitFactor = 10
from PhysicsTools.HeppyCore.framework.heppy import getHeppyOption
#Load all common analyzers
from CMGTools.VVResonances.analyzers.core_cff import *
#-------- SAMPLES AND TRIGGERS -----------
from CMGTools.VVResonances.samples.samples_13TeV_Spring15 import *
selectedComponents = [RSGravToWWToLNQQ_2000]
#-------- Analyzer
from CMGTools.VVResonances.analyzers.tree_cff import *
#-------- SEQUENCE
sequence = cfg.Sequence(coreSequence + [leptonTreeProducer])
#-------- HOW TO RUN
test = 1
if test == 1:
# test a single component, using a single thread.
comp = RSGravToWWToLNQQ_2000
comp.files = comp.files[:1]
selectedComponents = [comp]
comp.splitFactor = 1
elif test == 2:
# test all components (1 thread per component).
for comp in selectedComponents:
comp.splitFactor = 1
comp.files = comp.files[:1]
WJets,
data_Run2012A,
data_Run2012B,
data_Run2012C_v1,
data_Run2012C_v2,
]
sequence = cfg.Sequence([
# eventSelector,
jsonAna,
triggerAna, #No!
vertexAna,
TauMuAna, #Yes!
dyJetsFakeAna,
#JAN contact Colin when you want to deal with exclusive samples
# WNJetsAna,
higgsWeighter,
vbfAna, #Yes!
pileUpAna,
embedWeighter,
tauWeighter,
muonWeighter,
treeProducer, #Yes!
# treeProducerXCheck
])
if syncntuple:
sequence.append(treeProducerXCheck) #Yes!
test = 1
if test == 1:
comp = HiggsVBF125
particleAna = cfg.Analyzer('ParticleAnalyzer')
treeProducer = cfg.Analyzer('PFSimTreeProducer')
###############################################################################
testComp = cfg.MCComponent(name='Test',
files=glob.glob('Prod_211/Job*/sim.root'),
xSection=1,
nGenEvents=1,
triggers=[],
effCorrFactor=1)
testComp.splitFactor = len(testComp.files)
selectedComponents = [testComp]
sequence = cfg.Sequence([
# effAndSmear,
reader,
jetAna,
particleAna,
treeProducer
])
# creation of the processing configuration.
# we define here on which components to run, and
# what is the sequence of analyzers to run on each event.
config = cfg.Config(components=selectedComponents, sequence=sequence)
printComps(config.components, True)
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
stopper = cfg.Analyzer('Stopper')
sequence = cfg.Sequence([
fourJetAna,
createTreeProducer(fourJetAna),
])
Hig125 = cfg.MCComponent(
name='Hig125',
files=[
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_80.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_81.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_82.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_83.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_84.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_85.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_86.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_87.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_88.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_89.root',
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
comp.splitFactor = 14
config = cfg.Config(components=selectedComponents, sequence=sequence)
printComps(config.components, True)
)
#-------- SAMPLES
from CMGTools.TTHAnalysis.samples.samples_13TeV_CSA14 import *
selectedComponents = [TTJets_PU20bx25]
#-------- SEQUENCE
sequence = cfg.Sequence([
skimAnalyzer,
#eventSelector,
jsonAna,
pileUpAna,
ttHGenAna,
ttHVertexAna,
ttHLepAna,
ttHLepMCAna,
ttHJetAna,
ttHJetMCAna,
treeProducer,
])
#-------- HOW TO RUN
test = 1
if test == 1:
# test a single component, using a single thread.
# necessary to debug the code, until it doesn't crash anymore
comp = TTJets_PU20bx25
comp.files = comp.files[:1]
selectedComponents = [comp]
# selectedComponents = copy.copy(data_list_2012)
selectedComponents = data_list_2012
# selectedComponents.append( DYJets )
# selectedComponents = []
selectedComponents.extend([TTJets, DYJets])
selectedComponents.extend(mc_diboson)
VVgroup = [comp.name for comp in mc_diboson]
sequence = cfg.Sequence([
# eventSelector,
jsonAna,
vertexAna,
ZMuMuAna,
leg1Weighter,
leg2Weighter,
# WNJetsAna,
vbfAna,
treeProducer
])
test = 0
if test == 1:
comp = DYJets
comp.files = comp.files[:10]
selectedComponents = [comp]
comp.splitFactor = 1
elif test == 2:
selectedComponents = copy.copy(data_list_2012)
selectedComponents.extend(embed_list_2012)
selectedComponents = [WJetsSoup, TTJets, DYJets]
# selectedComponents = [WJets, W1Jets, W2Jets, W3Jets, W4Jets, TTJets, DYJets]
higgs = mc_higgs
selectedComponents.extend( higgs )
selectedComponents.extend( diboson_list )
selectedComponents.extend( data_list_2011 )
selectedComponents.extend( embed_list_2011 )
sequence = cfg.Sequence( [
# eventSelector,
jsonAna,
triggerAna,
vertexAna,
tauEleAna,
dyJetsFakeAna,
dyLLReweighterTauEle,
WNJetsAna,
higgsWeighter,
vbfAna,
embedWeighter,
tauWeighter,
eleWeighter,
treeProducer
] )
if syncntuple:
sequence.append( treeProducerXCheck)
test = 1
if test==1:
comp = HiggsVBF125
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')
sequence = cfg.Sequence([
httAna,
createTreeProducer(httAna),
])
Hig125 = cfg.MCComponent(
name='Hig125',
zzfiles=[
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_80.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_81.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_82.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_83.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_84.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_85.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_86.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_87.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_88.root',
'root://eoscms//eos/cms/store/cmst3/user/pjanot/LEP3/ZZ/cmgTuple_ZZ_89.root',
selectedComponents.extend(higgs)
selectedComponents.extend(mc_higgs_susy)
if not simulatedOnly:
selectedComponents.extend(embed_list)
sequence = cfg.Sequence([
# eventSelector,
jsonAna,
triggerAna,
vertexAna,
TauMuAna,
dyJetsFakeAna,
# WNJetsAna,
# WNJetsTreeAna,
NJetsAna,
higgsWeighter,
jetAna,
vbfSimpleAna,
pileUpAna,
embedWeighter,
tauDecayModeWeighter,
tauWeighter,
muonWeighter,
treeProducer,
# treeProducerXCheck
])
if syncntuple:
sequence.append(treeProducerXCheck) #Yes!
selectedComponents = [
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),
])
# 'MuHT' : triggers_1muHT,
# 'eleHT' : triggers_1eleHT
# }
)
#-------- SAMPLES AND TRIGGERS -----------
from CMGTools.TTHAnalysis.samples.samples_13TeV_CSA14 import *
selectedComponents = [
SingleMu, DoubleElectron, TTHToWW_PUS14, DYJetsToLL_M50_PU20bx25,
TTJets_PUS14
]
#-------- SEQUENCE
sequence = cfg.Sequence(susyCoreSequence + [
ttHEventAna,
treeProducer,
])
#-------- HOW TO RUN
test = 1
if test == 1:
# test a single component, using a single thread.
comp = TTJets_PUS14
#comp.files = ['root://eoscms//eos/cms/store/cmst3/group/susy/alobanov/MC/MiniAOD/13TeV_Gl_Gl_4q_Gl1400_LSP300_Chi315_MiniAOD.root']
# comp.files = comp.files[:1]
selectedComponents = [comp]
comp.splitFactor = 1
elif test == 2:
# test all components (1 thread per component).
for comp in selectedComponents:
comp.splitFactor = len(comp.files)
#########################################################################################
dummyAna = cfg.Analyzer(
'Analyzer'
)
selectedComponents = [DYJetsFall11, QCDMu]
DYJetsChamonix.splitFactor = splitFactor
DYJetsFall11.splitFactor = splitFactor
QCDMu.splitFactor = splitFactor
# QCDMu.files = QCDMu.files[:5]
test = False
if test:
sam = DYJetsFall11
sam.files = sam.files[:1]
selectedComponents = [sam]
sam.splitFactor = 1
sequence = cfg.Sequence( [jetAna, jetAnaU] )
config = cfg.Config( components = selectedComponents,
sequence = sequence )
from CMGTools.TTHAnalysis.samples.samples_13TeV_CSA14 import *
selectedComponents = [ SingleMu, DoubleElectron, TTHToWW_PUS14, DYJetsToLL_M50_PU20bx25, TTJets_PUS14 ]
#-------- SEQUENCE
sequence = cfg.Sequence([
skimAnalyzer,
#eventSelector,
jsonAna,
triggerAna,
pileUpAna,
ttHGenAna,
susyScanAna,
ttHVertexAna,
ttHLepAna,
ttHLepMCAna,
ttHTauAna,
ttHTauMCAna,
ttHJetAna,
ttHJetMCAna,
ttHCoreEventAna,
ttHEventAna,
treeProducer,
])
#-------- HOW TO RUN
test = 1
if test==1:
# test a single component, using a single thread.
from PhysicsTools.HeppyCore.framework.heppy import getHeppyOption
#Load all common analyzers
from CMGTools.VVResonances.analyzers.core_cff import *
#-------- SAMPLES AND TRIGGERS -----------
from CMGTools.VVResonances.samples.samples_13TeV_Spring15 import *
selectedComponents = mcSamples
#-------- Analyzer
from CMGTools.VVResonances.analyzers.tree_cff import *
#-------- SEQUENCE
sequence = cfg.Sequence(coreSequence + [vvSkimmer, vvTreeProducer])
#-------- HOW TO RUN
test = 0
if test == 1:
# test a single component, using a single thread.
comp = RSGravToWWToLNQQ_2000
comp.files = comp.files[:1]
selectedComponents = [comp]
comp.splitFactor = 1
elif test == 2:
# test all components (1 thread per component).
for comp in selectedComponents:
comp.splitFactor = 1
comp.files = comp.files[:1]