Differences

This shows you the differences between two versions of the page.

--- software:g4simulation:analysis [2019/03/07 11:15] – created himuser
+++ software:g4simulation:analysis [2022/01/14 18:44] (current) – external edit 127.0.0.1
@@ Line 1: / Line 1: @@
 ====== Analysis ======
-An example how to read the .root files can be found in rootmacros/test.C.
+There are already several files which help you to analyse root files from primasoft:
-This test.C file can be executed with:
-<code>
-/primasoft/geant4$ root -l rootmacros/test.C+
-</code>
-If executed, a histogram will be displayed with Root.
+  - PrimaSimAnalysis (can be found here: /geant4/src/PRIMASimAnalysis.cc, /geant4/include/PRIMASimAnalysis.hh)
+  - Config file for Analysis (examples can be found here: /geant4/sim_ana_utilities/ana_default.conf)
+  - Rootmacros (examples can be found here: /geant4/rootmacros)
-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",calInUse,taggerInUse)).
+This is what they do:
+== PrimaSimAnalysis ==
+  * 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.
+== Rootmacro ==
+  * an example how to read the .root files can be found in rootmacros/test.C
+  * execute a rootmacro (here: test.C): <code>
+/primasoft/geant4$ root -l rootmacros/test.C+
+</code>
+  * the root files include information about the primary particles, the deposited energy in the calorimeter and in the tagger.
+  * 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: <code> c.OpenTreeFiles("test_bwec.root","sim_ana_utilities/ana_default.conf");  </code>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 examples
+  * event_display.C the energy signal for each crystal for one event (usage: ''geant4$ root -l 'rootmacros/event_display.C("path/file.root",eventnumber)''')