10 Int_t particle_type=0;
14 void RICHdigit (Int_t evNumber1=0,Int_t evNumber2=0, Int_t merging)
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");
36 // Connect the Root Galice file containing Geometry, Kine and Hits
38 TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
39 if (file) file->Close();
40 file = new TFile("galice.root","UPDATE");
42 // Get AliRun object from file or create it if not on file
46 if (gClassTable->GetID("AliRun") < 0) {
47 gROOT->LoadMacro("loadlibs.C");
53 gAlice = (AliRun*)file->Get("gAlice");
54 if (gAlice) printf("AliRun object found on file\n");
55 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
58 gAlice = (AliRun*)file->Get("gAlice");
59 if (gAlice) printf("AliRun object found on file\n");
60 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
63 AliRICH *RICH = (AliRICH*) gAlice->GetDetector("RICH");
66 printf("Merging is ON\n");
68 printf("Merging is OFF\n");
70 // Creation of merger object
71 AliRICHMerger* merger = new AliRICHMerger();
74 merger->SetMode(merging);
75 merger->SetSignalEventNumber(0);
76 merger->SetBackgroundEventNumber(0);
77 merger->SetBackgroundFileName("bg.root");
80 RICH->SetMerger(merger);
86 for (int nev=0; nev<= evNumber2; nev++) {
87 Int_t nparticles = gAlice->GetEvent(nev);
88 cout <<endl<< "Processing event:" <<nev<<endl;
89 cout << "Particles :" <<nparticles<<endl;
90 if (nev < evNumber1) continue;
91 if (nparticles <= 0) return;
94 //gAlice->MakeTree("D");
95 //RICH->MakeBranch("D");
96 //RICH->Digitise(nev, particle_type);
97 //gAlice->SDigits2Digits("RICH");
98 //gAlice->Tree2Tree("D");
99 RICH->MakeBranch("D");
100 RICH->SDigits2Digits(nev, particle_type);
106 printf("\nEnd of Macro *************************************\n");