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 |
} |
} |