You are viewing an old version of this page. View the current version.

Compare with Current View Page History

« Previous Version 3 Next »

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 
    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

    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


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