Hi, I have added a new module to brat/ana: BrEffPackage.h, BrEffPackage.cxx which is a very sketchy and clumsily programmed "efficiency module" which reads raw data and geant events in parallel, assuming the following input: InputNode |---BrEventNode GeantDataNode (normal geant data) | |---BrEventNode RawDataNode (normal raw data) It handles digitization, event adding on the TPC sequence level and track reconstruction (normal chain and directly from Geant track/hits to generate "ideal" tracks.) It can handle single TPCs or spectrometer arms - only MRS properly implemented for the time being. It creates the following output: OutputNode |--BrEventNode GeantTrackNode (tracks direct from geant hits and tracks) | |--BrEventNode RawTrackNode (tracks from normal raw event) | |--BrEventNode SingleTrackNode (tracks from digitized geant) | |--BrEventNode AddedTrackNode (tracks from added event) each node with a set of track BrDataTables. A still rather messy application program is found on: ~trine/brat_test/effpackage_test/EffPackageTester.cxx A prototype macro reading the output file, checking overlap between geant tracks and reconstructed tracks and making a few histograms is found on: ~trine/brat_test/effpackage_test/ReadEffResults_MRS.C Truls and I will (hopefully) try to do some more work on this during the coming days - in particular Truls has some good ideas for developing the package and making it more elegant and user-friendly. Cheers, Trine
This archive was generated by hypermail 2b29 : Wed Jun 20 2001 - 12:48:07 EDT