00001 #include "MapsTrackManager.hh"
00002 #include "MapsTrack.hh"
00003 #include "MapsSensor.hh"
00004 #include "MapsException.hh"
00005 #include <cmath>
00006 #include <vector>
00007 #include "TRandom2.h"
00008 #include "TFile.h"
00009 #include "TH1F.h"
00010 int main(int argc, const char **argv) {
00011
00012 for(int d(0); d < argc; d++) {
00013 std::cout << argv[d] << "\n";
00014 }
00015 if(argc != 3) {
00016 std::cout << "Usage: Rewriter <input file> <outputfile>\n";
00017 return 0;
00018 }
00019
00020 MapsTrackManager mtm2;
00021 std::cout << "Recreating from root file...\n";
00022 char input[100];
00023 strcpy(input, argv[1]);
00024
00025 mtm2.recreateFromRootFile(input);
00026 std::cout << "Rewriting to new root file...\n";
00027 std::vector<MapsSensor*> s(mtm2.getSensors());
00028 for(std::vector<MapsSensor*>::iterator it = s.begin(); it != s.end(); ++it) {
00029 MapsSensor* sensor = *it;
00030 if(sensor->id() == 2 || sensor->id() == 8)
00031 sensor->setPhi(3.14159265358979323846);
00032
00033 }
00034
00035
00036
00037 std::vector<MapsTrack*> tracks(mtm2.getTracks());
00038 for(std::vector<MapsTrack*>::iterator it = tracks.begin(); it != tracks.end(); ++it)
00039 (*it)->eraseGlobalHits();
00040
00041 char output[100];
00042 strcpy(output, argv[2]);
00043 mtm2.exportToRootFile(output);
00044
00045
00046 std::cout << "Done.\n";
00047 return 0;
00048 }
00049
00050