def test_1(self):
"""Test FCC component creation"""
dset = Dataset(dataset_name_fccsw, dataset_pattern_fccsw,
cache=False)
comp = FCCComponent(dataset_name_fccsw,
xSection=dset.xsection())
self.assertListEqual(dset.list_of_good_files(),
comp.files)
import os
from fcc_datasets.fcc_component import FCCComponent
# import fcc_datasets.basedir as basedir
# basedir.basename = '/Users/cbernet/Datasets/FCC/fcc_ee_higgs/samples'
ZH = FCCComponent('heppy/ee_to_ZH_to_nunuww/CLD/Jun13/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=1)
ZZ = FCCComponent('heppy/ee_to_ZH_to_nunuww/CLD/Jun13/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000 * 200.,
uncertainty=1)
WW = FCCComponent('heppy/ee_to_ZH_to_nunuww/CLD/Jun13/ee_to_WW_Dec6_large',
xSection=1.64e-8,
nGenEvents=20000 * 200.,
uncertainty=1)
ffbar = FCCComponent(
'heppy/ee_to_ZH_to_nunuww/CLD/Jun13/ee_to_ffbar_Sep12_B_4',
xSection=7.9e-8,
nGenEvents=10000 * 491.,
uncertainty=1)
# mode = 'pythia/ee_to_ZH_Oct30'
mode = 'all'
nfiles = 1
from heppy.papas.detectors.CLIC import clic
from heppy.papas.detectors.CMS import cms
detector = clic
### definition of input samples
### from components.ZH_Znunu import components as cps
##from fcc_ee_higgs.components.all import load_components
##cps = load_components(mode='pythia')
from fcc_datasets.fcc_component import FCCComponent
zz = FCCComponent(
'pythia/ee_to_ZZ_Sep12_A_2',
splitFactor=1
)
zh = FCCComponent(
'pythia/ee_to_ZH_Oct30',
splitFactor=1
)
ffbar = FCCComponent(
'pythia/ee_to_ffbar_Sep12_B_4',
splitFactor=1
)
ww = FCCComponent(
'pythia/ee_to_WW_Dec6_large',
splitFactor=1
# import pdb; pdb.set_trace()
# mode = 'pythia/ee_to_ZH_Oct30'
mode = 'pythia/ee_to_ZZ_Sep12_A_2'
# mode = 'all'
nfiles = 1
# mode = 'test'
### definition of input samples
### from components.ZH_Znunu import components as cps
##from fcc_ee_higgs.components.all import load_components
##cps = load_components(mode='pythia')
from fcc_datasets.fcc_component import FCCComponent
zh = FCCComponent(
'pythia/ee_to_ZH_Oct30',
splitFactor=1
)
zz = FCCComponent(
'pythia/ee_to_ZZ_Sep12_A_2',
splitFactor=1
)
cpslist = [
zh,
zz
]
cps = dict( (c.name, c) for c in cpslist)
selectedComponents = cps.values()
jet_correction = True
mode = 'all'
nfiles = None
from heppy.papas.detectors.CLIC import clic
from heppy.papas.detectors.CMS import cms
detector = clic
### definition of input samples
### from components.ZH_Znunu import components as cps
##from fcc_ee_higgs.components.all import load_components
##cps = load_components(mode='pythia')
from fcc_datasets.fcc_component import FCCComponent
zz = FCCComponent('pythia/ee_to_ZZ_Sep12_A_2', splitFactor=1)
zh = FCCComponent('pythia/ee_to_ZH_Oct30', splitFactor=1)
wwh = FCCComponent('pythia/ee_WW_to_Hnunu_Sep25', splitFactor=1)
ffbar = FCCComponent('pythia/ee_to_ffbar_Sep12_B_4', splitFactor=1)
ww = FCCComponent('pythia/ee_to_WW_Dec6_large', splitFactor=1)
from fcc_ee_higgs.components.tools import get_components
selectedComponents = get_components(mode, [wwh], nfiles)
# read FCC EDM events from the input root file(s)
# do help(Reader) for more information
from heppy.analyzers.fcc.Reader import Reader
import os
from fcc_datasets.fcc_component import FCCComponent
# For ZH, ZZ, WW, there is a cut E>10 on reconstructed leptons!
ZH = FCCComponent(
'heppy/ee_to_ZH_to_qqwwll/CLD/Jun6/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000*100.,
uncertainty=0.2
)
ZZ = FCCComponent(
'heppy/ee_to_ZH_to_qqwwll/CLD/Jun6/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000*200.,
uncertainty=0.2
)
WW = FCCComponent(
'heppy/ee_to_ZH_to_qqwwll/CLD/Jun6/ee_to_WW_Dec6_large',
xSection=1.64e-8,
nGenEvents=20000*199.,
uncertainty=0.2
)
ll = FCCComponent(
# cut on lepton energy ok
'heppy/ee_to_ZH_to_qqwwll/CLD/Jun6/ee_to_2l_Mar8',
xSection=9.3e-9,
nGenEvents=20000*499,
uncertainty=0.2
import os
from fcc_datasets.fcc_component import FCCComponent
ZH = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CMS/Jan30/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 98.,
uncertainty=1)
ZZ = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CMS/Jan30/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000 * 198.,
uncertainty=1)
WW = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CMS/Jan30/ee_to_WW_Dec6_large',
xSection=1.64e-8,
nGenEvents=20000 * 197.,
uncertainty=1)
##ffbar = FCCComponent(
## 'heppy/ee_to_ZH_Z_to_nunu_H_to_bb/Oct23/pythia/ee_to_ffbar_Sep12_B_4',
## xSection=7.9e-8,
## nGenEvents=10000*500.
##)
import os
from fcc_datasets.fcc_component import FCCComponent
ZH = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CLIC_FCCee/Jan30/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=0.1)
ZZ = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CLIC_FCCee/Jan30/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000 * 200.,
uncertainty=0.05)
WW = FCCComponent(
'heppy/ee_to_ZH_Z_to_ll/CLIC_FCCee/Jan30/ee_to_WW_Dec6_large',
xSection=1.64e-8,
nGenEvents=20000 * 200.,
uncertainty=0.05)
##ffbar = FCCComponent(
## 'heppy/ee_to_ZH_Z_to_nunu_H_to_bb/Oct23/pythia/ee_to_ffbar_Sep12_B_4',
## xSection=7.9e-8,
## nGenEvents=10000*500.
##)
import os
from fcc_datasets.fcc_component import FCCComponent
# For ZH, ZZ, WW, there is a cut E>10 on reconstructed leptons!
ZH = FCCComponent('heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar7/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=0.1)
ZZ = FCCComponent('heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar7/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000 * 200.,
uncertainty=0.05)
WW = FCCComponent('heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar7/ee_to_WW_Dec6_large',
xSection=1.64e-8,
nGenEvents=20000 * 200.,
uncertainty=0.05)
ll = FCCComponent(
# cut on lepton energy ok
'heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar8/ee_to_2l_Mar8',
xSection=9.3e-9,
nGenEvents=20000 * 499,
uncertainty=0.05)
import os
from fcc_datasets.fcc_component import FCCComponent
# import fcc_datasets.basedir as basedir
# basedir.basename = '/Users/cbernet/Datasets/FCC/fcc_ee_higgs/samples'
ZH = FCCComponent('heppy/ee_to_ZH_Z_to_nunu_H_to_bb/CLIC_FCCee/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=1)
ZZ = FCCComponent(
'/heppy/ee_to_ZH_Z_to_nunu_H_to_bb/CLIC_FCCee/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000 * 200.,
uncertainty=1)
ffbar = FCCComponent(
'heppy/ee_to_ZH_Z_to_nunu_H_to_bb/CLIC_FCCee/ee_to_ffbar_Sep12_B_4',
xSection=7.9e-8,
nGenEvents=10000 * 497.,
uncertainty=1)
mode = 'all'
nfiles = sys.maxint
# nfiles = 1
# mode = 'test'
min_gen_z = 0
min_rec_z = 1
from heppy.papas.detectors.CLIC import clic
from heppy.papas.detectors.CMS import cms
detector = clic
from fcc_datasets.fcc_component import FCCComponent
zh = FCCComponent(
# 'pythia/ee_to_ZH_Oct30',
# 'pythia/ee_to_ZH_Z_to_mumu_noBS_Oct1',
'pythia/ee_to_ZH_BS_Oct2',
splitFactor=4)
zz = FCCComponent(
# 'pythia/ee_to_ZZ_Sep12_A_2',
'pythia/ee_to_ZZ_BS_Oct2',
splitFactor=1)
ww = FCCComponent('pythia/ee_to_WW_Dec6_large', splitFactor=1)
##ffbar = FCCComponent(
## 'pythia/ee_to_ffbar_Sep12_B_4',
## splitFactor=1
##)
import os
from fcc_datasets.fcc_component import FCCComponent
# mode = 'no_beam_smearing_new'
# mode = 'isofix'
mode = 'isofix_bs'
# For ZH, ZZ, WW, there is a cut E>10 on reconstructed leptons!
ZH = FCCComponent('heppy/ee_to_ZH_to_llww/CLD/May23/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=0.2)
ZH_bs = FCCComponent('heppy/ee_to_ZH_to_llww/CLD/ee_to_ZH_BS_je5_Oct2',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=0.2)
ZH_nobs_new = FCCComponent('heppy/ee_to_ZH_to_llww/CLD/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=0.2)
ZH_isofix = FCCComponent('heppy/ee_to_ZH_to_llww/CLD/IsoFix/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=0.2)
ZH_isofix_bs = FCCComponent(
'heppy/ee_to_ZH_to_llww/CLD/IsoFix/ee_to_ZH_BS_Oct2',
xSection=2.e-10,
import os
from fcc_datasets.fcc_component import FCCComponent
beam_smearing = False
# For ZH, ZZ, WW, there is a cut E>10 on reconstructed leptons!
ZH = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CLD/Mar28/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=0.2)
ZH_bs = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CLD/ee_to_ZH_BS_Oct2',
xSection=2.e-10,
nGenEvents=5000 * 97.,
uncertainty=0.2)
ZZ = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CLD/Mar28/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000 * 199.,
uncertainty=0.2)
ZZ_bs = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CLD/ee_to_ZZ_BS_Oct2',
xSection=1.35e-9,
nGenEvents=10000 * 199.,
uncertainty=0.2)
WW = FCCComponent('heppy/ee_to_ZH_Z_to_ll/CLD/Mar28/ee_to_WW_Dec6_large',
xSection=1.64e-8,
nGenEvents=20000 * 200.,
uncertainty=0.2)
import os
from fcc_datasets.fcc_component import FCCComponent
ZH = FCCComponent('heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar6/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=0.1)
ZZ = FCCComponent('heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar6/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000 * 200.,
uncertainty=0.05)
WW = FCCComponent('heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar6/ee_to_WW_Dec6_large',
xSection=1.64e-8,
nGenEvents=20000 * 200.,
uncertainty=0.05)
ffbar = FCCComponent(
'heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar6/ee_to_ffbar_Sep12_B_4',
xSection=7.9e-8,
nGenEvents=10000 * 491.,
uncertainty=0.05)
ll = FCCComponent('heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar6/ee_to_ffbar_2l_Mar6',
xSection=9.3e-9,
nGenEvents=10000 * 97,
uncertainty=0.05)
# import pdb; pdb.set_trace()
# mode = 'pythia/ee_to_ZH_Oct30'
# mode = 'pythia/ee_to_ZZ_Sep12_A_2'
mode = 'all'
nfiles = 4
# nfiles = 1
# mode = 'test'
min_gen_z = 0
min_rec_z = 1
from heppy.papas.detectors.CLIC import clic
from heppy.papas.detectors.CMS import cms
detector = clic
from fcc_datasets.fcc_component import FCCComponent
zh = FCCComponent('pythia/ee_to_ZH_Z_to_mumu_Oct30', splitFactor=4)
cpslist = [
zh,
# zz, ww, ffbar2l
]
from fcc_ee_higgs.components.tools import get_components
selectedComponents = get_components(mode, cpslist, nfiles)
##zh.files = 'ee_ZH_Zmumu_Htautau.root'
##zh.splitFactor = 1
# read FCC EDM events from the input root file(s)
# do help(Reader) for more information
from heppy.analyzers.fcc.Reader import Reader
# if collection is listed then print loop.event.papasevent will include the collections
Event.print_patterns=['zeds*', 'higgs*', 'jets*', 'bquarks', 'recoil*', 'collections']
# definition of the collider
# help(Collider) for more information
from heppy.configuration import Collider
Collider.BEAMS = 'ee'
Collider.SQRTS = 240.
# setting the base path for samples here, for testing
import fcc_datasets.basedir as basedir
from fcc_datasets.fcc_component import FCCComponent
comp = FCCComponent(
'papas/ee_to_ZZ_condor_A_703',
cache=False,
splitFactor=1
)
comp.splitFactor = len(comp.files)
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'
)
import os
from fcc_datasets.fcc_component import FCCComponent
# import fcc_datasets.basedir as basedir
# basedir.basename = '/Users/cbernet/Datasets/FCC/fcc_ee_higgs/samples'
ZHnunu = FCCComponent(
'heppy/ee_to_ZH_Z_to_nunu_H_to_bb/Oct23/pythia/ee_to_ZH_Z_to_nunu_Jun21_A_1',
xSection=1.3e-11 * 3,
nGenEvents=5000 * 100.)
ZZ = FCCComponent(
'heppy/ee_to_ZH_Z_to_nunu_H_to_bb/Oct23/pythia/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000 * 200.)
ffbar = FCCComponent(
'heppy/ee_to_ZH_Z_to_nunu_H_to_bb/Oct23/pythia/ee_to_ffbar_Sep12_B_4',
xSection=7.9e-8,
nGenEvents=10000 * 500.)
mode = 'pythia/ee_to_ZH_Oct30'
nfiles = 1
# mode = 'all'
# nfiles = None
from heppy.papas.detectors.CLIC import clic
from heppy.papas.detectors.CMS import cms
detector = clic
### definition of input samples
from fcc_datasets.fcc_component import FCCComponent
zh = FCCComponent('pythia/ee_to_ZH_Oct30', splitFactor=4)
zh_qqww = FCCComponent('pythia/ee_ZH_Zqq_HWW_Wll', splitFactor=4)
zh_qqtt = FCCComponent('pythia/ee_ZH_Zqq_Htt_tll', splitFactor=4)
zz = FCCComponent('pythia/ee_to_ZZ_Sep12_A_2', splitFactor=1)
ww = FCCComponent('pythia/ee_to_WW_Dec6_large', splitFactor=1)
ffbar2l = FCCComponent('pythia/ee_to_2l_Mar8', splitFactor=1)
from fcc_ee_higgs.components.tools import get_components
selectedComponents = get_components(mode, [zh], nfiles)
# read FCC EDM events from the input root file(s)
nfiles = 1
mode = 'all'
nfiles = None
# mode = 'pythia/ee_to_ZH_Oct30'
# nfiles = 1
from heppy.papas.detectors.CLIC import clic
from heppy.papas.detectors.CMS import cms
detector = clic
### definition of input samples
from fcc_datasets.fcc_component import FCCComponent
zh = FCCComponent('pythia/ee_to_ZH_Oct30', splitFactor=4)
# zh_qqww = FCCComponent(
# 'pythia/ee_ZH_Zqq_HWW_Wll',
# splitFactor=4
# )
# zh_qqtt = FCCComponent(
# 'pythia/ee_ZH_Zqq_Htt_tll',
# splitFactor=4
# )
zz = FCCComponent('pythia/ee_to_ZZ_Sep12_A_2', splitFactor=1)
ww = FCCComponent('pythia/ee_to_WW_Dec6_large', splitFactor=1)
import os
from fcc_datasets.fcc_component import FCCComponent
# import fcc_datasets.basedir as basedir
# basedir.basename = '/Users/cbernet/Datasets/FCC/fcc_ee_higgs/samples'
ZH = FCCComponent(
'heppy/ee_to_ZH_Z_to_nunu_H_to_bb/Dec12/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000*100.,
uncertainty=0.1
)
ZZ = FCCComponent(
'/heppy/ee_to_ZH_Z_to_nunu_H_to_bb/Dec12/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000*200.,
uncertainty=0.1
)
ffbar = FCCComponent(
'heppy/ee_to_ZH_Z_to_nunu_H_to_bb/Dec12/ee_to_ffbar_Sep12_B_4',
xSection=7.9e-8,
nGenEvents=10000*500.,
uncertainty=0.1
)
WW = FCCComponent(
'heppy/ee_to_ZH_Z_to_nunu_H_to_bb/Dec12/ee_to_WW_Dec6_large',
xSection=1.6e-8,
nGenEvents=20000*195.,
uncertainty=0.1
import os
from fcc_datasets.fcc_component import FCCComponent
# import fcc_datasets.basedir as basedir
# basedir.basename = '/Users/cbernet/Datasets/FCC/fcc_ee_higgs/samples'
##ZH = FCCComponent(
## 'heppy/ee_to_ZH_to_qqbb/May30/ee_to_ZH_Oct30',
## xSection=2.e-10,
## nGenEvents=5000*93,
## uncertainty=1
##)
ZH = FCCComponent(
'heppy/ee_to_ZH_to_qqbb/Sep25/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000*98,
uncertainty=1
)
ZZ = FCCComponent(
'heppy/ee_to_ZH_to_qqbb/May30/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000*178.,
uncertainty=1
)
WW = FCCComponent(
'heppy/ee_to_ZH_to_qqbb/May30/ee_to_WW_Dec6_large',
xSection=1.64e-8,
nGenEvents=20000*181.,
uncertainty=1
random.seed(0xdeadbeef)
# definition of the collider
from heppy.configuration import Collider
Collider.BEAMS = 'ee'
Collider.SQRTS = 240.
# mode = 'pythia/ee_to_ZH_Oct30'
mode = 'all'
nfiles = None
from fcc_datasets.fcc_component import FCCComponent
zh = FCCComponent('pythia/ee_to_ZH_Oct30', splitFactor=4)
zh = FCCComponent('pythia/ee_to_ZH_Oct30', splitFactor=1)
ffbar = FCCComponent('pythia/ee_to_ffbar_Sep12_B_4', splitFactor=1)
ww = FCCComponent('pythia/ee_to_WW_Dec6_large', splitFactor=1)
from fcc_ee_higgs.components.tools import get_components
selectedComponents = get_components(mode, [zh], nfiles)
# read FCC EDM events from the input root file(s)
# do help(Reader) for more information
from heppy.analyzers.fcc.Reader import Reader
import os
from fcc_datasets.fcc_component import FCCComponent
ZH = FCCComponent('heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar2/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000 * 100.,
uncertainty=0.1)
ZZ = FCCComponent('heppy/ee_to_ZH_Z_to_lltautau/CLD/Mar2/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000 * 199.,
uncertainty=0.05)
import os
from fcc_datasets.fcc_component import FCCComponent
# import fcc_datasets.basedir as basedir
# basedir.basename = '/Users/cbernet/Datasets/FCC/fcc_ee_higgs/samples'
ZH = FCCComponent(
'heppy/ee_to_ZH_to_nunubb/CLD/May17/ee_to_ZH_Oct30',
xSection=2.e-10,
nGenEvents=5000*99.,
uncertainty=1
)
WWH = FCCComponent(
'heppy/ee_to_ZH_to_nunubb/CLD/Sep25/ee_WW_to_Hnunu_Sep25',
xSection=6.379e-12,
nGenEvents=5000*100.,
uncertainty=1
)
ZZ = FCCComponent(
'heppy/ee_to_ZH_to_nunubb/CLD/May17/ee_to_ZZ_Sep12_A_2',
xSection=1.35e-9,
nGenEvents=10000*195.,
uncertainty=1
)
WW = FCCComponent(
'heppy/ee_to_ZH_to_nunubb/CLD/May17/ee_to_WW_Dec6_large',
xSection=1.64e-8,
