STEER/AliTracker.cxx

   1 /**************************************************************************
   2  * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
   3  *                                                                        *
   4  * Author: The ALICE Off-line Project.                                    *
   5  * Contributors are mentioned in the code where appropriate.              *
   6  *                                                                        *
   7  * Permission to use, copy, modify and distribute this software and its   *
   8  * documentation strictly for non-commercial purposes is hereby granted   *
   9  * without fee, provided that the above copyright notice appears in all   *
  10  * copies and that both the copyright notice and this permission notice   *
  11  * appear in the supporting documentation. The authors make no claims     *
  12  * about the suitability of this software for any purpose. It is          *
  13  * provided "as is" without express or implied warranty.                  *
  14  **************************************************************************/
  15
  16 /* $Id$ */
  17
  18 //-------------------------------------------------------------------------
  19 //               Implementation of the AliTracker class
  20 //  that is the base for AliTPCtracker, AliITStrackerV2 and AliTRDtracker
  21 //        Origin: Iouri Belikov, CERN, Jouri.Belikov@cern.ch
  22 //-------------------------------------------------------------------------
  23
  24 #include <TMath.h>
  25 #include <TH1F.h>
  26 #include <TGeoMatrix.h>
  27
  28 #include "AliTracker.h"
  29 #include "AliGeomManager.h"
  30 #include "AliCluster.h"
  31 #include "AliKalmanTrack.h"
  32 #include "AliGlobalQADataMaker.h"
  33
  34 Bool_t AliTracker::fFillResiduals=kFALSE;
  35 TObjArray **AliTracker::fResiduals=NULL;
  36 AliRecoParam::EventSpecie_t AliTracker::fEventSpecie=AliRecoParam::kDefault;
  37
  38 ClassImp(AliTracker)
  39
  40 AliTracker::AliTracker():
  41   AliTrackerBase(),
  42   fEventInfo(NULL)
  43 {
  44   //--------------------------------------------------------------------
  45   // The default constructor.
  46   //--------------------------------------------------------------------
  47 }
  48
  49 //__________________________________________________________________________
  50 AliTracker::AliTracker(const AliTracker &atr):
  51   AliTrackerBase(atr),
  52   fEventInfo(atr.fEventInfo)
  53 {
  54   //--------------------------------------------------------------------
  55   // The default constructor.
  56   //--------------------------------------------------------------------
  57 }
  58
  59 //__________________________________________________________________________
  60 void AliTracker::FillClusterArray(TObjArray* /*array*/) const
  61 {
  62   // Publishes all pointers to clusters known to the tracker into the
  63   // passed object array.
  64   // The ownership is not transfered - the caller is not expected to delete
  65   // the clusters.
  66
  67   AliWarning("should be overriden by a sub-class.");
  68 }
  69
  70 //__________________________________________________________________________
  71 void AliTracker::CookLabel(AliKalmanTrack *t, Float_t wrong) const {
  72   //--------------------------------------------------------------------
  73   //This function "cooks" a track label. If label<0, this track is fake.
  74   //--------------------------------------------------------------------
  75   Int_t noc=t->GetNumberOfClusters();
  76   if (noc<1) return;
  77   Int_t *lb=new Int_t[noc];
  78   Int_t *mx=new Int_t[noc];
  79   AliCluster **clusters=new AliCluster*[noc];
  80
  81   Int_t i;
  82   for (i=0; i<noc; i++) {
  83      lb[i]=mx[i]=0;
  84      Int_t index=t->GetClusterIndex(i);
  85      clusters[i]=GetCluster(index);
  86   }
  87
  88   Int_t lab=123456789;
  89   for (i=0; i<noc; i++) {
  90     AliCluster *c=clusters[i];
  91     lab=TMath::Abs(c->GetLabel(0));
  92     Int_t j;
  93     for (j=0; j<noc; j++) if (lb[j]==lab || mx[j]==0) break;
  94     lb[j]=lab;
  95     (mx[j])++;
  96   }
  97
  98   Int_t max=0;
  99   for (i=0; i<noc; i++) if (mx[i]>max) {max=mx[i]; lab=lb[i];}
 100
 101   for (i=0; i<noc; i++) {
 102     AliCluster *c=clusters[i];
 103     //if (TMath::Abs(c->GetLabel(1)) == lab ||
 104     //    TMath::Abs(c->GetLabel(2)) == lab ) max++;
 105     if (TMath::Abs(c->GetLabel(0)!=lab))
 106         if (TMath::Abs(c->GetLabel(1)) == lab ||
 107             TMath::Abs(c->GetLabel(2)) == lab ) max++;
 108   }
 109
 110   if ((1.- Float_t(max)/noc) > wrong) lab=-lab;
 111   t->SetFakeRatio((1.- Float_t(max)/noc));
 112   t->SetLabel(lab);
 113
 114   delete[] lb;
 115   delete[] mx;
 116   delete[] clusters;
 117 }
 118
 119 //____________________________________________________________________________
 120 void AliTracker::UseClusters(const AliKalmanTrack *t, Int_t from) const {
 121   //------------------------------------------------------------------
 122   //This function marks clusters associated with the track.
 123   //------------------------------------------------------------------
 124   Int_t noc=t->GetNumberOfClusters();
 125   for (Int_t i=from; i<noc; i++) {
 126      Int_t index=t->GetClusterIndex(i);
 127      AliCluster *c=GetCluster(index);
 128      c->Use();
 129   }
 130 }
 131
 132 void AliTracker::FillResiduals(const AliExternalTrackParam *t,
 133                               Double_t *p, Double_t *cov,
 134                               UShort_t id, Bool_t updated) {
 135   //
 136   // This function fills the histograms of residuals
 137   // The array of these histos is external for this AliTracker class.
 138   // Normally, this array belong to AliGlobalQADataMaker class.
 139   //
 140   if (!fFillResiduals) return;
 141   if (!fResiduals) return;
 142
 143   const Double_t *residuals=t->GetResiduals(p,cov,updated);
 144   if (!residuals) return;
 145
 146   TH1F *h=0;
 147   Int_t esIndex = AliRecoParam::AConvert(fEventSpecie) ;
 148   AliGeomManager::ELayerID layer=AliGeomManager::VolUIDToLayer(id);
 149   h=(TH1F*)fResiduals[esIndex]->At(2*layer-2);
 150   if (h) h->Fill(residuals[0]);
 151   h=(TH1F*)fResiduals[esIndex]->At(2*layer-1);
 152   if (h) h->Fill(residuals[1]);
 153
 154   if (layer==5) {
 155     if (p[1]<0) {  // SSD1 absolute residuals
 156       h = (TH1F*)fResiduals[esIndex]->At(40);
 157       if (h) h->Fill(t->GetY()-p[0]); //C side
 158       h = (TH1F*)fResiduals[esIndex]->At(41);
 159       if (h) h->Fill(t->GetZ()-p[1]);
 160     } else {
 161       h = (TH1F*)fResiduals[esIndex]->At(42);
 162       if (h) h->Fill(t->GetY()-p[0]); //A side
 163       h = (TH1F*)fResiduals[esIndex]->At(43);
 164       if (h) h->Fill(t->GetZ()-p[1]);
 165     }
 166   }
 167   if (layer==6) {  // SSD2 absolute residuals
 168     if (p[1]<0) {
 169       h = (TH1F*)fResiduals[esIndex]->At(44);
 170       if (h) h->Fill(t->GetY()-p[0]); //C side
 171       h = (TH1F*)fResiduals[esIndex]->At(45);
 172       if (h) h->Fill(t->GetZ()-p[1]);
 173     } else {
 174       h = (TH1F*)fResiduals[esIndex]->At(46);
 175       if (h) h->Fill(t->GetY()-p[0]); //A side
 176       h = (TH1F*)fResiduals[esIndex]->At(47);
 177       if (h) h->Fill(t->GetZ()-p[1]);
 178     }
 179   }
 180
 181 }
 182
 183 void AliTracker::FillResiduals(const AliExternalTrackParam *t,
 184                                const AliCluster *c, Bool_t /*updated*/) {
 185   //
 186   // This function fills the histograms of residuals
 187   // The array of these histos is external for this AliTracker class.
 188   // Normally, this array belong to AliGlobalQADataMaker class.
 189   //
 190   // For the moment, the residuals are absolute !
 191   //
 192
 193   if (!fFillResiduals) return;
 194   if (!fResiduals) return;
 195
 196   UShort_t id=c->GetVolumeId();
 197   const TGeoHMatrix *matrixT2L=AliGeomManager::GetTracking2LocalMatrix(id);
 198
 199   // Position of the cluster in the tracking c.s.
 200   Double_t clsTrk[3]={c->GetX(), c->GetY(), c->GetZ()};
 201   // Position of the cluster in the local module c.s.
 202   Double_t clsLoc[3]={0.,0.,0.};
 203   matrixT2L->LocalToMaster(clsTrk,clsLoc);
 204
 205
 206   // Position of the intersection point in the tracking c.s.
 207   Double_t trkTrk[3]={t->GetX(),t->GetY(),t->GetZ()};
 208   // Position of the intersection point in the local module c.s.
 209   Double_t trkLoc[3]={0.,0.,0.};
 210   matrixT2L->LocalToMaster(trkTrk,trkLoc);
 211
 212   Double_t residuals[2]={trkLoc[0]-clsLoc[0], trkLoc[2]-clsLoc[2]};
 213
 214   TH1F *h=0;
 215   Int_t esIndex = AliRecoParam::AConvert(fEventSpecie) ;
 216   AliGeomManager::ELayerID layer=AliGeomManager::VolUIDToLayer(id);
 217   h=(TH1F*)fResiduals[esIndex]->At(2*layer-2);
 218   if (h) h->Fill(residuals[0]);
 219   h=(TH1F*)fResiduals[esIndex]->At(2*layer-1);
 220   if (h) h->Fill(residuals[1]);
 221
 222 }
 223