]>
Commit | Line | Data |
---|---|---|
1 | #ifndef __CINT__ | |
2 | #include "alles.h" | |
3 | #include "AliRun.h" | |
4 | #include "AliRunLoader.h" | |
5 | #include "AliLoader.h" | |
6 | #include "AliMagF.h" | |
7 | #include "AliTPCtracker.h" | |
8 | #include "AliITS.h" | |
9 | #include "AliITSgeom.h" | |
10 | #include "AliITSRecPoint.h" | |
11 | #include "AliITSclusterV2.h" | |
12 | #include "AliITSsimulationFastPoints.h" | |
13 | #include "AliITStrackerV2.h" | |
14 | ||
15 | #endif | |
16 | Int_t AliTPCHits2SDigits(Int_t nevent=1) | |
17 | { | |
18 | ||
19 | // new version by J.Belikov | |
20 | ||
21 | // Connect the Root Galice file containing Geometry, Kine and Hits | |
22 | ||
23 | //it assures full cleaning of prevous session | |
24 | if (gAlice) | |
25 | { | |
26 | delete gAlice->GetRunLoader(); | |
27 | delete gAlice;//if everything was OK here it is already NULL | |
28 | gAlice = 0x0; | |
29 | } | |
30 | ||
31 | AliRunLoader *rl = AliRunLoader::Open("galice.root","Event","update"); | |
32 | if (!rl) | |
33 | { | |
34 | cerr<<"Can't load RunLoader from "<<inFile_new<<" !\n"; | |
35 | return 1; | |
36 | } | |
37 | ||
38 | // Get AliRun object from file or create it if not on file | |
39 | ||
40 | rl->LoadgAlice(); | |
41 | ||
42 | gAlice = rl->GetAliRun(); | |
43 | if (!gAlice) { | |
44 | cerr<<"AliTPCHits2Digits.C : AliRun object not found on file\n"; | |
45 | delete rl; | |
46 | return 2; | |
47 | } | |
48 | ||
49 | // gAlice->GetEvent(0); | |
50 | AliTPC *TPC = (AliTPC*)gAlice->GetDetector("TPC"); | |
51 | AliLoader * tpcl = rl->GetLoader("TPCLoader"); | |
52 | if ((TPC == 0x0) || (tpcl == 0x0)) | |
53 | { | |
54 | cerr<<"AliTPCHits2Digits.C : Can not find TPC or TPCLoader\n"; | |
55 | delete rl; | |
56 | return 3; | |
57 | } | |
58 | ||
59 | tpcl->LoadHits("READ"); | |
60 | tpcl->LoadSDigits("RECREATE"); | |
61 | ||
62 | TStopwatch timer; | |
63 | timer.Start(); | |
64 | ||
65 | // uncomment lines below to set active sectors | |
66 | //Int_t sec[10]={4,5,6,7,8,4+36,5+36,6+36,7+36,8+36}; | |
67 | //TPC->SetActiveSectors(sec,10); | |
68 | ||
69 | for(Int_t eventn =0;eventn<nevent;eventn++){ | |
70 | printf("Processing event %d \n",eventn); | |
71 | rl->GetEvent(eventn); | |
72 | TPC->SetTreeAddress(); | |
73 | TPC->SetActiveSectors(); // all sectors set active | |
74 | printf("\nActive sectors\n"); | |
75 | for (Int_t i=0;i<72;i++) if (TPC->IsSectorActive(i)) printf("%d\t",i); | |
76 | ||
77 | TPC->Hits2SDigits2(eventn); | |
78 | } | |
79 | ||
80 | delete rl; | |
81 | timer.Stop(); | |
82 | timer.Print(); | |
83 | ||
84 | return 0; | |
85 | }; | |
86 | ||
87 | ||
88 | ||
89 |