fRunLoader->UnloadKinematics();
fRunLoader->UnloadTrackRefs();
fRunLoader->UnloadTracks();
+ fMuonTrackRef->Delete();
delete fMuonTrackRef;
delete fMUONData;
}
iHitMin = 0;
isNewTrack = kTRUE;
- if (!trackRefs->GetEntries()) continue;
-
- while (isNewTrack) {
+ if (!trackRefs->GetEntries()) continue;
+ while (isNewTrack) {
+
for (Int_t iHit = iHitMin; iHit < trackRefs->GetEntries(); iHit++) {
trackReference = (AliTrackReference*)trackRefs->At(iHit);
hitForRec->SetZ(z);
hitForRec->SetBendingReso2(0.0);
hitForRec->SetNonBendingReso2(0.0);
- iChamber = ChamberNumber(z);
+ iChamber = AliMUONConstants::ChamberNumber(z);
hitForRec->SetChamberNumber(iChamber);
muonTrack->AddTrackParamAtHit(trackParam);
delete muonTrack;
delete trackParam;
delete hitForRec;
+ trackRefs->Delete();
delete trackRefs;
}
{
// Return TClonesArray of reconstructed tracks
+ GetMUONData()->ResetRecTracks();
GetMUONData()->SetTreeAddress("RT");
fTrackReco = GetMUONData()->RecTracks();
GetMUONData()->GetRecTracks();
hitForRec->SetNonBendingCoor(xRec);
hitForRec->SetBendingCoor(yRec);
hitForRec->SetZ(zRec);
- iChamber = ChamberNumber(zRec);
+ iChamber = AliMUONConstants::ChamberNumber(zRec);
hitForRec->SetChamberNumber(iChamber);
hitForRec->SetBendingReso2(0.0);
hitForRec->SetNonBendingReso2(0.0);
delete trackNew;
delete hitForRec;
delete trackParam;
+ newMuonTrackRef->Delete();
delete newMuonTrackRef;
}
isTrackOK = kTRUE;
for (Int_t ch = 0; ch < 10; ch++) {
if (!isChamberInTrack[ch]) isTrackOK = kFALSE;
- }
- if (isTrackOK) fReconstructibleTracks++;
+ } if (isTrackOK) fReconstructibleTracks++;
if (!isTrackOK) fMuonTrackRef->Remove(track); // remove non reconstructible tracks
}
fMuonTrackRef->Compress();
}
-//_____________________________________________________________________________
-Int_t AliMUONRecoCheck::ChamberNumber(Float_t z) const
-{
- // return chamber number according z position of hit. Should be taken from geometry ?
-
- Float_t dMaxChamber = AliMUONConstants::DzSlat() + AliMUONConstants::DzCh() + 0.25; // cm st 3 &4 & 5
- if ( z > (AliMUONConstants::DefaultChamberZ(4)+50.)) dMaxChamber = 7.; // cm stations 1 & 2
- Int_t iChamber;
-
- for (iChamber = 0; iChamber < 10; iChamber++) {
-
- if (TMath::Abs(z-AliMUONConstants::DefaultChamberZ(iChamber)) < dMaxChamber) {
- return iChamber;
- }
- }
- return -1;
-}