| 24 |
|
|
| 25 |
trkTrack = pamTrack->GetTrkTrack(); |
trkTrack = pamTrack->GetTrkTrack(); |
| 26 |
|
|
| 27 |
|
// Check if good clusters have meaningful energy releases |
| 28 |
|
for (unsigned int ip = 0; ip < 6; ip++) { |
| 29 |
|
if (trkTrack->XGood(ip) == 1 && trkTrack->GetDEDX(ip, 0) <= 0.) |
| 30 |
|
trkTrack->xgood[ip] = 0; // Set it to bad |
| 31 |
|
if (trkTrack->YGood(ip) == 1 && trkTrack->GetDEDX(ip, 1) <= 0.) |
| 32 |
|
trkTrack->ygood[ip] = 0; // Set it to bad |
| 33 |
|
} |
| 34 |
|
|
| 35 |
if (!(trkTrack->nstep < 100 && trkTrack->chi2 > 0.)) |
if (!(trkTrack->nstep < 100 && trkTrack->chi2 > 0.)) |
| 36 |
return 0; // check that the track is meaningful |
return 0; // check that the track is meaningful |
| 37 |
|
|
| 38 |
// int imatch=0; |
// identification of TRK-track TOF-index (this index is generally different from TRK-track TRK-index) |
|
|
|
|
// identification of TRK-track TOF-index (this index is generally different from TRK-track TRK-index) |
|
| 39 |
iTrackTof = -1000; |
iTrackTof = -1000; |
| 40 |
for (int i=0; i<event->GetToFLevel2()->ntrk(); i++) { |
for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) { |
| 41 |
if ( event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == trkTrack->GetSeqNo() ) { |
if (event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno == trkTrack->GetSeqNo()) { |
| 42 |
iTrackTof = i; |
iTrackTof = i; |
|
// cout << "trkseqno " << event->GetToFLevel2()->GetToFTrkVar(i)->trkseqno << " TOF track id " << i << endl; |
|
|
// imatch++; |
|
| 43 |
break; |
break; |
| 44 |
} |
} |
| 45 |
} |
} |