1 RICHdetect (Int_t evNumber1=0,Int_t evNumber2=0) {
2 // Dynamically link some shared libs
4 if (gClassTable->GetID("AliRun") < 0) {
5 gROOT->LoadMacro("loadlibs.C");
15 // Connect the Root Galice file containing Geometry, Kine and Hits
17 TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
18 if (file) file->Close();
19 file = new TFile("galice.root","UPDATE");
22 //printf ("I'm after Map \n");
24 // Get AliRun object from file or create it if not on file
27 gAlice = (AliRun*)file->Get("gAlice");
28 if (gAlice) printf("AliRun object found on file\n");
29 if (!gAlice) gAlice = new AliRun("gAlice","RICH reconstruction program");
32 gAlice = (AliRun*)file->Get("gAlice");
33 if (gAlice) printf("AliRun object found on file\n");
34 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
37 //printf ("I'm after gAlice \n");
39 AliRICH *RICH = (AliRICH*) gAlice->GetDetector("RICH");
41 // Create Recontruction algorithm object
42 AliRICHDetect *detect = new AliRICHDetect("RICH reconstruction algorithm","Reconstruction");
47 for (int nev=0; nev<= evNumber2; nev++) {
48 Int_t nparticles = gAlice->GetEvent(nev);
49 printf("\nProcessing event: %d\n",nev);
50 printf("\nParticles : %d\n",nparticles);
51 if (nev < evNumber1) continue;
52 if (nparticles <= 0) return;
53 if (RICH) detect->Detect();
55 sprintf(hname,"TreeR%d",nev);
56 gAlice->TreeR()->Write(hname);
57 gAlice->TreeR()->Reset();
61 printf("\nEnd of Macro *************************************\n");