--- PamCut/TrkCuts/TrkDedxHeCut/TrkDedxHeCut.cpp 2010/05/03 14:36:27 1.2 +++ PamCut/TrkCuts/TrkDedxHeCut/TrkDedxHeCut.cpp 2015/03/25 17:38:08 1.3 @@ -11,7 +11,11 @@ int TrkDedxHeCut::Check(PamLevel2 *event) { - TrkTrack *track = event->GetTrack(0)->GetTrkTrack(); +// TrkTrack *track = event->GetTrack(0)->GetTrkTrack(); + if(event->GetNTracks(_trkAlg)==0)return 0; + ExtTrack *track = event->GetTrack(0,_trkAlg)->GetExtTrack(); + + float R = track->GetRigidity(); //float dEdx = track->GetDEDX(); @@ -21,29 +25,6 @@ float totReleaseY = 0., maxReleaseY = 0.; int nX = 0, nY = 0; float dEdxView = 0; - /* - ToFLevel2 *tofL2 = event->GetToFLevel2(); - int trkSeqNo = 0; - for (int i = 0; i < event->GetToFLevel2()->ntrk(); i++) { - if (tofL2->GetToFTrkVar(i)->trkseqno == track->GetSeqNo()) { - trkSeqNo = i; - break; - } - } - if ((tofL2->GetdEdx(trkSeqNo, 0, 100) + tofL2->GetdEdx(trkSeqNo, 1, 100) + tofL2->GetdEdx(trkSeqNo, 2, 100) - + tofL2->GetdEdx(trkSeqNo, 3, 100)) / 4. < 2) { - cout << "Check" << endl; - for (int ip = 0; ip < 6; ip++) { - cout << "TOF" << ip << ": " << tofL2->GetdEdx(trkSeqNo, ip, 100) << endl; - dEdxView = track->GetDEDX(ip, 0); - if (dEdxView > 0 && track->XGood(ip)) - cout << "X" << ip << ": " << dEdxView << endl; - dEdxView = track->GetDEDX(ip, 1); - if (dEdxView > 0 && track->YGood(ip)) - cout << "Y" << ip << ": " << dEdxView << endl; - } - } - */ for (int ip = 0; ip < 6; ip++) { // X view dEdxView = track->GetDEDX(ip, 0);