#use tfile service to provide a single TFile to all modules where they
#can write any root object. If the name is 'outputfile' or the one specified in treeProducer
#also the treeProducer uses this file
from PhysicsTools.HeppyCore.framework.services.tfile import TFileService
output_service = cfg.Service(
TFileService,
'outputfile',
name="outputfile",
fname='tree.root',
option='recreate'
)
# the following two lines are just for automatic testing
# they are not needed for running on your own samples
from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles
testfiles=miniAodFiles()
sample = cfg.Component(
#specify the file you want to run on
# files = ["/scratch/arizzi/Hbb/CMSSW_7_2_2_patch2/src/VHbbAnalysis/Heppy/test/ZLL-8A345C56-6665-E411-9C25-1CC1DE04DF20.root"],
files = testfiles,
name="SingleSample", isMC=False,isEmbed=False
)
from PhysicsTools.Heppy.utils.cmsswPreprocessor import CmsswPreprocessor
preprocessor = CmsswPreprocessor("makeAK5Jets.py")
# the following is declared in case this cfg is used in input to the heppy.py script
from PhysicsTools.HeppyCore.framework.eventsfwlite import Events
selectedComponents = [sample]
config = cfg.Config( components = selectedComponents,
]
#use tfile service to provide a single TFile to all modules where they
#can write any root object. If the name is 'outputfile' or the one specified in treeProducer
#also the treeProducer uses this file
from PhysicsTools.HeppyCore.framework.services.tfile import TFileService
output_service = cfg.Service(TFileService,
'outputfile',
name="outputfile",
fname='tree.root',
option='recreate')
# the following two lines are just for automatic testing
# they are not needed for running on your own samples
from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles
testfiles = miniAodFiles()
sample = cfg.Component(
#specify the file you want to run on
# files = ["/scratch/arizzi/Hbb/CMSSW_7_2_2_patch2/src/VHbbAnalysis/Heppy/test/ZLL-8A345C56-6665-E411-9C25-1CC1DE04DF20.root"],
files=testfiles,
name="SingleSample",
isMC=False,
isEmbed=False)
from PhysicsTools.Heppy.utils.cmsswPreprocessor import CmsswPreprocessor
preprocessor = CmsswPreprocessor("makeAK5Jets.py")
# the following is declared in case this cfg is used in input to the heppy.py script
from PhysicsTools.HeppyCore.framework.eventsfwlite import Events
selectedComponents = [sample]
from PhysicsTools.Heppy.analyzers.objects.TauAnalyzer import TauAnalyzer
TauAna = TauAnalyzer.defaultConfig
from PhysicsTools.Heppy.analyzers.objects.JetAnalyzer import JetAnalyzer
JetAna = JetAnalyzer.defaultConfig
sequence = [VertexAna, LepAna, TauAna, PhoAna, JetAna, treeProducer]
from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles
sample = cfg.Component(
# files = "/scratch/arizzi/heppy/CMSSW_7_2_0_pre8/src/PhysicsTools/Heppy/test/E21AD523-E548-E411-8DF6-00261894388F.root",
files=miniAodFiles(),
name="ATEST",
isMC=False,
isEmbed=False)
# the following is declared in case this cfg is used in input to the heppy.py script
selectedComponents = [sample]
from PhysicsTools.HeppyCore.framework.eventsfwlite import Events
config = cfg.Config(components=selectedComponents,
sequence=sequence,
events_class=Events)
# and the following runs the process directly
if __name__ == '__main__':
from PhysicsTools.HeppyCore.framework.looper import Looper
from __future__ import print_function
import os
import PhysicsTools.HeppyCore.framework.config as cfg
from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles
# set to True if you want several parallel processes
multi_thread = False
# input component
# several input components can be declared,
# and added to the list of selected components
inputSample = cfg.MCComponent(
'test_component',
files=miniAodFiles(), # returns a list of file for this release
# a list of local or xrootd files can be specified by hand.
)
if 'RelValZMM' not in inputSample.files[0]:
print('''WARNING: this tutorial is supposed to run on Z->mumu events.
Do not expect meaningful results for this sample:''')
print(inputSample)
if multi_thread:
inputSample.splitFactor = len(inputSample.files)
selectedComponents = [inputSample]
# a very simple muon analyzer
# read miniaod muons and wrap them in python muons
from PhysicsTools.Heppy.analyzers.examples.SimpleMuonAnalyzer import SimpleMuonAnalyzer
muons = cfg.Analyzer(
import os
import PhysicsTools.HeppyCore.framework.config as cfg
from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles
# input component
# several input components can be declared,
# and added to the list of selected components
inputSample = cfg.Component(
'test_component',
files = miniAodFiles(),
)
inputSample.isMC = True
# inputSample.splitFactor = 2
selectedComponents = [inputSample]
from PhysicsTools.Heppy.analyzers.examples.SimpleJetAnalyzer import SimpleJetAnalyzer
jets = cfg.Analyzer(
SimpleJetAnalyzer,
'jets',
ptmin = 30.
)
from PhysicsTools.Heppy.analyzers.examples.SimpleTreeAnalyzer import SimpleTreeAnalyzer
tree = cfg.Analyzer(
SimpleTreeAnalyzer
)
# definition of a sequence of analyzers,
import os
import PhysicsTools.HeppyCore.framework.config as cfg
from PhysicsTools.Heppy.utils.miniAodFiles import miniAodFiles
# set to True if you want several parallel processes
multi_thread = False
# input component
# several input components can be declared,
# and added to the list of selected components
inputSample = cfg.MCComponent(
"test_component",
files=miniAodFiles(), # returns a list of file for this release
# a list of local or xrootd files can be specified by hand.
)
if "RelValZMM" not in inputSample.files[0]:
print """WARNING: this tutorial is supposed to run on Z->mumu events.
Do not expect meaningful results for this sample:"""
print inputSample
if multi_thread:
inputSample.splitFactor = len(inputSample.files)
selectedComponents = [inputSample]
# a very simple muon analyzer
# read miniaod muons and wrap them in python muons
from PhysicsTools.Heppy.analyzers.examples.SimpleMuonAnalyzer import SimpleMuonAnalyzer
muons = cfg.Analyzer(SimpleMuonAnalyzer, "muons")
