//-------------------------------------------------------------------------
#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"
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)
}
+//__________________________________________________________________________
+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];
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]);
+ }
+ }
+
+}
+