1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 ///////////////////////////////////////////////////////////////////////////////
20 // T0 (T-Zero) Detector //
21 // This class contains the base procedures for the T0 //
26 <img src="gif/AliT0Class.gif">
29 <font size=+2 color=red>
30 <p>The responsible person for this module is
31 <a href="mailto:Alla.Maevskaia@cern.ch">Alla Maevskaia</a>.
38 ///////////////////////////////////////////////////////////////////////////////
40 #include "TClonesArray.h"
42 #include "AliLoader.h"
48 #include "AliT0Digitizer.h"
49 #include "AliT0RawData.h"
50 #include "AliT0RecPoint.h"
51 #include "AliT0digit.h"
56 //static AliT0digit *digits;
58 //_____________________________________________________________________________
60 : AliDetector(), fIdSens(0), fDigits(NULL), fRecPoints(NULL)
63 // Default constructor for class AliT0
71 //_____________________________________________________________________________
72 AliT0::AliT0(const char *name, const char *title)
73 : AliDetector(name,title), fIdSens(0), fDigits(new AliT0digit()), fRecPoints(new AliT0RecPoint())
76 // Standard constructor for T0 Detector
81 // Initialise Hit array
82 fHits = new TClonesArray("AliT0hit", 405);
83 gAlice->GetMCApp()->AddHitList(fHits);
84 // fDigits = new AliT0digit();
85 // fRecPoints = new AliT0RecPoint();
88 //PH SetMarkerColor(kRed);
91 //_____________________________________________________________________________
103 cout<<" delete fDigits; "<<endl;
106 fRecPoints ->Delete();
108 cout<<" delete fRecPoints; "<<endl;
113 //_____________________________________________________________________________
114 void AliT0::AddHit(Int_t track, Int_t *vol, Float_t *hits)
119 TClonesArray &lhits = *fHits;
120 new(lhits[fNhits++]) AliT0hit(fIshunt,track,vol,hits);
124 //_____________________________________________________________________________
126 void AliT0::AddDigit(Int_t besttimeright, Int_t besttimeleft, Int_t meantime,
127 Int_t timediff, Int_t sumMult, Int_t refpoint,
128 TArrayI *timeCFD, TArrayI *qt0, TArrayI *timeLED, TArrayI *qt1)
131 // Add a T0 digit to the list.
135 fDigits = new AliT0digit();
137 fDigits-> SetTimeBestA(besttimeright);
138 fDigits->SetTimeBestC(besttimeleft);
139 fDigits-> SetMeanTime(meantime);
140 fDigits-> SetDiffTime(timediff);
141 fDigits-> SetSumMult(sumMult);
142 fDigits->SetTimeCFD(*timeCFD);
143 fDigits->SetTimeLED(*timeLED);
144 fDigits->SetQT0(*qt0);
145 fDigits->SetQT1(*qt1);
146 fDigits->SetRefPoint(refpoint);
149 //-------------------------------------------------------------------------
153 // Initialis the T0 after it has been built
156 if(AliLog::GetGlobalDebugLevel()>0) {
157 printf("\n%s: ",ClassName());
158 for(i=0;i<35;i++) printf("*");
160 for(i=0;i<35;i++) printf("*");
161 printf("\n%s: ",ClassName());
163 // Here the T0 initialisation code (if any!)
164 for(i=0;i<80;i++) printf("*");
169 //---------------------------------------------------------------------------
170 void AliT0::MakeBranch(Option_t* option)
173 // Create Tree branches for the T0.
177 // H Make a branch of TClonesArray of AliT0Hit's
178 // D Make a branch of TClonesArray of AliT0Digit's
180 // R Make a branch of AliT0RecPoints
183 sprintf(branchname,"%s",GetName());
185 const char *cH = strstr(option,"H");
186 const char *cD = strstr(option,"D");
187 const char *cR = strstr(option,"R");
188 const char *cS = strstr(option,"S");
190 if (cH && fLoader->TreeH())
192 if (fHits == 0x0) fHits = new TClonesArray("AliT0hit", 405);
193 AliDetector::MakeBranch(option);
197 if (cD && fLoader->TreeD())
199 if (fDigits == 0x0) fDigits = new AliT0digit();
200 // MakeBranchInTree(fLoader->TreeD(), branchname,
201 // &fDigits, 405, 0);
202 fLoader->TreeD()->Branch(branchname,"AliT0digit",&fDigits,405,1);
203 // fLoader->TreeD()->Print();
205 if (cR && fLoader->TreeR())
207 if (fRecPoints == 0x0) fRecPoints = new AliT0RecPoint();
208 MakeBranchInTree(fLoader->TreeR(), branchname,
209 &fRecPoints, 405, 0);
211 if (cS && fLoader->TreeS())
213 if (fDigits == 0x0) fDigits = new AliT0digit();
214 // MakeBranchInTree(fLoader->TreeD(), branchname,
215 // &fDigits, 405, 0);
216 fLoader->TreeS()->Branch(branchname,"AliT0digit",&fDigits,405,1);
217 // fLoader->TreeD()->Print();
222 //_____________________________________________________________________________
223 void AliT0::ResetHits()
228 AliDetector::ResetHits();
231 //____________________________________________________________________
232 void AliT0::ResetDigits()
235 // Reset number of digits and the digits array for this detector
237 if (fDigits) fDigits->Clear();
240 //_____________________________________________________________________________
241 void AliT0::SetTreeAddress()
244 TTree *treeH = fLoader->TreeH();
248 if (fHits == 0x0) fHits = new TClonesArray("AliT0hit", 405);
251 AliDetector::SetTreeAddress();
252 TTree *treeD = fLoader->TreeD();
254 if (fDigits == 0x0) fDigits = new AliT0digit();
255 TBranch* branch = treeD->GetBranch ("T0");
256 if (branch) branch->SetAddress(&fDigits);
259 TTree *treeR = fLoader->TreeR();
261 if (fRecPoints == 0x0) fRecPoints = new AliT0RecPoint() ;
262 TBranch* branch = treeR->GetBranch ("T0");
263 if (branch) branch->SetAddress(&fRecPoints);
265 // SDigitizer for Federico
266 TTree *treeS = fLoader->TreeS();
268 if (fDigits == 0x0) fDigits = new AliT0digit();
269 TBranch* branch = treeS->GetBranch ("T0");
270 if (branch) branch->SetAddress(&fDigits);
276 //_____________________________________________________________________________
277 void AliT0::MakeBranchInTreeD(TTree *treeD, const char *file)
280 // Create TreeD branches for the FMD
282 const Int_t kBufferSize = 4000;
284 sprintf(branchname,"%s",GetName());
287 MakeBranchInTree(treeD, branchname,&fDigits, kBufferSize, file);
291 //_____________________________________________________________________________
292 AliDigitizer* AliT0::CreateDigitizer(AliRunDigitizer* manager) const
294 return new AliT0Digitizer(manager);
296 //____________________________________________________________________________
297 void AliT0::Digits2Raw()
300 // Starting from the T0 digits, writes the Raw Data objects
302 // AliT0Loader* pStartLoader = (AliT0Loader*)fLoader;
303 fLoader ->LoadDigits("read");
304 TTree* treeD = fLoader->TreeD();
306 AliError("no digits tree");
309 if (fDigits == 0x0) fDigits = new AliT0digit();
311 TBranch *branch = treeD->GetBranch("T0");
313 branch->SetAddress(&fDigits);
315 AliError("Branch T0 DIGIT not found");
318 AliT0RawData rawWriter;
319 rawWriter.SetVerbose(0);
321 AliDebug(2,Form(" Formatting raw data for T0 "));
323 // rawWriter.RawDataT0(treeD->GetBranch("T0"));
324 rawWriter.RawDataT0(fDigits);
327 fLoader->UnloadDigits();
331 //____________________________________________________________________________
332 void AliT0::Raw2Digits(AliRawReader *rawReader,TTree* digitsTree)
335 //T0 raw data-> digits conversion
336 // reconstruct time information from raw data
337 // cout<<" AliT0::Raw2Digits(AliRawReader *rawReader,TTree* digitsTree) "<<
338 // rawReader<<" "<<digitsTree<<endl;
341 // AliT0RawReader myrawreader(rawReader,digitsTree);
342 AliT0RawReader myrawreader(rawReader);
343 if (!myrawreader.Next())
344 AliDebug(1,Form(" no raw data found!! %i", myrawreader.Next()));
345 Int_t allData[110][5];
346 for (Int_t i=0; i<110; i++) {
347 allData[i][0]=myrawreader.GetData(i,0);
351 fDigits = new AliT0digit();
352 digitsTree->Branch("T0","AliT0digit",&fDigits,405,1);
355 TArrayI *timeLED = new TArrayI(24);
356 TArrayI * timeCFD = new TArrayI(24);
357 TArrayI *chargeQT0 = new TArrayI(24);
358 TArrayI *chargeQT1 = new TArrayI(24);
360 for (Int_t in=0; in<24; in++)
362 timeLED->AddAt(allData[in+1][0],in);
363 timeCFD->AddAt(allData[in+25][0],in);
364 chargeQT0->AddAt(allData[in+55][0],in);
365 chargeQT1->AddAt(allData[in+79][0],in);
366 AliDebug(2, Form(" readed Raw %i %i %i %i %i", in, timeLED->At(in),timeCFD->At(in),chargeQT0->At(in),chargeQT1->At(in)));
369 fDigits->SetTimeCFD(*timeCFD);
370 fDigits->SetQT0(*chargeQT1);
372 fDigits->SetTimeLED(*timeLED);
373 fDigits->SetQT1(*chargeQT1);
375 fDigits->SetMeanTime(allData[49][0]);
376 fDigits->SetDiffTime(allData[50][0]);
377 fDigits->SetTimeBestA(allData[51][0]);
378 fDigits->SetTimeBestC(allData[52][0]);