]>
Commit | Line | Data |
---|---|---|
0df3ca52 | 1 | #include "iostream.h" |
2 | ||
3 | void MUONraw (Int_t evNumber1=0,Int_t evNumber2=0) | |
4 | { | |
5 | ////////////////////////////////////// | |
6 | // // | |
7 | // ROOT macro for ALICE Dimuon Arm: // | |
8 | // Clusterization of digits // | |
9 | // // | |
10 | ////////////////////////////////////// | |
11 | // | |
12 | // Adds the tree TR for raw clusters | |
13 | // to the ROOT file "galice.root" | |
14 | // containing the digits (tree TD). | |
15 | // | |
16 | // Arguments: | |
17 | // evNumber1 = first event number to act on in file "galice.root" | |
18 | // evNumber2 = last event number to act on in file "galice.root" | |
19 | // | |
20 | // Input/output file: | |
21 | // "galice.root" | |
22 | // | |
23 | //__________________________________________________________________________ | |
24 | ||
25 | // Dynamically link some shared libs | |
26 | ||
27 | if (gClassTable->GetID("AliRun") < 0) { | |
28 | gROOT->LoadMacro("loadlibs.C"); | |
29 | loadlibs(); | |
30 | } | |
31 | ||
32 | // Connect the Root Galice file containing Geometry, Kine and Hits | |
33 | ||
34 | TFile *file = (TFile*)gROOT->GetListOfFiles()->FindObject("galice.root"); | |
35 | if (!file) file = new TFile("galice.root","UPDATE"); | |
36 | ||
37 | // Get AliRun object from file or create it if not on file | |
38 | ||
39 | if (!gAlice) { | |
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"); | |
43 | } | |
44 | ||
45 | // Set reconstruction models | |
46 | // | |
47 | // Get pointers to Alice detectors and Digits containers | |
48 | AliMUON *MUON = (AliMUON*) gAlice->GetModule("MUON"); | |
49 | RecModel = new AliMUONClusterFinderAZ(0,1); //AZ | |
50 | //RecModel = new AliMUONClusterFinderAZ(1,0); //AZ | |
51 | for (Int_t i=0; i<10; i++) { | |
52 | MUON->SetReconstructionModel(i,(AliMUONClusterFinderVS*)RecModel); | |
53 | } | |
54 | // | |
55 | // Loop over events | |
56 | // | |
57 | Int_t Nh=0; | |
58 | Int_t Nh1=0; | |
59 | gAlice->RunReco("MUON", evNumber1, evNumber2); | |
60 | } | |
61 |