1 #include "AliITSdigit.h"
2 #include "AliITSLoader.h"
3 #include "AliRunLoader.h"
6 ///////////////////////////////////////////////////////////////////////////
8 // it manages the I/O for:
9 // raw clusters, primary vertices
11 // and tracks propagated to the origin
12 //////////////////////////////////////////////////////////////////////////
13 const TString AliITSLoader::fgkDefaultRawClustersContainerName = "TreeC";
14 const TString AliITSLoader::fgkDefaultBackTracksContainerName = "TreeB";
15 const TString AliITSLoader::fgkDefaultVerticesContainerName = "Vertex";
16 const TString AliITSLoader::fgkDefaultV0ContainerName = "V0";
17 const TString AliITSLoader::fgkDefaultCascadeContainerName = "Cascade";
18 ClassImp(AliITSLoader)
20 /**********************************************************************/
21 AliITSLoader::AliITSLoader():AliLoader(),
24 // Default constructor
26 /*********************************************************************/
27 AliITSLoader::AliITSLoader(const Char_t *name,const Char_t *topfoldername):
28 AliLoader(name,topfoldername),
32 AliDataLoader* rawClustersDataLoader = new AliDataLoader(
33 fDetectorName + ".RawCl.root",fgkDefaultRawClustersContainerName,
35 fDataLoaders->Add(rawClustersDataLoader);
36 rawClustersDataLoader->SetEventFolder(fEventFolder);
37 rawClustersDataLoader->SetFolder(GetDetectorDataFolder());
39 AliDataLoader* backTracksDataLoader = new AliDataLoader(
40 fDetectorName + ".BackTracks.root",fgkDefaultBackTracksContainerName,
41 "Back Propagated Tracks");
42 fDataLoaders->Add(backTracksDataLoader);
43 backTracksDataLoader->SetEventFolder(fEventFolder);
44 backTracksDataLoader->SetFolder(GetDetectorDataFolder());
46 AliDataLoader* vertexDataLoader = new AliDataLoader(
47 fDetectorName + ".Vertex.root",fgkDefaultVerticesContainerName,
48 "Primary Vertices","O");
49 fDataLoaders->Add(vertexDataLoader);
50 vertexDataLoader->SetEventFolder(fEventFolder);
51 vertexDataLoader->SetFolder(GetDetectorDataFolder());
53 AliDataLoader* v0DataLoader = new AliDataLoader(
54 fDetectorName + ".V0s.root",fgkDefaultV0ContainerName,"V0 Vertices");
55 fDataLoaders->Add(v0DataLoader);
56 v0DataLoader->SetEventFolder(fEventFolder);
57 v0DataLoader->SetFolder(GetDetectorDataFolder());
59 AliDataLoader* cascadeDataLoader = new AliDataLoader(
60 fDetectorName + ".Cascades.root",fgkDefaultCascadeContainerName,
62 fDataLoaders->Add(cascadeDataLoader);
63 cascadeDataLoader->SetEventFolder(fEventFolder);
64 cascadeDataLoader->SetFolder(GetDetectorDataFolder());
66 /**********************************************************************/
67 AliITSLoader::AliITSLoader(const Char_t *name,TFolder *topfolder):
68 AliLoader(name,topfolder),
72 AliDataLoader* rawClustersDataLoader = new AliDataLoader(
73 fDetectorName + ".RawCl.root",fgkDefaultRawClustersContainerName,
75 fDataLoaders->Add(rawClustersDataLoader);
76 rawClustersDataLoader->SetEventFolder(fEventFolder);
77 rawClustersDataLoader->SetFolder(GetDetectorDataFolder());
79 AliDataLoader* backTracksDataLoader = new AliDataLoader(
80 fDetectorName + ".BackTracks.root",fgkDefaultBackTracksContainerName,
81 "Back Propagated Tracks");
82 fDataLoaders->Add(backTracksDataLoader);
83 backTracksDataLoader->SetEventFolder(fEventFolder);
84 backTracksDataLoader->SetFolder(GetDetectorDataFolder());
86 AliDataLoader* vertexDataLoader = new AliDataLoader(
87 fDetectorName + ".Vertex.root",fgkDefaultVerticesContainerName,
88 "Primary Vertices","O");
89 fDataLoaders->Add(vertexDataLoader);
90 vertexDataLoader->SetEventFolder(fEventFolder);
91 vertexDataLoader->SetFolder(GetDetectorDataFolder());
93 AliDataLoader* v0DataLoader = new AliDataLoader(
94 fDetectorName + ".V0.root",fgkDefaultV0ContainerName,"V0 Vertices");
95 fDataLoaders->Add(v0DataLoader);
96 v0DataLoader->SetEventFolder(fEventFolder);
97 v0DataLoader->SetFolder(GetDetectorDataFolder());
99 AliDataLoader* cascadeDataLoader = new AliDataLoader(
100 fDetectorName + ".Cascade.root",fgkDefaultCascadeContainerName,
102 fDataLoaders->Add(cascadeDataLoader);
103 cascadeDataLoader->SetEventFolder(fEventFolder);
104 cascadeDataLoader->SetFolder(GetDetectorDataFolder());
108 /**********************************************************************/
109 AliITSLoader::~AliITSLoader(){
111 AliDataLoader* dl = 0;
113 dl = GetRawClLoader();
114 fDataLoaders->Remove(dl);
117 dl = GetBackTracksDataLoader();
118 fDataLoaders->Remove(dl);
121 dl = GetVertexDataLoader();
122 fDataLoaders->Remove(dl);
125 dl = GetV0DataLoader();
126 fDataLoaders->Remove(dl);
129 dl = GetCascadeDataLoader();
130 fDataLoaders->Remove(dl);
132 if(fGeom)delete fGeom;
136 //----------------------------------------------------------------------
137 AliITS* AliITSLoader::GetITS(){
138 // Returns the pointer to the ITS, kept on the file. A short cut metthod
144 // Returns a pointer to the ITS, if not found returns 0.
148 its = dynamic_cast<AliITS*> (gAlice->GetDetector(
149 GetDetectorName().Data()));
156 ar = rl->GetAliRun();
158 its = dynamic_cast<AliITS*> (ar->GetDetector(GetDetectorName().Data()));
161 //----------------------------------------------------------------------
162 void AliITSLoader::SetupDigits(AliITS *its){
163 // Sets up to store ITS Digits in side AliITS::fDtype TObjArray
165 // AliITS *its Pointer to the ITS
171 its->SetTreeAddressD(TreeD());
174 //----------------------------------------------------------------------
175 void AliITSLoader::SetupDigits(TObjArray *digPerDet,Int_t n,
176 const Char_t **digclass){
177 // Sets up digPerDet to store ITS Digits.
179 // TObjArray *digPerDet A pointer to a TObject Array size>=3.
180 // Int_t n The size of the TObjArray and digclass array
181 // Char_t **digclass Array of digit class names
183 // TObjArray *digPerDet Setup and linked to the tree of digits
187 TClonesArray *cl = 0;
191 const Char_t *det[3] = {"SPD","SDD","SSD"};
194 Error("SetUpDigits","TObject Array digPerDet does not exist");
197 m = digPerDet->GetSize();
199 Error("SetUpDigits","TObject Array digPerDet=%p must have a size"
200 " at least that of n=%d",digPerDet,n);
203 Error("SetUpDigits","TObject Array digPerDet=%p must have a size >2",
209 if(digPerDet->At(i)==0){ // set up TClones Array
210 digPerDet->AddAt(new TClonesArray(digclass[i],1000),i);
211 if(n==3) sprintf(branch,"ITSDigits%s",det[i]);
212 else sprintf(branch,"ITSDigits%d",i+1);
213 br = td->GetBranch(branch);
214 br->SetAddress(&((*digPerDet)[i]));
215 continue; // do next one.
217 cl = dynamic_cast<TClonesArray*> (digPerDet->At(i));
218 if(!cl && digPerDet->At(i)!=0){ // not a TClonesArray
219 Error("SetUpDigits","TObject Array digPerDet-At(%d)=%p must be "
220 "zeroed or filled with TClonesArrays",i,digPerDet);
223 if(!(cl->GetClass()->GetBaseClass(AliITSdigit::Class()))){
224 Error("SetUPDigits","TClones array at digPerDet[%d}=%p must be"
225 "derived from AliITSdigit",i,digPerDet->At(i));
228 if(n==3) sprintf(branch,"ITSDigits%s",det[i]);
229 else sprintf(branch,"ITSDigits%d",i+1);
230 br = td->GetBranch(branch);
231 br->SetAddress(&((*digPerDet)[i]));
235 //---------------------------------------------------------------------
236 AliITSdigit * AliITSLoader::GetDigit(TObjArray *digPerDet,Int_t module,
238 // Gets the digit for for a specific detector type and module.
239 // To be used in conjustion with Setupdigits(AliITS *its).
241 // TObjArray *digPereDet Pointer to the Array of digits
242 // Int_t module Module number
243 // Int_t digit Digit number
247 // returns the pointer to the digit. if zero then last digit for that
251 Error("GetDigit","digPerDet=%p, module=%d, digit=%d",
252 digPerDet,module,digit);
258 //---------------------------------------------------------------------
259 AliITSdigit * AliITSLoader::GetDigit(AliITS *its,Int_t module,Int_t digit){
260 // Gets the digit for for a specific detector type and module.
261 // To be used in conjustion with Setupdigits(AliITS *its).
263 // AliITS *its Pointer to the ITS
264 // Int_t module Module number
265 // Int_t digit Digit number
269 // returns the pointer to the digit. if zero then last digit for that
271 //AliITSDetType *idtype;
272 AliITSgeom *geom = its->GetITSgeom();
273 Int_t idet = geom->GetModuleType(module);
274 TClonesArray *digits;
277 TreeD()->GetEvent(module);
278 digits = its->DigitsAddress(idet);
279 if(digit>-1 && digit<digits->GetEntriesFast()){ // if in range.
280 return (AliITSdigit*) digits->At(digit);
285 //----------------------------------------------------------------------
286 void AliITSLoader::MakeTree(Option_t *opt){
287 // invokes AliLoader::MakeTree + specific ITS tree(s)
288 // Valid options: H,S,D,R,T and C (C=raw clusters)
289 AliLoader::MakeTree(opt);
290 const char *oC = strstr(opt,"C");
291 if (oC) MakeRawClustersContainer();
293 const char *oB = strstr(opt,"B");
294 if (oB) MakeBackTracksContainer();
296 const char *oV0 = strstr(opt,"V0");
297 if (oV0) MakeV0Container();
299 const char *oX = strstr(opt,"X");
300 if (oX) MakeCascadeContainer();
303 //----------------------------------------------------------------------
304 AliITSgeom* AliITSLoader::GetITSgeom(Bool_t force) {
305 // retrieves the ITS geometry from file
306 if(fGeom && !force)return fGeom;
311 AliRunLoader *runLoader = GetRunLoader();
312 if (!runLoader->GetAliRun()) runLoader->LoadgAlice();
313 if (!runLoader->GetAliRun()) {
314 Error("GetITSgeom", "couldn't get AliRun object");
318 TDirectory *curdir = gDirectory;
319 runLoader->CdGAFile();
320 fGeom = (AliITSgeom*)gDirectory->Get("AliITSgeom");
323 Error("GetITSgeom","no ITS geometry available");
326 AliDebug(1,"AliITSgeom object has been fetched from galice.root file");
329 //______________________________________________________________________
330 void AliITSLoader::SetITSgeom(AliITSgeom *geom){
331 // Replaces the AliITSgeom object read from file with the one
334 // AliITSgeom *geom The AliITSgeom object to replace the one
335 // read from the file
341 if(fGeom==geom) return; // Same do nothing