#include <TStopwatch.h>
-#include "AliL3StandardIncludes.h"
-#include "AliL3Logging.h"
-#include "AliL3Transform.h"
-#include "AliL3Hough.h"
+#include "AliHLTStandardIncludes.h"
+#include "AliHLTLogging.h"
+#include "AliHLTTransform.h"
+#include "AliHLTHough.h"
#include "AliLog.h"
-#include <AliL3ITSclusterer.h>
-#include <AliL3ITSVertexerZ.h>
-#include <AliL3ITStracker.h>
+#include <AliHLTITSclusterer.h>
+#include <AliHLTITSVertexerZ.h>
+#include <AliHLTITStracker.h>
#include "AliHoughFilter.h"
#include <AliMagFMaps.h>
#include <AliKalmanTrack.h>
#include <AliITSgeom.h>
+#include <AliESDVertex.h>
ClassImp(AliHoughFilter)
//_____________________________________________________________________________
-AliHoughFilter::AliHoughFilter()
+AliHoughFilter::AliHoughFilter():
+fPtmin(0),
+fITSgeom(NULL)
{
// default constructor
// Init debug level
- AliL3Log::fgLevel = AliL3Log::kError;
- if (AliDebugLevel() > 0) AliL3Log::fgLevel = AliL3Log::kWarning;
- if (AliDebugLevel() > 1) AliL3Log::fgLevel = AliL3Log::kInformational;
- if (AliDebugLevel() > 2) AliL3Log::fgLevel = AliL3Log::kDebug;
+ AliHLTLog::fgLevel = AliHLTLog::kError;
+ if (AliDebugLevel() > 0) AliHLTLog::fgLevel = AliHLTLog::kWarning;
+ if (AliDebugLevel() > 1) AliHLTLog::fgLevel = AliHLTLog::kInformational;
+ if (AliDebugLevel() > 2) AliHLTLog::fgLevel = AliHLTLog::kDebug;
// Init TPC HLT geometry
const char *path = gSystem->Getenv("ALICE_ROOT");
Char_t pathname[1024];
strcpy(pathname,path);
strcat(pathname,"/HLT/src");
- if (!AliL3Transform::Init(pathname, kFALSE))
+ if (!AliHLTTransform::Init(pathname, kFALSE))
AliError("HLT initialization failed!");
// Init magnetic field
AliMagF* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., AliMagFMaps::k5kG);
AliTracker::SetFieldMap(field,kTRUE);
- fPtmin = 0.1*AliL3Transform::GetSolenoidField();
+ fPtmin = 0.1*AliHLTTransform::GetSolenoidField();
// Init ITS geometry
fITSgeom = new AliITSgeom();
}
//_____________________________________________________________________________
-Bool_t AliHoughFilter::Filter(AliRawEvent* event, AliESD* esd)
+Bool_t AliHoughFilter::Filter(AliRawEvent* event, AliESDEvent* esd)
{
// Run fast online reconstruction
// based on the HLT tracking algorithms
TStopwatch timer;
timer.Start();
- if(!fITSgeom)
- AliError("ITS geometry not created!");
- AliL3ITSclusterer clusterer(fITSgeom);
+ AliHLTITSclusterer clusterer(0);
AliRawReader *itsrawreader=new AliRawReaderRoot(event);
clusterer.Digits2Clusters(itsrawreader,treeClusters);
delete itsrawreader;
//_____________________________________________________________________________
-void AliHoughFilter::RunITSvertexer(AliESD* esd, TTree *treeClusters)
+void AliHoughFilter::RunITSvertexer(AliESDEvent* esd, TTree *treeClusters)
{
// Run SPD vertexerZ
// Store the result in the ESD
TStopwatch timer;
timer.Start();
- AliL3ITSVertexerZ vertexer;
+ AliHLTITSVertexerZ vertexer;
AliESDVertex *vertex = vertexer.FindVertexForCurrentEvent(fITSgeom,treeClusters);
esd->SetVertex(vertex);
AliInfo(Form("ITS vertexer has finished in %f seconds\n",timer.RealTime()));
}
//_____________________________________________________________________________
-void AliHoughFilter::RunTPCtracking(AliRawEvent* event, AliESD* esd)
+void AliHoughFilter::RunTPCtracking(AliRawEvent* event, AliESDEvent* esd)
{
// Run hough transform tracking in TPC
// The z of the vertex is taken from the ESD
const AliESDVertex *vertex = esd->GetVertex();
Float_t zvertex = vertex->GetZv();
- AliL3Hough *hough1 = new AliL3Hough();
+ AliHLTHough *hough1 = new AliHLTHough();
hough1->SetThreshold(4);
hough1->CalcTransformerParams(fPtmin);
hough1->Init(100,4,event,zvertex);
hough1->SetAddHistograms();
- AliL3Hough *hough2 = new AliL3Hough();
+ AliHLTHough *hough2 = new AliHLTHough();
hough2->SetThreshold(4);
hough2->CalcTransformerParams(fPtmin);
hough2->StartProcessInThread(18,35);
if(hough1->WaitForThreadFinish())
- ::Fatal("AliL3Hough::WaitForThreadFinish"," Can not join the required thread! ");
+ ::Fatal("AliHLTHough::WaitForThreadFinish"," Can not join the required thread! ");
if(hough2->WaitForThreadFinish())
- ::Fatal("AliL3Hough::WaitForThreadFinish"," Can not join the required thread! ");
+ ::Fatal("AliHLTHough::WaitForThreadFinish"," Can not join the required thread! ");
/* In case we run HLT code in the main thread
for(Int_t slice=0; slice<=17; slice++)
}
//_____________________________________________________________________________
-void AliHoughFilter::RunITStracking(AliESD* esd, TTree *treeClusters)
+void AliHoughFilter::RunITStracking(AliESDEvent* esd, TTree *treeClusters)
{
// Run the ITS tracker
// The tracks from the HT TPC tracking are used as seeds
const AliESDVertex *vertex = esd->GetVertex();
vertex->GetXYZ(vtxPos);
- AliL3ITStracker itsTracker(fITSgeom);
+ AliHLTITStracker itsTracker(0);
itsTracker.SetVertex(vtxPos,vtxErr);
itsTracker.LoadClusters(treeClusters);