/[PAMELA software]/PamCut/TrkCuts/TrkGeomCut/TrkGeomCut.cpp
ViewVC logotype

Contents of /PamCut/TrkCuts/TrkGeomCut/TrkGeomCut.cpp

Parent Directory Parent Directory | Revision Log Revision Log


Revision 1.3 - (show annotations) (download)
Wed Mar 25 17:38:09 2015 UTC (9 years, 8 months ago) by pam-fi
Branch: MAIN
CVS Tags: HEAD
Changes since 1.2: +5 -1 lines
Port to 10th reduction.

Some cuts have been modified to use the extended track objects introduced with 10th reduction software. Some others have still to be ported.
Git ref: 4bd0c70baf73c56bddc08dca365cbff243ad8bc8
(Nicola)

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
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
20 Float_t xT, yT;
21
22 for (int igf = 0; igf < TrkParams::nGF; igf++) { // loop on GF planes
23
24 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 // check done for all 14 planes
34 if (!(TrkParams::xGF_min[igf] + xT < trkTrack->xGF[igf] && trkTrack->xGF[igf] < TrkParams::xGF_max[igf] - xT))
35 return 0;
36 if (!(TrkParams::yGF_min[igf] + yT < trkTrack->yGF[igf] && trkTrack->yGF[igf] < TrkParams::yGF_max[igf] - yT))
37 return 0;
38
39 } // end loop on GF planes
40
41 return CUTOK;
42 }

  ViewVC Help
Powered by ViewVC 1.1.23