/[PAMELA software]/PamelaLevel2/doc/examples/README.txt
ViewVC logotype

Diff of /PamelaLevel2/doc/examples/README.txt

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1 by pam-fi, Fri Jun 16 16:43:55 2006 UTC revision 1.4 by pam-fi, Wed Jan 3 13:28:50 2007 UTC
# Line 1  Line 1 
1  This directory contains some very simple macros, just to show how to use the PamLevel2 class.  This directory contains some very simple macros, just to show how to use the PamLevel2 class.
2  The macros need as input the name of a Level2 .root file produced by DarthVader.  The macros need as input the name of a Level2 .root file produced by DarthVader
3    (example1 and example2) or a list of files (see example3).
4    
5  Class implementation is under development and some problems still need to be fixed.  example1 - Example to integrate all detector info and solve the y-view tracks ambiguity
6  However this simple examples should allow to start working with Pamela data and to test the  example2 - Example to get pamela tracks and integrate the trajectory in the apparatus
7  general environment.  example3 - Example to loop over events, reading more than one file
8    Loop     - Example to loop over events, apply some selection cuts (My-Selection.cpp) and filling an
9               output file with histograms (My-Histos.cpp) and/or trees.
10    
11  In order to run the scripts:  In order to run the scripts:
12    
# Line 17  In order to run the scripts: Line 20  In order to run the scripts:
20    
21  [] .x load.C  [] .x load.C
22  [] .L example1.C  [] .L example1.C
23  [] example1("level2-data-path/level2-file.root")  [] example1("level2-data-path/level2-file, execute the command.root")
24  Tracker      : set branch address TrkLevel2  Tracker      : set branch address TrkLevel2
25  Calorimeter  : set branch address CaloLevel2  Calorimeter  : set branch address CaloLevel2
26  ToF          : set branch address ToFLevel2  ToF          : set branch address ToFLevel2
# Line 30  Anticounter  : set branch address AcLeve Line 33  Anticounter  : set branch address AcLeve
33    
34   Done   Done
35    
36  4) The macro creates some histograms, which you can list and plot as follows:  4) Macros example1 and example2 produce some histograms. Example3 just loop over events and tracks and
37    counts time. It is usefull to check processing time.
38    
39  [].ls  
40  TROOT*          Rint    The ROOT of EVERYTHING  5) To create Loop.exe execute the command (you need the PAMELA environment variables):
41   OBJ: TH1F      qtot    Total energy in the calorimeter : 0 at: 0xb3f7400  
42   OBJ: TH1F      npaddle Total number of hit ToF paddles : 0 at: 0xb3f8768  > make Loop.exe
43   OBJ: TH1F      ntr     Number of fitted tracks : 0 at: 0xb3f8b60  
44   OBJ: TH1F      rig     Particle rigidity : 0 at: 0xb3f8f30  The executable is create by compiling and linking the files:
45   OBJ: TH1F      resxs   Spatial residual (X) on the 1^ calo plane (sorted tracks) : 0 at: 0xb3f94a0  
46   OBJ: TH1F      resys   Spatial residual (Y) on the 1^ calo plane (sorted tracks) : 0 at: 0xb3f9a00  Loop.cpp
47   OBJ: TH1F      resxi   Spatial residual (X) on the 1^ calo plane (image tracks) : 0 at: 0xb3f9f70  My-Selection.cpp
48   OBJ: TH1F      resyi   Spatial residual (Y) on the 1^ calo plane (image tracks) : 0 at: 0xb3fa4e0  My-Histos.cpp
49  []rig->Draw()  
50    You can use My-Selection.cpp and My-Histos.cpp as a template to inplement your own
51    selection cuts and analysis.
52    More generally, you can compile and link Loop.cpp with your own source files (you need
53    to edit also the Makefile), where you should implement at least the following functions:
54    
55    bool Select(PamLevel2*);
56    void CreateHistos(TFile*);
57    void FillHistos(PamLevel2*);
58    void SaveHistos(TFile*);
59    
60    To know how to use Loop.exe, issue:
61    ./Loop.exe --help
62    
63    (If you execute "make Loop-0.exe" the old example code, previously posted on the cvs,
64    is compiled and linked. This is not a very general example, better use the new one)
Legend:
Removed from v.1.1  
changed lines
  Added in v.1.4
  ViewVC Help
Powered by ViewVC 1.1.23