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

Diff of /PamCut/TrkCuts/TrkGeomCut/TrkGeomCut.cpp

Parent Directory Parent Directory | Revision Log Revision Log | View Patch Patch

revision 1.1.1.1 by pam-fi, Wed May 27 13:30:09 2009 UTC revision 1.3 by pam-fi, Wed Mar 25 17:38:09 2015 UTC
# Line 2  Line 2 
2   * TrkGeomCut.cpp   * TrkGeomCut.cpp
3   *   *
4   *  Created on: 13-mar-2009   *  Created on: 13-mar-2009
5   *      Author: Sergio Ricciarini, Nicola Mori   *      Author: N. Mori, S. Ricciarini
6   */   */
7    
8  /*! @file TrkGeomCut.cpp The TrkGeomCut class implementation file */  /*! @file TrkGeomCut.cpp The TrkGeomCut class implementation file */
# Line 11  Line 11 
11    
12  int TrkGeomCut::Check(PamLevel2 *event) {  int TrkGeomCut::Check(PamLevel2 *event) {
13    
14    TrkTrack *trkTrack = event->GetTrack(0)->GetTrkTrack();  
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    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      // check done for all 14 planes
34      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))
35        return 0;        return 0;
36      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))
37        return 0;        return 0;
38    
39    } // end loop on GF planes    } // end loop on GF planes
Legend:
Removed from v.1.1.1.1  
changed lines
  Added in v.1.3
  ViewVC Help
Powered by ViewVC 1.1.23