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 // Alice segment manager object //
22 ///////////////////////////////////////////////////////////////////////////////
23 #include <Riostream.h>
29 #include "TClonesArray.h"
30 #include "TDirectory.h"
36 #include "AliSegmentID.h"
37 #include "AliSegmentArray.h"
38 #include "TObjString.h"
41 //_____________________________________________________________________________
42 ClassImp(AliSegmentArray)
44 AliSegmentArray::AliSegmentArray()
58 AliSegmentArray::AliSegmentArray(const char *classname, Int_t n)
63 // Create an array of objects of classname. The class must inherit from
64 // AliSegmentID . The second argument adjust number of entries in
72 SetName("SegmentArray");
73 SetTitle("SegmentArray");
76 if (MakeArray(n)==kFALSE){
77 Error("AliSegmentArray", "can't allocate %d segments in memory",n);
82 AliSegmentArray::AliSegmentArray(const AliSegmentArray &segment):
87 // to be later implemented
90 AliSegmentArray &AliSegmentArray::operator = (const AliSegmentArray & /*segment*/)
93 //to be later implemented
97 AliSegmentArray::~AliSegmentArray()
100 // default destructor
106 if (fTreeOwner) delete fTree;
109 if (fTreeIndex) delete fTreeIndex;
110 // if (fClass!=0) delete fClass;
114 Bool_t AliSegmentArray::SetClass(const char *classname)
117 //set class of stored object
119 //delete fClass; not ower of fClass
123 if (fTreeOwner) delete fTree;
138 ::Fatal("AliSegmentArray::AliSegmentArray", "ROOT system not initialized");
140 fClass = gROOT->GetClass(classname);
142 Error("AliSegmentArray", "%s is not a valid class name", classname);
145 if (!fClass->InheritsFrom(AliSegmentID::Class())) {
146 Error("AliSegmentArray", "%s does not inherit from AliSegmentID", classname);
153 AliSegmentID * AliSegmentArray::NewSegment()
156 //create object according class information
157 if (fClass==0) return 0;
158 AliSegmentID * segment = (AliSegmentID * )fClass->New();
159 if (segment == 0) return 0;
164 Bool_t AliSegmentArray::AddSegment(AliSegmentID *segment)
167 // add segment to array
169 if (segment==0) return kFALSE;
170 if (fSegment==0) return kFALSE;
171 if (fClass==0) return kFALSE;
172 if (!(segment->IsA()->InheritsFrom(fClass))){
173 Error("AliSegmentArray", "added class %s is not of proper type ",
174 segment->IsA()->GetName());
177 fSegment->AddAt(segment,segment->GetID());
178 fNSegment = fSegment->GetLast()+1;
182 AliSegmentID * AliSegmentArray::AddSegment(Int_t index)
185 // add segment to array
187 if (fSegment==0) return 0;
188 if (fClass==0) return 0;
189 AliSegmentID * segment = NewSegment();
190 if (segment == 0) return 0;
191 fSegment->AddAt(segment,index);
192 segment->SetID(index);
193 fNSegment = fSegment->GetLast()+1;
198 void AliSegmentArray::ClearSegment(Int_t index)
201 //remove segment from active memory
203 //PH if ((*fSegment)[index]){
204 if (fSegment->At(index)){
205 // (*fSegment)[index]->Delete(); //not working for TClonesArray
206 //PH delete (*fSegment)[index]; //because problem with deleting TClonesArray
207 //PH fSegment->RemoveAt(index);
208 delete fSegment->RemoveAt(index);
213 Bool_t AliSegmentArray::MakeArray(Int_t n)
216 //make array of pointers to Segments
222 fSegment = new TObjArray(n);
224 if (fSegment) return kTRUE;
227 void AliSegmentArray::MakeTree(TTree* tree)
229 //Make tree with the name
230 AliSegmentID * psegment = NewSegment();
232 //PH fBranch = fTree->Branch("Segment",psegment->IsA()->GetName(),&psegment,64000);
233 fBranch = fTree->Branch("Segment",psegment->IsA()->GetName(),&psegment,64000,99);
237 void AliSegmentArray::MakeTree(char *file)
239 // AliSegmentID segment;
240 AliSegmentID * psegment = NewSegment();
245 fTree = new TTree("Segment Tree","Tree with segments");
250 //PH fBranch = fTree->Branch("Segment",psegment->IsA()->GetName(),&psegment,64000);
251 fBranch = fTree->Branch("Segment",psegment->IsA()->GetName(),&psegment,64000,99);
254 TString outFile = gAlice->GetBaseFile();
255 outFile = outFile + "/" + file;
256 fBranch->SetFile(outFile.Data());
257 TDirectory *wd = gDirectory;
258 TBranch *b = fBranch;
259 TIter next( b->GetListOfBranches());
260 while ((b=(TBranch*)next())) {
261 b->SetFile(outFile.Data());
263 cout << "Diverting branch " << "Segment" << " to file " << outFile << endl;
270 Bool_t AliSegmentArray::MakeDictionary(Int_t size)
273 //create index table for tree
275 if (size<1) return kFALSE;
276 if (fTreeIndex) delete fTreeIndex;
277 fTreeIndex = new TArrayI();
278 fTreeIndex->Set(size);
280 AliSegmentID * psegment = NewSegment(); //MI change
281 fBranch->SetAddress(&psegment);
282 TBranch * brindix = fTree->GetBranch("fSegmentID");
283 Int_t nevent = (Int_t)fTree->GetEntries();
284 for (Int_t i = 0; i<nevent; i++){
285 brindix->GetEvent(i);
286 Int_t treeIndex=psegment->GetID();
287 if (fTreeIndex->fN<treeIndex) fTreeIndex->Set(Int_t(Float_t(treeIndex)*1.5)+1);
288 // Int_t index = segment.GetID();
289 (*fTreeIndex)[treeIndex]=i+1; //
291 if (psegment) delete psegment;
295 Bool_t AliSegmentArray::ConnectTree(TTree* tree)
298 if (fTree == 0) return kFALSE;
299 fBranch = fTree->GetBranch("Segment");
300 if (fBranch==0) return kFALSE;
301 MakeDictionary(TMath::Max(fNSegment,Int_t(fTree->GetEntries())));
302 MakeArray(fTreeIndex->fN);
307 Bool_t AliSegmentArray::ConnectTree(const char * treeName)
309 //connect tree from current directory
319 fTree =(TTree*)gDirectory->Get(treeName);
321 if (fTree == 0) return kFALSE;
322 fBranch = fTree->GetBranch("Segment");
323 if (fBranch==0) return kFALSE;
324 MakeDictionary(TMath::Max(fNSegment,Int_t(fTree->GetEntries())));
325 MakeArray(fTreeIndex->fN);
329 AliSegmentID *AliSegmentArray::LoadSegment(Int_t index)
332 //load segment with index to the memory
335 if (fTreeIndex ==0 ) MakeDictionary(3000);
336 //firstly try to load dictionary
337 if (fTreeIndex ==0 ) return 0;
338 if (fBranch==0) return 0;
339 if (index>fTreeIndex->fN) return 0;
340 //PH AliSegmentID *s = (AliSegmentID*)(*fSegment)[index];
341 AliSegmentID *s = (AliSegmentID*)fSegment->At(index);
342 if (s==0) s= NewSegment();
344 // new AliSegmentID(index);
347 Int_t treeIndex =(*fTreeIndex)[index];
348 if (treeIndex<1) return 0;
349 else treeIndex--; //I don't like it Int table I have index shifted by 1
350 fBranch->SetAddress(&s);
351 fTree->GetEvent(treeIndex);
352 //PH (*fSegment)[index] = (TObject*) s;
353 fSegment->AddAt((TObject*) s, index);
360 AliSegmentID *AliSegmentArray::LoadEntry(Int_t index)
363 //load segment at position inex in tree to the memory
366 if (fBranch==0) return 0;
367 if (index>fTree->GetEntries()) return 0;
368 AliSegmentID * s = NewSegment();
371 fBranch->SetAddress(&s);
372 fTree->GetEvent(index);
376 Int_t nindex = s->GetID();
377 ClearSegment(nindex);
378 //PH (*fSegment)[nindex] = (TObject*) s;
379 fSegment->AddAt((TObject*) s, nindex);
383 void AliSegmentArray::StoreSegment(Int_t index)
386 //make segment persistent
388 const AliSegmentID * ksegment = (*this)[index];
389 if (ksegment == 0 ) return;
390 if (fTree==0) MakeTree();
391 fBranch->SetAddress(&ksegment);
396 void AliSegmentArray::Streamer(TBuffer &R__b)
398 TObjString treeName, * ptreeName=&treeName;
399 if (R__b.IsReading()) {
400 Version_t R__v = R__b.ReadVersion(); if (R__v) { }
401 TNamed::Streamer(R__b);
403 if (fTree && fTreeOwner) delete fTree;
404 ConnectTree(ptreeName->String());
406 R__b.WriteVersion(AliSegmentArray::IsA());
407 TNamed::Streamer(R__b);
409 // sprintf(ch,"%s",fTrre->GetTitle());
410 treeName.String() = fTree->GetTitle();