]>
Commit | Line | Data |
---|---|---|
d19e1df9 | 1 | Int_t FindKrClustersRaw(const char *fileName="data.root"){ |
2 | ||
3 | //define tree | |
4 | TFile *hfile=new TFile("adc.root","RECREATE","ADC file"); | |
5 | // Create a ROOT Tree | |
6 | TTree *mytree = new TTree("Kr","Krypton cluster tree"); | |
7 | ||
8 | ||
9 | AliRawReader *reader = new AliRawReaderRoot(fileName); | |
10 | //AliRawReader *reader = new AliRawReaderDate(fileName); | |
11 | reader->Reset(); | |
12 | ||
13 | TStopwatch timer; | |
14 | timer.Start(); | |
15 | ||
16 | AliAltroRawStreamFast* stream = new AliAltroRawStreamFast(reader); | |
17 | stream->SelectRawData("TPC"); | |
18 | ||
19 | //one general output | |
20 | AliTPCclustererKr *clusters = new AliTPCclustererKr(); | |
21 | clusters->SetOutput(mytree); | |
22 | clusters->SetRecoParam(0); | |
23 | ||
24 | //only for geometry parameters loading - temporarly | |
25 | AliRunLoader* rl = AliRunLoader::Open("galice.root"); | |
26 | AliTPCParam *param=(AliTPCParamSR *)gDirectory->Get("75x40_100x60_150x60"); | |
27 | //if (!param) {cerr<<"TPC parameters have not been found !\n"; return 4;} | |
28 | clusters->SetParam(param); | |
29 | ||
30 | ||
31 | Int_t evtnr=0; | |
32 | while (reader->NextEvent()) { | |
33 | //output for each event | |
34 | // AliTPCclustererKr *clusters = new AliTPCclustererKr(); | |
35 | // clusters->SetOutput(mytree); | |
36 | ||
37 | //if(evtnr++<35)continue; | |
38 | cout<<"Evt = "<<evtnr<<endl; | |
39 | clusters->finderIO(reader); | |
40 | evtnr++; | |
41 | ||
42 | ||
43 | } | |
44 | ||
45 | //mytree->Print();//print rootuple summary | |
46 | // Save all objects in this file | |
47 | hfile->Write(); | |
48 | // Close the file | |
49 | hfile->Close(); | |
50 | ||
51 | timer.Stop(); | |
52 | timer.Print(); | |
53 | ||
54 | delete stream; | |
55 | ||
56 | return 0; | |
57 | } |