3 void MUONdigit (Int_t evNumber1=0,Int_t evNumber2=0,Int_t nsignal =25)
5 /////////////////////////////////////////////////////////////////////////
6 // This macro is a small example of a ROOT macro
7 // illustrating how to read the output of GALICE
8 // and do some analysis.
10 /////////////////////////////////////////////////////////////////////////
12 // Dynamically link some shared libs
14 if (gClassTable->GetID("AliRun") < 0) {
15 gROOT->LoadMacro("loadlibs.C");
20 // Connect the Root Galice file containing Geometry, Kine and Hits
22 TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root");
23 if (file) file->Close();
24 file = new TFile("galice.root","UPDATE");
27 printf ("I'm after Map \n");
29 // Get AliRun object from file or create it if not on file
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");
36 printf ("I'm after gAlice \n");
38 AliMUON *MUON = (AliMUON*) gAlice->GetModule("MUON");
44 for (int nev=0; nev<= evNumber2; nev++) {
45 Int_t nparticles = gAlice->GetEvent(nev);
46 cout << "nev " <<nev<<endl;
47 cout << "nparticles " <<nparticles<<endl;
48 if (nev < evNumber1) continue;
49 if (nparticles <= 0) return;
51 Int_t nbgr_ev=Int_t(nev/nsignal);
52 // printf("nbgr_ev %d\n",nbgr_ev);
53 //if (MUON) MUON->Digitise(nev,nbgr_ev,"Add"," ","galice_bgr.root");
54 if (MUON) MUON->Digitise(nev,nbgr_ev,"rien"," ","galice_bgr.root");
56 // sprintf(hname,"TreeD%d",nev);