/[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.58 by mocchiut, Wed Aug 27 14:44:37 2014 UTC revision 1.60 by pam-fi, Thu Jan 28 14:38:24 2016 UTC
# Line 381  Int_t TrkTrack::GetLeverArmXY(){ Line 381  Int_t TrkTrack::GetLeverArmXY(){
381      int first_plane = -1;      int first_plane = -1;
382      int last_plane  = -1;      int last_plane  = -1;
383      for(Int_t ip=0; ip<6; ip++){      for(Int_t ip=0; ip<6; ip++){
384          if( XGood(ip)*YGood(ip) && first_plane == -1 )first_plane = ip;          if( XGood(ip) && YGood(ip) && first_plane == -1 )first_plane = ip;
385          if( XGood(ip)*YGood(ip) && first_plane != -1 )last_plane = ip;          if( XGood(ip) && YGood(ip) && first_plane != -1 )last_plane = ip;
386      }      }
387      if( first_plane == -1 || last_plane == -1){      if( first_plane == -1 || last_plane == -1){
388          cout<< "Int_t TrkTrack::GetLeverArmXY() -- XGood(ip)*YGood(ip) always false ??? "<<endl;          cout<< "Int_t TrkTrack::GetLeverArmXY() -- XGood(ip)*YGood(ip) always false ??? "<<endl;
# Line 939  void TrkTrack::Fit(double pfixed, int& f Line 939  void TrkTrack::Fit(double pfixed, int& f
939          for(int i=0; i<5; i++) al[i]=al_ini[i];          for(int i=0; i<5; i++) al[i]=al_ini[i];
940      }      }
941    
942    
943        //ELENA 2015
944        int ngf = TrkParams::nGF;
945        float* zgf =  TrkParams::zGF;    
946        Trajectory tj = Trajectory(ngf,zgf);
947        tj.DoTrack(al,ZINI);
948        for(int i=0; i<14; i++){
949            xGF[i] = tj.x[i];
950            yGF[i] = tj.y[i];
951        }
952        if(ngf!=14){
953            cout << "TrkTrack::Fit("<<pfixed<<","<<fail<<","<<iprint<<","<<froml1<<"<<) "<<endl;;
954            cout << "TrkParams::nGF = "<<TrkParams::nGF<<" != 14 "<<endl;
955            cout << "back-incompatibility !?!?! acceptance check not reliable"<<endl;
956        }
957        
958        
959    
960  };  };
961  /**  /**
962   * Reset the fit parameters   * Reset the fit parameters
# Line 1584  void TrkLevel2::Set(){ Line 1602  void TrkLevel2::Set(){
1602  //  //
1603  //  //
1604  //--------------------------------------  //--------------------------------------
1605    void TrkLevel2::SetTrackArray(TClonesArray *track){
1606        if(track && strcmp(track->GetClass()->GetName(),"TrkTrack")==0){
1607            if(Track)Track->Clear("C");    
1608            Track = track;
1609        }
1610    }
1611    //--------------------------------------
1612    //
1613    //
1614    //--------------------------------------
1615  void TrkLevel2::Dump(){  void TrkLevel2::Dump(){
1616                    
1617          //          //
Legend:
Removed from v.1.58  
changed lines
  Added in v.1.60
  ViewVC Help
Powered by ViewVC 1.1.23