]>
Commit | Line | Data |
---|---|---|
2e9f335b | 1 | #if !defined(__CINT__) || defined(__MAKECINT__) |
2 | #include "TFile.h" | |
3 | #include "AliTPCBuffer.h" | |
4 | #endif | |
5 | ||
6 | ||
7 | void AliTPCDDL(char* FileName,Int_t eth=0){ | |
8 | //eth is a threshold. | |
a79660fb | 9 | //Digits stored into a file have an amplitude value greater than "eth" |
2e9f335b | 10 | TFile *cf=TFile::Open(FileName); |
11 | // old geometry (3.07) | |
12 | //AliTPCParamSR *param =(AliTPCParamSR *)cf->Get("75x40_100x60"); | |
13 | // if new geometry comment out the line above and uncomment the one below | |
14 | AliTPCParamSR *param =(AliTPCParamSR *)cf->Get("75x40_100x60_150x60"); | |
15 | AliTPCDigitsArray *digarr=new AliTPCDigitsArray; | |
16 | digarr->Setup(param); | |
17 | char cname[100]; | |
18 | //old geometry | |
19 | //sprintf(cname,"TreeD_75x40_100x60_%d",eventn); | |
20 | // if new geometry comment out the line above and uncomment the one below | |
21 | Int_t eventn=0; | |
22 | sprintf(cname,"TreeD_75x40_100x60_150x60_%d",eventn); | |
23 | digarr->ConnectTree(cname); | |
24 | AliTPCBuffer *b=new AliTPCBuffer("AliTPCDDL.dat"); | |
9f992f70 | 25 | |
26 | //Verbose level | |
27 | // 0: Silent | |
28 | // 1: cout messages | |
29 | // 2: txt files with digits | |
30 | //BE CAREFUL, verbose level 2 MUST be used only for debugging and | |
31 | //it is highly suggested to use this mode only for debugging digits files | |
32 | //reasonably small, because otherwise the size of the txt files can reach | |
33 | //quickly several MB wasting time and disk space. | |
34 | b->SetVerbose(0); | |
35 | ||
36 | ||
2e9f335b | 37 | Int_t nrows=Int_t(digarr->GetTree()->GetEntries()); |
38 | cout<<"Number of entries "<<nrows<<endl; | |
39 | Int_t PSector=-1; | |
40 | Int_t SubSec=0; | |
41 | for (Int_t n=0; n<nrows; n++) { | |
42 | AliSimDigits *digrow=(AliSimDigits*)digarr->LoadEntry(n); | |
43 | ||
44 | Int_t sec,row; // sector and row number (in the TPC) | |
45 | param->AdjustSectorRow(digrow->GetID(),sec,row); | |
46 | // cout<<sec<<" row "<<row<<endl; | |
47 | if(PSector!=sec){ | |
48 | SubSec=0; | |
49 | PSector=sec; | |
50 | }//end if | |
51 | ||
52 | if(sec<36){ | |
53 | //inner sector [0;35] | |
54 | if(row!=30) | |
55 | //the whole row is written into the output file | |
56 | b->WriteRowBinary(eth,digrow,0,0,0,sec,SubSec,row); | |
57 | else{ | |
58 | //only the pads in the range [37;48] are written into the output file | |
59 | b->WriteRowBinary(eth,digrow,37,48,1,sec,SubSec,row); | |
60 | SubSec=1; | |
61 | //only the pads outside the range [37;48] are written into the output file | |
62 | b->WriteRowBinary(eth,digrow,37,48,2,sec,SubSec,row); | |
63 | }//end else | |
64 | }//end if | |
65 | else{ | |
66 | //outer sector [36;71] | |
67 | if(row==54)SubSec=2; | |
68 | if((row!=27)&&(row!=76)) | |
69 | b->WriteRowBinary(eth,digrow,0,0,0,sec,SubSec,row); | |
70 | else{ | |
71 | if(row==27){ | |
72 | //only the pads outside the range [43;46] are written into the output file | |
73 | b->WriteRowBinary(eth,digrow,43,46,2,sec,SubSec,row); | |
74 | SubSec=1; | |
75 | //only the pads in the range [43;46] are written into the output file | |
76 | b->WriteRowBinary(eth,digrow,43,46,1,sec,SubSec,row); | |
77 | } | |
78 | if(row==76){ | |
79 | //only the pads outside the range [33;88] are written into the output file | |
80 | b->WriteRowBinary(eth,digrow,33,88,2,sec,SubSec,row); | |
81 | SubSec=3; | |
82 | //only the pads in the range [33;88] are written into the output file | |
83 | b->WriteRowBinary(eth,digrow,33,88,1,sec,SubSec,row); | |
84 | } | |
85 | } | |
86 | }//end else | |
87 | }//end for | |
88 | cf->Close(); | |
89 | cout<<"File created !"<<endl; | |
90 | cout<<"Total number of digits: "<<b->GetDigNumber()<<endl; | |
91 | delete b; | |
92 | return; | |
93 | }//end AliTPCDataChallenge |