- MatchTracks();
-
- AliInfo(Form("Number of matched tracks = %d",fnmatch));
-
- for (Int_t i=0; i<ntrk; i++) {
- AliESDtrack *t=event->GetTrack(i);
- //t->SetESDEvent(event);
- AliESDtrack *seed =(AliESDtrack*)fSeeds->At(i);
-
- if ( (seed->GetStatus()&AliESDtrack::kTOFin)!=0 ) {
- t->SetStatus(AliESDtrack::kTOFin);
- if ( (seed->GetStatus()&AliESDtrack::kTOFout)!=0 ) {
- t->SetStatus(AliESDtrack::kTOFout);
- //t->SetTOFclusterArray(seed->GetNTOFclusters(),seed->GetTOFclusterArray());
- t->SortTOFcluster();
-
- // Make attention, please:
- // AliESDtrack::fTOFInfo array does not be stored in the AliESDs.root file
- // it is there only for a check during the reconstruction step.
- Float_t info[10]; seed->GetTOFInfo(info);
- t->SetTOFInfo(info);
- AliDebug(3,Form(" distance=%f; residual in the pad reference frame: dX=%f, dZ=%f", info[0],info[1],info[2]));
-
- /*
- Double_t alphaA = (Double_t)t->GetAlpha();
- Double_t xA = (Double_t)t->GetX();
- Double_t yA = (Double_t)t->GetY();
- Double_t zA = (Double_t)t->GetZ();
- Double_t p1A = (Double_t)t->GetSnp();
- Double_t p2A = (Double_t)t->GetTgl();
- Double_t p3A = (Double_t)t->GetSigned1Pt();
- const Double_t *covA = (Double_t*)t->GetCovariance();
-
- // Check done:
- // by calling the AliESDtrack::UpdateTrackParams,
- // the current track parameters are changed
- // and it could cause refit problems.
- // We need to update only the following track parameters:
- // the track length and expected times.
- // Removed AliESDtrack::UpdateTrackParams call
- // Called AliESDtrack::SetIntegratedTimes(...) and
- // AliESDtrack::SetIntegratedLength() routines.
-
- AliTOFtrack *track = new AliTOFtrack(*seed);
- t->UpdateTrackParams(track,AliESDtrack::kTOFout); // to be checked - AdC
- delete track;
- Double_t time[AliPID::kSPECIESC]; t->GetIntegratedTimes(time);
- Double_t alphaB = (Double_t)t->GetAlpha();
- Double_t xB = (Double_t)t->GetX();
- Double_t yB = (Double_t)t->GetY();
- Double_t zB = (Double_t)t->GetZ();
- Double_t p1B = (Double_t)t->GetSnp();
- Double_t p2B = (Double_t)t->GetTgl();
- Double_t p3B = (Double_t)t->GetSigned1Pt();
- const Double_t *covB = (Double_t*)t->GetCovariance();
- AliDebug(2,"Track params -now(before)-:");
- AliDebug(2,Form(" X: %f(%f), Y: %f(%f), Z: %f(%f) --- alpha: %f(%f)",
- xB,xA,
- yB,yA,
- zB,zA,
- alphaB,alphaA));
- AliDebug(2,Form(" p1: %f(%f), p2: %f(%f), p3: %f(%f)",
- p1B,p1A,
- p2B,p2A,
- p3B,p3A));
- AliDebug(2,Form(" cov1: %f(%f), cov2: %f(%f), cov3: %f(%f)"
- " cov4: %f(%f), cov5: %f(%f), cov6: %f(%f)"
- " cov7: %f(%f), cov8: %f(%f), cov9: %f(%f)"
- " cov10: %f(%f), cov11: %f(%f), cov12: %f(%f)"
- " cov13: %f(%f), cov14: %f(%f), cov15: %f(%f)",
- covB[0],covA[0],
- covB[1],covA[1],
- covB[2],covA[2],
- covB[3],covA[3],
- covB[4],covA[4],
- covB[5],covA[5],
- covB[6],covA[6],
- covB[7],covA[7],
- covB[8],covA[8],
- covB[9],covA[9],
- covB[10],covA[10],
- covB[11],covA[11],
- covB[12],covA[12],
- covB[13],covA[13],
- covB[14],covA[14]
- ));
- */
- Double_t time[AliPID::kSPECIESC]; t->GetIntegratedTimes(time,AliPID::kSPECIESC);
- AliDebug(2,Form(" TOF params: %6d %f %f %f %f %f %6d %3d %f",
- i,
- t->GetTOFsignalRaw(),t->GetTOFsignal(),t->GetTOFsignalToT(),
- t->GetTOFsignalDz(),t->GetTOFsignalDx(),t->GetTOFCalChannel(),
- t->GetTOFcluster(),t->GetIntegratedLength()));
- AliDebug(2,Form(" %f %f %f %f %f %f %f %f %f",
- time[0], time[1], time[2], time[3], time[4], time[5], time[6], time[7], time[8]));
- }
- }
- }
-
- fSeeds->Clear();
- fTracks->Clear();