temporary workaround for circular dependency libSTEERbase <-> libRAWDatabase, see...

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

diff --git a/STEER/AliTracker.cxx b/STEER/AliTracker.cxx
index f9bbbb708516d6721d887279799babeb516b8af6..18d65ebc013fe4ae8b33d3717d29fe1c0e127847 100644 (file)
--- a/STEER/AliTracker.cxx
+++ b/STEER/AliTracker.cxx
@@ -24,6 +24,7 @@
 #include <TMath.h>
 #include <TH1F.h>
 #include <TGeoManager.h>
+#include <TGeoMatrix.h>
 
 #include "AliMagF.h"
 #include "AliTracker.h"
@@ -37,7 +38,8 @@ Bool_t AliTracker::fgUniformField=kTRUE;
 Double_t AliTracker::fgBz=kAlmost0Field;
 const AliMagF *AliTracker::fgkFieldMap=0;
 Bool_t AliTracker::fFillResiduals=kFALSE;
-TObjArray *AliTracker::fResiduals=0;
+TObjArray **AliTracker::fResiduals=NULL;
+AliRecoParam::EventSpecie_t AliTracker::fEventSpecie=AliRecoParam::kDefault;
 
 ClassImp(AliTracker)
 
@@ -48,7 +50,7 @@ AliTracker::AliTracker():
   fZ(0),
   fSigmaX(0.005),
   fSigmaY(0.005),
-  fSigmaZ(0.010)
+  fSigmaZ(0.010) 
 {
   //--------------------------------------------------------------------
   // The default constructor.
@@ -95,6 +97,17 @@ 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 {
   //--------------------------------------------------------------------
@@ -157,7 +170,7 @@ void AliTracker::UseClusters(const AliKalmanTrack *t, Int_t from) const {
   }
 }
 
-Double_t AliTracker::GetBz(Float_t *r) {
+Double_t AliTracker::GetBz(const Float_t *r) {
   //------------------------------------------------------------------
   // Returns Bz (kG) at the point "r" .
   //------------------------------------------------------------------
@@ -167,7 +180,7 @@ Double_t AliTracker::GetBz(Float_t *r) {
 }
 
 Double_t 
-AliTracker::MeanMaterialBudget(Double_t *start,Double_t *end,Double_t *mparam)
+AliTracker::MeanMaterialBudget(const Double_t *start, const Double_t *end, Double_t *mparam)
 {
   // 
   // Calculate mean material budget and material properties between 
@@ -388,9 +401,68 @@ void AliTracker::FillResiduals(const AliExternalTrackParam *t,
 
   TH1F *h=0;
   AliGeomManager::ELayerID layer=AliGeomManager::VolUIDToLayer(id);
-  h=(TH1F*)fResiduals->At(2*layer-2);
+  h=(TH1F*)fResiduals[fEventSpecie]->At(2*layer-2);
   h->Fill(residuals[0]);
-  h=(TH1F*)fResiduals->At(2*layer-1);
+  h=(TH1F*)fResiduals[fEventSpecie]->At(2*layer-1);
   h->Fill(residuals[1]);
+
+  if (layer==5) {
+    if (p[1]<0) {  // SSD1 absolute residuals
+       ((TH1F*)fResiduals[fEventSpecie]->At(40))->Fill(t->GetY()-p[0]); //C side
+       ((TH1F*)fResiduals[fEventSpecie]->At(41))->Fill(t->GetZ()-p[1]);
+    } else {             
+       ((TH1F*)fResiduals[fEventSpecie]->At(42))->Fill(t->GetY()-p[0]); //A side
+       ((TH1F*)fResiduals[fEventSpecie]->At(43))->Fill(t->GetZ()-p[1]);
+    }           
+  }
+  if (layer==6) {  // SSD2 absolute residuals
+    if (p[1]<0) {
+       ((TH1F*)fResiduals[fEventSpecie]->At(44))->Fill(t->GetY()-p[0]); //C side
+       ((TH1F*)fResiduals[fEventSpecie]->At(45))->Fill(t->GetZ()-p[1]);
+    } else {
+       ((TH1F*)fResiduals[fEventSpecie]->At(46))->Fill(t->GetY()-p[0]); //A side
+       ((TH1F*)fResiduals[fEventSpecie]->At(47))->Fill(t->GetZ()-p[1]);
+    }
+  }
+
 }
 
+void AliTracker::FillResiduals(const AliExternalTrackParam *t,
+                               const AliCluster *c, 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.  
+  // 
+  // For the moment, the residuals are absolute !
+  //
+
+  if (!fFillResiduals) return; 
+  if (!fResiduals) return; 
+
+  UShort_t id=c->GetVolumeId();
+  const TGeoHMatrix *matrixT2L=AliGeomManager::GetTracking2LocalMatrix(id);
+
+  // Position of the cluster in the tracking c.s.
+  Double_t clsTrk[3]={c->GetX(), c->GetY(), c->GetZ()};
+  // Position of the cluster in the local module c.s.
+  Double_t clsLoc[3]={0.,0.,0.};
+  matrixT2L->LocalToMaster(clsTrk,clsLoc);
+
+
+  // Position of the intersection point in the tracking c.s.
+  Double_t trkTrk[3]={t->GetX(),t->GetY(),t->GetZ()};
+  // Position of the intersection point in the local module c.s.
+  Double_t trkLoc[3]={0.,0.,0.};
+  matrixT2L->LocalToMaster(trkTrk,trkLoc);
+
+  Double_t residuals[2]={trkLoc[0]-clsLoc[0], trkLoc[2]-clsLoc[2]};
+
+  TH1F *h=0;
+  AliGeomManager::ELayerID layer=AliGeomManager::VolUIDToLayer(id);
+  h=(TH1F*)fResiduals[fEventSpecie]->At(2*layer-2);
+  h->Fill(residuals[0]);
+  h=(TH1F*)fResiduals[fEventSpecie]->At(2*layer-1);
+  h->Fill(residuals[1]);
+
+}