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	pam-fi	1.1	/*
2			* TrkGeomCut.cpp
3			*
4			* Created on: 13-mar-2009
5	pam-fi	1.2	* Author: N. Mori, S. Ricciarini
6	pam-fi	1.1	*/
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	pam-fi	1.2
16			Float_t xT, yT;
17	pam-fi	1.1
18			for (int igf = 0; igf < TrkParams::nGF; igf++) { // loop on GF planes
19
20	pam-fi	1.2	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	pam-fi	1.1	// check done for all 14 planes
30	pam-fi	1.2	if (!(TrkParams::xGF_min[igf] + xT < trkTrack->xGF[igf] && trkTrack->xGF[igf] < TrkParams::xGF_max[igf] - xT))
31	pam-fi	1.1	return 0;
32	pam-fi	1.2	if (!(TrkParams::yGF_min[igf] + yT < trkTrack->yGF[igf] && trkTrack->yGF[igf] < TrkParams::yGF_max[igf] - yT))
33	pam-fi	1.1	return 0;
34
35			} // end loop on GF planes
36
37			return CUTOK;
38			}