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 **************************************************************************/
17 /////////////////////////////////////////////////////////////////////////
19 // Class for ITS RecPoint reconstruction //
21 ////////////////////////////////////////////////////////////////////////
25 #include "AliRunLoader.h"
26 #include "AliITSDetTypeRec.h"
27 #include "AliITSLoader.h"
28 #include "AliITSreconstruction.h"
29 #include "AliITSgeom.h"
32 ClassImp(AliITSreconstruction)
34 //______________________________________________________________________
35 AliITSreconstruction::AliITSreconstruction():
44 // Default constructor.
50 // A zero-ed constructed AliITSreconstruction class.
51 fDet[0] = fDet[1] = fDet[2] = kTRUE;
53 //______________________________________________________________________
55 AliITSreconstruction::AliITSreconstruction(AliRunLoader *rl):
64 fDet[0] = fDet[1] = fDet[2] = kTRUE;
66 //______________________________________________________________________
67 AliITSreconstruction::AliITSreconstruction(const char* filename):
76 // Standard constructor.
78 // const char* filename filename containing the digits to be
79 // reconstructed. If filename = 0 (nil)
80 // then no file is opened but a file is
81 // assumed to already be opened. This
82 // already opened file will be used.
86 // A standardly constructed AliITSreconstruction class.
88 fDet[0] = fDet[1] = fDet[2] = kTRUE;
90 fRunLoader = AliRunLoader::Open(filename);
91 if (fRunLoader == 0x0)
93 Error("AliITSreconstruction","Can not load the session",filename);
96 fRunLoader->LoadgAlice();
97 gAlice = fRunLoader->GetAliRun();
100 Error("AliITSreconstruction","gAlice not found on file. Aborting.");
106 //______________________________________________________________________
107 AliITSreconstruction::AliITSreconstruction(const AliITSreconstruction &/*rec*/):TTask(/*rec*/){
110 Error("Copy constructor","Copy constructor not allowed");
113 //______________________________________________________________________
114 AliITSreconstruction& AliITSreconstruction::operator=(const AliITSreconstruction& /*source*/){
115 // Assignment operator. This is a function which is not allowed to be
117 Error("operator=","Assignment operator not allowed\n");
121 //______________________________________________________________________
122 AliITSreconstruction::~AliITSreconstruction(){
123 // A destroyed AliITSreconstruction class.
129 //______________________________________________________________________
130 Bool_t AliITSreconstruction::Init(){
137 // kTRUE if no errors initilizing this class occurse else kFALSE
139 if (fRunLoader == 0x0)
141 Error("Init","Run Loader is NULL");
144 fRunLoader->LoadgAlice();
145 fRunLoader->LoadHeader();
147 fLoader = (AliITSLoader*) fRunLoader->GetLoader("ITSLoader");
149 Error("Init","ITS loader not found");
153 // Now ready to init.
155 fRunLoader->CdGAFile();
156 fITSgeom = (AliITSgeom*)gDirectory->Get("AliITSgeom");
158 fDetTypeRec = new AliITSDetTypeRec();
159 fDetTypeRec->SetITSgeom(fITSgeom);
160 fDetTypeRec->SetDefaults();
161 fDet[0] = fDet[1] = fDet[2] = kTRUE;
164 //fEnt = gAlice->GetEventsPerRun();
165 fEnt = Int_t(fRunLoader->TreeE()->GetEntries());
167 fLoader->LoadDigits("read");
168 fLoader->LoadRecPoints("recreate");
169 fLoader->LoadRawClusters("recreate");
170 if (fLoader->TreeR() == 0x0) fLoader->MakeTree("R");
171 if (fLoader->TreeC() == 0x0) fLoader->MakeTree("C");
173 fDetTypeRec->SetLoader(fLoader);
174 fDetTypeRec->MakeBranchR(0);
175 fDetTypeRec->MakeBranchC();
176 fDetTypeRec->SetTreeAddress();
177 fDetTypeRec->SetTreeAddressR(fLoader->TreeR());
181 Info("Init"," Done\n\n\n");
185 //______________________________________________________________________
186 Bool_t AliITSreconstruction::InitRec(){
187 // Sets up Reconstruction part of AliITSDetType..
195 //AliITSDetType *idt;
196 fDetTypeRec->SetLoader(fLoader);
199 Info("InitRec","SPD");
200 //idt = fITS->DetType(kSPD);
201 AliITSsegmentationSPD *segSPD = (AliITSsegmentationSPD*)fDetTypeRec->GetSegmentationModel(0);
202 TClonesArray *digSPD = fDetTypeRec->DigitsAddress(kSPD);
203 TClonesArray *recpSPD = fDetTypeRec->ClustersAddress(kSPD);
204 Info("InitRec","idt = %#x; digSPD = %#x; recpSPD = %#x",fDetTypeRec,digSPD,recpSPD);
205 AliITSClusterFinderSPD *recSPD = new AliITSClusterFinderSPD(segSPD,digSPD,recpSPD);
206 fDetTypeRec->SetReconstructionModel(kSPD,recSPD);
207 } // end if fDet[kSPD].
210 Info("InitRec","SDD");
211 // idt = fITS->DetType(kSDD);
212 AliITSsegmentationSDD *segSDD = (AliITSsegmentationSDD*)
213 fDetTypeRec->GetSegmentationModel(1);
214 AliITSresponseSDD *resSDD = (AliITSresponseSDD*)
215 fDetTypeRec->GetCalibrationModel(fDetTypeRec->GetITSgeom()->GetStartSDD());
216 TClonesArray *digSDD = fDetTypeRec->DigitsAddress(kSDD);
217 TClonesArray *recpSDD = fDetTypeRec->ClustersAddress(kSDD);
218 AliITSClusterFinderSDD *recSDD =new AliITSClusterFinderSDD(segSDD,
221 fDetTypeRec->SetReconstructionModel(kSDD,recSDD);
222 } // end if fDet[kSDD]
225 Info("InitRec","SSD");
226 //idt = fITS->DetType(kSSD);
227 AliITSsegmentationSSD *segSSD = (AliITSsegmentationSSD*)
228 fDetTypeRec->GetSegmentationModel(2);
229 TClonesArray *digSSD = fDetTypeRec->DigitsAddress(kSSD);
230 AliITSClusterFinderSSD *recSSD =new AliITSClusterFinderSSD(segSSD,
232 recSSD->SetITSgeom(fDetTypeRec->GetITSgeom());
233 fDetTypeRec->SetReconstructionModel(kSSD,recSSD);
234 } // end if fDet[kSSD]
236 fDetTypeRec->SetDefaultClusterFinders();
237 Info("InitRec"," Done\n");
240 //______________________________________________________________________
241 void AliITSreconstruction::Exec(const Option_t *opt){
242 // Main reconstruction function.
244 // Option_t * opt list of subdetector to digitize. =0 all.
252 if(strstr(opt,"All")||strstr(opt,"ALL")||strstr(opt,"ITS")||opt==0){
253 fDet[0] = fDet[1] = fDet[2] = kTRUE;
256 fDet[0] = fDet[1] = fDet[2] = kFALSE;
257 if(strstr(opt,"SPD")) fDet[kSPD] = kTRUE;
258 if(strstr(opt,"SDD")) fDet[kSDD] = kTRUE;
259 if(strstr(opt,"SSD")) fDet[kSSD] = kTRUE;
260 if(fDet[kSPD] && fDet[kSDD] && fDet[kSSD]) lopt = "All";
262 } // end if strstr(opt,...)
265 cout << "Initilization Failed, Can't run Exec." << endl;
268 for(evnt=0;evnt<fEnt;evnt++)
271 Info("Exec","Processing Event %d",evnt);
274 fRunLoader->GetEvent(evnt);
275 fDetTypeRec->SetITSgeom(fITSgeom);
276 fDetTypeRec->SetLoader(fLoader);
277 if (fLoader->TreeR() == 0x0) fLoader->MakeTree("R");
278 fDetTypeRec->MakeBranchR(0);
279 if (fLoader->TreeC() == 0x0){
280 fDetTypeRec->MakeTreeC();
281 fDetTypeRec->MakeBranchC();
283 fDetTypeRec->SetTreeAddressR(fLoader->TreeR());
284 fDetTypeRec->SetTreeAddressD(fLoader->TreeD());
285 fDetTypeRec->DigitsToRecPoints(evnt,0,lopt);
288 //______________________________________________________________________
289 void AliITSreconstruction::SetOutputFile(TString filename){
290 // Set a new file name for recpoints.
291 // It must be called before Init()
292 if(!fLoader)fLoader = (AliITSLoader*) fRunLoader->GetLoader("ITSLoader");
294 Info("SetOutputFile","name for rec points is %s",filename.Data());
295 fLoader->SetRecPointsFileName(filename);
298 Error("SetOutputFile",
299 "ITS loader not available. Not possible to set name: %s",filename.Data());