TrkCuts/TrkGeomCut/TrkGeomCut.cpp

/*
 * TrkGeomCut.cpp
 *
 *  Created on: 13-mar-2009
 *      Author: N. Mori, S. Ricciarini
 */

/*! @file TrkGeomCut.cpp The TrkGeomCut class implementation file */

#include "TrkGeomCut.h"

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] + xT < trkTrack->xGF[igf] && trkTrack->xGF[igf] < TrkParams::xGF_max[igf] - xT))
      return 0;
    if (!(TrkParams::yGF_min[igf] + yT < trkTrack->yGF[igf] && trkTrack->yGF[igf] < TrkParams::yGF_max[igf] - yT))
      return 0;

  } // end loop on GF planes

  return CUTOK;
}
1	/*
2	* TrkGeomCut.cpp
3	*
4	* Created on: 13-mar-2009
5	* Author: N. Mori, S. Ricciarini
6	*/
7
8	/! @file TrkGeomCut.cpp The TrkGeomCut class implementation file /
9
10	#include "TrkGeomCut.h"
11
12	int TrkGeomCut::Check(PamLevel2 *event) {
13
14	TrkTrack *trkTrack = event->GetTrack(0)->GetTrkTrack();
15
16	Float_t xT, yT;
17
18	for (int igf = 0; igf < TrkParams::nGF; igf++) { // loop on GF planes
19
20	if (igf==5 \|\| igf==10) { // M1 and M2 planes
21	xT = _xMagnTol;
22	yT = _yMagnTol;
23	}
24	else {
25	xT = _xTol;
26	yT = _yTol;
27	}
28
29	// check done for all 14 planes
30	if (!(TrkParams::xGF_min[igf] + xT < trkTrack->xGF[igf] && trkTrack->xGF[igf] < TrkParams::xGF_max[igf] - xT))
31	return 0;
32	if (!(TrkParams::yGF_min[igf] + yT < trkTrack->yGF[igf] && trkTrack->yGF[igf] < TrkParams::yGF_max[igf] - yT))
33	return 0;
34
35	} // end loop on GF planes
36
37	return CUTOK;
38	}