/[PAMELA software]/DarthVader/TrackerLevel2/src/TrkLevel2.cpp
ViewVC logotype

Diff of /DarthVader/TrackerLevel2/src/TrkLevel2.cpp

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

revision 1.49 by pam-fi, Tue Nov 25 14:41:37 2008 UTC revision 1.50 by pam-fi, Fri Dec 5 08:26:46 2008 UTC
# Line 400  Int_t TrkTrack::GetLeverArmY(){ Line 400  Int_t TrkTrack::GetLeverArmY(){
400      return (last_plane-first_plane+1);      return (last_plane-first_plane+1);
401  }  }
402  /**  /**
403     * Returns the track "lever-arm" on the x+y view, defined as the distance (in planes) between
404     * the upper and lower x,y (couple) measurements (the maximum value of lever-arm is 6).
405     */
406    Int_t TrkTrack::GetLeverArmXY(){
407        int first_plane = -1;
408        int last_plane  = -1;
409        for(Int_t ip=0; ip<6; ip++){
410            if( XGood(ip)*YGood(ip) && first_plane == -1 )first_plane = ip;
411            if( XGood(ip)*YGood(ip) && first_plane != -1 )last_plane = ip;
412        }
413        if( first_plane == -1 || last_plane == -1){
414            cout<< "Int_t TrkTrack::GetLeverArmXY() -- XGood(ip)*YGood(ip) always false ??? "<<endl;
415            return 0;
416        }
417        return (last_plane-first_plane+1);
418    }
419    /**
420   * Returns the reduced chi-square of track x-projection   * Returns the reduced chi-square of track x-projection
421   */   */
422  Float_t  TrkTrack::GetChi2X(){  Float_t  TrkTrack::GetChi2X(){
# Line 724  void TrkTrack::FillMiniStruct(cMini2trac Line 741  void TrkTrack::FillMiniStruct(cMini2trac
741              if(       XGood(i) && !YGood(i) ){              if(       XGood(i) && !YGood(i) ){
742                  track.xm_a[i] = xm[i] - omega * segment;                  track.xm_a[i] = xm[i] - omega * segment;
743                  track.ym_a[i] = ym[i] + segment;                  track.ym_a[i] = ym[i] + segment;
744                  track.zm_a[i] = zm[i] + beta * segment;  //          track.zm_a[i] = zm[i] + beta * segment;//not used yet
745                  track.xm_b[i] = xm[i] + omega * segment;                  track.xm_b[i] = xm[i] + omega * segment;
746                  track.ym_b[i] = ym[i] - segment;                  track.ym_b[i] = ym[i] - segment;
747                  track.zm_b[i] = zm[i] - beta * segment;  //          track.zm_b[i] = zm[i] - beta * segment;//not used yet
748              }else if( !XGood(i) && YGood(i) ){              }else if( !XGood(i) && YGood(i) ){
749                  track.xm_a[i] = xm[i] + segment;                  track.xm_a[i] = xm[i] + segment;
750                  track.ym_a[i] = ym[i] + omega * segment;                  track.ym_a[i] = ym[i] + omega * segment;
751                  track.zm_a[i] = zm[i] - gamma * segment;  //          track.zm_a[i] = zm[i] - gamma * segment;//not used yet
752                  track.xm_b[i] = xm[i] - segment;                  track.xm_b[i] = xm[i] - segment;
753                  track.ym_b[i] = ym[i] - omega * segment;                  track.ym_b[i] = ym[i] - omega * segment;
754                  track.zm_b[i] = zm[i] + gamma * segment;  //          track.zm_b[i] = zm[i] + gamma * segment;//not used yet
755              }              }
756          }          }
757                    
# Line 1259  Float_t TrkTrack::GetResidual_max(int ip Line 1276  Float_t TrkTrack::GetResidual_max(int ip
1276    
1277  };  };
1278  /**  /**
1279   * \brief Give the average spatial residual   * \brief Give the anerage spatial residual
1280   */   */
1281  Float_t TrkTrack::GetResidual_av(int ip, int iv){  Float_t TrkTrack::GetResidual_av(int ip, int iv){
1282      //      //
# Line 1446  void TrkSinglet::Dump(){ Line 1463  void TrkSinglet::Dump(){
1463      cout << endl << "plane        : " << plane;      cout << endl << "plane        : " << plane;
1464      cout << endl << "coord[2]     : "; while( i<2 && cout << coord[i] << " ") i++;      cout << endl << "coord[2]     : "; while( i<2 && cout << coord[i] << " ") i++;
1465      cout << endl << "sgnl         : " << sgnl;      cout << endl << "sgnl         : " << sgnl;
1466      cout << endl << "max.strip    : " << GetCluster_MaxStrip();      cout << endl << "max.strip    : ";
1467      cout << endl << "multiplicity : " << GetCluster_Multiplicity();      cout << endl << "multiplicity : ";
1468  }  }
1469  //--------------------------------------  //--------------------------------------
1470  //  //
Legend:
Removed from v.1.49  
changed lines
  Added in v.1.50
  ViewVC Help
Powered by ViewVC 1.1.23