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) {


    if(event->GetNTracks(_trkAlg)==0)return 0;

//  TrkTrack *trkTrack = event->GetTrack(0)->GetTrkTrack();
    ExtTrack *trkTrack = event->GetTrack(0,_trkAlg)->GetExtTrack();
  
  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	pam-fi	1.3
15			if(event->GetNTracks(_trkAlg)==0)return 0;
16
17			// TrkTrack *trkTrack = event->GetTrack(0)->GetTrkTrack();
18			ExtTrack *trkTrack = event->GetTrack(0,_trkAlg)->GetExtTrack();
19	pam-fi	1.2
20			Float_t xT, yT;
21	pam-fi	1.1
22			for (int igf = 0; igf < TrkParams::nGF; igf++) { // loop on GF planes
23
24	pam-fi	1.2	if (igf==5 \|\| igf==10) { // M1 and M2 planes
25			xT = _xMagnTol;
26			yT = _yMagnTol;
27			}
28			else {
29			xT = _xTol;
30			yT = _yTol;
31			}
32
33	pam-fi	1.1	// check done for all 14 planes
34	pam-fi	1.2	if (!(TrkParams::xGF_min[igf] + xT < trkTrack->xGF[igf] && trkTrack->xGF[igf] < TrkParams::xGF_max[igf] - xT))
35	pam-fi	1.1	return 0;
36	pam-fi	1.2	if (!(TrkParams::yGF_min[igf] + yT < trkTrack->yGF[igf] && trkTrack->yGF[igf] < TrkParams::yGF_max[igf] - yT))
37	pam-fi	1.1	return 0;
38
39			} // end loop on GF planes
40
41			return CUTOK;
42			}