10 Int_t particle_type=0;
14 void RICHdigit (Int_t evNumber1=0,Int_t evNumber2=0)
16 /////////////////////////////////////////////////////////////////////////
17 // This macro is a small example of a ROOT macro
18 // illustrating how to read the output of GALICE
19 // and do some analysis.
21 /////////////////////////////////////////////////////////////////////////
23 // Dynamically link some shared libs
26 if (gClassTable->GetID("AliRun") < 0) {
27 gROOT->LoadMacro("loadlibs.C");
32 // Connect the Root Galice file containing Geometry, Kine and Hits
34 TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
35 if (file) file->Close();
36 file = new TFile("galice.root","UPDATE");
38 // Get AliRun object from file or create it if not on file
42 if (gClassTable->GetID("AliRun") < 0) {
43 gROOT->LoadMacro("loadlibs.C");
54 gAlice = (AliRun*)file->Get("gAlice");
55 if (gAlice) printf("AliRun object found on file\n");
56 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
59 gAlice = (AliRun*)file->Get("gAlice");
60 if (gAlice) printf("AliRun object found on file\n");
61 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
64 AliRICH *RICH = (AliRICH*) gAlice->GetDetector("RICH");
66 AliRICHChamber* iChamber;
68 printf("Generating tresholds...\n");
70 for(Int_t i=0;i<7;i++)
72 iChamber = &(RICH->Chamber(i));
73 iChamber->GenerateTresholds();
80 for (int nev=0; nev<= evNumber2; nev++) {
81 Int_t nparticles = gAlice->GetEvent(nev);
82 cout <<endl<< "Processing event:" <<nev<<endl;
83 cout << "Particles :" <<nparticles<<endl;
84 if (nev < evNumber1) continue;
85 if (nparticles <= 0) return;
86 if (RICH) RICH->Digitise(nev, particle_type);
88 //sprintf(hname,"TreeD%d",nev);
89 //gAlice->TreeD()->Write(hname);
90 //gAlice->TreeD()->Reset();
95 printf("\nEnd of Macro *************************************\n");