TRD/qaRec/macros/clusterResiduals.C

   1 #if !defined(__CINT__) || defined(__MAKECINT__)
   2 #include <TRD/AliTRDgeometry.h>
   3 #include <TRD/AliTRDcluster.h>
   4 #include <TRD/AliTRDseedV1.h>
   5 #include <TRD/AliTRDtrackV1.h>
   6 #endif
   7
   8 void clusterResiduals(const AliTRDtrackV1* track, Double_t* &res, Int_t& n)
   9 {
  10   if (!track)  return;
  11
  12   Int_t in = 0;
  13   n = track->GetNumberOfClusters();
  14   res = new Double_t[n];
  15   memset(res, 0, n*sizeof(Double_t));
  16
  17   AliTRDseedV1 *fTracklet = 0x0;
  18   for(Int_t ily=0; ily<AliTRDgeometry::kNlayer; ily++){
  19     if(!(fTracklet = track->GetTracklet(ily))) continue;
  20     if(!fTracklet->IsOK()) continue;
  21     if(!fTracklet->Fit(kFALSE)) continue;
  22
  23     AliTRDcluster *c = 0x0;
  24     for(Int_t ic=AliTRDseed::knTimebins-1; ic>=0; ic--){
  25       if(!(c = fTracklet->GetClusters(ic))) continue;
  26
  27       res[in++] = fTracklet->GetYat(c->GetX()) - c->GetY();
  28     }
  29   }
  30   return;
  31 }