)
from madminer.utils.interfaces.delphes import run_delphes
from madminer.utils.interfaces.delphes_root import parse_delphes_root_file
from madminer.utils.interfaces.hepmc import extract_weight_order
from madminer.utils.interfaces.lhe import parse_lhe_file, extract_nuisance_parameters_from_lhe_file
logger = logging.getLogger(__name__)
h5_file = sys.argv[1]
event_path = sys.argv[2]
input_file = sys.argv[3]
benchmark_file = sys.argv[4]
mg_dir = '/madminer/software/MG5_aMC_v2_6_2'
dp = DelphesReader(h5_file)
#### get benchmark name of the job
file = open(benchmark_file, 'r')
benchmark = str(file.read())
print(benchmark)
print(type(benchmark))
####
dp.add_sample(
lhe_filename=event_path +
'/unweighted_events.lhe.gz', #'/madminer/code/mg_processes/signal/Events/run_01/unweighted_events.lhe.gz'
hepmc_filename=event_path +
'/tag_1_pythia8_events.hepmc.gz', #'/madminer/code/mg_processes/signal/Events/run_01/tag_1_pythia8_events.hepmc.gz'
is_background=False,
sampled_from_benchmark=benchmark, #'sm',
miner.run_multiple(
sample_benchmarks=additional_benchmarks,
mg_directory=mg_dir,
temp_directory='./',
mg_process_directory=
'/data_CMS/cms/cortinovis/ewdim6/mg_processes_ew_1M_az/signal_pythia2',
proc_card_file='cards/proc_card_signal.dat',
param_card_template_file='cards/param_card_template.dat',
pythia8_card_file='cards/pythia8_card.dat',
run_card_files=['cards/run_card_signal_small.dat'],
log_directory='/data_CMS/cms/cortinovis/ewdim6/logs_ew_1M_az/signal2',
initial_command="PYTHONPATH=/usr/lib64/python",
python2_override=True,
)
delphes = DelphesReader(
'/data_CMS/cms/cortinovis/ewdim6/data_ew_1M_az/setup.h5')
delphes.add_sample(
lhe_filename=
'/data_CMS/cms/cortinovis/ewdim6/mg_processes_ew_1M_az/signal_pythia/Events/run_01/unweighted_events.lhe.gz',
hepmc_filename=
'/data_CMS/cms/cortinovis/ewdim6/mg_processes_ew_1M_az/signal_pythia/Events/run_01/tag_1_pythia8_events.hepmc.gz',
sampled_from_benchmark='sm',
is_background=False,
k_factor=1,
)
for i, benchmark in enumerate(additional_benchmarks):
delphes.add_sample(
lhe_filename=
'/data_CMS/cms/cortinovis/ewdim6/mg_processes_ew_1M_az/signal_pythia2/Events/run_0{}/unweighted_events.lhe.gz'
pythia8_card_file='cards/pythia8_card.dat',
param_card_template_file='cards/param_card_template.dat',
run_card_file='cards/run_card_background.dat',
log_directory='logs/background',
)
"""
# Finally, note that both `MadMiner.run()` and `MadMiner.run_multiple()` have a `only_create_script` keyword. If that is set to True, MadMiner will not start the event generation directly, but prepare folders with all the right settings and ready-to-run bash scripts. This might make it much easier to generate Events on a high-performance computing system.
# ## 2. Run Delphes
# The `madminer.delphes.DelphesReader` class wraps around Delphes, a popular fast detector simulation, to simulate the effects of the detector.
# In[9]:
delphes = DelphesReader('data/setup.h5')
# After creating the `DelphesReader` object, one can add a number of event samples (the output of running MadGraph and Pythia in step 1 above) with the `add_sample()` function.
#
# In addition, you have to provide the information which sample was generated from which benchmark with the `sampled_from_benchmark` keyword, and set `is_background=True` for all background samples.
# In[10]:
delphes.add_sample(
lhe_filename=
'mg_processes/signal_pythia_runIter{}/Events/run_01/unweighted_events.lhe.gz'
.format(runIteration),
hepmc_filename=
'mg_processes/signal_pythia_runIter{}/Events/run_01/tag_1_pythia8_events.hepmc.gz'
.format(runIteration),
sampled_from_benchmark='sm',
mg_dir = '/home/software/MG5_aMC_v2_6_2/'
# ## 1. Generate events
# Let's load our setup:
# In[4]:
miner = MadMiner()
miner.load("data/setup.h5")
delphes = DelphesReader('data/setup.h5')
# After creating the `DelphesReader` object, one can add a number of event samples (the output of running MadGraph and Pythia in step 1 above) with the `add_sample()` function.
#
# In addition, you have to provide the information which sample was generated from which benchmark with the `sampled_from_benchmark` keyword, and set `is_background=True` for all background samples.
# careful to keep same order of benchmarks as in the simulation file
benchmarks = ['sm', 'no-higgs','1.5pow4-higgs', '2pow4-higgs']
# Works up to 8 benchmarks
assert len(benchmarks) < 9
path = "./mg_processes/"
for i, benchmark in enumerate(benchmarks):
lhe_filename_list = [f for f in glob.glob(path + \