There are already several files which help you to analyse root files from primasoft:

1. PrimaSimAnalysis (can be found here: /geant4/src/PRIMASimAnalysis.cc, /geant4/include/PRIMASimAnalysis.hh)
2. Config file for Analysis (examples can be found here: /geant4/sim_ana_utilities/ana_default.conf)
3. Rootmacros (examples can be found here: /geant4/rootmacros)

This is what they do:

• Sorts hits into clusters
• Single crystal threshold is applied. The value van be set in the config file, default is 2 MeV. This is analogue the to hit detection threshold in a measurement.
• Different crystal-setups can be used. The setp can be chosen in the config file, default is to use all 786 crystals.
  • all crystals are used (768)
  • the inner 640 crystals are used symmetrically (640)
  • the inner 640 crystals are used in an oval shape (640_oval)
  • the inner 560 crystals are used symmetrically (560)
• ff
• 3 functions are left empty (BeginOfRun(), DoEvent(), EndOfRun()). These need to be defined in the rootmacro

• execute a rootmacro:

  root -l rootmacros/test.C+