{
"cells": [
{
"cell_type": "markdown",
"id": "30b2a263-4237-4909-8acd-2fcd0972be62",
"metadata": {},
"source": [
"# Event Generation Process\n",
"\n",
"The event generation process is a critical first step in our simulation pipeline. It involves creating simulated particle interactions that mimic real-world physics events. This process is typically carried out using specialized scripts, such as `gen_zee.py`, which are designed to generate events based on specific parameters and configurations.\n",
"\n",
"During event generation, we define various parameters, including the number of events to generate, the output file format, and the characteristics of the particles involved. The script utilizes a random number generator to simulate the inherent uncertainties in particle interactions, ensuring that the generated events are representative of actual experimental conditions.\n",
"\n",
"The generated events can then be used for further analysis, including shower propagation, digitalization, and event reconstruction. This allows researchers to study the properties of particles and their interactions in a controlled environment, facilitating a deeper understanding of high-energy physics phenomena.\n",
"\n",
"## Table of Contents:\n",
"\n",
"* [How to create $Z\\rightarrow ee$ events?](#zee)\n",
"* [How to create minumum bias events?](#mb)\n",
"* [How to create Jets?](#jets)\n",
"* [How to create single electrons?](#electron)\n",
"* [How to use the job configuration?](#jobs)\n",
"* [Let's run!](#run)\n"
]
},
{
"cell_type": "markdown",
"id": "9ef83ec5-21a8-484d-8456-8aca042d610a",
"metadata": {},
"source": [
"\n",
"## How to create $Z\\rightarrow ee$ events?\n",
"\n",
"Let's create a Z boson decay into an electron-positron pair using the `gen_zee.py` command. 🎉 This process is crucial for simulating the interactions that occur in high-energy physics experiments. By generating these events, we can study the properties of the Z boson and its decay channels in detail. \n",
"\n",
"To initiate the event generation, we specify parameters such as the average pileup, event number, and output file name. This allows us to customize the simulation according to our experimental needs."
]
},
{
"cell_type": "code",
"execution_count": 2,
"id": "496b8481-9920-4e06-8383-6bbdcd5a870d",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"\u001b[1;32mUsage:\u001b[0m \u001b[1;38;5;244mgen_zee.py\u001b[0m [\u001b[32m--run-number\u001b[0m \u001b[34mRUN_NUMBER\u001b[0m] [\u001b[32m--output-level\u001b[0m \u001b[34mOUTPUT_LEVEL\u001b[0m]\n",
" [\u001b[32m--eta-max\u001b[0m \u001b[34mETA_MAX\u001b[0m] [\u001b[32m--force-forward-electron\u001b[0m]\n",
" [\u001b[32m--zero-vertex-particles\u001b[0m] [\u001b[32m--zee-file\u001b[0m \u001b[34mZEE_FILE\u001b[0m]\n",
" [\u001b[32m-e\u001b[0m \u001b[34mEVENT_NUMBERS\u001b[0m] \u001b[32m-o\u001b[0m \u001b[34mOUTPUT_FILE\u001b[0m [\u001b[32m-s\u001b[0m \u001b[34mSEED\u001b[0m]\n",
" [\u001b[32m--nov\u001b[0m \u001b[34mNUMBER_OF_EVENTS\u001b[0m] [\u001b[32m--events-per-job\u001b[0m \u001b[34mEVENTS_PER_JOB\u001b[0m]\n",
" [\u001b[32m-nt\u001b[0m \u001b[34mNUMBER_OF_THREADS\u001b[0m] [\u001b[32m-m\u001b[0m] [\u001b[32m--overwrite\u001b[0m] [\u001b[32m--dry-run\u001b[0m]\n",
" [\u001b[32m--job-file\u001b[0m \u001b[34mJOB_FILE\u001b[0m]\n",
"\n",
"\u001b[1;32mOptions:\u001b[0m\n",
" \u001b[32m--run-number\u001b[0m \u001b[34mRUN_NUMBER\u001b[0m\n",
" \u001b[38;5;244mThe run number.\u001b[0m\n",
" \u001b[32m--output-level\u001b[0m \u001b[34mOUTPUT_LEVEL\u001b[0m\n",
" \u001b[38;5;244mThe output level messenger.\u001b[0m\n",
" \u001b[32m--eta-max\u001b[0m \u001b[34mETA_MAX\u001b[0m \u001b[38;5;244mThe eta max used in generator.\u001b[0m\n",
" \u001b[32m--force-forward-electron\u001b[0m\n",
" \u001b[38;5;244mForce at least one electron into forward region.\u001b[0m\n",
" \u001b[32m--zero-vertex-particles\u001b[0m\n",
" \u001b[38;5;244mFix the z vertex position in simulation to zero for\u001b[0m\n",
" \u001b[38;5;244mall selected particles. It is applied only at G4 step,\u001b[0m\n",
" \u001b[38;5;244mnot in generation.\u001b[0m\n",
" \u001b[32m--zee-file\u001b[0m \u001b[34mZEE_FILE\u001b[0m \u001b[38;5;244mThe pythia zee file configuration.\u001b[0m\n",
" \u001b[32m-e\u001b[0m, \u001b[32m--event-numbers\u001b[0m \u001b[34mEVENT_NUMBERS\u001b[0m\n",
" \u001b[38;5;244mThe event number list separated by ','. e.g.\u001b[0m\n",
" \u001b[32m--event-numbers\u001b[0m\u001b[38;5;244m '0,1,2,3'\u001b[0m\n",
" \u001b[32m-o\u001b[0m, \u001b[32m--output-file\u001b[0m \u001b[34mOUTPUT_FILE\u001b[0m\n",
" \u001b[38;5;244mThe output file.\u001b[0m\n",
" \u001b[32m-s\u001b[0m, \u001b[32m--seed\u001b[0m \u001b[34mSEED\u001b[0m \u001b[38;5;244mThe pythia seed (zero is the clock system)\u001b[0m\n",
" \u001b[32m--nov\u001b[0m, \u001b[32m--number-of-events\u001b[0m \u001b[34mNUMBER_OF_EVENTS\u001b[0m\n",
" \u001b[38;5;244mThe total number of events to run.\u001b[0m\n",
" \u001b[32m--events-per-job\u001b[0m \u001b[34mEVENTS_PER_JOB\u001b[0m\n",
" \u001b[38;5;244mThe number of events per job\u001b[0m\n",
" \u001b[32m-nt\u001b[0m, \u001b[32m--number-of-threads\u001b[0m \u001b[34mNUMBER_OF_THREADS\u001b[0m\n",
" \u001b[38;5;244mThe number of threads\u001b[0m\n",
" \u001b[32m-m\u001b[0m, \u001b[32m--merge\u001b[0m \u001b[38;5;244mMerge all files.\u001b[0m\n",
" \u001b[32m--overwrite\u001b[0m \u001b[38;5;244mRerun all jobs.\u001b[0m\n",
" \u001b[32m--dry-run\u001b[0m \u001b[38;5;244mPerform a dry run without executing jobs.\u001b[0m\n",
" \u001b[32m--job-file\u001b[0m \u001b[34mJOB_FILE\u001b[0m \u001b[38;5;244mThe JSON file used to configure this job automatically\u001b[0m\n",
" \u001b[38;5;244mby overwritting some arguments.\u001b[0m\n"
]
}
],
"source": [
"!gen_zee.py"
]
},
{
"cell_type": "markdown",
"id": "36bd50ae-9e41-4a79-b51f-614812e660ce",
"metadata": {},
"source": [
"This command will create 20 events, 5 event per job using 2 threads. The `dry_run` command is used here just to avoid the execution. This is a test."
]
},
{
"cell_type": "code",
"execution_count": 3,
"id": "cb307891-4b40-47aa-8f92-fd82df4a5c39",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"{'Zee.EVT.0.root': {'evt': [0, 1, 2, 3, 4], 'seed': 413654000},\n",
" 'Zee.EVT.1.root': {'evt': [5, 6, 7, 8, 9], 'seed': 813847340},\n",
" 'Zee.EVT.2.root': {'evt': [10, 11, 12, 13, 14], 'seed': 451585302},\n",
" 'Zee.EVT.3.root': {'evt': [15, 16, 17, 18, 19], 'seed': 43469774}}\n",
"{'eta_max': 3.2,\n",
" 'force_forward_electron': False,\n",
" 'logging_level': 'INFO',\n",
" 'run_number': 0,\n",
" 'zee_file': '/home/joao.pinto/git_repos/lorenzetti/build/data/evtgen/zee_config.cmnd',\n",
" 'zero_vertex_particles': False}\n"
]
}
],
"source": [
"!gen_zee.py -o Zee.EVT.root --nov 20 --events-per-job 5 -nt 2 --dry-run"
]
},
{
"cell_type": "markdown",
"id": "ff07d9ae-0027-4175-bdd0-bb5e578f91ee",
"metadata": {},
"source": [
"**NOTE**: To run the generation, just remove the `dry-run` argument."
]
},
{
"cell_type": "markdown",
"id": "e2b44be5-1b3c-43c8-b845-0c1600076fb0",
"metadata": {},
"source": [
"\n",
"## How to create minimum bias events?"
]
},
{
"cell_type": "code",
"execution_count": 4,
"id": "3b0c9fc2-0297-4407-ae06-ecb3c12658a3",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"\u001b[1;32mUsage:\u001b[0m \u001b[1;38;5;244mgen_minbias.py\u001b[0m [\u001b[32m--run-number\u001b[0m \u001b[34mRUN_NUMBER\u001b[0m] [\u001b[32m--output-level\u001b[0m \u001b[34mOUTPUT_LEVEL\u001b[0m]\n",
" [\u001b[32m--eta-max\u001b[0m \u001b[34mETA_MAX\u001b[0m]\n",
" [\u001b[32m--pileup-per-bunch-crossing\u001b[0m \u001b[34mPILEUP_PER_BUNCH_CROSSING\u001b[0m]\n",
" [\u001b[32m--bc-id-start\u001b[0m \u001b[34mBC_ID_START\u001b[0m] [\u001b[32m--bc-id-end\u001b[0m \u001b[34mBC_ID_END\u001b[0m]\n",
" [\u001b[32m--bc-duration\u001b[0m \u001b[34mBC_DURATION\u001b[0m] [\u001b[32m--pileup-file\u001b[0m \u001b[34mPILEUP_FILE\u001b[0m]\n",
" [\u001b[32m-e\u001b[0m \u001b[34mEVENT_NUMBERS\u001b[0m] \u001b[32m-o\u001b[0m \u001b[34mOUTPUT_FILE\u001b[0m [\u001b[32m-s\u001b[0m \u001b[34mSEED\u001b[0m]\n",
" [\u001b[32m--nov\u001b[0m \u001b[34mNUMBER_OF_EVENTS\u001b[0m]\n",
" [\u001b[32m--events-per-job\u001b[0m \u001b[34mEVENTS_PER_JOB\u001b[0m]\n",
" [\u001b[32m-nt\u001b[0m \u001b[34mNUMBER_OF_THREADS\u001b[0m] [\u001b[32m-m\u001b[0m] [\u001b[32m--overwrite\u001b[0m] [\u001b[32m--dry-run\u001b[0m]\n",
" [\u001b[32m--job-file\u001b[0m \u001b[34mJOB_FILE\u001b[0m]\n",
"\n",
"\u001b[1;32mOptions:\u001b[0m\n",
" \u001b[32m--run-number\u001b[0m \u001b[34mRUN_NUMBER\u001b[0m\n",
" \u001b[38;5;244mThe run number.\u001b[0m\n",
" \u001b[32m--output-level\u001b[0m \u001b[34mOUTPUT_LEVEL\u001b[0m\n",
" \u001b[38;5;244mThe output level messenger.\u001b[0m\n",
" \u001b[32m--eta-max\u001b[0m \u001b[34mETA_MAX\u001b[0m \u001b[38;5;244mMaximum Eta.\u001b[0m\n",
" \u001b[32m--pileup-per-bunch-crossing\u001b[0m \u001b[34mPILEUP_PER_BUNCH_CROSSING\u001b[0m\n",
" \u001b[38;5;244mThe fixed number of pileup for each bunch crossing.\u001b[0m\n",
" \u001b[32m--bc-id-start\u001b[0m \u001b[34mBC_ID_START\u001b[0m\n",
" \u001b[38;5;244mThe bunch crossing id start.\u001b[0m\n",
" \u001b[32m--bc-id-end\u001b[0m \u001b[34mBC_ID_END\u001b[0m\n",
" \u001b[38;5;244mThe bunch crossing id end.\u001b[0m\n",
" \u001b[32m--bc-duration\u001b[0m \u001b[34mBC_DURATION\u001b[0m\n",
" \u001b[38;5;244mThe bunch crossing duration (in nanoseconds).\u001b[0m\n",
" \u001b[32m--pileup-file\u001b[0m \u001b[34mPILEUP_FILE\u001b[0m\n",
" \u001b[38;5;244mThe pythia pileup file configuration.\u001b[0m\n",
" \u001b[32m-e\u001b[0m, \u001b[32m--event-numbers\u001b[0m \u001b[34mEVENT_NUMBERS\u001b[0m\n",
" \u001b[38;5;244mThe event number list separated by ','. e.g.\u001b[0m\n",
" \u001b[32m--event-numbers\u001b[0m\u001b[38;5;244m '0,1,2,3'\u001b[0m\n",
" \u001b[32m-o\u001b[0m, \u001b[32m--output-file\u001b[0m \u001b[34mOUTPUT_FILE\u001b[0m\n",
" \u001b[38;5;244mThe output file.\u001b[0m\n",
" \u001b[32m-s\u001b[0m, \u001b[32m--seed\u001b[0m \u001b[34mSEED\u001b[0m \u001b[38;5;244mThe pythia seed (zero is the clock system)\u001b[0m\n",
" \u001b[32m--nov\u001b[0m, \u001b[32m--number-of-events\u001b[0m \u001b[34mNUMBER_OF_EVENTS\u001b[0m\n",
" \u001b[38;5;244mThe total number of events to run.\u001b[0m\n",
" \u001b[32m--events-per-job\u001b[0m \u001b[34mEVENTS_PER_JOB\u001b[0m\n",
" \u001b[38;5;244mThe number of events per job\u001b[0m\n",
" \u001b[32m-nt\u001b[0m, \u001b[32m--number-of-threads\u001b[0m \u001b[34mNUMBER_OF_THREADS\u001b[0m\n",
" \u001b[38;5;244mThe number of threads\u001b[0m\n",
" \u001b[32m-m\u001b[0m, \u001b[32m--merge\u001b[0m \u001b[38;5;244mMerge all files.\u001b[0m\n",
" \u001b[32m--overwrite\u001b[0m \u001b[38;5;244mRerun all jobs.\u001b[0m\n",
" \u001b[32m--dry-run\u001b[0m \u001b[38;5;244mPerform a dry run without executing jobs.\u001b[0m\n",
" \u001b[32m--job-file\u001b[0m \u001b[34mJOB_FILE\u001b[0m \u001b[38;5;244mThe JSON file used to configure this job automatically\u001b[0m\n",
" \u001b[38;5;244mby overwritting some arguments.\u001b[0m\n"
]
}
],
"source": [
"!gen_minbias.py"
]
},
{
"cell_type": "markdown",
"id": "b566e5cf-db9f-4e38-868a-3421132896a2",
"metadata": {},
"source": [
"This command creates minimum bias files using the initial seed provided by the user. Each bunch crossing will be filled with pileup equal to the number specified by the `--pileup-per-bunch-crossing` parameter, ensuring that `--eta-max` ($\\eta_{max}\\leq2.5$). In total, 100 events will be generated, with 10 events per job, utilizing only one thread."
]
},
{
"cell_type": "code",
"execution_count": 5,
"id": "e3513bde-24a5-415a-8fb9-07b0d7c23919",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"{'MB.EVT.0.root': {'evt': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9], 'seed': 20387086},\n",
" 'MB.EVT.1.root': {'evt': [10, 11, 12, 13, 14, 15, 16, 17, 18, 19],\n",
" 'seed': 874807135},\n",
" 'MB.EVT.2.root': {'evt': [20, 21, 22, 23, 24, 25, 26, 27, 28, 29],\n",
" 'seed': 471648295},\n",
" 'MB.EVT.3.root': {'evt': [30, 31, 32, 33, 34, 35, 36, 37, 38, 39],\n",
" 'seed': 152360332},\n",
" 'MB.EVT.4.root': {'evt': [40, 41, 42, 43, 44, 45, 46, 47, 48, 49],\n",
" 'seed': 354838602},\n",
" 'MB.EVT.5.root': {'evt': [50, 51, 52, 53, 54, 55, 56, 57, 58, 59],\n",
" 'seed': 442384353},\n",
" 'MB.EVT.6.root': {'evt': [60, 61, 62, 63, 64, 65, 66, 67, 68, 69],\n",
" 'seed': 388525902},\n",
" 'MB.EVT.7.root': {'evt': [70, 71, 72, 73, 74, 75, 76, 77, 78, 79],\n",
" 'seed': 669594479},\n",
" 'MB.EVT.8.root': {'evt': [80, 81, 82, 83, 84, 85, 86, 87, 88, 89],\n",
" 'seed': 435449645},\n",
" 'MB.EVT.9.root': {'evt': [90, 91, 92, 93, 94, 95, 96, 97, 98, 99],\n",
" 'seed': 840205357}}\n",
"{'bc_id_end': 4,\n",
" 'bc_id_start': -21,\n",
" 'eta_max': 2.5,\n",
" 'logging_level': 'INFO',\n",
" 'mb_file': '/home/joao.pinto/git_repos/lorenzetti/build/data/evtgen/minbias_config.cmnd',\n",
" 'pileup_per_bunch_crossing': 1.0,\n",
" 'run_number': 0}\n"
]
}
],
"source": [
"!gen_minbias.py -o MB.EVT.root --seed 512 -nt 1 --pileup-per-bunch-crossing 1 --eta-max 2.5 --events-per-job 10 --nov 100 --dry-run"
]
},
{
"cell_type": "markdown",
"id": "b5c5d703-ac51-45db-805a-11098b1b7e29",
"metadata": {},
"source": [
"**NOTE**: To run the generation, just remove the `dry-run` argument."
]
},
{
"cell_type": "markdown",
"id": "15dfe31a-8a85-4f3f-a864-31a686fcde8c",
"metadata": {},
"source": [
"\n",
"## How to create Jets?"
]
},
{
"cell_type": "code",
"execution_count": 6,
"id": "7bdd76c1-903a-49c8-9cc3-64e1b79a8046",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"\u001b[1;32mUsage:\u001b[0m \u001b[1;38;5;244mgen_jets.py\u001b[0m [\u001b[32m--run-number\u001b[0m \u001b[34mRUN_NUMBER\u001b[0m] [\u001b[32m--output-level\u001b[0m \u001b[34mOUTPUT_LEVEL\u001b[0m]\n",
" [\u001b[32m--eta-min\u001b[0m \u001b[34mETA_MIN\u001b[0m] [\u001b[32m--eta-max\u001b[0m \u001b[34mETA_MAX\u001b[0m]\n",
" [\u001b[32m--energy-min\u001b[0m \u001b[34mENERGY_MIN\u001b[0m] [\u001b[32m--energy-max\u001b[0m \u001b[34mENERGY_MAX\u001b[0m]\n",
" [\u001b[32m--jf17-file\u001b[0m \u001b[34mJF17_FILE\u001b[0m] [\u001b[32m-e\u001b[0m \u001b[34mEVENT_NUMBERS\u001b[0m] \u001b[32m-o\u001b[0m \u001b[34mOUTPUT_FILE\u001b[0m\n",
" [\u001b[32m-s\u001b[0m \u001b[34mSEED\u001b[0m] [\u001b[32m--nov\u001b[0m \u001b[34mNUMBER_OF_EVENTS\u001b[0m]\n",
" [\u001b[32m--events-per-job\u001b[0m \u001b[34mEVENTS_PER_JOB\u001b[0m] [\u001b[32m-nt\u001b[0m \u001b[34mNUMBER_OF_THREADS\u001b[0m]\n",
" [\u001b[32m-m\u001b[0m] [\u001b[32m--overwrite\u001b[0m] [\u001b[32m--dry-run\u001b[0m] [\u001b[32m--job-file\u001b[0m \u001b[34mJOB_FILE\u001b[0m]\n",
"\n",
"\u001b[1;32mOptions:\u001b[0m\n",
" \u001b[32m--run-number\u001b[0m \u001b[34mRUN_NUMBER\u001b[0m\n",
" \u001b[38;5;244mThe run number.\u001b[0m\n",
" \u001b[32m--output-level\u001b[0m \u001b[34mOUTPUT_LEVEL\u001b[0m\n",
" \u001b[38;5;244mThe output level messenger.\u001b[0m\n",
" \u001b[32m--eta-min\u001b[0m \u001b[34mETA_MIN\u001b[0m \u001b[38;5;244mThe eta min used in generator.\u001b[0m\n",
" \u001b[32m--eta-max\u001b[0m \u001b[34mETA_MAX\u001b[0m \u001b[38;5;244mThe eta max used in generator.\u001b[0m\n",
" \u001b[32m--energy-min\u001b[0m \u001b[34mENERGY_MIN\u001b[0m\n",
" \u001b[38;5;244mMinimum energy\u001b[0m\n",
" \u001b[32m--energy-max\u001b[0m \u001b[34mENERGY_MAX\u001b[0m\n",
" \u001b[38;5;244mMaximum energy\u001b[0m\n",
" \u001b[32m--jf17-file\u001b[0m \u001b[34mJF17_FILE\u001b[0m\n",
" \u001b[38;5;244mThe pythia JF17 file configuration.\u001b[0m\n",
" \u001b[32m-e\u001b[0m, \u001b[32m--event-numbers\u001b[0m \u001b[34mEVENT_NUMBERS\u001b[0m\n",
" \u001b[38;5;244mThe event number list separated by ','. e.g.\u001b[0m\n",
" \u001b[32m--event-numbers\u001b[0m\u001b[38;5;244m '0,1,2,3'\u001b[0m\n",
" \u001b[32m-o\u001b[0m, \u001b[32m--output-file\u001b[0m \u001b[34mOUTPUT_FILE\u001b[0m\n",
" \u001b[38;5;244mThe output file.\u001b[0m\n",
" \u001b[32m-s\u001b[0m, \u001b[32m--seed\u001b[0m \u001b[34mSEED\u001b[0m \u001b[38;5;244mThe pythia seed (zero is the clock system)\u001b[0m\n",
" \u001b[32m--nov\u001b[0m, \u001b[32m--number-of-events\u001b[0m \u001b[34mNUMBER_OF_EVENTS\u001b[0m\n",
" \u001b[38;5;244mThe total number of events to run.\u001b[0m\n",
" \u001b[32m--events-per-job\u001b[0m \u001b[34mEVENTS_PER_JOB\u001b[0m\n",
" \u001b[38;5;244mThe number of events per job\u001b[0m\n",
" \u001b[32m-nt\u001b[0m, \u001b[32m--number-of-threads\u001b[0m \u001b[34mNUMBER_OF_THREADS\u001b[0m\n",
" \u001b[38;5;244mThe number of threads\u001b[0m\n",
" \u001b[32m-m\u001b[0m, \u001b[32m--merge\u001b[0m \u001b[38;5;244mMerge all files.\u001b[0m\n",
" \u001b[32m--overwrite\u001b[0m \u001b[38;5;244mRerun all jobs.\u001b[0m\n",
" \u001b[32m--dry-run\u001b[0m \u001b[38;5;244mPerform a dry run without executing jobs.\u001b[0m\n",
" \u001b[32m--job-file\u001b[0m \u001b[34mJOB_FILE\u001b[0m \u001b[38;5;244mThe JSON file used to configure this job automatically\u001b[0m\n",
" \u001b[38;5;244mby overwritting some arguments.\u001b[0m\n"
]
}
],
"source": [
"!gen_jets.py"
]
},
{
"cell_type": "markdown",
"id": "a4c08152-4c83-492c-9fc2-7eccd7575eb9",
"metadata": {},
"source": [
"This command will generate jet files, producing a total of 100 events, with 10 events per job, utilizing 10 threads. The initial seed is set to 512, and the run number corresponds to the current date provided by the user."
]
},
{
"cell_type": "code",
"execution_count": 7,
"id": "24c8edcf-0850-46c3-a169-953fa2658bf8",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"{'Jets.EVT.0.root': {'evt': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9], 'seed': 20387086},\n",
" 'Jets.EVT.1.root': {'evt': [10, 11, 12, 13, 14, 15, 16, 17, 18, 19],\n",
" 'seed': 874807135},\n",
" 'Jets.EVT.2.root': {'evt': [20, 21, 22, 23, 24, 25, 26, 27, 28, 29],\n",
" 'seed': 471648295},\n",
" 'Jets.EVT.3.root': {'evt': [30, 31, 32, 33, 34, 35, 36, 37, 38, 39],\n",
" 'seed': 152360332},\n",
" 'Jets.EVT.4.root': {'evt': [40, 41, 42, 43, 44, 45, 46, 47, 48, 49],\n",
" 'seed': 354838602},\n",
" 'Jets.EVT.5.root': {'evt': [50, 51, 52, 53, 54, 55, 56, 57, 58, 59],\n",
" 'seed': 442384353},\n",
" 'Jets.EVT.6.root': {'evt': [60, 61, 62, 63, 64, 65, 66, 67, 68, 69],\n",
" 'seed': 388525902},\n",
" 'Jets.EVT.7.root': {'evt': [70, 71, 72, 73, 74, 75, 76, 77, 78, 79],\n",
" 'seed': 669594479},\n",
" 'Jets.EVT.8.root': {'evt': [80, 81, 82, 83, 84, 85, 86, 87, 88, 89],\n",
" 'seed': 435449645},\n",
" 'Jets.EVT.9.root': {'evt': [90, 91, 92, 93, 94, 95, 96, 97, 98, 99],\n",
" 'seed': 840205357}}\n",
"{'energy_max': 6500,\n",
" 'energy_min': 17,\n",
" 'eta_max': 3.2,\n",
" 'eta_min': 0.0,\n",
" 'jf17_file': '/home/joao.pinto/git_repos/lorenzetti/build/data/evtgen/jet_config.cmnd',\n",
" 'logging_level': 'INFO',\n",
" 'run_number': 20251115}\n"
]
}
],
"source": [
"!gen_jets.py -o Jets.EVT.root --nov 100 --events-per-job 10 -nt 10 -m -s 512 --run-number 20251115 --dry-run"
]
},
{
"cell_type": "markdown",
"id": "a2ca6655-f191-491c-97a8-8b58dc842d3c",
"metadata": {},
"source": [
"\n",
"## How to create single electrons?"
]
},
{
"cell_type": "code",
"execution_count": 8,
"id": "5f96ec3e-e06d-4418-9e03-749a1656268a",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"\u001b[1;32mUsage:\u001b[0m \u001b[1;38;5;244mgen_single.py\u001b[0m [\u001b[32m--run-number\u001b[0m \u001b[34mRUN_NUMBER\u001b[0m] [\u001b[32m--output-level\u001b[0m \u001b[34mOUTPUT_LEVEL\u001b[0m]\n",
" \u001b[32m-p\u001b[0m \u001b[34mPARTICLE\u001b[0m [\u001b[32m--energy-min\u001b[0m \u001b[34mENERGY_MIN\u001b[0m]\n",
" [\u001b[32m--energy-max\u001b[0m \u001b[34mENERGY_MAX\u001b[0m] [\u001b[32m--energy\u001b[0m \u001b[34mENERGY\u001b[0m] [\u001b[32m--eta\u001b[0m \u001b[34mETA\u001b[0m]\n",
" [\u001b[32m--phi\u001b[0m \u001b[34mPHI\u001b[0m] [\u001b[32m--do-eta-ranged\u001b[0m \u001b[34mDO_ETA_RANGED\u001b[0m]\n",
" [\u001b[32m--eta-min\u001b[0m \u001b[34mETA_MIN\u001b[0m] [\u001b[32m--eta-max\u001b[0m \u001b[34mETA_MAX\u001b[0m]\n",
" [\u001b[32m--do-phi-ranged\u001b[0m \u001b[34mDO_PHI_RANGED\u001b[0m] [\u001b[32m--phi-min\u001b[0m \u001b[34mPHI_MIN\u001b[0m]\n",
" [\u001b[32m--phi-max\u001b[0m \u001b[34mPHI_MAX\u001b[0m] [\u001b[32m-e\u001b[0m \u001b[34mEVENT_NUMBERS\u001b[0m] \u001b[32m-o\u001b[0m \u001b[34mOUTPUT_FILE\u001b[0m\n",
" [\u001b[32m-s\u001b[0m \u001b[34mSEED\u001b[0m] [\u001b[32m--nov\u001b[0m \u001b[34mNUMBER_OF_EVENTS\u001b[0m]\n",
" [\u001b[32m--events-per-job\u001b[0m \u001b[34mEVENTS_PER_JOB\u001b[0m] [\u001b[32m-nt\u001b[0m \u001b[34mNUMBER_OF_THREADS\u001b[0m]\n",
" [\u001b[32m-m\u001b[0m] [\u001b[32m--overwrite\u001b[0m] [\u001b[32m--dry-run\u001b[0m] [\u001b[32m--job-file\u001b[0m \u001b[34mJOB_FILE\u001b[0m]\n",
"\n",
"\u001b[1;32mOptions:\u001b[0m\n",
" \u001b[32m--run-number\u001b[0m \u001b[34mRUN_NUMBER\u001b[0m\n",
" \u001b[38;5;244mThe run number.\u001b[0m\n",
" \u001b[32m--output-level\u001b[0m \u001b[34mOUTPUT_LEVEL\u001b[0m\n",
" \u001b[38;5;244mThe output level messenger.\u001b[0m\n",
" \u001b[32m-p\u001b[0m, \u001b[32m--particle\u001b[0m \u001b[34mPARTICLE\u001b[0m\n",
" \u001b[38;5;244mThe single particle type: ['Electron', 'Neutrino',\u001b[0m\n",
" \u001b[38;5;244m'Photon', 'PiZero', 'PiPlus']\u001b[0m\n",
" \u001b[32m--energy-min\u001b[0m \u001b[34mENERGY_MIN\u001b[0m\n",
" \u001b[38;5;244mEnergy min in GeV.\u001b[0m\n",
" \u001b[32m--energy-max\u001b[0m \u001b[34mENERGY_MAX\u001b[0m\n",
" \u001b[38;5;244mEnergy max in GeV.\u001b[0m\n",
" \u001b[32m--energy\u001b[0m \u001b[34mENERGY\u001b[0m \u001b[38;5;244mEnergy in GeV.\u001b[0m\n",
" \u001b[32m--eta\u001b[0m \u001b[34mETA\u001b[0m \u001b[38;5;244mEta position.\u001b[0m\n",
" \u001b[32m--phi\u001b[0m \u001b[34mPHI\u001b[0m \u001b[38;5;244mPhi position.\u001b[0m\n",
" \u001b[32m--do-eta-ranged\u001b[0m \u001b[34mDO_ETA_RANGED\u001b[0m\n",
" \u001b[38;5;244mEnable eta range.\u001b[0m\n",
" \u001b[32m--eta-min\u001b[0m \u001b[34mETA_MIN\u001b[0m \u001b[38;5;244mMinimum Eta.\u001b[0m\n",
" \u001b[32m--eta-max\u001b[0m \u001b[34mETA_MAX\u001b[0m \u001b[38;5;244mMaximum Eta.\u001b[0m\n",
" \u001b[32m--do-phi-ranged\u001b[0m \u001b[34mDO_PHI_RANGED\u001b[0m\n",
" \u001b[38;5;244mEnable phi range.\u001b[0m\n",
" \u001b[32m--phi-min\u001b[0m \u001b[34mPHI_MIN\u001b[0m \u001b[38;5;244mMinimum Phi.\u001b[0m\n",
" \u001b[32m--phi-max\u001b[0m \u001b[34mPHI_MAX\u001b[0m \u001b[38;5;244mMaximum Phi.\u001b[0m\n",
" \u001b[32m-e\u001b[0m, \u001b[32m--event-numbers\u001b[0m \u001b[34mEVENT_NUMBERS\u001b[0m\n",
" \u001b[38;5;244mThe event number list separated by ','. e.g.\u001b[0m\n",
" \u001b[32m--event-numbers\u001b[0m\u001b[38;5;244m '0,1,2,3'\u001b[0m\n",
" \u001b[32m-o\u001b[0m, \u001b[32m--output-file\u001b[0m \u001b[34mOUTPUT_FILE\u001b[0m\n",
" \u001b[38;5;244mThe output file.\u001b[0m\n",
" \u001b[32m-s\u001b[0m, \u001b[32m--seed\u001b[0m \u001b[34mSEED\u001b[0m \u001b[38;5;244mThe pythia seed (zero is the clock system)\u001b[0m\n",
" \u001b[32m--nov\u001b[0m, \u001b[32m--number-of-events\u001b[0m \u001b[34mNUMBER_OF_EVENTS\u001b[0m\n",
" \u001b[38;5;244mThe total number of events to run.\u001b[0m\n",
" \u001b[32m--events-per-job\u001b[0m \u001b[34mEVENTS_PER_JOB\u001b[0m\n",
" \u001b[38;5;244mThe number of events per job\u001b[0m\n",
" \u001b[32m-nt\u001b[0m, \u001b[32m--number-of-threads\u001b[0m \u001b[34mNUMBER_OF_THREADS\u001b[0m\n",
" \u001b[38;5;244mThe number of threads\u001b[0m\n",
" \u001b[32m-m\u001b[0m, \u001b[32m--merge\u001b[0m \u001b[38;5;244mMerge all files.\u001b[0m\n",
" \u001b[32m--overwrite\u001b[0m \u001b[38;5;244mRerun all jobs.\u001b[0m\n",
" \u001b[32m--dry-run\u001b[0m \u001b[38;5;244mPerform a dry run without executing jobs.\u001b[0m\n",
" \u001b[32m--job-file\u001b[0m \u001b[34mJOB_FILE\u001b[0m \u001b[38;5;244mThe JSON file used to configure this job automatically\u001b[0m\n",
" \u001b[38;5;244mby overwritting some arguments.\u001b[0m\n"
]
}
],
"source": [
"!gen_single.py"
]
},
{
"cell_type": "markdown",
"id": "86bbe783-a5a9-4819-9de6-8464a3c9cd12",
"metadata": {},
"source": [
"Let's create single electron particles between an specific $\\eta$ region in the detector."
]
},
{
"cell_type": "code",
"execution_count": 9,
"id": "c22295ba-c8f5-47c6-a92d-289525a6e4d1",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"{'Electron.EVT.0.root': {'evt': [0], 'seed': 413654000},\n",
" 'Electron.EVT.1.root': {'evt': [1], 'seed': 813847340},\n",
" 'Electron.EVT.2.root': {'evt': [2], 'seed': 451585302},\n",
" 'Electron.EVT.3.root': {'evt': [3], 'seed': 43469774},\n",
" 'Electron.EVT.4.root': {'evt': [4], 'seed': 278009743},\n",
" 'Electron.EVT.5.root': {'evt': [5], 'seed': 548977049},\n",
" 'Electron.EVT.6.root': {'evt': [6], 'seed': 521760890},\n",
" 'Electron.EVT.7.root': {'evt': [7], 'seed': 434794719},\n",
" 'Electron.EVT.8.root': {'evt': [8], 'seed': 841597327},\n",
" 'Electron.EVT.9.root': {'evt': [9], 'seed': 891047769}}\n",
"{'do_eta_ranged': False,\n",
" 'do_phi_ranged': False,\n",
" 'energy': -1,\n",
" 'energy_max': -1,\n",
" 'energy_min': -1,\n",
" 'eta': 0.0,\n",
" 'eta_max': 1.8,\n",
" 'eta_min': 0.0,\n",
" 'logging_level': 'INFO',\n",
" 'particle': 'Electron',\n",
" 'phi': 1.52,\n",
" 'phi_max': 3.14,\n",
" 'phi_min': -3.14,\n",
" 'run_number': 0}\n"
]
}
],
"source": [
"!gen_single.py -p Electron --eta-min 0 --eta-max 1.8 -nt 1 --events-per-job 1 --nov 10 -o Electron.EVT.root -m --dry-run"
]
},
{
"cell_type": "markdown",
"id": "25ac41ac-1fe8-462d-9c5c-1bdf9a77735a",
"metadata": {},
"source": [
"\n",
"## How to use the job configuration?\n",
"\n",
"The job configuration can be specified using the `--job-file` parameter, allowing the user to provide a JSON file that overrides the default parameters for the job execution."
]
},
{
"cell_type": "code",
"execution_count": 10,
"id": "fb5e9247-661e-4de7-b185-6369b51f3513",
"metadata": {},
"outputs": [],
"source": [
"import json\n",
"\n",
"with open('config.json', 'w') as f:\n",
" d = {\n",
" 'run_number': 20251116,\n",
" 'eta_min' : 0,\n",
" 'eta_max' : 2.5,\n",
" 'seed' : 512,\n",
" 'particle' : 'Electron',\n",
" 'energy_min': 50, # GeV\n",
" 'energy_max': 150, # GeV,\n",
" }\n",
" json.dump(d,f)"
]
},
{
"cell_type": "code",
"execution_count": 11,
"id": "ff3c80b1-8151-44b1-924d-e562e7735772",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"overwritting run_number with value 0 to new value 20251116\n",
"overwritting particle with value Electron to new value Electron\n",
"overwritting energy_min with value -1 to new value 50\n",
"overwritting energy_max with value -1 to new value 150\n",
"overwritting eta_min with value -2.5 to new value 0\n",
"overwritting eta_max with value 2.5 to new value 2.5\n",
"overwritting seed with value 0 to new value 512\n",
"{'Electron.EVT.0.root': {'evt': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9],\n",
" 'seed': 20387086},\n",
" 'Electron.EVT.1.root': {'evt': [10, 11, 12, 13, 14, 15, 16, 17, 18, 19],\n",
" 'seed': 874807135},\n",
" 'Electron.EVT.2.root': {'evt': [20, 21, 22, 23, 24, 25, 26, 27, 28, 29],\n",
" 'seed': 471648295},\n",
" 'Electron.EVT.3.root': {'evt': [30, 31, 32, 33, 34, 35, 36, 37, 38, 39],\n",
" 'seed': 152360332},\n",
" 'Electron.EVT.4.root': {'evt': [40, 41, 42, 43, 44, 45, 46, 47, 48, 49],\n",
" 'seed': 354838602},\n",
" 'Electron.EVT.5.root': {'evt': [50, 51, 52, 53, 54, 55, 56, 57, 58, 59],\n",
" 'seed': 442384353},\n",
" 'Electron.EVT.6.root': {'evt': [60, 61, 62, 63, 64, 65, 66, 67, 68, 69],\n",
" 'seed': 388525902},\n",
" 'Electron.EVT.7.root': {'evt': [70, 71, 72, 73, 74, 75, 76, 77, 78, 79],\n",
" 'seed': 669594479},\n",
" 'Electron.EVT.8.root': {'evt': [80, 81, 82, 83, 84, 85, 86, 87, 88, 89],\n",
" 'seed': 435449645},\n",
" 'Electron.EVT.9.root': {'evt': [90, 91, 92, 93, 94, 95, 96, 97, 98, 99],\n",
" 'seed': 840205357}}\n",
"{'do_eta_ranged': False,\n",
" 'do_phi_ranged': False,\n",
" 'energy': -1,\n",
" 'energy_max': 150,\n",
" 'energy_min': 50,\n",
" 'eta': 0.0,\n",
" 'eta_max': 2.5,\n",
" 'eta_min': 0,\n",
" 'logging_level': 'INFO',\n",
" 'particle': 'Electron',\n",
" 'phi': 1.52,\n",
" 'phi_max': 3.14,\n",
" 'phi_min': -3.14,\n",
" 'run_number': 20251116}\n"
]
}
],
"source": [
"!gen_single.py --job-file config.json --nov 100 --events-per-job 10 -nt 10 --dry-run -o Electron.EVT.root -p Electron"
]
},
{
"cell_type": "markdown",
"id": "21b99973-1e6b-4232-b22a-01e574df13ce",
"metadata": {},
"source": [
"\n",
"## Let's run!"
]
},
{
"cell_type": "code",
"execution_count": 12,
"id": "6d7e7aa9-dd9e-4471-a4b1-3eca9175dc72",
"metadata": {},
"outputs": [
{
"name": "stdout",
"output_type": "stream",
"text": [
"Module libc not found.\n",
"overwritting run_number with value 0 to new value 20251116\n",
"overwritting particle with value Electron to new value Electron\n",
"overwritting energy_min with value -1 to new value 50\n",
"overwritting energy_max with value -1 to new value 150\n",
"overwritting eta_min with value -2.5 to new value 0\n",
"overwritting eta_max with value 2.5 to new value 2.5\n",
"overwritting seed with value 0 to new value 512\n",
"{'Electron.EVT.0.root': {'evt': [0, 1, 2, 3, 4, 5, 6, 7, 8, 9],\n",
" 'seed': 20387086}}\n",
"{'do_eta_ranged': False,\n",
" 'do_phi_ranged': False,\n",
" 'energy': -1,\n",
" 'energy_max': 150,\n",
" 'energy_min': 50,\n",
" 'eta': 0.0,\n",
" 'eta_max': 2.5,\n",
" 'eta_min': 0,\n",
" 'logging_level': 'INFO',\n",
" 'particle': 'Electron',\n",
" 'phi': 1.52,\n",
" 'phi_max': 3.14,\n",
" 'phi_min': -3.14,\n",
" 'run_number': 20251116}\n",
"\n",
" *------------------------------------------------------------------------------------* \n",
" | | \n",
" | *------------------------------------------------------------------------------* | \n",
" | | | | \n",
" | | | | \n",
" | | PPP Y Y TTTTT H H III A Welcome to the Lund Monte Carlo! | | \n",
" | | P P Y Y T H H I A A This is PYTHIA version 8.313 | | \n",
" | | PPP Y T HHHHH I AAAAA Last date of change: 10 Jan 2025 | | \n",
" | | P Y T H H I A A | | \n",
" | | P Y T H H III A A Now is 04 Dec 2025 at 17:15:11 | | \n",
" | | | | \n",
" | | Program documentation and an archive of historic versions is found on: | | \n",
" | | | | \n",
" | | https://pythia.org/ | | \n",
" | | | | \n",
" | | PYTHIA is authored by a collaboration consisting of: | | \n",
" | | | | \n",
" | | Javira Altmann, Christian Bierlich, Naomi Cooke, Nishita Desai, | | \n",
" | | Ilkka Helenius, Philip Ilten, Leif Lonnblad, Stephen Mrenna, | | \n",
" | | Christian Preuss, Torbjorn Sjostrand, and Peter Skands. | | \n",
" | | | | \n",
" | | The complete list of authors, including contact information and | | \n",
" | | affiliations, can be found on https://pythia.org/. | | \n",
" | | Problems or bugs should be reported on email at authors@pythia.org. | | \n",
" | | | | \n",
" | | The main program reference is C. Bierlich et al, | | \n",
" | | 'A comprehensive guide to the physics and usage of Pythia 8.3', | | \n",
" | | SciPost Phys. Codebases 8-r8.3 (2022) [arXiv:2203.11601 [hep-ph]] | | \n",
" | | | | \n",
" | | PYTHIA is released under the GNU General Public Licence version 2 | | \n",
" | | or later. Please respect the MCnet Guidelines for Generator Authors | | \n",
" | | and Users. | | \n",
" | | | | \n",
" | | Disclaimer: this program comes without any guarantees. | | \n",
" | | Beware of errors and use common sense when interpreting results. | | \n",
" | | | | \n",
" | | Copyright (C) 2025 Torbjorn Sjostrand | | \n",
" | | | | \n",
" | | | | \n",
" | *------------------------------------------------------------------------------* | \n",
" | | \n",
" *------------------------------------------------------------------------------------* \n",
"\n",
"\u001b[0mCxx.EventTape INFO Initialize generator...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Initializing the SingleParticle...\u001b[0m\n",
"\u001b[0mCxx.Pythia8 INFO Initialize Pythia8 generator with seed 20387086...\u001b[0m\n",
"\n",
" *------- PYTHIA Flag + Mode + Parm + Word + FVec + MVec + PVec + WVec Settings (changes only) ------------------* \n",
" | | \n",
" | Name | Now | Default Min Max | \n",
" | | | | \n",
" | ProcessLevel:all | off | on | \n",
" | Random:seed | 20387086 | -1 900000000 | \n",
" | Random:setSeed | on | off | \n",
" | | \n",
" *------- End PYTHIA Flag + Mode + Parm + Word + FVec + MVec + PVec + WVec Settings -----------------------------* \n",
"\n",
" -------- PYTHIA Particle Data Table (changed only) ------------------------------------------------------------------------------\n",
" \n",
" id name antiName spn chg col m0 mWidth mMin mMax tau0 res dec ext vis wid\n",
" no onMode bRatio meMode products \n",
"\n",
" no particle data has been changed from its default value \n",
"\n",
" -------- End PYTHIA Particle Data Table -----------------------------------------------------------------------------------------\n",
"\n",
"\u001b[0mCxx.EventTape INFO Running event 0 with eventNumber 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 63012.1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 63.0121 mass = 0.000511 pp = 63.0121 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 3.19941 py = 62.9308 pz = 3.85838e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\n",
" -------- PYTHIA Event Listing (complete event) ---------------------------------------------------------------------------------\n",
" \n",
" no id name status mothers daughters colours p_x p_y p_z e m \n",
" 0 90 (system) -11 0 0 0 0 0 0 3.199 62.931 0.000 63.012 0.001\n",
" 1 11 e- 1 0 0 0 0 0 0 3.199 62.931 0.000 63.012 0.001\n",
" Charge sum: -1.000 Momentum sum: 3.199 62.931 0.000 63.012 0.001\n",
"\n",
" -------- End PYTHIA Event Listing -----------------------------------------------------------------------------------------------\n",
"\n",
" -------- PYTHIA Event Listing (complete event) ---------------------------------------------------------------------------------\n",
" \n",
" no id name status mothers daughters colours p_x p_y p_z e m \n",
" 0 90 (system) -11 0 0 0 0 0 0 3.199 62.931 0.000 63.012 0.001\n",
" 1 11 e- 1 0 0 0 0 0 0 3.199 62.931 0.000 63.012 0.001\n",
" Charge sum: -1.000 Momentum sum: 3.199 62.931 0.000 63.012 0.001\n",
"\n",
" -------- End PYTHIA Event Listing -----------------------------------------------------------------------------------------------\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 2.22045e-16 phi = 1.52 Et = 63.0121\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Running event 1 with eventNumber 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 103247\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 103.247 mass = 0.000511 pp = 103.247 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 5.24233 py = 103.114 pz = 6.32207e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 0 phi = 1.52 Et = 103.247\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Running event 2 with eventNumber 2\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 148115\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 148.115 mass = 0.000511 pp = 148.115 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 7.52048 py = 147.924 pz = 9.06944e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 0 phi = 1.52 Et = 148.115\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Running event 3 with eventNumber 3\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 61243.8\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 61.2438 mass = 0.000511 pp = 61.2438 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 3.10962 py = 61.1648 pz = 3.7501e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 2.22045e-16 phi = 1.52 Et = 61.2438\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Running event 4 with eventNumber 4\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 88009.9\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 88.0099 mass = 0.000511 pp = 88.0099 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 4.46866 py = 87.8964 pz = 5.38905e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 0 phi = 1.52 Et = 88.0099\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Running event 5 with eventNumber 5\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 95549.3\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 95.5493 mass = 0.000511 pp = 95.5493 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 4.85147 py = 95.4261 pz = 5.85071e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 0 phi = 1.52 Et = 95.5493\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Running event 6 with eventNumber 6\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 114742\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 114.742 mass = 0.000511 pp = 114.742 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 5.82597 py = 114.594 pz = 7.02593e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 0 phi = 1.52 Et = 114.742\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Running event 7 with eventNumber 7\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 51921.6\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 51.9216 mass = 0.000511 pp = 51.9216 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 2.63629 py = 51.8546 pz = 3.17928e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 0 phi = 1.52 Et = 51.9216\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Running event 8 with eventNumber 8\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 98189.7\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 98.1897 mass = 0.000511 pp = 98.1897 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 4.98553 py = 98.0631 pz = 6.01239e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 0 phi = 1.52 Et = 98.1897\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Running event 9 with eventNumber 9\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Start particle gun...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle: eta = 0 phi = 1.52 energy = 84175.3\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Fill particle with ID 11\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Storing particle into the event record...\u001b[0m\n",
"\u001b[0mCxx.Electron INFO id = 11 energy = 84.1753 mass = 0.000511 pp = 84.1753 eta = 0 phi = 1.52 at rest = 0 has lifetime = 0\u001b[0m\n",
"\u001b[0mCxx.Electron INFO px = 4.27396 py = 84.0668 pz = 5.15425e-15\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Particle stored with i = 1\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Adding particle: id = 11 eta = 0 phi = 1.52 Et = 84.1753\u001b[0m\n",
"\u001b[0mCxx.EventTape INFO Finalize the Event generator.\u001b[0m\n",
"\u001b[0mCxx.Electron INFO Finalize the SingleParticle Event.\u001b[0m\n",
"\u001b[0mCxx.Pythia8 INFO Finalize the Pythia8 Event.\u001b[0m\n",
"\n",
" *------- PYTHIA Error and Warning Messages Statistics ----------------------------------------------------------* \n",
" | | \n",
" | times message | \n",
" | | \n",
" | 0 no errors or warnings to report | \n",
" | | \n",
" *------- End PYTHIA Error and Warning Messages Statistics ------------------------------------------------------* \n",
"\u001b[0mCxx.StoreGate INFO Saving StoreGate...\u001b[0m\n",
"\u001b[0mCxx.StoreGate INFO Writing all root objects into the file and delete all decorators\u001b[0m\n",
"\u001b[0mCxx.StoreGate INFO Decorators will no be saved into the root file.\u001b[0m\n",
"hadd -f Electron.EVT.root Electron.EVT.0.root\n",
"Module libc not found.\n",
"hadd Target file: Electron.EVT.root\n",
"hadd compression setting for all output: 101\n",
"hadd Source file 1: Electron.EVT.0.root\n",
"hadd Target path: Electron.EVT.root:/\n",
"Warning in :\n",
"Warning in : Unchecked status codes encountered during the job\n",
"Warning in : Number of unchecked successes: 3\n",
"Warning in : To fail on an unchecked code, call StatusCode::enableFailure() at the job's start\n",
"Warning in :\n"
]
}
],
"source": [
"!gen_single.py --job-file config.json --nov 10 -nt 1 -o Electron.EVT.root -p Electron"
]
}
],
"metadata": {
"kernelspec": {
"display_name": "Python 3 (ipykernel)",
"language": "python",
"name": "python3"
},
"language_info": {
"codemirror_mode": {
"name": "ipython",
"version": 3
},
"file_extension": ".py",
"mimetype": "text/x-python",
"name": "python",
"nbconvert_exporter": "python",
"pygments_lexer": "ipython3",
"version": "3.10.12"
}
},
"nbformat": 4,
"nbformat_minor": 5
}