#include "AliRawReader.h"
#include "AliTPCclustererMI.h"
#include "AliTPCtrackerMI.h"
-#include "AliTPCpidESD.h"
+//#include "AliTPCpidESD.h"
#include "AliTPCParam.h"
#include "AliTPCParamSR.h"
#include "AliTPCcalibDB.h"
+#include "AliTracker.h"
+#include "AliMagF.h"
ClassImp(AliTPCReconstructor)
-AliTPCRecoParam * AliTPCReconstructor::fgkRecoParam =0; // reconstruction parameters
-Int_t AliTPCReconstructor::fgStreamLevel = 0; // stream (debug) level
+Int_t AliTPCReconstructor::fgStreamLevel = 1; // stream (debug) level
AliTPCReconstructor::AliTPCReconstructor():
//
// default constructor
//
- if (!fgkRecoParam) {
- AliError("The Reconstruction parameters nonitialized - Used default one");
- fgkRecoParam = AliTPCRecoParam::GetHighFluxParam();
- }
-
+ //
+ //
+ AliTPCcalibDB * calib = AliTPCcalibDB::Instance();
+ const AliMagF * field = (AliMagF*)TGeoGlobalMagField::Instance()->GetField();
+ calib->SetExBField(field);
AliTPCParam* param = GetTPCParam();
if (!param) {
AliWarning("Loading default TPC parameters !");
fClusterer = new AliTPCclustererMI(param);
}
+AliTPCReconstructor::AliTPCReconstructor(const AliTPCReconstructor& /*rec*/):
+AliReconstructor(),
+fClusterer(NULL)
+{
+ //
+ // Dummy copu constructor
+ //
+}
+
+AliTPCReconstructor& AliTPCReconstructor::operator=(const AliTPCReconstructor&){
+ //
+ // dummy operator
+ //
+ return *this;
+}
+
//_____________________________________________________________________________
AliTPCReconstructor::~AliTPCReconstructor()
{
- if (fgkRecoParam) delete fgkRecoParam;
if (fClusterer) delete fClusterer;
}
void AliTPCReconstructor::Reconstruct(AliRawReader* rawReader, TTree* clustersTree) const {
// single event local reconstruction
// of TPC data starting from raw data
- TString option = GetOption();
- if (option.Contains("OldRCUFormat"))
- fClusterer->SetOldRCUFormat(kTRUE);
fClusterer->SetOutput(clustersTree);
fClusterer->Digits2Clusters(rawReader);
//_____________________________________________________________________________
void AliTPCReconstructor::FillESD(TTree */*digitsTree*/, TTree */*clustersTree*/,
- AliESDEvent* esd) const
+ AliESDEvent* /*esd*/) const
{
// make PID
-
- Double_t parTPC[] = {47., 0.07, 5.};
+/* Now done in AliESDpid
+ Double_t parTPC[] = {50., 0.07, 5.}; // MIP nnormalized to channel 50 -MI
AliTPCpidESD tpcPID(parTPC);
tpcPID.MakePID(esd);
+*/
}