3 void MUONrawclusters (Int_t evNumber1=0,Int_t evNumber2=0)
5 //////////////////////////////////////
7 // ROOT macro for ALICE Dimuon Arm: //
8 // Clusterization of digits //
10 //////////////////////////////////////
12 // Adds the tree TR for raw clusters
13 // to the ROOT file "galice.root"
14 // containing the digits (tree TD).
17 // evNumber1 = first event number to act on in file "galice.root"
18 // evNumber2 = last event number to act on in file "galice.root"
23 //__________________________________________________________________________
25 // Dynamically link some shared libs
27 if (gClassTable->GetID("AliRun") < 0) {
28 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 = new TFile("galice.root","UPDATE");
37 // Get AliRun object from file or create it if not on file
40 gAlice = (AliRun*)file->Get("gAlice");
41 if (gAlice) printf("AliRun object found on file\n");
42 if (!gAlice) gAlice = new AliRun("gAlice","Alice test program");
45 // Set reconstruction models
47 // Get pointers to Alice detectors and Digits containers
48 AliMUON *MUON = (AliMUON*) gAlice->GetModule("MUON");
49 for (Int_t i=0; i<10; i++) {
50 RecModel = new AliMUONClusterFinderVS();
51 // RecModel->SetTracks(16,17);
52 // RecModel->SetTracks(266,267);
53 MUON->SetReconstructionModel(i,RecModel);
60 for (int nev=evNumber1; nev<= evNumber2; nev++) {
61 Int_t nparticles = gAlice->GetEvent(nev);
62 cout << "nev " << nev <<endl;
63 cout << "nparticles " << nparticles <<endl;
64 if (nev < evNumber1) continue;
65 if (nparticles <= 0) return;
67 TClonesArray *Particles = gAlice->Particles();
68 TTree *TD = gAlice->TreeD();
69 Int_t nent=TD->GetEntries();
71 MUON->FindClusters(nev,nent-2);