+ //__________________________________________________________________________
+void AliMUONTrackK::Print(const char* /*opt*/) const
+{
+/// Print out track information
+
+ cout << 1/(*fTrackPar)(4,0) << " " << fChi2 << " " << fNmbTrackHits << " " << fNSteps << endl;
+
+ AliMUONHitForRec *hit;
+ for (Int_t i = 0; i < fNmbTrackHits; ++i) {
+ hit = (AliMUONHitForRec*) ((*fTrackHits)[i]);
+ printf ("%5d", hit->GetChamberNumber());
+ }
+ cout << endl;
+ if (fgDebug > 0) {
+ for (Int_t i = 0; i < fNmbTrackHits; ++i) {
+ hit = (AliMUONHitForRec*) ((*fTrackHits)[i]);
+ //cout << ((AliMUONHitForRec*)((*fTrackHits)[i1]))->GetHitNumber() << " ";
+ //cout << ((AliMUONHitForRec*)((*fTrackHits)[i1]))->GetZ() << " ";
+ printf ("%5d", fgHitForRec->IndexOf(hit));
+ }
+ cout << endl;
+ }
+
+ // from raw clusters
+ /*
+ TClonesArray *rawclusters = 0x0;
+ AliMUONRawCluster *clus = 0x0;
+ for (Int_t i = 0; i < fNmbTrackHits; ++i) {
+ hit = (AliMUONHitForRec*) ((*fTrackHits)[i]);
+ if (hit->GetHitNumber() < 0) { // combined cluster / track finder
+ Int_t index = -hit->GetHitNumber() / 100000;
+ Int_t iPos = -hit->GetHitNumber() - index * 100000;
+ //clus = (AliMUONRawCluster*) fgCombi->DetElem(index-1)->RawClusters()->UncheckedAt(iPos);
+ } else {
+ rawclusters = fgTrackReconstructor->GetMUONData()->RawClusters(hit->GetChamberNumber());
+ clus = (AliMUONRawCluster*) rawclusters->UncheckedAt(hit->GetHitNumber());
+ }
+ printf ("%5d", clus->GetTrack(1)%10000000);
+ }
+ cout << endl;
+ for (Int_t i = 0; i < fNmbTrackHits; ++i) {
+ hit = (AliMUONHitForRec*) ((*fTrackHits)[i]);
+ if (hit->GetHitNumber() < 0) { // combined cluster / track finder
+ Int_t index = -hit->GetHitNumber() / 100000;
+ Int_t iPos = -hit->GetHitNumber() - index * 100000;
+ clus = (AliMUONRawCluster*) fgCombi->DetElem(index-1)->RawClusters()->UncheckedAt(iPos);
+ } else {
+ rawclusters = fgTrackReconstructor->GetMUONData()->RawClusters(hit->GetChamberNumber());
+ clus = (AliMUONRawCluster*) rawclusters->UncheckedAt(hit->GetHitNumber());
+ }
+ if (clus->GetTrack(2) != -1) printf ("%5d", clus->GetTrack(2)%10000000);
+ else printf ("%5s", " ");
+ }
+ cout << endl;
+ */
+
+ for (Int_t i = 0; i < fNmbTrackHits; ++i) {
+ hit = (AliMUONHitForRec*) ((*fTrackHits)[i]);
+ printf ("%5d", hit->GetNTrackHits());
+ }
+ cout << endl;
+ /*
+ for (Int_t i = 0; i < fNmbTrackHits; ++i) {
+ hit = (AliMUONHitForRec*) ((*fTrackHits)[i]);
+ cout << hit << " ";
+ }
+ cout << endl;
+ */
+
+ for (Int_t i = 0; i < fNmbTrackHits; ++i) {
+ //cout << ((AliMUONHitForRec*)((*fTrackHits)[i1]))->GetHitNumber() << " ";
+ cout << ((AliMUONHitForRec*)((*fTrackHits)[i]))->GetZ() << " ";
+ //cout << fgHitForRec->IndexOf(((AliMUONHitForRec*)((*fTrackHits)[i1]))) << " ";
+ }
+ cout << endl;
+ for (Int_t i = 0; i < fNmbTrackHits; ++i) printf("%8.4f", (*fChi2Smooth)[i]);
+ cout << endl;
+ cout << "---------------------------------------------------" << endl;
+}
+