Brief overview of the codes
gcrMake
- compiled code on noric (or SLAC farm)
- input: list of gcr and merit file locations in XROOTD
- output: root file with filled histograms and text file with event rates
- easy to implement any other analysis (e.g., CNO events through ACD tiles)
- same code runs on data AND simulations
gcrHistos
- interactive code based on a single class (GcrHistos), used by gcrMake
- input: ROOT and asci files produced by gcrMake
...
- output: histograms arranged in panel, fit results, calibration accuracy
- fast (loop on events is already done)
- easy to use: we can zoom on any particular channel, layer, tower etc.
- modular: we can add different data sets
Code Block root [0] .L gcrHistos.cxx root [1] g=new GcrHistos("titi/gcr_nomSciOps_noSkirtCno_noCno");// opens the file and read the histograms <><><><> opening filename=titi/gcr_nomSciOps_noSkirtCno_noCno.root <><><><> reading histograms <><><><> BuildRates <><><><> BuildEvents <><><><> BuildHits <><><><> BuildEngine4TopAcd <><><><> BuildMs <><><><> BuildPeaks (class GcrHistos*)0x97efa88 root [2] g-
...
>ShowEvents(); //shows event based distributions:nb of TKR triggering towers, OBF status word, nb of hits/event root [3] g->ShowHits(); //shows hit based distributions:path-length, raw and corrected energy spectra root [4] g->ShowEngine4TopAcd(); // shows the results on CNO events passing through ACD tiles root [5] g->ShowMs(); // shows distributions characterizing multiple scattering of protons root [6] g->SumPeaks(); // sums up log histograms to build layer, super-layer, tower, and whole CAL distribution root [7] g->FitAllPeaks(1); // fit all peaks with Landau distribution! root [
...
8] g->ShowPeakPanel(1, "CAL"); // displays carbon ("1") fit histogram for the whole CAL
Code Block root [
...
9] g->ShowPeakPanel(1, "LAY", 0); // displays carbon ("1") fit histogram for the 16 layers in super-layer 0 root [
...
10] g->ShowPeakPanel(1, "LAY", 7); // displays carbon ("1") fit histogram for the 16 layers in super-layer 7
Code Block root [
...
11] g->ShowPeakPanel(1, "TOW", 5);
...
...
...
// displays carbon ("1") fit histogram for the 16 layers in tower