def test_list(self):
event = Event(0)
event.the_list = range(10)
cfg_ana = cfg.Analyzer(Selector,
output='filtered',
input_objects='the_list',
filter_func=lambda x: x % 2 == 0)
cfg_comp = cfg.Component('test', files=[])
filter = Selector(cfg_ana, cfg_comp, self.outdir)
filter.process(event)
self.assertItemsEqual(event.filtered, [0, 2, 4, 6, 8])
def test_dict(self):
event = Event(0)
event.the_dict = dict([(x, x**2) for x in range(10)])
cfg_ana = cfg.Analyzer(Selector,
output='filtered',
input_objects='the_dict',
filter_func=lambda x: x == 9)
cfg_comp = cfg.Component('test', files=[])
filter = Selector(cfg_ana, cfg_comp, self.outdir)
filter.process(event)
self.assertDictEqual(event.filtered, {3: 9})
import heppy.framework.config as cfg
import logging
import imp
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sample = imp.load_source(
'heppylist',
'/afs/cern.ch/work/h/helsens/public/FCCDicts/HELHC_heppySampleList_helhc_v01.py'
)
comp = cfg.Component(
'example',
files=[
"/eos/experiment/fcc/helhc/generation/DelphesEvents/helhc_v01/mgp8_pp_tth0123j_5f_hbb/events_000838653.root"
])
selectedComponents = [comp]
# FCChhAnalyses specific
from heppy.FCChhAnalyses.analyzers.Reader import Reader
source = cfg.Analyzer(
Reader,
weights='mcEventWeights',
gen_particles='skimmedGenParticles',
met='met',
electrons='electrons',
electronITags='electronITags',
electronsToMC='electronsToMC',
import os, sys
import copy
import heppy.framework.config as cfg
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sys.path.append(os.path.join(os.getenv('FCCDICTSDIR', deffccdicts), ''))
comp = cfg.Component(
'example',
files=[
"/eos/experiment/fcc/ee/generation/DelphesEvents/fcc_v01/p8_ee_ZH_ecm240/events_086259737.root"
])
from FCCee_heppySampleList_fcc_v01 import *
selectedComponents = [p8_ee_ZH_ecm240, p8_ee_ZZ_ecm240, p8_ee_WW_ecm240]
p8_ee_ZH_ecm240.splitFactor = 100
p8_ee_ZZ_ecm240.splitFactor = 100
p8_ee_WW_ecm240.splitFactor = 100
#selectedComponents = [comp]
from FCCeeAnalyses.analyzers.eeReader import Reader
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
# setting the random seed for reproducible results
import heppy.statistics.rrandom as random
random.seed(0xdeadbeef)
# definition of the collider
from heppy.configuration import Collider
Collider.BEAMS = 'ee'
Collider.SQRTS = 240.
# input definition
comp = cfg.Component('ee_ZH_Z_Hbb', files=['ee_ZH_Z_Hbb.root'])
selectedComponents = [comp]
# read FCC EDM events from the input root file(s)
# do help(Reader) for more information
from heppy.analyzers.fcc.Reader import Reader
source = cfg.Analyzer(Reader,
gen_particles='GenParticle',
gen_vertices='GenVertex')
# the papas simulation and reconstruction sequence
from heppy.test.papas_cfg import papas_sequence, detector
# Use a Selector to select leptons from the output of papas simulation.
# Currently, we're treating electrons and muons transparently.
# we could use two different instances for the Selector module
import os
import copy
import heppy.framework.config as cfg
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
comp = cfg.Component('example', files=['top_mu.root'])
selectedComponents = [comp]
from heppy_fcc.analyzers.FCCReader import FCCReader
source = cfg.Analyzer(FCCReader)
from ROOT import gSystem
gSystem.Load("libdatamodel")
from eventstore import EventStore as Events
from heppy_fcc.analyzers.JetClusterizer import JetClusterizer
gen_jets = cfg.Analyzer(JetClusterizer,
instance_label='gen',
particles='gen_particles_stable')
from heppy_fcc.analyzers.PFSim import PFSim
from heppy_fcc.fastsim.detectors.CMS import CMS
import os
import copy
import heppy.framework.config as cfg
import logging
import imp
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sample=imp.load_source('heppylist', '/afs/cern.ch/work/h/helsens/public/FCCDicts/FCC_heppySampleList_fcc_v02.py')
comp = cfg.Component(
'example',
files = ["/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/p8_pp_jj_lo_tagger/events_001556815.root"]
#files = ["/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/p8_pp_Zprime_20TeV_ttbar_tagger/events_000260643.root"]
#files = ["/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/p8_pp_RSGraviton_20TeV_ww_tagger/events_002015707.root"]
)
selectedComponents = [
sample.p8_pp_Zprime_20TeV_ttbar_tagger,
sample.p8_pp_RSGraviton_20TeV_ww_tagger,
sample.p8_pp_jj_lo_tagger,
]
splitFac = 20
sample.p8_pp_RSGraviton_20TeV_ww_tagger.splitFactor = splitFac
sample.p8_pp_Zprime_20TeV_ttbar_tagger.splitFactor = splitFac
sample.p8_pp_jj_lo_tagger.splitFactor = splitFac
import heppy.framework.config as cfg
from getFiles import getFiles
gun_211_0_20 = cfg.Component(
'gun_211_0_20',
files=getFiles("/Gun_221_0_20_ptflat/743_v1/AODSIM", cache=True),
)
gun_211_0_10 = cfg.Component(
'gun_211_0_10',
files=getFiles("/Gun_211_0_10_ptflat/743_v1/AODSIM"),
)
gun_211_MatEff_0_20 = cfg.Component(
'gun_211_MatEff_0_20',
files=getFiles("/Gun_221_MatEff_0_20_ptflat/743_v1/AODSIM", cache=True),
)
gun_22_0_50 = cfg.Component(
'gun_22_0_50',
files=getFiles("/Gun_22_0_50/743_v2/RECOSIM"),
)
gun_22_0_50_eta3 = cfg.Component(
'gun_22_0_50_eta3',
files=getFiles("/Gun_22_0_50_eta3/743_v2/AODSIM"),
)
gun_130_0_50 = cfg.Component(
'gun_130_0_50',
import heppy.framework.config as cfg
import logging
import imp
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sample = imp.load_source(
'heppylist',
'/afs/cern.ch/work/h/helsens/public/FCCDicts/HELHC_heppySampleList_helhc_v01.py'
)
comp = cfg.Component(
'example',
files=[
"/eos/experiment/fcc/helhc/generation/DelphesEvents/helhc_v01/p8_pp_Zprime_10TeV_ttbar/events_198227905.root"
])
selectedComponents = [
sample.p8_pp_Zprime_2TeV_ttbar,
sample.p8_pp_Zprime_4TeV_ttbar,
sample.p8_pp_Zprime_5TeV_ttbar,
sample.p8_pp_Zprime_6TeV_ttbar,
sample.p8_pp_Zprime_8TeV_ttbar,
sample.p8_pp_Zprime_10TeV_ttbar,
sample.p8_pp_Zprime_12TeV_ttbar,
sample.p8_pp_Zprime_14TeV_ttbar,
sample.mgp8_pp_jj_5f_HT_500_1000,
sample.mgp8_pp_jj_5f_HT_1000_2000,
sample.mgp8_pp_jj_5f_HT_2000_5000,
import imp
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sample = imp.load_source(
'heppylist',
'/afs/cern.ch/work/h/helsens/public/FCCDicts/HELHC_heppySampleList_helhc_v01.py'
)
comp = cfg.Component(
'example',
#files = ["/eos/experiment/fcc/helhc/generation/DelphesEvents/helhc_v01/p8_pp_Zprime_10TeV_ttbar/events_198227905.root"]
files=[
"/eos/experiment/fcc/helhc/generation/DelphesEvents/helhc_v01/p8_pp_ZprimeSSM_4TeV_jj//events_087545417.root"
]
# files = ["/eos/experiment/fcc/helhc/generation/DelphesEvents/helhc_v01/mgp8_pp_jj_5f_HT_1000_2000/events_096768951.root"]
# files = ["/eos/experiment/fcc/helhc/generation/DelphesEvents/helhc_v01/mgp8_pp_jj_5f_HT_2000_5000/events_163309325.root"]
)
selectedComponents = [
sample.p8_pp_Zprime_2TeV_ttbar,
sample.p8_pp_Zprime_4TeV_ttbar,
sample.p8_pp_Zprime_5TeV_ttbar,
sample.p8_pp_Zprime_6TeV_ttbar,
sample.p8_pp_Zprime_8TeV_ttbar,
sample.p8_pp_Zprime_10TeV_ttbar,
sample.p8_pp_Zprime_12TeV_ttbar,
sample.p8_pp_Zprime_14TeV_ttbar,
sample.mgp8_pp_jj_5f_HT_500_1000,
import os
import heppy.framework.config as cfg
import logging
logging.basicConfig(level=logging.WARNING)
# input component
# several input components can be declared and added to the list of selected components
input_component = cfg.Component('Z2uu', files = ['/afs/cern.ch/work/a/ansemkiv/private/FCC/analysis/Z2uubar_100.root'])
selected_components = [input_component]
# analyzers
# analyzer for signal events
from heppy_fcc.analyzers.Z2UUAnalyzer import Z2UUAnalyzer
ana = cfg.Analyzer(Z2UUAnalyzer)
# definition of a sequence of analyzers, the analyzers will process each event in this order
sequence = cfg.Sequence([ana])
# finalization of the configuration object.
from ROOT import gSystem
gSystem.Load('libdatamodel')
from EventStore import EventStore as Events
config = cfg.Config(components = selected_components, sequence = sequence, services = [], events_class = Events)
import os
import copy
import heppy.framework.config as cfg
import logging
#import sys
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
comp = cfg.Component('example', files=['ee_ZH_Zmumu_Hbb.root'])
selectedComponents = [comp]
from heppy.analyzers.fcc.Reader import Reader
source = cfg.Analyzer(
Reader,
gen_particles='GenParticle',
)
# Use a Selector to select stable gen particles for simulation
# from the output of "source"
# help(Selector) for more information
from heppy.analyzers.Selector import Selector
gen_particles_stable = cfg.Analyzer(
Selector,
output='gen_particles_stable',
# output = 'particles',
input_objects='gen_particles',
filter_func=lambda x: x.status() == 1 and x.pdgid() not in [12, 14, 16] and
reload(logging)
logging.basicConfig(level=logging.WARNING)
# setting the random seed for reproducible results
import heppy.statistics.rrandom as random
random.seed(0xdeadbeef)
# definition of the collider
from heppy.configuration import Collider
Collider.BEAMS = 'ee'
Collider.SQRTS = 240.
# input definition
import glob
ZH = cfg.Component(
'ZH',
files=glob.glob(
'/afs/cern.ch/work/n/nsolomon/FCC/WorkDir/Mumu/ee_ZH_Z_Hmumu_0.root'))
ZH.splitFactor = len(ZH.files)
Bkg = cfg.Component(
'Bkg',
files=glob.glob(
'/afs/cern.ch/work/n/nsolomon/FCC/WorkDir/backgroundmumu_*.root'))
Bkg.splitFactor = len(Bkg.files)
selectedComponents = [ZH, Bkg]
# read FCC EDM events from the input root file(s)
# do help(Reader) for more information
from heppy.analyzers.fcc.Reader import Reader
source = cfg.Analyzer(
Reader,
gen_particles='GenParticle',
import os
import copy
import heppy.framework.config as cfg
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
comp = cfg.Component(
'example',
files = 'example.root'
)
selectedComponents = [comp]
from heppy_fcc.analyzers.FCCReader import FCCReader
source = cfg.Analyzer(
FCCReader,
mode = 'ee',
gen_particles = 'GenParticle'
)
from ROOT import gSystem
gSystem.Load("libdatamodelDict")
from EventStore import EventStore as Events
from heppy_fcc.analyzers.PFSim import PFSim
from heppy_fcc.fastsim.detectors.CMS import CMS
papas = cfg.Analyzer(
import os, sys
import copy
import heppy.framework.config as cfg
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sys.path.append('/afs/cern.ch/work/h/helsens/public/FCCDicts/')
comp = cfg.Component(
'example',
#files = ["/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v01/mgp8_pp_hh01j_5f/events0.root"]
files=[
"/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/mgp8_pp_hh_lambda100_5f_haa/events_000000127.root"
])
from FCC_heppySampleList_fcc_v02 import *
selectedComponents = [
mgp8_pp_jjaa_5f, mgp8_pp_hh_lambda090_5f_haa, mgp8_pp_hh_lambda095_5f_haa,
mgp8_pp_hh_lambda096_5f_haa, mgp8_pp_hh_lambda097_5f_haa,
mgp8_pp_hh_lambda098_5f_haa, mgp8_pp_hh_lambda099_5f_haa,
mgp8_pp_hh_lambda100_5f_haa, mgp8_pp_hh_lambda101_5f_haa,
mgp8_pp_hh_lambda102_5f_haa, mgp8_pp_hh_lambda103_5f_haa,
mgp8_pp_hh_lambda104_5f_haa, mgp8_pp_hh_lambda105_5f_haa,
mgp8_pp_hh_lambda110_5f_haa, mgp8_pp_tth01j_5f_haa, mgp8_pp_bbh_haa,
mgp8_pp_vh012j_5f_haa, mgp8_pp_h012j_5f_haa
]
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
# setting the random seed for reproducible results
import heppy.statistics.rrandom as random
random.seed(0xdeadbeef)
# definition of the collider
from heppy.configuration import Collider
Collider.BEAMS = 'ee'
Collider.SQRTS = 91.
# input definition
ee_Z_bbbar = cfg.Component('ee_Z_bbbar', files=['data/ee_Z_ddbar.root'])
selectedComponents = [ee_Z_bbbar]
# read FCC EDM events from the input root file(s)
# do help(Reader) for more information
from heppy.analyzers.fcc.Reader import Reader
source = cfg.Analyzer(Reader,
gen_particles='GenParticle',
gen_vertices='GenVertex')
from heppy.test.papas_cfg import gen_particles_stable, papas_sequence, detector, papas, papasdisplay, papasdisplaycompare
from heppy.test.papas_cfg import papasdisplaycompare as display
# Make jets
from heppy.analyzers.fcc.JetClusterizer import JetClusterizer
import os, sys
import copy
import heppy.framework.config as cfg
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sys.path.append('/afs/cern.ch/work/h/helsens/public/FCCDicts/')
comp = cfg.Component(
'example',
#files = ["root://eospublic.cern.ch///eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v01/pp_ll012j_5f_HT_15000_25000/events100.root"]
files = ["events100.root"]
)
from heppySampleList_fcc_v01 import *
#from heppySampleList_cms import *
selectedComponents = [
pp_ttv01j_5f,
pp_v0123j_5f,
]
pp_ttv01j_5f.splitFactor = 10
pp_v0123j_5f.splitFactor = 10
import os, sys
import copy
import heppy.framework.config as cfg
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sys.path.append('/afs/cern.ch/work/h/helsens/public/FCCDicts/')
comp = cfg.Component(
'example',
files=[
"root://eospublic.cern.ch//eos/fcc/hh/generation/DelphesEvents/decay/pp_h012j_5f_za/events997.root"
])
from heppySampleList import *
'''selectedComponents = [
pp_h012j_5f_HT_0_100_hlla,
pp_h012j_5f_HT_100_400_hlla,
pp_h012j_5f_HT_400_1000_hlla,
pp_h012j_5f_HT_1000_1900_hlla,
pp_h012j_5f_HT_1900_4400_hlla,
pp_h012j_5f_HT_4400_8500_hlla,
pp_h012j_5f_HT_8500_100000_hlla,
pp_h012j_5f_hlla,
pp_tth01j_5f_HT_0_1100_hlla,
pp_tth01j_5f_HT_1100_2700_hlla,
pp_tth01j_5f_HT_2700_4900_hlla,
import os, sys
import copy, math
import heppy.framework.config as cfg
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sys.path.append('/afs/cern.ch/work/h/helsens/public/FCCDicts/')
comp = cfg.Component(
'example',
files = ["/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/mgp8_pp_h012j_5f_haa/events_000001373.root"]
)
from FCC_heppySampleList_fcc_v02 import *
selectedComponents = [
mgp8_pp_h012j_5f_haa,
mgp8_pp_tth01j_5f_haa,
mgp8_pp_vbf_h01j_5f_haa,
mgp8_pp_vh012j_5f_haa,
#mgp8_gg_aa01j_mhcut_5f,
#mgp8_pp_aa012j_mhcut_5f
mgp8_pp_aa012j_mhcut_5f_HT_0_100,
mgp8_pp_aa012j_mhcut_5f_HT_100_300,
mgp8_pp_aa012j_mhcut_5f_HT_1100_100000,
mgp8_pp_aa012j_mhcut_5f_HT_300_500,
mgp8_pp_aa012j_mhcut_5f_HT_500_700,
import os
import heppy.framework.config as cfg
import logging
from ROOT import gSystem
from EventStore import EventStore as Events
from heppy_fcc.analyzers.BackgroundBs2DsDsKWithDs2PiPiPiPiAndDs2TauNuAnalyzer import BackgroundBs2DsDsKWithDs2PiPiPiPiAndDs2TauNuAnalyzer
logging.basicConfig(level=logging.WARNING)
# input component
# several input components can be declared and added to the list of selected components
input_component = cfg.Component(
'ILD-like',
files=[
'/afs/cern.ch/work/a/ansemkiv/private/FCC/analysis/background_Bs2DsDsK_with_Ds2PiPiPiPi_and_Ds2TauNu_100k.root'
])
selected_components = [input_component]
# analyzers
# analyzer for Bs -> Ds Ds K* events
bgana = cfg.Analyzer(
BackgroundBs2DsDsKWithDs2PiPiPiPiAndDs2TauNuAnalyzer,
smear_momentum=True,
momentum_x_resolution=0.01,
momentum_y_resolution=0.01,
momentum_z_resolution=0.01,
smear_pv=True,
import os
import copy
import heppy.framework.config as cfg
import sys
import logging
import imp
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sys.path.append('/afs/cern.ch/work/h/helsens/public/FCCDicts/')
comp = cfg.Component(
'example',
files = ["/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/p8_pp_ExcitedQ_30TeV_qq/events_113722907.root"]
)
sample=imp.load_source('heppylist', '/afs/cern.ch/work/h/helsens/public/FCCDicts/FCC_heppySampleList_fcc_v02.py')
selectedComponents = [
sample.p8_pp_ExcitedQ_15TeV_qq,
sample.p8_pp_ExcitedQ_20TeV_qq,
sample.p8_pp_ExcitedQ_25TeV_qq,
sample.p8_pp_ExcitedQ_30TeV_qq,
sample.p8_pp_ExcitedQ_35TeV_qq,
sample.p8_pp_ExcitedQ_40TeV_qq,
sample.p8_pp_ExcitedQ_45TeV_qq,
sample.p8_pp_ExcitedQ_50TeV_qq,
sample.mgp8_pp_jj_lo,
]
import os
import copy
import heppy.framework.config as cfg
import sys
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sys.path.append('/afs/cern.ch/work/h/helsens/public/FCCDicts/')
comp = cfg.Component(
'example',
files=[
"root://eospublic.cern.ch///eos/fcc/hh/generation/DelphesEvents/fcc_v01/pp_Zprime_5TeV_ttbar/events0.root"
])
from heppySampleList_fcc_v01 import *
selectedComponents = [
pp_jj012j_5f_HT_0_500,
pp_jj012j_5f_HT_500_1000,
pp_jj012j_5f_HT_1000_2000,
pp_jj012j_5f_HT_2000_4000,
pp_jj012j_5f_HT_4000_7200,
pp_jj012j_5f_HT_7200_15000,
pp_jj012j_5f_HT_15000_25000,
pp_jj012j_5f_HT_25000_35000,
pp_jj012j_5f_HT_35000_100000,
pp_Zprime_2TeV_ttbar,
pp_Zprime_5TeV_ttbar,
pp_Zprime_10TeV_ttbar,
import os
import copy
import heppy.framework.config as cfg
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
comp = cfg.Component('example', files=['FCCDelphesOutput.root'])
selectedComponents = [comp]
from heppy.analyzers.fcc.Reader import Reader
source = cfg.Analyzer(
Reader,
weights='mcEventWeights',
gen_particles='skimmedGenParticles',
electrons='electrons',
electronITags='electronITags',
electronsToMC='electronsToMC',
muons='muons',
muonITags='muonITags',
muonsToMC='muonsToMC',
)
from ROOT import gSystem
gSystem.Load("libdatamodelDict")
from EventStore import EventStore as Events
#############################
import heppy.framework.config as cfg
import logging
import imp
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sample = imp.load_source(
'heppylist',
os.path.join(os.getenv('FCCDICTSDIR', deffccdicts), '') +
'FCC_heppySampleList_fcc_v02.py')
comp = cfg.Component(
'example',
files=[
"/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/p8_pp_RSGraviton_10TeV_ww/events_193939464.root"
])
selectedComponents = [
sample.p8_pp_RSGraviton_2TeV_ww,
sample.p8_pp_RSGraviton_5TeV_ww,
sample.p8_pp_RSGraviton_10TeV_ww,
sample.p8_pp_RSGraviton_15TeV_ww,
sample.p8_pp_RSGraviton_20TeV_ww,
sample.p8_pp_RSGraviton_25TeV_ww,
sample.p8_pp_RSGraviton_30TeV_ww,
sample.p8_pp_RSGraviton_35TeV_ww,
sample.mgp8_pp_jj_5f_HT_500_1000,
sample.mgp8_pp_jj_5f_HT_1000_2000,
sample.mgp8_pp_jj_5f_HT_2000_5000,
import heppy.framework.config as cfg
import logging
import imp
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sample = imp.load_source(
'heppylist',
'/afs/cern.ch/work/h/helsens/public/FCCDicts/FCC_heppySampleList_fcc_v02.py'
)
comp = cfg.Component(
'example',
files=[
"/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/p8_pp_ZprimeSSM_20TeV_ll/events_135275618.root"
])
selectedComponents = [
sample.p8_pp_ZprimeSSM_2TeV_ll,
sample.p8_pp_ZprimeSSM_4TeV_ll,
sample.p8_pp_ZprimeSSM_5TeV_ll,
sample.p8_pp_ZprimeSSM_6TeV_ll,
sample.p8_pp_ZprimeSSM_8TeV_ll,
sample.p8_pp_ZprimeSSM_10TeV_ll,
sample.p8_pp_ZprimeSSM_12TeV_ll,
sample.p8_pp_ZprimeSSM_14TeV_ll,
sample.p8_pp_ZprimeSSM_15TeV_ll,
sample.p8_pp_ZprimeSSM_16TeV_ll,
sample.p8_pp_ZprimeSSM_18TeV_ll,
## Configuration script for signal (B0s -> tau+ tau-) analyzer
import os
import heppy.framework.config as cfg
import logging
logging.basicConfig(level=logging.WARNING)
# input component
# several input components can be declared and added to the list of selected components
input_component = cfg.Component('Bs2tautau', files = ['/afs/cern.ch/work/a/ansemkiv/private/FCC/analysis/Bs2tautau_10k.root'])
selected_components = [input_component]
# analyzers
# analyzer for signal events
from heppy_fcc.analyzers.Bs2TauTauAnalyzer import Bs2TauTauAnalyzer
sigana = cfg.Analyzer(Bs2TauTauAnalyzer,
smear_pv = True,
# IDL-like res
pv_x_resolution = 0.0025,
pv_y_resolution = 0.0025,
pv_z_resolution = 0.0025,
# progressive res
# pv_x_resolution = 0.001,
# pv_y_resolution = 0.001,
# pv_z_resolution = 0.001,
# smear_sv = True,
# IDL-like res
# sv_x_resolution = 0.007,
# sv_y_resolution = 0.007,
# comment the following line to see all the collections stored in the event
# if collection is listed then print loop.event.papasevent will include the collections
# Event.print_patterns=['rec_particles', 'gen_particles_stable']
Event.print_patterns = ['rec_particles', 'gen_particles_stable', 'sum_all*']
# definition of the collider
# help(Collider) for more information
from heppy.configuration import Collider
Collider.BEAMS = 'ee'
Collider.SQRTS = 240.
# definition of an input sample (also called a component)
comp = cfg.Component(
'ee_Z_ee',
files=[
# here we have a single input root file.
# the absolute path must be used to be able to run on the batch.
os.environ['HEPPY'] + '/test/ee_Z_ee.root'
])
# selecting the list of components to be processed. Here only one.
selectedComponents = [comp]
# read FCC EDM events from the input root file(s)
# do help(Reader) for more information
from heppy.analyzers.fcc.Reader import Reader
source = cfg.Analyzer(Reader,
gen_particles='GenParticle',
gen_vertices='GenVertex')
# importing the papas simulation and reconstruction sequence,
from heppy.framework.event import Event
# comment the following line to see all the collections stored in the event
# if collection is listed then print loop.event.papasevent will include the collections
# Event.print_patterns=['*']
# definition of the collider
# help(Collider) for more information
from heppy.configuration import Collider
Collider.BEAMS = 'ee'
Collider.SQRTS = 240.
# definition of an input sample (also called a component)
comp = cfg.Component(
'ee_ZH_Zmumu_Hbb',
files=[
# here we have a single input root file.
# the absolute path must be used to be able to run on the batch.
os.path.abspath('ee_ZH_Zmumu_Hbb.root')
])
# selecting the list of components to be processed. Here only one.
selectedComponents = [comp]
# read FCC EDM events from the input root file(s)
# do help(Reader) for more information
from heppy.analyzers.fcc.Reader import Reader
source = cfg.Analyzer(Reader,
gen_particles='GenParticle',
gen_vertices='GenVertex')
# Use a Selector to select leptons from the output of papas simulation.
import os, sys
import copy
import heppy.framework.config as cfg
import logging
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sys.path.append('/eos/user/s/selvaggi/FCCSW/dicts')
comp = cfg.Component(
'example',
files=[
"/eos/user/s/selvaggi/FCCSW/DelphesEvents/cms/pp_jjaa_5f_14TeV/events25.root"
])
from heppySampleList import *
selectedComponents = [pp_jjaa_5f_14TeV, pp_hh_14TeV]
pp_jjaa_5f_14TeV.splitFactor = 20
pp_hh_14TeV.splitFactor = 20
pp_jjja_5f_14TeV.splitFactor = 20
selectedComponents = [comp]
#from heppy.FCChhAnalyses.analyzers.OldReader import Reader
from heppy.FCChhAnalyses.analyzers.Reader import Reader
import logging
import imp
# next 2 lines necessary to deal with reimports from ipython
logging.shutdown()
reload(logging)
logging.basicConfig(level=logging.WARNING)
sample = imp.load_source(
'heppylist',
'/afs/cern.ch/work/h/helsens/public/FCCDicts/FCC_heppySampleList_fcc_v02.py'
)
comp = cfg.Component(
'example',
#files = ["/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/p8_pp_jj_lo_tagger/events_001556815.root"]
#files = ["/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/p8_pp_Zprime_20TeV_ttbar_tagger/events_000260643.root"]
files=[
"/eos/experiment/fcc/hh/generation/DelphesEvents/fcc_v02/mgp8_pp_jj_lo_5f/events_052659188.root"
])
selectedComponents = [
#sample.p8_pp_Zprime_5TeV_ttbar,
sample.p8_pp_Zprime_10TeV_ttbar,
sample.p8_pp_Zprime_15TeV_ttbar,
sample.p8_pp_Zprime_20TeV_ttbar,
sample.p8_pp_Zprime_25TeV_ttbar,
sample.p8_pp_Zprime_30TeV_ttbar,
sample.p8_pp_Zprime_35TeV_ttbar,
#sample.p8_pp_Zprime_40TeV_ttbar,
sample.mgp8_pp_jj_lo_5f,
sample.mgp8_pp_tt_lo,
