]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - MONITOR/AliMonitorHLTHough.cxx
new cuts and tasks added
[u/mrichter/AliRoot.git] / MONITOR / AliMonitorHLTHough.cxx
index f8e30a547fc8f67fd9fc0d5fb9162edf0412ddeb..4f0d8fe249a270762e059d5768a2f2c380d538a0 100644 (file)
 #include "AliMonitorHisto.h"
 #include "AliMonitorTrend.h"
 #include "AliTPCParam.h"
+#include "AliLog.h"
 #include <TFolder.h>
 #include <stdlib.h>
-#include <AliL3MemHandler.h>
-#include <AliL3TrackArray.h>
-#include <AliL3SpacePointData.h>
-#include <AliL3HoughTrack.h>
-#include <AliL3Transform.h>
+#include <AliHLTMemHandler.h>
+#include <AliHLTTrackArray.h>
+#include <AliHLTSpacePointData.h>
+#include <AliHLTHoughTrack.h>
+#include <AliHLTTransform.h>
 
 //_____________________________________________________________________________
-AliMonitorHLTHough::AliMonitorHLTHough(AliTPCParam* param)
+AliMonitorHLTHough::AliMonitorHLTHough(AliTPCParam* param):
+  AliMonitor(),
+  fParam(param),
+  fClustersCharge(NULL),
+  fNClustersVsRow(NULL),
+  fNClustersVsSector(NULL),
+  fNTracks(NULL),
+  fTrackPt(NULL),
+  fTrackEta(NULL),
+  fTrackPhi(NULL),
+  fTrackNHits(NULL),
+  fTrackDEdxVsP(NULL),
+  fTrackDEdx(NULL),
+  fTrackEtaVsPhi(NULL),
+  fPtEtaVsPhi(NULL)
 {
 // create a HLT monitor object with the given parameters
 
-  fParam = param;
 }
 
-//_____________________________________________________________________________
-AliMonitorHLTHough::AliMonitorHLTHough(const AliMonitorHLTHough& monitor) :
-  AliMonitor(monitor)
-{
-  Fatal("AliMonitorHLTHough", "copy constructor not implemented");
-}
-
-//_____________________________________________________________________________
-AliMonitorHLTHough& AliMonitorHLTHough::operator = (const AliMonitorHLTHough& 
-                                                   /*monitor*/)
-{
-  Fatal("operator =", "assignment operator not implemented");
-  return *this;
-}
-
-
 //_____________________________________________________________________________
 void AliMonitorHLTHough::CreateHistos(TFolder* folder)
 {
@@ -128,26 +126,26 @@ void AliMonitorHLTHough::CreateHistos(TFolder* folder)
 //_____________________________________________________________________________
 void AliMonitorHLTHough::FillHistos(AliRunLoader* /*runLoader*/, 
                                    AliRawReader* /*rawReader*/, 
-                                   AliESD* /*esd*/)
+                                   AliESDEvent* /*esd*/)
 {
 // fill the HLT Hough transform monitor histograms
 
-  AliL3MemHandler clusterHandler[36][6];
-  AliL3SpacePointData *clusters[36][6];
+  AliHLTMemHandler clusterHandler[36][6];
+  AliHLTSpacePointData *clusters[36][6];
   for (Int_t iSector = 0; iSector < fParam->GetNInnerSector(); iSector++) {
     for (Int_t iPatch = 0; iPatch < 6; iPatch++) {
       char fileName[256];
       sprintf(fileName, "hlt/fitter/points_%d_%d.raw", iSector,iPatch);
       if (!clusterHandler[iSector][iPatch].SetBinaryInput(fileName)) {
-       Warning("FillHistos", "could not open file %s", fileName);
+       AliWarning(Form("could not open file %s", fileName));
        continue;
       }
-      clusters[iSector][iPatch] = (AliL3SpacePointData*) clusterHandler[iSector][iPatch].Allocate();
+      clusters[iSector][iPatch] = (AliHLTSpacePointData*) clusterHandler[iSector][iPatch].Allocate();
       UInt_t nClusters = 0;
       clusterHandler[iSector][iPatch].Binary2Memory(nClusters, clusters[iSector][iPatch]);
 
       for (UInt_t iCluster = 0; iCluster < nClusters; iCluster++) {
-       AliL3SpacePointData& cluster = clusters[iSector][iPatch][iCluster];
+       AliHLTSpacePointData& cluster = clusters[iSector][iPatch][iCluster];
        fClustersCharge->Fill(cluster.fCharge);
        fNClustersVsRow->Fill(cluster.fPadRow);
        fNClustersVsSector->Fill(iSector);
@@ -160,21 +158,21 @@ void AliMonitorHLTHough::FillHistos(AliRunLoader* /*runLoader*/,
   fNClustersVsSector->ScaleErrorBy(10.);
 
 
-  AliL3MemHandler memHandler;
+  AliHLTMemHandler memHandler;
   Int_t nHoughTracks = 0;
 
   char fileName[256];
   sprintf(fileName, "hlt/fitter/tracks.raw");
   if (!memHandler.SetBinaryInput(fileName)) {
-    Warning("FillHistos", "could not open file hlt/fitter/tracks.raw");
+    AliWarning("could not open file hlt/fitter/tracks.raw");
     return;
   }
-  AliL3TrackArray* tracks = new AliL3TrackArray;
+  AliHLTTrackArray* tracks = new AliHLTTrackArray;
   memHandler.Binary2TrackArray(tracks);
 
   nHoughTracks += tracks->GetNTracks();
   for (Int_t iTrack = 0; iTrack < tracks->GetNTracks(); iTrack++) {
-    AliL3HoughTrack *track = (AliL3HoughTrack*)tracks->GetCheckedTrack(iTrack);
+    AliHLTHoughTrack *track = (AliHLTHoughTrack*)tracks->GetCheckedTrack(iTrack);
     if(!track) continue;
 
     track->CalculateHelix();
@@ -199,9 +197,9 @@ void AliMonitorHLTHough::FillHistos(AliRunLoader* /*runLoader*/,
       Int_t iPatch = (hitID>>22) & 0x7;
       UInt_t position = hitID&0x3fffff;
       UChar_t padrow = clusters[iSector][iPatch][position].fPadRow;
-      Float_t pWidth = AliL3Transform::GetPadPitchWidthLow();
+      Float_t pWidth = AliHLTTransform::GetPadPitchWidthLow();
       if (padrow>63)
-       pWidth = AliL3Transform::GetPadPitchWidthUp(); 
+       pWidth = AliHLTTransform::GetPadPitchWidthUp(); 
       Float_t corr=1.; if (padrow>63) corr=0.67;
       sampleDEdx[iHit] = clusters[iSector][iPatch][position].fCharge/pWidth*corr;
       Double_t crossingangle = track->GetCrossingAngle(padrow,iSector);