from CMGTools.HToZZ4L.analyzers.hzz4lCore_modules_cff import *
#-------- SAMPLES AND TRIGGERS -----------
#-------- SEQUENCE
from CMGTools.HToZZ4L.samples.samples_13TeV_Spring15 import *
selectedComponents = mcSamples + dataSamples
sequence = cfg.Sequence(hzz4lCoreSequence)
for comp in mcSamples:
comp.triggers = []
comp.vetoTriggers = []
from PhysicsTools.HeppyCore.framework.heppy import getHeppyOption
test = getHeppyOption('test')
if test == "1":
comp = QQHZZ4L
comp.files = comp.files[:1]
comp.splitFactor = 1
comp.fineSplitFactor = 1 if getHeppyOption('single') else 5
selectedComponents = [comp]
if getHeppyOption('events'):
eventSelector.toSelect = [
eval("(" + x.replace(":", ",") + ")")
for x in getHeppyOption('events').split(",")
]
sequence = cfg.Sequence([eventSelector] + hzz4lCoreSequence)
print "Will select events ", eventSelector.toSelect
elif test == '2':
for comp in selectedComponents:
#comp.files = comp.files[:1]
#comp.files = comp.files[57:58] # to process only file [57]
#comp.triggers = triggers_HT900 + triggers_HTMET # to apply trigger skimming
# ------------------------------------------------------------------------------------------- #
from PhysicsTools.HeppyCore.framework.services.tfile import TFileService
output_service = cfg.Service(
TFileService,
'outputfile',
name="outputfile",
fname='mt2.root',
option='recreate'
)
# the following is declared in case this cfg is used in input to the heppy.py script
from PhysicsTools.HeppyCore.framework.eventsfwlite import Events
from CMGTools.TTHAnalysis.tools.EOSEventsWithDownload import EOSEventsWithDownload
event_class = EOSEventsWithDownload
if getHeppyOption("nofetch"):
event_class = Events
config = cfg.Config( components = selectedComponents,
sequence = sequence,
services = [output_service],
events_class = event_class)
# events_class = Events)
#printComps(config.components, True)
#for comp in selectedComponents:
# comp.splitFactor = 1
# comp.fineSplitFactor = 40
#-------- HOW TO RUN -----------
from PhysicsTools.HeppyCore.framework.heppy import getHeppyOption
#test = getHeppyOption('test')
test = 'SR'
if test == '1':
comp = Monojet_M_1000_AV
monoJetSkim.metCut = 200
comp.files = comp.files[:1]
#comp.files = [ '/afs/cern.ch/work/a/avartak/public/dmVM10.root' ]
comp.splitFactor = 1
comp.fineSplitFactor = 1
if not getHeppyOption('single'):
comp.fineSplitFactor = 4
selectedComponents = [ comp ]
elif test == '2':
for comp in selectedComponents:
comp.files = comp.files[:1]
comp.splitFactor = 1
comp.fineSplitFactor = 1
elif test == 'EOS':
comp = DYJetsToLL_M50#TTJets
comp.files = comp.files[:1]
if getHeppyOption('Wigner'):
print "Will read from WIGNER"
comp.files = [ 'root://eoscms//eos/cms/store/mc/Phys14DR/DYJetsToLL_M-50_13TeV-madgraph-pythia8/MINIAODSIM/PU20bx25_PHYS14_25_V1-v1/00000/0432E62A-7A6C-E411-87BB-002590DB92A8.root' ]
else:
print "Will read from CERN Meyrin"
ttHFatJetAna,
ttHAlphaTAna,
ttHAlphaTControlAna,
monoJetVarAna,
MonoJetEventAna,
treeProducer,
])
# -- fine splitting, for some private MC samples with a single file
#for comp in selectedComponents:
# comp.splitFactor = 1
# comp.fineSplitFactor = 40
#-------- HOW TO RUN -----------
from PhysicsTools.HeppyCore.framework.heppy import getHeppyOption
test = getHeppyOption('test')
if test == '1':
monoJetSkim.metCut = 0
comp = DYJetsToLL_M50_HT100to200
comp.files = comp.files[:1]
comp.splitFactor = 1
comp.fineSplitFactor = 1
selectedComponents = [comp]
elif test == '2':
for comp in selectedComponents:
comp.files = comp.files[:1]
comp.splitFactor = 1
comp.fineSplitFactor = 1
elif test == 'EOS':
comp = DYJetsToLL_M50 #TTJets
comp.files = comp.files[:1]
comp.files = comp.files[:1]
## output histogram
outputService=[]
from PhysicsTools.HeppyCore.framework.services.tfile import TFileService
output_service = cfg.Service(
TFileService,
'outputfile',
name="outputfile",
fname='vvTreeProducer/tree.root',
option='recreate'
)
outputService.append(output_service)
from PhysicsTools.HeppyCore.framework.eventsfwlite import Events
from CMGTools.TTHAnalysis.tools.EOSEventsWithDownload import EOSEventsWithDownload
event_class = EOSEventsWithDownload
if getHeppyOption("nofetch"):
event_class = Events
config = cfg.Config( components = selectedComponents,
sequence = sequence,
services = [],
events_class = event_class)
from CMGTools.H2TauTau.proto.analyzers.TauTreeProducer import TauTreeProducer
from CMGTools.H2TauTau.proto.analyzers.TauGenTreeProducer import TauGenTreeProducer
# from CMGTools.H2TauTau.proto.samples.phys14.connector import httConnector
from CMGTools.TTHAnalysis.samples.samples_13TeV_PHYS14 import DYJetsToLL_M50, TTH
from CMGTools.TTHAnalysis.samples.ComponentCreator import ComponentCreator
# common configuration and sequence
from CMGTools.H2TauTau.htt_ntuple_base_cff import genAna, vertexAna
# Get all heppy options; set via "-o production" or "-o production=True"
# production = True run on batch, production = False (or unset) run locally
production = getHeppyOption('production')
treeProducer = cfg.Analyzer(
TauTreeProducer,
name='TauTreeProducer'
)
genTreeProducer = cfg.Analyzer(
TauGenTreeProducer,
name='TauGenTreeProducer'
)
tauAna = cfg.Analyzer(
TauAnalyzer,
name='TauAnalyzer',
#for comp in selectedComponents:
# comp.splitFactor = 1
# comp.fineSplitFactor = 40
#-------- SEQUENCE -----------
sequence = cfg.Sequence(susyCoreSequence + [
ttHJetTauAna,
ttHEventAna,
treeProducer,
])
#-------- HOW TO RUN -----------
from PhysicsTools.HeppyCore.framework.heppy import getHeppyOption
test = getHeppyOption('test')
if test == '1':
comp = TTH
if getHeppyOption('T1tttt'):
comp = SMS_T1tttt_2J_mGl1500_mLSP100
if getHeppyOption('H4L'):
comp = GGHZZ4L
comp.files = comp.files[:1]
comp.splitFactor = 1
if not getHeppyOption('single'):
comp.fineSplitFactor = 4
selectedComponents = [comp]
elif test == '2':
for comp in selectedComponents:
comp.files = comp.files[:1]
comp.splitFactor = 1
from CMGTools.H2TauTau.proto.analyzers.TauMuAnalyzer import TauMuAnalyzer
from CMGTools.H2TauTau.proto.analyzers.H2TauTauTreeProducerTauMu import H2TauTauTreeProducerTauMu
from CMGTools.H2TauTau.proto.analyzers.TauDecayModeWeighter import TauDecayModeWeighter
from CMGTools.H2TauTau.proto.analyzers.TauFakeRateWeighter import TauFakeRateWeighter
from CMGTools.H2TauTau.proto.analyzers.LeptonWeighter import LeptonWeighter
from CMGTools.H2TauTau.proto.analyzers.SVfitProducer import SVfitProducer
from CMGTools.H2TauTau.proto.samples.phys14.connector import httConnector
# common configuration and sequence
from CMGTools.H2TauTau.htt_ntuple_base_cff import commonSequence, genAna, dyJetsFakeAna, puFileData, puFileMC, eventSelector
# Get all heppy options; set via "-o production" or "-o production=True"
# production = True run on batch, production = False (or unset) run locally
production = getHeppyOption('production')
production = False
# mu-tau specific configuration settings
# 'Nom', 'Up', 'Down', or None
shift = None
syncntuple = True
computeSVfit = False
# When ready, include weights from CMGTools.H2TauTau.proto.weights.weighttable
# 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'
#-------- SAMPLES AND TRIGGERS -----------
#-------- SEQUENCE
from CMGTools.HToZZ4L.samples.samples_13TeV_PHYS14 import *
selectedComponents = []
sequence = cfg.Sequence(hzz4lCoreSequence)
for comp in mcSamples:
comp.triggers = triggers_multilep
comp.vetoTriggers = []
from PhysicsTools.HeppyCore.framework.heppy import getHeppyOption
test = getHeppyOption('test')
if test == "1":
comp = GGHZZ4L
comp.files = [ 'root://eoscms//eos/cms/store/mc/Phys14DR/GluGluToHToZZTo4L_M-125_13TeV-powheg-pythia6/MINIAODSIM/PU20bx25_tsg_PHYS14_25_V1-v1/00000/3295EF7C-2070-E411-89C4-7845C4FC35DB.root ' ]
comp.splitFactor = 1
comp.fineSplitFactor = 1 if getHeppyOption('single') else 5
selectedComponents = [ comp ]
elif test == '2':
for comp in selectedComponents:
comp.files = comp.files[:1]
comp.splitFactor = 1
comp.fineSplitFactor = 1
elif test == 'Debug':
comp = GGHZZ4L
comp.files = [ 'root://eoscms//eos/cms/store/mc/Phys14DR/GluGluToHToZZTo4L_M-125_13TeV-powheg-pythia6/MINIAODSIM/PU20bx25_tsg_PHYS14_25_V1-v1/00000/3295EF7C-2070-E411-89C4-7845C4FC35DB.root ' ]
comp.splitFactor = 1