Versions Compared

Key

  • This line was added.
  • This line was removed.
  • Formatting was changed.
  • Can GCRcalib help to calibrate ACD HI charge scale? (Fred)
    Following the presentation of last week, I analysed 50 runs of opssim2 data. The collection time is 70h (effective duration of 60h).
    For each event where GCRCalib found an even value for the atomic number in the range 6-14 (i.e. C, O, Ne, Mg , or Si), I loop on the AcdTrkHitPoca list in ACD recon. For those hits which have a suitable DOCA to Tkr1 track (doca>-50mm, i.e. Tkr1 track within 50mm of the tile), I get the number of MIPs from the AcdHit collection. This is a very slow analysis (8h30 to analyse 50 runs!). The following plot show the obtained MIP spectrum, combining all ACDtiles: the vertical lines

I modified the code I used to analyse the GCR events in opssim2 data, and added the reading of digi events (now it reads merit, gcr and digi files). I obtained the spectrum of each PM of each tile when GCRcalib is able to infer a value for the atomic number Z. The plots below show the low-range and high-range spectrum (top and medium plots, resp.) for all species, and the high-range spectrum for events tagged as carbons (bottom plots). This is for the PMT A of the central tile at the top of the ACD (left panel) and of the central tile (2nd row, 3rd column) of an ACD side (right panel).

Image Removed

...

  • red vertical lines stands for 36 (C), 64 (O), 100 (Ne), 144 (Mg) and 196 (Si).
    Image Added
    Only the carbon peak is visible. The rate of events above 20 MIPs is 0.06Hz, thus really lower than the ~1Hz of GCRCalib throughput.