Some extra lines to read PDC06 data.

[u/mrichter/AliRoot.git] / MONITOR / monitorGDC.cxx

diff --git a/MONITOR/monitorGDC.cxx b/MONITOR/monitorGDC.cxx
index 0b35ee0040162c5186bd7900803c1b3117711762..1c138ab0cbe41b9a779f2368f9a2d212bbba894b 100644 (file)
--- a/MONITOR/monitorGDC.cxx
+++ b/MONITOR/monitorGDC.cxx
@@ -27,31 +27,32 @@
 
 #include <TError.h>
 #include <TSysEvtHandler.h>
-#ifdef DATE_SYS
+#ifdef ALI_DATE
 #include <TROOT.h>
 #include <TSystem.h>
 #include <TDatime.h>
 #include "AliRawReaderDate.h"
 #include "event.h"
 #include "monitor.h"
-#include <AliL3StandardIncludes.h>
-#include <AliL3Transform.h>
-#include <AliL3MemHandler.h>
-#include <AliL3TrackArray.h>
-#include <AliL3HoughMaxFinder.h>
-#include <AliL3HoughBaseTransformer.h>
-#include <AliL3Hough.h>
-#include <AliL3Benchmark.h>
+#include <AliHLTStandardIncludes.h>
+#include <AliHLTTransform.h>
+#include <AliHLTMemHandler.h>
+#include <AliHLTTrackArray.h>
+#include <AliHLTHoughMaxFinder.h>
+#include <AliHLTHoughBaseTransformer.h>
+#include <AliHLTHough.h>
+#include <AliHLTBenchmark.h>
+#include "AliESDVertex.h"
 #include <AliKalmanTrack.h>
+#include "AliITSgeomTGeo.h"
 #include "AliITSgeom.h"
 #include "AliMagF.h"
 #include "AliMagFMaps.h"
-#include <AliL3ITSclusterer.h>
-#include <AliL3ITSVertexerZ.h>
-#include <AliL3ITStracker.h>
+#include <AliHLTITSclusterer.h>
+#include <AliHLTITSVertexerZ.h>
+#include <AliHLTITStracker.h>
 #endif
 
-
 //_____________________________________________________________________________
 class AliGDCInterruptHandler : public TSignalHandler {
 public:
@@ -71,7 +72,7 @@ AliGDCInterruptHandler::AliGDCInterruptHandler() :
 
 
 //_____________________________________________________________________________
-#ifdef DATE_SYS
+#ifdef ALI_DATE
 int main(int argc, char** argv)
 {
   // set ROOT in batch mode
@@ -101,18 +102,36 @@ int main(int argc, char** argv)
   if (status) ::Fatal("monitorDeclareMp", monitorDecodeError(status));
 
   // initialize HLT transformations
-  if (!AliL3Transform::Init("./", kFALSE)) {
-    ::Fatal("AliL3Transform::Init", "HLT initialization failed");
+  if (!AliHLTTransform::Init("./", kFALSE)) {
+    ::Fatal("AliHLTTransform::Init", "HLT initialization failed");
   }
-  AliESD *esd = new AliESD;
-  AliKalmanTrack::SetConvConst(
-     1000/0.299792458/AliL3Transform::GetSolenoidField()
-  );
+  AliESDEvent *esd = new AliESDEvent;
+  esd->CreateStdContent();
+  //  AliKalmanTrack::SetConvConst(
+  //     1000/0.299792458/AliHLTTransform::GetSolenoidField()
+  //  );
   AliITSgeom *geom = new AliITSgeom();
   geom->ReadNewFile("$ALICE_ROOT/ITS/ITSgeometry_vPPRasymmFMD.det");
   if (!geom) return 1;
-  AliMagF* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., AliMagFMaps::k5kG);
-  AliTracker::SetFieldMap(field);
+  Int_t sfield = 0;
+  switch ((Int_t)(AliHLTTransform::GetSolenoidField()+0.5)) {
+  case 2:
+    sfield = AliMagFMaps::k2kG;
+    break;
+  case 4:
+    sfield = AliMagFMaps::k4kG;
+    break;
+  case 5:
+    sfield = AliMagFMaps::k5kG;
+    break;
+  default:
+    ::Fatal("AliHLTTransform::GetSolenoidField", "Incorrect magnetic field");
+  }
+  AliMagF* field = new AliMagFMaps("Maps","Maps", 2, 1., 10., sfield);
+  AliTracker::SetFieldMap(field,kTRUE);
+
+  // Init PID
+  AliPID pid;
 
   // create the signal handler
   AliGDCInterruptHandler* handler = new AliGDCInterruptHandler;
@@ -148,18 +167,18 @@ int main(int argc, char** argv)
 
       } else {
 
-       AliL3Benchmark *fBenchmark = new AliL3Benchmark();
+       AliHLTBenchmark *fBenchmark = new AliHLTBenchmark();
        fBenchmark->Start("Overall timing");
 
        // ITS clusterer and vertexer
        fBenchmark->Start("ITS Clusterer");
-       AliL3ITSclusterer clusterer(geom);
+       AliHLTITSclusterer clusterer(0);
        AliRawReader *itsrawreader=new AliRawReaderDate(ptr);
        TTree* treeClusters = new TTree("TreeL3ITSclusters"," "); //make a tree
        clusterer.Digits2Clusters(itsrawreader,treeClusters);
        fBenchmark->Stop("ITS Clusterer");
        
-       AliL3ITSVertexerZ vertexer;
+       AliHLTITSVertexerZ vertexer;
        AliESDVertex *vertex = vertexer.FindVertexForCurrentEvent(geom,treeClusters);
        Double_t vtxPos[3];
        Double_t vtxErr[3]={0.005,0.005,0.010};
@@ -168,12 +187,12 @@ int main(int argc, char** argv)
        esd->SetVertex(vertex);
 
        // TPC Hough reconstruction
-       Float_t ptmin = 0.1*AliL3Transform::GetSolenoidField();
+       Float_t ptmin = 0.1*AliHLTTransform::GetSolenoidField();
        Float_t zvertex = vtxPos[2];
 
        // Run the Hough Transformer
        fBenchmark->Start("Init");
-       AliL3Hough *hough1 = new AliL3Hough();
+       AliHLTHough *hough1 = new AliHLTHough();
 
        hough1->SetThreshold(4);
        hough1->CalcTransformerParams(ptmin);
@@ -184,7 +203,7 @@ int main(int argc, char** argv)
        fBenchmark->Stop("Init");
 
        fBenchmark->Start("Init");
-       AliL3Hough *hough2 = new AliL3Hough();
+       AliHLTHough *hough2 = new AliHLTHough();
 
        hough2->SetThreshold(4);
        hough2->CalcTransformerParams(ptmin);
@@ -201,9 +220,9 @@ int main(int argc, char** argv)
 
        //      gSystem->Sleep(20000);
        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! ");
 
        gSystem->MakeDirectory("hough1");
        hough1->WriteTracks("./hough1");
@@ -226,7 +245,7 @@ int main(int argc, char** argv)
          hough1->AddTracks();
          fBenchmark->Stop("AddTracks");
 
-         //      AliL3TrackArray* tracks = (AliL3TrackArray*)hough1->GetTracks(0);
+         //      AliHLTTrackArray* tracks = (AliHLTTrackArray*)hough1->GetTracks(0);
          //      nglobaltracks += tracks->GetNTracks();
        }
        for(int slice=18; slice<=35; slice++)
@@ -244,7 +263,7 @@ int main(int argc, char** argv)
          hough2->AddTracks();
          fBenchmark->Stop("AddTracks");
 
-         //      AliL3TrackArray* tracks = (AliL3TrackArray*)hough2->GetTracks(0);
+         //      AliHLTTrackArray* tracks = (AliHLTTrackArray*)hough2->GetTracks(0);
          //      nglobaltracks += tracks->GetNTracks();
        }
 
@@ -252,7 +271,7 @@ int main(int argc, char** argv)
        nglobaltracks += hough2->FillESD(esd);
 
        // ITS tracker
-       AliL3ITStracker itsTracker(geom);
+       AliHLTITStracker itsTracker(0);
        itsTracker.SetVertex(vtxPos,vtxErr);
 
        itsTracker.LoadClusters(treeClusters);