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:

• The deposited energy in every crystal is stored: “calEdep”
• The signal energy in every crystal is simulated: “calSignal”
• Hits are sorted into clusters: “calClusters”, where the individual hits are stored and the total cluster energy is accessible
• The momentum direction and the energy of the initial particle is stored: “P”
• Single crystal threshold is applied. The value can 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 setup can be chosen in the config file, default is to use all 786 crystals. These are the options:
  • 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)
• Scintillation hits is still under developement …..
• 3 functions are left empty (BeginOfRun(), DoEvent(), EndOfRun()). These need to be defined in the rootmacro.

• execute a rootmacro (here: test.C):

  root -l rootmacros/test.C+

• in the rootmacros the functions BeginOfRun(), DoEvent() and EndOfRun() need to be defined, please have a look at the examples
• in OpenTreeFiles() you can define the root file which will be analysed and you can optionally choose a config file:

   c.OpenTreeFiles("test_bwec.root","sim_ana_utilities/ana_default.conf");  

  If the path is relative and not absolute for the root and config file, pay attention that the relative path changes when you execute the rootmacro from a different folder.

• test.C is a basic example
• energy_res.C and acceptance.C are more advanced exaples

An example how to read the .root files can be found in rootmacros/test.C. This test.C file can be executed with:

/primasoft/geant4$ root -l rootmacros/test.C+

If executed, a histogram will be displayed with Root.

For the analysis some pre-calculations are made in include/PRIMASimAnalysis and src/PRIMASimAnalysis. These are imported in test.C and the functions BeginOfRun(), DoEvent() and EndOfRun can be adjusted in MyEventAnalysis. The root files include information about the primary particles, the deposited energy in the calorimeter and in the tagger. For the calorimeter there are several functions implemented in SimAnalysis, for example a ClusterFinder.

When reading the root files with OpenTreeFiles(), flags for the calorimeter and the tagger can be set, which enable their use respectively (the syntax is: OpenTreeFiles(“filename.root”,“configfile.conf”)).