adding first version of common HLT track data struct for barrel tracks

[u/mrichter/AliRoot.git] / STEER / AliTracker.cxx

diff --git a/STEER/AliTracker.cxx b/STEER/AliTracker.cxx
index b096be58cfaec0d94b5f6f403f693bd601b0335e..89388ab9955a240c5dc5ff9f4cbc816decc7a3e6 100644 (file)
--- a/STEER/AliTracker.cxx
+++ b/STEER/AliTracker.cxx
@@ -22,10 +22,12 @@
 //-------------------------------------------------------------------------
 #include <TClass.h>
 #include <TMath.h>
+#include <TH1F.h>
 #include <TGeoManager.h>
 
 #include "AliMagF.h"
 #include "AliTracker.h"
+#include "AliGeomManager.h"
 #include "AliCluster.h"
 #include "AliKalmanTrack.h"
 
@@ -34,6 +36,8 @@ extern TGeoManager *gGeoManager;
 Bool_t AliTracker::fgUniformField=kTRUE;
 Double_t AliTracker::fgBz=kAlmost0Field;
 const AliMagF *AliTracker::fgkFieldMap=0;
+Bool_t AliTracker::fFillResiduals=kFALSE;
+TObjArray *AliTracker::fResiduals=0;
 
 ClassImp(AliTracker)
 
@@ -91,12 +95,24 @@ void AliTracker::SetFieldMap(const AliMagF* map, Bool_t uni) {
 
 }
 
+//__________________________________________________________________________
+void AliTracker::FillClusterArray(TObjArray* /*array*/) const
+{
+  // Publishes all pointers to clusters known to the tracker into the
+  // passed object array.
+  // The ownership is not transfered - the caller is not expected to delete
+  // the clusters.
+
+  AliWarning("should be overriden by a sub-class.");
+}
+
 //__________________________________________________________________________
 void AliTracker::CookLabel(AliKalmanTrack *t, Float_t wrong) const {
   //--------------------------------------------------------------------
   //This function "cooks" a track label. If label<0, this track is fake.
   //--------------------------------------------------------------------
   Int_t noc=t->GetNumberOfClusters();
+  if (noc<1) return;
   Int_t *lb=new Int_t[noc];
   Int_t *mx=new Int_t[noc];
   AliCluster **clusters=new AliCluster*[noc];
@@ -216,7 +232,8 @@ AliTracker::MeanMaterialBudget(Double_t *start,Double_t *end,Double_t *mparam)
   TGeoNode *startnode = gGeoManager->InitTrack(start, dir);
   //printf("%s length=%f\n",gGeoManager->GetPath(),length);
   if (!startnode) {
-    printf("ERROR: start point out of geometry\n");
+    AliErrorClass(Form("start point out of geometry: x %f, y %f, z %f",
+                 start[0],start[1],start[2]));
     return 0.0;
   }
   TGeoMaterial *material = startnode->GetVolume()->GetMedium()->GetMaterial();
@@ -366,3 +383,45 @@ Double_t mass, Double_t maxStep, Bool_t rotateTo, Double_t maxSnp){
   return kTRUE;
 }
 
+void AliTracker::FillResiduals(const AliExternalTrackParam *t,
+                             Double_t *p, Double_t *cov, 
+                              UShort_t id, Bool_t updated) {
+  //
+  // This function fills the histograms of residuals 
+  // The array of these histos is external for this AliTracker class.
+  // Normally, this array belong to AliGlobalQADataMaker class.  
+  // 
+  if (!fFillResiduals) return; 
+  if (!fResiduals) return; 
+
+  const Double_t *residuals=t->GetResiduals(p,cov,updated);
+  if (!residuals) return;
+
+  TH1F *h=0;
+  AliGeomManager::ELayerID layer=AliGeomManager::VolUIDToLayer(id);
+  h=(TH1F*)fResiduals->At(2*layer-2);
+  h->Fill(residuals[0]);
+  h=(TH1F*)fResiduals->At(2*layer-1);
+  h->Fill(residuals[1]);
+
+  if (layer==5) {
+    if (p[1]<0) {  // SSD1 absolute residuals
+       ((TH1F*)fResiduals->At(40))->Fill(t->GetY()-p[0]); //C side
+       ((TH1F*)fResiduals->At(41))->Fill(t->GetZ()-p[1]);
+    } else {             
+       ((TH1F*)fResiduals->At(42))->Fill(t->GetY()-p[0]); //A side
+       ((TH1F*)fResiduals->At(43))->Fill(t->GetZ()-p[1]);
+    }           
+  }
+  if (layer==6) {  // SSD2 absolute residuals
+    if (p[1]<0) {
+       ((TH1F*)fResiduals->At(44))->Fill(t->GetY()-p[0]); //C side
+       ((TH1F*)fResiduals->At(45))->Fill(t->GetZ()-p[1]);
+    } else {
+       ((TH1F*)fResiduals->At(46))->Fill(t->GetY()-p[0]); //A side
+       ((TH1F*)fResiduals->At(47))->Fill(t->GetZ()-p[1]);
+    }
+  }
+
+}
+