]> git.uio.no Git - u/mrichter/AliRoot.git/blame - START/AliSTARTRawReader.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Library split
[u/mrichter/AliRoot.git] / START / AliSTARTRawReader.cxx
CommitLineData
ee74718b 1#include "AliSTARTRawReader.h"
2#include "AliSTARTRawData.h"
db173afc 3#include "AliSTARTdigit.h"
ee74718b 4
5#include <Riostream.h>
6#include "TMath.h"
7#include "TH1F.h"
8#include "TArrayI.h"
9#include "AliLog.h"
10
11ClassImp(AliSTARTRawReader)
db173afc 12
13 AliSTARTRawReader::AliSTARTRawReader (AliRawReader *rawReader, TTree* tree)
14 : TTask("STARTRawReader","read raw START data"),
15 fTree(tree),
16 fRawReader(rawReader)
ee74718b 17{
18 //
19// create an object to read STARTraw digits
20 AliDebug(1,"Start ");
db173afc 21
ee74718b 22
23}
24 AliSTARTRawReader::~AliSTARTRawReader ()
25{
db173afc 26 //
ee74718b 27}
28//------------------------------------------------------------------------------------------------
29
30UInt_t AliSTARTRawReader::UnpackWord(UInt_t packedWord, Int_t startBit, Int_t stopBit)
31{
32 //This method unpacks a words of StopBit-StartBit+1 bits starting from "StopBits"
33 UInt_t word;
34 UInt_t offSet;
35 Int_t length;
36 length=stopBit-startBit+1;
37 offSet=(UInt_t)TMath::Power(2,length)-1;
38 offSet<<=startBit;
39 word=packedWord&offSet;
40 word>>=startBit;
41 return word;
42}
43//---------------------------------------------------------------------------------------
db173afc 44void AliSTARTRawReader::NextThing()
ee74718b 45{
46// read the next raw digit
47// returns kFALSE if there is no digit left
ee74718b 48 UInt_t word, unpackword;
db173afc 49 Int_t time, adc, pmt;
50 TArrayI *timeTDC1 = new TArrayI(24);
51 TArrayI * chargeTDC1 = new TArrayI(24);
52 TArrayI *timeTDC2 = new TArrayI(24);
53 TArrayI *chargeTDC2 = new TArrayI(24);
54 TArrayI *sumMult = new TArrayI(6);
55
a93ccc90 56 fRawReader->Reset();
ee74718b 57 fRawReader->Select(13);
58
db173afc 59 if (!fRawReader->ReadHeader()){
60 Error("ReadSTARTRaw","Couldn't read header");
61 return;
62 }
63 AliSTARTdigit *fDigits = new AliSTARTdigit();
64 fTree->Branch("START","AliSTARTdigit",&fDigits,400,1);
ee74718b 65
db173afc 66 fRawReader->ReadNextInt(word);
67 for (Int_t i=0; i<24; i++)
ee74718b 68 {
69 word=0;
70 unpackword=0;
71 fRawReader->ReadNextInt(word);
72 unpackword=UnpackWord(word,0,5);
db173afc 73 pmt=unpackword;
ee74718b 74 word=0;
75 unpackword=0;
ee74718b 76 fRawReader->ReadNextInt(word);
77 unpackword=UnpackWord(word,8,31);
db173afc 78 time=unpackword;
79 timeTDC1->AddAt(time,pmt);
80
ee74718b 81 word=0;
82 unpackword=0;
ee74718b 83 fRawReader->ReadNextInt(word);
84 unpackword=UnpackWord(word,0,5);
db173afc 85 pmt=unpackword;
ee74718b 86 word=0;
87 unpackword=0;
ee74718b 88 fRawReader->ReadNextInt(word);
ee74718b 89 unpackword=UnpackWord(word,8,31);
db173afc 90 time=unpackword;
91 timeTDC2->AddAt(time,pmt);
92
93 // QTC
ee74718b 94 word=0;
95 unpackword=0;
96 fRawReader->ReadNextInt(word);
97 unpackword=UnpackWord(word,0,5);
db173afc 98 pmt=unpackword;
ee74718b 99 word=0;
100 unpackword=0;
ee74718b 101 fRawReader->ReadNextInt(word);
102 unpackword= UnpackWord(word,8,31);
db173afc 103 time=unpackword; //T1
104
ee74718b 105 word=0;
106 unpackword=0;
107
108 fRawReader->ReadNextInt(word);
109 unpackword=UnpackWord(word,0,5);
db173afc 110 pmt=unpackword;
111 word=0;
112 unpackword=0;
113 fRawReader->ReadNextInt(word);
114 unpackword= UnpackWord(word,8,31);
115 adc=unpackword-time; // T1+ T2 (A)
116 chargeTDC1->AddAt(adc,pmt);
117 //QTC amplified
ee74718b 118
119 word=0;
120 unpackword=0;
121 fRawReader->ReadNextInt(word);
db173afc 122 unpackword=UnpackWord(word,0,5);
123 pmt=unpackword;
124 word=0;
125 unpackword=0;
126 fRawReader->ReadNextInt(word);
ee74718b 127 unpackword= UnpackWord(word,8,31);
db173afc 128 time=unpackword; //T1
129
130 word=0;
131 unpackword=0;
ee74718b 132
db173afc 133 fRawReader->ReadNextInt(word);
134 unpackword=UnpackWord(word,0,5);
135 pmt=unpackword;
136 word=0;
137 unpackword=0;
138 fRawReader->ReadNextInt(word);
139 unpackword= UnpackWord(word,8,31);
140 adc=unpackword-time; // T1+ T2 (A)
141 chargeTDC2->AddAt(adc,pmt);
ee74718b 142
db173afc 143
ee74718b 144 }
db173afc 145 fDigits->SetTime(*timeTDC1);
146 fDigits->SetADC(*chargeTDC1);
147
148
149 word=0;
150 unpackword=0;
151
152 fRawReader->ReadNextInt(word);
153 unpackword=UnpackWord(word,0,5);
154 pmt=unpackword;
155
156 word=0;
157 unpackword=0;
158
159 fRawReader->ReadNextInt(word);
160 unpackword=UnpackWord(word,8,31);
161 time=unpackword;
162 fDigits->SetMeanTime(time);
163
164 // Best time right &left
165 word=0;
166 unpackword=0;
167
168 fRawReader->ReadNextInt(word);
169 unpackword=UnpackWord(word,0,5);
170 pmt=unpackword;
171
172 word=0;
173 unpackword=0;
174
175 fRawReader->ReadNextInt(word);
176 unpackword=UnpackWord(word,8,31);
177 time=unpackword;
178 fDigits->SetTimeBestRight(time);
179
180
181 // best time left
182 word=0;
183 unpackword=0;
184
185 fRawReader->ReadNextInt(word);
186 unpackword=UnpackWord(word,0,5);
187 pmt=unpackword;
188
189 word=0;
190 unpackword=0;
191
192 fRawReader->ReadNextInt(word);
193 unpackword=UnpackWord(word,8,31);
194 time=unpackword;
195 fDigits->SetTimeBestLeft(time);
196
197
198 // best time differece
199 word=0;
200 unpackword=0;
201
202 fRawReader->ReadNextInt(word);
203 unpackword=UnpackWord(word,0,5);
204 pmt=unpackword;
205
206 word=0;
207 unpackword=0;
208
209 fRawReader->ReadNextInt(word);
210 unpackword=UnpackWord(word,8,31);
211 time=unpackword;
212 fDigits->SetDiffTime(time);
213
214 // multiplicity
215 for (Int_t im=0; im<6; im++)
216 {
217 word=0;
218 unpackword=0;
219 fRawReader->ReadNextInt(word);
220 unpackword=UnpackWord(word,0,5);
221 pmt=unpackword;
222 word=0;
223 unpackword=0;
224 fRawReader->ReadNextInt(word);
225 unpackword=UnpackWord(word,8,31);
226 time=unpackword;
227 sumMult->AddAt(time,im);
228 }
229 fDigits->SetSumMult(*sumMult);
230
231 fTree->Fill();
232
ee74718b 233}
db173afc 234