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 // class for TPC reconstruction //
22 ///////////////////////////////////////////////////////////////////////////////
25 #include "AliTPCReconstructor.h"
26 #include "AliRunLoader.h"
28 #include "AliRawReader.h"
29 #include "AliTPCclustererMI.h"
30 #include "AliTPCtrackerMI.h"
31 #include "AliTPCpidESD.h"
32 #include "AliTPCParam.h"
33 #include "AliTPCParamSR.h"
34 #include "AliTPCcalibDB.h"
36 ClassImp(AliTPCReconstructor)
39 AliTPCRecoParam * AliTPCReconstructor::fgkRecoParam =0; // reconstruction parameters
40 Int_t AliTPCReconstructor::fgStreamLevel = 0; // stream (debug) level
43 AliTPCReconstructor::AliTPCReconstructor():
48 // default constructor
53 AliTPCParam* param = GetTPCParam();
55 AliWarning("Loading default TPC parameters !");
56 param = new AliTPCParamSR;
58 fClusterer = new AliTPCclustererMI(param);
61 //_____________________________________________________________________________
62 AliTPCReconstructor::~AliTPCReconstructor()
64 if (fgkRecoParam) delete fgkRecoParam;
65 if (fClusterer) delete fClusterer;
68 //_____________________________________________________________________________
69 void AliTPCReconstructor::Reconstruct(TTree* digitsTree, TTree* clustersTree) const {
70 // single event local reconstruction
72 fClusterer->SetInput(digitsTree);
73 fClusterer->SetOutput(clustersTree);
74 fClusterer->Digits2Clusters();
77 //_____________________________________________________________________________
78 void AliTPCReconstructor::Reconstruct(AliRawReader* rawReader, TTree* clustersTree) const {
79 // single event local reconstruction
80 // of TPC data starting from raw data
81 TString option = GetOption();
82 if (option.Contains("OldRCUFormat"))
83 fClusterer->SetOldRCUFormat(kTRUE);
85 fClusterer->SetOutput(clustersTree);
86 fClusterer->Digits2Clusters(rawReader);
89 //_____________________________________________________________________________
90 AliTracker* AliTPCReconstructor::CreateTracker() const
92 // create a TPC tracker
94 AliTPCParam* param = GetTPCParam();
96 AliWarning("Loading default TPC parameters !");
97 param = new AliTPCParamSR;
99 param->ReadGeoMatrices();
100 return new AliTPCtrackerMI(param);
103 //_____________________________________________________________________________
104 void AliTPCReconstructor::FillESD(TTree */*digitsTree*/, TTree */*clustersTree*/,
105 AliESDEvent* esd) const
109 Double_t parTPC[] = {47., 0.07, 5.};
110 AliTPCpidESD tpcPID(parTPC);
115 //_____________________________________________________________________________
116 AliTPCParam* AliTPCReconstructor::GetTPCParam() const
118 // get the TPC parameters
120 AliTPCParam* param = AliTPCcalibDB::Instance()->GetParameters();
128 const AliTPCRecoParam* AliTPCReconstructor::GetRecoParam(){
130 // Get reconstruction parameters
135 // 1. try to get reco parameters from OCDB
137 fgkRecoParam = AliTPCcalibDB::Instance()->GetRecoParam(0);
138 //Info("","Reconstruction parameters from OCDB used");
140 // 2. If not initialized take default
143 fgkRecoParam = AliTPCRecoParam::GetHighFluxParam();
144 //Error("","Default reconstruction parameters used");