- printf ("I'm after Map \n");
-
-// Get AliRun object from file or create it if not on file
-
- if (!gAlice) {
- gAlice = (AliRun*)file->Get("gAlice");
- if (gAlice) printf("AliRun object found on file\n");
- if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
- }
- printf ("I'm after gAlice \n");
+ if (gClassTable->GetID("AliRun") < 0) {
+ gInterpreter->ExecuteMacro("loadlibs.C");
+ }
+ else {
+ if(gAlice){
+ delete AliRunLoader::Instance();
+ delete gAlice;
+ gAlice=0;
+ }
+ }
+
+
+ // Get geometry
+ AliGeomManager::LoadGeometry("geometry.root");
+
+ // Connect the Root Galice file containing Geometry, Kine and Hits
+ AliRunLoader* rl = AliRunLoader::Open(inFile.Data());
+ if (rl == 0x0)
+ {
+ Error("AliITSHits2FastRecPoints.C","Can not open session RL=NULL");
+ return;
+ }
+
+ Int_t retval = rl->LoadgAlice();
+ if (retval)
+ {
+ Error("AliITSHits2FastRecPoints.C","LoadgAlice returned error");
+ delete rl;
+ return;
+ }
+ gAlice=rl->GetAliRun();
+ rl->LoadHeader();
+ retval = rl->LoadKinematics();
+ if (retval)
+ {
+ Error("AliITSHits2FastRecPoints.C","LoadKinematics returned error");
+ delete rl;
+ return;
+ }
+
+ AliITSInitGeometry initgeom;
+ AliITSgeom *geom = initgeom.CreateAliITSgeom();
+ printf("Geometry name: %s \n",(initgeom.GetGeometryName()).Data());
+
+ AliITSLoader* gime = (AliITSLoader*)rl->GetLoader("ITSLoader");
+ if (gime == 0x0)
+ {
+ ::Error("AliITSHits2FastRecPoints.C","can not get ITS loader");
+ delete rl;
+ return;
+ }
+ gime->SetITSgeom(geom);
+ retval = gime->LoadHits("read");
+ if (retval)
+ {
+ ::Error("AliITSHits2FastRecPoints.C","LoadHits returned error");
+ delete rl;
+ return;
+ }
+ gime->SetRecPointsFileName("ITS.FastRecPoints.root");
+ retval = gime->LoadRecPoints("update");
+ if (retval)
+ {
+ ::Error("AliITSHits2FastRecPoints.C","LoadRecPoints returned error");
+ delete rl;
+ return;
+ }
+