layersWithTracklet = 0;
UInt_t stackHits[fgkTrdStacks];
- Double_t x[3];
+ Double_t x[3] = { 0. };
memset(stackHits, 0, fgkTrdStacks*sizeof(UInt_t));
#ifdef TRD_TM_DEBUG
Double_t r = fgkSaveInnerRadius;
while (r < fgkSaveOuterRadius){
- track->GetXYZAt(r, magFieldinKiloGauss, x);
- stack = EstimateStack(x);
- if (stack >= 0){
- stackHits[stack]++;
- if (stackHits[stack] > 16) // experimental
- break;
+ if (track->GetXYZAt(r, magFieldinKiloGauss, x)) {
+ stack = EstimateStack(x);
+ if (stack >= 0){
+ stackHits[stack]++;
+ if (stackHits[stack] > 16) // experimental
+ break;
#ifdef TRD_TM_DEBUG
- printf(" r=%.3fcm %.2f/%.2f - %d hits for stack %d S%02d-%d (mag=%.1f)\n",
- r, x[0], x[1], stackHits[stack], stack, stack/5, stack%5, magFieldinKiloGauss);
+ printf(" r=%.3fcm %.2f/%.2f - %d hits for stack %d S%02d-%d (mag=%.1f)\n",
+ r, x[0], x[1], stackHits[stack], stack, stack/5, stack%5, magFieldinKiloGauss);
#endif
+ }
}
r += 1.;
}
(trdTrack->GetLabel() != label))
continue;
+ if ((trdTrack->GetSector() < 0) || (trdTrack->GetSector() > 17) ||
+ (trdTrack->GetStack() < 0) || (trdTrack->GetStack() > 4))
+ continue;
+
stack = TrdSecSiLsi(trdTrack->GetSector(), trdTrack->GetStack());
trdPt = (esdEvent->GetMagneticField() > 0.) ? (-1.*trdTrack->Pt()) : trdTrack->Pt();
matchTrack = NULL;
}
if ((matchTrack) && (matchRating >= fMinMatchRating)){
+ AliDebug(1, Form("S%02d-%d trd %d - esd %d match! pt: %.2f %.2f",
+ trdTrack->GetSector(), trdTrack->GetStack(), iTrdTrack, matchEsdTrackIndexInStack,
+ trdPt, matchTrack->GetSignedPt()));
#ifdef TRD_TM_DEBUG
printf("#TRACKMATCHING S%02d-%d trd %d - esd %d match! pt: %.2f %.2f\n",
trdTrack->GetSector(), trdTrack->GetStack(), iTrdTrack, matchEsdTrackIndexInStack,