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()
60 AliSegmentArray::AliSegmentArray(const char *classname, Int_t n)
61 :TNamed("SegmentArray","SegmentArray"),
73 // Create an array of objects of classname. The class must inherit from
74 // AliSegmentID . The second argument adjust number of entries in
79 if (MakeArray(n)==kFALSE){
80 Error("AliSegmentArray", "can't allocate %d segments in memory",n);
85 AliSegmentArray::AliSegmentArray(const AliSegmentArray &segment)
98 // to be later implemented
101 AliSegmentArray &AliSegmentArray::operator = (const AliSegmentArray & /*segment*/)
103 //assignment operator
104 //to be later implemented
108 AliSegmentArray::~AliSegmentArray()
111 // default destructor
117 if (fTreeOwner) delete fTree;
118 else fTree->Reset();}
120 if (fTreeIndex) delete fTreeIndex;
121 // if (fClass!=0) delete fClass;
125 Bool_t AliSegmentArray::SetClass(const char *classname)
128 //set class of stored object
130 //delete fClass; not ower of fClass
134 if (fTreeOwner) delete fTree;
149 ::Fatal("AliSegmentArray::AliSegmentArray", "ROOT system not initialized");
151 fClass = gROOT->GetClass(classname);
153 Error("AliSegmentArray", "%s is not a valid class name", classname);
156 if (!fClass->InheritsFrom(AliSegmentID::Class())) {
157 Error("AliSegmentArray", "%s does not inherit from AliSegmentID", classname);
164 AliSegmentID * AliSegmentArray::NewSegment()
167 //create object according class information
168 if (fClass==0) return 0;
169 AliSegmentID * segment = (AliSegmentID * )fClass->New();
170 if (segment == 0) return 0;
175 Bool_t AliSegmentArray::AddSegment(AliSegmentID *segment)
178 // add segment to array
180 if (segment==0) return kFALSE;
181 if (fSegment==0) return kFALSE;
182 if (fClass==0) return kFALSE;
183 if (!(segment->IsA()->InheritsFrom(fClass))){
184 Error("AliSegmentArray", "added class %s is not of proper type ",
185 segment->IsA()->GetName());
188 fSegment->AddAt(segment,segment->GetID());
189 fNSegment = fSegment->GetLast()+1;
193 AliSegmentID * AliSegmentArray::AddSegment(Int_t index)
196 // add segment to array
198 if (fSegment==0) return 0;
199 if (fClass==0) return 0;
200 AliSegmentID * segment = NewSegment();
201 if (segment == 0) return 0;
202 fSegment->AddAt(segment,index);
203 segment->SetID(index);
204 fNSegment = fSegment->GetLast()+1;
209 void AliSegmentArray::ClearSegment(Int_t index)
212 //remove segment from active memory
214 //PH if ((*fSegment)[index]){
215 if (fSegment->At(index)){
216 // (*fSegment)[index]->Delete(); //not working for TClonesArray
217 //PH delete (*fSegment)[index]; //because problem with deleting TClonesArray
218 //PH fSegment->RemoveAt(index);
219 delete fSegment->RemoveAt(index);
224 Bool_t AliSegmentArray::MakeArray(Int_t n)
227 //make array of pointers to Segments
233 fSegment = new TObjArray(n);
235 if (fSegment) return kTRUE;
238 void AliSegmentArray::MakeTree(TTree* tree)
240 //Make tree with the name
241 AliSegmentID * psegment = NewSegment();
243 //PH fBranch = fTree->Branch("Segment",psegment->IsA()->GetName(),&psegment,64000);
244 fBranch = fTree->Branch("Segment",psegment->IsA()->GetName(),&psegment,64000,99);
248 void AliSegmentArray::MakeTree(char *file)
250 // AliSegmentID segment;
251 AliSegmentID * psegment = NewSegment();
256 fTree = new TTree("Segment Tree","Tree with segments");
258 else fTree->Reset();}
261 //PH fBranch = fTree->Branch("Segment",psegment->IsA()->GetName(),&psegment,64000);
262 fBranch = fTree->Branch("Segment",psegment->IsA()->GetName(),&psegment,64000,99);
265 TString outFile = gAlice->GetBaseFile();
266 outFile = outFile + "/" + file;
267 fBranch->SetFile(outFile.Data());
268 TDirectory *wd = gDirectory;
269 TBranch *b = fBranch;
270 TIter next( b->GetListOfBranches());
271 while ((b=(TBranch*)next())) {
272 b->SetFile(outFile.Data());
274 cout << "Diverting branch " << "Segment" << " to file " << outFile << endl;
281 Bool_t AliSegmentArray::MakeDictionary(Int_t size)
284 //create index table for tree
286 if (size<1) return kFALSE;
287 if (fTreeIndex) delete fTreeIndex;
288 fTreeIndex = new TArrayI();
289 fTreeIndex->Set(size);
291 AliSegmentID * psegment = NewSegment(); //MI change
292 fBranch->SetAddress(&psegment);
293 TBranch * brindix = fTree->GetBranch("fSegmentID");
294 Int_t nevent = (Int_t)fTree->GetEntries();
295 for (Int_t i = 0; i<nevent; i++){
296 brindix->GetEvent(i);
297 Int_t treeIndex=psegment->GetID();
298 if (fTreeIndex->fN<treeIndex) fTreeIndex->Set(Int_t(Float_t(treeIndex)*1.5)+1);
299 // Int_t index = segment.GetID();
300 (*fTreeIndex)[treeIndex]=i+1; //
302 if (psegment) delete psegment;
306 Bool_t AliSegmentArray::ConnectTree(TTree* tree)
309 if (fTree == 0) return kFALSE;
310 fBranch = fTree->GetBranch("Segment");
311 if (fBranch==0) return kFALSE;
312 MakeDictionary(TMath::Max(fNSegment,Int_t(fTree->GetEntries())));
313 MakeArray(fTreeIndex->fN);
318 Bool_t AliSegmentArray::ConnectTree(const char * treeName)
320 //connect tree from current directory
330 fTree =(TTree*)gDirectory->Get(treeName);
332 if (fTree == 0) return kFALSE;
333 fBranch = fTree->GetBranch("Segment");
334 if (fBranch==0) return kFALSE;
335 MakeDictionary(TMath::Max(fNSegment,Int_t(fTree->GetEntries())));
336 MakeArray(fTreeIndex->fN);
340 AliSegmentID *AliSegmentArray::LoadSegment(Int_t index)
343 //load segment with index to the memory
346 if (fTreeIndex ==0 ) MakeDictionary(3000);
347 //firstly try to load dictionary
348 if (fTreeIndex ==0 ) return 0;
349 if (fBranch==0) return 0;
350 if (index>fTreeIndex->fN) return 0;
351 //PH AliSegmentID *s = (AliSegmentID*)(*fSegment)[index];
352 AliSegmentID *s = (AliSegmentID*)fSegment->At(index);
353 if (s==0) s= NewSegment();
355 // new AliSegmentID(index);
358 Int_t treeIndex =(*fTreeIndex)[index];
359 if (treeIndex<1) return 0;
360 else treeIndex--; //I don't like it Int table I have index shifted by 1
361 fBranch->SetAddress(&s);
362 fTree->GetEvent(treeIndex);
363 //PH (*fSegment)[index] = (TObject*) s;
364 fSegment->AddAt((TObject*) s, index);
371 AliSegmentID *AliSegmentArray::LoadEntry(Int_t index)
374 //load segment at position inex in tree to the memory
377 if (fBranch==0) return 0;
378 if (index>fTree->GetEntries()) return 0;
379 AliSegmentID * s = NewSegment();
382 fBranch->SetAddress(&s);
383 fTree->GetEvent(index);
387 Int_t nindex = s->GetID();
388 ClearSegment(nindex);
389 //PH (*fSegment)[nindex] = (TObject*) s;
390 fSegment->AddAt((TObject*) s, nindex);
394 void AliSegmentArray::StoreSegment(Int_t index)
397 //make segment persistent
399 const AliSegmentID * ksegment = (*this)[index];
400 if (ksegment == 0 ) return;
401 if (fTree==0) MakeTree();
402 fBranch->SetAddress(&ksegment);
407 void AliSegmentArray::Streamer(TBuffer &R__b)
409 TObjString treeName, * ptreeName=&treeName;
410 if (R__b.IsReading()) {
411 Version_t R__v = R__b.ReadVersion(); if (R__v) { }
412 TNamed::Streamer(R__b);
414 if (fTree && fTreeOwner) delete fTree;
415 ConnectTree(ptreeName->String());
417 R__b.WriteVersion(AliSegmentArray::IsA());
418 TNamed::Streamer(R__b);
420 // sprintf(ch,"%s",fTrre->GetTitle());
421 treeName.String() = fTree->GetTitle();