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 //-----------------------------------------------------------------------------
19 /// \class AliMUONClusterStoreV1
21 /// Implementation of VClusterStore.
23 /// This one is a basic implementation, let's say "legacy" one, i.e.
24 /// compatible with what we stored in MUON.RecPoints.root files before
25 /// the switch to data stores.
27 /// \author Laurent Aphecetche, Subatech
29 //-----------------------------------------------------------------------------
31 #include "AliMUONClusterStoreV1.h"
34 #include "AliMUONRawCluster.h"
35 #include "AliMUONTOTCAStoreIterator.h"
36 #include "AliMUONTreeManager.h"
37 #include "AliMpConstants.h"
38 #include "AliMpDEManager.h"
39 #include <TClonesArray.h>
40 #include <TObjArray.h>
44 ClassImp(AliMUONClusterStoreV1)
47 //_____________________________________________________________________________
48 AliMUONClusterStoreV1::AliMUONClusterStoreV1()
49 : AliMUONVClusterStore(),
50 fClusters(new TObjArray(AliMpConstants::NofChambers()))
52 /// ctor. Set correct ownerships
53 fClusters->SetOwner(kTRUE);
54 for ( Int_t i = 0; i < fClusters->GetSize(); ++i )
56 TClonesArray* tca = new TClonesArray("AliMUONRawCluster",100);
58 fClusters->AddAt(tca,i);
63 //_____________________________________________________________________________
64 AliMUONClusterStoreV1::AliMUONClusterStoreV1(const AliMUONClusterStoreV1&)
65 : AliMUONVClusterStore(),
69 AliError("Please implement me");
72 //_____________________________________________________________________________
73 AliMUONClusterStoreV1&
74 AliMUONClusterStoreV1::operator=(const AliMUONClusterStoreV1&)
76 /// assignment operator
77 AliError("Please implement me");
81 //_____________________________________________________________________________
82 AliMUONClusterStoreV1::~AliMUONClusterStoreV1()
89 //_____________________________________________________________________________
91 AliMUONClusterStoreV1::Add(const AliMUONRawCluster& cluster)
93 /// Add a cluster to this store
94 Int_t iChamber = AliMpDEManager::GetChamberId(cluster.GetDetElemId());
95 TClonesArray* array = ChamberClusters(iChamber);
100 new((*array)[array->GetLast()+1]) AliMUONRawCluster(cluster);
104 //_____________________________________________________________________________
106 AliMUONClusterStoreV1::ChamberClusters(Int_t chamberId) const
108 /// Get the internal array of clusters for a given chamber
109 TClonesArray* array = static_cast<TClonesArray*>(fClusters->At(chamberId));
112 AliError(Form("Cannot get Clusters for chamberId=%d",chamberId));
118 //_____________________________________________________________________________
120 AliMUONClusterStoreV1::ChamberClustersPtr(Int_t chamberId) const
122 /// Get the internal array of clusters for a given chamber
123 TClonesArray* array = static_cast<TClonesArray*>(fClusters->At(chamberId));
126 AliError(Form("Cannot get Clusters for chamberId=%d",chamberId));
129 return fClusters->GetObjectRef(array);
132 //_____________________________________________________________________________
134 AliMUONClusterStoreV1::Connect(TTree& tree, Bool_t alone) const
136 /// Connect this to the tree, i.e. make the branches or set their addresses.
138 AliMUONTreeManager tman;
141 TBranch* b = tree.GetBranch("MUONRawClusters1");
143 Bool_t isMaking = (b == 0);
147 for ( Int_t i = 0; i < AliMpConstants::NofTrackingChambers(); ++i )
149 TString branchName(Form("MUONRawClusters%d",i+1));
150 ok = ok && tman.MakeBranch(tree,ClassName(),"TClonesArray",
151 branchName.Data(),ChamberClustersPtr(i));
157 if (alone) tman.UpdateBranchStatuses(tree,"MUONRawClusters");
158 for ( Int_t i = 0; i < AliMpConstants::NofTrackingChambers(); ++i )
160 TString branchName(Form("MUONRawClusters%d",i+1));
161 ok = ok && tman.SetAddress(tree,branchName.Data(),
162 ChamberClustersPtr(i));
168 //_____________________________________________________________________________
170 AliMUONClusterStoreV1::Remove(AliMUONRawCluster& cluster)
173 Int_t iChamber = AliMpDEManager::GetChamberId(cluster.GetDetElemId());
174 TClonesArray* array = ChamberClusters(iChamber);
175 TObject* o = array->Remove(&cluster);
180 return static_cast<AliMUONRawCluster*>(o);
183 //_____________________________________________________________________________
185 AliMUONClusterStoreV1::Clear(Option_t*)
187 /// Reset internal arrays
189 /// Reset the tclonesarray, but keep the tobjarray's size constant.
190 for ( Int_t i = 0; i < fClusters->GetSize(); ++i )
192 ChamberClusters(i)->Clear("C");
196 //_____________________________________________________________________________
198 AliMUONClusterStoreV1::CreateIterator() const
200 /// Return an iterator to loop over our clusters
201 return new AliMUONTOTCAStoreIterator(fClusters,0,AliMpConstants::NofTrackingChambers()-1);
204 //_____________________________________________________________________________
206 AliMUONClusterStoreV1::CreateChamberIterator(Int_t firstChamber, Int_t lastChamber) const
208 /// Return an iterator to loop over our clusters
209 return new AliMUONTOTCAStoreIterator(fClusters,firstChamber,lastChamber);
212 //_____________________________________________________________________________
214 AliMUONClusterStoreV1::GetSize() const
216 /// Return the number of clusters we hold
218 for ( Int_t i = 0; i < fClusters->GetSize(); ++i )
220 n += ChamberClusters(i)->GetLast()+1;