...
...
The macro doing the main job is EventFilter.C, which was written by W. Focke H. Kelly and A. W. Borgland. This macro creates a new digi.root OR recon.root OR root OR mc.root file with a sub-sample of the initial events, which are selected according to a cut in the "merit" or "svac" variables. This macro runs very nicely, but has two small drawbacks:
- If you want to process the 3 files (digi, recon and mc) which is usually the case, you have to run the macro 3 times
- You have to be sure you are loading the libraries from the right software package and release version before running it
The macro FilterRootFiles.C is a very simple macro which makes the event filtering a bit easier. With this macro, the user can specify a software package as argument (like "BeamtestRelease-v3r0907p2" or "EngineeringModel-v6r070329p16") and the macro loads the appropriate root libraries (using the macro LoadLibraries.C) and then executes EventFilter.C for the 3 types of files (in case they are specified as arguments) sequencially.
An example of running the macro FilterRootFiles.C is given by macro FilterEvents_BeamTest_CR_CERN322.C, which selects is the macro that defines the software package to be used, the input and output filenames and the selection cut to be applied to the events. In partituclar this macro uses the software package "BeamtestRelease-v3r0907p2" to select events from beam test run 700000322 according to the selection criteria:
char* cuts = "TMath::Abs((acos(Tkr1ZDir)*(180./3.14159))-180.) < 5 && CalEnergyRaw > 10";
You can easily modify this macro to fit your needs.
3 - Setting the FRED GLAST panel to visualize the desired data
...