--- PamCut/TrkCuts/TrkGeomCut/TrkGeomCut.cpp 2009/05/27 13:30:09 1.1 +++ PamCut/TrkCuts/TrkGeomCut/TrkGeomCut.cpp 2010/03/10 08:54:20 1.2 @@ -2,7 +2,7 @@ * TrkGeomCut.cpp * * Created on: 13-mar-2009 - * Author: Sergio Ricciarini, Nicola Mori + * Author: N. Mori, S. Ricciarini */ /*! @file TrkGeomCut.cpp The TrkGeomCut class implementation file */ @@ -12,13 +12,24 @@ int TrkGeomCut::Check(PamLevel2 *event) { TrkTrack *trkTrack = event->GetTrack(0)->GetTrkTrack(); + + Float_t xT, yT; for (int igf = 0; igf < TrkParams::nGF; igf++) { // loop on GF planes + if (igf==5 || igf==10) { // M1 and M2 planes + xT = _xMagnTol; + yT = _yMagnTol; + } + else { + xT = _xTol; + yT = _yTol; + } + // check done for all 14 planes - if (!(TrkParams::xGF_min[igf] + _xTol < trkTrack->xGF[igf] && trkTrack->xGF[igf] < TrkParams::xGF_max[igf] - _xTol)) + if (!(TrkParams::xGF_min[igf] + xT < trkTrack->xGF[igf] && trkTrack->xGF[igf] < TrkParams::xGF_max[igf] - xT)) return 0; - if (!(TrkParams::yGF_min[igf] + _yTol < trkTrack->yGF[igf] && trkTrack->yGF[igf] < TrkParams::yGF_max[igf] - _yTol)) + if (!(TrkParams::yGF_min[igf] + yT < trkTrack->yGF[igf] && trkTrack->yGF[igf] < TrkParams::yGF_max[igf] - yT)) return 0; } // end loop on GF planes