Align AOD creation with latest TOF (#51)

- Add more explanatory comments
- Add vertex spread to AOD

Author: njacazio

examples/aod/createO2tables.C

@@ -29,12 +29,13 @@ R__LOAD_LIBRARY(libDelphesO2)
 // Detector parameters
 const double Bz = 0.2; // [T]
 // TOF
-const double tof_radius = 100.; // [cm]
-const double tof_length = 200.; // [cm]
-const double tof_sigmat = 0.02; // [ns]
+const double tof_radius = 100.; // [cm] Radius of the TOF detector (used to compute acceptance)
+const double tof_length = 200.; // [cm] Length of the TOF detector (used to compute acceptance)
+const double tof_sigmat = 0.02; // [ns] Resolution of the TOF detector
+const double tof_sigmat0 = 0.2; // [ns] Time spread of the vertex
 // RICH
-const double rich_radius = 100.; // [cm]
-const double rich_length = 200.; // [cm]
+const double rich_radius = 100.; // [cm] Radius of the RICH detector (used to compute acceptance)
+const double rich_length = 200.; // [cm] Length of the RICH detector  (used to compute acceptance)
 const double rich_index = 1.03;
 const double rich_radiator_length = 2.;
 const double rich_efficiency = 0.4;
@@ -152,7 +153,7 @@ void createO2tables(const char* inputFile = "delphes.root",
 
   // TOF layer
   o2::delphes::TOFLayer toflayer;
-  toflayer.setup(tof_radius, tof_length, tof_sigmat);
+  toflayer.setup(tof_radius, tof_length, tof_sigmat, tof_sigmat0);
   // RICH layer
   o2::delphes::RICHdetector richdetector;
   richdetector.setup(rich_radius, rich_length);

examples/scripts/createO2tables.py

@@ -6,7 +6,6 @@
 
 import argparse
 import configparser
-import subprocess
 import os
 import shutil
 import multiprocessing
@@ -124,6 +123,7 @@ def opt(entry, require=True):
     # detector configuration
     bField = opt("bField")
     sigmaT = opt("sigmaT")
+    sigmaT0 = opt("sigmaT0")
     radius = opt("radius")
     length = opt("length")
 
@@ -175,10 +175,11 @@ def check_duplicate(option_name):
     msg("  nruns   = ", nruns)
     msg("  nevents = ", nevents)
     msg(" --- with detector configuration", color=bcolors.HEADER)
-    msg("  bField  = ", bField, " [kG] ")
-    msg("  sigmaT  = ", sigmaT, " [ns] ")
-    msg("  radius  = ", radius, " [cm] ")
-    msg("  length  = ", length, " [cm] ")
+    msg("  bField  = ", bField, "\t[kG]")
+    msg("  sigmaT  = ", sigmaT, "\t[ns]")
+    msg("  sigmaT0 = ", sigmaT0, "\t[ns]")
+    msg("  radius  = ", radius, "\t[cm]")
+    msg("  length  = ", length, "\t[cm]")
     msg("  LUT     = ", lut_tag)
     msg("  etaMax  = ", etaMax)
 
@@ -225,6 +226,8 @@ def set_config(config_file, config, value):
                f"{sigmaT}""e\-9/")
     set_config("createO2tables.C",
                "const double tof_sigmat =", f"{sigmaT}""\;/")
+    set_config("createO2tables.C",
+               "const double tof_sigmat0 =", f"{sigmaT0}""\;/")
 
     run_list = range(nruns)
 

examples/scripts/default_configfile.ini

@@ -7,8 +7,12 @@
 # magnetic field [kG]
 bfield = 5.
 
+## TOF
 # time resolution [ns]
 sigmat = 0.020
+# vertex spread [ns]
+sigmat0 = 0.20
+## RICH
 
 # radius [cm] used to compute the acceptance in eta (together with the detector length)
 radius = 100.