3 void MUONrawclusters (char* filename, 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("$ALICE_ROOT/macros/loadlibs.C");
32 // Creating Run Loader and openning file containing Hits
33 AliRunLoader * RunLoader = AliRunLoader::Open(filename,"MUONFolder","UPDATE");
34 if (RunLoader ==0x0) {
35 printf(">>> Error : Error Opening %s file \n",filename);
39 // Loading AliRun master
40 RunLoader->UnloadgAlice();
41 RunLoader->LoadgAlice();
42 gAlice = RunLoader->GetAliRun();
44 // Loading MUON subsystem
45 AliMUON * MUON = (AliMUON *) gAlice->GetDetector("MUON");
46 AliLoader * MUONLoader = RunLoader->GetLoader("MUONLoader");
48 Int_t ievent, nevents;
49 nevents = RunLoader->GetNumberOfEvents();
51 for (Int_t i=0; i<10; i++) {
52 RecModel = new AliMUONClusterFinderVS();
53 // RecModel->SetTracks(16,17);
54 // RecModel->SetTracks(266,267);
55 RecModel->SetGhostChi2Cut(10);
56 MUON->SetReconstructionModel(i,RecModel);
63 // gAlice->RunReco("MUON", evNumber1, evNumber2);
64 for(Int_t ievent=evNumber1; ievent<evNumber2; ievent++) {
65 printf("event %d\n",ievent);
66 RunLoader->GetEvent(ievent);
67 MUONLoader->LoadDigits("read");
68 if (MUONLoader->TreeR() == 0x0) MUONLoader->MakeTree("R");
69 MUON->MakeBranch("R");
70 MUON->SetTreeAddress();
72 MUONLoader->UnloadDigits();
73 MUONLoader->UnloadRecPoints();