Brief overview of the codes


  • compiled code – usually run 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


  • interactive code based on a single class, "GcrHistos" (the same used by gcrMake to fill histograms)
  • input: ROOT and asci files produced by gcrMake
  • output: histograms arranged in panel, fit results, calibration accuracy
  • fast (loop on events is already done by gcrMake)
  • easy to use: one can zoom in any particular channel, layer, tower etc.
  • modular: one can add different data sets
    root [0] .L gcrHistos.cxx                                      
    root [1] g=new GcrHistos("gcr_nomSciOps_noSkirtCno_noCno");// opens the file and reads 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 of the stufy on CNO events vs. ACD tiles
    root [5] g->ShowMs(); // shows the distributions characterizing the multiple scattering of protons
    root [6] g->SumPeaks(); // sums up log histograms to build layer, super-layer, tower, and whole CAL histograms 
    root [7] g->FitAllPeaks(1); // fits all carbon ("1") peaks with Landau distribution
    root [8] g->ShowPeakPanel(1, "CAL"); // displays carbon ("1") fit histogram for the whole CAL
    root [9] g->ShowPeakPanel(1, "LAY", 0); // displays carbon ("1") fit histograms for the 16 layers in super-layer 0
    root [10] g->ShowPeakPanel(1, "LAY", 7); // displays carbon ("1") fit histograms for the 16 layers in super-layer 7

    root [11] g->ShowPeakPanel(1, "TOW", 5); // displays carbon ("1") fit histograms for the 8 layers in tower 5
  • No labels