TODO for tracking project

[doloresgarcia]

In relation to the MC particles

Create a modular script with good documentation and one function per section

1. Energy:
   1. What is the energy of the particles actually hitting the DC (produce with more stats, add in the title the number of BX files used to produce it) + 1 plot for MCs with 1 hit + 1 plot for MCs > 1 hit.
   2. what is the energy of the particles actually hitting the DC - low values
   3. what is the energy of the particles actually hitting the DC when they come from photons, check if it's parent is a photon (https://github.com/BrieucF/FCC_scripts/blob/d379288807a3dad9b7a67d87a646bb25abf6f65e/background_studies/dc_occupancy.py#L365)
   4. what is the pt of the particles actually hitting the DC and with vertex radius below 10 cm
   5. [not produce for now] what is the de/DX of the background particles/muons (not clear for now)?
2. PDG
   1. what is the particle pdg/id hitting DC? (PDG ID of bkg particles hitting the DCH) (montecarlo guide to pdg, https://pdg.lbl.gov/2007/reviews/montecarlorpp.pdf) Bar plot
   3. most particles leaving signals in DCH are electrons or positrons (because photon kick out electron that then leave a signal), how many of this were produced by photons?
   4. are the particles hitting dc primaries or secondaries? Filled once per MCParticle, Is the particle a primary?
3. Location
   1. Where are the particle leading to hits in the DC coming from r (mm) vs z (mm)?
   2. where are the particles leading to hits in the DC coming from in x vs y?
4. Occupancy:
   1. percentage_of_fired_cells
   2. how many cells are fired by a bkg particle? (different cell IDs per MC)
   3. how many cells are fired by a bkg particle in log scale?
   4. Average occupancy per layer (number of layer vs %occupancy) check the example code script (if not cellID_unique_identifier in dict_cellID_nHits.keys(): # the cell was not fired yet) Example code https://github.com/BrieucF/FCC_scripts/blob/main/background_studies/dc_occupancy.py

In relation to the hit distribution

1. Average number of hits per event as a function of z and phi (Azimutal angle):
   1. mean for 1 background file
   2. Mean with the overlayed files and use the flag to detect background hits (separate between signal and background)
2. DCH simHit position (x vs y)
   Example (from Nate's talk)