Stefano Roberto Soleti, 25 May 2021
In this notebook we will analyze the charge distributions in Module0 data and in the simulation. The simulation sample was produced using CORSIKA as cosmic-ray producer, edep-sim for the simulation of the passage of the particles in the matter, and larnd-sim for the detector simulation.
Here we plot (1) the number of hits as a function of the track length and (2) the number of hits for segments of different lengths along the reconstructed tracks. The following cuts are applied:
Here we plot the dQ/dx for segments of different lengths along the reconstructed tracks. The following cuts are applied:
The charge has been corrected by the gain difference (4 mV/$10^3$e$^-$ for simulation and 4.08 mV/$10^3$e$^-$ for data) and by the electron lifetime ($2.2~\mathrm{\mu}$s for simulation and polynomial fit provided by CSU for data).
The $dQ/dx$ is plotted as a function of $\cos\theta$ and $\phi$ between the track and the anode plane. The charge distributions for each bin have been fitted with a Moyal distribution convoluted with a Gaussian.
The full $dQ/dx$ distributions for each segment length are plotted here and fitted with a Moyal distribution convoluted with a Gaussian. The histograms are normalized to their integral.
Since there are gaps in the detector due to disabled channels, here we group the tracks according to their relative angles and distance. If $\cos\theta > 0.94$ and the distance between two endpoints is $< 20~$cm then we group the two tracks together.
Then, groups of track that have an endpoint farther than $10~$cm from the detector borders are selected as stopping particle candidates.
Finally, we plot the distribution of the $dQ/dx$ as a function of the residual range for segments of $2~$cm. The theoretical curves for protons and muons correspond to the most probable values of the Landau-Vavilov-Bichsel distribution and are corrected by the Birks recombination factor, using the parameters obtained by ICARUS.
61283.837301145584
2.06861514705665
